python proj utm2wgs wgs2utm

最新推荐文章于 2024-10-08 18:16:12 发布

zolty

最新推荐文章于 2024-10-08 18:16:12 发布

阅读量1.2k

点赞数

分类专栏： python 文章标签： python 开发语言后端

本文链接：https://blog.csdn.net/qq_41854291/article/details/122612075

版权

python 专栏收录该内容

25 篇文章

订阅专栏

本文介绍如何使用Python进行WGS84坐标系与UTM坐标系之间的转换，包括wgs2utm和utm2wgs两个方向的转换，并提供具体的实现代码示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1.wgs2utm

zone 需要选择

2.utm2wgs

# zolty add v into vtd by udp
from pyproj import Proj, Transformer


# lon,lat = (119.9579946997, 30.5206528196)
def wgs2tm(lon, lat):
    pipeline_str = (
        '+proj=utm +zone=50 +datum=WGS84 +units=m +no_defs'
        '+proj=longlat +datum=WGS84 +no_defs'
    )
    pipe_trans = Transformer.from_pipeline(pipeline_str)
    xt, yt = pipe_trans.transform(lon, lat)
    print("utm no diff: ", xt, yt)
    #  diff
    xt = xt - 785883
    yt = yt - 3380679
    print("utm + diff : ", xt, yt)


def wgs2utm(lon, lat):
    p = Proj(proj='utm', zone=50, ellps='WGS84')
    x, y = p(lon, lat)
    print("utm no diff: ", x, y)


def utm2wgs(x, y):
    p = Proj(proj='utm', zone=50, ellps='WGS84')
    lon, lat = p(x, y, inverse=True)
    print("lon , lat  : ", lon, lat)


if __name__ == '__main__':
    wgs2tm(119.9579946997, 30.5206528196)
    wgs2utm(119.9579946997, 30.5206528196)
    utm2wgs(785883, 3380679)