Matlab绘制中国地图超全教程详解

 更新时间:2022年03月02日 14:53:18   作者:slandarer  
这篇文章主要介绍了如何利用Matlab绘制中国地图,文中的示例代码讲解详细,对我们学习Matlab有一定的帮助,感兴趣的小伙伴可以学习一下

依旧需要用到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)

标注省份名称

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)

% 读取首都和省级行政中心位置信息并标记
pnames=shaperead('res1_4m.shp','UseGeoCoords',true);
for i=1:numel(pnames)
    textm(pnames(i).Lat,pnames(i).Lon,pnames(i).NAME, ...
   'HorizontalAlignment','center','FontSize',5,'color','k')
end

单独绘制某个省份:

worldmap([42,55],[120,137]);

% 从中国省份数据文件中获取黑龙江省
P=shaperead('bou2_4p.shp','UseGeoCoords',true,...
  'Selector',{@(name) strcmp(name,'黑龙江省'),'NAME'});
patchm(P.Lat,P.Lon, [0.5 0.7 0.5])

中国公路交通图

borders=shaperead('bou1_4p.shp','UseGeoCoords',true);
roads=shaperead('roa_4m.shp','UseGeoCoords',true);
worldmap('China');

geoshow(borders,'FaceColor',[239,238,234]./255)
geoshow(roads,'Color',[250,227,158]./255,'LineWidth',1.5)

中国铁路交通图

borders=shaperead('bou1_4p.shp','UseGeoCoords',true);
rails=shaperead('rai_4m.shp','UseGeoCoords',true);
worldmap('China');

geoshow(borders,'FaceColor',[239,238,234]./255)
geoshow(rails,'Color',[.2,.2,.2],'LineWidth',1,'LineStyle','-.')

中国河流图

borders=shaperead('bou1_4p.shp','UseGeoCoords',true);
riversp=shaperead('hyd1_4p.shp','UseGeoCoords',true);
riversl=shaperead('hyd1_4l.shp','UseGeoCoords',true);
worldmap('China');

geoshow(borders,'FaceColor',[239,238,234]./255)
geoshow(riversp,'FaceColor',[127,141,181]./255,'EdgeColor',[127,141,181]./255)
geoshow(riversl,'Color',[127,141,181]./255,'LineWidth',1)

组合美化图

美化图一

ax=worldmap('China');
setm(ax,'FFaceColor',[153,179,204]./255)

% 绘制附近大陆
antarctica = shaperead('landareas.shp', 'UseGeoCoords', true,...
  'Selector',{@(name) strcmp(name,'Africa and Eurasia'), 'Name'});
patchm(antarctica.Lat, antarctica.Lon, [0.5 0.7 0.5])

% 绘制中国国界
bordersl=shaperead('bou1_4l.shp','UseGeoCoords',true);
geoshow(bordersl,'Color',[.3,.3,.3],'LineWidth',1)
bordersp=shaperead('bou1_4p.shp','UseGeoCoords',true);
patchm([bordersp.Lat],[bordersp.Lon],[239,238,234]./255);

% 绘制河流
riversp=shaperead('hyd1_4p.shp','UseGeoCoords',true);
riversl=shaperead('hyd1_4l.shp','UseGeoCoords',true);
geoshow(riversp,'FaceColor',[127,141,181]./255,'EdgeColor',[127,141,181]./255)
geoshow(riversl,'Color',[127,141,181]./255,'LineWidth',1)

% 绘制公路
roads=shaperead('roa_4m.shp','UseGeoCoords',true);
geoshow(roads,'Color',[250,227,158]./255,'LineWidth',1.5)

tightmap

美化图二

ax=worldmap('china');
setm(ax,'mapprojection','mercator');

% 绘制国界
bordersl=shaperead('bou1_4l.shp','UseGeoCoords',true);
geoshow(bordersl,'Color',[.3,.3,.3],'LineWidth',1)
% 绘制省份
colorList=[197,228,212;255,254,216;251,197,221;199,205,231]./255;
colorList=repmat(colorList,[250,1]);
provinces=shaperead('bou2_4p.shp','UseGeoCoords',true);
colorList=colorList(1:numel(provinces),:);
faceColors=makesymbolspec('Polygon',{'INDEX',[1 numel(provinces)],'FaceColor',colorList});
geoshow(provinces,'DisplayType','polygon','SymbolSpec',faceColors)

% 首都
pnames=shaperead('res1_4m.shp','UseGeoCoords',true);
geoshow(pnames(1).Lat,pnames(1).Lon,'DisplayType','point','Marker','p',...
    'MarkerEdgeColor',[.8,0,0],'MarkerFaceColor',[.8,0,0])

