(lat,lon,'ko-')

%% 
% That flight line is not spaced very evenly.  Perhaps you want equal spacingREPLACE_WITH_DASH_DASH
% let's get a path with a point every 500 m: 

[lati,loni] = pspath(lat,lon,500); 
plotps(lati,loni,'rx-')

%% 
% Maybe you prefer spline interpolation: 

[lati_spline,loni_spline] = pspath(lat,lon,500,'method','spline'); 
plotps(lati_spline,loni_spline,'b^-')

%% 
% We can zoom in for a closer look: 

axis([-1499557   -1495441    -468396    -465326])

%% 
% With even spacing, it's now po