摘要
首先介绍了偏振探测器的结构,建立了探测器的信号响应模型,并基于该响应模型制定了相应的参数定标方案,实现了对入射辐射偏振态的高精度计算和对偏振图像的非均匀性校正。实验数据表明,校正后计算得到的线偏振度是真实值的97.8%~101.5%,满足了偏振探测的精度要求;强度图像和线偏振度图像相对于未校正图像的非均匀性分别降低了93.64%和93.67%,有效提高了图像质量。
偏振是光的一个基本特
焦平面偏振探测器的非均匀性严重影响了偏振信息的计算和偏振图像的质量,因此必须对探测器参数进行标定和非均匀性校正。Powel
焦平面偏振探测器是将微偏振片阵列直接集成到探测器的焦平面上,每个像元对应着一个微偏振片,构成偏振像元。在偏振像元阵列中,将4个(2×2)相邻的像元组合成一个超像元,如

图1 相邻的四像元合成超像元
Fig.1 One super-pixel consists of four neighboring pixels
. | (1) |
得到Stokes矢量后,可通过式(2)计算得到线偏振度和偏振角。
. | (2) |
本文采用的探测器是国产自研的InGaAs近红外线列焦平面偏振探测器(0.9~1.7 μm),探测器包括红外像元和偏振像元,像元尺寸为30 μm×30 μm,探测器的排列结构如

图2 探测器的结构
Fig.2 Structure of the detector
由于微纳加工的工艺限制,在微偏振片阵列中,微偏振片的透射率、消光比和偏振角(微偏振片的偏振轴与参考方向的夹角)等存在一定的差异,影响了探测器的性能。为了降低因光衍射效应产生的像元串扰,偏振像元上集成的微偏振片并没有完全覆盖红外像元,微偏振片与红外像元的面积比为0.87:1,减弱了微偏振片的偏振调制效果。在不考虑像元间串扰的条件下,本文的校正方法可以消除微偏振片阵列引起的误差。
在偏振像元阵列中,入射辐射透过微偏振片,被焦平面上的像元吸收,如

图3 入射辐射进入偏振像元示意图
Fig.3 The incident radiation enters the polarization pixel diagram
对于单个偏振像元,其入射辐射用Stokes矢量表示,微偏振片的出射辐射强度 的表达式如式(3)所示。
, | (3) |
式中为微偏振片的透射率,为消光比,为偏振角。探测器的响应只与入射辐射强度相关,对辐射的偏振特性不敏感。探测器的响应与辐射强度的关系如式(4)所示。
, | (4) |
式中为探测器偏振像元的响应值,为像元的增益系数,为探测器的量子效率,A为探测器的接受面积, 为单位立体角,h为普朗克常数,v为辐射频率,b为像元的偏置系数。将式(3)代入式(4)中,可以得到:
, | (5) |
在式(5)中,像元的增益系数与微偏振片的透射率同时作用于探测器的响
, | (6) |
通过以上分析可知,在不考虑像元串扰等因素的情况下,偏振像元的响应会受到微偏振片的透射率、消光比、偏振角以及像元的增益系数和偏置系数的影响。对于传统的成像系统,探测器的非均匀性只与像元的增益系数和偏置系数有关;而偏振成像系统多了一个微偏振片阵列环节,在微偏振片阵列中,相同方向偏振轴的微偏振片因制造加工误差,导致微偏振片的透射率、消光比以及偏振角存在差异,进而使得偏振探测器的非均匀性问题更加严重。
偏振探测器定标和非均匀性校正的目的是:对于同一均匀的同一强度入射辐射(无偏光),偏振探测器像元的响应一致;对于同一均匀的偏振入射辐射,探测器计算出来的偏振态是一致
为了让偏振探测器的像元响应与接收到辐射强度成线性关系,需要先标定像元的累积增益和偏置系数。本文采用多点校正法标定了探测器累积增益和偏置系数。
定标采用高温黑体作为辐射源,黑体辐射可以视为均匀无偏光。将探测器正对着高温黑体,将黑体设置多个不同的温度,分别采集数据。根据斯蒂芬-玻尔兹曼定律如式(7)计算出不同温度下探测器响应波段的辐射出射度,采用多点校正法计算出探测器在不同的响应范围内的累积增益和偏置系数。
. | (7) |
校正得到偏振探测器像元的累积增益和偏置系数后,令 ,则有:
. | (8) |
对于微偏振片阵列中的微偏振片,其参数是不变的,因此可以输入多组不同的偏振态辐射来计算微偏振片的参数,如式(9)所示。利用一个高消光比的偏振片放置在偏振探测器和高温黑体之间,保持黑体的温度不变,旋转偏振片的角度,得到多组不同的入射偏振态。
. | (9) |
利用最小二乘法进行拟合,可以得到微偏振的参数,如式(10)所示,式中+表示伪逆矩阵。
. | (10) |
对于一个超像元,在不考虑瞬时视场误差(IFOV)的情况下,则有:
. | (11) |
在式(11)中,表示超像元的4个相邻微偏振片的偏振特性矩阵即超像元的分析矩阵,(i=0, 45,90,135)为超像元的响应值,和(i=0,45,90,135)为超像元4个相邻微偏振片的参数。
根据式(11),通过计算超像元的分析矩阵的伪逆矩阵,可以得到入射辐射的Stokes矢量,如式(12)所示,+表示伪逆矩阵。
. | (12) |
对于理想的焦平面偏振探测器,探测器的微偏振片参数和响应均是理想的,对于相同偏振角的偏振像元,其响应是一致的。校正后探测器像元的响应值可以利用理想的超像元参数和计算得到,即超像元校正,如式(14)所示。
, | (13) |
. | (14) |
探测器校正装置如

