能否利用Cartopy识别本地Shapefile的地理坐标系？

当然可以通过Cartopy轻松实现这个需求！Cartopy的shpreader.Reader对象不仅能读取几何数据，还能自动解析Shapefile的坐标系信息，下面是具体的实现方法：

方法一：直接读取Reader的crs属性 #

Cartopy会自动识别Shapefile配套的.prj文件，将其中的坐标系信息转换为Cartopy专属的CRS对象，你可以直接通过reader.crs访问：

import cartopy.io.shapereader as shpreader

# 初始化Shapefile读取器
reader = shpreader.Reader('/path_to_shapefile/myShapeFile.shp')

# 获取并打印坐标系信息
shape_crs = reader.crs
print("Shapefile的地理坐标系:", shape_crs)
# 若需要查看对应的Proj4格式字符串，可执行：
print("对应的Proj4字符串:", shape_crs.proj4_init)

方法二：手动解析.prj文件（兼容特殊格式） #

如果遇到.prj文件格式较特殊导致reader.crs无法正常识别的情况，你也可以手动读取文件内容再转换为CRS对象：

from cartopy import crs as ccrs
import pyproj

# 读取.prj文件的WKT格式内容
with open('/path_to_shapefile/myShapeFile.prj', 'r') as f:
    prj_wkt = f.read()

# 将WKT内容转换为Cartopy可用的CRS对象
proj_crs = pyproj.CRS.from_wkt(prj_wkt)
shape_crs = ccrs.Projection(proj_crs)

print("手动解析的坐标系:", shape_crs)

整合到你的绘图代码中 #

你可以把坐标系识别逻辑和原绘图代码结合，确保绘图时使用与Shapefile匹配的投影参数：

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

reader = shpreader.Reader('/path_to_shapefile/myShapeFile.shp')
shape_crs = reader.crs

# 使用识别出的坐标系创建绘图轴（替代固定的OSGB投影）
ax = plt.axes(projection=shape_crs)
ax.add_geometries(reader.geometries(), shape_crs, edgecolor='k', facecolor='none')

# 添加坐标系信息标注
plt.title(f"Shapefile坐标系: {shape_crs.name}")
plt.show()

这样既完成了坐标系的识别，也能避免因投影不匹配导致的绘图错误。

内容的提问来源于stack exchange，提问作者user1718097