改进暗通道先验和策略性融合的图像去雾算法.pdf

1、第 49卷 第 10期2023年 10月Computer Engineering 计算机工程改进暗通道先验和策略性融合的图像去雾算法王效灵，胡志杰，徐帅帅，黄浩如（浙江工商大学 信息与电子工程学院，杭州 310018）摘要：针对图像去雾领域的暗通道先验算法存在光晕效应、颜色失真及对天空区域处理不佳等问题，提出一种改进暗通道先验和策略性融合的图像去雾算法。依据暗通道先验原理得到多尺寸最小值滤波窗口下的透射率和对应大气光值，设计基于结构相似性指标的调节和侧重因子用于多窗口透射率拟合。利用非线性规划模型和偏度理论计算全局大气光值，引入并调整置信度对天空区域的透射率进行快速补偿，结合大气散射模型恢复去

2、雾图。将去雾图转换成 HSV 模型并对亮度 V空间进行增强，同时对增强前后的图像进行策略性融合。实验结果表明，相比于现有的暗通道先验去雾算法，该算法的可见边增率、信息熵和平均梯度分别提升了 79%131%、3%9%、39%81%，具有较好的去雾效果，避免了光晕效应和颜色失真现象，适用于处理包含不同尺寸天空区域的雾图。关键词：图像去雾；暗通道先验；结构相似性；非线性规划；置信度；策略性融合开放科学（资源服务）标志码（OSID）：中文引用格式：王效灵，胡志杰，徐帅帅，等.改进暗通道先验和策略性融合的图像去雾算法 J.计算机工程，2023，49（10）：212-221.英文引用格式：WANG X L

3、，HU Z J，XU S S，et al.Image dehazing algorithm using improved dark channel prior and strategic fusion J.Computer Engineering，2023，49（10）：212-221.Image Dehazing Algorithm Using Improved Dark Channel Prior and Strategic FusionWANG Xiaoling，HU Zhijie，XU Shuaishuai，HUANG Haoru（School of Information and E

4、lectronic Engineering，Zhejiang Gongshang University，Hangzhou 310018，China）【Abstract】To address the problems of halo effect，color distortion，and poor treatment of sky regions in the Dark Channel Prior（DCP）algorithm for image dehazing，a novel image dehazing algorithm using an improved DCP and strategi

5、c fusion is proposed.Based on the DCP principle，the transmission rate and corresponding atmospheric light value are obtained using a multi-scale minimum value filtering window.A design utilizing Structural Similarity（SSIM）index is employed to adjust and emphasize factors for multi-transmission rate

6、fitting.A nonlinear optimization model and the skewness theory are utilized to calculate the global atmospheric light value，and confidence is introduced and adjusted to rapidly compensate for the transmission rate in sky regions.Through combination with an atmospheric scattering model，the dehazed im

7、age is restored.The dehazed image is then transformed into the HSV color space to enhance the brightness in the V channel，followed by a strategic enhancement fusion.Experimental results show that compared with the existing DCP dehazing algorithms，the proposed algorithm achieves significant improveme

8、nts in visible edge increment rate，information entropy，and average gradient by 79%-131%，3%-9%，and 39%-81%，respectively.This demonstrates its superior dehazing performance by avoiding halo effects and color distortion phenomena，and it is suitable for handling hazy images containing sky regions of dif

9、ferent sizes.【Key words】image dehazing；Dark Channel Prior（DCP）；Structural Similarity（SSIM）；Non-Linear Programming（NLP）；confidence；strategic fusionDOI：10.19678/j.issn.1000-3428.00659080概述 雾霾是常见的大气污染现象，是由空气中的灰尘硫酸还有硝酸等颗粒物组成的溶胶系统造成的视觉障碍，导致雾霾天成像结果低质，不利于计算机视觉的后续研究，因此图像去雾具有重要的研究意义。基金项目：浙江省重点研发科技计划项目（2021C0

