解析几何:计算两条线段的交点

开发 前端
求两线段的交点,本质就是解方程,需要用到克莱姆法则,计算出来的交点是直线交点,不一定是线段交点,需要再判断点是否在线段范围内。

大家好,我是前端西瓜哥。

今天来实现计算两条线段的交点的解析几何算法。

我们要实现 getLineSegIntersection 方法:提供两条线段,计算它们的交点。

每条线段会用两个点坐标表示。

const getLineSegIntersection = (p1, p2, p3, p4) => {
  // 待实现
}

// 测试用例
getLineSegIntersection(
  { x: 1, y: 1 }, { x: 4, y: 4 },
  { x: 1, y: 4 }, { x: 4, y: 1 }
);
// 期望 { x: 2.5, y: 2.5 }

思路

思路很简单,就是解两条直线对应的一个二元一次方程组,求出 x 和 y。

如果无解或多解,说明直线平行,交点不存在。

如果有解,可拿到唯一交点,但也只能说明直线有交点,还需要判断线段是否有交点。

所以我们需要判断交点是否在线段的区间上。如果是,说明两线段有交点,返回交点。

克拉姆法则

解方程组需要用到 克拉姆法则。

对于:

可转换为矩阵形式表示:

然后计算主矩阵(最左边的矩阵)的行列式,对角相乘然后相减:

如果行列式为 0,说明没有唯一解。

如果不为 0,则有唯一解:

回到我们的两条直线,我们用两点式表示直线:

转换成 Ax+By=C 的格式,得到:

于是:

const a = y2 - y1;
const b = x1 - x2;
const c = x1 * y2 - x2 * y1;

第二条线段同理:

const d = y4 - y3;
const e = x3 - x4;
const f = x3 * y4 - x4 * y3;

算法实现

interface Point {
  x: number;
  y: number;
}

const getLineSegIntersection = (
  p1: Point,
  p2: Point,
  p3: Point,
  p4: Point
): Point | null => {
  const { x: x1, y: y1 } = p1;
  const { x: x2, y: y2 } = p2;
  const { x: x3, y: y3 } = p3;
  const { x: x4, y: y4 } = p4;

  const a = y2 - y1;
  const b = x1 - x2;
  const c = x1 * y2 - x2 * y1;

  const d = y4 - y3;
  const e = x3 - x4;
  const f = x3 * y4 - x4 * y3;

  // 计算分母
  const denominator = a * e - b * d;

  // 判断分母是否为 0(代表平行)
  if (Math.abs(denominator) < 0.000000001) {
    // 这里有个特殊的重叠但只有一个交点的情况,可以考虑处理一下
    return null;
  }

  const px = (c * e - f * b) / denominator;
  const py = (a * f - c * d) / denominator;

  // 判断交点是否在两个线段上
  if (
    px >= Math.min(x1, x2) &&
    px <= Math.max(x1, x2) &&
    py >= Math.min(y1, y2) &&
    py <= Math.max(y1, y2) &&
    px >= Math.min(x3, x4) &&
    px <= Math.max(x3, x4) &&
    py >= Math.min(y3, y4) &&
    py <= Math.max(y3, y4)
  ) {
    return { x: px, y: py };
  }

  return null;
};

变体

这个算法可以做一些变体,实现其他的算法。

变体1:两线段是否有交点。

返回值换成布尔值即可。

变体2:计算两直线的交点。

把判断直线交点是否在线段上的逻辑去掉,然后直接返回点坐标即可。

优化点

重叠但却只有一个交点的情况。

如果线段平行,有两种情况:

  • 没有重叠(0 个解)
  • 有部分重叠(多解)

如果部分重叠,可能有多个点,多个点的情况下也不知道拿哪个点作为交点好,这种情况下还是返回 null。

但有一个特殊的情况:重叠只有一个点(比如线段 a 的末点刚好是线段 b 的起点)。如果你的场景下判断比较严格,你可以选择返回这个点。要实现这部分也是有点点复杂的。

误差处理。线段的两个端点的距离非常小,计算出的结果也会非常小,可能会进入了 0 的绝对误差范围了,考虑改成相对误差。

溢出风险。数值很大时有溢出风险,可以考虑计算一个缩放值,缩小后计算,计算完再放大回去。

结尾

总结一下,求两线段的交点,本质就是解方程,需要用到克莱姆法则,计算出来的交点是直线交点,不一定是线段交点,需要再判断点是否在线段范围内。

不复杂,就是有一点点小细节。

责任编辑:姜华 来源: 前端西瓜哥
相关推荐

2023-08-30 07:34:41

2012-12-20 10:19:08

华为路由器接入设置

2023-04-05 14:31:19

Java计算移动

2019-04-04 13:36:25

云计算互联网云厂商

2011-06-07 11:21:34

路由负载

2009-06-19 15:25:13

ITSMNSM运维管理

2014-10-24 15:17:07

Android

2014-12-24 09:15:54

PaaS开源云服务

2012-05-11 13:15:12

戴尔虚拟化

2022-06-06 23:22:44

互联网产品模式

2019-11-06 15:16:12

16GB8GB内存

2011-06-21 15:24:20

HP高密度计算服务器

2021-11-10 23:44:21

笔记本触摸板技巧

2016-08-18 09:53:33

软件定义存储

2022-08-11 13:11:48

斯坦福大学英伟达VR 头显

2010-01-12 18:05:56

Linux Redha

2017-01-11 15:45:52

中国联通光缆亚欧5号

2018-07-30 11:03:08

4GB内存8GB

2009-03-06 12:17:24

IBMPowerSystemPower6

2021-12-08 09:00:25

LeetCode容器算法
点赞
收藏

51CTO技术栈公众号