Petr Sladek (sladep1) & Milan Kratochvil (kratom4) at CTU Prague
"crude" point-to-point with no interpolation
clear all; % load passpoints (input image only) load pta.mat; im1 = im2double(rgb2gray(imread('test.png')));
figure(1); clf; imshow(im1); axis on; hold on; % X Y component of ptAi correspondence plot(ptA(1,:),ptA(2,:),'+b','MarkerSize',10,'LineWidth',1); % forward collineation seed: A = [1.8 0.2 0.2;0.5 2.5 0;0.001 0.002 1]; % "crude" collineation [imo T1] = imcollin(A,im1); % get complete passpoints ptA->ptB eval;
1: min= 2.193420e+000 max= 1.637024e+002 2: min= 2.991027e+000 max= 2.307692e+002
figure(2); clf; imshow(imo); axis on; hold on; plot(ptB(1,:),ptB(2,:),'+b','MarkerSize',10,'LineWidth',1);