function M1 = perform_image_similitude(M,u,u1,v,v1)

% perform_image_similitude
%
%   M1 = perform_image_similitude(M,u,u1,v,v1)
%
% Compute the affine similitude that map u to u1
% and v to v1, and then resample the image M.
% p and p1 are assumed to be in [0,1]