网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

近红外偏振探测器的定标及非均匀性校正  PDF

  • 尹佳琪 1,2,3
  • 王世勇 1,3
  • 张瑞 1,2,3
  • 李范鸣 1,3
1.中国科学院上海技术物理研究所,上海 200083; 2.中国科学院大学,北京 100049; 3.中国科学院红外探测与成像技术重点实验室,上海 200083

中图分类号: TN215

最近更新:2020-04-23

DOI:10.11972/j.issn.1001-9014.2020.02.011

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

首先介绍了偏振探测器的结构,建立了探测器的信号响应模型,并基于该响应模型制定了相应的参数定标方案,实现了对入射辐射偏振态的高精度计算和对偏振图像的非均匀性校正。实验数据表明,校正后计算得到的线偏振度是真实值的97.8%~101.5%,满足了偏振探测的精度要求;强度图像和线偏振度图像相对于未校正图像的非均匀性分别降低了93.64%和93.67%,有效提高了图像质量。

引言

偏振是光的一个基本特[

1]。偏振图像可以从场景中获得目标场景的矢量信息,提供强度和光谱图像无法获取的目标信[2]。偏振探测技术在环境监[3]、图像去[4]、生物医[5]、遥[6]等领域有着广泛的应用。目前,偏振成像系统主要包括分时、分振幅、分孔径、焦平面等类型。焦平面偏振探测器具有实时好、紧凑的结构和高可靠性等优势,但是存在明显的非均匀性问[7]

焦平面偏振探测器的非均匀性严重影响了偏振信息的计算和偏振图像的质量,因此必须对探测器参数进行标定和非均匀性校正。Powell [

8]、Huang F[9]和Feng B[10]等采用最小二乘法,刘海[11]和彭[12]等采用矩阵法标定了偏振探测器的参数,提高了偏振探测器的探测精度和偏振图像的均匀性。但是探测器的响应与入射辐射不是完全呈线性关系,在不同的工作环境中探测器的参数是不同的,需要多次定标,大大增加了校正的工作量。本文基于近红外线列偏振探测器的特点,建立了偏振探测器像元的响应模型,并根据偏振探测器的响应模型制定了对探测器参数的定标方案,提出了响应非均匀性和微偏振片参数不一致性分别校正的方法。

1 偏振探测器模型

1.1 近红外线列偏振探测器

焦平面偏振探测器是将微偏振片阵列直接集成到探测器的焦平面上,每个像元对应着一个微偏振片,构成偏振像元。在偏振像元阵列中,将4个(2×2)相邻的像元组合成一个超像元,如图1所示,四个像元对应的微偏振片分别是0°、45°、90°、135°的金属光栅,用于测量偏振光在这四个角度的线偏振光分量,分别记为I0I45I90I135,并采用这四个强度值计算Stokes矢量S=S0S1S2T,如式(1)所示。

图1 相邻的四像元合成超像元

Fig.1 One super-pixel consists of four neighboring pixels

S0=0.5(I0+I90+I45+I135)S1=I0-I90S2=I45-I135
.
(1)

得到Stokes矢量后,可通过式(2)计算得到线偏振度和偏振角。

DOLP=S12+S22S0AOLP=12arctan(S2S1)
.
(2)

本文采用的探测器是国产自研的InGaAs近红外线列焦平面偏振探测器(0.9~1.7 μm),探测器包括红外像元和偏振像元,像元尺寸为30 μm×30 μm,探测器的排列结构如图2所示。探测器有270×4个像元,第一列和第四列为红外像元,第二列和第三列偏振像元。

图2 探测器的结构

Fig.2 Structure of the detector

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

1.2 偏振像元响应模型

在偏振像元阵列中,入射辐射透过微偏振片,被焦平面上的像元吸收,如图3所示。探测器将接受到的辐射进行光电转化,形成电信号,再通过A/D转换得到数字响应值。

图3 入射辐射进入偏振像元示意图

Fig.3 The incident radiation enters the polarization pixel diagram

对于单个偏振像元,其入射辐射用Stokes矢量Sin表示,微偏振片的出射辐射强度Lout 的表达式如式(3)所示。

Lout=ττε-1ε+1cos(2θ)τε-1ε+1sin(2θ)TSin
(3)

式中τ为微偏振片的透射率,ε为消光比,θ为偏振角。探测器的响应只与入射辐射强度相关,对辐射的偏振特性不敏感。探测器的响应与辐射强度的关系如式(4)所示。

I=gLout+b
(4)

