figure(3);
set(3,'Name','Calibration results','NumberTitle','off');
subplot(2,1,1);
ssmy=[-2.224834e+02 ;0.000000e+00 ;1.990932e-03; -4.004428e-06; 1.218167e-08];
ssit=[-2.169126e+02 ;0.000000e+00; 1.585854e-03 ;-3.183199e-07; 3.767708e-09];
x=0:0.02*370:370;
% plot(0:floor(374.956334),polyval([ssmy(end:-1:1)],[0:floor(374.956334)]),':',0:floor(369.620933),polyval([ssit(end:-1:1)],[0:floor(369.620933)]),'k.'); grid on; axis equal; 
mmh=plot(x,polyval([ssmy(end:-1:1)],[0:0.02*double(370):double(370)]),'*k',x,polyval([ssit(end:-1:1)],[0:0.02*double(370):double(370)]),'ok');legend({'Proposed','Urban'},'Location','NorthWest'); grid off; %axis equal;
% plot(0:floor(374.956334),polyval([ssmy(end:-1:1)],[0:floor(374.956334)]),'-.',0:floor(369.620933),polyval([ssit(end:-1:1)],[0:floor(369.620933)]),'--'); grid on; axis equal; 
xlabel('Distance ρ from the image center in pixels');
ylabel('f(ρ)');
title('Forward projection function');
%
% subplot(2,1,2);
% plot(0:floor(calib_data.ocam_model.width/2),180/pi*atan2(0:floor(calib_data.ocam_model.width/2),-polyval([ss(end:-1:1)],[0:floor(calib_data.ocam_model.width/2)]))-90); grid on;
% xlabel('Distance ''rho'' from the image center in pixels');
% ylabel('Degrees');
% title('Angle of optical ray as a function of distance from circle center (pixels)');