中 d/2个最高灰度值和d/2个最低灰度值，用gr(s,t)表示剩余的mn-d个像素。0

原始图片

f=imread('rice.png'); figure(1);imshow(f); title('原始图像');

g=imnoise(f,'salt& pepper',0.2); figure(2);imshow(g); title('椒盐噪声污染的的'); g1=double(g)/225;

j1=medfilt2(g1,'symmetric'); figure(3);imshow(j1); title('中值滤波图像');

j2=ordfilt2(g1,median(1:3*3),ones(3,3),'symmetric');

figure(4);imshow(j2); title('中点滤波图像'); j3=ordfilt2(g1,1,ones(3,3)); figure(5);imshow(j3); title('最小值滤波图像'); j4=ordfilt2(g1,9,ones(3,3)); figure(6);imshow(j4); title('最大值滤波图像');

椒盐噪声污染的图像中点滤波图像中值滤波图像最小值滤波图像最大值滤波图像

C=imread('tissue.png'); subplot(1,2,1); imshow(C);LEN=(30);

THETA=45;PSF=fspecial('motion',LEN,THETA); MF=imfilter(C,PSF,'circular','conv'); subplot(1,2,2), imshow(MF); imwrite(MF,tissue-MF.jpg');

F=checkerboard(8); figure(1);

imshow(F,[]);

PSF=fspecial('motion',7,45); MF=imfilter(F,PSF,'circular');

noise=imnoise(zeros(size(F)),'gaussian',0,0.001); MFN=MF+noise; figure(2); imshow(MFN,[]);

NSR=sum(noise(:).^2)/sum(MFN(:).^2); imshow(deconvwnr(MFN,PSF),[]); figure(4);

imshow(deconvwnr(MFN,PSF,NSR),[]);

五.实验分析

图像在获取过程中，由于成像系统的非线性、飞行器的姿态变化