边缘保持最优化滤波方法

1、边缘保持最优化滤波方法吴文瑾摘黄敬贺要：在图像处理中，去除噪声的同时常常会给边缘与细节带来不可恢复的损失。为了解决这一问题，讲解一种可以在去噪的同时，保留边缘细节的边缘保持最优化滤波方法。首先对该算法原理与实现方式进行了阐述，之后通过实验对该滤波方法的效果加以分析，最后对该算法的优缺点进行了总结。关键词：滤波；平滑；边缘保持；EPOS1引言在摄影时，由于分辨率不足、感光度设置偏高等原因，经常会遇到一些布满椒盐噪声或斑点噪声的照片，这些细小的噪声严重影响了照片效果。在这种情况下，通常可以使用均值滤波等滤波方法对照片进行平滑处理。平滑后，一般可以达到去除噪声的效果。但是一些边缘也随之变得模糊，导致

2、照片丢失很多细节。有没有一种方法可以达到既去除噪声的目的，又不破坏边缘和细节呢？答案是肯定的。边缘保持最优化（EPOS）滤波算法就是这样一种在去除噪声的同时可以最佳保持边缘和细节的滤波方法，下面将对这种方法的原理和实现加以详细介绍。2EPOS滤波原理2.1算法综述EPOS滤波是一种可以智能调整滤波窗口大小和形状的平滑算法。它通过首先使用噪声方差的大小对均匀区域和含有边缘或孤立点的区域加以判断区分，然后仅对每个像素邻域中的均匀部分进行均值处理，来达到在滤波的同时保证边缘不变模糊的目的。在确定滤波窗口的大小和形状时，该算法将窗口的邻域分为8个不相重合的部分（如图1），通过逐步剔除与中心像素不相同的

3、邻域，应用最后剩余区域中的像素灰度统计量估算中心像元最终的灰度值。2.2算法流程与分析EPOS执行滤波的流程如图2。通过分析其计算过程可知，EPOS滤波的执行流程可以分为如下几个部分：判断是否符合均值滤波条件，窗口区域分块并计算各块方差，选出方差最大的一个并进行分块剔除，重新计算剩下分块方差。若仍不符合条件，进行窗口缩小，并重复上述过程，直至符合条件或窗口缩小为33的为止。由于每个窗口都是三角形，需要分别计算，给程序编码带来一些麻烦。综合分析，可以设置3个循环，分别为窗口移动循环、窗口缩小循环和区域删减循环，进行EPOS滤波计算。662009.23图1EPOS算法将滤波窗口划分成的8个不相重合

