Underwater polarization imaging method based on active illumination modulation
-
摘要:
浑浊水体中的光学衰减和吸收效应常常导致成像质量的明显下降,对诸多基于水下视觉图像的应用产生了不可忽视的影响。介绍了一种水下散射抑制的偏振成像方法,其关键思想是通过主动照明调制实现后向散射光偏振程度最大化,从而最大程度地抑制后向散射光。通过目标场景穆勒矩阵的测量、后向散射光的斯托克斯矢量的给定、入射主动照明斯托克斯矢量的确定,以及偏振特征矢量的计算,最终结合这些信息对散射介质中的图像进行复原。实验结果表明,该方法在处理浑浊水体下的成像问题方面表现出色。与其他方法相比,它显著提升了图像质量,为水下环境中的成像应用提供了有力的解决方案。该方法的独特之处在于其主动偏振照明调制技术,通过主动照明控制,实现了后向散射光的偏振特性调控,从而有效地减少了散射引起的成像退化,为水下偏振成像领域带来了新的研究思路。
Abstract:The optical attenuation and absorption effects in turbid water often lead to a significant decrease in imaging quality, which has an undeniable impact on many applications based on underwater visual images. A polarization imaging method for underwater scattering suppression was proposed, the crucial idea of which was to maximize the degree of polarization of backscattered light through active illumination modulation, thus the backscattered light was suppressed to the greatest extent. Through the measurement of the Mueller matrix of the target scene, the given Stokes vector of the backscattered light, the determination of the Stokes vector of the incident active illumination, and the calculation of the polarization feature vectors, the images in the scattering medium were finally recovered by combining above information. Experimental results show that the method performs well in dealing with the imaging problems under turbid water. Compared with other methods, it significantly improves image quality and provides a powerful solution for imaging applications in underwater environments. The uniqueness of this method lies in its active illumination modulation technique, which achieves the modulation of the polarization characteristics of the backscattered light through active illumination control, thereby effectively reducing the imaging degradation caused by scattering, and bringing new research ideas to the field of underwater polarization imaging.
-
引言
飞机类目标作为现代战争的重要平台,广泛应用于侦察、攻击、运输等领域。基于此,其也成为红外探测系统、红外制导武器探测和打击的主要对象。准确获取来袭飞机的红外辐射特性,对于红外预警探测具有重要意义[1],因此,测量和研究飞机类目标的红外辐射特性是一项重要的基础工作,需长期进行。受保障条件限制,单独组织飞机类目标的红外辐射特性测量困难较大,一般会采取“搭车”方式测量。例如组织侦察探测实验时,利用已有的红外告警设备,结合实验开展目标红外辐射特性测量。由于红外告警设备并非专用的目标红外辐射特性测量设备,在具备原始图像采集的基础上,如何进行准确标定往往是测量前需要解决的关键问题[2]。
针对红外探测系统辐射标定的研究已有诸多报道,例如,吕原等提出了绝对辐射标定和基于多温度场的场地绝对辐射标定方法,实现了实时高精度绝对辐射标定[3]。蔡李靖等提出了一种无黑体式红外热像仪测温标定和温度补偿技术,并通过实验进行了验证[4]。凌玲等对傅里叶变换红外光谱仪开展了实验室标定,并分析了光谱仪测量过程中的不确定来源[5]。李宁等针对3 μm~5 μm红外焦平面阵列的标定拟合及评估,提出了利用剔除局外点的方法对拟合曲线进行优化处理,改善了辐射标定拟合的不确定度[6]。杨词银、李宪圣、李宁等针对大口径红外光电系统,分别提出了基于大面源黑体的辐射标定方法、用小口径照明光管标定大口径红外探测系统的方法和基于光栅单色仪和标准辐射计的相对光谱标定方法,并进行了实验验证及误差分析[7-10]。邹前进等采用不同的标定方法对焦距固定红外探测系统进行了标定实验,并对不同方法的适用情况进行了分析[11]。禄晓飞等基于红外标定参数受外场环境影响因素,提出了外场标定参数验证的方法流程,并进行了精度检验[12]。本文在此基础上,针对外场实际使用需求,分析验证了近距离直接扩展源法、远距离小源法对红外告警设备的标定误差,并给出了使用建议。
1 辐射标定及方法
标定一般采用红外探测系统测量标准红外辐射源的方式进行,目的是建立红外探测系统数字图像输出值(灰度)与其接收到的辐射亮度之间的对应关系,以便在目标特性测量过程中,准确得到目标的辐射亮度或强度。
标定过程首先将红外探测系统参数固定,并依次对不同温度的黑体进行探测,得到其在不同温度条件下的图像输出值(灰度)D,通过式(1)计算其辐射亮度响应度α和自身偏移值D0之间的对应关系,标定模型如下[13]:
$$ D = \alpha L + {D_0} $$ (1) 式中L为黑体在红外探测系统工作波段内的辐射亮度,不同温度对应的黑体辐射亮度通过式(2)计算得到[14]。
$$ L = \frac{\varepsilon }{\pi }\int\limits_{{\lambda _1}}^{{\lambda _2}} {\frac{{{{{c}}_1}}}{{{\lambda ^5}}} \cdot \frac{1}{{{{\mathrm{e}}^{{{{{{c}}_2}} / {\lambda T}}}} - 1}}{\rm{d}}\lambda } $$ (2) 式中:λ1~λ2为红外探测系统工作波段;ε为黑体发射率;c1、c2为辐射常数;T为黑体温度。
目前,针对红外探测系统的标定方法主要有近距离扩展源法、远距离小源法和平行光管法[15]。
近距离扩展源法一般采用大面源黑体作为标准红外辐射源,要求黑体辐射面覆盖红外探测系统光学镜头,并靠近红外探测系统光学镜头,以减小外界环境及背景的影响,该方法可以对红外探测系统探测器的每个像元进行标定。
当黑体辐射面不能覆盖红外探测系统光学镜头,或受现场条件限制,黑体无法靠近红外探测系统光学镜头时,可采用远距离小源法进行标定。远距离小源法一般将黑体置于距红外探测系统一定距离处,并保证红外探测系统能够对黑体清晰成像。
平行光管法需要使用大口径平行光管,标定过程首先将红外探测系统与平行光管对接,其次将黑体置于平行光管焦平面,形成无穷远辐射源。该方法由于标定系统较为复杂,一般不在外场使用。
本文主要探讨外场测量设备的标定,因此主要对近距离扩展源法和远距离小源法进行分析,不涉及平行光管法。
2 标定及结果
2.1 设备及参数
被标定设备为某长波红外告警设备,工作波段为8 μm~12 μm,像元尺寸为25 μm×25 μm,量化位数为16位。标准红外辐射源为以色列CI Systems公司生产的SR800扩展面黑体(见图1),参数见表1。
表 1 SR800扩展面黑体参数Table 1. Parameters of SR800 extended surface blackbody参 数 参数值 温度范围/℃ −40~170 温度精度/℃ ±0.02 辐射面积/mm2 300×300 发射率 0.97 标定时环境温度为28℃左右,为保证标定有效且标定温度分布在环境温度两侧,将标定温度设置在10℃~50℃之间。SR800扩展面黑体温度依次设置为10℃、20℃、30℃、40℃、50℃,利用长波红外告警设备对黑体进行测量,由式(1)、式(2)得到其辐射亮度响应度α和偏移值D0之间的对应关系。
2.2 近距离直接扩展源法标定结果
由于该长波红外告警设备采取旋转扫描方式对目标成像,标定过程中黑体无法紧贴其光学镜头,难以完全排除外界环境及背景的影响。近距离直接扩展源法的标定结果见图2,其中横坐标为黑体辐射亮度L,纵坐标为长波红外告警设备标定区域图像中心10×10个像元的输出灰度平均值。
对近距离直接扩展源法的标定结果进行线性拟合,得到长波红外告警设备的辐射亮度响应关系为D=225.07×L+
5973.7 。2.3 远距离小源法标定结果
标定时由于已将长波红外告警设备参数固定,距离太近无法对黑体清晰成像,距离太远大气及背景辐射对标定影响较大。因此,选取37 m和63 m两个能够清晰成像且距离较近的位置进行标定,37 m处黑体红外图像见图3。
标定结果如图4所示,其中横坐标为黑体辐射亮度L,纵坐标为长波红外告警设备标定区域图像中心10×10个像元的输出灰度平均值。
对远距离小源法的标定结果进行线性拟合,得到长波红外告警设备在37 m处的辐射亮度响应关系为D=198.63×L+
6031.9 ,在63 m处的辐射亮度响应关系为D=181.90×L+5701.2 。3 结果分析
从式(1)可以看出,在长波红外告警设备测量范围内,理论上2个不同的温度点就可以确定其辐射亮度响应度α和偏移值D0。表2~表4分别给出了由相邻温度点得到的近距离扩展源法、远距离小源法标定结果。
表 2 近距离直接扩展源法两点标定结果Table 2. Two points calibration results of short-range direct extended source method标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 232.3 231.5 218.4 222.2 偏移值D0 5873.8 5885.0 6112.5 6034.1 表 3 远距离小源法(37 m)两点标定结果Table 3. Two points calibration results of long-distance small source method (37 m)标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 193.0 206.3 194.4 199.2 偏移值D0 6100.5 5910.6 6117.7 6018.1 表 4 远距离小源法(63 m)两点标定结果Table 4. Two points calibration results of long-distance small source method (63 m)标定温度/℃ 10,20 20,30 30,40 40,50 响应度$\alpha $ 179.9 176.8 193.9 173.2 偏移值D0 5727.7 5772.1 5475.3 5907.5 表2中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.2节中拟合结果的相对误差分别在−2.96%~3.21%和−1.67%~2.32%之间。表3中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.3节中拟合结果的相对误差分别在−2.83%~3.86%和−2.01%~1.42%之间。表4中相邻温度点的辐射亮度响应度$\alpha $和偏移值D0与2.3节中拟合结果的相对误差分别在−4.78%~6.60%和−3.96%~3.62%之间。可以看出,近距离扩展源法的两点标定结果一致性较好,远距离小源法在37 m处的两点标定结果一致性也较好,但63 m处的两点标定结果一致性有所下降,说明随着距离的增加,外界环境及背景的影响有所增强。
在不考虑大气辐射传输的情况下,本次标定由设备带来的误差主要包括:SR800扩展面黑体光谱辐射出射度误差,长波红外告警设备光学系统透过率不确定度引起的误差,输出信号引起的误差,以及探测器响应不稳定引起的误差。SR800扩展面黑体温度精度为±0.02,结合文献[8,15-16]对黑体辐射出射度误差2%左右的分析,SR800扩展面黑体引起的测量误差${\sigma _1} $大约为2%。长波红外告警设备光学系统的光谱透过率实测误差${\sigma _2} $约为2%。根据当前长波红外告警设备技术状态,结合文献[15-16]的分析,其信号输出值D的不确定度${\sigma _3} $优于1%,探测器响应不稳定引起的误差${\sigma _4} $约为2%。
误差的合成标准不确定度可以表示为
$$ {\sigma }_{总}={\left({\sigma }_{1}^{2}+{\sigma }_{2}^{2}+{\sigma }_{3}^{2}+{\sigma }_{4}^{2}\right)}^{1/2} $$ (3) 经计算,由设备带来的总误差约为3.6%。
下面探讨不同标定方法带来的误差。为了验证2种标定结果的准确性,将SR800扩展面黑体作为目标布设于距长波红外告警设备140 m处进行辐射特性测量,利用式(4)对SR800扩展面黑体辐射亮度进行反演。
$$ {L_{\mathrm{F}}} = {{\left( {\frac{{D - {D_0}}}{\alpha } - {L_{{\mathrm{path}}}}} \right)} / {{\tau _\alpha }}} $$ (4) 式中:LF为反演的SR800扩展面黑体辐射亮度;τα为SR800扩展面黑体与长波红外告警设备之间的大气透过率;Lpath为SR800扩展面黑体与长波红外告警设备之间的大气程辐射,其结果通过大气辐射传输计算软件MODTRAN计算得到。测量期间,大气能见度为14.7 km,湿度为39%,温度为28℃,通过MODTRAN计算得到0 m、37 m、63 m距离的大气透过率分别为1、
0.9495 、0.9282 ;大气程辐射分别为0、0.5484 W∙m−2∙sr−1、0.7851 W∙m−2∙sr−1。表5~表7分别给出了采用近距离扩展源法、远距离小源法标定后的测量误差。表 5 近距离直接扩展源法测量误差Table 5. Measurement errors of short-range direct extended source method温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 10.97 −5.67 20 14.30 13.00 −9.09 30 17.36 15.36 −11.52 40 20.81 18.08 −13.12 50 24.67 20.79 −15.73 表 6 远距离小源法(37 m)测量误差Table 6. Measurement errors of long-distance small source method (37 m)温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 12.18 4.75 20 14.30 14.60 2.07 30 17.36 17.12 −1.38 40 20.81 20.19 −2.98 50 24.67 23.26 −5.72 表 7 远距离小源法(63 m)测量误差Table 7. Measurement errors of long-distance small source method (63 m)温度/℃ 黑体辐射亮度/
W∙m−2∙sr−1反演的辐射亮度/
W∙m−2∙sr−1相对
误差/%10 11.63 15.08 29.66 20 14.30 17.59 23.01 30 17.36 20.51 18.15 40 20.81 23.87 14.70 50 24.67 27.22 10.34 $$ \delta = \left( {\frac{{{L_{\mathrm{F}}} - {L_{\mathrm{H}}}}}{{{L_{\mathrm{H}}}}}} \right) \times 100{\text{%}} $$ (5) 式中LH为SR800扩展面黑体实际辐射亮度。
图5给出了采用近距离扩展源法、远距离小源法标定后的测量误差曲线。
1) 目标温度在10℃~50℃范围内,采用近距离直接扩展源法标定的测量误差在5.67%~15.73%之间。采用远距离小源法在37 m处标定的测量误差在1.38%~5.72%之间,在63 m处标定的测量误差在10.34%~29.66%之间。
2) 由于长波红外告警设备以周扫方式工作,无法对目标凝视成像,采用近距离直接扩展源法标定时,无法将黑体紧靠其光学系统,导致部分外界环境及背景辐射进入长波红外告警设备光学系统。另外,长波红外告警设备光学系统焦距较长,周扫过程中,由于角速度大对近距离放置的黑体探测稳定度较差。除3.6%的设备误差外,根据式(6),标定过程中环境及背景辐射、标定方法造成的误差可能在4.38%~15.31%之间。
$$ {\sigma }_{1}={\left({\sigma }_{总}^{2}\text-{\sigma }_{2}^{2}\right)}^{1/2} $$ (6) 3) 采用远距离小源法对长波红外告警设备进行标定,除3.6%的设备误差外,根据式(6),在37 m处,标定过程中环境及背景辐射、标定方法造成的误差优于4.45%,在63 m处,标定过程中环境及背景辐射、标定方法造成的误差可能在9.69%~29.44%之间。37 m处的标定误差明显优于63 m处的标定误差,由此可见,对于采取周扫方式工作的长波红外告警设备,远距离小源法是对其标定的有效手段。从验证结果可以看出,在保证长波红外告警设备能对黑体清晰成像的前提下,降低标定距离可有效减小环境及背景辐射对标定误差的影响。因此,在实际操作过程中,针对设备特点,找出其能对黑体清晰成像的最短距离,是远距离小源法标定的关键。
4 结论
本文利用长波红外告警设备、SR800扩展面黑体进行了外场标定实验,并验证了近距离直接扩展源法、远距离小源法的标定误差。结果表明,对于采取周扫方式工作的长波红外告警设备,近距离直接扩展源法误差较大,该方法适用性较差。远距离小源法是对周扫方式工作的长波红外告警设备进行标定的有效手段,但如何合理选择标定距离以提高其标定精度,是需要重点关注的问题。针对长波红外告警设备参数及特点,找出标定距离与标定误差的对应关系,总结规律并给出最优的标定方案,是未来研究的主要方向。
-
表 1 图6中光强图及不同算法客观评价指标
Table 1 Objective evaluation indicators for light intensity images and different algorithms in Fig. 6
指标 光强图 LRSD DCP Retinex Ours-1 Ours-2 EME 1.128 1.925 3.435 3.401 5.197 12.946 STD 0.017 0.040 0.048 0.033 0.041 0.067 Entropy 6.912 7.361 7.513 7.302 7.383 7.745 C 0.220 0.526 0.404 0.467 0.354 0.518 -
[1] JAFFE J S. Computer modeling and the design of optimal underwater imaging systems[J]. IEEE Journal of Oceanic Engineering, 1990, 15(2): 101-111. doi: 10.1109/48.50695
[2] HE D M, SEET G G L. Divergent-beam lidar imaging in turbid water[J]. Optics & Lasers in Engineering, 2004, 41(1): 217-231.
[3] HU H, YANG S, LI X, et al. Polarized image super-resolution via a deep convolutional neural network [J]. Optics Express, 31(5): 8535-8547.
[4] 梁健, 巨海娟, 张文飞, 等. 偏振光学成像去雾技术综述[J]. 光学学报, 2017, 37(4): 0400001. doi: 10.3788/AOS201737.0400001 LIANG Jian, JU Haijuan, ZHANG Wenfei, et al. Review of dehazing techniques for polarized optical imaging[J]. Acta Optica Sinica, 2017, 37(4): 0400001. doi: 10.3788/AOS201737.0400001
[5] LI X, XU J, ZHANG L, et al. Underwater image restoration via stokes decomposition[J]. Optics Letters, 2022, 47(11): 2854-2857. doi: 10.1364/OL.457964
[6] LI X B, HU H F, LIN Z, et al. Polarimetric image recovery method combining histogram stretching for underwater imaging[J]. Scientific Reports, 2018, 8(1): 1-10.
[7] LAND E H. The retinex theory of color vision[J]. Scientific American, 1977, 237(6): 108-128. doi: 10.1038/scientificamerican1277-108
[8] HE K, SUN J, TANG X. Single image haze removal using dark channel prior[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2010, 33(12): 2341-2353.
[9] TREIBITZ T, SCHECHNER Y Y. Active polarization descattering[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2009, 31(3): 385-399.
[10] TAN R T. Visibility in bad weather from a single image[C]//Proceedings of the 2008 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. Anchorage, AK, USA: IEEE, 2008: 1-8.
[11] CAI B, XU X, JIA K, et al. Dehazenet: an end-to-end system for single image haze removal[J]. IEEE Transactions on Image Processing, 2016, 25(11): 5187-5198.
[12] LI B, PENG X, WANG Z, et al. Aod-Net: all-in-one dehazing network [C]// Proceedings of the 2017 IEEE International Conference on Computer Vision (ICCV). Venice, Italy: IEEE, 2017: 4780-4788.
[13] ZHANG H, PATEL V M. Densely connected pyramid dehazing network[C]// Proceedings of the 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition. Salt Lake City, UT, USA: IEEE, 2018: 3194-3203.
[14] QI P, LI X, ZHAI J, et al. U2r-Pgan: unpaired underwater-image recovery with polarimetric generative adversarial network[J]. Optics and Lasers in Engineering, 2022, 157: 107112. doi: 10.1016/j.optlaseng.2022.107112
[15] NARASIMHAN S G, NAYAR S K. Vision and the atmosphere[J]. International Journal of Computer Vision, 2002, 48(3): 233-254. doi: 10.1023/A:1016328200723
[16] 周文舟, 范晨, 胡小平, 等. 多尺度奇异值分解的偏振图像融合去雾算法与实验 [J]. 中国光学, 2021, 14: 298-306. ZHOU Wenzhou , FAN Chen, HU Xiaoping, et al. Polarization image fusion dehazing algorithm and experiment based on multiscale singular value decomposition[J]. Chinese Optics, 2021, 14: 298-306.
[17] LIANG J, REN L, LIANG R. Low-pass filtering based polarimetric dehazing method for dense haze removal[J]. Optics Express, 2021, 29(18): 28178-28189. doi: 10.1364/OE.427629
[18] QU Y, ZOU Z. Non-sky polarization-based dehazing algorithm for non-specular objects using polarization difference and global scene feature[J]. Optics Express, 2017, 25(21): 25004-25022. doi: 10.1364/OE.25.025004
[19] SUN R, LIAO T, FAN Z, et al. Polarization dehazing method based on separating and iterative optimizing airlight from the frequency domain for different concentrations of haze[J]. Applied Optics, 2022, 61(35): 10362-10373.
[20] DONG Z M, ZHENG D F, HUANG Y T, et al. A polarization-based image restoration method for both haze and underwater scattering environment[J]. Scientific Reports, 2022, 12(1): 1836. doi: 10.1038/s41598-022-05852-1
[21] GOLDSTEIN D H. Polarized Light [M]. Boston: CRC Press Taylor&Francis Group, 2017.
[22] KIM J H, JANG W D, SIM J Y, et al. Optimized contrast enhancement for real-time image and video dehazing[J]. Journal of Visual Communication & Image Representation, 2013, 24(3): 410-425.
[23] ZUIDERVELD K. Contrast limited adaptive histogram equalization[J]. Graphics Gems, 1994, 4: 474-485.
[24] HE K, SUN J, TANG X. Single image haze removal using dark channel prior[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2011, 33: 2341-2353.
[25] WEI Y, HAN P L, LIU F, et al. Estimation and removal of backscattered light with nonuniform polarization information in underwater environments[J]. Optics Express, 2022, 30(22): 40208-40220. doi: 10.1364/OE.471337
[26] AGAIAN S S, PANETTA K, GRIGORYAN A M. Transform-based image enhancement algorithms with performance measure[J]. IEEE Transactions on Image Processing, 2001, 10(3): 367-382. doi: 10.1109/83.908502
[27] 卢彦飞, 张涛, 郑健, 等. 基于局部标准差与显著图的模糊图像质量评价方法[J]. 吉林大学学报(工学版), 2016, 46: 1337-1343. LU Yanfei, ZHANG Tao, ZHENG Jian, et al. No-reference blurring image quality assessment based on local standard deviation and saliency map[J]. Journal of Jilin University(Engineering and Technology Edition), 2016, 46: 1337-1343.
[28] 唐菊香, 宋凤丽. 基于图像信息熵优化参数的能见度检测方法[J]. 电子技术与软件工程, 2020(23): 111-113. TANG Juxiang, SONG Fengli. Visibility detection method based on image information entropy optimisation parameters[J]. Electronic Technology & Software Engineering, 2020(23): 111-113.
[29] 陈显, 陈兵. 灰度图像对比度提高方法及其评价[J]. 数字通信世界, 2018(10): 225-226. doi: 10.3969/J.ISSN.1672-7274.2018.10.188 CHEN Xian, CHEN Bing. Contrast enhancement method for grey scale images and its evaluation[J]. Digital Communication World, 2018(10): 225-226. doi: 10.3969/J.ISSN.1672-7274.2018.10.188