【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)

原文地址:https://www.cnblogs.com/ddzhen/p/14955389.html