%%PART#1**************************************************************
I = imread('ss4.jpg');
imshow(I)
pp=rgb2gray(I);
imshow(pp)
E=edge(pp,'sobel');
imshow(E)
d=0;
t=0;
[m,n]=size(E);
%%PART#2******************************************************************
for y= 1:n
for x= 1:m
  t=t+1;
  R(x,y)=E(t);
  if R(x,y) >0
      d=d+1;
      ax(d)=x;
      ay(d)=y;
end
end
end
%%PART#3***********************************************************************
   g=0;
for h = 1:d

    for i = 1 :  360     
     r=  34 ;  
     g=g+1;
             xi(g) = ax(h)+ r * cos(i * (pi / 180));
        yi(g) =  ay(h)+r * sin(i * (pi / 180));
      
    end        
end

%%PART#4***************************************************************************
i1=0;
for ty=1 :g
i1=i1+1;
for tx=1:g
if xi(tx) ==xi(ty)
 sx(i1)=sx(i1)+1
end
end
end

i2=0;
for ty=1 :g
i2=i2+1;
for tx=1:g
if yi(tx) ==yi(ty)
 sy(i2)=sy(i2)+1
end
end
end

Cx=max(sx);  
Cy=max(sy);
%%%%%%%%*************************************************************************************