osgeo和pyproj：经纬度坐标和高斯坐标互相转换

最新推荐文章于 2025-07-21 11:04:49 发布

原创

最新推荐文章于 2025-07-21 11:04:49 发布 · 2.1k 阅读

17 ·

CC 4.0 BY-SA版权

本文为博主原创文章，未经博主允许请勿转载！作者：ISmileLi 博客地址：http://blog.youkuaiyun.com/toby54king。

文章标签：

#gdal #经纬度转高斯 #高斯转经纬度 #坐标转换 #python

这篇博客介绍了在地图处理和自动驾驶中如何使用Python的osgeo和pyproj库进行经纬度到高斯坐标的转换。在gdal V3.0以后，转换操作需要设置转换策略以避免出现‘inf’结果。博主提供了经过测试的转换代码，并分享了在github上的下载链接。

一、前言
搞地图和自动驾驶的都知道，坐标转换是非常频繁的事情，有时候需要在各种坐标之间来回的转换，最近使用python代码处理地图数据，在使用osgeo库中的gdal时，发现了gdal v2和V3的一些不同之处，研究了一下，这里分享出来。

二、问题描述
经纬度转高斯的过程中，发现3.0一直出现的转换结果是’inf’，经过查看官方github上的issue，才知道，gdal V3.0以后，转换需要设置转换策略，具体看后面代码中的说明，现象截图如下：
在这里插入图片描述

三、解决后封装的代码
下面代码使用osgeo库和pyproj库分别实现，底层调用都是gdal，代码经过测试。

#!/usr/bin/python3

__author__ = 'ISmileLi'

from osgeo import gdal, ogr, osr
from pyproj import Transformer

'''
osgeo底层坐标转换使用的库还是proj,下面函数中的espg值需要根据自己的需求进行修改，
下文测试使用的是wgs84与中国区高斯-克吕格EPSG码为21460区的转换
'''

def lonLat_to_gauss(lon, lat, from_epsg=4326, to_epsg=21460):
    '''
    经纬度转高斯
    :param lon:
    :param lat:
    :param from_epsg:
    :param to_EPSG:
    :return:
    '''

    from_spa = osr.SpatialReference()
    '''
    gdal版本大于3.0以后必须设置转换策略才能正确显示结果，否则结果将会输出'inf'
    可以了解官方的这个issue说明：https://github.com/OSGeo/gdal/issues/1546
    '''
    if int(gdal.__version__[0]) >=