如何从与Polygon相交的Shapely LineString中提取相交线段端点

你不需要拆分原始LineString为多个小线段，可以直接利用Shapely的intersection方法计算线与多边形的交集，从交集结果中提取实际相交部分的端点。以下是适配你代码的解决方案：

修改后的完整代码 #

from shapely.geometry import LineString, Polygon
import matplotlib.pyplot as plt

intersections_list = []

x_line_patch =  [0, 4, -2, 8, 3, 10, 5]
y_line_patch =  [21, 17, 14, 11, 9, 6, 1]

x_line_patch2 =  [x + 20 for x in x_line_patch]
y_line_patch2 =  y_line_patch

# 创建多边形
start_x = 5
end_x = 8
start_y = 5
end_y = 8
polygon = Polygon([(start_x, start_y), (end_x, start_y), (end_x, end_y), (start_x, end_y)])
lines = [LineString(zip(x_line_patch, y_line_patch)), LineString(zip(x_line_patch2, y_line_patch2))]

# 处理相交线，提取相交线段的端点
for line in lines:
    if line.intersects(polygon):
        # 计算线与多边形的交集
        intersection = line.intersection(polygon)
        # 根据交集的几何类型提取端点
        if hasattr(intersection, 'geoms'):
            # 处理MultiLineString或MultiPoint
            for geom in intersection.geoms:
                if isinstance(geom, LineString):
                    # 提取LineString的两个端点
                    coords = list(geom.coords)
                    intersections_list.extend(coords)
                # 如果是Point类型，可根据需求决定是否加入
                # elif isinstance(geom, Point):
                #     intersections_list.append(list(geom.coords)[0])
        elif isinstance(intersection, LineString):
            # 单个LineString，直接取两端点
            coords = list(intersection.coords)
            intersections_list.extend(coords)
        elif intersection.geom_type == 'Point':
            # 线与多边形相切或端点在多边形内的情况，按需处理
            intersections_list.append(list(intersection.coords)[0])

# 可视化验证
x, y = polygon.exterior.xy
plt.plot(x, y)
plt.fill(x, y, alpha=0.3)

# 绘制原始线
x, y = lines[0].xy
x2, y2 = lines[1].xy
plt.plot(x, y)
plt.plot(x2, y2)

# 绘制提取的相交端点
if intersections_list:
    inter_x = [p[0] for p in intersections_list]
    inter_y = [p[1] for p in intersections_list]
    plt.scatter(inter_x, inter_y, color='red', s=50, zorder=10)

plt.tight_layout()
plt.show()

# 输出结果
print("相交线段端点列表：", intersections_list)

关键说明 #

1. 直接计算交集：使用line.intersection(polygon)可以直接得到线与多边形的重叠部分，避免了拆分原始线的额外复杂度。
2. 处理多种交集类型：
   • 当线穿过多边形时，交集是LineString，直接提取其两端点即可
   • 如果线与多边形有多个不连续的相交部分，会返回MultiLineString，需要遍历每个子线提取端点
   • 若线的端点在多边形内或线与多边形相切，会返回Point，可根据业务需求决定是否保留这类点
3. 适配你的数据结构：不需要获取原始的x/y坐标列表，直接对已有的LineString对象操作即可，完全符合你实际代码中的场景。

内容的提问来源于stack exchange，提问作者Chris Ze Third