10、3166）。作者简介：王效灵（1967），男，教授级高级工程师、博士，主研方向为图像处理、信号处理、智能监控；胡志杰、徐帅帅、黄浩如，硕士研究生。收稿日期：2022-10-04 修回日期：2022-11-23 Email：图形图像处理文章编号：1000-3428（2023）10-0212-10 文献标志码：A 中图分类号：TP391目前，图像去雾领域大致有基于图像增强、基于深度学习和基于图像复原三大研究方向。基于图像增强的方法有直方图均衡1、基于 Retinex 理论2-4等，虽然雾图的对比度等得到了增强，也改善了视觉效果，但无法真实地实现图像去雾。基于深度学习的方法是通过训练网络模型实现端到

11、端去雾，缺点是训练耗时，且目前合成的雾数据集多而自然雾数据集少。ZHANG 等5提出 DCPD 模型，结合语义分割的生成器和判别器实现去雾。CHEN 等6提出一种基于生成对抗网络的无监督学习训练去雾模型GCANet。REN 等7提出一种多尺度特征的 CNN 模型来训练雾图的透射率进行去雾。ZHAO 等8提出两阶段弱监督 RefineDNet模型，一阶段采用暗通道恢复可见性，二阶段使用未配对的雾状清晰图像进行对抗性学习提高复原真实性。YANG 等9提出可实现自增强的去雾网络框架，但对过亮区域的透射率估计偏高，引起深度网络对这些区域的深度预测值偏低。LI等10提出一种沙漏式结构且基于残差网络的 H

12、ourglass-DehazeNet，在多尺度上进行灵活的信息交互和聚合。基于图像复原的方法是使用特定物理成像模型实现雾气的退化，是目前流行的真实去雾方法。LIU 等11根据亮度和饱和度呈正相关的规律估计透射率，但对于景深处的去雾效果差。HE等12提出经典的暗通道先验（Dark Channel Prior，DCP）去雾算法，利用统计先验知识并结合大气散射模型高效复原无雾图，但出现了一定的光晕伪影、色彩偏移等问题，且在天空区域中会失效。RAIKWAR等13将透射率估计问题映射为模糊图像与无雾图像间最小值通道差的估计问题。JU 等14在大气散射模型中增加光吸收系数来提升恢复图像的可见度，但其拉伸策

13、略的全局性较差。何涛等15利用可自动调整的雾气保留参数和容差机制改进了暗通道先验算法，改善对天空区域的处理，但容差系数的固定选择无法适用于不同场景的雾图。LI等16和王硕等17运用暗通道置信度的理论，自适应修正了 DCP 不适用区域的透射率，但处理后图像的景深处去雾效果差。王国栋等18利用散射介质成像模型，基于低秩分解并结合像素置乱方法进行去雾，但对雾气分布不均匀的图像处理不佳。吴向平等19利用深度学习模型得到雾图的物理景深估计模型，基于景深信息进行自适应图像 Retinex操作实现去雾。本文提出基于结构相似性（Structural Similarity，SSIM）指标和侧重因子的稳定透射率拟

14、合算法和结合非线性规划（Non-Linear Programming，NLP）模型和偏度知识的新大气光值确定方法，结合并调整置信度解决 DCP 在天空区域的失效问题，同时对图像进行策略性的增强融合。1 算法基础1.1大气散射模型计算机视觉和图形图像领域一般利用大气散射模型来描述雾天图像降质的过程，可表示如下：I(x)=J(x)t(x)+A1-t(x)（1）其中：x表示图像像素；I(x)表示雾图；J(x)表示无雾图；A表示全局大气光值；t(x)表示透射率；A1-t(x)表示大气光幕20。在均匀介质中，t(x)可表示如下：t(x)=exp-d(x)（2）其中：表示介质散射系数；d(x)表示成像景深

