基于涡旋波片的空间调制偏振检测方法

1.本发明属于偏振测量的技术领域，具体涉及一种基于涡旋波片的空间调制偏振检测方法。


背景技术：

2.偏振测量用于表征光束的偏振态，光与样品相互作用会改变入射光的偏振态，样品对光的偏振态的改变可用stokes参量表示，stokes参量可以分为圆退偏参量、线退偏参量等与样品微结构密切相关的、具有实际物理意义的、可量化的偏振参数，可用于获取样品的偏振特性和结构参数。偏振测量作为光和样品偏振特性分析的重要工具，已经被广泛的应用于生物医学、量子通信、激光雷达等领域典型的偏振仪可分为时序偏振仪、同时偏振仪。时序偏振仪利用延迟器和偏振分析仪、液晶调制器或光弹调制器通过一系列分时测量来获取stokes参量，这种测量方法无法实现实时的偏振检测。同时测量方法包括分振幅法、分焦面法、分孔径法是一种并行测量方法，将带测光光束分成多个通道进行探测，这些通道可以在空间频域内通过分光器件获得，每个检测通道仅检测某种特定的偏振态。然而上述偏振测量技术需多次器件旋转或多路分光而无法精确测量弱光场偏振信息。
3.最近，基于涡旋波片的空间调制偏振测量法被提出，该方法具有简单、准确、鲁棒性强等优势。论文文献“spatially modulated polarimetry based on a vortex retarder and fourier analysis”使用涡旋二分之一波片将入射光偏振态映射到一幅空间变化的强度分布，通过数据处理算法解算出实验测得强度分布所表征的待测偏振态。然而，由于涡旋二分之一波片物理属性的限制，该偏振仪仅能测量完全偏振光无法测量部分偏振光。论文文献“single-shot measurement of polarization state at low light field using mueller-mapping star test polarimetry”提出基于涡旋四分之一波片空间调制的偏振仪解决了上述问题，并发展了弱光偏振测方法。由于该方法建立探测面的强度分布同空间调制器mueller矩阵之间的关系，因此理论模型的标定至少需要20次旋转波片，标定过程步骤复杂耗时。此外，该方法使用直接枚举法求解归一化的最小二乘问题，求解耗时且精度受限于stokes参量的采样间隔。因此，现有的基于涡旋波片的空间调制偏振测法精度和速度有待进一步提高。上述方法仅适用于偏振态的测量方法，并不适用于弱光场mueller矩阵测量。


技术实现要素：

