GDAL经纬度转84坐标系python
时间: 2023-09-16 16:15:17 浏览: 143
在Python中使用GDAL库进行经纬度坐标转换为84坐标系(WGS84)的示例代码如下:
```python
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`变量为你需要转换的经纬度坐标。
阅读全文