数字图像处理上机报告04.docVIP

一、图像复原 1.题目要求：已知p03-05-01为原图象，设计一个高斯低通滤波器,对其进行滤波,再对其施加椒盐噪声,其噪声参数和截止频率自行选择,采用维纳滤波方法实施图象复原。 2.方法：先读取图像，利用公式（其中D0为截止频率）对图像进行高斯低通滤波，再调用imnoise 对滤波后的图像加椒盐噪声，最后利用公式对图像进行维纳滤波。 程序代码： clear; picture=imread(p03-06-01.bmp); subplot(221); imshow(picture); title(原图); pfft=fft2(picture); pfft=fftshift(pfft); [M,N]=size(pfft); picture1=zeros(M,N); D=50; m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); h1=exp(-(d^2)/(2*(D^2))); picture1(i,j)=pfft(i,j)*h1; end end picture3=ifftshift(picture1); picture3=uint8(ifft2(picture3)); picture2=imnoise(picture3,salt pepper,0.02); subplot(222); imshow(picture2); title(处理后的图像); pfft2=fft2(picture2); pfft2=fftshift(pfft2); picture2=double(picture2); su=0; for i=1:M for j=1:N su=picture2(i,j)+su; end end aver=su/(M*N); for i=1:M for j=1:N s=(picture2(i,j)-aver)^2; end end picture4=zeros(M,N); s=s/(M*N); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); h1=exp(-(d^2)/(2*(D^2))); h=(h1^2/(h1^2+0.02/s))/h1; a=real(h); b=imag(h); c=real(pfft2(i,j)); d=imag(pfft2(i,j)); picture4(i,j)=complex(a*c-b*d,a*d+b*c); end end picture5=ifftshift(picture4); picture5=uint8(ifft2(picture5)); subplot(223); imshow(picture5); title(维纳滤波后的图像); 运行结果： 可以看出，维纳滤波对椒盐噪声的抑制能力并不强，滤波效果不怎么好。 几何校正 1.题目要求：将图象p03-06-02修正为三角菱形 2.方法：先读取图像，利用线性畸变求出函数h1(x,y),h2(x,y)，再用双线性内插的方法得到像素的值。 3.程序代码： clear; picture=imread(p03-10-01.tif); subplot(121); imshow(picture); [M,N]=size(picture); x=[385;755;35]; m1=[1,400,1;1,780,780;1,20,780]; a=(inv(m1))*x;%x变换的矩阵 y=[20;725;700]; b=(inv(m1))*y;%y变换的矩阵 picture1=zeros(M,N); for i=1:M for j=1:N x=[1,i,j]*a; y=[1,i,j]*b; c=x-fix(x); d=y-fix(y); picture1(i,j)=(1-c)*(1-d)*picture(fix(x),fix(y))+(1-c)*d*picture(fix(x),fix(y)+1)+c*(1-d)*picture(fix(x)+1,fix(y))+c*d*picture(fix(x)+1,fix(y)+1); end end subplot(122);