Python获取栅格数据的空间坐标系统及EPSG编号
1、代码如下
from osgeo import gdal,osr
#栅格文件存在的地址
f2 = 'C:\\Users\\zijiancao\\Desktop\\osgeopy-data\\osgeopy-data\\osgeopy-data\\Utah\\cache.tif'
ds2 = gdal.Open(f2)
sr2 = osr.SpatialReference(wkt=ds2.GetProjection())
#获得坐标系,若遇到无空间参考系的,无输出
print(sr2)
#获得ESPG,
原创
2021-03-16 19:02:41 ·
1516 阅读 ·
0 评论