完整的图像复原

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