4、的部分图2EPOS执行滤波的流程图3程序实现/灰度值计算min=65535;floatd0=0.5227,dmax=0.0;for(i=0;iheight;i+)/窗口移动，i,j为窗口中心for(j=0;j=m0)&(j=n0)&(iheight-m0)&(j1;m-)/窗口缩小循环n=m;floatu=uu(AA,i,j,m,n,width);/窗口均值intd1=0;for(k=i-m;k=i+m;k+)for(l=j-n;l=j+n;l+)d1+=(AAk*width+l-u)*(AAk*width+l-u);d1=d1/(2*m+1)/(2*n+1);/窗口方差if(d1d0*u*u

5、)PPi*width+j=int(u+0.5);break;/判断是否符合整个区域均值滤波条件floatu08=0,d8=0,d118=0;/累加数组置零/计算各区域在未删减之前的均值和方差for(k=i-m;ki;k+)for(l=j+k-i;lj;l+)u00+=AAk*width+l;u00=u00/(2*m+1)/(2*n+1)*8;for(k=i-m;ki;k+)for(l=j+k-i;lj;l+)d110+=(AAk*width+l-u00)*(AAk*width+l-u00);d110=d110/(2*m+1)/(2*n+1)*8;/for(k=i-m;ki;k+)for(l=j

6、-k+i;lj;l+)u01+=AAk*width+l;u01=u01/(2*m+1)/(2*n+1)*8;for(k=i-m;ki;k+)for(l=j-k+i;lj;l+)d111+=(AAk*width+l-u01)*(AAk*width+l-u01);d111=d111/(2*m+1)/(2*n+1)*8;/for(k=i-m;ki;k+)for(l=j-k+i;lj+m;l+)u02+=AAk*width+l;u02=u02/(2*m+1)/(2*n+1)*8;for(k=i-m;ki;k+)for(l=j-k+i;lj+m;l+)d112+=(AAk*width+l-u02)*(A

7、Ak*width+l-u02);d112=d112/(2*m+1)/(2*n+1)*8;/for(k=i;ki+m;k+)for(l=j+k-i;lj+m;l+)u03+=AAk*width+l;u03=u03/(2*m+1)/(2*n+1)*8;for(k=i;ki+m;k+)for(l=j+k-i;lj+m;l+)d113+=(AAk*width+l-u03)*(AAk*width+l-u03);d113=d113/(2*m+1)/(2*n+1)*8;/for(k=i;ki+m;k+)for(l=j;lj+k-i;l+)u04+=AAk*width+l;u04=u04/(2*m+1)/(2

8、*n+1)*8;for(k=i;ki+m;k+)for(l=j;lj+k-i;l+)d114+=(AAk*width+l-u04)*(AAk*width+l-u04);d114=d114/(2*m+1)/(2*n+1)*8;/for(k=i;ki+m;k+)for(l=j-k+i;lj+m;l+)u05+=AAk*width+l;u05=u05/(2*m+1)/(2*n+1)*8;for(k=i;ki+m;k+)for(l=j-k+i;lj+m;l+)d115+=(AAk*width+l-u05)*(AAk*width+l-u05);d115=d115/(2*m+1)/(2*n+1)*8;/f

9、or(k=i;ki+m;k+)for(l=j-m;lj-k+i;l+)u06+=AAk*width+l;u06=u06/(2*m+1)/(2*n+1)*8;for(k=i;ki+m;k+)for(l=j-m;lj-k+i;l+)d116+=(AAk*width+l-u06)*(AAk*width+l-u06);d116=d116/(2*m+1)/(2*n+1)*8;/for(k=i-m;ki;k+)for(l=j-m;lj+k+i;l+)u07+=AAk*width+l;u07=u07/(2*m+1)/(2*n+1)*8;2009.2367for(l=j-m;lj+k+i;l+)d117+=(

10、AAk*width+l-u07)*(AAk*width+l-u07);d117=d117/(2*m+1)/(2*n+1)*8;while(z7)/执行区域删减循环/计算方差最大最小的区域for(i1=0;i18;i1+)for(j1=0;j1dmax)&(i1!=Ij1)dmax=d11i1,Imax=i1;if(d11i1min)min=d11i1;if(min=d0*u*u)break;/符合缩小窗口条件，退出区域删减循环，进入缩小窗口循环/计算删减后总的保留区域的均值Iz+=Imax;u=0;intcount=0;for(i1=0;i1z;i1+)for(j1=0;j18;j1+)if(

11、j1!=Ii1)u+=u0i1;count+;u=u/count;/计算区域删减后总的保留区域的方差d0=0;for(k=i-m;ki;k+)for(l=j+k-i;lj;l+)d0+=(AAk*width+l-u)*(AAk*width+l-u);d0=d0/(2*m+1)/(2*n+1)*8;/d1=0;for(k=i-m;ki;k+)for(l=j-k+i;lj;l+)d1+=(AAk*width+l-u)*(AAk*width+l-u);d1=d1/(2*m+1)/(2*n+1)*8;/d2=0;for(k=i-m;ki;k+)for(l=j-k+i;lj+m;l+)d2+=(AAk*

12、width+l-u)*(AAk*width+l-u);d2=d2/(2*m+1)/(2*n+1)*8;/d3=0;for(k=i;ki+m;k+)for(l=j+k-i;lj+m;l+)d3+=(AAk*width+l-u)*(AAk*width+l-u);682009.23d3=d3/(2*m+1)/(2*n+1)*8;/d4=0;for(k=i;ki+m;k+)for(l=j;lj+k-i;l+)d4+=(AAk*width+l-u)*(AAk*width+l-u);d4=d4/(2*m+1)/(2*n+1)*8;/d5=0;for(k=i;ki+m;k+)for(l=j-k+i;lj;l

13、+)d5+=(AAk*width+l-u)*(AAk*width+l-u);d5=d5/(2*m+1)/(2*n+1)*8;/d6=0;for(k=i;ki+m;k+)for(l=j-m;lj-k+i;l+)d6+=(AAk*width+l-u)*(AAk*width+l-u);d6=d6/(2*m+1)/(2*n+1)*8;/d7=0;for(k=i-m;ki;k+)for(l=j-m;lj+k-i;l+)d7+=(AAk*width+l-u)*(AAk*width+l-u);d7=d7/(2*m+1)/(2*n+1)*8;/d1=0;for(i1=0;i1z;i1+)for(j1=0;j1

14、8;j1+)if(j1!=Ii1)d1+=di1;d1=d1/count;/将达到要求的保留区域均值作为该像素的灰度值if(d1d0)PPi*width+j=int(u+0.5);break;/while区域删减循环if(d1d0)break;/同时退出缩小窗口循环/缩小窗口循环/if/窗口移动4滤波实验与效果分析对一幅实验图像（如图3）使用99窗口进行了EPOS滤波处理，结果图如图4所示，并采用同样大小的窗口做均值平滑滤波（如图5）进行效果比对。从滤波后的结果可以看出，均值滤波对原图整体进行了平滑，平滑后噪声消除效果很好，但边缘丢失情况严重，一些较小的孤立点在平滑处理后消失掉了。而在EPOS

