function gao1(varargin)
file_name=varargin{1};
[con1,con2]=size(varargin');
II=imread(file_name);
I=abs(rgb2gray(II));
%I=II;
ma=max(max(I));
mi=min(min(I));
figure(1)
imshow(I,[]);
[M,N]=size(I);
if(con1>=2)
con=varargin{2};
else
con=0.75;
end
BW1 = edge(I,'canny');
%BW1 = bwmorph(I,'remove');
jmi=1024;jma=0;
for i=1:M
for j=1:N
if 0==BW1(i,j) %&& i>135 && i<243 && j>270 && j<378
BW1(i,j) =1;
else
BW1(i,j)=0;
end
end
end
s=0;
for i=1:M
for j=1:N
if 0==BW1(i,j) && i>459 && i<567 && j>243 && j<351
if j<jmi
jmi=j;
end
if j>jma
jma=j;
end
end
end
end
jmi,jma,
BW=BW1;
b(N)=0;
for i=2:M-1
for j=2:N-1
if 1==BW(i,j) && 0==BW(i-1,j) && i>459 && i<567 && j>243 && j<351
if 0==b(j) && j>jmi && j<jma
BW(i,j)=0;
else
b(j)=1;
end
if 0==BW(i+1,j)
b(j)=1;
end
end
end
end
for i=1:M
for j=1:N
if 0==BW(i,j) && i>459 && i<567 && j>243 && j<351
s=s+1;
end
end
end
s1=0;t=0;ji=1024;ja=0;
for i=1:M
for j=1:N
if 0==BW(i,j) && i>459 && i<567 && j>243 && j<351
s1=s1+1;
if j>ja
ja=j;
end
if j<ji
ji=j;
end
end
end
if s1 >=s/2
break;
end
end
s,s1,i0=i,ji,ja
if s1==s/2
i01=i0+0.5;
else
s5=0;
for i=1:i0-1
for j=1:N
if 0==BW(i,j) && i>459 && i<567 && j>243 && j<351
s5=s5+1;
end
end
end
s5
i01=i0+(s/2-s5)/(ja-ji+1)-1
end
s2=0;t=0;ii=1024;ia=0;
for j=1:N
for i=1:M
if 0==BW(i,j) && i>459 && i<567 && j>243 && j<351
s2=s2+1;
if i>ia
ia=i;
end
if i<ii
ii=i;
end
end
end
if s2 >=s/2
break;
end
end
s,s2,j0=j,ii,ia
if s2==s/2
j0=j0+0.5
else
s3=0;
for j=1:j0-1
for i=1:M
if 0==BW(i,j) && i>459 && i<567 && j>243 && j<351
s3=s3+1;
end
end
end
s3
j01=j0+(s/2-s3)/(ia-ii+1)-1
end
figure(2)
imshow(BW,[]);
figure(3)
imshow(BW1,[]);
i01,j01