基于阈值分割优化的暗原色天空区域去雾方法
【技术领域】
[0001] 本发明属于图像处理技术领域，更进一步的设及图像去雾方法，可用于图像增强 和视频监控。
【背景技术】
[0002] 图像去雾的目的在于去除图像中的雾信息，提高图像的对比度，提高图像的质量。 由于天气等因素的干扰，在图像的采集过程中可能会受到雾靈、烟雾等因素的影响从而使 得图像的对比度下降，降低图像的质量，影响图像的后续处理及其使用价值。因此，图像去 雾具有重要的意义。
[0003] 巧义斌，李鹏，孙锦涛等发表的论文"一种交互式雾天图像去雾方法"中提出一种 基于3日1:;[]1糾（1?日1:;[]1日1-。〇1'1日义-1:11日〇巧）理论的多尺度1?日1:;[]1糾（151〇去雾算法。1?日1:;[]16又理 论是基于色彩恒常的一种代表性计算理论，其全称是视网膜皮层理论(Retinal-cortex- theory)。Retinex理论认为人类感知到的物体颜色与物体表面的反射性质有着密切的关 系，而与投射到人眼的光的特性关系不大。该去雾方法的具体步骤是:首先，利用取对数的 方式将原图像的照射分量和反射分量进行分离，并与高斯滤波函数进行卷积，得到低通滤 波后的图像;之后用原图像减去低通滤波后的图像，得到高频增强的图像;最后对该图像进 行取反对数并进行对比度增强，得到最终的结果图像。该方法虽然能够在一定程度上将图 像的雾信息去除，但只是从图像增强上对图像进行去雾，没有考虑到图像中雾的本质，因此 不能将图像中的雾信息去除干净，效果不佳。
[0004] 何凯明等人在2009年发表的论文"Single Image Haze Removal Using Dark Channel Prior"提出基于暗原色先验(Dark Channel Prior,DCP)的单幅图像去雾方法。该 方法认为，户外清晰的图像总符合暗原色先验理论，该理论认为在一幅户外清晰图像上，总 存在一些像素点的值很小，甚至趋近于0。该方法的具体步骤如下:首先，求取图像的暗通道 图，之后根据暗通道图求取图像的透射率和大气光值，最后根据求得的透射率和大气光值 复原图像。该方法从图像中雾的本质出发，当图像中没有高亮的天空区域时，能够很好的去 除图像中的雾信息，去雾效果良好。但当图像有大量高亮天空区域时，暗原色先验理论不再 适用于运些区域，导致复原效果不佳。
【发明内容】
[0005] 本发明的目的在于针对上述现有技术的不足，提出一种基于阔值分割优化的暗原 色天空区域去雾方法，W在使用暗原色运一理论有效去除非天空区域的雾信息的同时，优 化处理天空区域的效果，获得复原效果好的去雾图像。
[0006] 为实现本发明目的的技术方案包括如下：
[0007] (1)输入任选的一幅含有天空区域的带雾图像I，依次计算该带雾图像的暗通道值 和大气光值A;
[000引（2)计算带雾图像每个像素点的透射率，并将各像素点的透射率组合为雾天图像 的透射图t;
[0009] (3)对透射图t进行导向滤波，得到优化后的透射图t' ；
[0010] (4)求优化后的透射图t'的灰度分布直方图，并根据该直方图分布规律选取一个 阔值T，将优化后的透射图t'中像素值小于T的区域标记为天空区域；
[0011] (5)按照下式对天空区域的透射率进行放大：
[0012] t'（x)=yt'（x)，
[0013] 其中，t'（x)为像素点X的透射率，μ值为大于1的数，取值为μ=1.5;
[0014] (6)根据上述所求得的大气光值AW及像素点X的透射率t'（x)，按照下式得到最终 的去雾图像：
[0015]
[0016] 其中，to为边界约束条件，取值为to = 0.1，I(x)为带雾图像I的像素点X的值，J(x) 为最终的去雾图像J的像素点X的值。
[0017] 本发明与现有技术相比具有W下优点：
[0018] 第一，本发明通过求优化后的透射图的灰度分布直方图，并根据该直方图的分布 规律，选取一个阔值，可W准确的对天空区域进行标记；
[0019] 第二，本发明通过对天空区域的透射率进行放大，优化处理天空区域的效果，避免 天空区域出现色彩失真、光晕等现象，获得复原效果好的去雾图像。
【附图说明】
[0020] 图1是本发明的实现流程图；
[0021] 图2是本发明仿真使用的一幅户外带有部分天空区域的有雾图像，
[0022] 图3是本发明仿真使用的一幅户外带有大片天空区域的有雾图像；
[0023] 图4是用本发明与现有技术对测试图像2的实验结果对比图；
[0024] 图5是用本发明与现有技术对测试图像3的实验结果对比图。
【具体实施方式】
[0025] 下面结合附图对本发明做进一步的描述。
[00%]参照图1，本发明的具体步骤如下：
[0027] 步骤1:输入任选的一幅含有天空区域的带雾图像I，依次计算该带雾图像的暗通 道值和大气光值A:
[0028] (la)从图2和图3中任选一幅带天空区域的有雾图像作为输入的一幅带雾图像I;
[0029] (化)对带雾图像I的R，G，邮P红蓝绿Ξ个通道分别进行9X9的最小值滤波；
[0030] (1C)对带雾图像的每个像素点，取（la)中滤波后的R，G，ΒΞ个通道在该点上的最 小值，作为带雾图像在该点的暗通道值，如下式：
[0031]
[0032] 其中，Idatk(x)表示带雾图像I在像素点X的暗通道值，C表示图像的通道，Ω(Χ)表 示Wx为中屯、的9X9矩形区域。
[0033] (Id)在求得的带雾图像的暗通道图上选取前0.1 %的最亮点；
[0034] (le)在带雾图像中寻找与（Id)中选取的点位置相对应的点，再从运些对应的点中 选取亮度最大的像素点值，作为大气光值。
[0035] 步骤2:计算带雾图像每个像素点的透射率，并将各像素点的透射率组合为雾天图 像的透射图。
[0036] (2a)将雾天图像的大气散射模型表示如下：
[0037] I(x)=J(x)t(x)+A(l-t(x)), <1>
[0038] 其中，I(x)为带雾图像I的像素点X的值，J(x)为最终的去雾图像J的像素点X的值， t(x)为像素点X的透射率，A为大气光值;对式<1〉中的雾天图像及去雾图像J求其暗通道 值：
[0039]
[0040] 其中，C表示图像的通道，y表示图像中的像素点，Ω (X)表示Wx为中屯、的9X9的矩 形区域，r(y)表示图像I在像素点y处的通道C的值，r(y)表示图像J在像素点y处的通道C的 值；
[0041] (2b)根据暗原色理论，可W认为式<2〉中的，因此将式<2〉化 简为：
[0042]
[0043] 对式<3〉进行移项求得雾天图像I的各像素点的透射率：
[0044]
[0045] 步骤3:对透射图t进行导向滤波，得到优化后的透射图t'。
[0046] 步骤4:求优化后的透射图t'的灰度分布直方图，
[0047] (4a)统计优化后的透射图的所有像素点在各灰度级上的分布个数；
[0048] (4b) W灰度级为横坐标，W像素点的个数为纵坐标绘制优化后的透射图的灰度分 布直方图。
[0049] 步骤5:根据该直方图分布规律选取一个阔值T，将优化后的透射图t'中像素值小 于T的区域标记为天空区域。
[0050] 因为天空区域的灰度值比较大，不存在某一像素点的暗通道值趋近于0的