MATLAB上实现不同窗口的中值滤波

clc picture=input('请输入任意N*N的图片文件名：','s'); I=imread(picture); J=imnoise(I,'salt & pepper',0.03); S=imnoise(I,'speckle',0.03); G=imnoise(I,'gaussian',0.03); B=input('J是椒盐噪声，S是随机噪声，G是高斯噪声请选择需噪声：'); Y=input('1是正方形，2是十字形，3是线性，4是圆形;请选择所需的窗形状：'); if Y<2 q=input('请选择你所需正方形的尺寸K1='); %X=Y(q,B); q1=mod(q,2); if q1<1 disp('输入尺寸错误，请输入奇数'); break; else X=fangxing(q,B); subplot(121),imshow(B); title('加方差为0.03椒盐噪声后的图象'); subplot(122),imshow(uint8(X)); title('你所选择窗形滤波后的图像'); imwrite(uint8(X),'zhangyang.bmp'); end elseif Y<2 q=input('请选择你所需正方形的尺寸K1='); %X=Y(q,B); q1=mod(q,2); if q1<1 disp('输入尺寸错误，请输入奇数'); break; else X=fangxing(q,B); subplot(121),imshow(B); title('加方差为0.03椒盐噪声后的图象'); subplot(122),imshow(uint8(X)); title('你所选择窗形滤波后的图像'); imwrite(uint8(X),'zhangyang.bmp'); end elseif Y<3 q=input('请选择你所需正方形的尺寸K1='); %X=Y(q,B); q1=mod(q,2); if q1<1 disp('输入尺寸错误，请输入奇数'); break; else X=shizixing(q,B); subplot(121),imshow(B); title('加方差为0.03椒盐噪声后的图象'); subplot(122),imshow(uint8(X)); title('你所选择窗形滤波后的图像'); imwrite(uint8(X),'zhangyang.bmp'); end elseif Y<4 q=input('请选择你所需正方形的尺寸K1='); %X=Y(q,B); q1=mod(q,2); if q1<1 disp('输入尺寸错误，请输入奇数'); break; else X=xianxing(q,B); subplot(121),imshow(B); title('加方差为0.03椒盐噪声后的图象'); subplot(122),imshow(uint8(X)); title('你所选择窗形滤波后的图像'); imwrite(uint8(X),'zhangyang.bmp'); end else X2=yuanxing(B); subplot(121),imshow(B); title('加你所选择的方差为0.03噪声后的图象'); subplot(122),imshow(uint8(X2)); title('你所选择窗形滤波后的图像'); imwrite(uint8(X2),'zhangyang.bmp'); end