15、。因此，根据大气散射模型，输入雾图、大气光值和透射率就能反解方程得到去雾图。1.2暗通道先验原理除室外无雾图像的非天空区域外，图像的其他部分至少有一个颜色通道强度趋近于 0，对于雾图I(x)，暗通道可表示如下：Jdark(x)=mini (x)minc(RGB)Ic(i)0（3）其中：Jdark(x)表示暗通道；(x)表示最小值滤波窗口；c表示 R、G、B颜色通道。结合式（1）和式（3），对于给定的全局大气光值A，计算出雾图的透射率t(x)：t(x)=1-Jdark(x)A（4）其中：表示雾气保留系数，一般取为 0.95，维持图像的少量雾气以保证真实性。根据大气散射模型，得出最终的去雾图像J(

16、x)可表示如下：J(x)=I(x)-Amaxt(x)t+A（5）其中：t设置为 0.1，避免分母为 0。1.3结构相似性结构相似性21从亮度、对比度、结构化 3 个层面，对A、B进行相似性的客观评估，可表示如下：SSIM(AB)=Lig(AB)Con(AB)Str(AB)（6）其中：Lig(AB)表示亮度相似度；Con(AB)表示对比度相似度；Str(AB)表示结构化相似度。A、B两图SSIM 的值域为0，1，SSIM(AB)越大，说明图像A、B整体越相似。2改进暗通道先验和策略性融合的图像去雾算法 改进暗通道先验和策略性融合的图像去雾算法流程如图 1 所示。首先，获取多窗口下 DCP 算法的

17、透 射 率ti(x)和 相 应 大 气 光 值Ai；接 着，使 用 基 于SSIM 的调节因子i和侧重因子i拟合透射率；然后，假设一个非线性规划模型，通过系数约束和偏度约束，从Ai组中计算A，继续引入并调整置信度补偿天空区域的透射率，并快速引导滤波细化；最后，反解大气散射模型，转换为 HSV 模型增强亮度 V，结合置信度特点进行增强前后图像的策略性融合，得到最终的去雾图。第 49卷 第 10期王效灵，胡志杰，徐帅帅，等：改进暗通道先验和策略性融合的图像去雾算法目前，图像去雾领域大致有基于图像增强、基于深度学习和基于图像复原三大研究方向。基于图像增强的方法有直方图均衡1、基于 Retinex 理

18、论2-4等，虽然雾图的对比度等得到了增强，也改善了视觉效果，但无法真实地实现图像去雾。基于深度学习的方法是通过训练网络模型实现端到端去雾，缺点是训练耗时，且目前合成的雾数据集多而自然雾数据集少。ZHANG 等5提出 DCPD 模型，结合语义分割的生成器和判别器实现去雾。CHEN 等6提出一种基于生成对抗网络的无监督学习训练去雾模型GCANet。REN 等7提出一种多尺度特征的 CNN 模型来训练雾图的透射率进行去雾。ZHAO 等8提出两阶段弱监督 RefineDNet模型，一阶段采用暗通道恢复可见性，二阶段使用未配对的雾状清晰图像进行对抗性学习提高复原真实性。YANG 等9提出可实现自增强的去

19、雾网络框架，但对过亮区域的透射率估计偏高，引起深度网络对这些区域的深度预测值偏低。LI等10提出一种沙漏式结构且基于残差网络的 Hourglass-DehazeNet，在多尺度上进行灵活的信息交互和聚合。基于图像复原的方法是使用特定物理成像模型实现雾气的退化，是目前流行的真实去雾方法。LIU 等11根据亮度和饱和度呈正相关的规律估计透射率，但对于景深处的去雾效果差。HE等12提出经典的暗通道先验（Dark Channel Prior，DCP）去雾算法，利用统计先验知识并结合大气散射模型高效复原无雾图，但出现了一定的光晕伪影、色彩偏移等问题，且在天空区域中会失效。RAIKWAR等13将透射率估计

20、问题映射为模糊图像与无雾图像间最小值通道差的估计问题。JU 等14在大气散射模型中增加光吸收系数来提升恢复图像的可见度，但其拉伸策略的全局性较差。何涛等15利用可自动调整的雾气保留参数和容差机制改进了暗通道先验算法，改善对天空区域的处理，但容差系数的固定选择无法适用于不同场景的雾图。LI等16和王硕等17运用暗通道置信度的理论，自适应修正了 DCP 不适用区域的透射率，但处理后图像的景深处去雾效果差。王国栋等18利用散射介质成像模型，基于低秩分解并结合像素置乱方法进行去雾，但对雾气分布不均匀的图像处理不佳。吴向平等19利用深度学习模型得到雾图的物理景深估计模型，基于景深信息进行自适应图像 Re

