www.gusucode.com > map 案例源码 matlab代码程序 > map/OverlayVectorOnOrthophotoExample.m

    %% Overlay Vector on Orthophoto
% Overlay Boston roads on an orthophoto. Note that |mapshow| draws a new
% layer in the axes rather than replacing its contents.
%%
% Display image.

% Copyright 2015 The MathWorks, Inc.

figure
mapshow boston.tif
axis image off manual
%%
% Convert Boston roads to units of survey feet and overlay on orthophoto.
S = shaperead('boston_roads.shp');
surveyFeetPerMeter = unitsratio('sf','meter');
x = surveyFeetPerMeter * [S.X];
y = surveyFeetPerMeter * [S.Y];
mapshow(x,y)