如何判定Shapely中被LineString拆分的几何左右侧？

要判断拆分后的两个多边形（Aside和Bside）对应splitter行走方向的左右侧，核心思路是通过多边形的内部点（比如质心）与splitter线段的相对位置来判定——利用向量叉乘的符号就能直接区分左右，完全符合你定义的“沿LineString节点顺序行走时的左右手侧”规则。

原理说明 #

在Shapely使用的笛卡尔坐标系（y轴向上）中，假设splitter的起点为P1(x1,y1)，终点为P2(x2,y2)，待判断点为P(x,y)：

• 计算向量叉乘：cross = (x2 - x1) * (y - y1) - (y2 - y1) * (x - x1)
  • 如果cross > 0：点P在splitter的左侧（沿P1→P2方向行走时）
  • 如果cross < 0：点P在splitter的右侧

因为Aside和Bside是被splitter完全分隔的两个多边形，它们的质心（centroid）必然分别在splitter的两侧，所以用质心作为判断点非常可靠。

实现get_LRsides函数 #

替换你原来的get_LRsides函数为以下代码：

def get_LRsides(Aside, Bside, splitter):
    """Determine the 'left' and 'right' sides of an already split geometry"""
    # 获取splitter的起点和终点坐标
    p1 = splitter.coords[0]
    p2 = splitter.coords[-1]
    
    # 计算两个多边形的质心
    centroid_a = Aside.centroid.coords[0]
    centroid_b = Bside.centroid.coords[0]
    
    # 定义叉乘计算函数
    def calculate_cross(p):
        return (p2[0] - p1[0]) * (p[1] - p1[1]) - (p2[1] - p1[1]) * (p[0] - p1[0])
    
    # 计算两个质心的叉乘值
    cross_a = calculate_cross(centroid_a)
    cross_b = calculate_cross(centroid_b)
    
    # 根据叉乘符号判断左右
    if cross_a > 0:
        Lside, Rside = Aside, Bside
    else:
        Lside, Rside = Bside, Aside
    
    return Lside, Rside

关键细节说明 #

• 为什么用质心？：多边形的质心是其几何中心，一定位于多边形内部，不会因为多边形的边缘顶点靠近splitter而导致误判，能稳定代表整个区域的位置。
• 坐标系兼容性：Shapely默认的笛卡尔坐标系（y轴向上）完全匹配叉乘的左右判断逻辑，不需要额外转换坐标系。
• 边界情况处理：因为splitter已经将minimum_rotated_rectangle拆分为两个独立部分，质心不可能落在splitter上（否则splitter会穿过质心，导致拆分后的多边形包含质心，但实际拆分后是两个分离的多边形），所以无需处理cross=0的极端情况。

验证逻辑 #

这个函数可以直接集成到你的splitLR函数中，替换原有的get_LRsides调用即可。测试时，你可以手动构造一个简单的旋转矩形和splitter，比如：

# 测试示例
rect = geo.Polygon([(0,0), (2,1), (3,3), (1,2)])
splitter = geo.LineString([(0,2), (3,0)])
left, right = splitLR(rect, splitter)
# 可视化验证：left应该是沿splitter从(0,2)到(3,0)行走时的左侧区域

内容的提问来源于stack exchange，提问作者Rick