15、滤波结果中，图像的大部分区域进行了平滑处理，边缘却被很好地保持了下来，孤立点在平滑后没有消失的现象，但部分较大的噪声在处理后非但没有被消除，反而得到了放大。图3原始实验图像图4EPOS滤波（99窗口）结果图5均值滤波（99）窗口结果通过对EPOS滤波算法的原理进行分析，可以知道，造成这种效果的原因是：EPOS滤波通过将平滑窗口分成8个区域，逐个剔除选择参与平滑的像元来进行边缘保持，8个区域的划分几乎可以包含了所有边缘方向，若边缘方向与8个区域中任一区域都无法匹配，该算法还设置通过缩小窗口来减小这种方向偏差为平滑带来的影响，因此具有边缘保持最优化的效果；与此同时，该算法设置当窗口缩小为33时，像

16、元中心像素不进行平滑处理而是保持原来灰度，已达到保持孤立点的作用，同理，很细很尖锐的边缘也会被很好地保持下来，但这一原理也保存了一些噪声点，导致部分噪声被锐化。5结语EPOS滤波在去噪时具有很好的边缘和细节保持效果，是较为理想的椒盐或斑点噪声去除方法。但在应用这一滤波方法时，一些较大的噪声会被当做孤立点反而得到放大。因此，这一方法主要适用于均匀细小的椒盐噪声和斑点噪声的去除，对于较大的噪声，将达不到理想的效果。(收稿日期：2009-10-12)!飞天诚信开通农行新K宝网银安全绿色通道日前，飞天诚信科技有限公司在总结了自身服务于五十余家金融机构的经验基础上，凭借十余年在智能身份认证领域的技术积累，倾力推出了中国农业银行新K宝技术支持中心，此平台将致力于为广大飞天新K宝的用户提供全新的服务体验和完善的技术支持。据了解，飞天诚信公司是中国农业银行网上银行推出的新K宝的供应商之一。其为农业银行提供的五彩K宝被誉为新颖时尚的、功能更强的USBKey产品。此次推出的新K宝技术支持中心本着以客户为中心的原则，在内容上包含了飞天诚信中国农业银行新K宝在程序安装、证书下载、用户应用以及工具使用方面的问题。此外为了更方便客户的