可以使用 Cartopy 识别 shapefile 的地理坐标系吗？

Question

我有一个已下载到本地硬盘上的 shapefile。我可以使用 Cartopy 读取 shapefile 并使用 Matplotlib 绘制地图，如下所示：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
reader = shpreader.Reader('/path_to_shapefile/myShapeFile.shp')
ax = plt.axes(projection=ccrs.OSGB())
ax.add_geometries(reader.geometries(),ccrs.OSGB(),edgecolor = 'k',facecolor = 'none')
plt.show()

但是，我希望能够识别 shapefile 的地理坐标系（例如 WGS84 或 OSGB），以便确保它与我想要绘制的某些点的坐标值相匹配在同一张地图上。阅读 shapefile 后，我假设 Cartopy 可以显示使用的坐标系，但我不知道如何做到这一点。 Cartopy 可以这样做还是我需要使用不同的包？

Answer 1

Cartopy 在后台使用 pyshp 来读取 shapefile。据我所知，该图书馆不提供您正在寻找的信息。 （CartoPy 无论如何都不会公开有关 PyShp 的这些信息。）

您最好的选择可能是查看Fiona，这是 GDAL 的不错的 Python 绑定。您可以使用它来手动打开 shapefile 并进行检查——我假设这个库将提供您正在寻找的信息，如果 Python 空间中有任何内容。