完整的圖像復原.doc

1 維納濾波復原 I=checkerboard(8);%產(chǎn)生8x8的棋盤矩陣圖像noise=0.1*randn(size(I);%產(chǎn)生隨機噪聲PSF=fspecial(motion,21,11);%生成模糊濾波器Blurred=imfilter(I,PSF,circular);%將原始圖像I按模糊濾波器進行濾波增強處理BlurredNoisy=im2uint8(Blurred+noise);%生成加噪模糊圖像NP=abs(fftn(noise).2;%計算NPOW=sum(NP(:)/prod(size(noise);NCORR=fftshift(real(ifftn(NP);IP=abs(fftn(I).2;IPOW=sum(IP(:)/prod(size(noise);ICORR=fftshift(real(ifftn(IP);ICORR1=ICORR(:,ceil(size(I,1)/2);NSR=NPOW/IPOW;subplot(221);imshow(BlurredNoisy,);title(模糊和噪聲圖像);subplot(222);imshow(deconvwnr(BlurredNoisy,PSF,NSR),);title(deconbwnr(A,PSF,NSR);subplot(223);imshow(deconvwnr(BlurredNoisy,PSF,NCORR,ICORR),);title(deconbwnr(A,PSF,NCORR,ICORR);subplot(224);imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),);title(deconbwnr(A,PSF,NPOW,ICORR_1_D);規(guī)則化濾波復原程序 I=checkerboard(8);PSF=fspecial(gaussian,7,10);%生成高斯低通濾波器V=.01;BlurredNoisy=imnoise(imfilter(I,PSF),gaussian,0,V);%對模糊圖像加高斯白噪聲，噪聲強度為0.01NOISEPOWER=V*prod(size(I);%計算噪聲強度J LAGRA=deconvreg(BlurredNoisy,PSF,NOISEPOWER);%使用平滑約束最小而成濾波算法對圖像去卷積subplot(221);imshow(BlurredNoisy);%顯示blurrednoisy圖像title(A=Blurred and Noisy);subplot(222);imshow(J);%顯示去模糊的復原圖像title(J LAGRA=deconvreg(A,PSF,NP);subplot(223);imshow(deconvreg(BlurredNoisy,PSF,LAGRA/10);%顯示使用小范圍搜索的拉普拉斯算子對圖像進行復原title(deconvreg(A,PSF,0.1*LAGRA);subplot(224);imshow(deconvreg(BlurredNoisy,PSF,LAGRA*10);%顯示使用大范圍搜索的拉普拉斯算子對圖像進行復原title(deconvreg(A,PSF,10*LAGRA);3Lucy-Richardson復原程序舉例(我們組丁旭已做好這個了)%Lucy-Richardson復原程序I=checkerboard(8);PSF=fspecial(gaussian,7,10);V=.0001;BlurredNoisy=imnoise(imfilter(I,PSF),gaussian,0,V);WT=zeros(size(I);WT(5:end-4,5:end-4)=1;J1=deconvlucy(BlurredNoisy,PSF);J2=deconvlucy(BlurredNoisy,PSF,20,sqrt(V);J3=deconvlucy(BlurredNoisy,PSF,20,sqrt(V),WT);subplot(221);imshow(BlurredNoisy);title(A=Blurred and Noisy);subplot(222);imshow(J1);title(deconvlucy(A,PSF);subplot(223);imshow(J2);title(deconvlucy(A,PSF,NI,DP);subplot(224);imshow(J3);title(deconvlucy(A,PSF,NI,DP,WT);4盲去卷積復原%盲去卷積復原I=checkerboard(8); 創(chuàng)建一個8*8棋盤圖像PSF=fspecial(gaussian,7,10);% 生成高斯低通濾波器（創(chuàng)建點擴散函數(shù)PSF）V=.0001;BlurredNoisy=imnoise(imfilter(I,PSF),gaussian,0,V); 對模糊圖像添加高斯噪聲WT=zeros(size(I);%生成大小為8x8的0矩陣WT(5:end-4,5:end-4)=1; 把矩陣里面的0用1替換INITPSF=ones(size(PSF);%生成大小與PSF相同的1矩陣FUN=inline(PSF+P1,PSF,P1);%建立函數(shù)J P=deconvblind(BlurredNoisy,INITPSF,20,10*sqrt(V),WT,FUN,0);%對blurrednoisy圖像進行盲圖像復原，迭代次數(shù)為20次，產(chǎn)生一個閥值為0.1偏移矩陣subplot(221);imshow(BlurredNoisy);%顯示模糊加噪的圖像title(A=Blurred and Noisy);subplot(222);imshow(PSF,);%生成真實的PSF圖像ti