图4 校正装置图
Fig.4 Correction device diagram
采集原始数据、处理数据和校正数据流程如下:
1) 将探测器正对着高温黑体(不经过偏振片),黑体设置成不同的温度(260~400℃,温度间隔为10℃),每个温度采集256次数据,数据深度为16位;
2) 为了降低探测器的随机噪声对定标参数的影响,将256次数据求平均,得到15组数据,根据式(7)计算出15个温度的辐射出射度。采用多点校正法计算出每个像元累积增益和偏置系数;
3) 将高温黑体的温度固定在380 ℃,偏振片固定在可精确控制旋转角度精密转台上,置于探测器和黑体之间。转台旋转5°采集一组数据(256次数据),将256次数据取平均,得到72组数据;
4) 步骤2)计算出的累积增益G和偏置系数b校正探测器采集的72组数据,偏振片可视为理想偏振片,计算出偏振片在不同角度时的入射Stokes矢量,再利用式(10)计算出每个像元的微偏振片参数;
5) 校正图像具体过程:对原始图像采用多点校正参数校正,除去每个像元的累积增益G和偏置系数b,再利用式(12)计算得到入射Stokes矢量估计;
6) 利用式(2)计算偏振度和偏振角。
探测器的响应曲线如

图5 探测器的响应曲线及不同方法的拟合曲线
Fig.5 Response curve of detector and fitting curve of different methods
为了评价不同拟合曲线的精度,采用平均误差(MAE)和均方误差(MSE)两种方法进行分析,MAE和MSE的计算公式如式(15)(16)所示。
, | (15) |
, | (16) |
式中M表示采集数据的组数,表示探测器的响应值,表示不同方法的拟合值。相对于原始数据,不同方法拟合曲线的MAE和MSE值如
在高温黑体温度为380 ℃的条件下,探测器和黑体之间的偏振片处在不同角度,角度间隔为5°,原始数据的四通道(超像元中四个不同角度的像元)响应如

图6 原始数据四通道响应
Fig.6 Raw data four-channel response

图7 响应校正后四通道响应
Fig.7 Four-channel response after non-uniform correction

图8 超像元校正后四通道响应
Fig.8 Four-channel response after super-pixel correction

图9 校正前和校正后计算出的线偏振度
Fig.9 DOLP calculated before and after calibration
为了验证校正的效果,定量评估校正模
, | (17) |
式(17)中表示图像的均值,表示偏振像元的总数,表示坏元和过热像元的数量,NU值越大表明图像不均匀性越严重,理想的NU值为0。探测器的偏振强度图像和线偏振度图像的NU计算结果如
未校正的偏振强度图像NU平均值为2.65%,线偏振度图像的NU平均值为13.60%,表明图像存在严重的不均匀性;在响应校正后,偏振强度图像和线偏振度图像的NU平均值分别为1.0083%和5.18%,图像不均匀性明显降低;偏振校正后,偏振强度图像和线偏振度图像的NU平均值分别为0.17 %和0.86%。
偏振校正后的偏振强度图像和线偏振度图像NU平均值相对于响应校正后的图像分别降低了83.30%和83.40%,说明微偏振片参数对偏振图像非均匀有着严重的影响;相对于未校正的图像分别降低了93.64%和93.67%,因此本文校正方法可以明显降低偏振图像的非均匀性。
将线列近红外偏振探测器固定在水平转台上,转台均速旋转(推扫式)即可采集得到原始图像数据如