21、tinex操作实现去雾。本文提出基于结构相似性（Structural Similarity，SSIM）指标和侧重因子的稳定透射率拟合算法和结合非线性规划（Non-Linear Programming，NLP）模型和偏度知识的新大气光值确定方法，结合并调整置信度解决 DCP 在天空区域的失效问题，同时对图像进行策略性的增强融合。1 算法基础1.1大气散射模型计算机视觉和图形图像领域一般利用大气散射模型来描述雾天图像降质的过程，可表示如下：I(x)=J(x)t(x)+A1-t(x)（1）其中：x表示图像像素；I(x)表示雾图；J(x)表示无雾图；A表示全局大气光值；t(x)表示透射率；A1-t(x

22、)表示大气光幕20。在均匀介质中，t(x)可表示如下：t(x)=exp-d(x)（2）其中：表示介质散射系数；d(x)表示成像景深。因此，根据大气散射模型，输入雾图、大气光值和透射率就能反解方程得到去雾图。1.2暗通道先验原理除室外无雾图像的非天空区域外，图像的其他部分至少有一个颜色通道强度趋近于 0，对于雾图I(x)，暗通道可表示如下：Jdark(x)=mini (x)minc(RGB)Ic(i)0（3）其中：Jdark(x)表示暗通道；(x)表示最小值滤波窗口；c表示 R、G、B颜色通道。结合式（1）和式（3），对于给定的全局大气光值A，计算出雾图的透射率t(x)：t(x)=1-Jdark

23、(x)A（4）其中：表示雾气保留系数，一般取为 0.95，维持图像的少量雾气以保证真实性。根据大气散射模型，得出最终的去雾图像J(x)可表示如下：J(x)=I(x)-Amaxt(x)t+A（5）其中：t设置为 0.1，避免分母为 0。1.3结构相似性结构相似性21从亮度、对比度、结构化 3 个层面，对A、B进行相似性的客观评估，可表示如下：SSIM(AB)=Lig(AB)Con(AB)Str(AB)（6）其中：Lig(AB)表示亮度相似度；Con(AB)表示对比度相似度；Str(AB)表示结构化相似度。A、B两图SSIM 的值域为0，1，SSIM(AB)越大，说明图像A、B整体越相似。2改进暗

24、通道先验和策略性融合的图像去雾算法 改进暗通道先验和策略性融合的图像去雾算法流程如图 1 所示。首先，获取多窗口下 DCP 算法的透 射 率ti(x)和 相 应 大 气 光 值Ai；接 着，使 用 基 于SSIM 的调节因子i和侧重因子i拟合透射率；然后，假设一个非线性规划模型，通过系数约束和偏度约束，从Ai组中计算A，继续引入并调整置信度补偿天空区域的透射率，并快速引导滤波细化；最后，反解大气散射模型，转换为 HSV 模型增强亮度 V，结合置信度特点进行增强前后图像的策略性融合，得到最终的去雾图。2132023年 10月 15日Computer Engineering 计算机工程2.1多窗口

25、透射率拟合暗通道先验算法去雾效果如图 2 所示（彩色效果见 计算机工程 官网 HTML 版）。选择 1515 像素的最小值滤波窗口对最小值通道进行处理，具有较好的去雾能力且复原度高，但透射率精细度较低、信息丢失导致去雾图存在光晕现象，如图 2（c）所示。通过逐渐缩小最小值滤波窗口的尺寸，虽然透射率越来越高，但复原效果越来越差，图像失真、过饱和严重，如图 2（i）所示。根据不同尺寸的最小值滤波窗口，基于暗通道先验算法可得到不同效果倾向的复原图。因此，本文提出一种基于 SSIM 指标的调节和侧重因子并用于拟合多窗口透射率的算法。首先，选取 33、66、99、1212、1515 像素的最小值滤波窗口

