動態(tài)模糊圖像復(fù)原MATLAB程序.doc

1、 退化程序clc;clear all;close all;I=imread(F:mmwB1圖2.jpg);%讀圖figure;subplot(4,3,1);imshow(I);title(原圖像);LEN=30;%運(yùn)動長度30THETA=30;%運(yùn)動角度30% LEN=60;% THETA=60;n=2;for i=1:3 for j=1:3PSF=fspecial(motion,LEN*i,THETA*j);%退化并研究運(yùn)動角度和長度對圖片模糊程度的影響 PSF=fspecial(motion,LEN,THETA);Blurred=imfilter(I,PSF,circular,conv);subplot(4,3,n);imshow(uint8(Blurred);title(模糊化);hold onn=n+1; endend%imwrite(Blurred,模糊60長60.png);%保存圖2、 運(yùn)動角度的求解%求解模糊運(yùn)動角度matlab代碼close all;clc;clear all;im=imread(F:mmwB1模糊60長60.png);img_gray=rgb2gray(im);%灰度化img_fft=fftshift(fft2(img_gray);N=abs(img_fft);P=(N-min(min(N)/(max(max(N)-min(min(N)*225;figure;imshow(P);title(頻譜圖（運(yùn)動角度與光斑方向垂直）);len=35;theta=0;PSF=fspecial(motion,len,theta);B=imfilter(img_gray,PSF,circular,conv);subplot(121);imshow(B);%模糊圖像B1=fft2(double(B); B2=mat2gray(log(abs(B1); subplot(122);imshow(B2);%模糊圖像的頻譜圖C=sum(B2,1);%對頻譜圖求列和m,n=size(C); x=0:1:n-1; y=C;figure,plot(x,y);title(頻譜列和曲線圖1)%繪制頻譜列和曲線圖3、 運(yùn)動長度算法%求解模糊運(yùn)動長度matlab代碼：im=imread(F:mmwB1模糊60長60.png);img_gray=rgb2gray(im);%灰度化h=fspecial(sobel);%sobel邊緣檢測img_double=double(img_gray);J=conv2(img_double,h,same);IP=abs(fft2(J);S=fftshift(real(ifft2(IP);figure;plot(S);title(模糊運(yùn)動長度);4、 噪聲分析%噪聲分析clc;clear allim=imread(F:mmwB1圖1.png);m,n,h=size(im);f11=ones(192,162,3);f22=ones(130,130,3);f33=ones(100,100,3);f44=ones(70,70,3);for i=1:190for j=1:162for k=1:3f11(i,j,k)=im(i,j,k);endendendfor i=1:130for j=501:630fork=1:3;f22(i,j-500,k)=im(i,j,k);endendfor i=721:870for j=11:170for k=1:3f33(i-720,j-10,k)=im(i,j,k);endendendfor i=761:830for j=561:630for k=1:3f33(i-760,j-560,k)=im(i,j,k);endendendfigure;subplot(221),hist(f11,100);subplot(222),hist(f22,100);subplot(223),hist(f33,100);subplot(224),hist(f44,100);title(噪聲分析2);5、 去噪還原clc;clear all;close all;I=imread(F:mmwB1圖1.png);%讀圖Len=60;Theta=60;PSF=fspecial(motion,Len,Theta); %模糊化BlurredA=imfilter(I,PSF,circular,conv);wnr1=deconvwnr(BlurredA,PSF);%維納濾波BlurredD=imfilter(I,PSF,circ,conv);INITPSF=ones(size(PSF);K DePSF=deconvblind(BlurredD,INITPSF,30);%盲去卷積法BlurredB=imfilter(I,PSF,conv);v=0.02;Blurred_I_Noisy=imnoise(BlurredB,gaussian,0,v);NP=v*prod(size(I);J=deconvreg(Blurred_I_Noisy,PSF,NP);%最小二乘法BlurredC=imfilter(I,PSF,symmetric,conv);v=0.002;BlurredNoisy=imnoise(BlurredC,gaussian,0,v);Luc=deconvlucy(BlurredNoisy,PSF,5);%L_Rl濾波subplot(221);imshow(I);title(原圖);subplot(222);imshow(BlurredA);title(模糊化);%subplot(233);imshow(wnr1);title(維納濾波);subplot(223);imshow(J);title(最小二乘法);imwrite(J,min_recover1.png);subplot(224);imshow(Luc);title(L_R法);imwrite(Luc,LR_recover1.png); 6、 截取部分圖片進(jìn)行對比程序clear all;clc;a=imread(F:mmwB1模糊60長60.png); %未處理質(zhì)量較差圖像b=a(64:120,67:126);a=imread(F:mmwmin_recover1.png);%算法處理后質(zhì)量較好圖象c=a(64:120,67:126);%從eyechart3中截取測試參考圖象，截取部分需要進(jìn)行縮放-%使之與eyechart1,eyechart2截取部分大小匹配-a=imread(F:mmwB1圖2.jpg);%高清晰參考圖象d=a(64:120,67:126);e=imresize(d,length(b(:,1),length(b(1,:),bicubic);%調(diào)整imwrite(b,area_模糊60長60.png);imwrite(c,area_最小二乘法復(fù)原圖.png);imwrite(e,area_圖2.png);subplot(1,3,1);imshow(e);title(模糊60長60截取參考);hold on;subplot(1,3,2);imshow(b);title(eyechart1截取部分);hold on;subplot(1,3,3);imshow(c);title(eyechart2截取部分);7、 模糊系數(shù)、質(zhì)量指數(shù)、PSNR的計(jì)算（評價(jià)）clc;clear;PSNRenable=1; %PSNR計(jì)算使能，為0不計(jì)算，為1，計(jì)算KBlurenable=1; %模糊系數(shù)KBlur計(jì)算使能，為0不計(jì)算，為1，計(jì)算Qenable=1; %質(zhì)量指數(shù)Q計(jì)算使能，為0不計(jì)算，為1，計(jì)算for m=1:2imsrcnamehead=area_模糊60長60; %源圖象文件名頭imsrcnameext=png; %源圖象文件名擴(kuò)展if m=1 %以area_eyechart1.bmp為測試圖象imdstname=strcat(area_圖2,.,imsrcnameext);%污染圖象文件名，可修改elseif m=2%以area_eyechart2.bmp為測試圖象imdstname=strcat(area_最小二乘法復(fù)原圖,.,imsrcnameext);%污染圖象文件名，可修改end%-iminfo=imfinfo(strcat(imsrcnamehead,.,imsrcnameext);%源圖象信息讀取 imsrc=imread(strcat(imsrcnamehead,.,imsrcnameext); %源圖象讀取imdst=imread(imdstname,imsrcnameext); %污染圖象讀取doubleimsrc=double(imsrc); %轉(zhuǎn)換為浮點(diǎn)類型doubleimdst=double(imdst); %轉(zhuǎn)換為浮點(diǎn)類型%-源圖象和污染圖象讀取W=iminfo.Width; %圖象寬H=iminfo.Height; %圖象高%-PSNR計(jì)算-if PSNRenable=1PSNR=0.0; %PSNR賦初值for j=1:Hfor i=1:WPSNR=PSNR+double(doubleimsrc(j,i)-doubleimdst(j,i)*(doubleimsrc(j,i)-doubleimdst(j,i);endendPSNR=PSNR/W/H;PSNR=10*log10(255*255/PSNR);%-PSNR計(jì)算完畢-end%-模糊系數(shù)KBlur計(jì)算-if KBlurenable=1Sin=0.0; %Sin賦初值Sout=0.0;for j=2:H-1for i=2:W-1t=doubleimsrc(j-1,i+1)+doubleimsrc(j+1,i-1)-doubleimsrc(j-1,i-1)-doubleimsrc(j+1,i+1);if t0t=-t;endSin=Sin+t; %源圖象鄰域邊緣能量計(jì)算t=doubleimdst(j-1,i+1)+doubleimdst(j+1,i-1)-doubleimdst(j-1,i-1)-doubleimdst(j+1,i+1);if t0t=-t;endendendSout=Sout+t; %污染圖象鄰域邊緣能量計(jì)算 KBlur=Sout/Sin;end%-KBlur計(jì)算完畢-%-質(zhì)量指數(shù)Q計(jì)算-if Qenable=1Q=0.0; %Q賦初值Qnum=0; %圖象以7X7塊大小計(jì)算每塊的Q，逐象素的移動塊窗口，這里Qnum為塊數(shù)量的計(jì)數(shù)for j=4:H-3for i=4:W-3midsrc=0.0;middst=0.0;varsrc=0.0;vardst=0.0; %源圖象和污染圖象塊內(nèi)的平均值和方差賦初值varsrcdst=0.0;%源圖象和污染圖象塊內(nèi)的協(xié)方差賦初值for n=-3:3for m=-3:3 midsrc=midsrc+doubleimsrc(j+n,i+m);middst=middst+doubleimdst(j+n,i+m);endendmidsrc=midsrc/49;middst=middst/49;%源圖象和污染圖象塊內(nèi)的平均值計(jì)算-for n=-3:3for m=-3:3varsrc=varsrc+(doubleimsrc(j+n,i+m)-midsrc)*(doubleimsrc(j+n,i+m)-midsrc);vardst=vardst+(doubleimdst(j+n,i+m)-middst)*(doubleimdst(j+n,i+m)-middst);varsrcdst=varsrcdst+(doubleimsrc(j+n,i+m)-midsrc)*(doubleimdst(j+n,i+m)-middst);endendvarsrc=va