I = imread('mmm.bmp');
imshow(I)
pp=rgb2gray(I);
imshow(pp)
E=edge(pp,'canny');
imshow(E)
d1=0;
t1=0;
d2=0;
t2=0;
[m,n]=size(E);
%.................................
for y1= 1:m
for x1= 1:n
  t1=t1+1;
  R1(y1,x1)=E(y1,x1);
  if R1(y1,x1) >0
      d1=d1+1;
      ax(d1)=x1;
     
end
end
end
%.....................................
for x2= 1:n
for y2= 1:m
  t2=t2+1;
  R2(y2,x2)=E(y2,x2);
  if R2(y2,x2) >0
      d2=d2+1;
      ay(d2)=y2;
     
end
end
end
%................
   r=54;
   g1=0;
i1=0;
 for h1 = 1:d1
if h1 >1
     g1=g1+1;
        xi(g1)  = ax(h1)-ax(h1-1) ;
 
        if sqrt(xi(g1)*xi(g1))>(2*r)-1 &sqrt(xi(g1)*xi(g1))<(2*r)+1
            i1=i1+1;
        ux(i1)=xi(g1);
          
            sx(i1)=ax(h1-1)+((ax(h1)-ax(h1-1))/2);
      
        end
    end   
end
%..................................
  g2=0;
i2=0;
 for h2 = 1:d2
if h2 >1
     g2=g2+1;
        yi(g2)  = ay(h2)-ay(h2-1) ;
 
        if sqrt(yi(g2)*yi(g2))>(2*r)-3 &sqrt(yi(g2)*yi(g2))<(2*r)+3
            i2=i2+1;
        uy(i2)=yi(g2);
                        
            sy(i2)=ay(h2-1)+ ((ay(h2)-ay(h2-1))/2) ;
            
                
        end
    end   
end
///////////////////////////////////
if i1 >1
ix =find(abs(mean(sx)-sx)<=10);
Cx=mean(sx(ix));
else
    Cx=sx;
end

if i2 >1
iy =find(abs(mean(sy)-sy)<=10);
Cy=mean(sy(iy));
else
    Cy=sy;
end
%%%%%%%%%%%%%%%%%%%%