在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的分辨率保存结果图片