CN110470399A - Full Stokes self-calibration measuring method based on polarization camera intensity collection twice - Google Patents
Full Stokes self-calibration measuring method based on polarization camera intensity collection twice Download PDFInfo
- Publication number
- CN110470399A CN110470399A CN201910745979.8A CN201910745979A CN110470399A CN 110470399 A CN110470399 A CN 110470399A CN 201910745979 A CN201910745979 A CN 201910745979A CN 110470399 A CN110470399 A CN 110470399A
- Authority
- CN
- China
- Prior art keywords
- wave plate
- stokes vector
- degree
- twice
- full
- 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.)
- Pending
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 8
- 238000003491 array Methods 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000000711 polarimetry Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Studio Devices (AREA)
Abstract
The invention discloses a kind of full Stokes vector self-calibration measuring methods based on polarization camera intensity collection twice, step 1, according to luminous intensity measurement variance, the corresponding calculation matrix of analyzer, it derives based on the functional relation of point focal plane polarization camera full the Stokes vector estimate variance and PSA calculation matrix of luminous intensity measurement twice, obtains the estimation variances sigma of 4 component elementsm 2；Step 2 makes optimal rotation angle twice that should meet 45 degree of difference；Step 3, according to optimal solution, intensity collection twice is obtained to 0 degree and 45 degree by rotating wave plate, and the corresponding calculation matrix W of analyzer PSA is calculated according to optimization8×4, full Stokes vector is calculated by light intensity value, step 4, according to the self-calibration of the relational implementation retardation of wave plate between measurement light intensity value, each component of full Stokes vector and rotatable retardation of wave plate.The method of the present invention is easy to operate, and estimate variance is small, greatly reduces the cost of polarization measurement system.
Description
Technical field
The invention belongs to polarimetry technical fields, more particularly to a kind of full Stokes self-calibration measuring method.
Background technique
Basic physical message one of of the polarization information as light wave, it is possible to provide other light-wave informations cannot provide tested
Object information, therefore the measurement of polarization information has very extensive application in many fields.Stokes vector describes light wave
Polarization state, contains most basic polarization information, thus the measurement of Stokes vector also becomes the main side in polarimetry field
One of to.
Point focal plane polarization camera is extensive because it is capable of real-time acquisition the polarization information obtained under four different polarization states
Applied to polarimetry and polarization imaging field.By placing rotatable wave plate before point focal plane polarizes camera, it can be achieved that complete
The measurement of Stokes vector.The selection of rotation angle determines the size of full Stokes vector estimate variance.On the other hand, but by
The influence of environment is subjected to when the phase-delay quantity of wave plate, it is therefore desirable to demarcate in advance, this brings to practical operation and application
Inconvenience.
Summary of the invention
It is an object of the invention to overcome deficiencies of the prior art, one kind is provided and is based on polarization camera light twice
The full Stokes self-calibration measuring method acquired by force is realized and is surveyed based on a point full Stokes vector estimation for focal plane polarization camera
The phase-delay quantity self-calibration of amount and rotating wave plate.
A kind of full Stokes vector self-calibration measuring method based on polarization camera intensity collection twice of the invention, the party
Method specific implementation the following steps are included:
Step 1, according to luminous intensity measurement variance, the corresponding calculation matrix of analyzer, derive based on a point focal plane polarization phase
The functional relation of machine full the Stokes vector estimate variance and PSA calculation matrix of luminous intensity measurement twice, obtains 4 component elements
Estimate variances sigmam 2, expression formula are as follows:
Wherein, σ2For full Stokes vector estimate variance, m is that the component elements of full Stokes vector are numbered；
Step 2 makes optimal rotation angle twice that should meet 45 degree of difference；
Step 3, according to optimal solution, intensity collection twice is obtained to 0 degree and 45 degree by rotating wave plate, and count according to optimization
Calculation obtains the corresponding calculation matrix W of analyzer PSA8×4, expression formula are as follows:
Wherein, θ1=0 °, θ2=45 ° respectively indicate the angle of wave plate rotation twice；The phase-delay quantity of δ expression wave plate；
Full Stokes vector, expression formula are calculated by light intensity value are as follows:
Wherein, S=[s0,s1,s2,s3] indicate full Stokes vector to be measured, I0~I3When indicating that wave plate is rotated to 0 degree,
The corresponding light intensity value obtained of four micro- polarizing films, I4~I7When indicating that wave plate is rotated to 45 degree, four micro- polarizing films are corresponding to be obtained
Light intensity value；
Step 4, according to measurement light intensity value, the pass between each component of full Stokes vector and rotatable retardation of wave plate
It is the analytical expression obtained from calibration algorithm are as follows:
Wherein, δ indicates phase-delay quantity；
To realize the self-calibration of retardation of wave plate.
Beneficial effects of the present invention and advantage are:
1, by rotating wave plate to 0 degree and 45 degree acquisition, intensity collection can be realized full Stokes vector and estimates the present invention twice
The self-calibration of meter and retardation of wave plate, it is easy to operate quick, it is easy to accomplish；
2, the gaussian additive noise of the invention not only for environment, and it is same for the system for mixing other types noise
With theoretic optimization feasibility；
3, the method for the present invention no longer needs to carry out the rotating wave plate in system individually additional staking-out work, operation letter
Single, estimate variance is small, greatly reduces the cost of polarization measurement system.
Detailed description of the invention
Fig. 1 is of the invention based on the full Stokes vector measurement method overall flow for dividing focal plane to polarize camera signal
Figure；
Fig. 2 is that the realization based on the full Stokes vector measurement method for dividing focal plane to polarize camera of the embodiment of the present invention fills
Set schematic diagram；
Fig. 3 is a point micro- polarizing film array distribution figure for focal plane polarization camera.
Appended drawing reference:
1, rotatable wave plate, 2, micro- polarization chip arrays, 3, light intensity detector (CCD), 4, divide focal plane polarize camera, 5,0
Spend the micro- polarizing film in direction, 6,45 degree of micro- polarizing films in direction, 7,135 degree of micro- polarizing films in direction, 8,90 degree of micro- polarizing films in direction.
Specific embodiment
Technical solution of the present invention is described in detail with reference to the accompanying drawings and examples.
Theoretical foundation of the invention:
Divide each unit pixel before dividing focal plane polarization camera to the micro- polarization etched as 2x2 by micro-nano technology technology
Chip arrays, angleRespectively 0 degree, 45 degree, 135 degree, 90 degree, it can be achieved that under four kinds of different polarization states light intensity it is real-time same
Pacing amount, so that one-shot measurement realizes Stokes vector (s0,s1,s2) measurement estimation.By this divide focal plane polarize camera
One rotatable wave plate of preceding placement and rotate can realize full Stokes vector (s twice0,s1,s2,s3) measurement estimation.It obtains
Light intensity matrix I, expression formula are as follows:
I=W8×4·S (1)
Wherein, W8×4Expression measures corresponding analyzer (PSA) calculation matrix, S=[s twice0,s1,s2,s3] indicate to be measured
Full Stokes vector.
By corresponding four polarization angles of micro- polarization chip arraysWith rotatable wave plate (phase-delay quantity δ)
Rotation angle, θ twicei, i ∈ [1,2] calculate light intensity matrix I in each element Iij, expression formula are as follows:
Wherein, i is rotation angle number, and j is polarization angle number；
Light intensity in anti-solution formula (1), obtains Stokes vector estimator to be measuredExpression formula are as follows:
Wherein,Indicate calculation matrix W8×4Pseudo inverse matrix.
Assuming that ambient noise is Gaussian, i.e. it is < I that each component of detection light intensity, which obeys mean value,ij>, variance σ2
Gaussian Profile, the measurement of each component is mutually independent of each other, then it is σ that the covariance matrix of light intensity, which is main diagonal element,2It is diagonal
Battle array.And then obtain the estimation variance of 4 component elements of Stokes vector, it may be assumed that
By optimizing calculation matrix, that is, wave plate rotation angle, each component elements of reduction Stokes vector are corresponding to be estimated
Variance is calculated, to realize the raising of the Estimation and Measurement precision of Stokes vector.
Use the variance of each element and the standard variance V as estimation precision:
Then optimization problem is to find optimal calculation matrix W8×4, so that variance V is minimum.It is theoretically available above-mentioned to ask
The minimum value of topic is 5.5 σ2, the wave plate of corresponding optimization rotate twice angle satisfaction | θ1-θ2|=45 °.
For the sake of simplicity, θ can be chosen1=0 °, θ2=45 ° are used as optimal rotation angle.In such situation, according to each light intensity meter
Calculate Stokes vector to be measured, expression formula are as follows:
Wherein, I0~I3When indicating that wave plate is rotated to 0 degree, the corresponding light intensity value obtained of four micro- polarizing films, I4~I7It indicates
When wave plate is rotated to 45 degree, the corresponding light intensity value obtained of four micro- polarizing films.
Meanwhile according to the pass between measurement light intensity value, each component of full Stokes vector and the mutually rotatable position retardation of wave plate
System obtains the analytical expression that the related light intensity of retardation of wave plate indicates:
To realize the synchronization self-calibration of retardation of wave plate.
As shown in Figure 1, for a kind of full Stokes vector self-calibration based on polarization camera intensity collection twice of the invention
Measurement method overall flow schematic diagram.This method is directed to the measuring system under gaussian additive noise environment, by that will be placed on point
The angle of rotatable wave plate in front of the polarization camera of focal plane is rotated respectively to 0 degree and 45 degree, realizes full Stokes arrow for the first time
The minimum variance estimate of amount parsing self-calibration synchronous with rotating wave plate phase-delay quantity.In practical applications, it is no longer necessary to right
Rotatable wave plate in system carries out individual phase-delay quantity staking-out work, and simple and quick, estimate variance is small, greatly drops
The low cost of measuring system.Specific implementation process the following steps are included:
Step 1, according to luminous intensity measurement variance, the corresponding calculation matrix of analyzer (PSA), derive based on point focal plane it is inclined
The full Stokes vector estimate variance of luminous intensity measurement and the (rotation depending on rotatable wave plate of PSA calculation matrix twice of vibration camera
Angle) functional relation；
Step 2 obtains optimal wave plate rotation angle by minimizing estimation method, verifies optimal rotation angle twice
45 degree of difference should be met；
Step 3, according to optimal solution, intensity collection twice is obtained to 0 degree and 45 degree by rotating wave plate, and count according to optimization
Calculation obtains the corresponding calculation matrix of PSA, and being calculated by light intensity value indicates full Stokes vector；
Step 4, according to measurement light intensity value, the relationship between each component of Stokes vector and rotatable retardation of wave plate
The self-calibration for realizing retardation of wave plate, obtained from the analytical expression of calibration algorithm.
The specific embodiment of the invention is described as follows:
In actual measurement, intensity collection totally 8 kinds of different polarization states twice are obtained to 0 degree and 45 degree by rotating wave plate
Under light intensity value, by above-mentioned formula (6) it is synchronous with (7) realize full Stokes vector estimation and retardation of wave plate from
Calibration.
As shown in Fig. 2, being the embodiment of the present invention based on the full Stokes vector measurement method for dividing focal plane to polarize camera
Realization device schematic diagram, selected light intensity detection device is a point focal plane polarization camera 4, by micro- 2 He of polarization chip arrays
3 two parts of ccd detector are constituted.Light with full Stokes vector to be measured enters a point focal plane after rotatable wave plate 1 and polarizes
Camera 4.The rotatable degree of wave plate 1 to 0 and 45 degree are rotated respectively, and then realize that a point measurement twice for focal plane polarization camera 4 obtains
Light intensity measurement (I under eight kinds of difference PSA polarization states0-I7).The device and optimization after rotatable wave plate rotation angle and
PSA calculation matrix can synchronize the low variance evaluation and wave for realizing full Stokes vector compared with traditional measuring device and method
The self-calibration of piece phase-delay quantity, it is simple and convenient, greatly reduce the cost of measuring system.
Claims (1)
1. a kind of full Stokes vector self-calibration measuring method based on polarization camera intensity collection twice, which is characterized in that should
Method specific implementation the following steps are included:
Step 1, according to luminous intensity measurement variance, the corresponding calculation matrix of analyzer, derive based on be divided to focal plane polarize camera two
The functional relation of full the Stokes vector estimate variance and PSA calculation matrix of secondary luminous intensity measurement, obtains the estimation of 4 component elements
Variances sigmam 2, expression formula are as follows:
Wherein, σ2For full Stokes vector estimate variance, m is that the component elements of full Stokes vector are numbered；
Step 2 makes optimal rotation angle twice that should meet 45 degree of difference；
Step 3, according to optimal solution, intensity collection twice is obtained to 0 degree and 45 degree by rotating wave plate, and calculate according to optimization
To the corresponding calculation matrix W of analyzer PSA8×4, expression formula are as follows:
Wherein, θ1=0 °, θ2=45 ° respectively indicate the angle of wave plate rotation twice；The phase-delay quantity of δ expression wave plate；
Full Stokes vector, expression formula are calculated by light intensity value are as follows:
Wherein, S=[s0,s1,s2,s3] indicate full Stokes vector to be measured, I0~I3When indicating that wave plate is rotated to 0 degree, four
The corresponding light intensity value obtained of micro- polarizing film, I4~I7When indicating that wave plate is rotated to 45 degree, the corresponding light intensity obtained of four micro- polarizing films
Value；
Step 4 is obtained according to the relationship between measurement light intensity value, each component of full Stokes vector and rotatable retardation of wave plate
Analytical expression derived from calibration algorithm are as follows:
Wherein, δ indicates phase-delay quantity；
To realize the self-calibration of retardation of wave plate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745979.8A CN110470399A (en) | 2019-08-13 | 2019-08-13 | Full Stokes self-calibration measuring method based on polarization camera intensity collection twice |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745979.8A CN110470399A (en) | 2019-08-13 | 2019-08-13 | Full Stokes self-calibration measuring method based on polarization camera intensity collection twice |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110470399A true CN110470399A (en) | 2019-11-19 |
Family
ID=68511803
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910745979.8A Pending CN110470399A (en) | 2019-08-13 | 2019-08-13 | Full Stokes self-calibration measuring method based on polarization camera intensity collection twice |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110470399A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114324185A (en) * | 2022-01-04 | 2022-04-12 | 浙江大学 | Underwater polarization detection device based on Stokes vector |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104034426A (en) * | 2014-06-11 | 2014-09-10 | 中国科学技术大学 | Real-time polarization state and phase measurement method based on pixel polarizing film array |
| CN108731810A (en) * | 2018-06-07 | 2018-11-02 | 合肥工业大学 | Divide focal plane polarization imaging measurement method and device |
| CN108918425A (en) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | A kind of Muller matrix measuring system and method |
-
2019
- 2019-08-13 CN CN201910745979.8A patent/CN110470399A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104034426A (en) * | 2014-06-11 | 2014-09-10 | 中国科学技术大学 | Real-time polarization state and phase measurement method based on pixel polarizing film array |
| CN108918425A (en) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | A kind of Muller matrix measuring system and method |
| CN108731810A (en) * | 2018-06-07 | 2018-11-02 | 合肥工业大学 | Divide focal plane polarization imaging measurement method and device |
Non-Patent Citations (2)
| Title |
|---|
| S HUHEI S HIBATA等: "Robust full Stokes imaging polarimeter with dynamic calibration", 《OPTICS LETTERS》 * |
| S TÉPHANE R OUSSEL等: "On the optimal ways to perform full Stokes measurements with a linear division-of-focal-plane polarimetric imager and a retarder", 《OPTICS LETTERS》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114324185A (en) * | 2022-01-04 | 2022-04-12 | 浙江大学 | Underwater polarization detection device based on Stokes vector |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hurley-Walker et al. | The Murchison Widefield Array Commissioning Survey: A low-frequency catalogue of 14,110 compact radio sources over 6,100 square degrees | |
| EP2961231A1 (en) | Angle determining system and method | |
| CN103363949B (en) | Mixed measurement analysis method for satellite antenna | |
| CN110132420B (en) | Polarization measuring device, polarization measuring method, and optical alignment method | |
| CN104897060A (en) | large Large field of view global measurement method using coordinates tracking control board | |
| CN101419064B (en) | Relative rotation measuring and rotating positioning system and method for small gapping place | |
| CN110470399A (en) | Full Stokes self-calibration measuring method based on polarization camera intensity collection twice | |
| CN103439001A (en) | Method and device for measuring and evaluating inhomogeneous vector polarized light | |
| CN104236857A (en) | Liquid crystal optical phase shift distribution detection system and method based on quarter-wave plate method | |
| CN106872936A (en) | Near field sources L-type acoustic vector-sensor array row ambiguity solution Multiple Parameter Estimation Methods | |
| CN110487229A (en) | Circumference equal-dividing hole minimum position degree calculation method | |
| CN104980237B (en) | Channel phases postpone automatic calibrating method and signal processor, RSU | |
| CN112240957B (en) | Method for correcting amplitude-phase characteristics of antenna in satellite navigation interference direction finding | |
| CN109490867A (en) | Waterborne target polarization remote sensing detectivity evaluation method | |
| CN110332994A (en) | Full Stokes vector measurement method based on polarization camera and optimization wave plate | |
| CN105872403B (en) | The dynamic rage extension method of Mueller matrix imaging | |
| CN104833425B (en) | The combined method of polaroid angle in a kind of polarized imaging system | |
| CN103115683B (en) | Method of measuring monochromatic light polarization state | |
| CN109459138A (en) | Mueller type ellipsometer ellipsometric parameter measurement method and device based on four luminous intensity measurements | |
| CN206959992U (en) | A kind of laser polarization state measuring instrument | |
| CN103630336A (en) | Dynamic interference measuring method based on random fast axis azimuth delay array | |
| CN109781317A (en) | Optical glass stress detection system and detection method | |
| CN201257719Y (en) | Viscose glue machine of X ray orientation device | |
| CN104931234A (en) | Device for precisely calibrating inclined angle between wave plate and polaroid, and calibration method | |
| CN107271041A (en) | A kind of laser polarization state measuring instrument |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191119 |
|
| WD01 | Invention patent application deemed withdrawn after publication |