多种窗口中值滤波

1、数字图像处理报告非线性平滑一一中值滤波一. 实验目的1. 了解统计排序滤波器在数字图象处理中的作用;2. 理解中值滤波的原理、特点、适用对象；3. 掌握中值滤波的各种方法。二. 实验原理中值滤波是一种非线性信号处理方法，与其对应的中值滤波器是一种统计排序滤波器，用象素邻域内各灰度级排序的中值来代替该象素的灰度级（用一个奇数点的移动窗口，将窗口中心点的值用窗口内的中值代替），一定条件下，可以克服线性滤波器（如邻域平滑滤波等）所带来的图象细节模糊， 而且对滤除脉冲干扰及图象扫描噪声最为有效（中值滤波的输出与输入噪声的密度分布有关，对随机噪声的抑制能力，中值滤波比平均值滤波要差一些；但对脉冲干扰，特

2、别是脉冲宽度小于窗口一半、相距较远的窄脉冲干扰，中值滤 波的效果较好）设有一个一维序列f1, f2,fn，取窗口长度（点数）为 m （m为奇数），对其进行中值 滤波，就是从输入序列中相继抽出m个数fi-v, fi-1, fi, ；+（其中fi为窗口中心点值，v=（m-1）/2）,再将这m个点按其数值大小排序，取其序号为中心点的那个数作为滤波输 出。用数学公式表示为：m1y厂 Medf-, , fj, fv i N,v二 2 （ 4i）这样处理的特点是：周期小于m/2 （窗口一半）的脉冲会被抑制，而周期大于m/2 （窗口一半） 的脉冲得到保留，在窗口内单调增加或单调减少的序列，中值滤波输出信号仍

3、保持输入信号不变（利用这个特点，可以使中值滤波既能去除图象中的噪声，又能保持图象中一些物体的边缘）。同理，二维中值滤波可由下式表示：yij = Med f j （ 42）A式中：A为窗口； fj为二维数据序列。这样处理的特点是：当在窗口内，被噪声污染的 象素不超过滤波器区域的一半时，噪声会被抑制；当在窗口内，被噪声污染的象素超过滤波器区域的一半时，噪声会得到保留。二维中值滤波的窗口形状和尺寸对滤波效果影响较大,不同的图象内容和不同的应用要求，往往采用不同的窗口形状和尺寸。常用的二维中值滤波窗口有线状、十字形、方形 等。窗口尺寸一般先用 3 X 3,再取5X 5逐渐增大，直到滤波效果满意为止。就

4、一般经验 来讲，对于有缓变的较长轮廓线物体的图象，采用方形窗口为宜。对于包含有尖顶物体的图象，用十字形窗口；而窗口大小则以不超过图象中最小有效物体的尺寸为宜。如果图象中点、线、尖角细节较多，则不宜采用中值滤波。对一些内容复杂的图象，可以使用复合 型中值滤波，如中值滤波的线性组合，即将几种窗口尺寸大小和形状不同的中值滤波器复合使用，只要各窗口都与中心对称，滤波输出可保持几个方向上的边缘跳变，而且跳变幅度可调节；其线性组合方程如下：NYjj 二 ak Mfd fj(4 3)kdAk式中：ak为不同中值滤波的系数，Ak为窗口。三. 实验内容1. 分别对含有高斯噪声、瑞利噪声、椒盐噪声的图象进行中值滤

