一种基于模型与直方图及灰度增强的图像雾霾去除方法与流程

本发明属于图像雾霾信号与信息处理及遥感图像处理技术领域，具体涉及一种基于模型与直方图及灰度增强的图像雾霾去除方法。

背景技术

随着工业化、城镇化和现代化的发展，雾霾天气的明显增加，持续时间增长。雾霾天气的浓度不同，对光学成像的影响效果也是不同的。在轻雾条件下，获取的光学图像，有一种朦胧感觉。如果受到重度霾天气的影响，获取的遥感图像往往不清晰，明显的偏白，特别是远处的部分。这不仅严重影响图像的质量，而且给遥感图像的处理、应用和解译带来较大的困难，因此，去除遥感图像雾霾影响并对其进行增强处理是遥感图像预处理的重要内容。遥感图像雾霾的去除，本质上就是图像的增强或者图像恢复，目的是提高图像的清晰度和质量。在雾霾天气下获得的遥感图像，通常是一种灰蒙蒙的感觉，细节信息变得模糊，甚至看不到细节信息，整幅图像呈现一片偏白的现象。原因是雾霾中存在一些乳胶状的颗粒，影响太阳光的照射和物体的反射，对于地物来说，相当于披了一件外套，使到达地物的光的饱和度降低，产生了模糊现象。雾霾越重，这种现象越严重。

众多学者针对图像雾霾去除方法和理论进行了探讨，有学者从物理模型、光照原理和散射原理等方面进行研究，分别提出了典型的retinex去雾霾算法、暗原色先验原理去雾霾算法和同态滤波去雾霾算法。虽然它们能在一定程度上减少雾霾的影响，但每种方法的不足也很明显，例如retinex算法不仅原理复杂，运算量相对较大，而且在复杂情况下会出现光晕现象；暗原色原理法参数估计的准确性和调整对去雾霾结果有较大影响；同态滤波是在频率域内对部分频率进行处理，以达到增强效果，但这个频率阈值经验性比较强。还有一种典型的传统方法是直方图均衡化处理，这种方法简单，速度快，对于轻雾的去除，效果很好，但是对于复杂的情况，效果不太理想。最近出现了一些新的理论和方法，例如retinex算法和暗原色算法的结合，小波分解法和retinex算法的结合，稀疏表示法和暗原色法的结合，以及机器学习的方法等。这些方法对于处理薄雾，通常能获得好的结果，但是对于复杂情况下，即浓雾或浓霾条件下，结果通常不太理想。所以，这些方法的普适性差，应用范围受限；同时不同方法的组合以及机器学习的方法都是以时间为代价来获得去雾结果。

技术实现要素：

为了克服上述现有技术的缺点，本发明的目的在于提供了一种基于模型与直方图及灰度增强的图像雾霾去除方法，能够有效地去除图像中的雾霾影响，对于没有受到雾霾影响的图像，能有效地增强或改善图像质量。

为了达到上述目的，本发明采取的技术方案为：

一种基于模型与直方图及灰度增强的图像雾霾去除方法，包括以下步骤：

步骤1：输入光学遥感图像；

步骤2：判断输入光学遥感图像的类型，是灰度图像还是rgb彩色图像；

步骤3：对于全色的灰度图像，需转换成rgb图像；

步骤4：把rgb图像转换成hsi图像，提取hsi图像的三个分量，即h分量图、s分量图和i分量图；

步骤5：对i分量图进行直方图规定化处理；

步骤6：对i分量图进行灰度对比度范围调整；

步骤7：对i分量图进行雾霾去除处理；

步骤8：对s分量图进行雾霾去除处理，与步骤7相同，只是把s分量图代替i分量图；

步骤9：把hsi图像转换成rgb图像；

步骤10：如果处理的图像是rgb图像，则处理结束，直接输出结果；

步骤11：如果处理的图像是全色灰度图像，则要由rgb图像转换成灰度图像，把其作为处理后的结果，然后直接输出。

所述的步骤1中光学遥感图像受到雾霾影响，或没有受到雾霾的干扰；同时，是彩色rgb图像，或是全色的灰度图像。

所述的步骤2中对输入的光学遥感图像类型进行判断的依据为：获取输入光学遥感图像的三个参数，即图像纵向像素个数m、横向像素个数n和图像维数k，如果图像维数k大于1，表明该图像是rgb彩色图像；否则该图像就是全色灰度图像。

