shapely相关操作

Sheyueyu

已于 2023-06-13 15:26:37 修改

阅读量270

点赞数

分类专栏：交通 shapely 文章标签： python

于 2023-06-13 15:24:09 首次发布

本文链接：https://blog.csdn.net/sheyueyu/article/details/131189269

版权

交通同时被 2 个专栏收录

8 篇文章 0 订阅

订阅专栏

shapely

2 篇文章 0 订阅

订阅专栏

shapely相关操作

1. 将为折线的linestring分为多段直线

line = LineString([(114.18453629400005, 30.507459477000054),
 (114.18823357300005, 30.50530951700017),
 (114.19034520000002, 30.50882250000012)])
line

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-uSeAJGKV-1686640879316)(./shapely相关操作.assets/image-20230613113435281.png)]

multi_list = []
for i in range(len(line.coords[:-1])):
    start, end = line.coords[i], line.coords[i+1]
    ls = LineString([start,end])
    multi_list.append(ls)
multi_list

[<shapely.geometry.linestring.LineString at 0x1f68a3f07f0>,
 <shapely.geometry.linestring.LineString at 0x1f68a0829e0>]

2. 延伸linestring

延伸前：

在这里插入图片描述

延伸后：

在这里插入图片描述

函数:(注意投影)

def extend_line_in_end(line, distance=10.2):
    """
    延长线段的两端
    ------
    line(linestring):投影过的线段
    distance(float):在两端延长的的长度
    """
    if not isinstance(line, LineString):
        raise ValueError("Input should be a LineString object")

    # 获取起点和终点
    start, end = line.coords[0], line.coords[-1]
    
    # 计算方向
    dx = end[0] - start[0]
    dy = end[1] - start[1]
    angle = math.atan2(dy, dx)

    # 计算新的起点和终点
    new_start = (start[0] - distance*math.cos(angle), start[1] - distance*math.sin(angle))
    new_end = (end[0] + distance*math.cos(angle), end[1] + distance*math.sin(angle))

    # 创建新的LineString
    new_line = LineString([new_start] + list(line.coords)[1:-1] + [new_end])
    
    return new_line