获取shp文件的空间参考信息

 针对带有.prj的shp文件，读取相关的空间参考信息。

# -*- coding:gb2312 -*-
from osgeo import osr,ogr

shp_path=r"C:\Users\skfzh\Documents\python地理空间分析指南\第五章\5.4\NYC_MUSEUMS_LAMBERT\NYC_MUSEUMS_GEO.shp"
ds=ogr.Open(shp_path)
utm_sr=ds.GetLayer().GetSpatialRef()
# print(utm_sr)
"""
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"]]
"""
print(utm_sr.GetAttrValue("GEOGCS"))
"""
WGS 84
"""