所述的步骤3中对全色灰度图像转换成rgb彩色图像处理为；灰度图像分别赋予每个通道，这样就得到一幅看似像灰度图像的rgb彩色图像，如果改变赋值大小，就会得到不同颜色组合的彩色图像；按式(1)进行处理，得到彩色图像，

式中i的表示灰度图像，r、g和b分别表示彩色图像的三个通道图像。

所述的步骤4中把rgb图像转换成hsi图像的具体步骤为：

步骤4.1：获取hsi图像的h分量图，按式(2)来计算和获取，

其中按式(3)计算，

步骤4.2：hsi图像的s分量图按式(4)来获得，

步骤4.3：hsi图像的i分量图按式(5)来获得，

所述的步骤5中对i分量图进行直方图规定化处理为：直方图反映不同灰度级出现的频率或次数，用公式(6)表示，

其中，p(rk)表示灰度值为rk的第k级灰度级出现的频率，nk表示灰度值为rk的像素点出现的个数的总和，n表示图像中所有像素的总和，即总的像素的个数；通过对图像直方图的调整，能够调整图像灰度级的空间分布情况；采用直方图规定化，而且规定化的灰度等级设为1024。

所述的步骤6中对i分量图进行灰度对比度范围调整为：采用线性的调整策略，用下式表示，

式中f(m,n)表示变化前的灰度值，且f(m,n)∈[a,b]，g(m,n)表示变化后的灰度值，且g(m,n)∈[c,d]；这里处理的图像都进行了归一化处理，关键是参数a、b、c和d的取值，对于没有受到雾霾影响或者只受到轻雾影响的图像，这四个参数分别设为a＝0.1，b＝0.9，c＝0.3，d＝1.0；对于受到重度雾霾影响的图像，则它们的设置分别为a＝0.1，b＝0.6，c＝0.5，d＝1.0。

步骤7中所述的对i分量图进行雾霾去除处理为：

图像雾霾的去除有模型方法和非模型方法，用来描述大气散射中成像过程的物理模型如式(8)所示，

i(m,n)＝j(m,n)·t+(1-t)·a(8)

式中的(m,n)表示每个像素的空间位置，i(m,n)表示最终进入成像系统的光照强度，即该像素点获取能量的幅度的大小，也即像素点的亮度，也是实际获取的图像；j(m,n)表示没有雾霾影响的条件获得的图像，t表示透射率，a表示大气散射中光的强度；其中透射率由下式计算，

t＝e^-βd(x)(9)

这里的参数β表示大气中介质的消光系数，在均匀介质中其值通常为常数，d(x)表示相机与景物之间的距离；

对式(8)进行调整，得到式(10)，

从式(10)中看出，只要知道知道参数t和a，就能够获得去雾后的图像j(m,n)；理想情况是t→1，即优的天气条件下获取图像，此时j(m,n)≈i(m,n)；因为对图像进行了归一化处理，参数t近似的用max{i(m,n)}值进行代替，能够显著提高运算速度。

所述的步骤9中把hsi图像转换成rgb图像的具体步骤为：

步骤9.1：图像都进行了归一化处理，所以各图像的灰度值取值范围是[0,1]；由h、s、i这三个分量来获取r、g、b三个分量计算公式如式(11)、式(12)或式(13)；

步骤9.2：若h在大于等于0°到小于120°的范围内，则r、g、b分量由式(11)获取；

步骤9.3：若h在大于等于120°到小于240°的范围内，则r、g、b分量由式(12)进行计算，

步骤9.4：若h在大于等于240°到小于360°的范围内，则r、g、b分量由式(13)进行计算，

所述的步骤11中由rgb图像转换成灰度图像为：直接把rgb图像的b分量作为处理后的灰度图像。

本发明与现有技术相比的优点在于：

(1)处理的对象广泛，既可以对灰度图像进行处理，也可以彩色图像进行处理；既能去除图像雾霾降低其影响程度，又能对没有受到雾霾影响的图像进行改善，提高其视觉效果和清晰度，表明本方法有好的应用范围和好的普适性；在进行雾霾处理时，简化了暗原色先验算法处理的过程，降低了处理时间，提高了处理速度。对于灰度图像的处理，最后由rgb图像转换成灰度图像时，不是通过平均计算，而是直接提取b分量作为最终的灰度图像，这也同样简化了计算。

