在matlab中实现遥感影像和shp文件的结合显示

clc;close all;clear;
road=shaperead('boston_roads.shp');  %读取shape文件
figure,
mapshow('boston.tif');  %读取tif遥感影像并显示
axis image manual off; %关闭本地坐标系统
surveyFeetPerMeter = unitsratio('sf','meter'); %将米级单位转换成surveyFeet单位
colorTab='rgbcymkw'; %颜色代码表，不同等级的公路显示不同颜色
for i=:numel(road)
    x = surveyFeetPerMeter * [road(i).X]; %将shp文件坐转换为tif文件坐标
    y = surveyFeetPerMeter * [road(i).Y];
    mapshow(x,y,'Color',colorTab(road(i).CLASS)); %根据公路类别来进行不同颜色显示
end
print('-dpng','-r600','d:/boston.png');  %以600的分辨率保存结果图片