利用python绘制中国地图（含省界、河流等）

更新时间：2020年09月21日 15:40:34 作者：王琛NUIST

这篇文章主要介绍了利用python绘制中国地图，帮助大家更好的理解和使用python绘制图形，感兴趣的朋友可以了解下

我们可以使用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()

以上就是利用python绘制中国地图的详细内容，更多关于python 绘制地图的资料请关注脚本之家其它相关文章！

您可能感兴趣的文章:

Python入门教程(三十七)Python中的删除文件
这篇文章主要介绍了Python入门教程(三十七)Python中的删除文件,在Python中,如果需要删除文件,必须导入OS模块,使用os模块的remove()方法,需要的朋友可以参考下
2023-05-05
Python日志模块logging的使用方法总结
这篇文章主要分享的是Python日志模块logging的使用方法总结，ogging模块默认级别是WARNING，意味着只会追踪该级别以上的事件，除非更改日志配置，想了解更多相关资料的小伙伴可以参考下面文章内容
2022-05-05
深入理解Python虚拟机中调试器实现原理与源码分析
本文主要给大家介绍python中调试器的实现原理，通过了解一个语言的调试器的实现原理我们可以更加深入的理解整个语言的运行机制，可以帮助我们更好的理解程序的执行，感兴趣的可以了解一下
2023-04-04
python机器学习基础K近邻算法详解KNN
这篇文章主要为大家介绍了python机器学习基础K近邻算法详解有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2021-11-11
关于django 1.10 CSRF验证失败的解决方法
今天小编就为大家分享一篇关于django 1.10 CSRF验证失败的解决方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2019-08-08
在脚本中单独使用django的ORM模型详解
这篇文章主要介绍了在脚本中单独使用django的ORM模型详解，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2020-04-04
记一次django内存异常排查及解决方法
这篇文章主要给大家介绍了关于一次django内存异常排查记解决方法，文中通过示例代码介绍的非常详细，对大家学习或者使用django具有一定的参考学习价值，需要的朋友们下面来一起学习学习吧
2020-08-08
python中的迭代和可迭代对象代码示例
这篇文章主要介绍了python中的迭代和可迭代对象代码示例，具有一定借鉴价值,需要的朋友可以参考下
2017-12-12
Python基于gevent实现文件字符串查找器
这篇文章主要介绍了Python基于gevent实现文件字符串查找器,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
2020-08-08
使用python实现kmean算法
这篇文章主要介绍了使用python实现kmean算法,kmean 是无监督学习的一种算法，主要是用来进行聚类分析的，他会在数据集中算出几个点作为簇中心，求这些数据集与这些簇中心的距离，并将距离同一个簇中心距离最近的数据归为一类,需要的朋友可以参考下
2023-04-04