《图像复原程序》doc版.doc

面是一个人写的盲去卷积算法实现图像去模糊，其中有些代码不是很能理解，求解释！尤其每次都做这个循环：for i=1:x_size %进行乘(-1)(i+j)处理for j=1:y_sizeflag=(-1)(i+j);tmp1(i,j)=tmp1(i,j)*flag;endend到底是干嘛？最好能对整段代码加些注释再贴到回复里面。谢谢！%程序内容：盲迭代反卷积方法恢复图象%编程时间：2007年12月16日%clear all;gray_S=imread(cameraman.jpg); %读入图像figure;imshow(gray_S);title(original image) %显示原始图x_size,y_size=size(gray_S);%求取PSF尺寸大小（圆形）p_size=5; %定义psf尺寸大小p_support=zeros(x_size,y_size); %psf掩模p_support(1:p_size,1:p_size)=1;beta=0.9; %参数设置g=im2double(gray_S); %转换为double型tmp1=g; flag=1; for i=1:x_size %进行乘(-1)(i+j)处理for j=1:y_sizeflag=(-1)(i+j);tmp1(i,j)=tmp1(i,j)*flag;endendG=fft2(tmp1);F0=fft2(tmp1);F1=zeros(x_size,y_size);%-循环迭代for k=1:40 %确定迭代次数 H=G./F0; %确定H，并且把小于0的部分置为0h=ifft2(H);tmp1=h;for i=1:x_sizefor j=1:y_sizeflag=(-1)(i+j);tmp1(i,j)=tmp1(i,j)*flag;endendtmp2=real(tmp1); %取实部tmp3=find(tmp2=abs(G(:);a2=find(abs(H(:)abs(G(:); F1(a1(:)=(1-beta)*F0(a1(:)+beta*(G(a1(:)./H(a1(:);F1(a2(:)=(1-beta)*F0(a2(:).(-1)+beta*H(a2(:)./G(a2(:).(-1);f1=ifft2(F1); %把图像空域小于0的部分置0tmp1=f1;for i=1:x_size %进行乘(-1)(i+j)处理for j=1:y_sizeflag=(-1)(i+j);tmp1(i,j)=tmp1(i,j)*flag;endendtmp2=real(tmp1);tmp3=find(tmp20);tmp2(tmp3(:)=0;f1=tmp2;tmp1=f1;for i=1:x_size %进行乘(-1)(i+j)处理for j=1:y_sizeflag=(-1)(i+j);tmp1(i,j)=tmp1(i,j)*flag;endendF0=fft2(tmp1); endfigure;imshow(uint8(real(f1)*255/max(f1(:);title(result image) %显示结果图imwrite(uint8(real(f1)*255/max(f1(:),best.bmp);%mesh(h(1:p_size,1:p_size)m=size(I,1);n=size(I,2);T(m,n)=0;min=5000000;%灰度值的绝对值之和第一次判断界值degree=0;%初始运动角度设为0for k=-pi/2:1*pi/180:pi/2 for i=3:m-3for j=3:n-3ii=i+2*sin(k);jj=j+2*cos(k);rii=round(ii);rjj=round(jj);x=ii-rjj;y=jj-rjj;T(i,j)=(1-x-y+x*y)*J(rii,rjj)+(y-x*y)*J(rii,rjj+1)+(x-x*y)*J(rii+1,rjj)+x*y*J(rii+1,rjj+1);T(i,j)=abs(T(i,j)-J(i,j);%双线性插值求的k角度下的微分图像求的a角度下的微分图像endenda=0;for i=1:mfor j=1:na=a+T(i,j);endendif aminmin=a;degree=-k*180/pi;endendI=imread(IMG_20130426_182723525.jpg); m=size(I,1); n=size(I,2); T(m,n)=0; min=5000000;% degree=0;%0 for k=-pi/2:1*pi/180:pi/2 for i=3:m-3 for j=3:n-3 ii=i+2*sin(k); jj=j+2*cos(k); rii=floor(ii); rjj=floor(jj); x=ii-rjj; y=jj-rjj; I(i,j)=(1-x-y+x*y)*J(rii,rjj)+(y-x*y)*J(rii,rjj+1)+(x-x*y)*J(rii+1,rjj)+x*y*J(rii+1,rjj+1); T(i,j)=abs(T(i,j)-J(i,j);