专利摘要

本发明涉及图像处理技术领域，提供一种雾霾图像中大气光值的提取方法及系统，方法包括：对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A，实现对大气光值的快速计算，适应不同大小和天空区域的图像，为雾霾图像快速去雾处理提供条件。

说明书

技术领域

本发明属于图像处理技术领域，尤其涉及一种雾霾图像中大气光值的提取方法及系统。

背景技术

雾霾场景中，大气粒子的散射作用会导致光学传感器捕获到的信息严重降质，图像中对比度和色彩保真度方面都会发生不同程度的衰减，直接影响到人类视觉的感受和机器视觉系统的正常工作，因此研究图像去雾方法具有重要意义。基于图像复原的方法针对雾天图像降质的物理过程建立雾天图像退化模型，通过反演运算补偿有雾图像中丢失的信息。该方法针对性强，得到的去雾效果自然，因而得到了大量关注并成为研究的热点。

根据大气散射理论，景物在雾天成像时的模型主要包括两部分：场景中物体反射的光在到达传感器的过程中受大气粒子的散射而衰减的过程和太阳光被大气中的悬浮粒子散射后到达传感器的过程。因此，在计算机视觉和计算机图形学中，描述有雾图像的散射模型可表示为：

I(x)＝J(x)t(x)+A(1-t(x))；

x为空间坐标，I(x)为有雾图像，J(x)为场景辐射照度或清晰无雾的图像，A为整体大气光值，t(x)为大气透射率函数(参数)。一旦求出透射率和大气光值，则可以恢复场景深度：

但是，目前对于雾霾图像的处理过程中，大气光值的估计过程较为复杂，从而影响了整个图像去雾的效率。

发明内容

本发明的目的在于提供一种雾霾图像中大气光值的提取方法，旨在解决现有技术中大气光值的估计过程较为复杂，从而影响了整个图像去雾的效率的问题。

本发明是这样实现的，一种雾霾图像中大气光值的提取方法，所述方法包括下述步骤：

对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；

对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；

对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；

提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A。

作为一种改进的方案，所述对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)的步骤具体包括下述步骤：

对获取到的最小值滤波图像M(x,y)进行水平投影，获取一列所有像素点的灰度值的累加参数H(y)；

计算水平投影的宽度b_H，并以此确定天空区域的自适应宽度值为2b_H+1；

对投影在2b_H+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像K(x,y)。

作为一种改进的方案，所述对获取到的最小值滤波图像M(x,y)进行水平投影，获取一列所有像素点的灰度值的累加参数H(y)的步骤之后，所述计算水平投影的宽度b_H的步骤之前还包括下述步骤：

对累加参数H(y)进行平滑滤波，消除图像中存在的尖峰值和奇异值。

作为一种改进的方案，所述对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)的步骤具体包括下述步骤：

对最大值区域图像K(x,y)进行垂直投影，获取一行所有像素点的灰度值的累加参数V(x)；

计算垂直投影的宽度b_V，并以此确定天空区域的自适应宽度值为2b_V+1；

对投影在2b_V+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像R(x,y)。

作为一种改进的方案，所述对最大值区域图像K(x,y)进行垂直投影，获取一行所有像素点的灰度值的累加参数V(x)的步骤之后，所述计算垂直投影的宽度b_V的步骤之前还包括下述步骤：

对累加参数V(x)进行平滑滤波，消除图像中存在的尖峰值和奇异值。

本发明的另一目的在于提供一种雾霾图像中大气光值的提取系统，所述系统包括：

最小值滤波模块，用于对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；

第一最大值区域图像获取模块，用于对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；

第二最大值区域图像获取模块，用于对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；

大气光值确定模块，用于提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A。

作为一种改进的方案，所述第一最大值区域图像获取模块具体包括：

水平投影模块，用于对获取到的最小值滤波图像M(x,y)进行水平投影，获取一列所有像素点的灰度值的累加参数H(y)；

