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

时间:2021-05-22

我们可以使用Basemap这个工具包来实现中国地图的绘制

首先需要加载一些包:

import numpy as npimport matplotlib.pyplot as pltfrom 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 npimport matplotlib.pyplot as pltfrom mpl_toolkits.basemap import Basemapplt.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 绘制地图的资料请关注其它相关文章!

声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。

相关文章