技术领域
[0001] 本
发明涉及数字
图像处理技术领域，具体涉及一种基于融合通道且对天空免疫的图像去雾方法。
背景技术
[0002] 近年来，雾霾天气较为常见，由于空气中大量
气溶胶粒子、微小
水滴等颗粒对光线的吸收和散射作用导致户外成像
质量严重退化，图像
颜色失真偏灰白，图像模糊细节不清晰，
对比度下降。不仅影响视觉效果，还直接限制了户外成像系统的效用的发挥，给后期图像处理和分析带来了困难。因此，对雾天图像进行清晰化处理十分必要。
[0003] 目前，图像去雾技术已经成为研究热点，对单幅图像去雾方法有图像增强方法和图像复原方法两种。图像增强方法主要是借助图像处理等手段，如直方图均衡化、频域滤波、Retinex等，通过削弱或去除一些不需要的信息，突出图像中的某些信息，在一定程度上能提高图像的对比度，改善视觉效果，但是没有从雾天图像的模糊退化机理入手，并不是实质上的去雾。
[0004] 图像复原方法去雾是从大气散射理论出发，考虑图像退化原因，建立雾天成像的散射模型，实现场景复原，能获得较为理想的去雾效果，其中以暗原色先验的去雾方法为代表。暗原色先验是一种对无雾图像的统计规律，首先粗估计出介质传播函数，然后利用软抠图原理或者引导
滤波器对透过率图像进行精细优化，该
算法取得令人叹服的去雾效果。
[0005] 然而在场景目标的
亮度与大气光相似的情况下，图像中存在大面积的天空区域基本找不到
像素值很小的点，此时暗原色假设失效。对于户外图像一般都会存在一定面积的天空图像，然而，以上无论是图像增强算法还是图像复原算法，这些算法有一个共同的缺点，就是会造成去雾后天空区域色彩失真，出现光晕现象和
块效应，严重影响去雾效果。相应的改进算法需要提取天空梯度信息等特征，进行天空识别，从而对天空和非天空区域进行分割，然后采取不同的透过率估计方法进行去雾处理，识别过程繁琐，计算量大。
发明内容
[0006] 本发明为解决现有图像去雾方法去雾后存在图像上天空区域色彩失真，出现光晕现象和块效应，导致去雾效果差以及现有方法存在识别过程繁琐计算量大等问题，提供一种基于融合通道且对天空免疫的图像去雾方法。
[0007] 基于融合通道且对天空免疫的图像去雾方法，该方法由以下步骤实现：
[0008] 步骤一、对原始有雾图像I(x)分别进行最小值滤波和最大值滤波，分别获得暗原色通道图像和亮原色通道图像；
[0009] 步骤二、取暗原色通道图像中亮度值最高的0.1％的像素，在原始有雾图像I(x)中寻找对应像素值的点，并取最大的像素值作为对大气光A的粗略估计；
[0010] 步骤三、引入修正因子ω，并根据步骤二获得的大气光A的粗略估计，获得暗通道透过率的粗略估计 和亮通道透过率的粗略估计 用下式表示为：
[0011]
[0012]
[0013] 式中，c为R、G、B三通道中的一个颜色通道，取暗原色通道图像的修正因子ωD＝0.8，亮原色通道图像的修正因子
[0014] 步骤四、根据步骤三获得暗通道透过率的粗略估计 和亮通道透过率的粗略估计 采用引导滤波分别对暗通道透过率的粗略估计 和亮通道透过率的粗略估计进行精细优化，分别得到暗通道透过率的精细估计 和亮通道透过率的精细估计[0015] 步骤五、判断暗通道透过率的精细估计 是否大于亮通道透过率的精细估计如果是，执行步骤六；如果否，则暗原色通道图像不是天空区域，融合透过率估计[0016] 步骤六、所述暗原色通道图像为天空区域，融合透过率估计 统计符合亮原色先验理论的天空区域对应的有雾图像部分的像素值，并将所述像素均值作为对大气光A的自适应准确估计；
[0017] 步骤七、根据步骤五、步骤六中综合得到的融合透过率 步骤六中获得的大气光A的自适应准确估计，采用大气散射成像模型获得去雾清晰化图像。
[0018] 本发明的有益效果：本发明基于大气散射模型，在暗原色先验的理论
基础上，从暗、亮双通道入手进行分析，利用亮原色先验弥补暗原色先验失效的情况，对天空部分图像有自适应的免疫性，并能对大气光A进行精准的自适应估计，有效提高了暗原色先验算法的普适性和鲁棒性。
附图说明
[0019] 图1为本发明的户外无雾图像暗原色图像与亮原色图像；其中，图1a为户外无雾原图像，图1b为暗原色图像，图1c为亮原色图像；
[0020] 图2为本发明的一种基于融合通道且对天空免疫的图像去雾方法
流程图；
[0021] 图3中分别为按照暗通道、亮通道、及融合通道对透过率图像进行估计的对比效果图；其中，图3a为原图，图3b为暗原色通道精细估计透过率效果图，3c为亮原色通道估计精细透过率效果图，图3d为融合通道估计透过率效果图；
[0022] 图4为采用本发明所述的基于融合通道且对天空免疫的图像去雾方法与现有直方图均衡化算法、暗通道先验算法对具有小面积天空雾图的去雾效果对比效果图；其中，图4a为原图，图4b为采用直方图均衡化算法对具有小面积天空雾图的去雾效果图；图4c为采用暗通道先验算法对具有小面积天空雾图的去雾效果图，图4d为采用本发明所述的基于融合通道且对天空免疫的图像去雾方法对具有小面积天空雾图的去雾效果图；
[0023] 图5为采用本发明所述的基于融合通道且对天空免疫的图像去雾方法与现有直方图均衡化算法、暗通道先验算法对具有大面积天空雾图的去雾效果对比图；其中，图5a为原图，图5b为采用直方图均衡化算法对具有大面积天空雾图的去雾效果图；图5c为采用暗通道先验算法对具有大面积天空雾图的去雾效果图，图5d为采用本发明所述的基于融合通道且对天空免疫的图像去雾方法对具有大面积天空雾图的去雾效果图；
[0024] 图6基于本算法与直方图均衡化算法、暗通道先验算法对具有大面积且天空与物体重叠的雾天图像去雾效果图；其中，图6a为原图，图6b为采用直方图均衡化算法对具有大面积且天空与物体重叠的雾天图像去雾效果图；图6c为采用暗通道先验算法对具有大面积且天空与物体重叠的雾天图像去雾效果图，图6d为采用本发明所述的基于融合通道且对天空免疫的图像去雾方法对具有大面积且天空与物体重叠的雾天图像去雾效果图。
具体实施方式
[0025] 具体实施方式一、结合图1至图6说明本实施方式，基于融合通道且对天空免疫的图像去雾方法，该方法由以下步骤实现：
[0026] 步骤一、对原始有雾图像I(x)分别进行最小值、最大值滤波，分别求取暗原色、亮原色通道的图像。
[0027] 步骤二、取暗通道图像按照像素值大小排序的前0.1％像素，并在原始有雾图像中寻找对应像素值点，取最大的像素值作为对大气光A的粗略估计。
[0028] 步骤三、由步骤二对大气光A的粗略估计，根据公式(5)、公式(10)分别得到暗通道的透过率 和亮通道的透过率 的粗略估计；
[0029] 其具体过程为：
[0030] 在
计算机视觉领域，大气散射成像模型得到广泛应用，如下式所示：
[0031] I(x)＝J(x)t(x)+A(1-t(x)) (1)
[0032] 其中，J(x)为要恢复的无雾图像，t(x)为透过率，c指R、G、B三通道中的一个颜色通道。将式(1)稍作处理，
变形为下式：
[0033]
[0034] 由暗原色先验理论，首先假设在每一个窗口内透过率t(x)为常数，定义为 A为常数，对式(2)两边求两次最小值运算，得到下式：
[0035]
[0036] 根据暗通道先验可知：
[0037]
[0038] 为了让人感到景深的存在，有必要在去雾的时有针对性的保留一部分
覆盖遥远景物的雾，引入一个在[0,1]之间的修正因子ω，最终可得到暗通道透过率的粗略估计[0039]
[0040] 其中，暗原色通道图像的修正因子ωD取值为：ωD＝0.8。
[0041] 对于户外场景中无雾图像统计发现，多数户外无雾图像的任意局部小块中，也存在一些像素，它们的某一个或几个颜色通道的强度值很高，甚至接近于255饱和值，在此我们称之为亮原色。
[0042] 经分析，户外无雾局部图像同时满足暗原色和亮原色，但也有局部图像不能同时满足这两种先验知识，如图1所示。显然，天空区域局部图像只满足亮原色先验，而不满足暗原色先验，那么采用暗通道将会出现错误，导致天空部分图像颜色失真、出现光晕现象和块效应。
[0043] 我们仿照暗通道推倒过程，对亮通道进行如下推倒：对式(2)两边求取两次最大值运算，得到式子：
[0044]
[0045] 由亮通道先验可知：
[0046]
[0047] 由(6)式可得：
[0048]
[0050]
[0051] 仿照暗通道理论，对图像进行归一化处理，并且加入修正因子，得式(10)，采用这个公式可得到利用亮通道透过率的粗略估计
[0052]
[0053] 其中，亮原色通道图像的修正因子ωL取值为：
[0054] 由式(5)可知： 是 减函数，说明当 取值越大，即透过率图越亮， 就越小，那么该局部图像符合暗通道先验规律概率越大，反之，透过率图越暗，该局部图像符合暗通道先验规律概率越小。
[0055] 由式(10)可知： 是 的增函数，说明当 取值越大，即透过率图越亮， 就越大，那么该局部图像符合亮通道先验规律概率越大，反之，透过率图越暗，该局部图像符合暗通道先验规律概率越小。
[0056] 步骤四、由步骤三对暗、亮通道的透过率的粗略估计 采用引导滤波分别对 进行精细优化估计，分别得到暗通道透过率的精细估计 和亮通道透过率的精细估计 对图3a中的雾天图像的原图进行精细透过率估计，分别得到暗原色通道精细估计透过率图如图3b，亮原色通道精细估计透过率图如图3c，图3d为本发明中融合通道估计透过率图。
[0057] 步骤五、判断 关系，来获取更加精确的融合透过率估计图 选取暗通道透过率的精细估计 和亮通道透过率的精细估计 中相对较大值(如式(11)所示)，得到最终的透过率估计 这样，可以使得暗原色和亮原色理论相互弥补，来获取更加精确的透过率估计图
[0058]
[0059] 如果 说明是天空部分，那么改进的融合透过率 反之，则不是天空部分，改进的融合透过率 最终得到改进的融合透过率图 如
图3d所示。
[0060] 步骤六、取步骤五中， 不符合暗原色先验理论的天空部分对应有雾图像部分的像素值，进行统计并求取均值，作为对大气光A的自适应准确估计。
[0061] 步骤七：取步骤五、步骤六中融合透过率图 大气光A的精细自适应估计值，根据大气散射成像模型得到去雾清晰化图像。
[0062] 本实施方式中，选取三种存在天空的有雾图像，一种是存在小面积天空的雾天图像，如图4所示；第二种是具有大面积天空的雾天图像，如图5所示；第三种是具有大面积天空且天空与物体重合的雾天图像，如图6所示；
[0063] 选取图像增强的典型性算法-直方图均衡化算法、图像复原中的典型去雾算法-暗通道先验算法，和本发明的一种基于融合通道且对天空免疫的图像去雾算法对雾天的图像进行清晰化处理，并对去雾效果进行对比，处理结果参见图4-图6。
[0064] 对比去雾效果表明：在雾天图像中，无论存在天空面积大小与否，直方图均衡化算法、暗通道先验算法对天空部分图像的处理结果都会存在色彩失真、光晕现象和块效应，参见图4-图6中的4b、4c、5b、5c、6b和6c，严重影响图像的去雾效果。
[0065] 采用本实施方式所述基于融合通道且对天空免疫的图像去雾算法处理得到的图像，无论天空部分图像面积大小，有效避免了以上算法对雾天图像天空部分的失真，得到整体和谐清晰的图像，说明了本发明算法对雾天图像处理的天空部分的免疫性。