水平投影宽度计算模块，用于计算水平投影的宽度b_H，并以此确定天空区域的自适应宽度值为2b_H+1；

第一图像确定模块，用于对投影在2b_H+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像K(x,y)。

作为一种改进的方案，所述系统还包括：

第一平滑滤波模块，用于对累加参数H(y)进行平滑滤波，消除图像中存在的尖峰值和奇异值。

作为一种改进的方案，所述第二最大值区域图像获取模块具体包括：

垂直投影模块，用于对最大值区域图像K(x,y)进行垂直投影，获取一行所有像素点的灰度值的累加参数V(x)；

垂直投影宽度计算模块，用于计算垂直投影的宽度b_V，并以此确定天空区域的自适应宽度值为2b_V+1；

第二图像确定模块，用于对投影在2b_V+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像R(x,y)。

作为一种改进的方案，所述系统还包括：

第二平滑滤波模块，用于对累加参数V(x)进行平滑滤波，消除图像中存在的尖峰值和奇异值。

在本发明实施例中，对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A，实现对大气光值的快速计算，适应不同大小和天空区域的图像，为雾霾图像快速去雾处理提供条件。

附图说明

图1是本发明提供的雾霾图像中大气光值的提取方法的实现流程图；

图2是本发明提供的对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)的实现流程图；

图3是本发明提供的对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)的实现流程图；

图4是本发明提供的雾霾图像中大气光值的提取系统的结构框图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

图1示出了本发明提供的雾霾图像中大气光值的提取方法的实现流程图，其具体包括下述步骤：

在步骤S101中，对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)。

其中， 其中，(x，y)为图像中的像素点，c为彩色图像的三色通道。

在步骤S102中，对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)。

在步骤S103中，对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)。

在步骤S104中，提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A。

在该步骤中，将属于填空区域的像素值提取出来，然后对该填空区域的像素值进行降序排列，选取亮度大于一定数值的像素的平均灰度值作为大气光值，该一定数量可以是10％，即A＝mean(max^0.1R(x,y))

图2示出了本发明提供的对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)的实现流程图，其具体包括下述步骤：

在步骤S201中，对获取到的最小值滤波图像M(x,y)进行水平投影，获取一列所有像素点的灰度值的累加参数H(y)。

其中，水平投影的公式表达为：

其中，1≤x≤n。

在步骤S202中，计算水平投影的宽度b_H，并以此确定天空区域的自适应宽度值为2b_H+1。

其中，该水平投影的宽度b_H的计算方式有多种，下述给出其中一种，在此不再赘述。

在步骤S203中，对投影在2b_H+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像K(x,y)。

在该步骤中，最大值区域： 其中，b_H+1≤t≤n-b_H；裁剪得到的最大值区域图像保存为K(x,y)，其大小为m×(2b_H+1)像素。

其中，窗口宽度为2b_H+1且为奇数，t为2b_H+1宽度内的任一一点。

其中，H(y′)为平滑滤波后的累加参数，即：

对步骤S201得到的累加参数H(y)进行平滑滤波，消除图像中存在的尖峰值和奇异值，其中：

其中，y＝σ+1,σ+2,...,m-σ；

其平滑效果随着滤波器宽度的增大而变的更加平滑，一般设置σ＝2，即滤波器的宽度为2σ+1，共5个像素。

进一步地，由于输入图像的尺寸不同，而且图像中天空区域在整幅图像中所占的比例也不同，如果天空区域所占的比例较大，则H_max/H_m所占的比例较小，可以选择较宽的水平投影宽度，反之则选择较小固定值，为了适应不同尺寸和天空区域的图像，确定水平投影的自适应宽度值为2b_H+1，其计算方式为：

