物探论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

查看: 2327|回复: 0

matlab m_map绘制中国行政区划 political region

[复制链接]
发表于 2013-4-23 20:08:34 | 显示全部楼层 |阅读模式
Question:
我想画三个省的行政边界,然后在这个范围内做等值线,我现在需要怎么呢?
Answer:
把m_proj中的范围设置为关注区域,先绘行政区划,hold on,然后再m_contourf

其实很想学mapping工具箱,这些功能似乎都是附带的,但是没找到机会,所以还是m_map吧
绘制行政区划一直不会,四处搜索,找到一个解决办法,但并不完美
这个办法优点:可以完整绘制全国各级行政区划以及行政单位
缺点,无法单独挑选单个省,也不能上色,只能画线

具体说来
行政区划资料来源: 国家基础地理信息系统 http://nfgis.nsdi.gov.cn/nfgis/chinese/c_xz.htm
登记姓名、单位、工作、邮箱之后就可以开始下载,感觉没什么验证过程,难道胡乱登记也可以么?
既然登记了,就把能下载的都下了,我们需要的省级区划是bou2_4m.shp
程序来源...是个论坛,忘记是哪位了,程序第一行有作者

依然是使用m_map工具箱,
如果想关注三个省,可以在m_proj命令中设置特定区域,当然区域内的省外界线还是会显示
shp文件是shape文件,用matlab的shaperead命令就可以读取,
以本程序为例,读出一组1785×1的构架数组map_path,每组代表了一小段界线,
如果能确定每个省的是哪部分,倒是可以把省勾出来,测试了一下,没发现什么规律,作罢。

% dacheng@NUIST
% draw map of China
map_path = shaperead('ch_map/bou2_4m/bou2_4l.shp'); % set path of .shp files
map_X = [map_path(.X];   % read x of the contour of province
map_Y = [map_path(.Y];   % read y of the contour of province
m_proj('Miller','lon',[70, 137],'lat',[0,60])   % set projection, and range
m_plot(map_X,map_Y) % draw political line
m_grid % draw grid
1.jpg
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|物探论坛 ( 鄂ICP备12002012号 微信号:iwutan )

GMT+8, 2024-4-19 23:07 , Processed in 0.079825 second(s), 22 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表