L1 = imread('portrait.jpg');
figure;imshow(L1);
[x,y,k] = ginput(1)
[x1,y1,k] = ginput(1)
while k~=3
 L2 = circshift (L1, [round(y1-y), round(x1-x)]);
 figure;imshow(L2);
 [x,y,k] = ginput(1);
 [x1,y1,k] = ginput(1);
end