(2)为了改善图像的质量，不影响图像的色调，对图像进行了模型转换处理，即rgb模型与hsi模型之间的相互转换，达到实现对图像去噪、去雾霾或提高质量处理的目的，但原图像的颜色不失真。

(3)采用直方图特征和灰度对比特征对图像进行双重增强处理，充分考虑和利用像素的灰度空间分布关系。

附图说明

图1为本发明方法的流程图。

图2为灰度图像的去雾处理实验结果，图a为去雾前图像；图b为去雾后图像。

图3为轻雾图像的不同方法比较实验结果，图(a)为原始图像；图(b)为暗原色法；图(c)为retinex法；图(d)为本发明方法。

图4为浓雾图像的不同方法比较实验结果，图(a)为原始图像；图(b)为暗原色法；图(c)为retinex法；图(d)为本发明方法。

具体实施方式

下面结合附图和实施例对本发明作详细描述。

参照图1，一种基于模型与直方图及灰度增强的图像雾霾去除方法，包括以下步骤：

步骤1：输入光学遥感图像，光学遥感图像受到雾霾影响，或没有受到雾霾的干扰；同时，是彩色rgb图像，或是全色的灰度图像，它们都不会影响本方法的性能；

步骤2：判断输入光学遥感图像的类型，是灰度图像还是rgb彩色图像，判断的依据为：获取输入光学遥感图像的三个参数，即图像纵向像素个数m、横向像素个数n和图像维数k，如果图像维数k大于1，表明该图像是rgb彩色图像；否则该图像就是全色灰度图像；

步骤3：对于全色的灰度图像，需转换成rgb图像；具体为：rgb彩色图像包含有三个通道，这三个通道分别赋予红色、绿色和蓝色，它们再组合成彩色图像；同样的原理，当灰度图像转换成彩色图像时，灰度图像分别赋予每个通道，这样就得到一幅看似像灰度图像的rgb彩色图像，如果改变赋值大小，就会得到不同颜色组合的彩色图像；本发明中，为了后续处理和恢复，按式(1)进行处理，得到彩色图像，

式中i的表示灰度图像，r、g和b分别表示彩色图像的三个通道图像；

步骤4：把rgb图像转换成hsi图像，提取hsi图像的三个分量，即h分量图、s分量图和i分量图；具体步骤为：

步骤4.1：获取hsi图像的h分量图，按式(2)来计算和获取，

其中按式(3)计算，

步骤4.2：hsi图像的s分量图按式(4)来获得，

步骤4.3：hsi图像的i分量图按式(5)来获得，

步骤5：对i分量图进行直方图规定化处理；直方图是数字图像的最基本统计特性，反映不同灰度级出现的频率或次数，用公式(6)表示，

其中，p(rk)表示灰度值为rk的第k级灰度级出现的频率，nk表示灰度值为rk的像素点出现的个数的总和，n表示图像中所有像素的总和，即总的像素的个数；通过对图像直方图的调整，能够调整图像灰度级的空间分布情况；图像直方图的处理通常包含两种方法：均衡化处理和规定化处理。均衡化是一种全局的处理方法，规定化是一种局部或全局的方法，因此，均衡化是规定化的一种特例；在本发明中采用直方图规定化，而且规定化的灰度等级设为1024；

步骤6：对i分量图进行灰度对比度范围调整：对灰度对比度进行调整处理，是本发明的一个重要步骤，因为不同的调整范围和不同的方式，都会导致不同的结果，可以是线性的调整，或是非线性的调整，在本发明中，采用线性的调整策略，用下式表示，

式中f(m,n)表示变化前的灰度值，且f(m,n)∈[a,b]，g(m,n)表示变化后的灰度值，且g(m,n)∈[c,d]；这里处理的图像都进行了归一化处理，目的是减少光照强度不均造成的影响；关键是参数a、b、c和d的取值，对于没有受到雾霾影响或者只受到轻雾影响的图像，这四个参数分别设为a＝0.1，b＝0.9，c＝0.3，d＝1.0；对于受到重度雾霾影响的图像，则它们的设置分别为a＝0.1，b＝0.6，c＝0.5，d＝1.0；

