function KellerEntranceHomography
im1 = imread('keller_left.png');
im2 = imread('keller_right.png');

im_warped = zeros(2000,4000,3);

u1 = [2806 1004;
      2456 753;
      1677 1234;
      2325 1474];
  
u2 = [1483 1541;
      1948 997;
      860 843;
      587 1316];

u1 = [u1 ones(4,1)];
u2 = [u2 ones(4,1)];

H1 = ComputeHomography(u2, u1);
H2 = ComputeHomography(u1, u2);

im_warped1 = ImageWarping(im1, H1);
im_warped2 = ImageWarping(im2, H2);

im_1 = 0.5*im_warped1 + 0.5*im2;
im_2 = 0.5*im_warped2 + 0.5*im1;

figure(1)
clf;
imshow(im_1);

figure(2)
clf;
imshow(im_2);

imwrite(im_warped1, 'keller_homography_left2right1.png');
imwrite(im_warped2, 'keller_homography_right2left1.png');

imwrite(im_1, 'keller_homography_left2right.png');
imwrite(im_2, 'keller_homography_right2left.png');

% figure(2)
% clf;
% imshow(im1);
% hold on
% plot(u1(:,1), u1(:,2), 'rx');
% 
% figure(3)
% clf;
% imshow(im2);
% hold on
% plot(u2(:,1), u2(:,2), 'rx');
return;
 


