demonstration of "crude" collineation

Petr Sladek (sladep1) & Milan Kratochvil (kratom4) at CTU Prague

Contents

first approximation of collineation

"crude" point-to-point with no interpolation

clear all;
% load passpoints (input image only)
load pta.mat;

im1 = im2double(rgb2gray(imread('test.png')));

input image

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 

output image

figure(2); clf;
imshow(imo);
axis on;
hold on;

plot(ptB(1,:),ptB(2,:),'+b','MarkerSize',10,'LineWidth',1);