Realizing the combination display of remote sensing image and shp file in MATLAB

clc;close all;clear;road=shaperead ('boston_roads.shp'); %read the shape file Figure,mapshow ('boston.tif'); %read the TIF remote sensing image and display axis image manual off;%close the local coordinate system Surveyfeetpermeter= Unitsratio ('SF','Meter'); %convert meter units to Surveyfeet units Colortab='rgbcymkw'; %color code table, different grades of highways show different colors forI=1: Numel (road) x= Surveyfeetpermeter * [Road (i). X]; %convert SHP files to tif file coordinates y= Surveyfeetpermeter *[Road (i).    Y]; Mapshow (x, Y,'Color', Colortab (Road (i). CLASS)); %different color display according to the road category EndPrint ('-dpng','-r600','D:/boston.png'); % Save the resulting picture at 600 resolution

Realizing the combination display of remote sensing image and shp file in MATLAB