$b_H=\left\{\begin{array}{cc}\left(H_{\max }/H_m\right)\×\left(w/40\right)& \left(H_{\max }/H_m\right)\≤1.5\\ w/40& \left(H_{\max }/H_m\right)\≥1.5\end{array}\right.;$

H_m＝mean(H′(y))，H_max＝max(H′(y))，其中，H_m为H′(y)序列的平均值，mean为求平均运算，H_max是H′(y)序列的最大值，max为求最大值运算。

通过该方法能够更好的处理不同大小和不同天空面积的图像。

图3示出了本发明提供的对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)的实现流程图，其具体包括下述步骤：

在步骤S301中，对最大值区域图像K(x,y)进行垂直投影，获取一行所有像素点的灰度值的累加参数V(x)。

其中，1≤y≤m。

在步骤S302中，计算垂直投影的宽度b_V，并以此确定天空区域的自适应宽度值为2b_V+1。

该计算方式如上述方式相同，在此不再赘述。

在步骤S303中，对投影在2b_V+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像R(x,y)。

该最大值区域： 其中，b_V+1≤t≤n-b_V；

其中，图像R(x,y)的大小为(2b_V+1)×(2b_H+1)。

在本发明实施例中，在对图像进行投影后，需要计算一定区域内的数据和，相邻的区域存在重复现象，为了提高计算精度，可以采用移位求和的加速算法，具体为：

假设数列i(x)，1≤x≤n，区域宽度为(2b+1)，则区域内的求和为：

b+1≤t≤n-b，t为透射率；

由于相邻区域Sum(t)和Sum(t-1)存在{i(t-b+1),i(t-b+2),...,i(t+b-1)}共2b-1个重复点，因此，可以直接利用前面区域进行计算，即：

Sum(t)＝Sum(t-1)-i(t-b)+i(t+b)；

通过该加速，可以将2b+1次求和运算减小到3次运算，大大提高了运算速度。

图4是本发明提供的雾霾图像中大气光值的提取系统的结构框图，为了便于说明，图中仅给出了与本发明实施例相关的部分。

最小值滤波模块11对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；第一最大值区域图像获取模块12对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；第二最大值区域图像获取模块13对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；大气光值确定模块14提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A。

其中，第一最大值区域图像获取模块12具体包括：

水平投影模块15对获取到的最小值滤波图像M(x,y)进行水平投影，获取一列所有像素点的灰度值的累加参数H(y)；

水平投影宽度计算模块16计算水平投影的宽度b_H，并以此确定天空区域的自适应宽度值为2b_H+1；

第一图像确定模块17对投影在2b_H+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像K(x,y)。

在该实施例中，第一平滑滤波模块18对累加参数H(y)进行平滑滤波，消除图像中存在的尖峰值和奇异值。

在本发明实施例中，第二最大值区域图像获取模块13具体包括：

垂直投影模块19对最大值区域图像K(x,y)进行垂直投影，获取一行所有像素点的灰度值的累加参数V(x)；

垂直投影宽度计算模块20计算垂直投影的宽度b_V，并以此确定天空区域的自适应宽度值为2b_V+1；

第二图像确定模块21对投影在2b_V+1宽度区域内像素值进行求和运算，并筛选出最大值区域，确定筛选出的最大值区域所对应的图像R(x,y)。

其中，第二平滑滤波模块22对累加参数V(x)进行平滑滤波，消除图像中存在的尖峰值和奇异值.

其中，上述各个模块的具体实现如上述方法实施例所记载，在此不再赘述。

在本发明实施例中，对输入的图像I(x,y)中的任意像素三通道进行最小值滤波，得到最小值滤波图像M(x,y)；对获取到的最小值滤波图像M(x,y)进行水平投影，并获取水平投影后的最大值区域图像K(x,y)；对所述最大值区域图像K(x,y)进行垂直投影，并获取垂直投影后的最大值区域图像R(x,y)；提取最大值区域图像R(x,y)中的像素值，选取一定数量的像素值的平均值作为大气光值A，实现对大气光值的快速计算，适应不同大小和天空区域的图像，为雾霾图像快速去雾处理提供条件。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

一种雾霾图像中大气光值的提取方法及系统专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。