ults in figure show

 OUTPUTS
  H        - estimated homography to transform I into Iref
  Ip       - tranformed version of I (slow to compute)

 EXAMPLE
  Iref = double(imread('cameraman.tif'))/255;
  H0 = [eye(2)+randn(2)*.1 randn(2,1)*10; randn(1,2)*1e-3 1];
  I = imtransform2(Iref,H0^-1,'pad','replicate');
  o=50; P=ones(o)*1; I(150:149+o,150:149+o)=P;
  prmAlign={'outThr',.1,'resample',.5,'type',1:8,'show'};
  [H,Ip]=imagesAlign(I,Iref,prmAlign{:},1);
  tic, for i=1:30, H=imagesAlign(I,Ire