GDAL经纬度转84坐标系

在Python中使用GDAL库进行经纬度坐标转换为84坐标系（WGS84）的示例代码如下：

from osgeo import ogr, osr

# 定义转换目标和源空间参考
source = osr.SpatialReference()
source.ImportFromEPSG(4326)  # EPSG代码4326表示经纬度坐标系

target = osr.SpatialReference()
target.ImportFromEPSG(4326)  # EPSG代码4326表示WGS84坐标系

# 创建转换对象
transform = osr.CoordinateTransformation(source, target)

# 定义经纬度坐标点
lon = 116.3975  # 经度
lat = 39.9087  # 纬度

# 创建点对象
point = ogr.Geometry(ogr.wkbPoint)
point.AddPoint(lon, lat)

# 进行坐标转换
point.Transform(transform)

# 获取转换后的坐标
lon84 = point.GetX()
lat84 = point.GetY()

print(f"转换后的坐标为: ({lon84}, {lat84})")

这段代码使用GDAL库中的`osr`模块来进行坐标转换，首先定义了目标空间参考（WGS84）和源空间参考（经纬度坐标系），然后创建了一个`CoordinateTransformation`对象用于进行坐标转换。接下来，定义了一个经纬度坐标点，并使用`Transform`方法进行转换，最后通过`GetX`和`GetY`方法获取转换后的坐标。请注意替换`lon`和`lat`变量为你需要转换的经纬度坐标。