5、波，采用多种窗口形状(线形、十字形、方形)，窗口尺寸可变(先用3X 3,再取5X 5逐渐增大);2. 比较中值滤波对含有不同噪声的图象的平滑效果，讨论中值滤波最适用于平滑哪种噪声；3. 比较不同的窗口形状和窗口尺寸对整幅图象以及局部图象平滑效果的影响；4. 比较中值滤波和低通掩模法在平滑效果以及轮廓模糊程度上的差异；5. 尝试多种窗口线形组合的中值滤波。四. 实验程序清单主程序clear alll=(imread(camerama n.tif);J1= im noise(l,salt & pepper,0.02);J2=i mno ise(l,gaussia n,0.02);J3=raylrn

6、d(im2double(l);domain3ju=1 1 1%domain需为n*n的窗口，形状可自己定义1 1 11 1 1;K3j1=ord(J1,5,domain3ju); %ord函数中需传递窗口，及窗口采中值num的位置，此处为5,K3j2=ord(J2,5,domain3ju); % 矩 形 窗 num=(n*n+1)/2, 十 字 窗 num=4*(n-1)/2)*(n-1)/2)+nK3j3=ord(J3,5,domain3ju); %一字型窗 num=2*(n-1)+(n+1)/2, 其余自己推算 K301=LP(J1);x31=uint8(real(K301);K3=x31

7、K401=LP(J2); x32=uint8(real( K401);K4=x32K5=LP(J3);figure(1) subplot(221),imshow(I)title(原图像)subplot(222),imshow(J1)title(加椒盐噪声)subplot(223),imshow(J2)title(加高斯噪声)subplot(224),imshow(J3)title(加瑞利噪声) figure(2)subplot(321),imshow(uint8( K3j1)title(加椒盐后滤波)subplot(322),imshow(uint8( K3j2)title(加高斯后滤波)su

8、bplot(323),imshow( K3j3)title(加瑞利后滤波)subplot(324),imshow( K3)title( 加椒盐 butterworth 低通滤波 )subplot(325),imshow( K4)title( 加高斯 butterworth 低通滤波 )subplot(326),imshow( K5)title( 加瑞利 butterworth 低通滤波 )domain3shi=0 1 01 1 10 1 0;K3s1=ord(J1,7,domain3shi);K3s2=ord(J2,7,domain3shi);K3s3=ord(J3,7,domain3shi)

9、;figure(3)subplot(321),imshow(uint8( K3s1) title( 加椒盐后滤波 )subplot(322),imshow(uint8( K3s2) title( 加高斯后滤波 ) subplot(323),imshow( K3s3) title( 加瑞利后滤波 ) subplot(324),imshow( K3) title( 加椒盐 butterworth 低通滤波 ) subplot(325),imshow( K4) title( 加高斯 butterworth 低通滤波 ) subplot(326),imshow( K5) title( 加瑞利 butt

10、erworth 低通滤波 )domain3yi=0 0 01 1 10 0 0;K3y1=ord(J1,8,domain3yi);K3y2=ord(J2,8,domain3yi);K3y3=ord(J3,8,domain3yi);figure(4)subplot(321),imshow(uint8( K3y1)title( 加椒盐后滤波 ) subplot(322),imshow(uint8( K3y2)title( 加高斯后滤波 )subplot(323),imshow( K3y3)title( 加瑞利后滤波 )subplot(324),imshow( K3)title( 加椒盐 butte

11、rworth 低通滤波 )subplot(325),imshow( K4)title( 加高斯 butterworth 低通滤波 )subplot(326),imshow( K5)title( 加瑞利 butterworth 低通滤波 )domain5ju=1 1 1 1 11 1 1 1 11 1 1 1 11 1 1 1 11 1 1 1 1 ;K5j1=ord(J1,13,domain5ju);K5j2=ord(J2,13,domain5ju);K5j3=ord(J3,13,domain5ju);figure(5)subplot(321),imshow(uint8( K5j1) titl

12、e( 加椒盐后滤波 )subplot(322),imshow(uint8( K5j2) title( 加高斯后滤波 )subplot(323),imshow( K5j3)title( 加瑞利后滤波 )subplot(324),imshow( K3)title( 加椒盐 butterworth 低通滤波 )subplot(325),imshow( K4)title( 加高斯 butterworth 低通滤波 )subplot(326),imshow( K5)title( 加瑞利 butterworth 低通滤波 )domain5shi=0 0 1 0 00 0 1 0 01 1 1 1 10 0

13、 1 0 00 0 1 0 0 ;K5s1=ord(J1,21,domain5shi);K5s2=ord(J2,21,domain5shi);K5s3=ord(J3,21,domain5shi);figure(6)subplot(321),imshow(uint8( K5s1)title( 加椒盐后滤波 )subplot(322),imshow(uint8( K5s2)title( 加高斯后滤波 )subplot(323),imshow( K5s3)title( 加瑞利后滤波 )subplot(324),imshow( K3)title( 加椒盐 butterworth 低通滤波 )subpl

14、ot(325),imshow( K4)title( 加高斯 butterworth 低通滤波 )subplot(326),imshow( K5)title( 加瑞利 butterworth 低通滤波 )domain5yi=0 0 0 0 00 0 0 0 01 1 1 1 10 0 0 0 00 0 0 0 0 ;K5y1=ord(J1,23,domain5yi);K5y2=ord(J2,23,domain5yi);K5y3=ord(J3,23,domain5yi);figure(7)subplot(321),imshow(uint8( K5y1)title( 加椒盐后滤波 ) subplot

15、(322),imshow(uint8( K5y2)title( 加高斯后滤波 )subplot(323),imshow( K5y3)title( 加瑞利后滤波 )subplot(324),imshow( K3)title( 加椒盐 butterworth 低通滤波 )subplot(325),imshow( K4)title( 加高斯 butterworth 低通滤波 )subplot(326),imshow( K5)title( 加瑞利 butterworth 低通滤波 )中值滤波function w=ord(J,num,domain) c1,k1=size(J);n,k2=size(dom

16、ain); A=zeros(c1+n-1,k1+n-1);for i=1:c1%补零for j=1:k1A(i+(n-1)/2,j+(n-1)/2)=J(i,j);endendm=(n+1)/2; m0=(n-1)/2;B=zeros(n,n);for k=m:(c1+m-1)%窗口中心移动位置%根据窗口产矩阵for l=m:(k1+m-1) for p1=-m0:m0 for q1=-m0:m0B(p1+m,q1+m)=A(k+p1,l+q1); end endC=B.*domain;n2=n*n;D=reshape(C,1,n2) ;for x=1:n2-1%冒泡排序 for y=1:n2

17、-x if D(y)D(y+1) temp=D(y);D(y)=D(y+1);D(y+1)=temp; end end end w(k-(n-1)/2,l-(n-1)/2)=D(num); end end低通滤波 function x2=LP(I2) f=double(I2); g=fft2(f);g=fftshift(g); N1,N2=size(g); n=2; d0=50;n1=fix(N1/2); n2=fix(N2/2); fori=1:N1 for j=2:N2d=sqrt(i-n1F2+(j-n2)T); h=1心+0.414*(d/d0F(2* n); s1(i,j)=h*g(

18、i,j); if(g(i,j)50);s2(i,j)=0; else s2(i,j)=g(i,j);end end end s1=ifftshift(s1); s2=ifftshift(s2); x2=ifft2(s1);实验结果加怖斯噪声原图及加噪图像力口葛甘斤后i、;(t漩力口椒盐 gttprworth(氐1淀滅加喘刊后漩涙*1 口 禹 昼斤 bu tt g rvwo rth |t, i. 応-戒州口羽tMJlte rwort h -fHt iR-懑一3*3十字窗口3*#十字窗口3*3矩形窗口力口瑞利后滤波加萬Kfr buttsrworth j氐通澹液加高肪后滤独3*#十字窗口力口柳盐后谑迫加趟盐but怕rvrth低逋蔗激加瑞f-