式中I为探测器偏振像元的响应值,g=QeAΩhv为像元的增益系数,Qe为探测器的量子效率,A为探测器的接受面积,Ω 为单位立体角,h为普朗克常数,v为辐射频率,b为像元的偏置系数。将式(3)代入式(4)中,可以得到:

I=gττε-1ε+1cos(2θ)τε-1ε+1sin(2θ)TSin+b
(5)

在式(5)中,像元的增益系数与微偏振片的透射率同时作用于探测器的响[

13],因此可以简单地认为是偏振像元的累积增益G(像元的增益系数和微偏振片透射率的乘积即G=gτ)。令p=ε-1ε+1cos(2θ)q=ε-1ε+1sin(2θ)则:

I=G1pqSin+b
(6)

通过以上分析可知,在不考虑像元串扰等因素的情况下,偏振像元的响应会受到微偏振片的透射率、消光比、偏振角以及像元的增益系数和偏置系数的影响。对于传统的成像系统,探测器的非均匀性只与像元的增益系数和偏置系数有关;而偏振成像系统多了一个微偏振片阵列环节,在微偏振片阵列中,相同方向偏振轴的微偏振片因制造加工误差,导致微偏振片的透射率、消光比以及偏振角存在差异,进而使得偏振探测器的非均匀性问题更加严重。

2 偏振探测器的校正方法

偏振探测器定标和非均匀性校正的目的是:对于同一均匀的同一强度入射辐射(无偏光),偏振探测器像元的响应一致;对于同一均匀的偏振入射辐射,探测器计算出来的偏振态是一致[

11]。在不同的积分时间或不同的入射辐射强度情况下,偏振探测器像元的增益系数和偏置系数是变化的,而微偏振片的透射率、消光比和偏振角等参数保持不变,因此偏振探测器的定标包括对偏振像元的累积增益G、偏置系数b的定标和微偏振片的消光比和偏振角等参数的定标,探测器的非均匀性校正分别为响应校正(探测器响应的非均匀性校正)和偏振校正(定标的微偏振片参数校正计算出入射偏振态)。

2.1 偏振探测器的定标和非均匀性校正

为了让偏振探测器的像元响应与接收到辐射强度成线性关系,需要先标定像元的累积增益G和偏置系数b。本文采用多点校正法标定了探测器累积增益G和偏置系数b

定标采用高温黑体作为辐射源,黑体辐射可以视为均匀无偏光。将探测器正对着高温黑体,将黑体设置多个不同的温度,分别采集数据。根据斯蒂芬-玻尔兹曼定律如式(7)计算出不同温度下探测器响应波段的辐射出射度,采用多点校正法计算出探测器在不同的响应范围内的累积增益G和偏置系数b

M=λ1λ2c1λ5[exp(c2λT)-1]dλ
.
(7)

校正得到偏振探测器像元的累积增益G和偏置系数b后,令Y=(I-b)G ,则有:

Y=1pqSin
.
(8)

对于微偏振片阵列中的微偏振片,其参数是不变的,因此可以输入多组不同的偏振态辐射来计算微偏振片的参数,如式(9)所示。利用一个高消光比的偏振片放置在偏振探测器和高温黑体之间,保持黑体的温度不变,旋转偏振片的角度,得到多组不同的入射偏振态。

Y1Y2YN=S1inTS2inTSNinT1pq
.
(9)

利用最小二乘法进行拟合,可以得到微偏振的参数,如式(10)所示,式中+表示伪逆矩阵。

1pq=S1inTS2inTSNinT+Y1Y2YN
.
(10)

对于一个超像元,在不考虑瞬时视场误差(IFOV)的情况下,则有:

Y=Y0Y45Y90Y135=1p0q01p45q451p90q901p135q135Sin=ASin
.
(11)

在式(11)中,ARn×4表示超像元的4个相邻微偏振片的偏振特性矩阵即超像元的分析矩阵,Yii=0, 45,90,135)为超像元的响应值,piqii=0,45,90,135)为超像元4个相邻微偏振片的参数。

根据式(11),通过计算超像元的分析矩阵A的伪逆矩阵,可以得到入射辐射的Stokes矢量Sin,如式(12)所示,+表示伪逆矩阵。

Sin=A+Y=1p0q01p45q451p90q901p135q135+Y0Y45Y90Y135
.
(12)

对于理想的焦平面偏振探测器,探测器的微偏振片参数和响应均是理想的,对于相同偏振角的偏振像元,其响应是一致的。校正后探测器像元的响应值Yc可以利用理想的超像元参数AidealSin计算得到,即超像元校正,如式(14)所示。

Aideal=12×1101011-1010-1
(13)
Yc=Yc0Yc45Yc90Yc135=AidealSin
.
(14)

