基于MATLAB的改进型中值滤波源代码设计.docVIP

基于MATLAB的改进型中值滤波源代码设计 ?%by?? 周波? 2011,10,27.4meiren=imread(meiren.bmp);jiazao=imnoise(meiren,salt pepper,0.7) ;%椒盐噪声[row,col]=size(meiren)? ;%返回矩阵的行和列。 ??????? ?%下面是扩充矩阵的部分?kc=kuochong(jiazao);?image=kc;?kc=kuochong(image);?image=kc;? kc=kuochong(image);爄mage=kc;?牋?output=image;?牋 for i=4:1:row+3牋牋牋 for j=4:1:col+3牋牋牋牋牋 ps=pans(image,i,j);牋牋牋牋牋 if(ps==1)牋牋牋牋牋牋牋 output(i,j)=image(i,j);牋牋牋牋牋 else牋牋牋牋牋牋牋 out=nchu3(image,i,j);牋牋牋牋牋牋牋 out=double(out);牋牋牋牋牋牋牋 牋牋牋牋牋 end牋牋牋牋牋牋牋 output(i,j)=ceil(out);牋牋?end牋牋牋 end牋牋? 牋牋牋 for i=1:1:row牋?牋牋牋?for j=1:1:col牋牋牋牋牋?s(i,j)=output(i+3,j+3);牋牋牋?end牋?end?shuchu=uint8(s);subplot(231),imshow(meiren);title(原图像);subplot(232),imshow(jiazao);title(加椒盐噪声后图像 p=0.7);subplot(233),imshow(shuchu);title(经过中值滤波后的图像);xitong3=medfilt2(jiazao,[3 3]);?? %调用系统函数进行中值滤波 subplot(234),imshow(xitong3),title(用Matlab系统函数进行3X3中值滤波之后的结果); xitong5=medfilt2(jiazao,[5 5]);?? %调用系统函数进行中值滤波 subplot(235),imshow(xitong5),title(用Matlab系统函数进行5X5中值滤波之后的结果);??????? function kc=kuochong(im) [row,col]=size(im); image=zeros(row+2,col+2); for i=1:1:row for j=1:1:col image(i+1,j+1)=im(i,j); end end for i=2:1:row+1 image(i,1)=im(i-1,1); end for j=2:1:col+1 image(1,j)=im(1,j-1); end for i=2:1:row+1 image(i,col+2)=im(i-1,col); end for j=2:1:col+1 image(row+2,j)=im(row,j-1); end image(1,1)=im(1,1); image(row+2,col+2)=im(row,col); image(row+2,1)=im(row,1); image(1,col+2)=im(1,col); kc=image; %by 周波¨ 2011,10,28.2 function out=nchu3(image,i1,j1) ima=[image(i1-1,j1),image(i1,j1-1),image(i1,j1),image(i1,j1+1),image(i1+1,j1)]; mi=min(ima); ma=max(ima); sz1=zeros(1,5); v=1; for k=1:1:5 if((ima(k)mi)(ima(k)ma)) sz1(v)=ima(k); v=v+1; end end if