function Image1 = Circledraw( x,y,r,color)
%note: 1012 is the number of horizontal pixels, change if the pixel numbers
%change
th = 0:pi/50:2*pi;
xunit = r*22 * cos(th) + x*22;
yunit = r*22* sin(th) +(858-y*22);
plot(xunit, yunit,'Color',color);

end