图10 实景图
Fig.10 Actual scene pictures
本文采用了国产自研的InGaAs近红外线列焦平面偏振探测器,针对探测器的特点,建立了探测器偏振像元的响应模型,根据响应模型采用多点校正法和带角度约束的最小二乘法分别标定了探测器的响应参数和微偏振片阵列参数,修正了探测器的响应非均匀性、微偏振片阵列透射率与消光比不均匀和微偏振片角度误差。利用探测器的定标参数计算得到的线偏振度值是真实线偏振度值的97.8%~101.5%,提高了探测器的偏振探测精度;同时校正了偏振强度图像和线偏振度图像的非均匀性,校正后的图像偏振强度图像和线偏振度图像的NU平均值分别为0.17 %和0.86%,相对于未校正图像分别降低了93.64%和93.67%,对真实场景进行成像,图像质量得到了有效地提升。实验数据表明本文方法能够有效提高探测器的偏振探测精度,降低偏振图像的非均匀性,提升偏振成像质量。
References
Tyo J S, Goldstein D L, Chenault D B, et al. Review of passive imaging polarimetry for remote sensing applications[J]. Appl Opt, 2006, 45(22):5453-5469. [百度学术]
Meng L, Kerekes J P. An analytical model for optical polarimetric imaging systems[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10):6615-6626. [百度学术]
Bostater C R, Mertikas S P, Neyt X, et al. Remote sensing of the ocean, sea ice, coastal waters, and large water regions 2016 - Infrared polarimetric sensing of oil on water[J]. SPIE Proceedings,2016, 9999:99990D. [百度学术]
Fang S, Xia X S, Huo X, et al. Image dehazing using polarization effects of objects and airlight[J]. Optics Express, 2014, 22(16):19523. [百度学术]
Qi J, Elson D S. Mueller polarimetric imaging for surgical and diagnostic applications: a review[J]. Journal of Biophotonics, 2017, 10: 950-982. [百度学术]
Snik F, Craven-Jones J, Escuti M, et al. An overview of polarimetric sensing techniques and technology with applications to different research fields[C]. Polarization: Measurement, Analysis, and Remote Sensing XI. International Society for Optics and Photonics, 2014. [百度学术]
WANG Xia, XIA Run-Qiu, JIN Wei-Qi, et al. Technology process of infrared polarization detection[J]. Infrared and Laser Engineering (王霞, 夏润秋, 金伟其, 等. 红外偏振成像探测技术进展. 红外与激光工程), 2014, 43(10):3175-3182. [百度学术]
Powell S B, Gruev V. Calibration methods for division-of-focal-plane polarimeters[J]. Optics Express, 2013, 21(18):21039. [百度学术]
Huang F, Fan-Ming L, Wei-Cong C, et al. Calibration method for division of focal plane polarimeters[J]. Applied Optics, 2018, 57(18):4992-. [百度学术]
Feng B, Shi Z, Liu H, et al. Polarized-pixel performance model for DoFP polarimeter [J]. Journal of Optics, 2018, 20(6):065703. [百度学术]
LIU Hai-Zhen, SHI Ze-Lin, FENG Bin. Mechanism and correction of infrared imaging non-uniformity in micro-polarizer arrays[J]. Optics and Precision Engineering(刘海峥, 史泽林, 冯斌. 微偏振片阵列红外成像非均匀性产生机理及其校正. 光学精密工程) 2018, 26(2): 480-491. [百度学术]
PENG Yong, FENG Bin, SHI Ze-lin, et al. Non-uniformity correction in polarization imaging obtained with integrated microgrid polarimeters integrated microgrid polarimeters [J]. Infrared and Laser Engineering (彭勇, 冯斌, 史泽林, 等. 微偏振片阵列成像的非均匀校正研究. 红外与激光工程), 2017,46(4):89-96. [百度学术]
Bowers, David L. Unpolarized calibration and nonuniformity correction for long-wave infrared microgrid imaging polarimeters[J]. Optical Engineering, 2008, 47(4):046403. [百度学术]