[X,map]=imread('Lighthouse','jpg');
% image(X);
%colormap(map), axis off
image(X); colormap(map);
for i=1:30
    [x,y]=ginput(1)
end