2.2 原始数据采集与校正流程

探测器校正装置如图4所示,包括高温黑体、可精确控制旋转角度精密转台、高消光比的金属线栅偏振片(Thorlabs公司提供的WP50L-UB,探测器响应波段的消光比>1000)、光学镜头和近红外线列偏振探测器。

图4 校正装置图

Fig.4 Correction device diagram

采集原始数据、处理数据和校正数据流程如下:

1) 将探测器正对着高温黑体(不经过偏振片),黑体设置成不同的温度(260~400℃,温度间隔为10℃),每个温度采集256次数据,数据深度为16位;

2) 为了降低探测器的随机噪声对定标参数的影响,将256次数据求平均,得到15组数据,根据式(7)计算出15个温度的辐射出射度。采用多点校正法计算出每个像元累积增益G和偏置系数b

3) 将高温黑体的温度固定在380 ℃,偏振片固定在可精确控制旋转角度精密转台上,置于探测器和黑体之间。转台旋转5°采集一组数据(256次数据),将256次数据取平均,得到72组数据;

4) 步骤2)计算出的累积增益G和偏置系数b校正探测器采集的72组数据,偏振片可视为理想偏振片,计算出偏振片在不同角度时的入射Stokes矢量,再利用式(10)计算出每个像元的微偏振片参数;

5) 校正图像具体过程:对原始图像采用多点校正参数校正,除去每个像元的累积增益G和偏置系数b,再利用式(12)计算得到入射Stokes矢量估计;

6) 利用式(2)计算偏振度和偏振角。

3 实验结果与分析

3.1 实验数据分析

探测器的响应曲线如图5所示。从图中可以看出探测器的响应值与入射辐射值并不完全呈线性关系,采用两点校正参数、多点校正参数和带有二次项的非线性校正参数进行拟合探测器的响应曲线,可以看出两点校正参数拟合的曲线误差最大,多点校正参数和带有二次项的非线性校正参数拟合的曲线误差较小。

图5 探测器的响应曲线及不同方法的拟合曲线

Fig.5 Response curve of detector and fitting curve of different methods

为了评价不同拟合曲线的精度,采用平均误差(MAE)和均方误差(MSE)两种方法进行分析,MAEMSE的计算公式如式(15)(16)所示。

MAE=1M1iMO(i)-F(i)
(15)
MSE=1M1iMO(i)-F(i)2
(16)

式中M表示采集数据的组数,O(i)表示探测器的响应值,F(i)表示不同方法的拟合值。相对于原始数据,不同方法拟合曲线的MAE和MSE值如表1所示。从表中数据可以看出,多点校正参数拟合值的平均误差和均方误差值均明显小于非线性校正参数和两点校正参数的拟合值,因此探测器的响应校正应采用线性的多点校正法。

在高温黑体温度为380 ℃的条件下,探测器和黑体之间的偏振片处在不同角度,角度间隔为5°,原始数据的四通道(超像元中四个不同角度的像元)响应如图6所示,可以看出探测器四通道的响应差异很大,存在严重的非均匀性;图7是探测器响应校正后的四通道响应,利用多点校正法校正原始数据,可以看出探测器四个通道的响应差异明显减小,但是峰谷和峰顶与探测器表面的微偏振片角度存在偏差;图8是在响应校正的基础上,利用超像元校正计算出的四通道响应,校正后的数据具有较好的均匀性,且校正了角度误差。

图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是探测器的响应校正和偏振校正后的偏振度值。在响应校正后,计算出的偏振度平均值为44.6%,是理论值的43.3%~46.6%,与实际值差别很大;偏振校正后,计算的偏振度值平均值为100.0%,偏振度的值在97.8%~101.5%之间波动。

图9 校正前和校正后计算出的线偏振度

Fig.9 DOLP calculated before and after calibration

为了验证校正的效果,定量评估校正模[

10]。在均匀的入射辐射条件下,改变偏振片的线偏振角,捕获原始偏振图像数据,计算强度图像和偏振度图像的均匀性。计算如式(17)所示。

NU=1R¯1M-Ndk=1MR(k)-R¯2
(17)

式(17)中R¯表示图像的均值,M表示偏振像元的总数,Nd表示坏元和过热像元的数量,NU值越大表明图像不均匀性越严重,理想的NU值为0。探测器的偏振强度图像和线偏振度图像的NU计算结果如表2表3所示。

表2表3列出了不同偏振态辐射下的强度图像和线偏振度图像的NU值,表中的Ixx=0, 30, 60, 90,120,150)表示置于探测器前方偏振片角度为x,线偏振度为1的偏振图像。

