[glat,glon] = antbounds_data('gl'); 
tracklat = [-75.4491  -75.4657  -75.3771  -75.1145];
tracklon = [-94.4371  -95.4879  -97.8617 -101.0862]; 

%% 
% For context, let's plot these data: 

% Initialize a 350 km map centererd on PIG: 
mapzoomps('pine island glacier','mapwidth',350,'inset','northeast')

% Plot a MOA image: 
modismoaps

% Plot Mouginot's grounding line: 
antbounds('gl','color','green','linewidth',2)

% Plot the track: 
plotps(tracklat,tracklon,'ro-','linewidth',2)

%% 
% Now find the intersect