CN115272041A - 起偏器与检偏器配置优化方法及起偏检偏系统 - Google Patents
起偏器与检偏器配置优化方法及起偏检偏系统 Download PDFInfo
- Publication number
- CN115272041A CN115272041A CN202210840465.2A CN202210840465A CN115272041A CN 115272041 A CN115272041 A CN 115272041A CN 202210840465 A CN202210840465 A CN 202210840465A CN 115272041 A CN115272041 A CN 115272041A
- Authority
- CN
- China
- Prior art keywords
- polarizer
- analyzer
- measurement
- configuration
- wave plate
- 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 167
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 160
- 238000005259 measurement Methods 0.000 claims description 200
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 16
- 230000002068 genetic effect Effects 0.000 claims description 12
- 230000008901 benefit Effects 0.000 claims description 10
- 239000013598 vector Substances 0.000 description 22
- 230000008859 change Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 241000962146 Alsophila tricolor Species 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
- G06T1/0007—Image acquisition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/60—Rotation of whole images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10141—Special mode during image acquisition
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Genetics & Genomics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种起偏器与检偏器配置优化方法及起偏检偏系统,所述起偏器与检偏器配置优化方法包括如下步骤:调整起偏器的仪器矩阵W和检偏器的仪器矩阵A,使起偏器和检偏器仪器矩阵的等权重方差EWV最小,以实现针对高斯噪声的优化;且起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,使泊松噪声所引起的估计方差独立于样本,且估计方差达到最小值。本发明能够最大程度地抑制噪声,并使噪声的规律与样本独立,无论测量何种样本,噪声的分布规律都是相同的。
Description
技术领域
本发明涉及偏振光学成像技术领域,特别是涉及一种起偏器与检偏器配置优化方法及起偏检偏系统。
背景技术
偏振成像技术具有非侵入、无损伤、在体原位、数据量大等优点,广泛应用于生物医学、海洋科学、大气遥感等领域。根据最终获得偏振信息形式的不同,偏振成像可以分为斯托克斯向量测量和缪勒矩阵测量。斯托克斯向量主要用于描述光的偏振特性,而缪勒矩阵则是表征样本的偏振特征,并可通过缪勒矩阵进一步提取出介质的微观结构信息,这也是偏振测量相比于传统的光学测量最大的优势所在。
在进行缪勒矩阵测量时,需要对入射光进行偏振调制,同时对出射光进行偏振属性的检测,从而获悉样本对光的偏振特性产生了何种影响与改变。通过多次调制与检测,就可以获得样本完备的偏振特性,亦即样本的缪勒矩阵。在这个过程中,调制入射光偏振态的器件被称为起偏器(PSG),而对光的偏振属性进行检测的器件称为检偏器(PSA)。起偏器与检偏器的结构类似,其核心都是对光的偏振态进行调制,起偏器是将自然光调制为特定偏振态的偏振光;检偏器是将入射偏振光反向调制,从而得到该入射偏振光在某一特定偏振态分量上的大小,并通过计算得到光的偏振属性。
斯托克斯矢量是一种描述光的偏振属性的一种方法,其表示为S=[S0 S1 S2 S3]T,是一个四维向量,S0表示光的强度;S1=I0-I90,是光在0度偏振方向的光强分量和90度偏振方向的光强分量之差;S2=I45-I135,是光在45度偏振方向的光强分量和135度偏振方向的光强分量之差;S4=IR-IL,是光在右旋偏振方向的光强分量和左旋偏振方向的光强分量之差。通常情况下,我们只关心光的偏振属性,因此对S0进行归一化处理,即S0=1,其余三个分量也按相同比例进行归一化处理。
邦加球是用来图形化描述偏振态斯托克斯矢量的一个单位球,任何一种偏振态都可以映射到邦加球上的一点,将斯托克斯矢量的S1、S2、S3分别作为x、y、z坐标绘制在笛卡尔坐标系中,即为该斯托克斯矢量的邦加球表示。
而缪勒矩阵是一个变换矩阵,体现一束光的斯托克斯矢量在散射前后发生的变化,及
Sout=M×Sin
式中,Sout是出射光的斯托克斯矢量,Sin是入射光的斯托克斯矢量,M是缪勒矩阵,是一个4×4的矩阵。由于CCD(电荷耦合器件)无法接受到偏振信息,仅能接收到光强信息。因此在实际测量中,需要进行至少四次互相独立的起偏及检偏。起偏即为用一束已知偏振态的偏振光入射,检偏即为通过测量计算得到出射光的偏振态Sout。由于斯托克斯矢量有四个分量,检偏时需要至少测量四个投影分量的光强大小,才能得到一束光的斯托克斯矢量。这样通过测量多种不同偏振态的偏振光照明后的出射偏振态,即可计算出斯托克斯矢量的变化矩阵,即缪勒矩阵。
目前偏振调制主要通过一个偏振片,再加一个或一系列相位延迟器件实现的,借助相位延迟器件不同的快轴角以及相位延迟大小之间的彼此配合得到多个不同偏振态。因此需要首先对快轴角及相位延迟对偏振态的影响进行说明,一个相位延迟器件的缪勒矩阵Mδ,θ为:
其中δ为线性相位延迟,θ为快轴角方向。
假设此时入射光的偏振态为Sin=[1 S1 S2 S3]T,其中T表示矩阵的转置,即任意的完全偏振光时,出射光的偏振态可用公式Sout=Sin×Mδ,θ表示。这样在同一确定的入射光(通常为线偏振片产生的0度线偏振光)从一相位延迟器出射时,我们可以得到一个确定的偏振态,通过改变相位延迟器件的快轴角θ或相位延迟δ,出射偏振态也会发生相应的改变。这样通过多次改变相位延迟器件,就可以得到多个不同且已知的偏振态,组成起偏器的仪器矩阵W或检偏器的仪器矩阵A。
相位延迟器件又可分为固定相位延迟(通过旋转改变快轴角的波片)以及电控式起偏元件(相位延迟可变的)。前者的主要优势在于精度高,缺点则是有运动部件,需要旋转电机带动波片旋转才能实现多种不同偏振态的调制,因而调制速度受电机速度影响,往往比较耗时。而后者主要有液晶相位延迟器(LCVR)、电光相位延迟器(EO)及光弹调制器(PEM)等通过改变施加在其两端的电压从而改变相位延迟的器件。其主要优势在于响应速度快,可达到毫秒量级。但因环境温度或电压波动等会造成一定的系统误差,且成本较高。此外,由于电控式起偏元件需要线缆进行控制,往往也不适合配合电机转动,因此其快轴角在测量中是无法自由旋转的。
为了使缪勒矩阵的测量更为精确,抑制噪声表现,可以对起偏器及检偏器进行设计与优化。
发明内容
为了弥补上述背景技术的不足,本发明提出一种起偏器与检偏器配置优化方法及起偏检偏系统,以解决缪勒矩阵测量中精度有限、无法很好地抑制噪声的问题。
本发明的技术问题通过以下的技术方案予以解决:
本发明公开了一种起偏器与检偏器配置优化方法,包括如下步骤:
调整起偏器的仪器矩阵W和检偏器的仪器矩阵A,使起偏器和检偏器仪器矩阵的等权重方差EWV最小,以实现针对高斯噪声的优化;且
起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,使泊松噪声所引起的估计方差独立于样本,且估计方差达到最小值。
在一些实施例中,还可能包括以下技术方案:
配置测量系统的起偏器以及检偏器,使其偏振态两两正交;随后搜索起偏器和检偏器仪器矩阵的最小等权重方差EWV;当起偏器和检偏器的仪器矩阵同时满足行和为0及最小EWV条件时,得到最优配置。
当起偏器和检偏器分别为一个旋转偏振片及一个旋转1/4波片时,使起偏器和检偏器偏振态两两正交的配置条件满足以下关系:
式中,θP为第一次测量配置的偏振片透光轴方向,θR为第一次测量配置的1/4波片快轴角方向,θP'为第二次测量配置的偏振片透光轴方向,θR'为第二次测量配置的1/4波片快轴角方向;所述方法通过偏振片调整起偏的主通光方向,再通过1/4波片实现特定偏振态的调制。
当起偏器和检偏器系统包括一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片时,使起偏器和检偏器偏振态两两正交的配置条件满足以下关系:
式中,θH为第一次测量配置的1/2波片快轴角方向,θR为第一次测量配置的1/4波片快轴角方向,θH'为第二次测量配置的1/2波片快轴角方向,θR'为第二次测量配置的1/4波片快轴角方向;
所述方法的线偏振光方向的调制从单偏振片调制变成了固定偏振片加一个1/2波片,能够避免入射光不是理想的自然光时,转动偏振片导致的不同方向的线偏振光光强不一致的情况。
当起偏器和检偏器系统包括一个固定偏振片,两个全波延迟相位调制器件时,使起偏器和检偏器偏振态两两正交的配置条件满足以下关系:
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟。所述方法首先通过第一个全波延迟相位调制器件对偏振态在邦加球上的S1OS3平面的圆上调制,随后借助第二个全波延迟相位调制器件得到目标偏振态。
当起偏器和检偏器分别为一个旋转偏振片及一个旋转1/4波片时,使起偏器和检偏器仪器最优的四点测量配置满足以下关系:
式中,θP为第一次测量配置的偏振片透光轴方向,θR为第一次测量配置的1/4波片快轴角方向,θP'为第二次测量配置的偏振片透光轴方向,θR'为第二次测量配置的1/4波片快轴角方向,θP”为第三次测量配置的偏振片透光轴方向,θR”为第三次测量配置的1/4波片快轴角方向,θP”'为第四次测量配置的偏振片透光轴方向,θR”'为第四次测量配置的1/4波片快轴角方向。
当起偏器和检偏器系统包括一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片时,使起偏器和检偏器最优的四点测量配置满足以下关系:
式中,θH为第一次测量配置的1/2波片快轴角方向,θR为第一次测量配置的1/4波片快轴角方向,θH'为第二次测量配置的1/2波片快轴角方向,θR'为第二次测量配置的1/4波片快轴角方向,θH”为第三次测量配置的1/2波片快轴角方向,θR”为第三次测量配置的1/4波片快轴角方向,θH”'为第四次测量配置的1/2波片快轴角方向,θR”'为第四次测量配置的1/4波片快轴角方向。
当起偏器和检偏器系统包括一个固定偏振片,两个全波延迟相位调制器件时,使起偏器和检偏器最优的四点测量配置满足以下关系:
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″是第三次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″是第三次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″′是第四次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″′是第四次测量配置的第二个全波延迟相位调制器件的相位延迟。
当起偏器和检偏器系统包括一个固定偏振片,两个半波延迟相位调制器件时,使起偏器和检偏器最优的四点测量配置满足以下关系:
式中,δH1、δH2分别为第一次测量配置的两个半波延迟相位调制器件的线性相位延迟,θH1、θH2则分别为其快轴角,δH1′、δH2′分别为第二次测量配置的两个半波延迟相位调制器件的线性相位延迟,δH1″、δH2″分别为第三次测量配置的两个半波延迟相位调制器件的线性相位延迟,δH1″′、δH2″′分别为第四次测量配置的两个半波延迟相位调制器件的线性相位延迟。
进一步地,用实际配置表示的偏振态组成一个仪器矩阵,并利用遗传算法或类似的寻优算法对该仪器矩阵进行最小EWV优化,即,对每一组设置一个变量,并根据该变量得到一组中剩余的能够由该变量明确表示的偏振态;之后,这些偏振态将组成一个有数个未知变量的仪器矩阵;通过遗传算法对仪器矩阵进行最小EWV寻优,即可计算出当这数个未知变量为何值时,仪器矩阵的EWV能达到最小值。
进一步地,起偏器和检偏器均满足全偏振态的调制。
本发明还公开了一种起偏检偏系统,包括一个旋转偏振片及一个旋转1/4波片,配置如下:
,该配置满足起偏器和检偏器的仪器矩阵W与A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优,以优化缪勒测量系统抗高斯-泊松混合噪声的性能。
本发明还公开了一种起偏检偏系统,包括一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片,配置如下:
,该配置满足起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优,以优化缪勒测量系统抗高斯-泊松混合噪声的性能。
本发明还公开了一种起偏检偏系统,包括一个固定偏振片,两个全波延迟相位调制器件,配置如下:
,该配置满足起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优,以优化缪勒测量系统抗高斯-泊松混合噪声的性能。
本发明还公开了一种起偏检偏系统,包括一个固定偏振片及两个半波延迟相位调制器件,配置如下:
,该配置满足起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优,以优化缪勒测量系统抗高斯-泊松混合噪声的性能。
本发明与现有技术对比的有益效果包括:本发明提供的起偏器与检偏器配置优化方法,通过调整起偏器的仪器矩阵W和检偏器的仪器矩阵A,使起偏器和检偏器仪器矩阵的等权重方差EWV最小,以实现针对高斯噪声的优化,且起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,使泊松噪声所引起的估计方差独立于样本,且估计方差达到最小值,从而实现最大程度地抑制噪声,并使噪声的规律与样本独立,无论测量何种样本,噪声的分布规律都是相同的。
在一些实施例中,本发明将测量配置与高斯-泊松混合噪声的优化相结合,通过采用多种不同器件下的起偏器及检偏器的最优测量配置,使起偏器和检偏器仪器矩阵的等权重方差EWV最小,进一步提升偏振测量,特别是缪勒矩阵测量中的性能精度,实现将缪勒矩阵测量中的误差传递影响降至最小,且使得泊松噪声表现与样本无关。同时,本方法具备普适性,可以对各种不同结构的测量系统进行优化。
附图说明
图1a至1b为本发明实施例中的旋转式相位延迟器件形成的一组正交偏振态图。
图2a至2b为本发明实施例中的相位可变延迟器件形成的一组正交偏振态图。
图3为本发明实施例中的一个旋转偏振片及一个旋转1/4波片的仪器示意图。
图4为本发明实施例中的一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片的仪器示意图。
图5为本发明实施例中的一个固定偏振片,两个全波延迟相位调制器件的仪器示意图。
图6为本发明实施例中的一个固定偏振片及两个半波延迟相位调制器件的仪器示意图。
图7a至7b为本发明实施例中的一个固定偏振片及两个半波延迟相位调制器件的仪器下可调制出的所有偏振态及最优的四点测量框架图。
图8为本发明实施例的起偏器与检偏器配置优化方法流程图。
具体实施方式
下面对照附图并结合优选的实施方式对本发明作进一步说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
需要说明的是,本实施例中的左、右、上、下、顶、底等方位用语,仅是互为相对概念,或是以产品的正常使用状态为参考的,而不应该认为是具有限制性的。
经研究发现,合理选择起偏器产生的入射偏振态(即起偏器的仪器矩阵W)以及检偏器的几个偏振态投影通道(即检偏器的仪器矩阵A)可以大幅提高缪勒矩阵测量系统的性能。因此本发明实施例的优化即为选取更加合适的W和A。而仪器矩阵W和A在何种情况下达到最优由以下推导给出:
在缪勒矩阵测量系统中往往存在高斯-泊松混合噪声,评估噪声对测缪勒矩阵估计方差的影响需要借助协方差分析。在缪勒矩阵测量系统中,直接测量得到的光强与仪器矩阵和样本缪勒矩阵的关系可以表示为:
I=ATMW
其中I为CCD测得的光强图像,M为样本的缪勒矩阵,W和A分别表示起偏器(PSG)和检偏器(PSA)的仪器矩阵,T代表矩阵的转置形式。为方便计算,借助克罗内克积将光强矩阵和样本的缪勒矩阵展开为向量形式:
其中VI和VM分别为上述光强矩阵I和样本的缪勒矩阵M的列向量形式,-1代表矩阵的逆或伪逆。根据VM和VI的关系,其协方差矩阵之间的关系可以定义为:
其中QA=(ATA)-1QW=(WTW)-1(这是算式的一种简写形式),EWVW和EWVA分别代表起偏器和检偏器仪器矩阵的等权重方差。由上式可以看出,针对高斯噪声的优化,仅需起偏器和检偏器仪器矩阵的等权重方差(EWV)最小,即可达到最优状态。
归一化后的VM其第一项永远是1。显然,在上式中,第一项是独立于样本的,而第二项与样本的后15个阵元相关,会导致随着样本的改变,泊松噪声所引起的估计方差会随之改变。幸运的是,当起偏器和检偏器仪器矩阵的EWV是最优时,为恒定的常数。因此,当满足下式时,CPoission(C泊松)的第二项会被置零,进而使得泊松噪声对缪勒矩阵所造成的估计方差独立于样本:
上式的等价条件为:起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0。在满足此条件且起偏器和检偏器仪器矩阵的EWV是最优时,泊松噪声所引起的估计方差独立于样本,且估计方差达到了最小值。此时最优化的CPoission可以表示为:
综上,当起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优时,缪勒矩阵测量系统抗高斯-泊松混合噪声的性能最优,系统总体的估计方差可以表示为:
至此,我们已经了解仪器矩阵W与A在何种情况下最优。但此时的仪器矩阵与实际的测量系统之间没有直接的联系,只是得到了仪器矩阵在满足什么条件时最优。若要使实际的测量系统能够实现最优,即使测量系统的仪器矩阵满足上述条件,仍需要进一步地优化。由于这一步优化是为了将实际的测量与理论联系起来,下面我们以几个实际的测量系统作为范例,展现优化的方法与过程。
考虑到旋转式与相位可变式的延迟器有着不同的特点,旋转式器件由高精度旋转电机带动,整体测量的误差较小,但由于运动部件响应速度较慢,测量一般比较耗时。而电控式相位延迟可变延迟器则由于温度漂移等的影响,精度会稍差些,但响应速度在毫秒量级,可以实现非常快速的测量,但成本也较高。因此两者往往不会混用,我们分别对这两类相位延迟器件组成的系统进行优化。旋转式相位延迟器件形成的一组正交偏振态如图1a至1b所示,相位可变延迟器件形成的一组正交偏振态如图2a至2b所示。
本发明实施例通过采用多种不同器件下的起偏器及检偏器的最优配置,进一步提升偏振测量,特别是缪勒矩阵测量中的性能精度。实现将缪勒矩阵测量中的误差传递影响降至最小,且使得泊松噪声表现与样品无关的最优配置。
前文提到,当起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优时,缪勒矩阵测量系统抗高斯-泊松混合噪声的性能最优。因此,如果使测量系统的W与A满足该条件,自然能使测量系统最优,实现噪声传递最小且泊松噪声样本无关。但是测量系统的W与A是由测量系统的起偏器与检偏器所产生的偏振态决定的,而产生何种偏振态是由测量系统的配置(各个元件的快轴角、相位延迟等)决定的。因此,优化W与A的问题就转换成对测量系统配置的优化问题。
显然,当仪器矩阵的偏振态两两正交时,两个正交偏振态的斯托克斯矢量的后三个分量互为相反数,整个仪器矩阵的行和自然为零,也就是满足了最优仪器矩阵的第一个要求。而在邦加球中,两个关于圆心对称的偏振态就是一组正交的偏振态。而仪器矩阵本身是由多个偏振态组成的,仪器矩阵的N个偏振态可以在邦加球上形成一个多面体。如果该多面体任意一个顶点都有一个关于圆心对称的顶点,那么每一个偏振态也都有着一个对应的正交偏振态。此时问题的关键就变成了如何配置测量系统的起偏器以及检偏器,使其偏振态两两正交。而根据相位延迟器件的缪勒矩阵Mδ,θ:
可以看出虽然我们可以很方便地计算出任意偏振光经过相位延迟器件后出射光的偏振态,但在进行偏振态调制时,我们更加注重的是什么样的快轴角及相位延迟关系能满足将入射光调制为目标偏振态。而这种情况下,常规的正向思路并不能实现我们的要求,往往需要将快轴角及相位延迟作为输入量,通过不断调整输入得到对应的输出偏振态,而不是根据需要的偏振态去设计能满足需求的相位延迟器件。
因此本发明实施例将首先从直观,而非代数的层面出发,考虑相位延迟器件的性质对偏振态变化的影响。其规律可以总结为:
a)相位延迟器件在固定快轴角变化线性相位延迟时,偏振态在邦加球上变化的轨迹会在S1、S2平面投影成一段直线;
c)相位延迟器件在固定快轴角变化线性相位延迟时,偏振态在邦加球上变化的轨迹会形成一段圆弧。当线性相位延迟从0变化到2π时,圆弧将变为一个完整的圆。同时,圆弧的方向为从入射偏振态起,沿与S2方向成夹角(当S3<0时,为反向)作弧线。
再将相位延迟器件具体为固定相位延迟的波片与相位延迟可变的电控式元件(此处以液晶可变相位延迟器LCVR为例)。波片的特点是结构相对简单,无需配套的控制设备,可以很方便地与旋转电机组合以改变快轴角,但波片的相位延迟却无法自由地改变,通常是一个由波片本身材料决定的一个固定值。LCVR的特点与波片恰恰相反,它有着复杂的控制器及电路限制,无法与运动部件配合来使用,因此其快轴角通常在搭建完光路后便不再改变,但其相位延迟可以很方便地借助电压进行控制,且电控式导致其响应速度远快于电机旋转。
结合波片与LCVR的特性我们可以发现,LCVR器件在实际使用中是先确定一个快轴角后,根据需要可以人为地控制改变相位延迟,更加适合进行偏振态的调制。而波片是固定相位延迟,通过电机控制快轴角,在进行偏振态调制时,较难使其达到一个能满足需求的偏振态。同时,由于波片(旋转器件)在精度上的优势及LCVR在响应速度上的优势,往往不会将两者混合使用。同时考虑不同的成本及系统需求,我们以下列几种不同结构的测量光路为例,给出当两次测量时,各个元件的配置(偏振片、波片的角度,电控式相位可变延迟器的相位延迟)满足何种关系时,可以使产生的两个偏振态正交:
a)一个旋转偏振片及一个旋转1/4波片。
式中,θP为第一次测量配置的偏振片透光轴方向,θR为第一次测量配置的1/4波片快轴角方向,θP'为第二次测量配置的偏振片透光轴方向,θR'为第二次测量配置的1/4波片快轴角方向。
该结构主要通过偏振片调整起偏的主通光方向,再通过1/4波片实现特定偏振态的调制。一个旋转偏振片及一个旋转1/4波片的仪器示意图如图3所示。
b)一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片。
式中,θH为第一次测量配置的1/2波片快轴角方向,θR为第一次测量配置的1/4波片快轴角方向,θH'为第二次测量配置的1/2波片快轴角方向,θR'为第二次测量配置的1/4波片快轴角方向。
该结构与上一结构类似,不过线偏振光方向的调制从单偏振片调制变成了固定偏振片加一个1/2波片,可以避免入射光不是理想的自然光时,转动偏振片导致的不同方向的线偏振光光强不一致的情况。一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片的仪器示意图如图4所示。
c)一个固定偏振片,两个全波延迟相位调制器件。
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟。该结构首先通过第一个全波延迟相位调制器件对偏振态在邦加球上的S1OS3平面的圆上调制,随后借助第二个全波延迟相位调制器件得到目标偏振态。一个固定偏振片,两个全波延迟相位调制器件的仪器示意图如图5所示。
当两次测量配置满足以上关系时,出射偏振态在邦加球上呈对称的关系,及该两次测量的S1、S2、S3分量的和分别为零,即可满足最优仪器矩阵的第一个条件。
但由于四次测量时,该方法调制出的四种偏振态两两正交,这样两组调制出的四个偏振态将形成一个条件数无限大的仪器矩阵,使缪勒矩阵无解。因此,两两正交的关系仅对四点以上,且测量次数为偶次的仪器矩阵有效。四点仅存在两种最优的仪器矩阵:
这两种仪器矩阵在邦加球上恰好形成一个正四面体,且三视图均为正方形,在对行和为零进行限制时,将不再采取两两正交的方式,而是通过正四面体的特性进行限制。正四面体的三个视图均为正方形,因此如果对三视图进行约束,使测量系统产生的四个偏振态能够在邦加球上形成一个三视图均为矩形的四点测量相比多点测量来说,有着绝对的测量速度上的优势,此处给出四点测量下的最优测量框架配置:
a)一个旋转偏振片及一个旋转1/4波片。
式中,θP为第一次测量配置的偏振片透光轴方向,θR为第一次测量配置的1/4波片快轴角方向,θP'为第二次测量配置的偏振片透光轴方向,θR'为第二次测量配置的1/4波片快轴角方向,θP”为第三次测量配置的偏振片透光轴方向,θR”为第三次测量配置的1/4波片快轴角方向,θP”'为第四次测量配置的偏振片透光轴方向,θR”'为第四次测量配置的1/4波片快轴角方向。当四次测量的配置满足该式给出的关系时,可使四面体的三视图均为矩形。
b)一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片。
式中,θH为第一次测量配置的1/2波片快轴角方向,θR为第一次测量配置的1/4波片快轴角方向,θH'为第二次测量配置的1/2波片快轴角方向,θR'为第二次测量配置的1/4波片快轴角方向,θH”为第三次测量配置的1/2波片快轴角方向,θR”为第三次测量配置的1/4波片快轴角方向,θH”'为第四次测量配置的1/2波片快轴角方向,θR”'为第四次测量配置的1/4波片快轴角方向。
c)一个固定偏振片,两个全波延迟相位调制器件。
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″是第三次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″是第三次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″′是第四次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″′是第四次测量配置的第二个全波延迟相位调制器件的相位延迟。
d)一个固定偏振片及两个半波延迟相位调制器件。
式中,δH1、δH2分别为第一次测量配置的两个半波延迟相位调制器件的线性相位延迟,θH1、θH2则分别为其快轴角,δH1′、δH2′分别为第二次测量配置的两个半波延迟相位调制器件的线性相位延迟,δH1″、δH2″分别为第三次测量配置的两个半波延迟相位调制器件的线性相位延迟,δH1″′、δH2″′分别为第四次测量配置的两个半波延迟相位调制器件的线性相位延迟。
该配置无法满足全偏振态的调制,但可以满足四点的最优框架。一个固定偏振片及两个半波延迟相位调制器件的仪器示意图如图6所示,该仪器下可调制出的所有偏振态及最优的四点测量框架图如图7a至7b所示,图7a至7b中r深色部分为一个22.5度偏振片及两个半波延迟相位调制器件(67.5度、90度)结构的起偏器/检偏器所能产生的所有偏振态。这种结构能产生的偏振态无法遍历整个邦加球,但恰好可以覆盖到四点理论最优仪器矩阵的四个偏振态(如浅色正四面体的四个顶点)。
当满足以上条件时,可以将四点配置的自由度缩减到一个或两个,大幅减少配置优化的难度,借助该关系,可以快速得出最优框架的实际配置。
通过上述约束,无论是对四次测量还是多次测量,均可以得到泊松噪声与样本无关的优化配置。此时我们获得了N/2组正交偏振态(N次起偏/检偏测量)或一组满足一定关系的4个未知偏振态(4次起偏/检偏测量)。但仅知道每一组内部的关系,而实际的偏振态与配置仍是未知的。此时我们需要对每一组设置一系列变量,并根据该变量得到一组中剩余的可以由该变量明确表示的偏振态。之后,这些偏振态将组成一个有数个未知变量的仪器矩阵,仪器矩阵本身可以视为由多个偏振态组成,即如果有四个斯托克斯矢量S1,S2,S3,S4;则这四个列向量组成的仪器矩阵就为W=[S1 S2 S3 S4],无论这四个斯托克斯矢量为完全已知,还是由未知变量表示,前者可以得到一个完全已知的仪器矩阵W,后者可以得到一个由变量表示的仪器矩阵W。通过遗传算法对仪器矩阵进行最小EWV寻优,即可计算出当这数个未知变量为何值时,仪器矩阵的EWV能达到最小值。此时的仪器矩阵以及测量配置就是明确且已知的,未知变量都已经通过该算法计算出来。
在一些实施例中,当起偏器和检偏器分别为一个旋转偏振片及一个旋转1/4波片时,根据其对应的几何约束条件,设置一系列变量并针对变量进行EMW寻优的过程如下:
首先确定第一组正交偏振态:设第一组正交偏振态中,生成第一个偏振态时,偏振片的通光轴角度为θ1P;四分之一波片(以下简称波片)的快轴角为θ1R。则生成与其正交的偏振态时,根据给出的几何约束的关系,偏振片的角度应为θ1P'=θ1P+90;波片的角度应为θ1R'=θ1R。如果要生成八次测量的仪器矩阵,需要同理生成另外3组正交偏振态,其偏振片与波片角度分别为:
θ2P、θ2R;θ2P+90、θ2R
θ3P、θ3R;θ3P+90、θ3R
θ4P、θ4R;θ4P+90、θ4R
这样就获得了由θ1P、θ1R;θ2P、θ2R;θ3P、θ3R;θ4P、θ4R八个变量表示的四组正交偏振态。每一组中的两个偏振态都可以由每组对应的变量表示。之后将这八个变量产生的一系列偏振态组成仪器矩阵,就可以得到由八个变量表示的EWV。然后再对得到的EWV借助遗传算法进行寻优。本发明实施例的结果是由遗传算法计算得到,但其它的不同寻优算法也能实现这一步,即对一个由多个可变量计算出的结果进行最小化寻优这一问题。
%遗传算法示例如下:
%%该语句直接调用matlab自带的遗传算法函数ga,下面对该函数的使用进行说明:
[OutputVariables,EWV]=ga(@(UnknowVariables)fitfunc(UnknowVariables),8,[],[],[],[],[0 0 0 0 0 0 0 0],[180 180 180 180 180 180 180 180]);
%OutputVariables,EWV为该遗传算法的最终输出,OutputVariables为优化出的配置设置,EWV为优化出的配置所形成的仪器矩阵的等权重方差。
%ga函数的第一个输入fitfunc为自己定义的目标函数,其内容及定义方法在代码部分40-75行详细介绍。
%ga函数的第二个输入8为变量的数量。在本例中,有8个变量,因此输入设置为8。
%ga函数的第三到六个输入为对遗传算法本身的一些参数设定,此处选为默认设定,及[]即可。
%ga函数的第七个输入为变量的最小值取值范围,此处8个全部给为0,即[0 0 00 0 0 0 0]。
%ga函数的第八个输入为变量的最大值取值范围,此处8个全部给为180,即[180180 180 180 180 180 180 180]。
%使八个变量的取值范围均为0到180度。
%根据得到的变量与几何约束的关系,计算出其余几个由变量表示的配置:
theat_R1=OutputVariables(1);
theat_P1=OutputVariables(2);
theat_P1_r=theat_P1+90;
theat_R1_r=theat_R1;
theat_R2=OutputVariables(3);
theat_P2=OutputVariables(4);
theat_P2_r=theat_P2+90;
theat_R2_r=theat_R2;
theat_R3=OutputVariables(5);
theat_P3=OutputVariables(6);
theat_P3_r=theat_P3+90;
theat_R3_r=theat_R3;
theat_R4=OutputVariables(7);
theat_P4=OutputVariables(8);
theat_P4_r=theat_P4+90;
theat_R4_r=theat_R4;
%将得到的配置与EWV显示出来。
disp(strcat('八次测量中偏振片的角度分别为',num2str(theat_P1),',',num2str(theat_P1_r),',',num2str(theat_P2),',',num2str(theat_P2_r),',',num2str(theat_P3),',',num2str(theat_P3_r),',',num2str(theat_P4),',',num2str(theat_P4_r)));
disp(strcat('八次测量中波片的角度分别为',num2str(theat_R1),',',num2str(theat_R1_r),',',num2str(theat_R2),',',num2str(theat_R2_r),',',num2str(theat_R3),',',num2str(theat_R3_r),',',num2str(theat_R4),',',num2str(theat_R4_r)));
disp(strcat('最终得到的仪器矩阵的EWV为',num2str(EWV)));
%该函数为自己定义的目标函数,输入为变量,输出为优化目标的值:
function Fitness=fitfunc(variables)
%首先对一些基础值进行定义,包括偏振片缪勒矩阵,波片缪勒矩阵,入射自然光,方便计算Stokes:
MLR=@(x,y)[1,0,0,0;
0,cosd(2*x)^2+sind(2*x)^2*cosd(y),cosd(2*x)*sind(2*x)*(1-cosd(y)),-sind(2*x)*sind(y);
0,cosd(2*x)*sind(2*x)*(1-cosd(y)),sind(2*x)^2+cosd(2*x)^2*cosd(y),cosd(2*x)*sind(y);
0,sind(2*x)*sind(y),-cosd(2*x)*sind(y),cosd(y)];%波片穆勒矩阵公式,x为快轴角,y为相位延迟;
MLD=@(x)0.5*[1,cosd(2*x),sind(2*x),0;
cosd(2*x),cosd(2*x)^2,-cosd(2*x)*sind(2*x),0;
sind(2*x),sind(2*x)*cosd(2*x),sind(2*x)^2,0;
0,0,0,0];%偏振片穆勒矩阵公式,x为主通光轴角度;
inputlight=[1;0;0;0];%输入为一束没有偏振效应的自然光。
%之后计算出一系列Stokes。
%第一组Stokes:
Stokes1=MLR(variables(1),90)*MLD(variables(2))*inputlight;%得到由第1和第2个变量表示的Stokes1。
Stokes1_r=MLR(variables(1),90)*MLD(variables(2)+90)*inputlight;%得到由第1和第2个变量表示的且与Stokes1正交的Stokes1_r,下同,不再赘述。
%第二组Stokes:
Stokes2=MLR(variables(3),90)*MLD(variables(4))*inputlight;
Stokes2_r=MLR(variables(3),90)*MLD(variables(4)+90)*inputlight;
%第三组Stokes:
Stokes3=MLR(variables(5),90)*MLD(variables(6))*inputlight;
Stokes3_r=MLR(variables(5),90)*MLD(variables(6)+90)*inputlight;
%第四组Stokes:
Stokes4=MLR(variables(7),90)*MLD(variables(8))*inputlight;
Stokes4_r=MLR(variables(7),90)*MLD(variables(8)+90)*inputlight;
%得到所有Stokes后,将其组成一个仪器矩阵:
InstrumentMatrix=[Stokes1 Stokes1_r Stokes2 Stokes2_r Stokes3Stokes3_r Stokes4 Stokes4_r];
%计算该仪器矩阵的EWV:
ewv=trace((pinv(InstrumentMatrix))'*pinv(InstrumentMatrix));
%将优化目标确定为EWV。
Fitness=ewv;
end
综上,本发明实施例的起偏器与检偏器配置优化方法流程图如图8所示。
本发明实施例系统地提出一种对偏振测量系统(起偏器/检偏器)进行优化的方法,该方法关键是对能满足缪勒矩阵测量时,16个阵元的估计泊松噪声方差与样本无关的仪器矩阵给出现实可实现的、具体的仪器配置。现有的测量大部分未考虑到对泊松噪声的抑制或者仅从测量原理出发,说明在什么仪器矩阵下可以抑制噪声,并未给出可以参考或可借助现有器件实际实现的配置。本发明实施例通过上述关系进行几何约束,从而计算出具体测量时的仪器配置,以及满足上述关系,由该方法得到对应的配置变形。
下面给出本发明实施例的实施流程及最终结果。
首先确定所要采用的测量系统结构,除已在上文给出的一种特殊结构外,该系统无论是作为起偏器还是检偏器,均需要满足全偏振态的调制。其次确定对四次采集测量还是多次采集测量进行优化,针对这两种情况确定对每组偏振态的约束条件,借助不同偏振元件改变偏振态的规律,将每组偏振态用测量系统的实际配置进行约束。之后将经过约束后,用实际配置表示的偏振态组成一个仪器矩阵,并利用遗传算法或类似的寻优算法对该仪器矩阵进行最小EWV优化。即可得到最优的实际配置。
经过本发明实施例的起偏器与检偏器配置优化方法给出的部分起偏检偏系统的配置如下:
a)一个旋转偏振片及一个旋转1/4波片。
b)一个固定偏振片、一个旋转1/2波片及一个旋转1/4波片。
c)一个固定偏振片,两个全波延迟相位调制器件。
d)一个固定偏振片及两个半波延迟相位调制器件。
上述a)至d)的配置均满足起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,且起偏器和检偏器仪器矩阵的EWV是最优,以优化缪勒测量系统抗高斯-泊松混合噪声的性能。
当采用相应的结构进行测量时,需要把各个元件的配置按表中给出参数进行设置即可。例如,当采用旋转偏振片和旋转波片进行测量时,可以采用a)表所给出的第一种四点配置,即,首先使偏振片旋转到4.87度,波片旋转到22.5度进行第一次起偏;再使偏振片旋转到85.13度,波片旋转到67.5度进行第二次起偏;然后使偏振片旋转到94.87度,波片旋转到112.5度进行第三次起偏;最后使偏振片旋转到175.13度,波片旋转到157.5度进行第四次起偏。同时,每次起偏也需要进行一次检偏,检偏的过程与起偏类似,根据使用的检偏器的结构,选择结构对应的表中给出的具体配置,进行四次或八次测量检偏。即,如果起偏、检偏都采用四次测量的模式,总共需要进行16次测量。
与传统的测量配置相比,本发明实施例的起偏器与检偏器配置优化方法系统性地将测量配置与高斯-泊松混合噪声的优化相结合,能够对实际测量产生指导作用,使测量系统达到最优。传统的没有采用最佳测量配置的测量系统,虽然依然能够进行测量,但是其精度有限,无法很好地抑制系统中的噪声,而且噪声规律也并不明确,难以进行后期的计算降噪等工作。该方法优化出的测量配置则改善了这一问题,一是使系统中的噪声传递最小,最大程度地抑制了噪声,二是使噪声的规律与样本独立,无论测量何种样本,噪声的分布规律都是相同的。同时,该方法也具备普适性,可以对各种不同结构的测量系统进行优化。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的技术人员来说,在不脱离本发明构思的前提下,还可以做出若干等同替代或明显变型,而且性能或用途相同,都应当视为属于本发明的保护范围。
Claims (15)
1.一种起偏器与检偏器配置优化方法,其特征在于,包括如下步骤:
调整起偏器的仪器矩阵W和检偏器的仪器矩阵A,使起偏器和检偏器仪器矩阵的等权重方差EWV最小,以实现针对高斯噪声的优化;且
起偏器的仪器矩阵W和检偏器的仪器矩阵A每一行的和为0,使泊松噪声所引起的估计方差独立于样本,且估计方差达到最小值。
2.如权利要求1所述的起偏器与检偏器配置优化方法,其特征在于,配置测量系统的起偏器以及检偏器,使其偏振态两两正交;随后搜索起偏器和检偏器仪器矩阵的最小等权重方差EWV;当起偏器和检偏器的仪器矩阵同时满足行和为0及最小EWV条件时,得到最优配置。
5.如权利要求2所述的起偏器与检偏器配置优化方法,其特征在于,当起偏器和检偏器系统包括一个固定偏振片,两个全波延迟相位调制器件时,使起偏器和检偏器偏振态两两正交的配置条件满足以下关系:
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟。所述方法首先通过第一个全波延迟相位调制器件对偏振态在邦加球上的S1OS3平面的圆上调制,随后借助第二个全波延迟相位调制器件得到目标偏振态。
8.如权利要求1所述的起偏器与检偏器配置优化方法,其特征在于,当起偏器和检偏器系统包括一个固定偏振片,两个全波延迟相位调制器件时,使起偏器和检偏器最优的四点测量配置满足以下关系:
式中,δF1是第一次测量配置的第一个全波延迟相位调制器件的相位延迟,θF1是其快轴角,δF2是第一次测量配置的第二个全波延迟相位调制器件的相位延迟,θF2是其快轴角,δF1′是第二次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2′是第二次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″是第三次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″是第三次测量配置的第二个全波延迟相位调制器件的相位延迟,δF1″′是第四次测量配置的第一个全波延迟相位调制器件的相位延迟,δF2″′是第四次测量配置的第二个全波延迟相位调制器件的相位延迟。
10.如权利要求3-9中任一项所述的起偏器与检偏器配置优化方法,其特征在于,用实际配置表示的偏振态组成一个仪器矩阵,并利用遗传算法或类似的寻优算法对该仪器矩阵进行最小EWV优化,即,对每一组设置一个变量,并根据该变量得到一组中剩余的能够由该变量明确表示的偏振态;之后,这些偏振态将组成一个有数个未知变量的仪器矩阵;通过遗传算法对仪器矩阵进行最小EWV寻优,即可计算出当这数个未知变量为何值时,仪器矩阵的EWV能达到最小值。
11.如权利要求3-8中任一项所述的起偏器与检偏器配置行优化方法,其特征在于,起偏器和检偏器均满足全偏振态的调制。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2022/106351 WO2023001121A1 (zh) | 2021-07-19 | 2022-07-19 | 起偏器与检偏器配置优化方法及起偏检偏系统 |
US18/415,388 US20240167941A1 (en) | 2021-07-19 | 2024-01-17 | Polarizer and analyzer configuration optimization method and polarizing and analyzing system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210707222 | 2022-06-21 | ||
CN2022107072221 | 2022-06-21 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115272041A true CN115272041A (zh) | 2022-11-01 |
Family
ID=83767587
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210840465.2A Pending CN115272041A (zh) | 2021-07-19 | 2022-07-18 | 起偏器与检偏器配置优化方法及起偏检偏系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115272041A (zh) |
-
2022
- 2022-07-18 CN CN202210840465.2A patent/CN115272041A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9976906B2 (en) | Light polarization state modulation and detection apparatuses and detection method | |
US20240167941A1 (en) | Polarizer and analyzer configuration optimization method and polarizing and analyzing system | |
CN107367329A (zh) | 一种图像、光谱、偏振态一体化获取装置及探测方法 | |
CN110346655A (zh) | 一种天线的极化参数测量装置及方法 | |
WO2020120942A1 (en) | Single-shot mueller matrix polarimeter | |
Duivenvoorden et al. | Full-sky beam convolution for cosmic microwave background applications | |
CN111307286B (zh) | 一种基于介质柱结构的大角度偏振探测超表面 | |
CN112629662A (zh) | 一种可重构分时偏振分析系统及探测方法 | |
CN114136894B (zh) | 一种基于涡旋波片的偏振检测系统的误差校准方法及装置 | |
JP2003337068A (ja) | 偏光計の校正方法 | |
US11656484B2 (en) | Voltage-tunable polarizer | |
CN115272041A (zh) | 起偏器与检偏器配置优化方法及起偏检偏系统 | |
de Juan Ovelar et al. | Modeling the instrumental polarization of the VLT and E-ELT telescopes with the M&m's code | |
CN108534993A (zh) | 一种液晶可变相位延迟器偏振特性检测方法及系统 | |
CN113624690B (zh) | 一种能够实现任意偏振态调制的起偏器及其优化方法 | |
Harrington et al. | Polarization Modeling and Predictions for DKIST Part 4: Calibration Accuracy Over Field of View, Retardance Spatial Uniformity And Achromat Design Sensitivity | |
CN114152578A (zh) | 基于涡旋波片的空间调制偏振检测方法 | |
Vaughn | Bandwidth and noise in spatiotemporally modulated Mueller matrix polarimeters | |
Wang et al. | Equivalent Mueller matrix method for 3-D axial error analysis in 2-D SoP measurement | |
CN114018830B (zh) | 一种基于液晶偏振光栅的线偏振方向检测方法 | |
CN111238634B (zh) | 一种基于干涉的光束轨道角动量谱分析方法与装置 | |
Sharma et al. | Mahakala: a Python-based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Space-times | |
CN114488623B (zh) | 一种级联式液晶偏振光栅衍射效率的计算方法 | |
Michielsen et al. | Covariance operators, Green functions, and canonical stochastic electromagnetic fields | |
CN105136298B (zh) | 一种差分快照式成像方法 |
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 |