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=
1: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的分辨率保存结果图片
转载于:https://www.cnblogs.com/denny402/p/4679057.html
相关资源：遥感图像几何纠正matlab编程
