使用python的GDAL库将坐标系的wkt转epsg

ylfmsn

已于 2023-06-15 10:22:02 修改

阅读量521

点赞数 1

分类专栏： gdal 文章标签： python

于 2023-06-15 08:16:50 首次发布

本文链接：https://blog.csdn.net/ylfmsn/article/details/131219976

版权

gdal 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

在Python中使用GDAL库通过WKT格式的坐标系数据获取坐标系的EPSG代码：

from osgeo import osr


def get_epsg_from_wkt(wkt):
    spatial_ref = osr.SpatialReference()
    spatial_ref.ImportFromWkt(wkt)
    epsg = spatial_ref.GetAuthorityCode(None)
    return int(epsg)


# WKT格式的坐标系数据
wkt_data = 'PROJCS["WGS 84 / UTM zone 31N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],AUTHORITY["EPSG","32631"],AXIS["Easting",EAST],AXIS["Northing",NORTH]]'

# 获取EPSG代码
epsg_code = get_epsg_from_wkt(wkt_data)
print(epsg_code)

gdal版本：3.2.1