Python 省和省会词典 python 绘制中国省界图

import shapefile
import numpy as np
from mpl_toolkits.basemap import Basemap

file=shapefile.Reader('gadm36_CHN_shp/gadm36_CHN_1.shp')  #读取省界.shp文件
shapes=file.shapes()   #获取point 
records=file.records() #获取省名称

pro_points=[]  #建立省边界列表
pro_names=[] #建立省名称列表

for i in range(len(shapes)):
    points=shapes[i].points #h获取经纬度数据    
    pro_name=file.records()[i][3]  #获取省名称
    lon =[]  
    lat  =[] 
    #将每个tuple的lon和lat组合起来
    [lon.append(points[i][0]) for i in range(len(points))]  
    [lat.append(points[i][1]) for i in range(len(points))]   
 
    lon=np.array(lon).reshape(-1,1)
    lat=np.array(lat).reshape(-1,1)
    loc=np.concatenate((lon,lat),axis=1)
    
    pro_points.append(loc)
    pro_names.append(pro_name)

lat_min,lat_max=37,55
lon_min,lon_max=115,135

fig = plt.figure(figsize=(18,10))
ax1 = fig.add_axes([0.1,0.1,0.8,0.8])

m=Basemap(projection='cyl',llcrnrlat=lat_min,llcrnrlon=lon_min,
       urcrnrlat=lat_max,urcrnrlon=lon_max,resolution='l',ax=ax1)
       
m.drawcoastlines()
m.drawcountries()

name_list=['Liaoning','Jilin','Beijing','Hebei','Heilongjiang']  #画出这几个省的边界

for name,point in zip(pro_names,pro_points):
    
    if name in name_list:
        
        lon=point[:,0]
        lat=point[:,1]
        
        plt.scatter(lon,lat,marker='.',c='k',s=0.5)

import h5py

file_name='province_boundary_lon_lat.hdf'  #保存的目标文件名和路径
f=h5py.File(file_name,mode='w')  #创建一个hdf文件

for pro_name,points in zip(pro_names,pro_points):
   #每一个省份创建一个群组 
    a=f.create_group(pro_name) 
    a['longitude']=points[:,0]
    a['latitude']=points[:,1]
 
f.close() #写入完成后，记得关掉。