26、，使用式（3）和式（4）处理雾图，分别得到对应的透射率t0(x)、t1(x)、t2(x)、t3(x)、t4(x)和大气光值A0、A1、A2、A3、A4，其中，t0(x)精细度相对较高但复原度相对较低，t4(x)反之。接着，结合式（6）定义调节因子i：i=maxSSSIMti(x)t0(x)SSSIMti(x)t4(x)minSSSIMti(x)t0(x)SSSIMti(x)t4(x)（7）其 中：SSSIMti(x)t0(x)，i 04为 相 对 精 细 程 度；SSSIMti(x)t4(x)为相对复原程度；i=1表示第i幅雾图的透射率在样本组中没有相对信息提供，i越大，说明能提供越多的相对信

27、息。然后，计算全局侧重数：=i=0NSSSIMti(x)t0(x)i=0NSSSIMti(x)t4(x)（8）其中：N=4；1表示在拟合过程中侧重于相对精细程度大的透射率，反之侧重于相对复原效果好的透射率。侧重因子i可表示如下：i=SSSIMti(x)t0(x)SSSIMti(x)t4(x)且 11SSSIMti(x)t0(x)SSSIMti(x)t4(x)且 11，其他（9）其中：若全局侧重数1，当透射率ti(x)对应的相对精细程度大于相对复原程度，则令其侧重因子i=，否则保持为 1；若全局侧重数1，当ti(x)对应的相对复原程度大于相对精细程度，则令侧重因子i=1，否则保持为 1。结合式（

28、7）式（9），拟合得到新的透射率t(x)：t(x)=i=0Nii ti(x)i=0Nii（10）表 1 给出了原始雾图的各透射率，可以看出通过i和i拟合出的透射率相对信息最多，稳定性更好。图 2暗通道先验算法去雾效果Fig.2Dehazing effect of dark channel prior algorithm图 1改进暗通道先验和策略性融合的图像去雾算法流程Fig.1Procedure of image dehazing algorithm using improved dark channel prior and strategic fusion214第 49卷 第 10期王效灵，

29、胡志杰，徐帅帅，等：改进暗通道先验和策略性融合的图像去雾算法2.2基于 NLP和偏度的大气光值计算利 用 暗 通 道 先 验 算 法 得 到 大 气 光 值A0、A1、A2、A3和A4，因为理论稳定或者波动幅度不大，所以称(A0A1A2A3A4)为大气光值组。当它们全部相等时，则认为得到稳定可靠的全局大气光值；当它们不完全相等时，则必然存在最大值、最小值导致大气光值误差。因此，基于 NLP 模型求取雾图的全局大气光值A，使A在一定约束条件下保持理论合理性，假设：max Z=Aavg-AminAmax-Amin x1+Amax-AavgAmax-Amin x2s.t.x1AminAavg x2A

30、avgAmax|PZ|0（11）其中：Z是待求取符合约束条件的最大值，避免了过小 的 不 准 确 大 气 光 值；Amin、Aavg和Amax分 别 是(A0A1A2A3A4)中的最小值、平均值和最大值；通过差值步长比Aavg-AminAmax-Amin、Amax-AavgAmax-Amin控制求和权重进行系数约束，增加异常Amin或者Amax情况下的容错性；条件约束是计算偏度PZ并迫使其尽可能靠近 0。亮通道原理22指出：雾图的像素点的亮通道强度趋近于其无雾时的大气光值强度。雾图亮通道的所有强度值可视作一个大气光值组。从目前规模较大的单幅图像去雾数据集 RESIDE23中随机取 4 000幅

31、雾图，根据亮通道原理22求取各像素的亮通道强度组成大气光值组，利用 MATLAB 库函数 Skewness计算各组的偏度值，4 000幅随机雾图大气光值组的偏度统计分布如图 3所示。图 34 000幅雾图大气光值组的偏度统计分布Fig.3Skewness statistical distribution of atmospheric light group in 4 000 hazy images增加样本继续进行统计，可以发现：对于不完全相等的随机大气光值组而言，理论成立的样本光越多，其组元素值表现为相对均匀地分布在均值两侧，即偏度趋近于 0。因此，使用Z对(A0A1A2A3A4)随机组进行偏

