资源描述：
《实验四图像复原》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验四图像复原一、实验目的利用反向滤波和维纳滤波进行图像复原。二、实验内容给定一个图像及其降质图像，分别为peppers.tiff和peppers_degraded.tiff 。降质模型已给出。与降值图象相关的参数是：N=256，n=5，1．利用反向滤波方法进行图像复原。2．利用维纳滤波方法进行图像复原，K＝17。对每一种方法都要求计算复原出来的图像的峰值信躁比。3．将原始图像、降质图、恢复后的图像同时显示出来，以便比较。三、源程序：%-------------------------------------------
2、-------------------------functionvarargout=pushbutton2_Callback(h,eventdata,handles,varargin)%复原clear;f1=zeros(256,256);f2=zeros(256,256);gf=zeros(256,256);h=zeros(256,256);h2=zeros(256,256);f=zeros(256,256);ff2=zeros(256,256);p1=0;p2=0;g=im2double(imread('pepper
3、s_degraded.tif'));f=im2double(imread('peppers.tif'));gf=fft2(g);fori=1:256forj=1:256h(i,j)=(1/(1+(((i-129)^2+(j-129)^2)/36^2)^5));endend%反向滤波4fori=1:256forj=1:256ff(i,j)=1/256^2*gf(i,j)/h(i,j);endendf1=ifft2((ff))/max(max(ifft2(ff)));%维纳滤波gf=fftshift(gf);fori=1:2
4、56forj=1:256ff2(i,j)=gf(i,j)*(1/h(i,j)*(abs(h(i,j))^2/((abs(h(i,j))^2+17))));endendf2=ifft2(fftshift(ff2))/max(max(ifft2(fftshift(ff2))));h2=ff2./gf;%求峰值信噪比fori=1:256forj=1:256p1=((g(i,j))-f1(i,j))^2+p1;p2=((g(i,j))-f2(i,j))^2+p2;endendpmse1=-10*log10(1/256^2*p1/
5、double((max(max(g))))^2)pmse2=-10*log10(1/256^2*p2/double((max(max(g))))^2)subplot(2,2,1)imshow(f)subplot(2,2,2)imshow(g)subplot(2,2,3)imshow(abs(f1))subplot(2,2,4)imshow(abs(f2))%--------------------------------------------------------------------functionvarargo
6、ut=pushbutton4_Callback(h,eventdata,handles,varargin)%边沿检测clear;f=double(imread('lena.bmp'));f1=zeros(512,512);f2=zeros(512,512);f3=zeros(512,512);4fori=2:511forj=2:511f1(i,j)=(sqrt((f(i,j)-f(i-1,j-1))^2+(f(i-1,j)-f(i,j-1))^2));f3(i,j)=(abs(f(i+1,j)+f(i-1,j)+f(i,
7、j+1)+f(i,j-1)-4*f(i,j)));f2(i,j)=(sqrt((f(i-1,j+1)+2*f(i,j+1)+f(i+1,j+1)-(f(i-1,j-1)+2*f(i,j-1)+f(i+1,j-1)))^2+(f(i-1,j-1)+2*f(i-1,j)+f(i-1,j+1)-(f(i+1,j-1)+2*f(i+1,j)+f(i+1,j+1))^2)));endendf1=uint8(f1);f2=uint8(f2);f3=uint8(f3);f=uint8(f);figureimshow(f)title('
8、原图')figureimshow(f1)title('robert');figureimshow((f2))title('sobol')figureimshow(f3)title('laplacian')%--------------------------------------------------------
