资源描述：
《数字图像处理上机报告04》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、一、图像复原1.题目要求：已知P03-05-01为原图象，设计一个高斯低通滤波器，对其进行滤波,再对其施加椒盐噪声，其噪芦参数和截止频率口行选择，采用维纳滤波方法实施图彖复原。2.方法：先读取图像，利用公式HM=e5（其中do为截止频率）对图像进行高斯低通滤波，再调用inmoise对滤波后的图像加椒盐噪声，最后利用公式F(w,v)=GM对图像进行维纳滤波。
2、H(S)2
3、+3.程序代码：clear;picture=imread('pO3-06-01.bmp1);subplot(221);imshow(picture);tit
4、le(*原图*);pfft=fft2(picture);pfft=fftshift(pfft);[MzN]=size(pfft);picturel=zeros(M,N);D=50;m=fix(M/2);n=fix(N/2);fori=l:Mforj=l:Nd=sqrt((i-m)八2+(j-n)八2);hl=exp(-(dA2)/(2*(DA2)));picturel(i,j)=pfft(i,j)*hl;endendpicture3=ifftshift(picturel);picture3=uint8(ifft2(pict
5、ure3));picture2=imnoise(picture3,1salt&pepperfz0.02);subplot(222);imshow(picture2);titled处理后的图像pfft2=fft2(picture2);pfft2=fftshift(pfft2);picture2=double(picture2);su=0;fori=l:Mforj=l:Nsu=picture2(i,j)+su;endendaver=su/(M*N);fori=l:Mforj=l:Ns=(picture2(i,j)-aver)八
6、2;endendpicture4=zeros(M,N);s=s/(M*N);fori=l:Mforj=l:Nd=sqrt((i-m)八2+(j-n)八2);hl=exp(-(dA2)/(2*(DA2)));h=(hl^2/(h"2+0.02/s))/hl;a=real(h);b=imag(h);c=real(pfft2(izj));d=imag(pfft2(i,j));pictured(izj)=complex(a*c-bF?a*d+b*c);endendpicture5=ifftshift(picture4);pictu
7、re5=uint8(ifft2(picture5));subplot(223);imshow(picture5);title(，维纳滤波后的图像J;1.运行结果：町以看出，维纳滤波对椒盐噪声的抑制能力并不强，滤波效果不怎么好。二、几何校正1•题口要求：将图象P03-06-02修正为三角菱形2方法:先读取图像,利用线性畸变求出函数hl(x/y)/h2(x/y),再用双线性内插的方法得到像素的值。3•程序代码：clear;picture=imread(fp03-10-01.tiff);subplot(121);imshow(p
8、icture);[MfN]=size(picture);x=[385;755;35];ml=[l,400zl;lz780,780;1,20,780];a=(inv(ml))*x;%x变换的矩阵y=[20;725;700];b=(inv(ml))*y;%y变换的矩阵picturel=zeros(M’N);fori=l:Mforj=l:NX=[l,i,j]*a;y=[lzirj]*b;c=x-fix(x);d=y-fix(y);picturel(i,j)=(1-c)*(1-d)*picture(fix(x),fix(y))+(
9、1-c)*d*picture(fix(x)rfix(y)+1)+c*(1-d)^picture(fix(x)+1zfix(y))+c*d*picture(fix(x)+1zfix(y)+l);endendsubplot(122);imshow(picturel);4•运行结果：运行结果不太理想，口J能是出于处标变换的函数没冇求准确。三、图象重建1•题目要求：000000000000000000000000000000001111111111000001111111111000001100000000000001100000
10、00000000111111111000001111111111000001111111110000000000000110000000000000110000011111111110000011111111110000000000000000000000000000000002•方法：对该图像进行[0
