资源描述：
《数字图像处理上机报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、数字图像处理上机报告(3)姓名班级学号指导教师傅华明、上机题目一、图象复原3.己知p03-03-01为原图象，pO3-O3-O2为经过高通滤波以后的退化图象，请采川逆滤波方法，假设高通滤波器为（1）高斯滤波器（2）巴特沃茨滤波器。每组滤波器采用3组高通截止频率。并计算均方误差最小，说明哪一组参数复原的阁象复原效果最好？二、几何校正1.将阁象P03-06-01修正为菱形三、图象重建（根据傅立叶反变换实施图象重建）2.00000000000000000000000000000000111111111110000111I11111110000000000
2、1111100000000011111000000000111110000000001111100000000°I111000000000111110000000001111100000000001111111111100001111111111100000000000000000000000000000000二、/A*:复原1.题目要求1.己知p03-03-01为原图象，pO3-O3-O2为经过高通滤波以后的退化图象，请采川逆滤波方法,假设高通滤波器为（1）高斯滤波器（2）巴特沃茨滤波器。每组滤波器采用3组高通截止频率。并计算均方误差最小，说明哪
3、一组参数复原的阁象复原效果最好？2.原理分析先进行傅立叶变换后低通滤波，再傅立叶反变换得到退化阁象，再傅立叶变换，最后逆滤波。逆滤波:厂(“，v)=G(u,V)H(w,v)i/H（u，v）称为逆滤波器。对上式再进行傅立叶反变换可得到.加，但实际上碰到的问题都是有噪声的，因而只能求F（W，0的估计值：然后再作傅立叶逆变换，得N(u,v、H(u,v)/(AV)=/(^y)+J#一切)dndv3.程序代码11:imread(’p03-03-01.bmp’)；subplot(3,3,1)，imshow(H),title(’原始图像’);12二imreadC
4、p03-03-02.bmp’)；subplot(3,3,4),imshow(12),ti11e(’退化图像’)；doubled”:double(I2):I二imreadCp03-03-01.bmp’)；k二I;I=fftshift(fft2(I));J二I;[m,n]:size⑴;x二fix(m/2):y二fix(n/2);%图象复原高斯滤波d=1d0:1;fori二1:mforj二1:nd=sqrt((i-xr2+(j-yK2);ifcf二0;h(iJ)=1./(1-1*exp((-d^2)/(2*d0^2))):I(i，j)二h(i，j).*1
5、(iJ)；elseI(i，j):I(i，j);endendend;I二real(ifft2(ifftshift(I))):subplot(3,3,2):imshowd,[]);title(’高斯滤波d=1’);11:imread(’p03-03-01.bmp’);I=double(I1)-double(I);D1:std2(I);DI%图象复原高斯滤波d=2dO二2;I:J;fori=1:mforj二1:nd二sqrt((i-x：T2+(j-yr2);if(T:0;h(i,j)=1./(1-1*exp((-(T2)/(2*d(T2)));I(i,j
6、)=h(i,j).*1(iJ)；elseI(i,j)=I(i，j);endendend:I=real(ifft2(ifftshift(I))):subplot(3,3,5)；imshow(I,口)；ti11e(’高斯滤波d=2’);II二imreadCp03-03-01.bmp’)；I=double(I1)-double(I):std2(I);D2:std2(I);D2%图象复原高斯滤波d=3d0:3;I=J；fori=1:mforj=1:nd=sqrt((i-xr2+(j-yr2);ifd~二0;h(i，j)=1./(1-1*exp((-d'2)
7、/(2*d(T2)));I(i，j)二h(i，j).*1(i，j);elseI(i,j)=I(i,j);endendend:I=real(ifft2(ifftshift(I)));subplot(3,3,8):imshow(I,[]);title(’高斯滤波d=3’)；II二imreadCp03-03-01.bmp’)；I=double(I1)-double(I):std2(I);D3二std2(I);D3I:J;%图象复原butterworth滤波D=1J=I;D0:1;foru二1:m;forv=1:n;D=sqrt((u-x)"2+(v-y)
8、"2):ifD〜=0;H1(u，v):(1+D0/Dr(2*1.5);J(u,v)=J(u,v)*H1(u，v);else
