如何在不使用第三方模块的情况下求两条线段的交点（若存在）

你提到的两个Stack Overflow问题分别是《如何计算两条直线的交点？》和《如何判断两条线段是否相交？》，但现有答案要么无法稳定生效，要么依赖第三方库，以下是纯原生实现的解决方案，无需任何外部依赖。

核心思路 #

用向量叉积判断线段的相对位置关系，结合参数方程计算交点坐标，全程仅使用基础数学运算，兼容所有线段相交场景。

关键概念 #

叉积：用于判断点与线段的位置关系，以及两条线段是否平行。对于点A(xa, ya)、B(xb, yb)、C(xc, yc)，叉积公式为：
cross = (xb - xa) * (yc - ya) - (yb - ya) * (xc - xa)

• 若cross > 0：点C在线段AB的左侧
• 若cross < 0：点C在线段AB的右侧
• 若cross = 0：点C在AB所在的直线上

计算步骤 #

1. 计算两条线段的方向向量，通过叉积判断是否平行或重合
2. 若不平行，计算交点对应的参数值，验证参数是否在[0,1]区间内（确保交点在线段上）
3. 处理特殊场景：端点相交、线段部分/完全重叠

代码实现 #

def segment_intersection(seg1, seg2):
    # 解析线段端点：seg1 = ((x1,y1), (x2,y2)), seg2 = ((x3,y3), (x4,y4))
    (x1, y1), (x2, y2) = seg1
    (x3, y3), (x4, y4) = seg2

    # 计算方向向量
    dx1 = x2 - x1
    dy1 = y2 - y1
    dx2 = x4 - x3
    dy2 = y4 - y3

    # 分母：判断两条线段是否平行
    denominator = dx1 * dy2 - dy1 * dx2

    # 处理平行或重合的情况
    if denominator == 0:
        # 辅助函数：判断点是否在线段上
        def point_on_segment(px, py, sx1, sy1, sx2, sy2):
            cross = (sx2 - sx1) * (py - sy1) - (sy2 - sy1) * (px - sx1)
            if cross != 0:
                return False
            return (min(sx1, sx2) <= px <= max(sx1, sx2)) and (min(sy1, sy2) <= py <= max(sy1, sy2))
        
        # 收集所有重合的端点
        intersection_points = []
        if point_on_segment(x3, y3, x1, y1, x2, y2):
            intersection_points.append((x3, y3))
        if point_on_segment(x4, y4, x1, y1, x2, y2):
            intersection_points.append((x4, y4))
        if point_on_segment(x1, y1, x3, y3, x4, y4):
            intersection_points.append((x1, y1))
        if point_on_segment(x2, y2, x3, y3, x4, y4):
            intersection_points.append((x2, y2))
        
        # 去重并返回结果
        if not intersection_points:
            return None  # 平行无交点
        unique_points = list(set(intersection_points))
        if len(unique_points) == 1:
            return unique_points[0]
        else:
            unique_points.sort()
            return tuple(unique_points)
    else:
        # 计算参数t和u，对应两条线段上的位置
        t_numerator = (x3 - x1) * dy2 - (y3 - y1) * dx2
        u_numerator = (x3 - x1) * dy1 - (y3 - y1) * dx1
        t = t_numerator / denominator
        u = u_numerator / denominator

        # 参数在[0,1]区间内，说明交点在线段上
        if 0 <= t <= 1 and 0 <= u <= 1:
            ix = x1 + t * dx1
            iy = y1 + t * dy1
            return (ix, iy)
        else:
            return None  # 交点在线段延长线上，无有效交点

使用示例 #

# 示例1：线段相交于内部点
seg1 = ((0,0), (2,2))
seg2 = ((0,2), (2,0))
print(segment_intersection(seg1, seg2))  # 输出 (1.0, 1.0)

# 示例2：线段交于端点
seg1 = ((0,0), (2,2))
seg2 = ((2,2), (3,3))
print(segment_intersection(seg1, seg2))  # 输出 (2, 2)

# 示例3：平行无交点
seg1 = ((0,0), (1,1))
seg2 = ((0,1), (1,2))
print(segment_intersection(seg1, seg2))  # 输出 None

# 示例4：线段部分重叠
seg1 = ((0,0), (3,3))
seg2 = ((1,1), (2,2))
print(segment_intersection(seg1, seg2))  # 输出 ((1, 1), (2, 2))

结果说明 #

函数会返回三种类型的结果：

• None：两条线段无交点（包括平行不重叠、交点在延长线上）
• 单个坐标元组：线段交于单点（内部点或端点）
• 包含两个坐标元组的元组：两条线段部分/完全重叠，返回重叠区间的端点

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