python proj utm2wgs wgs2utm

Posted 2022-06-23 zolty

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)