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)