32、度校正。偏度PZ可表示如下：PZ=1N+1i=0N+1()AZi-AZavgZ3Z=1N+1i=0N+1(AZi-AZavg)2（12）其 中：AZi表 示 新 大 气 光 值 组(A0A1A2A3A4Z)元 素；Z表 示 该 组 的 标 准 差；AZavg表 示 该 组 的均值。最后，将式（12）代入式（11）求解 NLP 模型符合条件的最优Z值，并取当前组的AZavg为全局大气光值A。2.3天空区域透射率优化暗通道先验算法不能处理天空区域，引入置信度17并调整，以补偿这些区域的透射率。首先，结合二维高斯函数，利用 HSV模型的 V、S空间信息，计算置信度C(x)：C(x)=max0.51-

33、exp-(Vmin(x)-1)2c-Smax(x)2c（13）其中：Vmin(x)、Smax(x)是分别对 V、S 分量进行最小值和最大值滤波的结果；文献 17 令 V、S 分量滤波窗口的大小等于其计算暗通道时的窗口大小，由于采用多窗口拟合透射率，为了稳定尺寸，选取第 2.1 节的中间窗口大小，即 99 像素；C(x)设定阈值为 0.5，防止透射率过补偿；c是自适应调节参数。c可表示如下：c=max0.05min 0.01(1-A255)0.2（14）其中：A为全局大气光值，可自动调节c以取得更好的算法效果。包含不同尺寸的天空区域的雾图和对应图如图 4 所示（彩色效果见 计算机工程 官网 HT

34、ML版）。观察图 4（d）并结合像素点值可得一般区域的置信度非常大并接近 1，利用该特点进行透射率的快速补偿 如式（15）所示。采用快速引导滤波器进一步细化得到最终的透射率，反解大气散射模型得到去雾图。对比图 4（b）和图 4（e）可以看出，在对置信度调整后，既能补偿天空区域的透射率，又能维持一般区域的稳定性，复原效果更佳。t*(x)=min1t(x)+1-C(x)（15）表 1原始雾图的透射性能 Table 1Transmission performance of the original hazy images窗口大小/像素33669912121515拟合值相对精细程度1.000 00.7

35、11 30.590 70.549 10.534 60.766 0相对复原程度0.534 60.643 60.740 30.867 71.000 00.836 3总计1.534 61.354 91.331 01.416 81.534 61.602 32152023年 10月 15日Computer Engineering 计算机工程2.4图像融合增强基于暗通道先验算法的复原图像往往整体偏暗，因此将复原图转换成 HSV模型，对亮度 V空间进行增强。在 V 空间进行对比度受限的自适应直方图均衡化（Contrast Limited Adaptive Histogram Equalization，CLA

36、HE），将对比度限制为 0.01防止图像过增强，使用 22像素的滑块进行局部均衡。如图 5 所示（彩色效果见 计算机工程 官网HTML 版），对比图像增强前后，非天空场景的视觉效果得到显著提升，但天空区域出现了一定的伪影和纹理，CLAHE 造成天空失真，反观增强前的天空场景保留完好，因此设计一种图像融合策略。首先，将置信图的天空部分（背景色）和非天空部分（前景色）天然粗分割，可直接用最大类间方差法快速进行二值分割。对于置信图C(x)，将非天空区域和天空背景的阈值记作k，属于非天空区域的像素点数占C(x)总像素数的比值为0，平均像素值为0。同理，天空背景有1和1。类间方差g可表示如下：g=0 1

37、(0-1)2（16）遍历C(x)的像素值，得到使g最大情况下的最优 阈 值k，分 割 置 信 图 得 到 二 值 图K(x)，如 图 6所示。令 CLAHE前后的去雾图分别为Jbef(x)和Jaft(x)，为了减弱K(x)的 0-1 突变边缘造成最终图像J(x)在拼接处的明显痕迹以及违和感，采用各向同性的高斯平滑核增加标准差来模糊置信度二值图，可表示如下：Kgaus(x)=K(x)12 exp(-x222)（17）设置方形高斯平滑核的尺寸如下：s=2 2+1（18）越小，拼接区域越突兀；越大，Jbef(x)天空区域的权值越不稳定，且会弱化图像整体增强效果，因此一般取 15。最终图像J(x)可表