步骤7：对i分量图进行雾霾去除处理：雾霾、灰尘和烟等因素对图像的获取产生较大的影响，会使获得的图像质量下降；主要原因是大气中的水汽、浮尘和烟尘颗粒对太阳光波起到散射作用，从而使在雾霾天气中获取的光学遥感图像都比较模糊，尤其是较远的部分；由于雾霾的存在，所以光在大气中的穿透率减少，到达目标的能量小；但在目标周围，由于雾霾存在，增强了散射作用，因此一幅包含雾霾的图像比无雾霾的图像会显得亮一些，即偏白一些，但模糊度会显著增强；因此，图像雾霾的去除有模型方法和非模型方法，用来描述大气散射中成像过程的物理模型如式(8)所示，

i(m,n)＝j(m,n)·t+(1-t)·a(8)

式中的(m,n)表示每个像素的空间位置，i(m,n)表示最终进入成像系统的光照强度，即该像素点获取能量的幅度的大小，也即通常说的像素点的亮度，也是实际获取的图像；j(m,n)表示没有雾霾影响的条件获得的图像，t表示透射率，a表示大气散射中光的强度；其中透射率由下式计算，

t＝e^-βd(x)(9)

这里的参数β表示大气中介质的消光系数，在均匀介质中其值通常为常数，d(x)表示相机与景物之间的距离；

对式(8)进行调整，得到式(10)，

从式(10)中看出，只要知道知道参数t和a，就能够获得去雾后的图像j(m,n)；理想情况是t→1，即优的天气条件下获取图像，此时j(m,n)≈i(m,n)；因为对图像进行了归一化处理，参数t近似的用max{i(m,n)}进行代替，能够显著提高运算速度，这是本发明的一个重要改进的地方；

步骤8：对s分量图进行雾霾去除处理，与步骤7相同，只是把s分量图代替i分量图；

步骤9：把hsi图像转换成rgb图像，具体步骤为：

步骤9.1：为了便于处理，同时减少光照不均的影响，图像都进行了归一化处理，所以各图像的灰度值取值范围是[0,1]；由h、s、i这三个分量来获取r、g、b三个分量计算公式如式(11)、式(12)或式(13)，这样就实现了hsi图像到rgb图像的转换；

步骤9.2：若h在大于等于0°到小于120°的范围内，则r、g、b分量由式(11)获取；

步骤9.3：若h在大于等于120°到小于240°的范围内，则r、g、b分量由式(12)进行计算，

步骤9.4：若h在大于等于240°到小于360°的范围内，则r、g、b分量由式(13)进行计算，

步骤10：如果处理的图像是rgb图像，则处理结束，直接输出结果；

步骤11：如果处理的图像是全色灰度图像，则要由rgb图像转换成灰度图像，把其作为处理后的结果，然后直接输出；

由于对于灰度图像进行处理，是先转换成rgb图像进行处理的，处理完成后，又需要把rgb图像转换成灰度图像，这里的转换过程比较简单，不需要计算，直接把rgb图像的b分量作为处理后的灰度图像。

雾霾天气对遥感图像的影响实质上是一种典型的遥感图像质量下降现象，所以去除雾霾影响本质上是图像的恢复或复原，是一种全域的处理方式。因此，本发明从图像增强原理和雾霾产生影响机理的角度提出地一种新的图像雾霾去除算法以及图像增强算法。如果图像是全色的灰度图像，则先把灰度图像转换成rgb图像；如果已经是rgb图像，则可以直接进入下一步处理。然后把rgb图像转换成hsi图像，接着提取h、s和i分量图像，对i分量图像进行直方图规定化处理，完后对其灰度值对比范围进行调整与增强处理，最后对s分量和i分量图像分别进行雾霾的去除处理，把处理后的hsi图像再转变成rgb图像，对于灰度图像还需要提取rgb模型中b分量，作为最后的结果。

参照图2，根据图像去雾前后的对比，看出去雾后图像更清晰。

参照图3和图4，通过本发明方法和原始图、暗原色法、retinex法进行对比，方法性能评价参数值参见表1，

表1

能够得出本发明不仅能有效去除薄雾霾，而且对浓雾霾也能降低其影响。对于薄雾或没雾的图像能起到增强图像的清晰程度和视觉效果。用实测数据进行了比较实验，不但取得了好的实验效果，而且有较宽的普适性，相比与典型的传统方法，运算时间明显减少，具有极大的应用潜力和前景。