【shapely】【投影变换】

# shapely.ops.transform做投影变换，可以跳过geopandas

from shapely.ops import transform

from functools import partial 

import pyproj

 

line_pnts = [113.428126, 34.808386], [113.540136, 34.808702]

proj = partial(pyproj.transform, pyproj.Proj(init="EPSG:4326"),pyproj.Proj(init="EPSG:4527"))

lpts = transform(proj,LineString(line_pnts))

print(lpts.wkt)

 

posted @ 2021-06-30 17:04  ddzhen  阅读(349)  评论(0)  编辑  收藏  举报