資源描述：
《數(shù)字圖像處理上機報告》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在學術(shù)論文-天天文庫。
1、數(shù)字圖像處理上機報告(3)姓名班級學號指導教師傅華明、上機題目一、圖象復原3.己知p03-03-01為原圖象，pO3-O3-O2為經(jīng)過高通濾波以后的退化圖象，請采川逆濾波方法，假設(shè)高通濾波器為（1）高斯濾波器（2）巴特沃茨濾波器。每組濾波器采用3組高通截止頻率。并計算均方誤差最小，說明哪一組參數(shù)復原的閣象復原效果最好？二、幾何校正1.將閣象P03-06-01修正為菱形三、圖象重建（根據(jù)傅立葉反變換實施圖象重建）2.00000000000000000000000000000000111111111110000111I11111110000000000
2、1111100000000011111000000000111110000000001111100000000°I111000000000111110000000001111100000000001111111111100001111111111100000000000000000000000000000000二、/A*:復原1.題目要求1.己知p03-03-01為原圖象，pO3-O3-O2為經(jīng)過高通濾波以后的退化圖象，請采川逆濾波方法,假設(shè)高通濾波器為（1）高斯濾波器（2）巴特沃茨濾波器。每組濾波器采用3組高通截止頻率。并計算均方誤差最小，說明哪
3、一組參數(shù)復原的閣象復原效果最好？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
