使用gdal的TransformPoints()方法可以将地理坐标系和投影坐标系下的坐标点相互转换,如WGS84转到UTM投影,最近在实际使用的时候发现,同样使用该方法,gdal2.0和3.0得到的结果坐标顺序不一致。测试如下
from osgeo import osr
#point1 at UTM zone 50N
PointULX=699960
PointULY=4200000
#point2 at UTM zone 50N
PointBRX=809760
PointBRY=4090200
source