[Matlab] 3D plot with map

S = shaperead('landareas.shp') 

%mapshow(S); 

plot([S.X],[S.Y],'k') 

xLim = [115, 130];

yLim = [20, 35];

set(gca,'XLim', xLim, 'YLim', yLim); 

rotate3d 

hold on

plot3(lon,lat,alt,'.','color',[0.0 1.0 1.0],'MarkerSize',5)