script1

clear
close all

load 2018_02_14_SaddlePoints.mat

vps = zeros(8,2);
ndx = [1:2];
for n=1:4
    figure
    rgb = imread(files{n});
    vp = vpcalc(rgb,imgPoints(:,:,n),wpts);
    vps(ndx,:) = vp;
    ndx = ndx+2;
end

Warning: Image is too big to fit on screen; displaying at 33% 
Warning: Image is too big to fit on screen; displaying at 33% 
Warning: Image is too big to fit on screen; displaying at 33% 
Warning: Image is too big to fit on screen; displaying at 33%

sz = size(rgb);
% leave off the far-away vp
[ pm, k, qrng, epts ] = line_fit( vps(2:end,1), vps(2:end,2) );
figure
plot(epts(:,1),epts(:,2),'k','LineWidth',1.5);
hold on
plot(vps(2:end,1),vps(2:end,2),'ko','MarkerSize',9,'LineWidth',1.5);
% imshow(rgb);
vpcalc(rgb,imgPoints(:,:,4),wpts);
hold on
plot([1 sz(2) sz(2) 1 1],[1 1 sz(1) sz(1) 1],'b','LineWidth',2);
hold off
grid
axis equal
axis ij

Warning: Image is too big to fit on screen; displaying at 33%