38、示如下：J(x)=Jbef(x)1-Kgaus(x)+Jaft(x)Kgaus(x)（19）融合图像的天空背景真实准确，非天空场景得到很好的视觉性增强且融合效果不突兀。3 实验结果与分析为了验证所提算法的有效性，搭建实验平台进行验证。硬件平台配置为 Intel Core i7-6700HQ2.60 GHz，12 GB 内存，Windows 10 专业版，软件平台为 MATLAB R2019a。实验图片来自去雾领域的图 5图像增强前后对比Fig.5Comparison of images before and after enhancement图 4包含不同尺寸的天空区域的雾图和对应图Fig.4

39、Hazy images with different sizes of sky regions and related images图 6置信图的二值分割Fig.6Binary segmentation of confidence image216第 49卷 第 10期王效灵，胡志杰，徐帅帅，等：改进暗通道先验和策略性融合的图像去雾算法相关论文和单幅图像去雾数据集 RESIDE23，混合主观测试子集 HSTS包含了 10张真实雾图以及 10张对应的已知样本真值图片（合成），从主客观角度与其他算法进行比较分析。3.1主观评价通过主观评价能直接快速对复原效果进行分析，所提算法与其他去雾算法的复原效

40、果比较如图 7所示（彩色效果见 计算机工程 官网 HTML版）。如 7 图所示，image 1、2、3 是不包含天空区域的雾图，image 4、5、6雾图分别存在不同尺寸的天空场景。由图 7 可以看出：文献 6 算法处理的 image 4、5、6天空区域效果不稳定，其中，image 4的天空严重失真，处理的 image 1、2、3 存在一定的颜色偏移，视觉效果较差；文献 12 算法由于较低的透射率导致image 3、5、6存在明显光晕效应，因为图像中高亮像素的影响导致 image 1 的颜色偏移、image 2、3 图像过暗，image 4、5、6 的天空失真严重；文献 15 算法引入容差机制

41、但由于容差系数过大，导致 image 1、2、3、4的去雾效果不理想；文献 17 算法用自适应校正透射率算法进行去雾，观察 image 1、3、5存在景深突 变 区 域 雾 气 残 留 问 题，image 2、6 颜 色 过 暗，image 4 天空轻微失真。经过大量雾图的测试实验可知，所提算法对包含不同尺寸天空区域的雾图表现良好，能够复原并增强场景细节和色彩，同时不失去图像去雾真实性。3.2客观评价相较于主观评价，客观评价可以更好地分析算法的性能表现。使用 HAUTIRE 等24提出的可见边增率（l）、信息熵（H）、平均梯度（G）3 个性能指标进行客观评价。l 计算图像的边缘信息变化程度，值

42、越大，图像复原效果越好25，可表示如下：l=lr-lolo（20）其 中：lr和lo分 别 为 恢 复 图 像 和 原 始 图 像 的 可 见边数。H 计算图像的信息量程度，值越大，包含的信息越多，可表示如下：H=-i=1MR(ai)logaR(ai)（21）图 7不同算法的复原效果比较 1Fig.7Comparison 1 of restoration effects of the different algorithms2172023年 10月 15日Computer Engineering 计算机工程其中：M是不同像素点值ai的总数；R(ai)表示每个ai值在所有像素中出现的概率。G 计

43、算图像的纹理复杂程度，值越大，图像越清晰，分辨率越高，可表示如下：G=x=1m-1y=1n-1G(xy)-G(x+1y)2+G(xy)-G(xy+1)22(m-1)(n-1)（22）其中：m和n分别是图像的总行数和总列数；G(xy)为图像第x行和第y列的灰度值。根据图 7，各算法处理得到的图像与各雾图之间的 l 指标如表 2 所示，其中，最优指标值用加粗字体标示。由表 2 可以看出，所提算法对不同场景下的雾图均有较高的 l值，较其他算法提升约 79%131%，在保护雾图强边缘细节的同时，能够恢复出大量被雾气影响的边缘信息。根据图 7，各算法处理得到的图像与各雾图之间的 H、G 指标如表 3、表

