利用无人机对一片区域进行测绘前,我们会先在地图上框选一个区域,然后再规划飞行的路线,而需要测绘的这片区域往往是一个多边形。在 Mesh iOS 即将推出的新版本中,我们加入了多边形区域的编辑功能,其中就涉及判断用户所编辑出来的多边形是否合法的问题。

    首先我们要确定一个标准:怎么样才算一个不合法的多边形 ?我们可以简单地通过下面这幅图来解释一下:

    image.png
    我们可以看出前面两个分别是凹多边形和凸多边形,而最后一张则是我们所说的不合法多边形,可以看出这个不合法的多边形的特征就是:它存在某条边与另外一条边相交的情况

    那么要判断一个多边形是否合法,我们只要判断组成多边形的所有线段是否存在相交的情况即可,当然,我们这里所说的相交是 规范相交 ,即 交点不在线段的端点上

    好了,那么现在的问题可以简化成:如何判断两条线段是否规范相交

    这里我们需要借助 向量的叉积 来进行判断。

    叉积,又称向量积,是对三维空间中的两个向量的二元运算。

    这里推荐 3Blue1Brown 的 视频 来快速回顾一下叉积的概念(下面的两幅截图来自此视频)。我们只需知道叉积的结果是有正负的,比如我们以向量 如何判断一个多边形是否合法 - 图2 为标准,如下图,向量 如何判断一个多边形是否合法 - 图3如何判断一个多边形是否合法 - 图4顺时针方向,那么 如何判断一个多边形是否合法 - 图5

    image.png
    如何判断一个多边形是否合法 - 图7

    如果向量 如何判断一个多边形是否合法 - 图8如何判断一个多边形是否合法 - 图9逆时针方向,那么 如何判断一个多边形是否合法 - 图10
    image.png
    如何判断一个多边形是否合法 - 图12

    我们如何利用叉积的特性运用到判断线段是否相交上呢?

    我们先看下面最直接的一个线段相交的情况:

    image.png
    线段 如何判断一个多边形是否合法 - 图14和 线段 如何判断一个多边形是否合法 - 图15 明显存在一个交点,从上面这张图我们可以做一个简单的结论:如果一条的线段的两个端点在另外一条线段两侧,那么这两条线段可能相交,注意这里说的是可能相交,稍后会讲到另外一种情况。

    我们可以将上面的图转换为向量的情况来看:

    image.png

    是不是觉得似曾相识,这跟上面提到的叉积的情况是不是很类似?
    向量 如何判断一个多边形是否合法 - 图17如何判断一个多边形是否合法 - 图18 的逆时针方向,那么:如何判断一个多边形是否合法 - 图19
    向量 如何判断一个多边形是否合法 - 图20如何判断一个多边形是否合法 - 图21 的顺时针方向,那么:如何判断一个多边形是否合法 - 图22

    用 A 表示 如何判断一个多边形是否合法 - 图23 的叉积结果,用 B 表示 如何判断一个多边形是否合法 - 图24 的叉积结果,那么 一条的线段的两个端点在另外一条线段两侧 这个几何现象可以用这个公式表示 :A*B < 0

    我们前面提到 如果一条的线段的两个端点在另外一条线段两侧,那么这两条线段可能相交 ,为什么是可能相交呢?如果我们将 线段 如何判断一个多边形是否合法 - 图25 往右边移动一下,会存在下面这种情况:
    image.png
    从上图可以看出,线段 如何判断一个多边形是否合法 - 图27 的两个端点在线段 如何判断一个多边形是否合法 - 图28 两侧,但是它们并没有相交。

    那么如何排除这种情况呢?其实很简单,我们之前都是以线段 如何判断一个多边形是否合法 - 图29 作为主视角,如果将主视角换成线段 如何判断一个多边形是否合法 - 图30,那么我们很容易看出 线段 如何判断一个多边形是否合法 - 图31 的两个端点并没有在 线段如何判断一个多边形是否合法 - 图32 的两侧。所以我们再次看回上面相交的那幅图,为了能够充分的判断两条线段相交,这次以 如何判断一个多边形是否合法 - 图33 为主视角看待这个问题,求叉积:
    image.png
    向量 如何判断一个多边形是否合法 - 图35如何判断一个多边形是否合法 - 图36 的逆时针方向,那么:如何判断一个多边形是否合法 - 图37
    向量 如何判断一个多边形是否合法 - 图38如何判断一个多边形是否合法 - 图39 的顺时针方向,那么:如何判断一个多边形是否合法 - 图40

    综上,我们可以得出:
    如何判断一个多边形是否合法 - 图41
    如何判断一个多边形是否合法 - 图42
    如何判断一个多边形是否合法 - 图43
    如何判断一个多边形是否合法 - 图44

    A*B < 0 && C*D < 0 的时候,两条线段规范相交。
    至于向量的叉积如何运算,这里就不细写了,给出一张计算草稿给大家过目一下:
    image.png
    根据计算草稿的内容,我们就很容易通过代码来实现了:

    1. private func isIntersect(line1: (CGPoint, CGPoint), line2: (CGPoint, CGPoint)) -> Bool {
    2. let p1 = line1.0
    3. let p2 = line1.1
    4. let q1 = line2.0
    5. let q2 = line2.1
    6. let a1 = (p2.x - p1.x) * (q1.y - p1.y) - (q1.x - p1.x) * (p2.y - p1.y)
    7. let a2 = (p2.x - p1.x) * (q2.y - p1.y) - (q2.x - p1.x) * (p2.y - p1.y)
    8. let b1 = (q2.x - q1.x) * (p1.y - q1.y) - (p1.x - q1.x) * (q2.y - q1.y)
    9. let b2 = (q2.x - q1.x) * (p2.y - q1.y) - (p2.x - q1.x) * (q2.y - q1.y)
    10. if a1 * a2 < 0 && b1 * b2 < 0 {
    11. return true
    12. }
    13. return false
    14. }

    由于笔者能力有限,文中如有错误还请各位读者不吝赐教。