-
摘要:
空间中子是影响航天器和航天员安全的重要辐射要素之一。优化中子探测器，提高测量效率，提升反演精度是中子测量的难点。中国空间站将搭载一种基于新型中子探测材料Cs2LiYCl6：Ce（CLYC）闪烁体的中子探测器，该探测器具有同时测量热中子和快中子，以及探测效率高等特点。针对该新型探测器的中子能谱反演，分析了不同能量中子在该探测器中的响应特点，分析了中子反演常用的概率迭代法和非负最小二乘（NNLS）法的优缺点，考虑到这2种方法在CLYC探测器反演应用中的不足，提出了基于增广矩阵的非负最小二乘（AM-NNLS）法。数值实验结果表明：AM-NNLS法具有反演运算效率高和反演相对误差小的特点，验证了所提方法的有效性。
Abstract:Neutron is one of the important radiation factors that affect the safety of spacecraft and astronauts. It is difficult for neutron measurement to optimize the neutron detector, improve the measurement efficiency, and improve the inversion accuracy. A neutron detector based on a new type of neutron detection material Cs2LiYCl6:Ce (CLYC) scintillator will be installed on the space station in China. This detector has the characteristics of measuring thermal neutrons and fast neutrons simultaneously, and has high detection efficiency, etc. For the neutron spectrum inversion of this detector, the response characteristics of different energy neutrons in the detetor are analyzed, and the advantages and disadvantages of the probabilistic iterative method and Non-Negative Least Square (NNLS) method commonly used in neutron spectrum inversion are analyzed, considering the disadvantages of these two methods in the inversion application of CLYC detector, a Non-Negative Least Square method based on the Augmented Matrix (AM-NNLS) is proposed. The numerical experiment results show that the AM-NNLS method has the characteristics of high inversion operation efficiency and small inversion relative error, which verifies the effectiveness of the method.
-
雾、霾的形成是由于大气中含有大量的悬浮介质导致光线产生散射或吸收作用，使得视觉能见度降低，如细小水滴、气溶胶、沙尘等颗粒。这种天气现象会使得一些图像及视频采集设备呈现低对比度、色彩偏移、细节丢失等问题，对计算机视觉领域的相关任务造成严重的影响。因此，雾天场景的清晰化复原具有重要的实际意义和研究价值。
目前去雾领域大致可分为图像增强、图像复原和神经网络这3类方法[1]。基于图像增强[2-4]是较为早期的方法，此类方法虽能增强或削弱图像中需要与不需要的信息，但未考虑图像退化的本质，使得最终的去雾效果无法达到预期要求，代表性方法有直方图均衡化[2]、同态滤波[3]和Retinex理论[4]等。
利用图像复原[5-8]是去雾领域中较为传统的一种方法，该系列方法主要是依据大气散射模型进行相应的去雾处理，从图像的退化本质出发，通过对先验知识进行假设来达到估计参数的目的，从而完成对雾图的清晰化处理。文献[5]提出暗通道先验方法，通过对大量雾图研究发现大多数雾图的暗通道值趋于0，间接利用消元的方法对大气散射模型进行化简，以得出参数，但该方法并未考虑到含有天空区域的雾天场景，由于天空区域的暗通道值并不为0，导致暗通道先验失效，使得含天空区域雾图复原结果颜色失真严重；文献[6]提出一种颜色衰减先验方法，利用图像的亮度与饱和度信息建立一种与场景深度相关的线性关系，进而得出无雾图像，但其大气散射系数选取较为单一，透射率估计过大，容易产生去雾不彻底的现象；文献[7]设计了一种分段函数，利用其代替最小值滤波操作，避免了复原结果出现光晕效应，但该方法在一定程度上存在去雾不彻底的问题；文献[8]提出一种基于伽马校正先验的去雾方法，该方法从原始雾图和均匀虚拟转换中提取深度比，再采用全局策略和视觉指示器来恢复场景反照率，但该方法存在估计误差，会导致复原结果出现光晕伪影。
以上方法均是基于传统大气散射模型对雾图进行复原，虽然能在一定程度上达到去雾的效果，但由于传统的大气散射模型未考虑到场景入射光照的影响，因此，在复原过程中，对光照不均匀区域产生的阴影部分并不能得到一个很好的复原效果。
此外，近年来由于神经网络在一些任务上取得了不俗的进展，相关的研究也逐渐运用到了去雾领域上[9-13]。该类方法大致可分为端到端与非端到端2类。文献[9]提出一体化多合一去雾网络，该网络将大气散射模型中的多个未知参数合并为一个进行训练，从而得到无雾图像；文献[10]建立了多个不同尺度的网络模型，利用这种多尺度网络学习雾图与透射率之间所存在的映射关系；文献[11]提出一种仿射双边网络，通过多引导双边学习高清图像的方法，仿射图像的高质量特征，以此达到去雾的目的；文献[12]通过直接学习原图与雾图之间的残差，提出一种门控语境聚合网络；文献[13]提出一种基于渐进特征融合的U-Net网络，可以不依赖大气散射模型直接学习有雾图像和无雾图像之间的高度非线性变换函数。
由于深度学习类方法处理有雾图像时缺少大量真实的场景图进行训练，使用的大多是合成图像，因此，此类方法去雾效果不稳定，在部分真实场景下易产生去雾不彻底的情况。
本文针对传统大气散射模型未考虑场景入射光，导致阴影部分难以复原的问题，提出一种改进型大气散射模型的双光幕边界约束的雾天图像复原算法，并引入均值不等关系对大气光幕进行上下边界约束，使得复原结果适应多种不同场景条件。实验结果表明，本文算法不仅去雾彻底，且对光照不均匀区域的阴影部分也能实现较好的复原效果。
1. 相关背景
1.1 大气散射模型
在数字图像处理及计算机视觉领域中，常用大气散射模型来阐述雾天图像的成像机理，即
I(x,y)=J(x,y)t(x,y)+A(x,y)[1−t(x,y)]=A(x,y)ρ(x,y)t(x,y)+A(x,y)[1−t(x,y)] (1) 式中：$ (x,y) $为像素坐标；$ {{\boldsymbol{I}}}(x,y) $为有雾图像；$ {{\boldsymbol{J}}}(x,y) $为无雾图像；$ {{\boldsymbol{A}}}(x,y) $为大气光值；$ {{\boldsymbol{\rho }}}(x,y) $为场景反照率；$ {{\boldsymbol{t}}}(x,y) $为透射率，当场景中的大气悬浮粒子分布较为均匀时，透射率$ {{\boldsymbol{t}}}(x,y) $可表示为
t(x,y)=exp(−βd(x,y)) (2) 式中：$ \beta $为大气散射系数；$ {{\boldsymbol{d}}}(x,y) $为景深。
模型中$ {{\boldsymbol{A}}}(x,y) $与$ {{\boldsymbol{t}}}(x,y) $为未知参量，$ {{\boldsymbol{J}}}(x,y) $或$ {{\boldsymbol{\rho}} }(x,y) $为待估计参量。在大气散射模型中$ {{\boldsymbol{J}}}(x,y){{\boldsymbol{t}}}(x,y) $为直接传输项，表示大气中悬浮的散射介质粒子经场景辐射后的衰减过程。$ {{\boldsymbol{A}}}(x,y)[1 - {{\boldsymbol{t}}}(x,y)] $为大气光幕，表示大气中的散射光，定义为$ {{\boldsymbol{V}}}(x,y) $，且满足：
{V(x,y)⩽W(x,y)=min[Ic(x,y)]V(x,y)>0 (3) 式中：$ {{\boldsymbol{W}}}(x,y) $为有雾图像的最小通道；$ c\in\{{\mathrm{red}}, {\mathrm{green}}, {\mathrm{blue}}\} $代表RGB三通道[14]。
1.2 大气散射模型的不足
根据传统的大气散射模型对直接传输项与大气光幕的定义，图像的退化原因一般具有3个因素，分别为雾霾对照射目标过程的光源的光谱散射引起的成像质量的原始退化、目标反射光到成像系统在传输路径上的光谱的散射衰减退化、直接进入成像系统的雾霾散射光噪声引起的图像退化[15]。考虑到这些因素，模型中将场景中所有的入射光照定义为大气光这一全局常量其实并不合理，当场景光照不均匀时，大气光幕项中由于所描述的光源未经衰减，所以此时该项传播路径中的粒子大多仍暴露在大气中，可认为全局常量，满足大气光的定义，但在直接传输项中的光源信息应当包括各种反射光与折射光，并非常量。因此，传统大气散射模型中将直接传输项中的光源信息设为恒定的大气光是不合适的。
因此，本文算法拟从场景入射光入手改进传统大气散射模型对直接传输项中光源信息的设定，使得改进后的大气散射模型不仅可以复原一般雾霾场景，针对光源照射不均匀的阴暗区域也可得到较好的复原效果。
2. 本文算法
本文通过引入场景入射光对大气散射模型进行改进，使得该模型在处理光照不均匀图像时，不仅可以较大程度的复原有雾图像的细节信息，同时可以提升图像阴暗区域的可见度。并将雾图分为雾度层和纹理层2部分，避免额外信息对复原结果产生影响。针对未知参数的估计，本文结合高斯衰减函数与均值不等关系设定不同的上下边界对大气光幕进行约束与估计，并利用雾图中通道改进局部大气光的获取，本文算法流程如图1所示。
2.1 改进型大气散射模型
文献[16]通过增加噪声光对传统雾天成像的大气散射模型进行改进，但仅考虑了雾霾场景下的人造光源信息，并未考虑到太阳光由于反射或折射造成的场景光照不均匀区域的信息复原。为完善这类问题，本文算法通过分析成像设备所接收光源的主要途径，结合Retinex理论引入场景入射光进一步改进大气散射模型。
Retinex理论的思想是认为人眼视觉感知到的色彩、亮度主要来自于一般场景中的物体反射，而非受光照影响[4]。根据该理论，一幅图像可表示为
S(x,y)=R(x,y)L(x,y) (4) 式中：$ {{\boldsymbol{S}}}(x,y) $为原图像；$ {{\boldsymbol{R}}}(x,y) $为反射图像；$ {{\boldsymbol{L}}}(x,y) $为场景入射光，入射光的大小决定了一幅图像其像素所能达到的动态范围，一般用原图像与环绕函数的卷积对其进行估计，表达式为
L(x,y)=S(x,y)⊗F (5) 式中：$ \otimes $表示卷积操作；$ F $为中心环绕函数。
而Retinex理论中对场景入射光的描述恰好满足直接传输项中对光源信息的要求，根据这一推论，本文对传统大气散射模型进行改进。将直接传输项中的大气光$ {{\boldsymbol{A}}}(x,y) $改为场景入射光$ {{\boldsymbol{L}}}(x,y) $，因此，式(1)可改写为
I(x,y)=L(x,y)ρ(x,y)t(x,y)+A(x,y)[1−t(x,y)]=L(x,y)ρ(x,y)[1−1A(x,y)V(x,y)]+V(x,y) (6) 由于Retinex理论在复原有雾图像时考虑到了场景入射光对复原效果的影响，基于这一点，将式(6)改为Retinex的形式，即
Snew(x,y)=Rnew(x,y)Lnew(x,y) (7) 式(7)中各个参数的表达式为
{Snew(x,y)=I(x,y)−V(x,y)Rnew(x,y)=[1−1A(x,y)V(x,y)]ρ(x,y)Lnew(x,y)=L(x,y) (8) 式中：$ {{{\boldsymbol{S}}}_{{\text{new}}}}(x,y) $为新模型下的原图像；$ {{{\boldsymbol{R}}}_{{\text{new}}}}(x,y) $为新模型下的反射图像；$ {{{\boldsymbol{L}}}_{{\text{new}}}}(x,y) $为场景入射光。改进前后的大气散射模型对比，如图2所示。
图 2 大气散射模型物理退化过程Figure 2. Physical degradation process of atmospheric scattering model此外，为避免在图像复原过程中，可能存在的纹理、噪声等信息对复原结果产生影响，使主观效果不佳。本文将有雾图像$ {{\boldsymbol{I}}}(x,y) $分为雾度层$ {{\boldsymbol{C}}}(x,y) $和纹理层$ {{\boldsymbol{T}}}(x,y) $，雾度层为有雾图像的雾霾与背景等信息，纹理层表示有雾图像的纹理细节、噪声等信息，表达式为
I(x,y)=C(x,y)+T(x,y) (9) 在图像复原过程中，先去除有雾图像纹理层，再按改进型大气散射模型估计各参数，得到初级无雾图像，最终将提取出的纹理层图像与初级无雾图像进行融合得到最终的复原结果。这样即可避免纹理层信息对复原结果造成的干扰，又能使得复原结果更加贴近真实场景，具体参数求取在后续章节进行。
图3为分别基于改进前后2种大气散射模型下的复原结果对比。可以看出，基于改进后的大气散射模型下的复原图像较传统大气散射模型更能复原出有雾图像在光照不均匀处产生的阴影区域，且细节信息更为丰富。
图 3 2种大气散射模型的复原效果对比Figure 3. Comparison of restoration effects of two atmospheric scattering models2.2 双光幕边界约束下的大气光幕估计
考虑到大气光幕中同时含有待估计参数透射率与大气光，因此，本文算法根据式(3)中大气散射模型对大气光幕不等关系的定义，设定上下边界对其进行约束，以求取较为准确的大气光幕。
2.2.1 大气光幕上边界
根据式(3)可以看出，大气光幕总是小于或等于有雾图像的最小通道且大于0，因此，选定雾图的最小通道$ {{\boldsymbol{W}}}(x,y) $为大气光幕的上边界。此外，随着场景深度的叠加，雾浓度升高，视觉可见度会降低，使得图像在远景区域的纹理信息不可见。为了使得所求大气光幕更能准确地表示雾浓度分布的情况，本文对有雾图像的上边界，即最小通道图采用引导滤波[17]处理，以消除纹理与部分无用边缘信息使其更为平滑，以逼近更为真实的大气光幕。
V1(x,y)=FguideW(x,y)) (10) 式中：$ {\boldsymbol{V}}_{1}(x,y) $表示大气光幕上边界；$F_ {\rm{guide}} $表示引导滤波操作。
2.2.2 大气光幕下边界
大气散射模型中大气光幕$ {{\boldsymbol{V}}}(x,y) $被定义为
V(x,y)=A(x,y)[1−exp(−βd(x,y))] (11) 假设大气光$ {{\boldsymbol{A}}}(x,y) $已知，由式（11）中可以看出透射率$ {{\boldsymbol{t}}}(x,y) $与场景深度$ {{\boldsymbol{d}}}(x,y) $呈衰减关系，而大气光幕的大小与透射率或景深的变化有关。因此，本文算法从景深与透射率出发求取大气光幕下边界。
根据文献[6]中的颜色衰减先验可知，随着场景深度的叠加，亮度升高，饱和度降低；即景深$ {{\boldsymbol{d}}}(x,y) $与亮度$ {{\boldsymbol{v}}}(x,y) $呈递增关系，与饱和度$ {{\boldsymbol{s}}}(x,y) $呈衰减关系。引入亮度和饱和度的定义式为
{v(x,y)=max[Ic(x,y)]s(x,y)=1−[min[Ic(x,y)]/max[Ic(x,y)]] (12) 对饱和度进行反转，得到反饱和度$ {{{\boldsymbol{s}}}'}(x,y) $为
s′(x,y)=1−s(x,y)=min[Ic(x,y)]/max[Ic(x,y)] (13) 根据式(12)和式(13)可以看出饱和度与反饱和度呈衰减关系，结合景深与饱和度呈衰减关系，可以得出景深与反饱和度呈递增关系。
结合透射率$ {{\boldsymbol{t}}}(x,y) $与场景深度$ {{\boldsymbol{d}}}(x,y) $呈衰减关系，场景深度$ {{\boldsymbol{d}}}(x,y) $与亮度$ {{\boldsymbol{v}}}(x,y) $、反饱和度$ {{{\boldsymbol{s}}}'}(x,y) $呈递增关系，可得出透射率$ {{\boldsymbol{t}}}(x,y) $与亮度$ {{\boldsymbol{v}}}(x,y) $、反饱和度$ {{{\boldsymbol{s}}}'}(x,y) $呈衰减关系。本文通过设定衰减因子$ p $来表示这种衰减关系。
{t1(x,y)=p1⋅[v(x,y)]t2(x,y)=p2⋅[s′(x,y)] (14) 以景深分布较为均匀的雾图为例，分析其亮度、饱和度、反饱和度的特征分布及其灰度性能曲线，如图4所示。根据图4(e)的亮度、饱和度、反饱和度随场景深度变化的灰度分布形式可以看出，随着场景深度的叠加，当亮度递增时，饱和度降低，反饱和度增高，进一步印证了本文假设的可行性。
鉴于高斯函数具有较好的衰减特性，本文设定一种均值为1，标准差为$ \sigma $，且幅值进行归一化处理后的高斯衰减函数来拟合前述的衰减关系，由此得到式(14)的近似估计为
{t1(x,y)=exp{−[1−v(x,y)]2/σ12}t2(x,y)=exp{−[1−s′(x,y)]2/σ22} (15) 考虑到有雾图像在远景处亮度值与反饱和度值偏高；而在近景处则出现亮度值与反饱和值偏低的现象。根据上述分析，场景深度与亮度和反饱和度呈递增关系，将二者进行乘积，场景深度与亮度和反饱和度的乘积形式亦满足递增关系，即
d(x,y)∝v(x,y)∝s′(x,y)∝v(x,y)⋅s′(x,y) (16) 联合式(12)和式(13)，可得出景深与雾图的最大、最小通道均呈递增关系为
d(x,y)∝max[Ic(x,y)]∝min[Ic(x,y)] (17) 结合所得这一结论，本文以有雾图像的灰度均值$ {\text{mean}}[{{{\boldsymbol{I}}}^c}(x,y)] $为阈值将图像划分为远景和近景区域，并设定$ a $和$ {b} $分别作为2个区域的特征值，其中$ a, b \in(0,1) $。$ a $和$ b $可表示为
{a=max[Ic(x,y)]−mean[Ic(x,y)]2b=mean[Ic(x,y)]−min[Ic(x,y)]2 (18) 获取特征值后，本文引入均值不等关系对式(15)所估计的透射率进行逼近，这样可以使算法在不涉及大规模优化迭代计算的前提下，最大程度上求取标准差区域内的最优解，使得算法具有自适应的功能。首先，根据数学均值不等关系有
2aba+b⩽√ab⩽a+b2 (19) 结合式(15)，对$ \sigma $进行定义：
{σ1=2aba+bσ2=a+b2 (20) 假设大气光已知，利用均值不等式可以得到大气光幕下边界的粗略估计$ {\boldsymbol{V}}'(x,y） $为
V′(x,y)=A(x,y)[1−√t1(x,y)⋅t2(x,y)] (21) 由于所估计的粗略大气光幕是采用边界逼近的方法，可能会损失部分细节信息，为获取更为准确的大气光幕信息，利用梯度约束[18]结合梯度算子求取边缘代价对粗略大气光幕下边界进行补偿：
E(x,y)=∑x,y[V1(x,y)−V′(x,y）]2+η{α[∂V1(x,y)/∂x]2+β[∂V1(x,y)/∂y]2} (22) 式中：$ {{\boldsymbol{E}}}(x,y) $为边缘代价；η为惩罚项系数，便于最小化全局的最优代价函数；$ {\boldsymbol{V}}_{1}(x,y)-{\boldsymbol{V}}'(x,y） $表示细化初始大气光幕上边界与粗略大气光幕下边界的纹理信息之差；等式的右项为惩罚项，用于标记图像的边缘信息；$ \alpha $、$ \beta $为边缘系数，本文采用对数压缩处理，使其具有自适应的能力。
{α=[∂ln|V1(x,y)|/∂x+ε]−1β=[∂ln|V1(x,y)|/∂y+ε]−1 (23) 为了防止分母为0，取$ \varepsilon $=0.1。由此可得到修正后的大气光幕下边界$ {{{\boldsymbol{V}}}_2}(x,y) $的表达式为
V2(x,y)=V′(x,y）+E(x,y) (24) 2.2.3 大气光幕估计
取上下边界之和的均值作为最终大气光幕$ {\boldsymbol{V}}(x, y) $，表达式为
V(x,y)=V1(x,y)+V2(x,y)2 (25) 图5为大气光幕估计和复原结果，可以看出所估计出的大气光幕不仅能表示雾浓度的分布特征，且细节保证良好，复原结果色彩鲜明。
图 5 大气光幕估计和复原结果示意图Figure 5. Schematic diagram of atmospheric light curtain estimation and restoration results2.3 场景入射光估计
根据2.1节分析可知，传统大气散射模型中对大气光的定位是不合理的，并没有考虑到场景反射光、背景折射光等光照因素的影响，从而导致雾图的阴暗区域无法得到很好的复原效果。考虑到Retinex理论具有较强的动态压缩能力，能够以较大程度复原图像阴影区域的细节信息，因此，结合式(8)中的入射光图像$ {{\boldsymbol{L}}}(x,y) $，新原图像$ {{{\boldsymbol{S}}}_{{\text{new}}}}(x,y) $，利用Retinex理论引入场景入射光，并得到其粗略估计$ {\boldsymbol{L}}'(x,y) $，表达式为
L′(x,y)=(I(x,y)−V(x,y))⊗F (26) 式中：$ F $选用标准差为0.5，卷积窗口尺寸$ h $由雾图$ {{\boldsymbol{I}}}(x,y) $中尺寸的最小值$ {m} $决定的高斯函数作为中心环绕函数，$ h $的表达式为
h=floor(m/50) (27) 式中：$ { {\mathrm{floor}}(\cdot) } $表示向负无穷方向取整。
为避免所得场景入射光在处理光照均匀场景时，对目标场景出现过度增强的影响；同时在处理光照不均匀场景时，为保证所得复原场景的真实性，保留些许光照不均匀区域的阴暗程度，对所求场景入射光进行限制：
L1(x,y)=max(L′(x,y),L0) (28) 式中：$ {{\boldsymbol{L}}_1}(x,y) $为经过限制约束处理后的初始场景入射光；$ {{{L}}_0} $为初始场景入射光的下限值（本文中$ {L_0}{\text{ = 0}}{\text{.1}} $）。
由于所求场景入射光是由Retinex理论得来，但基于图像增强的去雾方法并未考虑到图像退化的本质，可能会导致最终复原结果出现颜色失真的现象。因此，为了能基于图像退化的本质得到较好的去雾效果，本文引入亮通道先验[19]对所求场景入射光进行补偿，根据亮通道先验所述，在有雾图像的RGB三通道中至少存在一个通道具有较高的像素值，这说明在部分局部区域内这些像素具有很强的场景光反射能力和较高的亮度。所以，亮通道图可以对雾图中受光照影响产生阴影较少的区域起到很好的修补作用，如图像中的远景区域。使得所求场景入射光能覆盖到雾图的各个区域。通常情况下，亮通道被定义为
L2(x,y)=maxy∈Ω(x,y)[max(Ic(x,y)] (29) 式中：$ {{{\boldsymbol{L}}}_2}(x,y) $为亮通道值；$ {{\boldsymbol{\varOmega}} } $为一个以$ (x, y) $为中心的局部区域。
将所求粗略场景入射光与有雾图像的亮通道按照一定比例进行融合，从而达到对场景入射光进行补偿的目的，表达式为
L(x,y)=λ1⋅L1(x,y)+λ2⋅L2(x,y) (30) 式中：$ {\lambda _1} $和$ {\lambda _2} $为权重参数。由于$ {{{\boldsymbol{L}}}_2}(x,y) $在一定程度上也反映了场景光照的能力，也可作为场景入射光的一种近似估计，与$ {{{\boldsymbol{L}}}_1}(x,y) $在一定程度上具有相似性。考虑到亮通道图在一定程度上反应了图像的亮度信息，因此，本文利用图像存在的亮度特征，对$ {{{\boldsymbol{L}}}_1}(x,y) $和$ {{{\boldsymbol{L}}}_2}(x,y) $的相似性进行度量，从而求取相关系数。在雾天条件下获得的图像，随着雾浓度升高，图像的亮度值也逐渐增大。本文所得大气光幕可反应雾图的雾浓度分布情况，其灰度均值的大小恰可以反映图像的亮度变化，可以将其作为$ {{{\boldsymbol{L}}}_2}(x,y) $的前缀系数。考虑到$ {{{\boldsymbol{L}}}_1}(x,y) $和$ {{{\boldsymbol{L}}}_2}(x,y) $分别是从图像增强与复原两方面对场景入射光进行估计的，二者具有一定的互补性。因此，将二者进行加权融合后，相关系数可表示为
{λ2=mean[V(x,y)]λ1=1−λ2 (31) 本文所求场景入射光如图6所示。根据图中可以看出，虽然亮通道图可以反映图像的整体的光照水平，但并不能将其视为真正的场景入射光。因为场景入射光需要覆盖到图像的各个区域，但亮通道图属于单通道图，无法估算出由于光照不均匀导致反射能力较差处（如阴暗区域）的场景入射光值。但亮通道图作为场景入射光的补偿值是合适的，经过亮通道对粗略场景入射光的补偿后可以看出，补偿后的场景入射光图，不仅对场景的覆盖能力更广，且在一定程度上反应了景深信息，同时避免了由于场景入射光不足导致复原结果颜色失真的问题。
2.4 大气光估计
大气光$ {{\boldsymbol{A}}}(x,y) $是去雾问题中另一重要参数，大气光的选取决定了恢复图像的视觉效果。当大气光估计偏高时，结果偏暗；大气光估计偏低时，结果偏亮。文献[20]从图像的局部信息出发，以雾图的三通道最大值为出发点求取大气光，从而减弱了大多数算法从全局信息进行大气光估计带来的影响。但利用三通道最大值会使得大气光的估计偏高，影响复原效果。
为了解决该问题，并减少有雾图像中含有高亮区域对复原结果产生影响，本文提出一种基于中通道的局部大气光求取方法。参照式(19)的均值不等式，将雾图的最大通道和最小通道分别代入不等式的左右两边，可得到有雾图像的中通道表达式为
{max[Ic(x,y)]=a1+b12min[Ic(x,y)]=2a1b1a1+b1⇒mid[Ic(x,y)]=√a1b1 (32) 为避免高亮噪声干扰，对所求中通道进行形态学闭操作，并利用交叉双边滤波对其进行平滑处理，以此得到最终的大气光：
A(x,y)=FBF(imclose(mid[Ic(x,y)],se)) (33) 式中：$ F_{\mathrm{BF}}(\cdot) $表示交叉双边滤波操作；$ {{\mathrm{imclose}}(\cdot) } $表示形态学闭操作；$s_{\mathrm{e}}$表示结构元素对象。复原效果对比如图7所示。
图 7 大气光图及其复原效果对比Figure 7. Comparison of atmospheric light maps and their restoration effects2.5 图像复原
综合所求大气光幕，场景入射光，大气光，并结合式(7)和式(8)可得到最终的场景反照率为
ρ(x,y)=A(x,y)[I(x,y)−V(x,y)]L(x,y)[A(x,y)−V(x,y)] (34) 最后，将场景反照率限制在[0,1]范围内，可得到无雾图像$ {{\boldsymbol{O}}}(x,y) $，其表达式为
O(x,y)=min[max(ρ(x,y),0),1] (35) 此外，根据式(9)及其描述，本文引入纹理层图像$ {{\boldsymbol{T}}}(x,y) $改善最终复原效果。考虑到形态学开运算可以对目标图像进行去噪处理，因此，去噪后的图像可近似看作为雾度层图像$ {{\boldsymbol{C}}}(x,y) $，利用有雾图像$ {{\boldsymbol{I}}}(x,y) $减去雾度层图像$ {{\boldsymbol{C}}}(x,y) $即可得到纹理层图像$ {{\boldsymbol{T}}}(x,y) $。具体表达式为
C(x,y)=imopen(I(x,y),se) (36) T(x,y)=I(x,y)−C(x,y) (37) 式中：$ { {\mathrm{imopen}}(\cdot) } $表示形态学开运算。因此，最终的雾图复原结果$ {{\boldsymbol{J}}}(x,y) $可表示为
J(x,y)=O(x,y)+T(x,y) (38) 各参量即纹理层添加前后的复原效果对比如图8所示。从图8(c)和图8(d)可以看出，纹理层添加后“飞机机身”较添加前较清晰，添加的前机身有部分灰色杂质，且尾翼处的标志不够清晰，地面处颜色较暗。
图 8 纹理层图像及添加前后复原效果对比Figure 8. Comparison of texture layer images and restoration effects before and after adding算法复原结果及各个未知量估计如图9所示。
3. 实验结果及分析
本文选取了几组不同场景下的真实雾图进行算法验证。并从定性、定量的主客观两方面与一些经典算法进行对比。另外，选取了4幅深度学习算法常用的RESIDE测试集雾图验证对合成图像处理上的可行性。实验环境为MATLAB（R2018a），硬件环境为Intel(R) Core(TM) i7-10510U CPU @ 1.80 GHz，2.30 GHz，运行系统为Windows11。
3.1 去雾质量定性分析
本文选取4组8幅不同场景环境下的真实雾图，分别为光照不均匀图像、颜色变化剧烈图像、远景含天空图像、近景含人物图像。为了验证算法的普适性，在图像复原和深度学习2类方法上，本文各选3种算法进行实验对比，分别与文献[5]算法、文献[6]算法、文献[7]算法、文献[10]算法、文献[12]算法、文献[13]算法的复原结果进行对比。此外，另选取4幅合成雾图与不同算法进行实验对比。
图10为光照不均匀组图像，从图中可以看出文献[5]算法在远景处的复原效果较好，但由于暗通道先验中使用最小滤波操作导致远近景交替处出现光晕效应；文献[6]算法由于透射率估计偏大导致复原结果出现去雾不彻底的问题；文献[7]算法在近景去雾效果较好，在远则存留有稍许残雾，且复原结果整体呈现锐化过度的问题；文献[10]算法出现严重的颜色失真，尤其在近景区域；文献[12]算法与文献[13]算法由于缺乏对真实雾图的训练导致复原结果出现严重的去雾断层现象，特别是文献[13]算法尤为明显。且以上算法在光照不均匀的阴暗区域均无法复原出较为清晰的细节信息，相较于其他算法，本文算法不仅去雾彻底，色彩均衡，且对近景处由于光照不均导致的阴暗区域也能得到较好的复原结果。
图 10 不同算法下光照不均匀图像的复原结果Figure 10. Restoration results of images with uneven light under different algorithms图11为颜色变化剧烈组图像的复原结果，可以看出，文献[5]算法与文献[10]算法此时出现不同程度的颜色失真现象，文献[10]算法的偏色程度相较于文献[5]算法较严重；文献[6]算法与文献[7]算法则出现了不同程度的去雾不彻底现象，文献[6]算法的残雾存留程度相较于文献[7]算法更为严重；文献[12]算法与文献[13]算法同样因为训练不足的原因导致去雾效果不稳定。
图 11 不同算法下颜色变化剧烈图像的复原结果Figure 11. Restoration results of images with drastic color changes under different algorithms图12为远景含天空组图像的复原结果，可以看出文献[5]算法由于对暗通道的估计未考虑到远景天空处，导致透射率严重偏小，使得复原结果的远景天空处出现偏色问题；而文献[6]算法由于是对场景深度进行估计的算法，所以在该组图像中表现较前2组效果较好；文献[7]算法此时在近景区域则出现色失真的问题；文献[10]算法则在远景区域出现了失真的问题；文献[12]算法在该组去雾较为彻底，但色彩偏移较为严重；文献[13]算法则几乎无法去雾。
图 12 不同算法下远景含天空图像的复原结果Figure 12. Restoration results of distant view with sky images under different algorithms图13为近景含人物组图像的复原结果，此时各算法复原效果与图11组图像近景处较为相似。文献[5]算法、文献[7]算法、文献[10]算法、文献[12]算法此时呈现不同程度的颜色失真现象，并伴随些许色彩偏移的问题；文献[6]算法与文献[13]算法则出现去雾不彻底的现象，且该现象在文献[13]算法处理该组图像的人物周围尤其明显。
图 13 不同算法下近景含人物图像的复原结果Figure 13. Restoration results of close shot images with people under different algorithms图14为RESIDE测试集图像的复原结果，可以看出，在处理合成图像时，文献[5]算法、文献[6]算法、文献[7]算法同样出现了上述问题；文献[10]算法、文献[12]算法、文献[13]算法的复原效果较好，更接近真实图像，但文献[12]算法有轻微偏色的问题；而本文算法相较于其他算法，不仅贴合真实图像，且对光照阴影区域也能呈现较好的处理效果。从复原结果与真实场景的对比中可以看出，针对该组图像的第1幅与第4幅图像，图像的阴影区域信息本文算法可以较大程度的对其进行复原。
图 14 不同算法下RESIDE测试集图像的复原结果Figure 14. Restoration results of RESIDE test set images under different algorithms3.2 去雾质量定量分析
仅依照定性的主观评价的判断较为片面，因此，本文从定量的客观方面进行侧面验证。将新增可见边$ e $、平均梯度$r$、饱和像素点占比$ \delta $及运行时间$ t $作为评价指标[21]。其中，可见边增加率与平均梯度越大代表实验指标效果越好，饱和像素点占比与运行时间越小代表实验指标越好。实验对象为主观评价中所选8幅真实雾图的各算法与本文算法复原结果，各指标的数学形态如式(39)～式(41)所示：
e=nr−n0n0 (39) r=exp(1nr∑pi∈Xlog2ri) (40) δ=nsMN (41) 式中：$ n_{{\mathrm{r}}} $和$ n_{0} $分别为算法处理前后雾图和无雾图的可见边数量；$ r_{i} $为$ p_{i} $处去雾前后图像的平均梯度比；$ X $为无雾图像的可见边集合；$ n_{{\mathrm{s}}} $为饱和点像素个数；$ MN $为图像宽与高的尺寸乘积。各算法的定量评价指标如图15所示，图中横轴中色块对应的算法从左至右分别为文献[5]算法、文献[6]算法、文献[7]算法、文献[10]算法、文献[12]算法、文献[13]算法和本文算法。
图 15 各个算法的定量分析指标对比Figure 15. Comparison of quantitative analysis indexes of each algorithms针对测试集图像采用峰值信噪比和结构相似性来测试各算法对合成图像复原效果的优劣程度。两者数值越大表明效果越好。指标对比如表1所示（表中数据为均值）。
表 1 测试集评测指标Table 1. Test set evaluation index通过实验数据指标可以看出，各个算法在客观指标数据下均体现了各自的优势，文献[12]算法在部分图像复原中在新增可见边上存在一定优势；文献[6]算法在运行时间上要优于其他算法；但相较于其他2个指标，本文算法具有一定的优势。在测试集的指标验证方面，本文算法的峰值信噪比优于其他算法，虽然文献[13]算法在结构相似性上取得了较高的指标，但本文算法整体主观效果较好，且优于其他算法。综合以上定性、定量的主客观分析，本文算法具有一定的可行性和应用价值。
4. 结 论
1） 本文算法有效地解决了传统大气散射模型在处理光照不均匀图像时阴暗区域难以复原的问题。
2） 本文算法通过引入场景入射光、细化后的雾图最小通道和亮度、饱和度与透射率的关系结合均值不等关系和特征值估计设定衰减函数等方法，对大气光幕进行估计，并利用亮通道图补偿粗略的场景入射光，并结合中通道理念改进局部大气光。
3） 本文算法可以有效地复原光照不均匀的有雾图像的阴暗区域，并且去雾效果良好，边缘信息保持完整。
-
图 3 热中子以及5 MeV的单能中子在CLYC闪烁体中的沉积能谱
Figure 3. Deposited energy spectrum of thermal neutron and 5 MeV monoenergetic neutron in CLYC scintillator
图 4 CLYC闪烁体中对热中子和快中子(0.025 eV~100 MeV)的响应函数
Figure 4. Response functions of thermal neutron and fast neutron (0.025 eV-100 MeV) in CLYC scintillator
图 7 基于概率迭代法和NNLS法的中子微分谱
Figure 7. Neutron differential spectrum based on probabilisticiterative method and NNLS
图 8 NNLS法与AM-NNLS法的能谱对比
Figure 8. Comparative energy spectrum of non-negative least squares method and AM-NNLS
图 11 测量总计数不同时2种反演方法得到的能谱与真实入射谱之间的相对误差
Figure 11. Relative error between energy spectrum obtained from two inversion methods and actual input spectrum when measured values are different
反应 Q/MeV 能量区间 6Li+n→3H+ɑ +4.78 热中子及快中子 35Cl+n→35S+P +0.615 快中子 
下载: 导出CSV
表 2 概率迭代法与AM-NNLS法反演能谱所需时间对比
Table 2. Time comparison between probabilistic iterative method and AM-NNLS for inversion of energy spectrum
测量总计数 耗时/s 概率迭代法 AM-NNLS 102 138 5 103 143 4 104 156 5 105 150 5
下载: 导出CSV
-
[1] MOUNTFORD P J, TEMPERTON D H.Recommendations of the international commission on radiological protection(ICRP)1990[J].European Journal of Nuclear Medicine, 1992, 19(2):77-79. [2] MAZUR J E, CRAIN W R, LOOPER M D, et al.New measurements of total ionizing dose in the lunar environment[J].Space Weather-The International Journal of Research & Applications, 2011, 9(7):1-12. [3] KOSHIISHI H, MATSUMOTO H, CHISHIKI A, et al.Evaluation of the neutron radiation environment inside the International Space Station based on the bonner ball neutron detector experiment[J].Radiation Measurements, 2007, 42(9):1510-1520. doi: 10.1016/j.radmeas.2007.02.072 [4] HASSLER D M, ZEITLIN C, WIMMER-SCHWEINGRUBER R F, et al.The radiation assessment detector (RAD) investigation[J].Space Science Reviews, 2012, 170(1-4):503-558. doi: 10.1007/s11214-012-9913-1 [5] TANIGUCHI T, UEDA N, NAKAZAWA M, et al.Systematic study on spectral effects in the adjustment calculations using the NEUPAC-83 code[M].Berlin:Springer, 1985:685-691. [6] KOHLER J, EHRESMANN B, MARTIN C, et al.Inversion of neutron/gamma spectra from scintillator measurements[J].Nuclear Instruments and Methods in Physics Research B, 2011, 269:2641-2648. doi: 10.1016/j.nimb.2011.07.021 [7] 王冬, 何彬, 张全虎.用遗传算法求解中子能谱[J].原子能科学技术, 2010, 44(10):1270-1275.WANG D, HE B, ZHANG Q H.Unfolding neutron spectrum using genertic algorithm[J].Atomic Energy Science and Technology, 2010, 44(10):1270-1275(in Chinese). [8] 杨鑫, 李润东, 刘汉刚, 等.基于概率迭代的NDP反演方法[J].计算物理, 2012, 29(6):891-900. doi: 10.3969/j.issn.1001-246X.2012.06.014YANG X, LI R D, LIU H G, et al.An unfolding method of NDP based on probability iteration[J].Chinese Journal of Computational Physics, 2012, 29(6):891-900(in Chinese). doi: 10.3969/j.issn.1001-246X.2012.06.014 [9] D'OLYMPIA N W, NATHAN W.Development of novel neutron and gamma-ray scintillators: Cs2LiYCl6: Ce and CeBr3[D].Boston: University of Massachusetts, 2013: 8-40. [10] 胡圣荣, 戴纳新.病态线性方程组的新解法:增广方程组法[J].华南农业大学学报, 2009, 30(1):119-121.HU S R, DAI N X.A novel method for solving Ⅲ-conditioned liner system:Augmented system method[J].Journal of South China Agricultural University, 2009, 30(1):119-121(in Chinese). [11] 李肖.反符合探测杯在空间粒子探测中的应用[D].北京: 中国科学院, 2015: 5-23.LI X.The application of anti-coincidence detective cup technology in space particle detection[D].Beijing: University of Chinese Academy of Science, 2015: 5-23(in Chinese). [12] WHITNEY C M, SOUNDARA-PANDIAN L, JOHNSON E B, et al.Gamma-neutron imaging system utilizing pulse shape discrimination with CLYC[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2015, 784:346-351. doi: 10.1016/j.nima.2014.09.022 [13] LEE D W, STONEHILL L C, KLIMENKO A, et al.Pulse-shape analysis of Cs2LiYCl6:Ce scintillator for neutron and gamma-ray discrimination[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2012, 664(1):1-5. [14] WANG Q B, TUO X G, DENG C, et al.Characterization of a Cs2LiYCl6:Ce3+ scintillator coupled with two silicon photomultiplier arrays of different sizes[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2019, 942:1-5. [15] D'OLYMPIA N W, CHOWDHURY P, LISTER C J, et al.Pulse-shape analysis of CLYC for thermal neutrons, fast neutrons, and gamma-rays[J].Nuclear Instruments and Methods in Physics Research Section A:Accereators, Spectrometers, Detectors and Associated Equipment, 2013, 714:121-127. [16] BROWN D, CHADWICK M B, CAPOTE R, et al.ENDF/B-Ⅷ.0:The 8th major release of the nuclear reaction data library with CIELO-project cross sections, new stantards and thermal scattering data[J].Nuclear Data Sheets, 2018, 148(2):1-142. [17] AGOSTINELLI S, ALLISON J, AMAKO K, et al.GEANT4-A simulation toolkit[J].Nuclear Instruments and Methods in Physics Research Section A:Accelerators, Spectrometers, Detectors and Associated Equipment, 2003, 506(3):250-303. [18] ZHANG S Y, ZHANG X G, WANG C Q, et al.The geometric factor of high energy protons detector on FY-3 satellite[J].Science China Earth Sciences, 2014, 57(10):2558-2566. doi: 10.1007/s11430-014-4853-0 [19] KHARYTONOV A, BÖHM, E, WIMMER-SCHWEINGRUBER R F.Regularization methods used in error analysis of solar particle spectra measured on SOHO/EPHIN[J].Astronomy & Astrophysics, 2009, 495(2):663-675. [20] GUO J N, BANJAC S, ROSTEL L, et al.Implementation and validation of the GEANT4/AtRIS code to model the radiation environment at Mars[J].Journal of Space Weather and Space Climate, 2019, 9(3):1-30. [21] NOWAK R D, KOLACZYK E D.A statistical multiscale framework for Poisson inverse problems[J].IEEE Transactions on Information Theory, 2000, 46(5):1811-1825. doi: 10.1109/18.857793 [22] WILLIAM H P, BRIAN P F, SAUL A, et al.Numerical recipes:The art of scientific computing[M].Cambridge:Cambridge University Press, 2007:809-816. [23] LAWSON C L, HANSON R J.Solving least squares problems[M].Philadelphia:Society for lndustrial and Applied Mathematics, 1987:1-5. [24] BÖHM E, KHARYTONOV A, WIMMER-SCHWEINGRUBER R F.Solar energetic particle spectra from the SOHO-EPHIN sensor by application of regularization methods[J].Astronomy & Astrophysics, 2007, 473(2):673-682. [25] MORALE J L, NOCEDAL J.Remark on 'algorithm 778:L-BFGS-B:Fortran subroutines for large-scale bound constrained optimization'[J].ACM Transactions on Mathematical Software, 2011, 38(1):1-4. 期刊类型引用(4)
1. GuoHong Shen，DongHui Hou，Yuan Chang，XianGuo Zhang，HuanXin Zhang，Bin Yuan，BinQuan Zhang，Ying Sun. Neutron observations from the energetic particle detector on China's Space Station. Earth and Planetary Physics. 2025(02): 460-466 . 
必应学术2. 白超平，张珅毅，张鑫，孙越强，张帅，王子婷. 基于卷积神经网络的星载粒子鉴别平台及应用. 北京航空航天大学学报. 2025(04): 1313-1323 . 
本站查看3. 夏庆中，蒋勇，夏广新，何铁. 强γ抑制的~6LiF/ZnS(Ag)中子探测器研制及性能测试. 原子能科学技术. 2023(11): 2218-2224 . 百度学术4. 李雅男，李桃生，洪兵，郭浩民. 能量拓展型长中子计数器响应特性研究. 核技术. 2022(06): 53-60 . 百度学术其他类型引用(8)
-
下载:
下载:
必应学术
点击查看大图
计量
- 文章访问数: 1504
- HTML全文浏览量: 178
- PDF下载量: 428
- 被引次数: 12
