start_15_point.m

This is a computer vision project implemented in matlab to remove radial distortion from the image

ic1 = [311,144,1;
317,246,1;
330,453,1;
346,656,1;
363,849,1;
       
 997,150,1;
1003,232,1;
1009,398,1;
1015,563,1;
1018,724,1;

510,142,1;
610,239,1;
703,330,1;
619,431,1;
535,535,1;
630,620,1;
796,597,1];

ic3 = [143,210,1;
144,291,1;
143,452,1;
144,611,1;
143,768,1;
759,228,1;
760,301,1;
762,449,1;
762,597,1;
762,743,1;
306,214,1;
384,294,1;
462,372,1;
384,448,1;
304,529,1;
384,606,1;
542,602,1];

% im1 = imread('new\I1.jpg');
% figure;imshow(im1);
% hold on;
% scatter(ic1(:,1), ic1(:,2),'+');
% 
% im2 = imread('new\I_T.jpg');
% figure;imshow(im2);
% hold on;
% scatter(ic3(:,1), ic3(:,2),'+');

F = getF15pt2D(ic1(1:15,:),ic3(1:15,:));
F = F/F(4,4)
%% conventional F matrix
Fp = F(2:end, 2:end)
exp = null(Fp)
eyp = null(Fp')