如何用python在地图上显示shapefile？

Question

我目前正在使用python，用basemap绘制越南地图，就像这样。越南地图

我的形状文件: 请在此输入图片描述

我想用我已经有的shapefile来显示所有的岛群或小岛，但我不知道如何读取shapefile。我在网上试过一些代码，但是没有用。希望你们能帮助我。我真的很感激。非常感谢。这是我的代码。

import cartopy.crs as ccrs
from cartopy.feature import ShapelyFeature
from cartopy.io.shapereader import Reader
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
plt.figure(figsize=(12,7))
m = Basemap(projection='mill',
       llcrnrlat = 7,
       urcrnrlat = 25,
       llcrnrlon = 90,
       urcrnrlon = 120,
       resolution = 'l')
m.drawcoastlines()
m.drawcountries()
m.drawstates()
m.fillcontinents(color='None',lake_color='#FFFFFF')
m.drawmapboundary(color='k', linewidth=1.0, fill_color=None, zorder=None, ax=None)
parallels = np.arange(0.,81,10.)
# labels = [left,right,top,bottom]
m.drawparallels(parallels,labels=[False,True,True,False])
meridians = np.arange(10.,351.,10.)
m.drawmeridians(meridians,labels=[True,False,False,True])

fig, ax = plt.subplots(subplot_kw=dict(projection=ccrs.PlateCarree()))
ax.set_extent([90, 120, 25, 7])
# Set map extent
In_st_shp = 'E:\koppen-master\haidao\haidao.shp'
# provide path of shapefile
state_feature = ShapelyFeature(Reader(In_st_shp).geometries(),
                               crs=ccrs.PlateCarree(), edgecolor='k')
ax.add_feature(state_feature, facecolor="none")

plt.show()

Answer 1

这是使用cartopy的工作代码.请在代码中检查地图范围和设置shapefile的路径。

import shapefile
import cartopy.crs as ccrs
from cartopy.feature import ShapelyFeature
from cartopy.io.shapereader import Reader
import matplotlib.pyplot as plt
fig, ax = plt.subplots(subplot_kw=dict(projection=ccrs.PlateCarree()))
ax.set_extent([7, 23, 102, 111])
# Set map extent
In_st_shp = 'path_to_shape.shp'
# provide path of shapefile
state_feature = ShapelyFeature(Reader(In_st_shp).geometries(),
                               crs=ccrs.PlateCarree(), edgecolor='k')
ax.add_feature(state_feature, facecolor="none")