基于pyproj地理坐标和投影坐标转换

复制代码
from shapely.geometry import Point, LineString
from pyproj import Transformer
from pyproj import CRS


class wgs842proj(): def __init__(self, fcrs=None, tcrs=None): if fcrs is None: fcrs = CRS("EPSG:4326") if tcrs is None: tcrs = CRS("EPSG:4527") self.transformer = Transformer.from_crs(crs_4326, crs_proj) def pointProj(self, p): # p Point(113, 34) lng, lat = p.x, p.y y, x = self.transformer.transform(lat, lng) return Point([x, y]) def lineProj(self, line): # line LineString() lngs, lats = line.xy pts = [self.transformer.transform(lat, lng)for lat, lng in zip(lats, lngs)] pts = [_[::-1]for _ in pts] return LineString(pts) if __name__ == '__main__': p0 = [113.0, 34.0] p1 = [113.5, 34.0] ls = LineString([p0, p1]) p0 = Point(p0) p1 = Point(p1) op = wgs842proj() p0_ = op.pointProj(p0) p1_ = op.pointProj(p1) ls_ = op.lineProj(ls)
复制代码

 

posted @   ddzhen  阅读(627)  评论(0编辑  收藏  举报
编辑推荐:
· 开发者必知的日志记录最佳实践
· SQL Server 2025 AI相关能力初探
· Linux系列:如何用 C#调用 C方法造成内存泄露
· AI与.NET技术实操系列(二):开始使用ML.NET
· 记一次.NET内存居高不下排查解决与启示
阅读排行:
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· 没有Manus邀请码?试试免邀请码的MGX或者开源的OpenManus吧
· 【自荐】一款简洁、开源的在线白板工具 Drawnix
· 园子的第一款AI主题卫衣上架——"HELLO! HOW CAN I ASSIST YOU TODAY
· Docker 太简单,K8s 太复杂?w7panel 让容器管理更轻松!
点击右上角即可分享
微信分享提示