Infrared Image Non-uniformity Correction Algorithm Based on Full Convolutional Network
-
摘要: 针对红外成像系统在经过两点校正后,随时间漂移仍然会出现的非均匀性噪声,提出一种基于全卷积深度学习网络的红外图像非均匀性校正算法,使用子网络与主网络相结合的方式进行非均匀性校正。该算法设计了非均匀性等级估计子网络,将含有非均匀性噪声的红外图像输入子网络后,输出非均匀性等级估计图,并和待校正红外图像一并输入校正主网络。子网络生成的非均匀性等级估计图作为一个参数输入校正主网络,避免了网络只针对同一等级非均匀性产生过拟合。经过实验验证,该算法克服了传统的基于场景的算法所产生的边缘模糊问题,对条纹状非均匀性噪声校正效果较好,经过校正后的红外图像清晰度高、细节丰富、边缘清晰、图像质量良好。
-
关键词:
- 非均匀性等级估计子网络 /
- 红外图像 /
- 非均匀性校正 /
- 深度学习
Abstract: The infrared imaging system will still exhibit non-uniform noise after two-point correction. An infrared image non-uniformity correction algorithm based on a fully convolutional deep learning network was proposed in response to this problem. This algorithm combines the subnetwork and main network for non-uniformity correction. The network contains a nonuniformity-level estimation subnetwork. After inputting the infrared image with non-uniform noise into the non-uniformity level estimation subnetwork, the outputted non-uniformity level estimation map is input into the main network together with the original noise infrared image. The non-uniformity level estimated map generated by the subnetwork prevents the network from overfitting only for the non-uniformity of the same grade. After experimental verification, the algorithm overcomes the problem of edge blur generated by the scene-based algorithm. The algorithm will not appear blurred, the images have high definition and rich details, and the quality of images is good. -
0. 引言
近年来,随着红外制冷探测器朝尺寸小、重量轻、功耗低、成本低即低SWaP-C(size, weight and power, cost)方向快速发展,集成小型制冷机的该类中波红外探测器国外已大量应用于武器热瞄镜、便携式手持热像仪、小型无人机、无人车、遥控狙击手和遥控武器站、导弹导引头等空间受限的红外系统[1]。针对该类小型探测器设计一款仅由4片透镜组成的尺寸小、重量轻、成本低的中波红外连续变焦光学系统进而生产结构紧凑、低功耗和低成本的红外变焦热像仪将在手持观瞄具、边防监视系统、小型无人系统等平台得到广泛应用。
目前采用四片式的红外变焦光学多为双视场变焦系统。文献[2]为四片式非制冷双视场,文献[3]为四片制冷型长波双视场系统,文献[4]为4片制冷型中波双视场系统,目前仅采用4片透镜实现连续变焦的制冷型红外光学系统未见报道。
从变焦理论分析,一般常用的机械补偿变焦系统由典型的前固定组、变倍组、补偿组、后固定组4组透镜组成。变倍组一般是负透镜,而补偿组可以是正透镜组也可以是负透镜组,前者为正组补偿系统,后者称为负组补偿系统。从光学系统像差校正难易程度、减少透镜数量、降低光学透镜成本考虑,本文采用无后固定组正组补偿变焦系统,即变焦部分由会聚目标光线的前固定组正透镜、变倍负透镜、补偿正透镜构成。为压缩前固定组物镜口径并满足100%冷屏效率,系统采用二次成像方案,利用单片正光焦度透镜将一次像再次中继成像到探测器焦平面。为压缩轴向尺寸,采用二片平面反射镜将光路U型折转,最终实现仅由4片透镜构成的轻小型中波红外连续变焦光学系统。
1. 连续变焦理论计算模型
机械补偿连续变焦光学系统参数求解就是确定变焦系统在满足像面稳定和焦距在一定范围内连续变化的条件下系统中各组元的焦距、间隔、位移量等参数。通过建立数学模型能方便地计算和分析变焦过程、确定变焦系统高斯光学参数[5]。二组元正组补偿连续变焦系统运动方式如图 1所示。
变焦系统由于只有运动组才产生像面位移,只需抽出变倍组、补偿组加以分析。
因变倍组f2′的移动,引起整个运动组分的像面移动为m32(1-m22)dq,因补偿组f3′的移动,引起整个运动组分的像面移动为(1-m32)dΔ。为达到像面稳定,两个运动组像面移动量的代数和必为零。
$$ m_3^2\left(1-m_2^2\right) \mathrm{d} q+\left(1-m_3^2\right) \mathrm{d} \varDelta=0$$ (1) 而变倍组f2′、补偿组f3′微分移动量dq、dΔ与其倍率变化dm2、dm3之间的关系为:
$$ \mathrm{d} q = \frac{{{f_2}'}}{{m_2^2}}\mathrm{d}{m_2} $$ (2) $$ \mathrm{d} \varDelta=f_3{ }^{\prime} \mathrm{d} m_3 $$ (3) 将(2)、(3)代入(1),经整理得到二组元连续变焦微分方程如下:
$$ \frac{{1 - m_2^2}}{{m_2^2}}{f_2}'\mathrm{d}{m_2} + \frac{{1 - m_3^2}}{{m_3^2}}{f_3}'\mathrm{d}{m_3} = 0 $$ (4) 式(4)是多变量全微分型微分方程,设U(m2, m3)为原函数,则有dU(m2, m3)=0。
其通解为:
$$ U({m_2}, {m_3}) = {f_2}'(\frac{1}{{{m_2}}} + {m_2}) + {f_3}'(\frac{1}{{{m_3}}} + {m_3}) = C $$ (5) 式中:C为常量;设变倍组f2′、补偿组f3′初始状态都处于系统长焦位置,则:
$$ m_{2}=m_{2l}\text{;}m_{3}=m_{3l} $$ (6) 同样有:
$$ {f_2}'(\frac{1}{{{m_{2l}}}} + {m_{2l}}) + {f_3}'(\frac{1}{{{m_{3l}}}} + {m_{3l}}) = C $$ (7) 消去常量C,得到方程的特解:
$$ {f_2}'(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + {f_3}'(\frac{1}{{{m_3}}} - \frac{1}{{{m_{3l}}}} + {m_3} - {m_{3l}}) = 0 $$ (8) 将(8)式整理得到补偿组f3′的倍率m3构成的二次方程:
$$ m_3^2 - b{m_3} + 1 = 0 $$ (9) 其中
$$ b = - \frac{{{f_2}'}}{{{f_3}'}}(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + (\frac{1}{{{m_{3l}}}} + {m_{3l}}) $$ (10) 解得m3的两根为:
$$ {m_{31}} = \frac{{b + \sqrt {{b^2} - 4} }}{2} $$ (11) $$ {m_{32}} = \frac{{b - \sqrt {{b^2} - 4} }}{2} $$ (12) 系统参数求解过程如下:
1)将(2)式积分并整理得到变倍组f2′的倍率m2;
$$ {m_2} = \frac{1}{{\frac{1}{{{m_{2l}}}} + \frac{{{q_2}}}{{{f_2}'}}}} $$ (13) 2)根据求得的m2按照公式(10)求出系数b,再由(11)、(12)式解得补偿组f3′满足运动方程的两个解m31和m32。
3)根据补偿组的两个解求出满足运动方程补偿像面位移的移动量Δ1和Δ2:
$$ \varDelta_{1}=f_{3}′(m_{31}-m_{3l})$$ (14) $$ \varDelta_{2}=f_{3}′(m_{32}-m_{3l}) $$ (15) 4)求出系统的总变倍比
$$ {\varGamma _1} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{31}}}} $$ (16) $$ {\varGamma _2} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{32}}}} $$ (17) 变倍组f2′每移动q2对应着Γ1和Γ2,变倍组和补偿组一起同步运动直到预定的总变倍比为止,到达变倍比的要求的最终状态为系统短焦位置。根据上述连续变焦微分模型,利用(10)~(17)式,可解得各组分光焦度分配及间隔位置关系。
2. 连续变焦系统设计流程
针对复杂的变焦光学系统,建立模型是至关重要的环节,尽管根据(10)~(17)式可推导求解变焦系统高斯参数,但是当光学系统初始参数选择不合适, 会使得系统光焦度分配不合理导致系统像差校正难度大、间隔不合适导致各组件运动中相互碰撞等情况发生。鉴于连续变焦系统的复杂性,系统建模难度大,因此根据连续变焦理论模型编制连续变焦参数计算程序,辅助建立理想光学模型。建模后,设计工作的重点将放在选型以及评价函数的设置和动态修改上,使得设计系统快速收敛[6]。连续变焦光学系统设计流程如图 2所示。
首先,根据设计指标要求需要建立理想光学模型,确定每个组元的参数;再合理选型选材、设定评价函数,进入优化和全局优化;最后根据设计结果评价成像质量。其中函数优化和像质评价环节反复多次迭代,直至达到设计技术指标要求。
3. 连续变焦光学系统设计
3.1 设计指标
中波制冷型连续变焦光学系统采用昆明物理研究所研制生产的长度方向仅119 mm的小型化制冷型中波红外640×512焦平面探测器组件,该探测器具体参数如表 1所示。连续变焦光学系统主要设计指标见表 2。
表 1 探测器参数Table 1. Parameters of detectorDetector HgCdTe Array size 640×512 Pixel dimension/μm 15 NETD/mK ≤22 Spectral response/μm 3.7−4.8 μm Weight ≤380 g 表 2 光学系统设计指标Table 2. Parameters of optical system designWorking waveband/μm 3.7 to 4.8 Zoom 10:01 Field of view 20°×16° to 2.0°×1.6° F# 4 Focal length/mm 27.0~275 Working temperature/℃ -40 to 60 3.2 设计过程
首先,按照系统指标要求,根据连续变焦理论模型编制连续变焦正组补偿参数计算程序求解变焦光学系统各组元光焦度分配及位置间隔关系。在系统初始参数取值上主要考虑以下几点:
1)系统采用二次成像,既可以压缩前固定组直径又可以满足100%冷屏效率。中继组初始倍率取为-1×,前端变焦核心按照系统实际变焦范围进行取值,无需缩放;
2)系统无后固定组,在理论求解中将后固定组倍率取为1×即将后固定组定为无光焦度的虚拟面;
3)为压缩系统变焦过程中变倍组、补偿组位移量,变倍组长焦初始倍率取较大的倍率值,以符合最速变焦理论;
4)考虑光路U型折转,补偿组和中继组之间需较大的空间安置两片平面反射镜,则补偿组取较大的焦距值;
5)系统处于短焦位置时前固定组与变倍组应留出足够的间隔,使两组透镜不至于相碰,初始值设为0.55,补偿组与无光焦度后固定组虚拟面距离初始设为0.55。
利用计算程序,通过反复调整系统初始参数,观察组元间隔、光焦度分配是否合适,最终确定系统初始值为:
$$ \begin{gathered} f_{2}′=-1、f_{3}′=1.62、m_{2l}=-1.45、\\m_{3l}=-1.34、d_{12d}=0.55、d_{34d}=0.55 \end{gathered} $$ 表 3为连续变焦光学系统参数简易计算程序按照上述初始值计算得到的5个视场位置的变焦间隔参数分配结果。
表 3 变焦系统初始间隔参数Table 3. Initial spacing parameters of optical systemFocal length/mm 275 215 150 78.6 25.9 f1/f2 spacing/mm 83.49 80.06 74.56 59.44 18.90 f2/f3 spacing/mm 13.02 24.01 38.57 68.69 125.94 f3/f4 spacing/mm 67.23 59.67 50.61 35.61 18.90 其次,将上述程序计算的变焦光学系统各组元光焦度及间隔位置参数输入ZEMAX光学设计程序,得到系统近轴光学结构,如图 3所示。通过近轴光学结构分析,系统在各位置的焦距值与计算值吻合、总长一致、间隔布局合理、变焦曲线连续,验证程序计算结果正确,可进入下一步选型工作。
再次,考虑各组元材料与调焦镜的选取,由于系统需满足-40℃~+60℃工作环境下成像清晰要求,按照光学系统无热化设计理论,系统的光焦度分配、材料选取、元件间隔都要满足光焦度、消色差、消热差3个方程[7]。
$$ \sum\limits_{i = 1}^j {{h_i}{\varphi _i}} = \varphi $$ (18) $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\theta _i}} = 0 $$ (19) $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\chi _i}} = \sum\limits_{i = 1}^j {{\alpha _i}{L_i}} $$ (20) 式中:hi、ϕi、θi、χi分别为系统各透镜组近轴光线高度、光焦度、色差系数及热差系数;h1为第一透镜近轴光线高;ϕ为系统总光焦度;αi为各透镜间隔镜筒材料的线膨胀系数;Li为各间隔镜筒长度。
该系统采用分步设计技术实现系统连续变焦无热化。首先选取满足上述(18)(19)公式的光焦度和材料分配,实现系统常温状态连续变焦成像清晰及高低温情况下离焦量的线性变化,系统光学材料主要选用硅单晶、锗单晶及硒化锌。其次利用主动补偿技术使系统满足(20)式要求,由于系统透镜数量少,变倍组、补偿组采用凸轮轨道变焦,中继组的轴向移动会产生变倍效果即系统视场焦距发生变化,因此采用前固定组轴向移动来进行主动调焦消热。通过上述无热化设计方法,光学系统在-40℃~+60℃温度范围内保持其性能基本不变。
最后,设置多重结构,合理设计像差评价函数,不断调整优化。为提升连续变焦系统各视场成像清晰度要求,通过设置多个高次非球面和二元衍射面,以提供更多的自由度,有利于球差、色差、像散等各类像差的校正。
3.3 设计结果
轻小型中波红外连续变焦光学系统最终设计结果如图 4所示,从上到下依次为长焦275 mm、中焦100 mm、短焦27 mm的系统图。
系统前固定组采用硅单晶材料用于会聚目标景物光线、压缩变倍组透镜尺寸;变倍组采用锗单晶材料,利用锗单晶高折射率、高色散的特性实现大倍率的变焦;补偿组采用硒化锌材料主要利用其较低的温度折射率系数使其在高低温工作环境只产生较小的离焦量以便实时补偿;中继组选择具有低的温度折射率系数、低色散的硅单晶材料平衡前端变焦核的残留色差。
整个光学系统由4片透镜、两个平面反射镜组成,其中最大透镜为第一透镜其加工直径为71 mm,平面反射镜U型折转后光学系统轴向尺寸长度172 mm,横向尺寸宽度108 mm,光学零件总重量为64 g。该系统光学透镜数量少、重量轻,光路紧凑体积小,冷屏效率100%,适配小型化制冷探测器符合连续变焦光学系统轻小型设计理念。
3.4 二元衍射面分析
由于二元衍射面在消热差、消色差方面的优异特性,系统采用了两个二元衍射面用于减少透镜数量、简化系统设计、提高系统成像质量。
在补偿组透镜硒化锌材料上引入的二元衍射面参数为Norm Radius=15 mm,A1=-33.416,A2=4.580。经计算得到二元衍射面环带深度随透镜径向的变化如图 5所示。硒化锌二元面环带数为4,最大环带深度2.918 μm,最小环带间隔宽度为2.06 mm。该面型环带间隔宽、加工环带数量少,易于单点金刚石车削加工。
在中继组硅单晶材料引入的二元衍射面参数为Norm Radius=16 mm,A1=-60.402,A2=-1.254。其环带深度随透镜径向的变化如图 6所示。硅透镜二元面环带数为9,最大环带深度1.72 μm,最小环带间隔宽度为0.86 mm。该透镜由于材料硬、环带多相对加工难度大,目前昆明物理所光学中心采用单点金刚石车削加工工艺,能制造出满足指标要求的硅基底二元光学元件。
3.5 凸轮曲线计算
本文应用动态光学理论[8],根据像移补偿公式计算补偿组运动曲线。由于变焦组和补偿组均为沿光轴的一维移动,稳像方程为:
$$ {\beta _{2m}}{\beta _2}(1 - {\beta _{1m}}{\beta _1}){q_1} + (1 - {\beta _{2m}}{\beta _2}){q_2} = 0 $$ (21) 式中:β1为变倍组初始位置的垂轴放大率;β1m为变倍组运动后的垂轴放大率;β2为补偿组初始位置的垂轴放大率;β2m补偿组运动后的垂轴放大率;q1为变倍组沿光轴位移量;q2为补偿组沿光轴位移量。
式中:
$$ {\beta _{1m}} = \frac{{{\beta _1}{f_1}'}}{{{f_1}' - {\beta _1}{q_1}}} $$ (22) $$ {\beta _{2m}} = \frac{{{\beta _2}{f_2}'}}{{{f_2}' + (1 - {\beta _1}{\beta _{1m}}){\beta _2}{q_2} - {\beta _2}{q_2}}} $$ (23) 由(21)式~(23)式可得出q1和q2的运动关系,即:
$$ Aq_{2}^{2}+Bq_{2}+C=0$$ (24) 式中:
$$ A=(f_{1}′-β_{1}q_{1})β_{2}\text{;} $$ $$ \begin{align} B=&β_{1}β_{2}q_{1}^{2}+[f_{2}′(1-β_{2}^{2})β_{1}-f_{1}′(1-β_{1}^{2})β_{2}]q_{1}- \\&f_{1}′f_{2}′(1-β_{1}^{2}) \end{align}$$ $$ C=β_{2}^{2}f_{2}′[β_{1}q_{1}-f_{1}′(1-β_{1}^{2})]q_{1} $$ 得到:
$$ {q_2} = \frac{{ - B \pm \sqrt {{B^2} - 4AC} }}{{2A}} $$ (25) 根据上述求解公式计算该变焦系统凸轮曲线如图 7所示。变倍组最大行程为28.6 mm、补偿组最大行程35 mm;补偿组曲线变化平滑,有利于凸轮轨道加工。
4. 系统像质评价
4.1 光学传递函数
系统3个焦距状态下的光学传递函数(MTF)如图 8所示。在3个焦距状态下的系统MTF满足使用要求,光学系统成像质量清晰。
4.2 点列图
系统3个焦距状态下的点列图如图 9所示。各视场RMS(root mean square)均小于一个像素,最大弥散斑RMS半径为12.6 μm,小于像元尺寸;最大弥散斑几何半径为24.5 μm,与系统艾里斑半径20.9 μm相当。系统成像质量良好,满足使用要求。
4.3 畸变
系统畸变情况如图 10所示,在长焦端小视场位置时,最大畸变量为1.2%,在短焦端大视场位置时的最大畸变量为2.4%,该变焦系统畸变对连续成像无明显影响。
4.4 系统高低温成像质量
在高低温工作环境中,系统采用轴向移动前固定组进行主动调焦消热。系统长焦端受环境温度变化,成像质量影响较大,本文主要分析长焦275 mm在高低温下经补偿后的系统成像质量。图 11为系统在高低温下长焦经补偿后的系统调制传递函数。图 12为系统在高低温下长焦经补偿后的系统点列图。从高低温传函图及点列图中看出系统在-40℃~60℃范围内成像质量良好,满足使用要求。
5. 结论
基于小型化制冷中波640×512、像元间距15 μm的焦平面探测器,设计一款具有SWaP-C特征的正组补偿连续变焦光学系统。系统由4片透镜两片平面反射镜组成,F#为4、视场变化范围为20°×16°~2.0°×1.6°,变倍比为10×、U型折叠后系统包络尺寸为172 mm×108 mm、最大物镜口径71 mm、光学零件总重量64 g、零件加工工艺成熟,变焦凸轮曲线平滑,在-40℃~60℃范围内保持较好的成像质量。该轻小型中波红外连续变焦光学系统在导航、搜索、跟踪、警戒、侦察等领域具有广阔的市场前景。
-
表 1 子网络参数
Table 1 Parameters of the subnetwork
Layer Filters Input Output Conv 3×3×32 256×256×1 256×256×32 Conv 3×3×32 256×256×32 256×256×32 Conv 3×3×32 256×256×32 256×256×32 Conv 3×3×1 256×256×32 256×256×1 表 2 校正主网络参数
Table 2 Parameters of main network
Layer Filters Input Output Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×64 256×256×64 256×256×64 Conv 3×3×1 256×256×1 256×256×1 表 3 各算法平均PSNR和SSIM
Table 3 PSNR and SSIM of each algorithm
Algorithm BFTH MHE DLS DMRN Ours PSNR/dB 32.93 33.92 34.39 34.86 35.90 SSIM 0.828 0.858 0.887 0.903 0.928 表 4 各算法平均粗糙度
Table 4 Roughness of each algorithm
Algorithm BFTH MHE DLS DMRN Ours ρ 0.107 0.072 0.064 0.058 0.049 -
[1] ZHOU Huixin, LI Qing, LIU Shangqian, et al. Nonuniformity and its correction principle of infrared focal plane arrays[J]. Laser & Infrared, 2003, 3(6): 46-48. http://www.researchgate.net/publication/293264376_Nonuniformity_and_its_correction_principle_of_infrared_focal_plane_arrays
[2] Scribner D A, Sarkady K A, Kruer M R, et al. Adaptive nonuniformity correction for IR focal-plane arrays using neural networks[C]// International Society for Optics and Photonics, 1991: 100-109.
[3] ZUO C, CHEN Q, GU G, et al. New temporal high-pass filter nonuniformity correction based on bilateral filter[J]. Optical Review, 2011, 18(2): 197-202. DOI: 10.1007/s10043-011-0042-y
[4] QIAN W, CHEN Q, GU G. Space low-pass and temporal high-pass nonuniformity correction algorithm[J]. Optical Review, 2010, 17(1): 24-29. DOI: 10.1007/s10043-010-0005-8
[5] Harris J G, Chiang Y M. Nonuniformity correction using the constant-statistics constraint: analog and digital implementations[C]// Proceedings of SPIE - The International Society for Optical Engineering, 1997, 3061: 895-905.
[6] KUANG Xiaodong, SUI Xiubao, CHEN Qian, et al. Single infrared image stripe noise removal using deep convolutional networks[J]. IEEE Photonics Journal, 2017, 9(4): 1-13. http://www.onacademic.com/detail/journal_1000039958065210_b903.html
[7] 赵春晖, 刘振龙. 改进的红外图像神经网络非均匀性校正算法[J]. 红外与激光工程, 2013, 42(4): 1079-1083. DOI: 10.3969/j.issn.1007-2276.2013.04.044 ZHAO Chunhui, LIU Zhenlong. Improved infrared image neural network non-uniformity correction algorithm[J]. Infrared and Laser Engineering, 2013, 42(4): 1079-1083. DOI: 10.3969/j.issn.1007-2276.2013.04.044
[8] 张龙, 董峰, 傅雨田. 基于神经网络的红外图像非均匀性校正[J]. 红外技术, 2018, 40(2): 164-169. http://hwjs.nvir.cn/article/id/hwjs201802011 ZHANG Long, DONG Feng, FU Yutian. Non-uniformity correction for infrared image using neural networks[J]. Infrared Technology, 2018, 40(2): 164-169. http://hwjs.nvir.cn/article/id/hwjs201802011
[9] MOU Xingang, LU Junjie, ZHOU Xiao, et al. Single frame infrared image adaptive correction algorithm based on residual network[C]//The 11th International Symposium on Photonics and Optoelectronics(SOPO), 2018: 17-23.
[10] 牟新刚, 陆俊杰, 周晓. 基于残差编解码网络的红外图像自适应校正算法[J]. 红外技术, 2020, 42(9): 833-839. http://hwjs.nvir.cn/article/id/hwjs202009004 MOU Xingang, LU Junjie, ZHOU Xiao. Adaptive correction algorithm of infrared image based on encoding and decoding residual network[J]. Infrared Technology, 2020, 42(9): 833-839. http://hwjs.nvir.cn/article/id/hwjs202009004
[11] HE Zewei, CAO Yanpeng, DONG Yafei, et al. Single-image-based nonuniformity correction of uncooled long-wave infrared detectors: a deep-learning approach[J]. Applied Optics, 2018, 57(18): 155-164 DOI: 10.1364/AO.57.00D155
[12] ZHANG Kai, ZUO Wangmeng, ZHANG Lei. Ffdnet: Toward a fast and flexible solution for CNN based image denoising[C]//IEEE Transactions on Image Processing, 2017: 4608-4622.
[13] GUO Shi, YAN Zifei, ZHANG Kai, et al. Toward convolutional blind denoising of real photographs[C]//IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 2019: 1712-1722.
[14] Ioffe S, Szegedy C. Batch normalization: Accelerating deep network training by reducing internal covariate shift[C]//International Conference on Machine Learning, 2015: 448-456.
[15] QIAN W, CHEN Q, GU G. Space low-pass and temporal high-pass nonuniformity correction algorithm[J]. Optical Review, 2010, 17(1): 24-29. DOI: 10.1007/s10043-010-0005-8
[16] CHANG Y, YAN L, LIU L, et al. Infrared Aerothermal nonuniform correction via deep multiscale residual network[J]. IEEE Geoscience and Remote Sensing Letters, 2019, 16(7): 1120-1124. DOI: 10.1109/LGRS.2019.2893519