点多边形 - 故障算法

时间:2015-05-04 18:59:23

标签: c polygon point point-in-polygon

我有一个alogrithm的C代码,检查给定的点是否在多边形内。它应该是正确的,我也在不同的地方看到这个代码。然而,当我使用它并不能完美地工作时 - 大约20%的答案是错误的。

int pnpoly(int nvert, double *vertx, double *verty, double testx, double 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;
}

我的主要功能可能有问题。可以给我一个检查这个算法的主要功能吗?

这是我的主要功能

int main(){

    double vertx[4] = {10, 10, 0, 0};
    double verty[4] = {10, 0, 10, 0};

    // for those two it returns "Inside"
    double testx = 6;
    double testy = 4;

    /* for those two it returns "Outside"
    double testx = 5;
    double testy = 4;
    */

    int result = pnpoly(4, vertx, verty, testx, testy);

    if (result) {
        printf("\nInside\n");
    }

    else {
        printf("\nOutside\n");
    }

    return 0;
}

1 个答案:

答案 0 :(得分:1)

您的多边形是自相交的。 (5,4)是“外部”是正常的

我认为你认为你的多边形是一个正方形,即使是自相交的多边形,算法也能很好地工作。