數(shù)字圖像處理-試驗(yàn)報(bào)告書(shū)七圖像復(fù)原

電子信息學(xué)院實(shí)驗(yàn)報(bào)告書(shū)課程名： 《數(shù)字圖像處理》題目：實(shí)驗(yàn)七圖像復(fù)原實(shí)驗(yàn)實(shí)驗(yàn)類(lèi)別： 【設(shè)計(jì)】 班級(jí)：BX0901學(xué)號(hào)：111103020103姓名： 竇中鋒評(píng)語(yǔ)：學(xué)習(xí)態(tài)度:【很好】【一般】【較差】程序編寫(xiě):【完整】【部分完整】【不完整】得出結(jié)論:【正確】【部分正確】【不正確】報(bào)告書(shū)寫(xiě):【規(guī)范】【一般】【不規(guī)范】成績(jī)： 指導(dǎo)教師： 范光宇批閱時(shí)間：2012年月日《數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告1、實(shí)驗(yàn)?zāi)康恼莆請(qǐng)D像復(fù)原的原理和方法；掌握利用MATLAB程序進(jìn)行圖像復(fù)原的過(guò)程。2、實(shí)驗(yàn)原理.圖像復(fù)原原理圖像復(fù)原：試圖利用退化過(guò)程的先驗(yàn)知識(shí)，去恢復(fù)已被退化圖像的本來(lái)面目。成像系統(tǒng)受各種因素的影響，導(dǎo)致了圖像質(zhì)量的降低，或者說(shuō)是退化。由于獲得圖像的方法不同(光學(xué)、光電子或電子等)，有多種退化形式，都使成像的分辨率和對(duì)比度退化，例如：傳感器噪聲攝像機(jī)聚焦不佳物體與攝像機(jī)之間的相對(duì)移動(dòng)光學(xué)系統(tǒng)的象差成像光源和射線的散射圖像退化的主要表現(xiàn)形式：圖像模糊圖像受到干擾由于成像系統(tǒng)造成圖像退化的典型現(xiàn)象是模糊，所以圖像復(fù)原的一個(gè)基本任務(wù)就是去模糊。圖像復(fù)原的基本思路：先建立退化的數(shù)學(xué)模型，然后根據(jù)該模型對(duì)退化圖像進(jìn)行擬合。圖像復(fù)原模型可以用連續(xù)數(shù)學(xué)和離散數(shù)學(xué)處理，處理項(xiàng)的實(shí)現(xiàn)可在空間域卷積，或在頻域相乘。圖像復(fù)原的方法：空間域?yàn)V波和頻域?yàn)V波。.使用MATLAB對(duì)圖像添加模糊與噪聲，再進(jìn)行圖像復(fù)原若不存在噪聲，圖像的模糊結(jié)果就完全來(lái)源于點(diǎn)擴(kuò)散函數(shù)PSF的作用，通過(guò)精確描述失真的PSF對(duì)模糊圖像進(jìn)行卷積操作，就可實(shí)現(xiàn)去模糊。分別創(chuàng)建一個(gè)仿真運(yùn)動(dòng)模糊PSF和一個(gè)均值濾波PSF并模糊原圖像%使用函數(shù)fspecial創(chuàng)建點(diǎn)擴(kuò)散函數(shù)PSFI=imread('girl.bmp');%I=I(10+[1:256],100+[1:256],:); %裁剪圖像subplot(1,2,1);imshow(I);title('originalimage')LEN=40; %指定運(yùn)動(dòng)位移為40個(gè)像素THETA=30; %運(yùn)動(dòng)角度為30度PSF=fspecial('motion',LEN,THETA);Blurred=imfilter(I,PSF,'circular','conv');subplot(1,2,2);imshow(Blurred);title('Blurredimage')H=fspecial('disk',10);%均值濾波PSFblurredfilter=imfilter(I,H);figure,imshow(blurredfilter);title('FilterBlurredimage')《數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告 -2-3、實(shí)驗(yàn)步驟(1)打開(kāi)計(jì)算機(jī)，啟動(dòng)MATLAB程序，設(shè)置MATLAB的程序組中的當(dāng)前活動(dòng)文件夾；(2)找到待處理的圖像文件；(3)根據(jù)實(shí)驗(yàn)內(nèi)容和實(shí)驗(yàn)要求進(jìn)行實(shí)驗(yàn)；(4)記錄和整理實(shí)驗(yàn)報(bào)告。4、實(shí)驗(yàn)內(nèi)容(1)編寫(xiě)MATLAB程序，對(duì)一圖像添加不同噪聲，并觀察其效果和直方圖。如例6.1所示。(2)編寫(xiě)MATLAB程序，用均值濾波方法對(duì)含噪圖像進(jìn)行濾波，觀察其效果。如例6.2所示。(3)編寫(xiě)MATLAB程序，用標(biāo)準(zhǔn)的均值、中值濾波方法對(duì)含噪圖像進(jìn)行濾波，觀察其效果。如例6.3所示。5、實(shí)驗(yàn)程序與結(jié)果(結(jié)果可以抓圖粘貼)(1)編寫(xiě)MATLAB程序，對(duì)一圖像添加不同噪聲，并觀察其效果和直方圖。如例6.1所示。程序：I=imread('Lena.bmp')figure,imshow(I),figure,hist(double(I),10)J=imnoise(I,'gaussian',0.05);figure,imshow(J),figure,hist(double(J),10)J=imnoise(I,'speckle',0.05);figure,imshow(J),figure,hist(double(J),10)J=imnoise(I,'salt&pepper,,0.05);figure,imshow(J),figure,hist(double(J),10)I=imread('Lena.bmp')figure,imshow(I),figure,hist(double(I),10)J=imnoise(I,'gaussian',0.05);figure,imshow(J),figure,hist(double(J),10)J=imnoise(I,'speckle',0.05);figure,imshow(J),figure,hist(double(J),10)J=imnoise(I,'salt&pepper,,0.05);figure,imshow(J),figure,hist(double(J),10)截圖：數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告-l°|x|JiltEdiILitInsertI-xls”邛」小I4]j□LSH昌上A/7西R?二Fil白Bd^tVlmIncAftToolsYlindovHalpU□?Pm*A/ |WK600 1 1 1 1 1數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告 -4-數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告 -5-數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告PileEdit山fIniM-l：InnisV:ndovHelp(2)編寫(xiě)MATLAB程序，用均值濾波方法對(duì)含噪圖像進(jìn)行濾波，觀察其效果。如例6.2所示。程序：img=imread('Lena.bmp');figure;imshow(img);img_noise=double(imnoise(img,'gaussian',0.06));figure,imshow(img_noise,[]);img_mean=imfilter(img_noise,fspecial('average',3));figure;imshow(imgmean,[]);title('de-noisebymeanfilter,);img_mean=exp(imfilter(log(img_noise),fspecial('average',3)));數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告figure;imshow(img_mean,[]);title('de-noisebymeanfilter,);Q=-1.5;img_mean=imfilter(img_noise「(Q+1),fspecial('average',3))./imfilter(img_noise「Q,fspecial('average',3));figure;imshow(img_mean,[]);title('de-noisebymeanfilter,);Q=1.5;img_mean=imfilter(img_noise「(Q+1),fspecial('average',3))./imfilter(img_noise「Q,fspecial('average',3));figure;imshow(img_mean,[]);title('de-noisebymeanfilter,);截圖:數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告10-10-數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告觀察其效果。(3)編寫(xiě)MATLAB程序，用標(biāo)準(zhǔn)的均值、中值濾波方法對(duì)含噪圖像進(jìn)行濾波，如例6.3所示。觀察其效果。程序：img=imread('Lena.bmp');figure;imshow(img);img_noise=double(imnoise(img,'salt&pepper',0.06));figure,imshow(img_noise,口)；img_mean=imfilter(img_noise,fspecial('average',5));figure;imshow(img_mean,口)；title('de-noisebymeanfilter,);img_median=medfilt2(img_noise);figure;imshow(img_median,口)；title('de-noisebymedianfilter,);img_median2=medfilt2(img_median);figure;imshow(img_median2,口)；title('de-noisebymedianfilter,);截圖：-11--11-數(shù)字圖像處理》實(shí)驗(yàn)報(bào)告二二212sE，iIeEdit^3n-Lai-irt工diItl￡a