未校正的偏振强度图像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%,因此本文校正方法可以明显降低偏振图像的非均匀性。

3.2 成像校正结果

将线列近红外偏振探测器固定在水平转台上,转台均速旋转(推扫式)即可采集得到原始图像数据如图10(a)所示,图10(a)的图像包含偏振四个通道的数据,可以利用超像元计算图像的偏振信息。超像元计算图像的偏振信息会导致图像的行数减至为原始图像的一半,大大降低了图像的分辨率。为了不降低图像的分辨率,本文采用偏振像元复用的方式,如图2所示,第一行和第二行的偏振像元构成第一个超像元,而第二行和第三行的偏振像元可以构成第二个超像元,如此类推,计算得到偏振图像的行数不会减半,如图10(b)、 (d)和(f)是强度图像S0图10(c)、 (e)和(g)是线偏振度图像DOLP。

图10 实景图

Fig.10 Actual scene pictures

图10(b)和(c)是未校正的原始数据利用式(1)和式(2)计算得到的偏振强度图像和线偏振度图像,可以看出强度图像和线偏振度图像均存在明显的横条纹,严重影响了图像的质量;图10(d)和(e)是响应校正后的图像,与未校正的图像相比,偏振强度图像和线偏振度图像均匀性得到了明显的改善,图像质量有了显著地提升,但是线偏振度图像还是存在一定横条纹;图10(f)和(g)是校正微偏振片参数后得到的图像,偏振强度图像和线偏振度图像的均匀性得到了进一步的改善,无肉眼可见的横条纹,图像质量进一步改善。

4 结论

本文采用了国产自研的InGaAs近红外线列焦平面偏振探测器,针对探测器的特点,建立了探测器偏振像元的响应模型,根据响应模型采用多点校正法和带角度约束的最小二乘法分别标定了探测器的响应参数和微偏振片阵列参数,修正了探测器的响应非均匀性、微偏振片阵列透射率与消光比不均匀和微偏振片角度误差。利用探测器的定标参数计算得到的线偏振度值是真实线偏振度值的97.8%~101.5%,提高了探测器的偏振探测精度;同时校正了偏振强度图像和线偏振度图像的非均匀性,校正后的图像偏振强度图像和线偏振度图像的NU平均值分别为0.17 %和0.86%,相对于未校正图像分别降低了93.64%和93.67%,对真实场景进行成像,图像质量得到了有效地提升。实验数据表明本文方法能够有效提高探测器的偏振探测精度,降低偏振图像的非均匀性,提升偏振成像质量。

References

1

Tyo J S, Goldstein D L, Chenault D B, et al. Review of passive imaging polarimetry for remote sensing applications[J]. Appl Opt, 2006, 4522):5453-5469. [百度学术

2

Meng L, Kerekes J P. An analytical model for optical polarimetric imaging systems[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 5210):6615-6626. [百度学术

3

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 Proceedings2016, 9999:99990D. [百度学术

4

Fang S, Xia X S, Huo X, et al. Image dehazing using polarization effects of objects and airlight[J]. Optics Express, 2014, 2216):19523. [百度学术

5

Qi J, Elson D S. Mueller polarimetric imaging for surgical and diagnostic applications: a review[J]. Journal of Biophotonics, 2017, 10: 950-982. [百度学术

6

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. [百度学术

7

WANG Xia, XIA Run-Qiu, JIN Wei-Qi, et al. Technology process of infrared polarization detection[J]. Infrared and Laser Engineering (王霞, 夏润秋, 金伟其, . 红外偏振成像探测技术进展. 红外与激光工程), 2014, 4310):3175-3182. [百度学术

8

Powell S B, Gruev V. Calibration methods for division-of-focal-plane polarimeters[J]. Optics Express, 2013, 2118):21039. [百度学术

9

Huang F, Fan-Ming L, Wei-Cong C, et al. Calibration method for division of focal plane polarimeters[J]. Applied Optics, 2018, 5718):4992-. [百度学术

10

Feng B, Shi Z, Liu H, et al. Polarized-pixel performance model for DoFP polarimeter [J]. Journal of Optics, 2018, 206):065703. [百度学术

11

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, 262): 480-491. [百度学术

12

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 (彭勇, 冯斌, 史泽林, . 微偏振片阵列成像的非均匀校正研究. 红外与激光工程), 2017464):89-96. [百度学术

13

Bowers, David L. Unpolarized calibration and nonuniformity correction for long-wave infrared microgrid imaging polarimeters[J]. Optical Engineering, 2008, 474):046403. [百度学术