两点截取线【shapely】【clipline】

复制代码
from shapely.geometry import Point, LineString
from shapely.ops import snap, split, nearest_points
def clipLine(ls, p1, p2):
    p1 = nearest_points(ls, p1)[0]
    geom = snap(ls, p1, 1)
    parts = split(geom, p1).geoms  # 新版本的shapely对该函数的返回值做了改动
    if len(parts) == 1:
        fp = parts[0]
    else:
        fp = parts[0] if parts[0].distance(p2) < parts[1].distance(p2) else parts[1]

    p2 = nearest_points(fp, p2)[0]
    geom = snap(fp, p2, 0.1)
    parts = split(geom, p2).geoms #

    if len(parts) == 1:
        res = parts[0]
    else:
        res = parts[0] if parts[0].distance(p1) < parts[1].distance(p1) else parts[1]
    return res


if __name__ == '__main__':

    ls = LineString([(0, 0), (5, 5)])
    p1 = Point((1, 1))
    p2 = Point([3, 3])
    res = clipLine(ls, p1, p2)
    print(res.wkt) # LINESTRING (1 1, 3 3)
复制代码

 

posted @   ddzhen  阅读(598)  评论(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 让容器管理更轻松!
点击右上角即可分享
微信分享提示