4.有鉴于此，本发明提供了一种基于涡旋波片的空间调制偏振检测方法，能够实现强/弱光场下的光子偏振态、光子偏振几率及mueller矩阵的检测，解决了现有空间调制偏振检测技术测量和标定速度慢，弱光偏振检测精度低的问题。
5.实现本发明的技术方案如下：
6.基于涡旋波片的空间调制stokes检测方法，建立任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系，通过四次旋转偏振器件(旋转3次单
快轴波片和1次偏振片)获取所述四种特定偏振光探测面的强度分布，实现检测系统的加工和对准误差的标定；
7.建立光子偏振态及偏振几率的理论模型，将实验获得的待测光场强度分布输入所述理论模型中，利用混合梯度下降算法求解最大的光子偏振几率及所对应的stokes参量，所述混合梯度下降算法即为枚举算法—梯度下降算法—枚举算法的结合算法。
8.进一步地，所述四种特定偏振光探测面为水平、+45度线偏振态和左、右圆偏振态入射时探测面。
9.进一步地，所述任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系为：
[0010][0011]
其中，i表示任意偏振光入射时探测面的强度分布，ih、i
p
、i
l
、ir分别表示水平、+45度线偏振态和左、右圆偏振态入射时探测面的强度分布，s1,s2,s3表示入射光归一化的偏振态。
[0012]
进一步地，所述方法对于待测弱光场的强度分布，使用阈值滤波对弱光图像进行去噪处理：
[0013][0014]
其中，ii是实验测得的弱光场下的待测偏振光强度分布在i像素点的强度值，n表示总像素数，λ表示可调参数，snr表示弱光强度分布的信噪比。
[0015]
进一步地，所述检测方法中，采用了基于涡旋四分之一波片的空间调制偏振仪。
[0016]
进一步地，所述混合梯度下降算法具体为：
[0017]
步骤1、使用直接枚举算法搜索初始的stokes参量的最优解，为步骤2梯度下降算法提供初始值；
[0018]
步骤2、将步骤1所得的初始值代入梯度下降算法迭代搜索stokes参量的最优解，每次搜索stokes参量的最优解需要更新偏振几率，当偏振几率的梯度为正时梯度下降的方向被设置为1，反之梯度下降的方向被设置为-1；
[0019]
步骤3、利用步骤2的结果使用直接枚举算法获得测量的stokes参量及最大的光子偏振几率。
[0020]
基于涡旋波片的空间调制mueller检测方法，建立任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系，通过四次旋转偏振器件(旋转3次单快轴波片和1次偏振片)获取所述四种特定偏振光探测面的强度分布，实现检测系统的加工和对准误差的标定；
[0021]
用多快轴1/4波片替换单快轴波片，用于产生入射样品的4种线性无关的偏振态，经过样品后产生4种出射偏振态；每一种所述出射偏振态依据以下检测方法求解；
[0022]
建立光子偏振态及偏振几率的理论模型，将每一种出射偏振态输入所述理论模型中，利用混合梯度下降算法求解最大的光子偏振几率及所对应的stokes参量，所述混合梯
度下降算法即为枚举算法—梯度下降算法—枚举算法的结合算法；
[0023]
进而获得样品的mueller矩阵。
[0024]
有益效果：
[0025]
第一，本发明选用基于涡旋四分之一波片的空间调制偏振仪，其可同时实现完全偏振光和部分偏振光偏振信息的检测，有效解决了现有基于涡旋二分之一波片的空间调制偏振仪无法测量部分偏振态的问题，结合阈值选通的低通傅立叶滤波技术，所述方法可以实现更低信噪比下的偏振检测，有望应用于活体生物偏振检测。
[0026]
第二，本发明建立的正向理论模型即任意偏振光的强度分布同四种特殊偏振态入射的强度分布之间的关系，该模型标定时仅需要四次旋转偏振器件获取四幅探测面的强度分布，简化了标定的步骤，节省了标定的时间，提高了系统的稳定性。
[0027]
第三，本发明建立了光子偏振态及偏振几率模型，不仅能解算待测光的偏振态而且根据偏振几率能够提高先验信息。
[0028]
第四，在数据处理端口，本发明提出了混合梯度下降算法来求解最大的偏振几率，结合了直接枚举算法和梯度下降算法的各自的优势，显著提高了计算待测光偏振态的速度和精度，有望实现实时偏振检测。
[0029]
第五，本发明在进行待测样品的mueller矩阵的测量过程中，使用水平偏振片和多快轴四分之一波片作为偏振态产生器，可同时产生四组线性无关的入射偏振态，经过涡旋四分之一波片空间调制后，相机单词拍摄可解算出出射样品的四种偏振态，加快了样品mueller矩阵测量的测量速度。
[0030]
第六，本发明系统集成度高，初始系统校准方便快捷，性能稳定，相机单次拍摄就能够获取待测光束的偏振信息；成本较低，可广泛用于快速、高精度的检测待测光偏振态。
附图说明
[0031]
图1为本发明偏振态及偏振几率检测实验系统光路图。
[0032]
图2为本发明样品mueller矩阵检测实验系统光路图。
[0033]
其中，101-光源，102-准直透镜，103-偏振态产生器，104-涡旋四分之一波片，105-偏振态分析器，106-聚焦透镜，107-探测面，108-电脑，109-待测样品。
具体实施方式
[0034]
下面结合附图并举实施例，对本发明进行详细描述。
[0035]
本发明提供了一种基于涡旋波片的空间调制偏振测量方法，具体过程包括：
[0036]
1、建立偏振检测理论模型即任意偏振态入射时探测面的强度分布同四种特殊偏振态入射的探测面的强度分布之间的关系。测量四种特殊的入射偏振态(水平、+45
°
线偏振态和左右圆偏振态)的探测面的强度分布，将上述实验获得的四幅强度分布带入所建立的理论模型中，以标定器件加工和对准所引起的测量误差。
[0037]
2、本发明建立了弱光场下图像的去噪方法，以实现高精度弱光场偏振检测。
[0038]
3、本发明建立光子偏振态及偏振几率模型用于高精度检测强/弱光场的偏振态及偏振几率。
[0039]
4、本发明提出了混合梯度下降算法求解偏振几率最大值及沿着梯度下降的方向
搜索最优的stokes参量。
[0040]
5、使用偏振片和多快轴四分之一波片构成偏振态产生器，可同时产生四种线性无关的入射偏振态，通过本发明方法可实现相机单次拍摄解算待测mueller矩阵。
[0041]
本发明建立理论模型即任意偏振态入射时探测面的强度分布同四种特殊偏振态入射时探测面的强度分布之间的关系，具体步骤为：
[0042]
步骤一，搭建偏振检测系统，635nmled光源由光纤导出后由准直透镜扩束，经偏振态产生器产生待测偏振态(用于核验所述方法测量stokes参量的准确性)，待测光束先后经过涡旋四分之一波片、偏振态分析器(由水平线偏振片和快轴沿水平方向放置的四分之一波片组成)构成的空间调制偏振检测系统，最后被聚焦透镜将光束聚焦到探测面，探测面的强度分布经过电脑数据处理输出待测光stokes参量的测量结果。
[0043]
步骤二，建立理论模型即任意偏振态入射时探测面的强度分布同四种特殊偏振态入射的探测面的强度分布之间的关系，待测光经过涡旋四分之一波片、偏振分析器、透镜聚焦后到探测面的强度分布可表示为：
[0044][0045]
这里i
out
表示探测面的强度分布，m
l
表示透镜将光束聚焦到探测器的mueller矩阵，m
p
表示偏振态分析器偏振片的mueller矩阵，mq表示偏振态分析器波片的mueller矩阵，m
zvqr
表示涡旋四分之一波片的mueller矩阵，表示待测光归一化的stokes参量。上述四个偏振器件的mueller矩阵可一般化为：
[0046][0047]
进一步的公式(1)可表示为：
[0048]
i＝m
00
+m
01
·
s1+m
02
·
s2+m
03
·
s3(3)
[0049]
四种特殊强度分布同mueller矩阵m关系为：
[0050]
ih＝m
00
+m
01
(4)
[0051]ip
＝m
00
+m
02
(5)
[0052]
ir＝m
00
+m
03
(6)
[0053]il
＝m
00-m
03
(7)
[0054]
这里ih、i
p
、i
l
、ir分别表示水平、+45度线偏振态和左右圆偏振态入射时探测面的强度分布，将上述公式带入公式(3)可得：
[0055][0056]
步骤三，将实验中测量的四幅探测面的强度分布带入公式(8)标定器件加工和对准引起的测量误差。
[0057]
此外，本发明建立了阈值选通的傅里叶低通滤波去噪算法，具体过程如下：
[0058]
在探测器获取了待测弱光场强度分布后，使用阈值滤波对弱光图像进行去噪处理：
[0059][0060]
其中，ii是实验测得的弱光场下的待测偏振光强度分布在i像素点的强度值，n表示总像素数，λ表示可调参数，snr表示弱光强度分布的信噪比，表示为：
[0061][0062]
这里er表示弱光强度分布的均值，eb表示背景噪声的均值，
σb
表示背景噪声的方差。经过阈值去噪后的强度分布再经过傅里叶低通滤波去噪后解算出偏振态及偏振几率。
[0063]
本发明建立了光子偏振态及偏振几率模型，具体步骤为：
[0064]
给定入射光偏振态单个光子落入探测面的第i个像素点的几率为：
[0065][0066]ii
表示入射光偏振态为时，所获取的探测面的强度分布在第i个像素的强度值。当探测面每一个像素接受的强度分别为i1,i2...ii,...in,每一像素点落入的光子数为k1,k2,...ki,...kn,因此探测面形成一幅强度分布对应入射偏振态为的几率为:
[0067][0068]
其中对于均匀极化的待测光源，ω为常数，表示取得偏振态的概率，p(ii)表示光子会打到第i像素点的概率，ki表示待测偏振光在探测面的强度分布在第i个像素点的光字数。表示获得探测面的一幅强度分布，入射光为的偏振几率。
[0069][0070]
q表示偏振几率，当偏振几率取得最大值时，所对应的stokes参量作为偏振测量结果。
[0071]
本发明建立了求解最大偏振几率的混合梯度下降算法，具体步骤为：
[0072]
步骤一，使用直接枚举算法(设置大的stokes采样间隔和较少的样本点)搜索初始的stokes参量的最优解，为步骤二梯度下降算法提供初始值。
[0073]
步骤二，进一步使用梯度下降算法迭代搜索stokes参量的最优解，每次搜索stokes参量的最优解需要更新偏振几率，表示为:
[0074]
[0075][0076][0077]
这里是stokes参量sj在第k迭代中对应的偏振几率，表示在第k次迭代中的stokes参量，在k次迭代中stokes参量被表示为：
[0078][0079]
其中step表示迭代步长，表示偏振几率梯度下降的方向，被设置为-1或1。进一步的偏振几率的梯度可表示为：
[0080][0081]
这里，表示偏振几率的梯度对于在第k次迭代中的stokes参量sj，进一步的梯度方向被设置为：
[0082][0083]
当偏振几率的梯度为正时梯度下降的方向被设置为1，反之梯度下降的方向被设置为-1。
[0084]
步骤三，使用直接枚举算法(设置小的stokes采样间隔和较少的样本点)获得最终的stokes参量及最大的光子偏振几率。
[0085]
在此基础上本发明建立了单拍mueller矩阵测量方法，具体步骤为：
[0086]
使用水平偏振片和多快轴四分之一波片作为偏振态产生器，多快轴四分之一波片由四个不同快轴方向的四分之一波片组成，可同时产生四种线性无关的偏振态入射。
[0087][0088]
表示由偏振态产生器所产生的四组线性无关的入射偏振态，
[0089]
表示相应入射偏振态经过样品后的出射偏振态，因此样品的mueller矩阵表示为：
[0090][0091]
结合所述光子偏振态测量方法，实验中可同时测得四种出射样品的偏振态，即可计算出样品的mueller矩阵m。
[0092]
本发明实施例提供一种快速、高精度的空间调制stokes测量方法，如图1所示，具体过程为：
[0093]
步骤一，搭建stokes偏振检测系统如图1所示，625nm led(带宽17nm)光源101由光纤导出后由准直透镜102扩束，经偏振态产生器103产生待测偏振态(用于核验所述方法测量stokes参量的准确性)，待测光束先后经过涡旋四分之一波片104、偏振态分析器105(由水平线偏振片和快轴沿水平放置的四分之一波片组成)构成的空间调制偏振检测系统，最
后被聚焦透镜106将光束聚焦到探测面107，然经过电脑108数据处理。
[0094]
步骤二，测量四种特殊偏振态入射的探测面的强度分布，并带入理论模型，用于标定器件加工和对准引起的测量误差。
[0095][0096]
这里，i表示任意偏振光入射时探测面的强度分布，ih、i
p
、i
l
、ir分别表示水平、+45度线偏振态和左右圆偏振态入射时探测面的强度分布。
[0097]
将实验中测量的四幅探测面的强度分布带入上述公式标定器件加工和对准引起的测量误差。
[0098]
步骤三，解算待测光stokes参量需要求偏振几率模型，如下所示；
[0099]
给定入射光偏振态单个光子落入探测面的第i个像素点的几率为：
[0100][0101]ii
表示入射光偏振态为时，所获取的探测面的强度分布在第i个像素的强度值。当探测面每一个像素接受的强度分别为i1,i2...ii,...in,每一像素点落入的光子数为k1,k2,...ki,...kn,因此探测面形成一幅强度分布对应入射偏振态为的几率为:
[0102][0103]
其中对于均匀极化的待测光源，ω为常数，表示取得偏振态的概率，p(ii)表示光子会打到第i像素点的概率。表示获得探测面的一幅强度分布，入射光为的偏振几率。
[0104][0105]
q表示偏振几率，当偏振几率取得最大值时，所对应的stokes参量作为偏振测量结果。
[0106]
步骤四，使用直接枚举算法(大的stokes采样间隔和较少的样本点)搜索初始的stokes参量的最优解，为步骤五梯度下降算法提供初始值。
[0107]
步骤五，进一步使用梯度下降算法迭代搜索stokes参量的最优解，每次搜索stokes参量的最优解需要更新偏振几率，表示为:
[0108][0109][0110]
[0111]
这里是stokes参量sj在第k迭代中对应的偏振几率，表示在第k次迭代中的stokes参量，在k次迭代中stokes参量被表示为：
[0112][0113]
其中step表示迭代步长，表示偏振几率梯度下降的方向，被设置为-1或1.进一步的偏振几率的梯度可表示为：
[0114][0115]
这里，表示偏振几率的梯度对于在第k次迭代中的stokes参量sj，进一步的梯度方向被设置为：
[0116][0117]
当偏振几率的梯度为正时梯度下降的方向被设置为1，反之梯度下降的方向被设置为-1。
[0118]
步骤六、利用步骤五的结果使用直接枚举算法获得测量的stokes参量及最大的光子偏振几率。
[0119]
本发明实施例提供一种快速、高精度的强/弱光场mueller矩阵检测方法，搭建mueller偏振检测系统如图2所示，625nm led(带宽17nm)光源101由光纤导出后由准直透镜102扩束，经偏振态产生器103产生待测偏振态(透光轴为水平方向的偏振片和拼接的多快轴四分之一波片共同组成偏振态产生器)，产生四组已知的入射光偏振态入射到待测生物样品109，从样品出射的光束经过涡旋四分之一波片104、偏振态分析器105(由水平线偏振片和快轴沿水平放置的四分之一波片组成)构，最后被聚焦透镜106将光束聚焦到探测面107，然经过电脑108数据处理，由此构成mueller矩阵偏振检测系统。
[0120]
强/弱光场mueller矩阵测量具体实施案例：
[0121]
步骤一，偏振态产生器的波片由偏振片和多快轴四分之一波片构成，多快轴波片会同时产生四种线性无关的入射偏振态表示为：
[0122]
步骤二，入射光经过待测样品、涡旋四分之一波片和偏振态分析器后，多快轴四分之一波片的拼缝会将探测面的强度分布分成四个不同形态的强度图，分别对四个强度图单独数据处理，数据处理的方法是如上述求解最大偏振几率和偏振态的过程，进而求解出样品出射的四种偏振态
[0123]
步骤三，将入射样品的四种偏振态和对应的出射样品的四种偏振态带入公式(21)求解出样品的mueller矩阵。
[0124]
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.基于涡旋波片的空间调制stokes检测方法，其特征在于，建立任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系，通过四次旋转偏振器件获取所述四种特定偏振光探测面的强度分布，实现检测系统的加工和对准误差的标定；所述四次旋转偏振器为旋转3次单快轴波片和1次偏振片；建立光子偏振态及偏振几率的理论模型，将实验获得的待测光场强度分布输入所述理论模型中，利用混合梯度下降算法求解最大的光子偏振几率及所对应的stokes参量，所述混合梯度下降算法即为枚举算法—梯度下降算法—枚举算法的结合算法。2.如权利要求1所述的基于涡旋波片的空间调制stokes检测方法，其特征在于，所述四种特定偏振光探测面为水平、+45度线偏振态和左、右圆偏振态入射时探测面。3.如权利要求1或2所述的基于涡旋波片的空间调制stokes检测方法，其特征在于，所述任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系为：其中，i表示任意偏振光入射时探测面的强度分布，i
h
、i
p
、i
l
、i
r
分别表示水平、+45度线偏振态和左、右圆偏振态入射时探测面的强度分布，s1,s2,s3表示入射光归一化的偏振态。4.如权利要求1所述的基于涡旋波片的空间调制stokes检测方法，其特征在于，所述方法对于待测弱光场的强度分布，使用阈值滤波对弱光图像进行去噪处理：其中，i
i
是实验测得的弱光场下的待测偏振光强度分布在i像素点的强度值，n表示总像素数，λ表示可调参数，snr表示弱光强度分布的信噪比。5.如权利要求1所述的基于涡旋波片的空间调制stokes检测方法，其特征在于，所述检测方法中，采用了基于涡旋四分之一波片的空间调制偏振仪。6.如权利要求1所述的基于涡旋波片的空间调制stokes检测方法，其特征在于，所述混合梯度下降算法具体为：步骤1、使用直接枚举算法搜索初始的stokes参量的最优解，为步骤2梯度下降算法提供初始值；步骤2、将步骤1所得的初始值代入梯度下降算法迭代搜索stokes参量的最优解，每次搜索stokes参量的最优解需要更新偏振几率，当偏振几率的梯度为正时梯度下降的方向被设置为1，反之梯度下降的方向被设置为-1；步骤3、利用步骤2的结果使用直接枚举算法获得测量的stokes参量及最大的光子偏振几率。7.基于涡旋波片的空间调制mueller检测方法，其特征在于，建立任意偏振光入射时探测面的强度分布同四种特定偏振光探测面的强度分布之间的关系，通过四次旋转偏振器件获取所述四种特定偏振光探测面的强度分布，实现检测系统的加工和对准误差的标定；所
述四次旋转偏振器为旋转3次单快轴波片和1次偏振片；用多快轴1/4波片替换单快轴波片，产生入射样品的4种线性无关的偏振态，经过样品后产生4种出射偏振态；每一种所述出射偏振态依据以下检测方法求解；建立光子偏振态及偏振几率的理论模型，将每一种出射偏振态输入所述理论模型中，利用混合梯度下降算法求解最大的光子偏振几率及所对应的stokes参量，所述混合梯度下降算法即为枚举算法—梯度下降算法—枚举算法的结合算法；进而获得样品的mueller矩阵。

技术总结
本发明公开了基于涡旋波片的空间调制偏振检测方法，建立任意偏振光入射时探测面的强度分布同四种特殊偏振光入射时探测面的强度分布之间的关系，所建立理论模型获取四幅探测面的强度分布可实现系统加工和对准误差的标定，简化了空间调制偏振仪的标定步骤，提高了偏振检测效率；提出了混合梯度下降算法来求解最大的光子偏振几率，该方法可实现高精度、快速的偏振态及偏振几率检测，提高了空间调制偏振仪的数据处理效率；结合阈值选通的傅里叶低通滤波技术，本发明方法通过相机单次拍摄可实现强/弱光场下的光子偏振态、光子偏振几率及穆勒矩阵的快速、高精度检测。高精度检测。高精度检测。


技术研发人员：李艳秋 宁天磊
受保护的技术使用者：北京理工大学
技术研发日：2021.11.26
技术公布日：2022/3/8