如何从GeoDataFrame.crs返回结果中提取EPSG编码

你可以直接调用CRS对象的内置方法拿到EPSG编码，有两种常用实现方式：

• 直接获取数值型EPSG编码
  调用to_epsg()方法，直接返回整数形式的编码，无对应EPSG时返回None：

  epsg_code = gdf.crs.to_epsg()
  # 示例输出：25833
  

  如果你需要降低匹配置信度阈值来获取最接近的EPSG，可以调整min_confidence参数，取值范围0~100，默认值为70：

  epsg_code = gdf.crs.to_epsg(min_confidence=20)
  

• 同时获取权威机构名称和编码
  调用to_authority()方法，会返回格式为(权威机构名, 编码)的元组，示例场景下返回('EPSG', '25833')，提取编码的写法如下：

  _, code = gdf.crs.to_authority()
  epsg_code = int(code)
  

注意：如果当前使用的是自定义坐标系，没有对应的官方EPSG编号，上述两个方法会返回None，属于正常现象。

内容的提问来源于stack exchange，提问作者four-eyes