Question

因此，对于我的计算机图形学类，我的任务是做一个多边形填充器，目前我的软件渲染器是用Python编码的。现在，我想测试一下我在How can I determine whether a 2D Point is within a Polygon?上找到的pointInPolygon代码，以便以后可以根据自己的方法创建自己的方法。

代码是：

int pnpoly(int nvert, float *vertx, float *verty, float testx, float testy)
{
  int i, j, c = 0;
  for (i = 0, j = nvert-1; i < nvert; j = i++) {
    if ( ((verty[i]>testy) != (verty[j]>testy)) &&
     (testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]) )
       c = !c;
  }
  return c;
}

我在Python中重新创建它的尝试如下：

def pointInPolygon(self, nvert, vertx, verty, testx, testy):
    c = 0
    j = nvert-1
    for i in range(nvert):
        if(((verty[i]>testy) != (verty[j]>testy)) and (testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i] + vertx[i]))):
           c = not c

        j += 1

    return c

但是显然，这将在第二次迭代中返回超出范围的索引，因为j = nvert并将崩溃。

谢谢。

Answer 1

您正在错误地读取棘手的C代码。 df.groupby(['date', 'column-3']).size()的意义是使j = i++递增一个，并且将 old 值分配给i。在循环结束时，类似的python代码也会执行j：

j = i

我们的想法是，对于j = nvert - 1 for i in range(nvert): ... j = i，这些值会出现

nvert == 3

另一种实现方法是j | i ---+--- 2 | 0 0 | 1 1 | 2等于j，

(i - 1) % nvert

即它滞后一个，并且索引形成一个环（就像顶点一样）

更多的Python代码将使用for i in range(nvert): j = (i - 1) % nvert ...并遍历坐标本身。您将有一个称为 vertices 的对（元组）列表，以及两个迭代器，其中一个是一个顶点在另一个顶点之前，并且由于以下原因循环回到起点itertools，类似：

itertools.cycle

如何在Python的多边形代码中实现这一点？

1 个答案: