hausdorff1.m

对于做模式识别的朋友有很大的帮助~HAUSDORFF距离！！！！！

i=imread('111.bmp')
i1=imread('222.bmp')
save2_filename='bg_frame';
save2_filename=strcat(save2_filename,'.bmp');
[m,n]=size(i)
index=1;
for t1=1:m
    for t2=1:n
        if i(t1,t2)>0;
            x(index)=t1;
            y(index)=t2;
            index=index+1;
            
        end
    end
end
p=1;
w=zeros(index,2);
for t1=1:m;
    for t2=1:n;
        if i(t1,t2)>0;
             w(p,:)=[t1 t2];
             p=p+1;
         end
     end
 end
x1=max(x);
x2=min(x);
y1=max(y);
y2=min(y);
%for k=x2:x1;
  %  i(k,y2)=255;
   % i(k,y1)=255;
   %end
%for k=y2:y1;
 %  i(x2,k)=255;
  %  i(x1,k)=255;
  %end
i3=size(x1-x2,y1-y2);

for x=x2:x1
    for y=y2:y1
        if i(x,y)>0;
        i3(x-x2+1,y-y2+1)=255;
    end
    end
end
subplot(1,2,1);
figure(1);
imshow(i3)
imwrite(i3,save2_filename,'bmp'); 
%**************************************************************
[mm,nn]=size(i1)
index1=1;
for tt1=1:mm
    for tt2=1:nn
        if i1(tt1,tt2)>0;
            xx(index1)=tt1;
            yy(index1)=tt2;
            index1=index1+1;
        end
    end
end
pp=1;
ww=zeros(index1,2);
for tt1=1:mm;
    for tt2=1:nn;
        if i1(tt1,tt2)>0;
             ww(pp,:)=[tt1 tt2];
             pp=pp+1;
         end
     end
 end
xx1=max(xx)+2;
xx2=min(xx)-2;
yy1=max(yy)+2;
yy2=min(yy)-2;
for k=xx2:xx1;
    i1(k,yy2)=255;
    i1(k,yy1)=255;
end
for k=yy2:yy1;
    i1(xx2,k)=255;
    i1(xx1,k)=255;
end
subplot(1,2,2);
figure(1);
imshow(i1)
%*******************************************************************
b=zeros(1,p);
% min_a=sqrt((w(1,1)-ww(1,1)).^2+(w(1,2)-ww(1,2)).^2);
for i=1:p
   min_a=sqrt((w(i,1)-ww(1,1)).^2+(w(i,2)-ww(1,2)).^2);
   for j=2:pp
        
        if sqrt((w(i,1)-ww(j,1)).^2+(w(i,2)-ww(j,2)).^2)<min_a
           
            min_a=sqrt((w(i,1)-ww(j,1)).^2+(w(i,2)-ww(j,2)).^2);
        end
        %  b(1,i)=min_a;
    end
     b(1,i)=min_a;        
 end
g=sort(b);
max_a=max(g)
%******************
c=zeros(1,pp);
% min_a=sqrt((w(1,1)-ww(1,1)).^2+(w(1,2)-ww(1,2)).^2);
for i=1:pp
   min_b=sqrt((ww(i,1)-w(1,1)).^2+(ww(i,2)-w(1,2)).^2);
   for j=2:p
        
        if sqrt((ww(i,1)-w(j,1)).^2+(ww(i,2)-w(j,2)).^2)<min_b
           
            min_b=sqrt((ww(i,1)-w(j,1)).^2+(ww(i,2)-w(j,2)).^2);
        end
        %  b(1,i)=min_a;
    end
     c(1,i)=min_b;        
 end
g1=sort(c);
max_b=max(g1)