Python 基于GDAL包的一些地理坐标投影转换 包括GDAL的安装以及环境配置设置

1.安装GDAL包

建议https://www.lfd.uci.edu/~gohlke/pythonlibs/#_gdal找对应的python版本下载之后安装

pip install GDAL‑3.4.2‑cp37‑cp37m‑win_amd64.whl

然后就是要设置一下环境变量

重启计算机之后再进行一些函数的调用就可以了

2.投影以及坐标转换等

def getSRSPair(dataset):
    '''
    获得给定数据的投影参考系和地理参考系
    :param dataset: GDAL地理数据
    :return: 投影参考系和地理参考系
    '''
    prosrs = osr.SpatialReference()
    prosrs.ImportFromWkt(dataset.GetProjection())
    geosrs = prosrs.CloneGeogCS()
    return prosrs, geosrs

def geo2lonlat(dataset, x, y):
    '''
    将投影坐标转为经纬度坐标（具体的投影坐标系由给定数据确定）
    :param dataset: GDAL地理数据
    :param x: 投影坐标x
    :param y: 投影坐标y
    :return: 投影坐标(x, y)对应的经纬度坐标(lon, lat)
    '''
    prosrs, geosrs = getSRSPair(dataset)
    ct = osr.CoordinateTransformation(prosrs, geosrs)
    coords = ct.TransformPoint(x, y)
    return coords[:2]

def lonlat2geo(datas