依旧需要 Mapping ToolBox 不会安装的请看我上一篇 虽然我们只读取shp文件,但需要保证文件夹里还有shx文件及dbf文件 各省边界线绘图provinces=shaperead('bou2_4l.shp','UseGeoCoords',true);
% 绘图
worldmap('China');
geoshow(provinces)省份填色图provinces=shaperead('bou2_4p.shp','UseGeoCoords',true);
% 绘图
worldmap('China');
faceColors=makesymbolspec('Polygon',{'INDEX',[1 numel(provinces)],'FaceColor',polcmap(numel(provinces))});
geoshow(provinces,'DisplayType','polygon','SymbolSpec',faceColors)改个颜色: provinces=shaperead('bou2_4p.shp','UseGeoCoords',true);
% 插值定义颜色
C=[222,238,209;126,190,174;144,213,220;
33,118,155;30,69,128;20,49,127]./255;
C1(:,1)=interp1(0:5,C(:,1),linspace(0,5,numel(provinces)),'linear')';
C1(:,2)=interp1(0:5,C(:,2),linspace(0,5,numel(provinces)),'linear')';
C1(:,3)=interp1(0:5,C(:,3),linspace(0,5,numel(provinces)),'linear')';
% 绘图
worldmap('China');
faceColors=makesymbolspec('Polygon',{'INDEX',[1 numel(provinces)],'FaceColor',C1});
geoshow(provinces,'DisplayType','polygon','SymbolSpec',faceColors)
|