function img(IMAGE)


%Construct a full path name
str=[IMAGE, '.tiff'];
IMG=imread(str);

%Show the image
figure
imshow(IMG)

% Read a couple of pixel positions; 
% The user must click several times in the image and then 
% press enter or middle click
[c,r,p]=impixel

%Crop 20x20 regions from the image centered at the position of the clicks
A=IMG( r(1)-10: r(1)+10 , c(1)-10: c(1)+10 );
B=IMG( r(2)-10: r(2)+10 , c(2)-10: c(2)+10 );
C=IMG( r(3)-10: r(3)+10 , c(3)-10: c(3)+10 );


%Show a mosaic of the first three regions
blank=zeros(21,21);

figure 
imshow([A, blank, B, blank, C])
