CN108956097B - Light field polarization state measuring method and device, computer equipment and storage medium - Google Patents
Light field polarization state measuring method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN108956097B CN108956097B CN201810749870.7A CN201810749870A CN108956097B CN 108956097 B CN108956097 B CN 108956097B CN 201810749870 A CN201810749870 A CN 201810749870A CN 108956097 B CN108956097 B CN 108956097B
- Authority
- CN
- China
- Prior art keywords
- polarization state
- target
- probability density
- density distribution
- mueller matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0292—Testing optical properties of objectives by measuring the optical modulation transfer function
Abstract
The application relates to a method and a device for measuring the polarization state of a light field, computer equipment and a storage medium. The method comprises the following steps: acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to a maximum probability value in a target polarization state probability density distribution function; and representing the polarization state of the light field to be measured by the target Stokes parameters. By adopting the method, the measurement precision of the polarization state of the light field can be improved under weak light, and the measurement speed of the polarization state of the light field can be improved under non-weak light.
Description
Technical Field
The present application relates to the field of optical technologies, and in particular, to a method and an apparatus for measuring a polarization state of an optical field, a computer device, and a storage medium.
Background
With the development of the optical field polarization measurement technology, a theoretical technology for detecting probability density of the polarization state of the optical field appears. The technology utilizes a phase modulation device in spatial distribution to generate a Point Spread Function (PSF) image related to polarization, processes the PSF image under a Bayesian statistical theory framework, and calculates to obtain a polarization state probability density distribution function of a light field to be measured.
However, in the conventional method, although a theoretical relationship between the point spread function of the optical field polarization state probability density detection and the pinhole imaging is established, the theoretical jones matrix of the specific phase modulator in the detection system is not shown in the theoretical relationship, and the theoretical relationship can be established only under the condition that the specific phase modulator selected by experimental measurement is ideal and error-free. Since the jones matrix of the phase modulator cannot be directly measured, it is usually necessary to use the ideal jones matrix of the designed phase modulation device in the test data processing. Therefore, the processing error of the phase modulation device in the detection system and the error of the detection system cannot be calibrated through experiments, and the measurement precision of the polarization state of the light field is not high enough.
Disclosure of Invention
In view of the above, it is necessary to provide a light field polarization state measurement method, a light field polarization state measurement device, a computer device, and a storage medium, which can improve the precision of light field polarization state measurement.
A method of light field polarization state measurement, the method comprising: acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; the calibration incident light is output through a phase modulator to be tested to obtain calibration incident light; calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function; and expressing the polarization state of the light field to be measured through the target Stokes parameter.
In one embodiment, the calculating the target mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data includes: acquiring first Mueller matrices of a plurality of groups of polarization state generators and second Mueller matrices of a polarization state analyzer and target light intensities corresponding to each group of the first Mueller matrices and the second Mueller matrices according to the incident polarization state data and the emergent polarization state data; calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting to obtain a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
In one embodiment, the inputting the target mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target mueller matrix includes: acquiring a point spread function light intensity distribution image corresponding to the light field to be detected; extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, the plurality of spatial position information and photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model; and converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
In one embodiment, the obtaining of the target stokes parameter corresponding to the maximum probability value in the target probability density distribution function includes: taking logarithm of the target polarization state probability density distribution function to obtain a three-dimensional function on a Ponga spherical coordinate system; selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on the Bonga spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter; and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
In one embodiment, the model of probability density distribution of polarization state of the target includes:
wherein the content of the first and second substances,U1=∫u1(X)d2X,and the target Mueller matrix of the phase modulator to be tested isXiIndicates the ith spatial position information (common I)iA spatial location); n is a radical ofiRepresenting and i-th spatial position information XiCorresponding photon number information;and representing the target Stokes parameters corresponding to the light field to be measured.
An optical field polarization state measurement apparatus, the apparatus comprising: the polarization state generator is used for generating incident polarization state data of the calibration incident light; the polarization state analyzer is used for generating emergent polarization state data of the calibration emergent light; the calibration incident light is output through a phase modulator to be tested to obtain calibration incident light; the processor is used for calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function; and expressing the polarization state of the light field to be measured through the target Stokes parameter.
In one embodiment, the polarization state generator and the polarization state analyzer are further configured to generate a first mueller matrix of the plurality of sets of polarization state generators and a second mueller matrix of the polarization state analyzer, and the apparatus further includes a detector configured to obtain a target light intensity corresponding to each set of the first mueller matrices and the second mueller matrices; the processor is further configured to calculate an instrument matrix corresponding to the light field to be measured according to the first mueller matrix and the second mueller matrix, and convert the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
In one embodiment, the polarization analyzer is used as a left-handed analyzer to modulate the light field to be measured; the detector is also used for acquiring a point spread function light intensity distribution image corresponding to the light field to be detected; the processor is further used for extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, the plurality of spatial position information and photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model; and performing conversion processing on the input polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
A computer device comprising a memory storing a computer program and a processor implementing the steps of the light field polarization state measurement method described in the various embodiments above when the computer program is executed.
A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps of the light field polarization state measurement method described in the various embodiments above.
According to the light field polarization state measuring method, the light field polarization state measuring device, the computer equipment and the storage medium, the target Mueller matrix of the phase modulator to be measured is obtained through obtaining the incident polarization state data of the calibration incident light and the emergent polarization state data of the calibration emergent light, and according to the incident polarization state data and the emergent polarization state data, the on-line calibration of the phase modulator to be measured is achieved. By providing a pre-constructed polarization state probability density distribution model, a target polarization state probability density distribution function reflecting Stokes parameters and corresponding probability values can be obtained according to the calibrated actual target Mueller matrix of the phase modulator to be measured. And representing the polarization state of the light field to be measured by obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function. Through on-line calibration of the phase modulator to be measured, the actual target Mueller matrix of the phase modulator to be measured in the light field to be measured can be measured, and through providing a pre-constructed polarization state probability density distribution model, the problems that the Jones matrix of a phase modulator device is hidden in an existing theoretical model, and the system error cannot be calibrated due to the fact that the Jones matrix of the device cannot be measured are solved, and therefore the precision of light field polarization state measurement is improved.
Drawings
FIG. 1 is a diagram illustrating an exemplary application of a method for measuring polarization state of a light field;
FIG. 2 is a schematic flow chart of a method for measuring polarization state of a light field according to an embodiment;
FIG. 3 is a schematic flow chart of a method for measuring polarization state of light field in another embodiment;
FIG. 4 shows an exemplary embodiment of sixteen point spread function light intensity distribution images collected experimentally;
FIG. 5 is a simulation image of sixteen normalized light intensity distributions in one embodiment;
FIG. 6 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The method for measuring the polarization state of the light field provided by the application can be applied to the application environment shown in fig. 1. Referring to fig. 1, the application environment includes a light source 102, a spatial filter 104, a collimating lens 106, a polarization state generator 108, a phase modulator 110 to be tested, a polarization state analyzer 112, a detector 114, and a terminal 116. Wherein, the light source 102 may be a he — ne laser. The spatial filter 104 may be constituted by a microscope objective and an aperture. The polarization state generator 108 may be composed of a polarizing plate P1And lambda/4 wave plate Q1The polarization analyzer 112 may be composed of a polarizing plate P2And lambda/4 wave plate Q2Is composed ofThe polarization state generator 108 is structurally symmetrical to the polarization state analyzer 112. The terminal 116 is in communication with the polarization state generator 108, the polarization state analyzer 112, and the detector 114, respectively, via a network. The terminal 114 may be, but is not limited to, various personal computers, notebook computers, smart phones, tablet computers, and portable wearable devices.
After the light source 102 generates the initial light beam, the light beam is converged to the aperture through the microscope objective in the spatial filter 104, and a spherical wave is diffracted from the aperture, so as to filter out high-frequency components and ensure high energy utilization rate. The collimating lens 106 collimates the spherical wave diffracted by the aperture into a parallel beam for further measurement. The parallel light beam passes through the polarization state generator 108 to obtain the incident polarization state data of the calibration incident light. The outgoing polarization state data of the calibration outgoing light can be obtained after the calibration incoming light passes through the phase modulator 110 to be tested and the polarization state analyzer 112. The terminal 116 may obtain the incident polarization state data and the emergent polarization state data, and calculate to obtain a target mueller matrix of the phase modulator 110 to be measured. When the terminal 116 inputs the target mueller matrix into the pre-constructed polarization state probability density distribution model, a target polarization state probability density distribution function corresponding to the target mueller matrix is obtained through output. The terminal 116 obtains a target stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function, and represents the polarization state of the light field to be measured in the application environment through the target stokes parameter.
In one embodiment, as shown in fig. 2, a method for measuring polarization state of optical field is provided, which is exemplified by the application of the method to the terminal 116 in fig. 1, and includes the following steps:
step 202, acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; and outputting the calibration incident light through the phase modulator to be tested to obtain calibration incident light.
The calibration incident light refers to a light beam for calibrating the phase modulator to be measured, and may be a light beam which passes through the polarization state generator and enters the phase modulator to be measured. The incident polarization state data refers to polarization state data of the nominal incident light. The calibration light is a light beam obtained by outputting the calibration incident light through the phase modulator to be measured and the polarization analyzer. The emergent polarization state data refers to polarization state data of the calibrated emergent light. The incident polarization state data includes, but is not limited to, a mueller matrix of the polarization state generator, and the emergent polarization state data includes, but is not limited to, a mueller matrix of the polarization state analyzer. The parallel light beam is converted into a nominal incident light having a specific polarization state after passing through the polarization state generator. After the calibration incident light passes through the phase modulator to be tested, calibration emergent light of another polarization state can be obtained at the polarization state analyzer.
The phase modulator to be measured refers to a phase modulator which needs to be calibrated in an actual Mueller matrix. A phase modulator refers to a device that can phase modulate a light beam. The phase modulator under test may be a specific phase modulator, such as a triple stress suppressing device SEO, or may be an arbitrary phase modulator PM. The phase modulator to be measured can be obtained by processing optical glass serving as a raw material by a method of applying triple symmetric pressure to the optical glass, wherein the triple symmetric pressure can enable the interior of the optical glass to generate triple symmetric stress birefringence.
And 204, calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data.
In one embodiment, the calculating the target mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data includes: acquiring a first Mueller matrix of a plurality of groups of polarization state generators and a second Mueller matrix of a polarization state analyzer according to the incident polarization state data and the emergent polarization state data, and acquiring target light intensity corresponding to each group of the first Mueller matrices and the second Mueller matrices; calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
Assuming that the number of measurements is Q, in the Q-th (Q-1, 2, …, Q) -th measurement, the illumination beam S0Transmission polarization state generationAfter the detector, the phase modulator to be detected and the polarization state analyzer reach the target light intensity I on the detection surface of the detectorqComprises the following steps:
in the above formula, P2Andrespectively represent a polarizing plate P2And lambda/4 wave plate Q2Due to the polarizer P2And lambda/4 wave plate Q2Form a polarization state analyzer, then P2Anda second mueller matrix collectively forming a polarization analyzer; p1Andrespectively represent a polarizing plate P1And lambda/4 wave plate Q1Due to the polarizer P1And lambda/4 wave plate Q1Form a polarization state generator, then P1Anda first mueller matrix collectively forming a polarization state generator; mSEOThe target mueller matrix of the phase modulator to be measured is obtained. Order:
the following can be obtained:
rewriting the above formula to a matrix form:
target Mueller matrix MSEOCorresponding column vector M'SEO=[m00m01… m33]TQ sets of target intensities I may be used after a total of Q measurementsqConversion to light intensity vector I:
wherein, WMAnd measuring the instrument matrix of the Mueller matrix of the phase modulator to be measured for the Mueller detection system. The instrument matrix is obtained by converting a first Mueller matrix of the Q-group polarization state generator and a second Mueller matrix of the polarization state analyzer. From the matrix calculation, M 'was obtained by the following equation'SEO:
M'SEO=WM +·I
In the above formula, WM +Is an instrument matrix WMWhen W is a generalized inverse matrix ofMWhen column is full rank, WM +From WMUniquely determined and can be derived from a least squares estimate:
WM +=(WM TWM)-1WM T
if necessary, measure M'SEOThe first mueller matrices of at least 16 sets of polarization state generators and the second mueller matrices of the polarization state analyzer, and the target light intensities corresponding to each set of the first mueller matrices and the second mueller matrices can be obtained through at least 16 measurements. For example, by selecting a/4 wave plate Q1And Q2The azimuth angle (the included angle between the fast axis and the x axis) of the array makes the instrument matrix WMIs equal to 16. Obtained when measured to give M'SEOWhen all 16 elements are contained, then M 'can be obtained'SEOAnd determining a target Mueller matrix of the phase modulator to be tested.
In one embodiment, the selection can be performed bySixteen sets of lambda/4 wave plates Q1And Q2Such that the instrument matrix WMIs 16. The experimental sixteen groups of wave plates have four Q azimuth angles1Azimuth and four Q2The azimuth angles are combined into a group two by two, and the sixteen groups of wave plates are combined with the polarizing plate P1、P2The sixteen polarization states produced are linearly independent of each other. For example, a first mueller matrix of the 16 sets of polarization state generators and a second mueller matrix of the polarization state analyzer, and a target light intensity corresponding to each set of the first mueller matrix and the second mueller matrix may be obtained by the sixteen sets of wave plate azimuth angles shown in table 1 below:
TABLE 1
Step 206, inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to the light field to be measured.
The left-handed analyzer (L HC) is composed of a lambda/4 wave plate and a linear polarizer, and is used for passing left-handed rotation and blocking right-handed rotation.
The pre-constructed polarization state probability density distribution model comprises a target Mueller matrix of the phase modulator to be tested, space position information in a point spread function light intensity distribution image and corresponding photon number information. The point spread function light intensity distribution image is generated by emergent light to be measured after light beams corresponding to a light field to be measured pass through the phase modulator to be measured and the left-handed polarization analyzer. And substituting specific spatial position information and corresponding photon number information of a point spread function light intensity distribution image measured by a detector into a pre-constructed polarization state probability density distribution model to obtain a target polarization state probability density distribution function. The target polarization state probability density distribution function reflects the mapping relation between the Stokes parameters and the corresponding probability values. The probability value refers to the probability of occurrence of the polarization state of the light field to be measured reflected by the Stokes parameter. The polarization state of the light field to be measured can be determined more accurately by calibrating the phase modulator to be measured and inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model.
Step 208, obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function; and representing the polarization state of the light field to be measured by the target Stokes parameters.
When the detector detects a point spread function light intensity distribution image, that is, when the spatial position information of N photons is detected, there are many possibilities for the N photons to have such an incident polarization state. In other words, it is theoretically possible for the different incident polarization states to cause the N photons to fall on the detection plane X at such specific N spatial points XnThe above. However, the probability that a photon will behave in this way with a different incident polarization state is different. The actually occurring event is the event with the highest probability among all possible events, and the observed event can be regarded as the event with the highest probability density, as defined in quantum theory. Therefore, the stokes parameter corresponding to the maximum value in the polarization state probability density distribution function can be used as the target stokes parameter of the light field to be detected corresponding to the detected point spread function light intensity distribution image. And the target Stokes parameters can be used for representing the polarization state of the light field to be measured.
In the light field polarization state method, the target Mueller matrix of the phase modulator to be measured is obtained by obtaining the incident polarization state data of the calibration incident light and the emergent polarization state data of the calibration emergent light and calculating according to the incident polarization state data and the emergent polarization state data, so that the on-line calibration of the phase modulator to be measured is realized. By providing a pre-constructed polarization state probability density distribution model, a target polarization state probability density distribution function reflecting Stokes parameters and corresponding probability values can be obtained according to the calibrated actual target Mueller matrix of the phase modulator to be measured. And representing the polarization state of the light field to be measured by obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function. Through on-line calibration of the phase modulator to be measured, the actual target Mueller matrix of the phase modulator to be measured in the light field to be measured can be measured, and through providing a pre-constructed polarization state probability density distribution model, the problems that the Jones matrix of a phase modulator device is hidden in an existing theoretical model, and the system error cannot be calibrated due to the fact that the Jones matrix of the device cannot be measured are solved, and therefore the precision of light field polarization state measurement is improved.
In one embodiment, inputting a target mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting a target polarization state probability density distribution function corresponding to the target mueller matrix, includes: acquiring a point spread function light intensity distribution image corresponding to a light field to be detected; extracting a plurality of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each spatial position information; inputting the target Mueller matrix, a plurality of pieces of spatial position information and photon number information corresponding to each piece of spatial position information into a pre-constructed polarization state probability density distribution model; and converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
The method includes the steps of analyzing a point spread function light intensity distribution image to obtain a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information, for example, generating a light intensity distribution simulation image according to the point spread function light intensity distribution image, and extracting the plurality of pieces of spatial position information in the light intensity distribution simulation image and the photon number information corresponding to each piece of spatial position information through a Mat L ab program.
In one embodiment, the model of the probability density distribution of the polarization state of the object comprises:
wherein the content of the first and second substances,U1=∫u1(X)d2X,and the target Mueller matrix of the phase modulator to be tested isXiIndicates the ith spatial position information (common I)iA spatial location); n is a radical ofiRepresenting and i-th spatial position information XiCorresponding photon number information;and representing the target Stokes parameters corresponding to the light field to be measured.
The general expression of the target Mueller matrix of the phase modulator to be tested isAfter the polarization state generator generates the light field to be measured, the polarization state of the light field to be measured can be represented by the standardized stokes parameters:the stokes parameter of the initial emergent light obtained after the light beam corresponding to the light field to be detected passes through the phase modulator to be detected can be expressed as follows:
the Mueller matrix of the wave plate of the left-handed polarization analyzer and the linear polarizer is respectively as follows:
the Mueller matrix of the left-handed polarization analyzer is the product of the wave plate and the Mueller matrix of the linear polarizer:
after the initial emergent light passes through the left-handed polarization analyzer, the Stokes parameters of the emergent light to be measured of the light field to be measured can be obtainedComprises the following steps:
the first component of the stokes parameter is the total light intensity of the light field to be measured. The distribution of the light intensity over the detection surface of the detector is then:
the light intensity distribution model in the conventional method is only established for the special phase modulator SEO, but is established for the actual polarization matrix of the phase modulator caused by the defects of the phase modulator, the errors of the experimental system and the likeThe change of (1) can only obtain the model parameter I represented by the light intensity distribution image corresponding to the left-handed and right-handed circularly polarized incident lights through experimental measurement in the traditional modeR、ILFor model IuAnd IR、ILThe relationship between them can still only be modeled in the ideal case. The light intensity distribution in the above embodiment may change with the change of the actual target mueller matrix of the phase modulator to be measured, so as to obtain a light intensity distribution model conforming to the actual light intensity distribution.
The above formula can be rewritten as:
wherein the content of the first and second substances,and (m)01-m31)、 (m02-m32)、(m03-m33) Is not equal to zero; x represents spatial position information.
The probability that a photon falls on any point X of the detection surface is as follows:
The model of the probability density distribution function of the polarization state at this time is as follows:
for a point spread function light intensity distribution image measured by experimentExtracting a plurality of spatial position information X in a point spread function light intensity distribution imagei(assuming common I)iIndividual spatial point) and the ith spatial position information XiCorresponding photon number information Ni(ii) a According to a plurality of spatial position information XiAnd with each spatial position information XiCorresponding photon number information NiThe probability density distribution function for the target polarization state can be obtained by:
in one embodiment, obtaining a target stokes parameter corresponding to a maximum probability value in a target polarization state probability density distribution function includes: taking logarithm of the target polarization state probability density distribution function to obtain a three-dimensional function on a Ponga spherical coordinate system; selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on a Pongall spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter; and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
Due to the number of spatial positions I in the actual measurementiPhoton number information N corresponding to each spatial positioniVery large, resulting in that the value of several in the probability density distribution function of the target polarization state will approach zero, which is not beneficial for subsequent data processing. Thus, the probability density distribution function for the target polarization state can be logarithmized instead of the target polarization stateAnd (3) carrying out data processing analysis:
since the log function is a monotonically increasing function, it does not affectDistribution trend and maximum value point value. The maximum value point(s) of the logarithmic target polarization state probability density distribution function can be obtained by calculating and processing the logarithmic target polarization state probability density distribution functionmax,logPmax). Wherein s ismaxRepresenting the target Stokes parameter, logP, corresponding to the maximum value of the severalmaxRepresenting the logarithm of the maximum value of the several degrees.
Is a three-dimensional function on a poincare spherical coordinate system. The Poincare sphere is a graphical representation of either polarization state. Since any elliptically polarized light only needs two azimuth angles to completely determine the polarization state, and the two azimuth angles can be represented by longitude and latitude on a spherical surface, one point on the spherical surface can represent one polarization state, and the combination of all the points on the Poincall represents all the possible polarization states. For example, at the Poincare sphere center, completely unpolarized light is represented; on the Poincare sphere, fully polarized light is represented; at any point within the Poincare sphere, partially polarized light is represented. The three-dimensional functions on the Pongar ball coordinate system can be input according to the test Stokes parameters corresponding to the coordinate points at every preset position on the selected Pongar ball, or the three-dimensional functions on the Pongar ball coordinate system with the test Stokes parameters exceeding the preset number can be generated through a random algorithm. By flexibly controlling the precision of the determined test Stokes parameters, the accuracy and the efficiency of the maximum probability value of the three-dimensional function on the calculated Bonga spherical coordinate system can be improved.
In one embodiment, preset point coordinate data, and first straight line data of a first straight line, second straight line data of a second straight line, and third straight line data of a third straight line corresponding to the point coordinate data may also be taken from a poincare coordinate system; the first straight line is parallel to the x axis of the Poincare coordinate system, the second straight line is parallel to the y axis of the Poincare coordinate system, and the second straight line is parallel to the z axis of the Poincare coordinate system; and inputting the first straight line data, the second straight line data and the third straight line data into the three-dimensional function respectively, and outputting to obtain a corresponding first target polarization state probability density distribution function, a corresponding second target polarization state probability density distribution function and a corresponding third target polarization state probability density distribution function. The maximum point on the three-dimensional function on the Poincare sphere can be observed more intuitively through the first target polarization state probability density distribution function, the second target polarization state probability density distribution function and the third target polarization state probability density distribution function.
The light field polarization state measurement method and the provided polarization state probability density distribution model in the above embodiments can also be applied to weak light polarization measurement and quantum entanglement measurement.
In one embodiment, as shown in fig. 3, another optical field polarization state measurement method is provided, which is described by taking the method as an example applied to the terminal 116 in fig. 1, and includes the following steps:
And step 304, acquiring first Mueller matrices of the multiple groups of polarization state generators and second Mueller matrices of the polarization state analyzers according to the incident polarization state data and the emergent polarization state data, and acquiring target light intensities corresponding to the first Mueller matrices and the second Mueller matrices of each group.
And step 306, calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting to obtain a generalized inverse matrix of the instrument matrix.
And step 308, converting the acquired multiple groups of target light intensities into light intensity vectors.
And 310, multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
And step 312, acquiring a point spread function light intensity distribution image corresponding to the light field to be detected.
As shown in fig. 4, sixteen point spread function light intensity distribution images were collected for the experiment. As shown in fig. 5, the mueller matrix M 'of the phase modulator to be measured is obtained after data processing is performed on the collected sixteen images'SEOImages using normalizationThe light intensity distribution simulation image shows that white and black represent positive and negative values, respectively.
In step 314, a plurality of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each spatial position information are extracted.
And step 316, inputting the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model.
The model of probability density distribution of target polarization state comprises:
wherein the content of the first and second substances,U1=∫u1(X)d2X,and the target Mueller matrix of the phase modulator to be tested isXiIndicates the ith spatial position information (common I)iA spatial location); n is a radical ofiRepresenting and i-th spatial position information XiCorresponding photon number information;and representing the target Stokes parameters corresponding to the light field to be measured.
And step 318, performing conversion processing on the polarization state probability density distribution model after the target mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target mueller matrix.
The abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to the light field to be measured.
And step 320, taking logarithm of the probability density distribution function of the target polarization state to obtain a three-dimensional function on the Poincare spherical coordinate system.
And 322, selecting a plurality of testing Stokes parameters according to preset precision, inputting the testing Stokes parameters into a three-dimensional function on a Boncard spherical coordinate system, and obtaining a testing probability value corresponding to each testing Stokes parameter.
Step 324, using the test stokes parameter corresponding to the obtained maximum test rate value as a target stokes parameter; and representing the polarization state of the light field to be measured by the target Stokes parameters.
In the light field polarization state method, the target Mueller matrix of the phase modulator to be measured is obtained by obtaining the incident polarization state data of the calibration incident light and the emergent polarization state data of the calibration emergent light and calculating according to the incident polarization state data and the emergent polarization state data, so that the on-line calibration of the phase modulator to be measured is realized. By providing a pre-constructed polarization state probability density distribution model, a target polarization state probability density distribution function reflecting Stokes parameters and corresponding probability values can be obtained according to the calibrated actual target Mueller matrix of the phase modulator to be measured. And representing the polarization state of the light field to be measured by obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function. Through online calibration of the phase modulator to be measured, the actual target Mueller matrix of the phase modulator to be measured in the light field to be measured can be measured, and through providing a pre-constructed polarization state probability density distribution model, the problems that the Jones matrix of a phase modulator device is hidden by an existing theoretical model, and the system error cannot be calibrated due to the fact that the Jones matrix of the device cannot be measured are solved. And the accuracy and the efficiency of the maximum probability value of the three-dimensional function on the calculated Ponga spherical coordinate system can be improved by flexibly controlling the precision of the determined test Stokes parameters, so that the precision of the measurement of the light field polarization state is improved.
It should be understood that although the steps in the flowcharts of fig. 2 and 3 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 2 and 3 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performing the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.
In one embodiment, there is provided an optical field polarization state measuring apparatus, comprising: the polarization state generator is used for generating incident polarization state data of the calibration incident light; the polarization state analyzer is used for generating emergent polarization state data of the calibration emergent light; outputting the calibration incident light through the phase modulator to be tested to obtain calibration incident light; the processor is used for calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to a maximum probability value in a target polarization state probability density distribution function; and representing the polarization state of the light field to be measured by the target Stokes parameters.
In one embodiment, the polarization state generator and the polarization state analyzer are further configured to generate a first mueller matrix of the plurality of sets of polarization state generators and a second mueller matrix of the polarization state analyzer, and the apparatus further includes a detector configured to obtain a target light intensity corresponding to each set of the first mueller matrices and the second mueller matrices; the processor is also used for calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
In one embodiment, the polarization analyzer is used as a left-handed analyzer for modulating the light field to be measured; the detector is also used for acquiring a point spread function light intensity distribution image corresponding to the light field to be detected; the processor is also used for extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, a plurality of pieces of spatial position information and photon number information corresponding to each piece of spatial position information into a pre-constructed polarization state probability density distribution model; and performing conversion processing on the input polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
In one embodiment, the processor is further configured to logarithm the target polarization state probability density distribution function to obtain a three-dimensional function on the poincare sphere coordinate system; selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on a Pongall spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter; and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
For specific limitations of the optical field polarization state measurement apparatus, reference may be made to the above limitations of the optical field polarization state measurement method, which are not described herein again. The modules in the light field polarization state measuring device can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a terminal, such as terminal 116 in FIG. 1, whose internal structure diagram may be as shown in FIG. 6. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a light field polarization state measurement method. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 6 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, there is provided a computer device comprising a memory storing a computer program and a processor implementing the following steps when the processor executes the computer program: acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; outputting the calibration incident light through the phase modulator to be tested to obtain calibration incident light; calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to a maximum probability value in a target polarization state probability density distribution function; and representing the polarization state of the light field to be measured by the target Stokes parameters.
In one embodiment, the step of calculating the target mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data when the processor executes the computer program includes the following steps: acquiring a first Mueller matrix of a plurality of groups of polarization state generators and a second Mueller matrix of a polarization state analyzer according to the incident polarization state data and the emergent polarization state data, and acquiring target light intensity corresponding to each group of the first Mueller matrices and the second Mueller matrices; calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
In one embodiment, when the processor executes the computer program, the step of inputting the target mueller matrix into a pre-constructed polarization state probability density distribution model and outputting a target polarization state probability density distribution function corresponding to the target mueller matrix includes the following steps: acquiring a point spread function light intensity distribution image corresponding to a light field to be detected; extracting a plurality of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each spatial position information; inputting the target Mueller matrix, a plurality of pieces of spatial position information and photon number information corresponding to each piece of spatial position information into a pre-constructed polarization state probability density distribution model; and converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
In one embodiment, the step of obtaining the target stokes parameter corresponding to the maximum probability value in the target probability density distribution function when the processor executes the computer program includes the following steps: taking logarithm of the target polarization state probability density distribution function to obtain a three-dimensional function on a Ponga spherical coordinate system; selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on a Pongall spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter; and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
In one embodiment, the computer program, when executed by a processor, implements a model of probability density distribution of a target polarization state comprising:
wherein the content of the first and second substances,U1=∫u1(X)d2X,and the target Mueller matrix of the phase modulator to be tested isXiIndicates the ith spatial position information (common I)iA spatial location); n is a radical ofiRepresenting and i-th spatial position information XiCorresponding photon number information;and representing the target Stokes parameters corresponding to the light field to be measured.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of: acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; outputting the calibration incident light through the phase modulator to be tested to obtain calibration incident light; calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to a maximum probability value in a target polarization state probability density distribution function; and representing the polarization state of the light field to be measured by the target Stokes parameters.
In one embodiment, the step of calculating the target mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data when the computer program is executed by the processor includes the following steps: acquiring a first Mueller matrix of a plurality of groups of polarization state generators and a second Mueller matrix of a polarization state analyzer according to the incident polarization state data and the emergent polarization state data, and acquiring target light intensity corresponding to each group of the first Mueller matrices and the second Mueller matrices; calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
In one embodiment, the step of inputting a target mueller matrix into a pre-constructed polarization state probability density distribution model and outputting a target polarization state probability density distribution function corresponding to the target mueller matrix when the computer program is executed by the processor includes the following steps: acquiring a point spread function light intensity distribution image corresponding to a light field to be detected; extracting a plurality of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each spatial position information; inputting the target Mueller matrix, a plurality of pieces of spatial position information and photon number information corresponding to each piece of spatial position information into a pre-constructed polarization state probability density distribution model; and converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
In one embodiment, the step of obtaining the target stokes parameter corresponding to the maximum probability value in the target probability density distribution function when the computer program is executed by the processor includes the following steps: taking logarithm of the target polarization state probability density distribution function to obtain a three-dimensional function on a Ponga spherical coordinate system; selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on a Pongall spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter; and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
In one embodiment, the computer program, when executed by the processor, implements a model of probability density distribution of a target polarization state comprising:
wherein the content of the first and second substances,U1=∫u1(X)d2X,and the target Mueller matrix of the phase modulator to be tested isXiIndicates the ith spatial position information (common I)iA spatial location); n is a radical ofiRepresenting and i-th spatial position information XiCorresponding photon number information;and representing the target Stokes parameters corresponding to the light field to be measured.
It will be understood by those of ordinary skill in the art that all or a portion of the processes of the methods of the embodiments described above may be implemented by a computer program that may be stored on a non-volatile computer-readable storage medium, which when executed, may include the processes of the embodiments of the methods described above, wherein any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (8)
1. A method of light field polarization state measurement, the method comprising:
acquiring incident polarization state data of calibration incident light and emergent polarization state data of calibration emergent light; the calibration incident light is output through a phase modulator to be tested to obtain calibration incident light;
calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data;
inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix, wherein the specific steps are as follows:
acquiring a point spread function light intensity distribution image corresponding to the light field to be detected; extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, the plurality of spatial position information and photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model; converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix;
the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected;
obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function;
and expressing the polarization state of the light field to be measured through the target Stokes parameter.
2. The method according to claim 1, wherein the calculating a target mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data includes:
acquiring first Mueller matrices of a plurality of groups of polarization state generators and second Mueller matrices of a polarization state analyzer and target light intensities corresponding to each group of the first Mueller matrices and the second Mueller matrices according to the incident polarization state data and the emergent polarization state data;
calculating according to the first Mueller matrix and the second Mueller matrix to obtain an instrument matrix corresponding to the light field to be measured, and converting to obtain a generalized inverse matrix of the instrument matrix;
converting the obtained multiple groups of target light intensity into light intensity vectors;
and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
3. The method according to claim 1, wherein the obtaining of the target stokes parameter corresponding to the maximum probability value in the target probability density distribution function of the polarization state comprises:
taking logarithm of the target polarization state probability density distribution function to obtain a three-dimensional function on a Ponga spherical coordinate system;
selecting a plurality of testing Stokes parameters according to preset precision and inputting the testing Stokes parameters into a three-dimensional function on the Bonga spherical coordinate system to obtain a testing probability value corresponding to each testing Stokes parameter;
and taking the test Stokes parameter corresponding to the obtained maximum test rate value as the target Stokes parameter.
4. An optical field polarization state measurement apparatus, the apparatus comprising:
the polarization state generator is used for generating incident polarization state data of the calibration incident light;
the polarization state analyzer is used for generating emergent polarization state data of the calibration emergent light; the calibration incident light is output through a phase modulator to be tested to obtain calibration incident light;
the processor is used for calculating to obtain a target Mueller matrix of the phase modulator to be measured according to the incident polarization state data and the emergent polarization state data; inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix; the abscissa and the ordinate of the target polarization state probability density distribution function are respectively a stokes parameter and a corresponding probability value; the target polarization state probability density distribution function corresponds to a light field to be detected; obtaining a target Stokes parameter corresponding to the maximum probability value in the target polarization state probability density distribution function; expressing the polarization state of the light field to be detected through the target Stokes parameters;
inputting the target Mueller matrix into a pre-constructed polarization state probability density distribution model, and outputting a target polarization state probability density distribution function corresponding to the target Mueller matrix, wherein the method specifically comprises the following steps:
acquiring a point spread function light intensity distribution image corresponding to the light field to be detected; extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, the plurality of spatial position information and photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model; and converting the polarization state probability density distribution model after the target Mueller matrix, the plurality of spatial position information and the photon number information corresponding to each spatial position information are input, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
5. The apparatus of claim 4, wherein the polarization state generator and the polarization state analyzer are further configured to generate a first Mueller matrix of a multi-group polarization state generator and a second Mueller matrix of a polarization state analyzer,
the device further comprises a detector, a first light source and a second light source, wherein the detector is used for acquiring target light intensity corresponding to each group of the first Mueller matrix and the second Mueller matrix;
the processor is further configured to calculate an instrument matrix corresponding to the light field to be measured according to the first mueller matrix and the second mueller matrix, and convert the instrument matrix into a generalized inverse matrix of the instrument matrix; converting the obtained multiple groups of target light intensity into light intensity vectors; and multiplying the generalized inverse matrix point by the light intensity vector to obtain a target Mueller matrix of the phase modulator to be measured.
6. The apparatus of claim 5, wherein the polarization analyzer is configured to act as a left-handed analyzer for modulating the light field to be measured;
the detector is also used for acquiring a point spread function light intensity distribution image corresponding to the light field to be detected;
the processor is further used for extracting a plurality of pieces of spatial position information in the point spread function light intensity distribution image and photon number information corresponding to each piece of spatial position information; inputting the target Mueller matrix, the plurality of spatial position information and photon number information corresponding to each spatial position information into a pre-constructed polarization state probability density distribution model; and performing conversion processing on the input polarization state probability density distribution model, and outputting to obtain a target polarization state probability density distribution function corresponding to the target Mueller matrix.
7. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor realizes the steps of the method of any one of claims 1 to 3 when executing the computer program.
8. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810749870.7A CN108956097B (en) | 2018-07-10 | 2018-07-10 | Light field polarization state measuring method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810749870.7A CN108956097B (en) | 2018-07-10 | 2018-07-10 | Light field polarization state measuring method and device, computer equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108956097A CN108956097A (en) | 2018-12-07 |
CN108956097B true CN108956097B (en) | 2020-07-31 |
Family
ID=64483725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810749870.7A Active CN108956097B (en) | 2018-07-10 | 2018-07-10 | Light field polarization state measuring method and device, computer equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108956097B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975998B (en) * | 2019-02-02 | 2023-03-31 | 中国科学院国家授时中心 | Electro-optical polarization modulation system and electro-optical polarization modulation method |
CN110514598B (en) * | 2019-08-20 | 2020-07-14 | 北京理工大学 | Spectral polarization detection system and detection method for frequency domain modulation |
CN110806266A (en) * | 2019-11-11 | 2020-02-18 | 北京理工大学 | Selection method of polarization state analyzer in polarization state detection system |
CN111595330B (en) * | 2020-05-29 | 2021-09-14 | 北京航空航天大学 | Night polarization course calculation method based on probability density function estimation |
CN112035794B (en) * | 2020-07-13 | 2024-02-23 | 西安电子科技大学 | Polarization identification method and device based on Mueller matrix and electronic equipment |
CN112134618B (en) * | 2020-09-23 | 2022-04-15 | 京东方科技集团股份有限公司 | Endoscope-based polarized light decoding communication method |
CN112146760B (en) * | 2020-09-28 | 2021-08-31 | 浙江大学 | Polarization uniformity measuring method, application, device, electronic equipment and storage medium of completely non-uniform polarized light |
CN113984346B (en) * | 2021-10-15 | 2023-09-12 | 西安工业大学 | Device and method for characterizing optical transmission matrix of polarization imaging system |
CN114152578B (en) * | 2021-11-26 | 2023-10-20 | 北京理工大学 | Spatial modulation polarization detection method based on vortex wave plate |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201187993Y (en) * | 2008-01-29 | 2009-01-28 | 北京理工大学 | Device for large distance light parallel regulation |
CN104457995A (en) * | 2014-12-15 | 2015-03-25 | 清华大学深圳研究生院 | Fast polarization detector and detecting method |
CN104535192A (en) * | 2015-01-09 | 2015-04-22 | 北京理工大学 | Stokes polarization meter error calibration and compensation method based on wave plate rotating method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7701561B2 (en) * | 2006-02-13 | 2010-04-20 | Boston Applied Technologies Incorporated | Polarization imaging apparatus |
-
2018
- 2018-07-10 CN CN201810749870.7A patent/CN108956097B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201187993Y (en) * | 2008-01-29 | 2009-01-28 | 北京理工大学 | Device for large distance light parallel regulation |
CN104457995A (en) * | 2014-12-15 | 2015-03-25 | 清华大学深圳研究生院 | Fast polarization detector and detecting method |
CN104535192A (en) * | 2015-01-09 | 2015-04-22 | 北京理工大学 | Stokes polarization meter error calibration and compensation method based on wave plate rotating method |
Non-Patent Citations (1)
Title |
---|
Mueller 矩阵成像偏振仪的误差标定和补偿研究;李建慧等;《激光与光电子学进展》;20161231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108956097A (en) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108956097B (en) | Light field polarization state measuring method and device, computer equipment and storage medium | |
Altepeter et al. | Photonic state tomography | |
CN101238350B (en) | Method and apparatus for determining liquid crystal cell parameters from full mueller matrix measurements | |
Ghosh et al. | Estimating specular roughness and anisotropy from second order spherical gradient illumination | |
CN110806266A (en) | Selection method of polarization state analyzer in polarization state detection system | |
CN111443045B (en) | Spectral imaging analysis system and spectral imaging analysis method | |
CN105675455A (en) | Method and device for reducing random system noise in particle size analyzer | |
CN111562223A (en) | Polarizing imaging device and method based on micro-polarizer array | |
US20230304921A1 (en) | Infrared measurement method and apparatus, computer device and storage medium | |
WO2022179386A1 (en) | Distributed quantum imaging method, apparatus and system, and computer-readable storage medium | |
CN113538381A (en) | Method and system for rapidly detecting Mueller matrix of sample in weak light field | |
Duivenvoorden et al. | Full-sky beam convolution for cosmic microwave background applications | |
del Hoyo et al. | High precision calibration method for a four-axis Mueller matrix polarimeter | |
CN114136894B (en) | Error calibration method and device of polarization detection system based on vortex wave plate | |
Ning et al. | Single-shot measurement of polarization state at low light field using Mueller-mapping star test polarimetry | |
JP2017535779A (en) | System and method for sensing a remote object | |
CN116802465A (en) | Polarization measurement method | |
CN112035794A (en) | Polarization identification method and device based on Mueller matrix and electronic equipment | |
CN114152578B (en) | Spatial modulation polarization detection method based on vortex wave plate | |
CN105973467B (en) | Acousto-optic signal processing system based on compressed sensing | |
Chipman | Polarization diversity active imaging | |
CN112880987B (en) | Polarization performance detection method and system of optical element | |
Staes et al. | Optimized Stokes imaging for highly resolved optical speckle fields, Part I: optimized experimental setup | |
CN113988304A (en) | Magnetic flux crosstalk calibration method and device, computer equipment and storage medium | |
TW201804140A (en) | Stress analysis method based on temporal phase unwrapping |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |