使用GDAL进行简单的坐标系转换

使用GDAL进行简单的坐标系转换

使用python GDAL进行简单的坐标系转换，暂时不考虑不同基准坐标系转换的精度问题。

安装环境

使用Ubuntu+Anaconda python 环境
conda install gdal

定义坐标系

from osgeo import gdal
from osgeo import osr

srs_wgs84 = osr.SpatialReference()
srs_wgs84.ImportFromEPSG(4326)


srs_beijing = osr.SpatialReference()
srs_beijing.ImportFromEPSG(2418)

print(srs_wgs84)

print(srs_beijing)

定义坐标转换Transform

xmin=42386600
ymin=4663300
xmax=42393800
ymax=4670900

ct = osr.CoordinateTransformation(srs_beijing, srs_wgs84)
min_coords = ct.TransformPoint(ymin, xmin)
print(min_coords)

max_coords = ct.TransformPoint(ymax, xmax)
print(max_coords)

参考

1、https://blog.csdn.net/mrbaolong/article/details/104212918