利用python绘制中国地图(含省界、河流等)
作者:王琛NUIST 时间:2021-12-19 03:19:59
我们可以使用Basemap这个工具包来实现中国地图的绘制
首先需要加载一些包:
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
Basemap包就是气象画图的利器,现在我们就可以愉快的画图了!
plt.figure(1)
map=Basemap()
map.drawcoastlines()
plt.title(r'$World\ Map$',fontsize=24)
plt.show()
第2行创建一个地图,第3行添加海岸线,这样一个世界地图就出来了,怎么样,很简单吧。(plt.show()这行代码是用来显示图片的)
我们发现这只是海岸线图,那么怎么将国界线添加上去呢?很简单,只要添加一行代码就可以了。
map.drawcountries()
那么怎么添加河流呢?可能有些同学已经猜到了,就是drawrivers()
map.drawrivers(color='blue',linewidth=0.3)
好了,现在我们可以开始画中国地图了!
其实只要在创建地图时指定一下范围就可以了,查阅资料发现,中国的经纬度范围是东经135度2分30秒-东经73度40分,北纬3度52分-北纬53度33分。
map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54)
好了,一个中国地图就出来了!但是我们发现,好像少了点什么,没错就是省界。我们可以在https://gadm.org/download_country_v3.html下载中国大陆和台湾省的行政区域的shape文件,下载后解压,然后加入下面的代码。
CHN='G:\python_material\MapOfChina'
CHN的值就是解压后的地图文件所在的地址。
下面我们就可以加入省界了!
map.readshapefile(CHN+'\gadm36_CHN_shp\gadm36_CHN_1',
'states',drawbounds=True)
别忘了把台湾省加上去
map.readshapefile(CHN+'\gadm36_TWN_shp\gadm36_TWN_1',
'taiwan',drawbounds=True)
还可以在地图上加上经纬度,比如我们要画5条经纬线,可以这么做:
parallels = np.linspace(3,55,5)
map.drawparallels(parallels,labels=[True,False,False,False])
meridians = np.linspace(70,140,5)
map.drawmeridians(meridians,labels=[False,False,False,True])
大功告成!但是,emmm,我们发现好像有点歪?我们可以在创建地图时选择投影参数。
map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=137,urcrnrlat=54,
projection = 'lcc', lat_1 = 33, lat_2 = 45, lon_0 = 100)
这回正式的完成了!
附上所有代码:
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
plt.figure(1)
map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54,
projection = 'lcc', lat_1 = 33, lat_2 = 45, lon_0 = 100)
map.drawcoastlines()
map.drawcountries()
map.drawrivers(color='blue',linewidth=0.3)
CHN='G:\python_material\MapOfChina'
map.readshapefile(CHN+'\gadm36_CHN_shp\gadm36_CHN_1',
'states',drawbounds=True)
map.readshapefile(CHN+'\gadm36_TWN_shp\gadm36_TWN_1',
'taiwan',drawbounds=True)
parallels = np.linspace(3,55,5)
map.drawparallels(parallels,labels=[True,False,False,False])
meridians = np.linspace(70,140,5)
map.drawmeridians(meridians,labels=[False,False,False,True])
plt.title(r'$China\ Map$',fontsize=24)
plt.show()
来源:https://www.jianshu.com/p/9d0100f455bb
标签:python,绘制,中国,地图
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
去掉前面的0的sql语句(前导零,零前缀)
2011-09-30 11:28:19
Python中输入和输出(打印)数据实例方法
2021-03-25 22:32:11
![](https://img.aspxhome.com/file/2023/7/127617_0s.png)
如何使用ASP来读写注册表
2007-09-20 13:08:00
用python将pdf转化为有声读物
2021-02-02 10:06:46
![](https://img.aspxhome.com/file/2023/0/105030_0s.jpg)
oracle 日期函数集合(集中版本)第1/2页
2009-06-19 17:23:00
PyQt5固定窗口大小的方法
2021-05-18 12:34:43
![](https://img.aspxhome.com/file/2023/0/99040_0s.jpg)
运行(runCode)复制(copyCode)保存(saveCode)代码框方法
2007-10-21 08:41:00
举例讲解Python的lambda语句声明匿名函数的用法
2021-05-29 13:02:53
详解用python -m http.server搭一个简易的本地局域网
2023-09-29 13:03:34
![](https://img.aspxhome.com/file/2023/7/110587_0s.png)
Python分析彩票记录并预测中奖号码过程详解
2023-07-20 04:49:18
![](https://img.aspxhome.com/file/2023/0/66000_0s.png)
Python给定一个句子倒序输出单词以及字母的方法
2022-11-21 00:36:50
Ajax的错误处理机制探讨
2007-09-07 09:53:00
让XML在ASP中发挥其长处
2008-01-16 19:07:00
PyQt5 关于Qt Designer的初步应用和打包过程详解
2023-07-31 20:26:32
![](https://img.aspxhome.com/file/2023/2/72132_0s.png)
python中星号变量的几种特殊用法
2021-03-20 10:57:50
对TensorFlow的assign赋值用法详解
2023-03-18 22:52:56
浅谈tensorflow语义分割api的使用(deeplab训练cityscapes)
2022-03-16 23:24:01
PHP PDOStatement::bindValue讲解
2023-06-14 11:16:02
深入认识javascript中的eval函数
2008-08-03 16:44:00
用python统计代码行的示例(包括空行和注释)
2022-06-28 02:15:30