在cartopy地图顶部绘制shapefile城市边界

import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
import matplotlib.pyplot as plt

fig = plt.figure(figsize=(10, 10))
ax = plt.axes(projection=ccrs.PlateCarree())

shapefile_path = 'path_to_shapefile'  # 替换为shapefile文件的路径
shapefile = shpreader.Reader(shapefile_path)

for record in shapefile.records():
    geometry = record.geometry
    ax.add_geometries([geometry], ccrs.PlateCarree(), edgecolor='black', facecolor='none')

plt.show()