1、1 維納濾波復原 I=checkerboard(8);%產生8x8的棋盤矩陣圖像noise=0.1*randn(size(I);%產生隨機噪聲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(rea

2、l(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(22

3、3);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)');規則化濾波復原程序 I=checkerboard(8);PSF=fspecial('gaussian',7,10);%生成高斯低通濾波器V=.01;BlurredNo

4、isy=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);%顯示去模糊的復原圖像tit

5、le('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-Rich

6、ardson復原程序舉例(我們組丁旭已做好這個了)%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(Blu

7、rredNoisy,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盲去卷積復原%

8、盲去卷積復原I=checkerboard(8); 創建一個8*8棋盤圖像PSF=fspecial('gaussian',7,10);% 生成高斯低通濾波器（創建點擴散函數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');%建立函數J P=deconvblind(BlurredNoisy,INITPSF,20,10*sqrt(V),WT,FUN,0);%對blurrednoisy圖像進行盲圖像復原，迭代次數為20次，產生一個閥值為0.1偏移矩陣subplot(221);imshow(BlurredNoisy);%顯示模糊加噪的圖像title('A=Blurred and Noisy');subplot(222);imshow(PSF,);%生成真實的PSF圖像title('True PS