tv2_woman.m

这是图像处理中的tv去噪模型的matlab程序

x=imread('barbara.bmp');
x=double(x);
x=imcrop(x,[257 1 255 255]);
I=x;
[M N]=size(x);
x=imnoise(x(:,:),'gaussian',0,0.02);
f(2:M+1,2:N+1)=I;
f(:,1)=f(:,2);f(:,N+2)=f(:,N+1);
f(1,:)=f(2,:);f(M+2,:)=f(M+1,:);
i=[];
i(2:M+1,2:N+1)=x;
i(:,1)=i(:,2);
i(:,N+2)=i(:,N+1);
i(1,:)=i(2,:);
i(N+2,:)=i(N+1,:);
u0=i;
u1=i;
L=0.1; miu=0.02;
h=1;eps=0.000001;
 u00=u0;
SNR11=[];PSNR11=[];MSE11=[];NOR2=[];PNOR2=[];
u2=u0;
for n=1:100
    z1=0.2;
    u3=u2;
    wn=1./(((u2(2:M+1,3:N+2)-u2(2:M+1,2:N+1)).^2+((u2(3:M+2,3:N+2)-u2(1:M,3:N+2)).*0.5).^2).^0.5+0.1);
    ws=1./(((u2(2:M+1,2:N+1)-u2(2:M+1,1:N)).^2+((u2(3:M+2,1:N)-u2(1:M,1:N)).*0.5).^2).^0.5+0.1);
  we=1./(((u2(3:M+2,2:N+1)-u2(2:M+1,2:N+1)).^2+((u2(3:M+2,3:N+2)-u2(3:M+2,1:N)).*0.5).^2).^0.5+0.1);
    ww=1./(((u2(2:M+1,2:N+1)-u2(1:M,2:N+1)).^2+((u2(1:M,3:N+2)-u2(1:M,1:N)).*0.5).^2).^0.5+0.1);
   hn=wn./(wn+ws+we+ww+z1);ho=z1./(wn+ws+we+ww+z1);hs=ws./(wn+ws+we+ww+z1);he=we./(wn+ws+we+ww+z1);hw=ww./(wn+ws+we+ww+z1);
    u2(2:M+1,2:N+1)=hn.*u2(2:M+1,3:N+2)+hs.*u2(2:M+1,1:N)+he.*u2(3:M+2,2:N+1)+hw.*u2(1:M,2:N+1)+ho.*u00(2:M+1,2:N+1);
    u2(:,1)=u2(:,2);
    u2(:,N+2)=u2(:,N+1);
   u2(1,:)=u2(2,:);
   u2(M+2,:)=u2(M+1,:);
   [NOR1,PNOR1]=normal1(u3,u2);
   NOR2(n)=NOR1;
   PNOR2(n)=PNOR1;
end
%figure;imagesc(u2); axis image; axis off; colormap(gray);
%[SNR1,PSNR1,MSE1]=imagesnr(f,u2);
%SNR11(n1)=SNR1;PSNR11(n1)=PSNR1;MSE11(n1)=MSE1;