44、 4 所示，其中，最优指标值用加粗字体标示。由表 3、表 4 可以看出：其他算法的处理导致一些雾图的失真或是本该平滑的区域（如天空）出现过多纹理，存在拉高了 H 和 G 指标的情况；所提算法的 H、G 表现优势明显，较其他算法分别提升 3%9%、39%81%，同时保留了平滑区域，符合主观评价中图像信息丰富、分辨率高、天空自然的复原效果。通过综合比较分析可知，所提算法处理不同场景下的雾图时去雾效果真实、细节丰富、图像清晰，取得了很好的去雾效果。3.3定性与定量比较及消融实验图 8 给出了文献 8，10，19 以及所提去雾算法在 HSTS 测试集上的定性结果（彩色效果见 计算机工程 官网 HTML

45、 版）。对照各样本真值图 8（f）可以看出：文献 8 算法在天空区域有云团和图像阴暗处的处理效果不佳；文献 10 算法的沙漏式多尺度操作会引起一些场景细节的丢失；文献 19算法的整体色彩还原度欠佳，采用双线性插值映射法减弱块与块间的不连续边缘能力有限。SSIM 和峰值信噪比（Peak Signal-to-Noise Ratio，PSNR）是两种广泛使用的评价去雾图像质量的指标。各算法处理得到的图像与各雾图样本真值之间的 SSIM、PSNR 指标定量比较如表 5 所示，其中，最优指标值用加粗字体标示。由表 5 可以看出：文献8，10算法在若干图上取得了较高的 SSIM 或PSNR指标，该结果是符

46、合预期的，因为文献 8 算法使用场景数据集 ITS 的全部样本进行训练，在一定程度上增强了模型对测试子集 HSTS 的有效性，文献 10 算法也是类似的；所提算法虽然基于暗通道先验理论，但是对于基于深度学习去雾常用的开源数据集，也具有一定的性能优势。为了验证改进模块的有效性，进行消融实验。消融实验进行如下分组：组 1（DCP），组 2（DCP+透射率拟合模块），组 3（DCP+天空区域优化模块），组 4（DCP+透射率拟合模块+天空区域优化模块），组 5（DCP+透射率拟合模块+天空区域优化模块+策略性融合模块，即所提算法）。消融实验结果如表 6 所示，可以看出：设置对照组 1、2和 3，透射

47、率拟合模块表 2不同算法处理各雾图的 l指标比较 Table 2Comparison of l indexes for processing various hazy images using different algorithmsimage123456文献6算法-0.065 90.048 60.136 40.016 6-0.013 7-0.460 1文献12算法0.013 30.097 70.122 80.133 4-0.022 8-0.001 8文献15算法0.004 50.167 80.011 8-0.014 8-0.206 1-0.282 6文献17算法0.014 3-0.100 4

48、0.171 70.090 7-0.096 1-0.253 6所提算法0.034 70.288 00.229 40.213 40.036 00.231 8表 3不同算法处理各雾图的 H指标比较 Table 3Comparison of H indexes for processing various hazy images using different algorithmsimage123456文献6算法7.414 07.387 67.002 87.510 07.445 06.277 5文献12算法7.400 97.039 96.707 06.925 97.190 66.487 2文献15算法

49、6.968 37.069 57.251 17.277 17.044 26.336 5文献17算法7.133 97.272 46.647 37.306 67.291 86.429 4所提算法7.796 17.497 17.505 07.786 17.670 37.022 5表 4不同算法处理各雾图的 G指标比较 Table 4Comparison of G indexes for processing various hazy images using different algorithmsimage123456文献6算法0.030 00.027 00.010 40.035 50.015 00

50、.011 7文献12算法0.032 60.026 60.013 50.034 80.013 00.013 9文献15算法0.024 10.015 60.011 80.027 80.009 50.010 4文献17算法0.031 60.022 00.008 50.033 00.009 80.012 3所提算法0.053 60.044 30.040 00.041 30.030 90.032 7218第 49卷 第 10期王效灵，胡志杰，徐帅帅，等：改进暗通道先验和策略性融合的图像去雾算法（增加透射率的相对信息）和天空区域优化模块（解决 DCP 不适用于天空、明亮区域的问题）均能提升DCP 算法的性