% 指北针
northarrow('latitude',50,'longitude',80,'scaleratio',.08,'FaceColor',[.4,.4,.4])

% 比例尺
scaleruler on
setm(handlem('scaleruler1'), ...
    'XLoc',-3.2e6,'YLoc',1.8e6, ...
    'MajorTick',0:500:2000,...
    'MinorTick',0:40:200,'FontSize',7)
scaleruler('units','nm')
setm(handlem('scaleruler2'), ...
    'XLoc',-3.15e6, ...
    'YLoc',1.7e6, ...
    'TickDir','down', ...
    'MajorTick',0:250:1000,...
    'MinorTick',0:40:200,...
    'MajorTickLength',km2nm(25),...
    'MinorTickLength',km2nm(12.5),'FontSize',7)

% 小地图
h2=axes('pos',[0.67 0.15 0.13 0.2]);
worldmap([1 22],[106 122])
setm(h2,'FFaceColor',[197,213,236]./255,'FEdgeColor',[.3,.3,.3],'FLineWidth',1.5)
insert1=shaperead('bou2_4l.shp','UseGeoCoords',true);
geoshow(insert1,'Color',[.2,.2,.2],'LineWidth',1)
insert2=shaperead('bou2_4p.shp','UseGeoCoords',true);
geoshow(insert2,'FaceColor','w')
mlabel off 
plabel off
gridm off

代码和数据:

链接:https://pan.baidu.com/s/1ZvMTLVElEWgQ3TFOmDXyCg  提取码:1234

到此这篇关于Matlab绘制中国地图超全教程详解的文章就介绍到这了,更多相关Matlab绘制地图内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • 深入探索C++ string的底层实现

    深入探索C++ string的底层实现

    C语言中的字符串是以字符数组的形式存储的,每个字符占用一个字节的内存空间,本文我们将和大家一起深入探讨一下string的底层实现,感兴趣的小伙伴快来和小编一起吧
    2023-08-08
  • C++类中变量也可以是引用的代码实例

    C++类中变量也可以是引用的代码实例

    今天小编就为大家分享一篇关于C++类中变量也可以是引用的代码实例,小编觉得内容挺不错的,现在分享给大家,具有很好的参考价值,需要的朋友一起跟随小编来看看吧
    2019-04-04
  • 基于make命令与makefile文件详解

    基于make命令与makefile文件详解

    下面小编就为大家分享一篇基于make命令与makefile文件详解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2018-01-01
  • C++实现LeetCode(81.在旋转有序数组中搜索之二)

    C++实现LeetCode(81.在旋转有序数组中搜索之二)

    这篇文章主要介绍了C++实现LeetCode(81.在旋转有序数组中搜索之二),本篇文章通过简要的案例,讲解了该项技术的了解与使用,以下就是详细内容,需要的朋友可以参考下
    2021-07-07
  • 一文总结C++中的异常

    一文总结C++中的异常

    异常是一种处理错误的方式,当一个函数发现自己无法处理的错误时就可以抛出异常,让函数的直接或间接调用者处理这个错误,本文给大家总结了C++中的异常,需要的朋友可以参考下
    2023-10-10
  • Qt使用windeployqt工具实现程序打包发布方法

    Qt使用windeployqt工具实现程序打包发布方法

    本文主要介绍了Qt使用windeployqt工具实现程序打包发布方法,文中通过示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-11-11
  • C数据结构循环链表实现约瑟夫环

    C数据结构循环链表实现约瑟夫环

    这篇文章主要介绍了C数据结构循环链表实现约瑟夫环的相关资料,需要的朋友可以参考下
    2017-05-05
  • C语言入门篇--充分理解操作符

    C语言入门篇--充分理解操作符

    本篇文章是基础篇,适合c语言刚入门的朋友,本文主要介绍了c语言的操作符基础理论,希望可以帮助大家快速入门c语言的世界,更好的理解c语言
    2021-08-08
  • FFmpeg中AVIOContext的使用方法详解

    FFmpeg中AVIOContext的使用方法详解

    AVIOContext是FFMPEG管理输入输出数据的结构体,这篇文章主要为大家详细介绍了这个结构体的具体使用,文中的示例代码讲解详细,需要的可以参考一下
    2023-08-08
  • Matlab实现绘制雷达图(蜘蛛图)

    Matlab实现绘制雷达图(蜘蛛图)

    这篇文章主要为大家详细介绍了如何利用Matlab实现雷达图(蜘蛛图)的绘制,文中的示例代码讲解详细,对我们学习Matlab有一定帮助,需要的可以参考一下
    2022-09-09

最新评论