发明名称 一种智能的无人艇视觉系统视频快速海雾去除方法
摘要 本发明属于视频图像增强领域,具体涉及一种智能的无人艇视觉系统视频快速海雾去除方法。本发明包括:获取海上视频图像;判别海上场景是否含雾:去除单幅图像海雾;去除快速的视频海雾。经过本发明所提方法的处理,能够极大的提高无人艇视频去雾的智能性和实时性,进而使得无人艇或其他水面航行器对于不同水面环境的适应能力。
申请公布号 CN104933680A 申请公布日期 2015.09.23
申请号 CN201510112352.0 申请日期 2015.03.13
申请人 哈尔滨工程大学 发明人 马忠丽;刘权勇;郝亮亮;文杰;顾衍明;张永奇
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 一种智能的无人艇视觉系统视频快速海雾去除方法,其特征在于:(1)获取海上视频图像:在海雾下,通过安装在无人艇或其他水面航行器上的可见光摄像头采集海雾视频图像;(2)判别海上场景是否含雾:通过实际运行过程中采集的首帧图像来识别当前场景是否含雾:(2.1)建立海上图像库:包括水面有雾图像和清晰图像样本库;(2.2)提取图像特征:提取图像改进后的均值、图像的能见度和暗通道图像强度3个特征作为水面有雾图像和清晰图像的分类特征:(2.2.1)改进的均值:图像均值表示为:<img file="FDA0000681843160000011.GIF" wi="212" he="143" />其中p<sub>i</sub>同样表示在一幅图像中灰度值为i的像素点出现的概率,改进的均值特征公式:μ′=μ/(M‑m)其中μ为所求的均值,M为灰度图像的最大灰度值,m为灰度图像的最小灰度值;(2.2.2)图像能见度:图像能见度:对于一幅灰度图像f的任意一个像素点x,若在以x为中心的4邻域内满足:min[f(x),f(y)]≤s<max[f(x),f(y)]则称图像f中的像素对(x,y)能被阈值s分割,其中y为以x为中心的4邻域像素点,即y∈V<sub>4</sub>(x);F(s)为图像f中所有能被阈值s分割的像素对(x,y)所组成的集合,对于每一个分割阈值s∈[0,255]都对应一个集合F(s),对于集合F(s)中的任意像素对(x,y),图像空间对比度为:<img file="FDA0000681843160000012.GIF" wi="867" he="153" />则集合F(s)的平均对比度为:<img file="FDA0000681843160000021.GIF" wi="636" he="136" />其中cardF(s)表示集合F(s)的像素对总数,得到最佳的分割阈值s<sub>0</sub>:<img file="FDA0000681843160000022.GIF" wi="399" he="93" />即将平均对比度取最大值时对应的阈值s<sub>0</sub>作为图像的分割阈值,在原灰度图像f中满足空间对比度2C<sub>x</sub>(s<sub>0</sub>)>0.05的像素点则为可见边缘像素点,得到图像能见度:<img file="FDA0000681843160000023.GIF" wi="412" he="134" />其中n为可见边缘像素点总数,m为图像中满足条件C<sub>x</sub>(s<sub>0</sub>)>0的所有像素点总数,<img file="FDA0000681843160000026.GIF" wi="45" he="60" />为可见边缘像素点组成的集合;(2.2.3)暗通道图像强度:根据前一节所求的可见边缘图像求得对应原图像中的可见边缘暗通道图像,统计暗通道图像中像素值小于阈值T的像素点总数以及可见边缘图像的像素总数,得到图像强度为:<img file="FDA0000681843160000024.GIF" wi="446" he="150" />其中I(i)为由可见边缘图像得到的暗通道灰度图像中像素值为i的像素点数量;其中暗通道图像为:<img file="FDA0000681843160000025.GIF" wi="619" he="97" />J<sup>c</sup>(x)表示彩色图像J的RGB三通道中的某一通道,Ω(x)表示以x为中心的方形窗口;得到的暗通道图像I:I(x)=J<sup>dark</sup>(x)·I<sub>e</sub>(x)其中I<sub>e</sub>为前面求取能见度时所求的可见边缘图像,该边缘图像在图像目标边缘处的值为1,显示为白色,非边缘处的值为0;(2.3)识别BP神经网络雾图:提取水面图像样本库中所有有雾图像和清晰图像的三个特征,利用BP神经网络进行分类器训练,得到当前周围水面场景是否含雾的识别结果,若周围场景图像的识别结果为含雾则执行步骤(3);(3)去除单幅图像海雾:(3.1)估计大气光值:若当周围场景有雾图像为I,首先对图像I进行最小值滤波,得到其暗通道图像J<sup>dark</sup>;然后取经过最小值滤波后的图像像素值最大的0.01%个像素值作为候选大气光值估计区域;对候选像素点进行8联通区域分析;再选取连通分量像素个数最多的连通分量作为大气光值估计的精确候选区域;最后在区域中选取最大的亮度值作为天空亮度即大气光值A的估计;(3.2)透射率的估计和优化:利用得到的暗通道图像J<sup>dark</sup>以及大气光值A,初步估计当前无人艇场景图像的透射率:<img file="FDA0000681843160000031.GIF" wi="420" he="83" />将原有雾图像的灰度图像作为引导滤波的引导图像,得到原有雾图像的灰度图像:I<sup>gray</sup>(x)=0.3×I<sup>r</sup>(x)+0.59×I<sup>g</sup>(x)+0.11×I<sup>b</sup>(x)I<sup>r</sup>,I<sup>g</sup>,I<sup>b</sup>分别表示原无人艇采集的海雾图像的RGB三个通道分量图像,利用引导滤波对透射率进行优化:<img file="FDA0000681843160000032.GIF" wi="386" he="95" /><img file="FDA0000681843160000033.GIF" wi="41" he="62" />为初步估计的透射率;t表示引导滤波优化后的透射率;W为引导滤波核函数:<img file="FDA0000681843160000034.GIF" wi="1045" he="176" />其中μ<sub>k</sub>和σ<sub>k</sub>分别为引导图像I<sup>gray</sup>中局部小区域ω<sub>k</sub>内的均值和方差,而|ω|表示该局部小区域内的像素总数;ε是调整因子;(3.3)复原单幅海雾图像:求解大气散射物理模型,得到海雾去除后的图像:<img file="FDA0000681843160000035.GIF" wi="501" he="142" />A为前面改进的大气光值估计算法所求的大气光值,t为引导滤波优化后的透射率;t<sub>0</sub>为防止分母为0的因子,J为经过本发明单幅图像海雾去除算法处理后的清晰图像,I为原海雾图像;(4)去除快速的视频海雾:(4.1)求取雾气分量:对于无人艇采集的海雾视频,首先对海雾视频的第i帧海雾图像进行单幅图像海雾去除清晰化处理,然后求取该视频第i帧的初始雾气分量:F<sub>i</sub>=I<sub>i</sub>‑J<sub>i</sub>I<sub>i</sub>表示当前海雾视频下的第i帧海雾图像,J<sub>i</sub>是对I<sub>i</sub>进行处理后得到的第i帧去雾清晰化图像,F<sub>i</sub>为视频第i帧的雾气分量;(4.2)引导滤波的雾气分量估计:估计海雾视频背景变化量:(4.2.1)将视频第i帧和第j帧图像进行做差运算,并进行二值化操作:<img file="FDA0000681843160000041.GIF" wi="1055" he="202" />上式中I<sup>gray</sup>(x,y,t<sub>j</sub>)表示海雾下视频第j帧图像在坐标(x,y)处的像素灰度值,I<sup>gray</sup>(x,y,t<sub>i</sub>)表示的是视频第i帧海雾图像在坐标(x,y)位置的灰度值,t<sub>j</sub>、t<sub>i</sub>分别表示视频的第j帧和第i帧,Th表示二值化操作的阈值,得到的二值图像D中目标运动或背景变化大的区域为白色,值为255,而变化不大的区域则为黑色,值为0;(4.2.2)求取海雾视频第i帧和第j帧之间的变化量:<img file="FDA0000681843160000043.GIF" wi="230" he="154" /><img file="FDA0000681843160000047.GIF" wi="105" he="84" />是视频帧灰度图像的像素总数,即采集的图片的尺寸,N<sub>D</sub>得到的二值图D中值为255的像素总数,即运动目标区域或背景变化大的区域;(4.2.3)若T<sub>D</sub>>T,则可认为视频第j帧相对于第i帧场景变化大,即场景中的雾气浓度变化较大;否则认为背景变化不大,即场景中雾气浓度变化不大,T为视频去雾时雾气分量的更新阈值;(4.3)去除视频海雾:若海雾视频第j帧与第i帧背景变化大,对第j帧进行单幅图像海雾去除处理;否则表明场景雾气浓度变化不大,利用引导滤波得到的第i帧雾气分量进行滤波,得到第j帧海雾视频的雾气分量<img file="FDA0000681843160000044.GIF" wi="86" he="90" />此时引导滤波的引导图像为当前要去雾的第j帧视频海雾图像;用引导滤波得到的第j帧视频海雾图像的雾气分量<img file="FDA0000681843160000045.GIF" wi="85" he="94" />得到第j帧视频海雾图像的去雾恢复结果:<img file="FDA0000681843160000046.GIF" wi="249" he="92" />I<sub>j</sub>表示的是原始第j帧海雾视频图像,J<sub>j</sub>表示对第j帧海雾视频图像去雾恢复后的图像。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室