基于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)
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 开发者必知的日志记录最佳实践
· 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 让容器管理更轻松!