多边形光线投射算法中的点答案

【问题标题】：Point in polygon ray casting algorithm多边形光线投射算法中的点
【发布时间】：2014-12-09 21:01:46
【问题描述】：

给定一个大小为 M、N 的矩阵和一个点列表形成一个多边形，我需要返回一个二进制矩阵，它标记多边形内的所有点。代码如下：使用水平线扫描，与多边形相交并将相交之间的所有点标记为真。

import cv2
import math
import numpy as np

def getABC(x1, y1, x2, y2):
    A = y2 - y1
    B = x1 - x2
    C = A*x1 + B*y1
    return (A, B, C)

def polygon(M, N, poly):
    """
    Return a mask matrix
    Where points inside the polygon is 1, outside is 0
    """
    out = np.zeros((M, N)). astype(bool)

    n = len(poly)
    for i in range(M):

        # horizontal scanning
        intersection_x = i
        intersection_y = []

        # check through all edges
        for edge in range(n + 1):
            v1_x, v1_y = poly[edge % n]
            v2_x, v2_y = poly[(edge + 1) % n]

        v1_x = int(v1_x)
        v1_y = int(v1_y)
        v2_x = int(v2_x)
        v2_y = int(v2_y)

            assert (v1_x <= M and v2_x <= M and v1_y <= N and v2_y <= N)

            A1, B1, C1 = getABC(v1_x, v1_y, v2_x, v2_y)
            A2 = 1
            B2 = 0
            C2 = i

            # find intersection
            if intersection_x >= min(v1_x, v2_x) and intersection_x <= max(v1_x, v2_x):
                det = A1*B2 - A2*B1
                if (det != 0):
                    tmp = (A1 * C2 - A2 * C1)/det
                    intersection_y.append(math.ceil(tmp))

        intersection_y = sorted(list(set(intersection_y)))



        if len(intersection_y) > 1:
            for k in range(1, len(intersection_y), 2):
                out[intersection_y[k - 1]:intersection_y[k], intersection_x] = True

    return out

if __name__ == "__main__":
out = polygon(512, 512,         
    [[
      321.929203539823, 
      414.4070796460177
    ], 
    [
      164.74414597854175, 
      512.0
    ], 
    [
      0.0, 
      509.9846825337252
    ], 
    [
      0.0, 
      221.4867256637168
    ], 
    [
      170.60176991150445, 
      2.902654867256672
    ], 
    [
      320.1592920353982, 
      91.39823008849561
    ], 
    [
      271.4867256637168, 
      201.1327433628319
    ], 
    [
      348.4778761061947, 
      228.56637168141594
    ], 
    [
      359.98230088495575, 
      302.9026548672567
    ], 
    [
      220.15929203539824, 
      329.4513274336283
    ]])

    cv2.imwrite("polygon.png", out * 255)

在拐角的情况下执行错误，例如上面的输入会生成：

我无法解决这个问题。根据交叉点的数量，我将分配 l0-l1 l2-l3 l4-l5 或 l1-l2 l3-l4 l5-l6 为真。

【问题讨论】：

标签： python algorithm computer-science corner-detection

【解决方案1】：

您可以通过应用以下算法进行修复：
1. 如果一条线与它内部的一条边相交（不是在端点之一），请始终保持相交点。
2.如果一条线在端点与一条边相交，当且仅当它在边的左端（x坐标较小）时才保留它。
3. 如果边缘是垂直的，则忽略它。

之后，只需对点进行排序而不删除重复项。

【讨论】：