Advertisement

osgeo 和 pyproj:处理经纬度与高斯坐标的转换

阅读量:

一、前言
在地图处理和自动驾驶领域中,
无论是从事地图处理还是自动驾驶领域的人士都能清楚知道,
频繁进行坐标转换是一项非常常见的操作,
特别是在同一项目中可能需要在同一时间进行不同坐标系之间的转换。
最近在使用Python代码处理地理数据时,
发现了一项值得注意的功能特性,
即在利用osgeo库中的GDAL模块进行数据处理时,
发现了GDAL v2和v3版本之间的差异后,
并对其进行了深入研究,
在此分享了相关经验。

在进行经纬度转高斯的过程中,遇到了一个问题,在GDAL V3.0版本之后的操作中发现的问题是,在GDAL官方GitHub上的相关Issue中指出,在进行特定操作时会遇到计算异常的情况。具体操作时需要设置相应的转换策略以避免出现无效值(NaN/Inf)。此外,请参考附图中的现象截图以更好地理解问题所在。

在这里插入图片描述

该段代码基于osgeo库和pyproj库分别进行开发,并且所有底层操作都依赖于gdal库的支持,在测试阶段已经实现了稳定运行

复制代码
    #!/usr/bin/python3
    
    __author__ = 'ISmileLi'

全部评论 (0)

还没有任何评论哟~