% get image(s) to rectify

I = imread('calib1.jpg');
imshow(I)

% specify four points in the source image

[x1, y1] = ginput(4)