《暗原色先验去雾算法的一些认识》由会员分享，可在线阅读，更多相关《暗原色先验去雾算法的一些认识（7页珍藏版）》请在金锄头文库上搜索。
1、暗原色先验去雾算法的一些认识摘要：本文首先简单介绍基于暗原色先验去雾算法，然后总结出各变量的求解，最后利用已 知程序检验算法的效果。引言：在雾天时，由于空气中大量悬浮粒子的散射和吸收作用，使得光学成像对比度降低， 颜色偏灰白色，严重影响视觉效果，这直接限制室外各种视觉系统的作用。因此对图像去雾 有着重要的现实意义。暗原色先验是基于户外无雾图像数据库得到的一条简单有效的图像先 验规律。该理论的核心是：绝大多数户外无雾图像的每个局部区域里，至少存在一个颜色通 道的强度值很低。称之为在暗像素。雾天中，这些暗像素会因为雾光的充斥而强度值变大。因此这些暗像素能够用来评估雾光的透射情况。根据McCarne
2、y的大气散射模型，根据光在雾天传输的物理特性可以得到较好的去雾原像。1. 雾天图像去雾模型McCarney 的大气散射模型(如下)被广泛应用于计算机视觉和图形学领域中I(x) = t(x)J(x) + (l-t(x)A (1)I是输入图像(观测到的)光强度，t指光线透射率，A是大气光成分，J是无雾时景物的光 线强度。图像去雾目标就是由已知的I求得未知参数J,A,t。由于方程的个数少于未知量的个 数，需增加约束条件来求解。暗原色先验即是其中的一种约束。2. 暗原色先验暗原色先验是基于大量户外无雾图像观察到的一条统计规律：在绝大数户外无雾图像的每个局部区域至少存在某个颜色通道的强度值很低。对户外无
3、雾图像 J 进行分块，对每个像素 块定义暗原色如下：J dark (x)二 min (min (J (y) (2)ccer,g,b yeQ(x)式中：Q (x)是以x为中心的正方形邻域，J为J三原色的一个通道，Jdark (x)即为图像Jc在这个邻域的暗原色。观察统计表明Jdarl于零。在带雾图像中，这些暗原色的值会升高。3. 模型求解对于t的估计需要以下的假设：每一个像素块的光线透射率是相同的。因此像素块的大小直 接影响透射率的估计。像素块较小则t的细节多，层次感好。但会增加运算量，且可能增加 局部错误；快较多时，得到的暗原色图分布会很单一，图片的层次感不明显，会使对比度下 降。因此分块时可
4、尝试不同的值，以期得到较好的值。假设大气光A已知(事实上t的估计 不需要A)，首先在模型(1)两端在每个像素块求取最小值：min( min (1 (y) = t(x)min( min (Jc(y) + (1-1(x) (3)yeQ( x) Ac yeQ( x) Acc根据暗原色先验无雾图像的暗原色趋于零，故由上式可得I ( y )t(x) = 1 一min( min ) (4)c yeQ( x) Ac而右式第二项即为带雾图像在x邻域的暗原色值，这是可以通过matlab得到的，因此局部 区域的t值可以求得。从而可以得到整幅图的透射率t(x)。但当t(x)趋于零时，图像趋于包含噪音，因此常常需要设
5、置一个透射率下限10 (通常设置为0.1)。大气光的估计：在一般的单一图像去雾中，通常整幅图最大强度的像素作为大气光的估测 最后得到 J 的求解公式为I (x) AJ(x) + A (5)max(t(x),t )4.实验结果与分析 利用已知的程序来验证算法的有效性，可以看到原本模糊的图像经处理后变得更加清晰，色 彩更鲜明，视觉上感觉更舒适。原图像灰度图土雯三加底園Original ImageHaze-Free Image:总结:暗原色先验使得单一图像去雾变得简洁且有效。但由于本身是一条统计规律，因此在 某些极端情况下（如全白色的图景）效果就会受到很大的影响。在实际应用中，暗原色先验 原理还可以
6、推广到音影材料的去噪中等等，这都可以进一步探索。参考文献：1.He Kai-Ming,SUN JIAN,TANG Xiao-ou.Single image haze removed using dark channel priorA.Proceedings of IEEEConference on computer Vision and Pattern RecognitionC.Mami,2009,1956-1963.2嵇晓强，戴明，孙丽娜，郎小龙，王洪.暗原色先验图像去雾算法研究.光电子，1005-0086 （2011）06-0926-053.胡伟，袁国栋，董朝，疏学明.基于暗通道优先的单幅
7、去雾新方法.计算机研究与发展，2010附录：(已知程序) close all clear allclcblockSize=8;%每个 block 为 8 个像素w0=0.6;t0=0.1;% A=200;I=imread(C:UsersasusDesktop 图片 17.jpg); h = figure;subplot(321) imshow(I);title(Original Image);subplot(323); grayI=rgb2gray(I);imshow(grayI,);title(原图像灰度图)原图像的灰度图subplot(324); imhist(grayI,64);%统计0
8、.02)error(This image need not Haze-Free-Proprocessing.); else if(percent0.01) w=0.3;else w=0.45;endendendh,w,s=size(I);min_I=zeros(h,w);for i=1:hfor j=1:wdark_I(i,j)=min(I(i,j,:);%取每个点的像素为RGB分量中最低的那个通道的值 endend Max_dark_channel=double(max(max(dark_I) dark_channel=double(dark_I);t=1-w0*(dark_channel/
9、Max_dark_channel);T=uint8(t*255);t=max(t,t0);%I1=double(I);J(:,:,1) = uint8(I1(:,:,1) - (1-t)*Max_dark_channel)./t);J(:,:,2) = uint8(I1(:,:,2) - (1-t)*Max_dark_channel)./t);J(:,:,3) =uint8(I1(:,:,3) - (1-t)*Max_dark_channel)./t);subplot(322)imshow(J);imwrite(J,tj2.jpg);title(Haze-Free Image:);subplot(325);grayJ=rgb2gray(J);imshow(grayJ,);title(去雾后灰度图)
