(1)首先绘制一张全球地图:
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
plt.figure(figsize=(16,8))
map1 = Basemap()
map1.drawcoastlines()
map1.drawcountries(linewidth=1.5)
plt.show()
(2)绘制中国地图
由于basemap中未包含中国行政区划,需自己下载。下载地址:https://gadm.org/download_country_v3.html,其中包含全世界任何一个国家的行政区划shape文件。
plt.figure(figsize=(16,8))
map2 = Basemap(llcrnrlon=77,llcrnrlat=14,urcrnrlon=140,urcrnrlat=51,projection='lcc', lat_1=33, lat_2=45, lon_0=100)
map2.drawcoastlines()
map2.drawcountries(linewidth=1.5)
map2.readshapefile('gadm36_CHN_shp/gadm36_CHN_1','states',drawbounds=True)
plt.show()