获取相应坐标系
import gdal
import numpy as np
from osgeo import osr
dataset = gdal.Open("C:\\Users\\Nihil\\Desktop\\Arcgisdata\\HPDEM.tif")
adfGeoTransform = dataset.GetGeoTransform()
get_ref = dataset.GetSpatialRef()
print(get_ref)
PROJCS["WGS 84 / UTM zone 48N",
GEOGCS[