GDAL坐标转换——TransformPoint
生活随笔
收集整理的這篇文章主要介紹了
GDAL坐标转换——TransformPoint
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
值得注意的是GDAL 3.x 改變了坐標策略,需要使用SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)來糾正
from osgeo import osr, __version__x = 359019 y = 3924437 prj = osr.SpatialReference() prj.ImportFromEPSG(32654) if __version__ >= "3.0.0":prj.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER) geo = osr.SpatialReference() geo.ImportFromEPSG(4326) if __version__ >= "3.0.0":geo.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)trans = osr.CoordinateTransformation(prj, geo) lng, lat, _ = trans.TransformPoint(x, y) print(lng, lat)trans1 = osr.CoordinateTransformation(geo, prj) x, y, _ = trans1.TransformPoint(lng, lat) print(x, y)總結
以上是生活随笔為你收集整理的GDAL坐标转换——TransformPoint的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: C++ 模板的显式具体化
- 下一篇: MCE公司:新型靶向RAS突变的共价抑制