matlab-code.rar_degraded image _图像复原_图像复原matlab_图像退化

f=imread('D:/picture/MRIBrain_10.bmp'); subplot(2,4,1); imshow(f); %f0=rgb2gray(f);%把三维图像f改成二维图像 f0=fspecial('motion',7,45);%1产生的移动退化滤波器,作为imfilter的滤波模板 imwrite(f0,'D:/mypicture/MRIBrain_10_1.bmp'); subplot(2,4,2); imshow(f0); f1=imfilter(f,f0,'corr','circular','same')%2为图像MRIBrain_10.bmp添加一个无能移动退化 imwrite(f1,'D:/mypicture/MRIBrain_10_2.bmp'); subplot(2,4,3); imshow(f1); f2=imnoise(f1,'gaussian',0,0.001)%3将高斯噪声叠加到已经产生移动退化的图像MRIBrain_10.bmp上 imwrite(f2,'D:/mypicture/MRIBrain_10_3.bmp'); subplot(2,4,4); imshow(f2); f3=DECONVWNR(f1,f0); imwrite(f3,'D:/mypicture/MRIBrain_10_4.bmp'); subplot(2,4,5); imshow(f3); f4=DECONVWNR(f1,f0,0.5); imwrite(f4,'D:/mypicture/MRIBrain_10_5.bmp'); subplot(2,4,6); imshow(f4); f5=DECONVWNR(f1,f0,1); imwrite(f5,'D:/mypicture/MRIBrain_10_6.bmp'); subplot(2,4,7); imshow(f5); f6=DECONVWNR(f1,f0,2); imwrite(f6,'D:/mypicture/MRIBrain_10_7.bmp'); subplot(2,4,8); imshow(f6); %省略DECONVWNR中参数NSR时，其复原效果较好；当NSR的取值越大时，其复原效果越差。