WGS84 与 UTM 互转（Python代码版）

妖狐艹你老母 2022-08-30 09:10 610阅读 0赞

## 前言 ##

本文不讲原理，只分享实践思路和代码。基于Python语言，使用pyproj库进行WGS 和UTM的转换。

**目录**

一、安装pyproj

二、WGS84 转 UTM-示例代码

三、UTM 转 WGS84-示例代码

四、算法效率

三、查询坐标系WKID

四、pyproj官网文档学习

--------------------

## 一、安装pyproj ##

pip install pyproj

## 二、WGS84 转 UTM-示例代码 ##

'''
    WGS84的经纬度 转 UTM的x,y
    '''
    from pyproj import Transformer
    
    
    # 参数1：WGS84地理坐标系统 对应 4326 
    # 参数2：坐标系WKID 广州市 WGS_1984_UTM_Zone_49N 对应 32649
    transformer = Transformer.from_crs("epsg:4326", "epsg:32649") 
    
    lat = 22.744435950
    lon = 113.595417400
    x, y = transformer.transform(lat, lon)
    print("x:", x, "y:", y)

输出结果：

x: 766544.7801555358  y: 2517564.4969607797

对比结果精准性，Apollo6.0代码转换的结果：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70][]

## 三、UTM 转 WGS84-示例代码 ##

'''
    UTM的x,y 转 WGS84的经纬度
    '''
    from pyproj import Transformer
    
    # 参数1：坐标系WKID 广州市 WGS_1984_UTM_Zone_49N 对应 32649
    # 参数2：WGS84地理坐标系统 对应 4326
    transformer = Transformer.from_crs("epsg:32649", "epsg:4326")
    
    x = 766544.7801555358
    y = 2517564.4969607797
    lat, lon = transformer.transform(x, y)
    print("x:", x, "y:", y)
    print("lat:", lat, "lon:", lon)

输出：x: 766544.7801555358           y: 2517564.4969607797  
          lat: 22.744435950000007        lon: 113.5954174

## 四、算法效率 ##

这里计算1s内，能运行多少次转换，以UTM 转 WGS84为例子。

'''
    计算1s内，能运行多少次转换（UTM的x,y 转 WGS84的经纬度）
    '''
    from pyproj import Transformer
    import time
    # 参数1：坐标系WKID 广州市 WGS_1984_UTM_Zone_49N 对应 32649
    # 参数2：WGS84地理坐标系统 对应 4326
    transformer = Transformer.from_crs("epsg:32649", "epsg:4326")
    
    sum_t=0.0            #花费的总时间
    i = 0
    
    # 计算1s内，能运行多少次转换。
    while sum_t < 1:
        time_start = time.time()  # 开始计时
        print("i:", i)
        i = i + 1
        ###########
        x = 766544.7801555358
        y = 2517564.4969607797
        lat, lon = transformer.transform(x, y)
        print("x:", x, "y:", y)
        print("lat:", lat, "lon:", lon)
        ###########
        time_end = time.time()                # 结束计时
        sum_t=(time_end - time_start)+sum_t   #运行所花时间
        print('time cost', sum_t, 's')
        print("\n")

这是我电脑的结果，1s能转换2489次左右。![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 1][]

## 三、查询坐标系WKID ##

可以从以下两个网站进行查询：

地理坐标系WKID：[https://developers.arcgis.com/javascript/3/jshelp/gcs.htm][https_developers.arcgis.com_javascript_3_jshelp_gcs.htm]

投影坐标系WKID：[https://developers.arcgis.com/javascript/3/jshelp/pcs.htm][https_developers.arcgis.com_javascript_3_jshelp_pcs.htm]

比如，查询广州市的，先查广州市的UTM区号WGS\_1984\_UTM\_Zone\_49N；然后在官方网页，用Ctrl + F搜索一下，能找到WGS\_1984\_UTM\_Zone\_49N 对应 32649

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 2][]

## 四、pyproj官网文档学习 ##

[https://www.osgeo.cn/pyproj/examples.html][https_www.osgeo.cn_pyproj_examples.html]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 3][]

参考：[https://blog.csdn.net/tap880507/article/details/111608529][https_blog.csdn.net_tap880507_article_details_111608529]

[UTM坐标和WGS84坐标（如何转换？][UTM_WGS84]

[https://blog.csdn.net/ywcpig/article/details/107516253][https_blog.csdn.net_ywcpig_article_details_107516253]

本文只提供参考学习，谢谢。

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70]: /images/20220829/3364b0f97ab6404192b4d333ba17d055.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 1]: /images/20220829/1f39eb63bd55485690502f362a536e00.png
[https_developers.arcgis.com_javascript_3_jshelp_gcs.htm]: https://developers.arcgis.com/javascript/3/jshelp/gcs.htm
[https_developers.arcgis.com_javascript_3_jshelp_pcs.htm]: https://developers.arcgis.com/javascript/3/jshelp/pcs.htm
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 2]: /images/20220829/9d8c00c586cf47c4bae541256920f733.png
[https_www.osgeo.cn_pyproj_examples.html]: https://www.osgeo.cn/pyproj/examples.html
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMjA0NDY0_size_16_color_FFFFFF_t_70 3]: /images/20220829/ac5b36b5f0dc46c08baa195e27e06dec.png
[https_blog.csdn.net_tap880507_article_details_111608529]: https://blog.csdn.net/tap880507/article/details/111608529
[UTM_WGS84]: https://blog.csdn.net/qq_41204464/article/details/104543636?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522162668837016780274133032%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fblog.%2522%257D&request_id=162668837016780274133032&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~first_rank_v2~rank_v29-1-104543636.pc_v2_rank_blog_default&utm_term=UTM&spm=1018.2226.3001.4450
[https_blog.csdn.net_ywcpig_article_details_107516253]: https://blog.csdn.net/ywcpig/article/details/107516253