4

ポイントがポリゴン内にあるかどうかを確認したいと思います。もちろん、私はグーグルで検索して、この質問が以前に回答されているかどうかを調べ、次のアルゴリズムを見つけました: http://www.ecse.rpi.edu/Homepages/wrf/Research/Short_Notes/pnpoly.html ポリゴンが部分的に開いていない限り、これは正常に機能します。 . 例えば： 例

AE は正常に検出されますが、B ポリゴンの開いた部分も閉じていると見なされます。このサンプル コードを実行すると、私の言いたいことがわかるでしょう。

#include <stdio.h>

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;
}

int main(int argc, char *argv[])
{
        // 1 closed [A]
        float x1[] = { 0, 1, 1, 0, 0 };
        float y1[] = { 0, 0, 1, 1, 0 };
        printf("1: %d (1 expected)\n", pnpoly(5, x1, y1, 0.8, 0.8));
        printf("1: %d (0 expected)\n", pnpoly(5, x1, y1, -0.8, -0.8));

        // 1 closed [B] with a partial open
        // please note that the vertex between [0,-1] and [0,0] is missing
        float x2[] = { 0, 1, 1, 0, 0, -1, -1, 0 };
        float y2[] = { 0, 0, 1, 1, 0,  0, -1, -1 };
        printf("2: %d (1 expected)\n", pnpoly(8, x2, y2, 0.8, 0.8));
        printf("2: %d (0 expected)\n", pnpoly(8, x2, y2, -0.8, -0.8)); // <- fails

        // 2 closed [C/D/E]
        float x3[] = { 0, 1, 1, 0, 0, -1, -1, 0, 0 };
        float y3[] = { 0, 0, 1, 1, 0,  0, -1, -1, 0 };
        printf("3: %d (1 expected)\n", pnpoly(9, x3, y3, 0.8, 0.8));
        printf("3: %d (1 expected)\n", pnpoly(9, x3, y3, -0.8, -0.8));

        return 0;
}

x2/y2 ポリゴンは、部分的に開いたブロックに接続された閉じたブロックで構成されます。pnpoly 関数は、開いたブロックの "内" の点を多角形内にあると見なします。

私の質問は次のとおりです。どうすればこの問題を解決できますか? それとも私は何かを見落としていますか？

前もって感謝します。

4

1 に答える 1

1

返事が遅くなるかもしれませんが、とにかく貢献したいと思います。

このようなアルゴリズムは最後のポイントを最初のポイントに接続するように扱うため、ポリゴン配列データで前処理を行う必要があります。特にこのアルゴリズムの説明でわかるように、ポリゴンを扱いますが、開いた領域はポリゴンではないため、このアルゴリズムでは扱うことができません。

次のアイデアを使用して、ポリゴンで前処理を行うことをお勧めします。

「多角形を表す可能性のある点のすべてのセットについて、開始点が繰り返されているのを見つけるたびに、サブ多角形に分割する必要があります。見つかったら、この多角形を保存し、この点から始めて、新しい多角形を作成し、セットの最後に到達するまで処理を続けます. したがって, 開始点で終わらないすべてのサブポリゴンは処理で無視されるべきであり, したがってアルゴリズムに提出されません. この前処理アプローチを使用することにより,終点が始点と同じサブポリゴンのみをアルゴリズムに送信します。」

お役に立てますように - それでも必要な場合は...

于 2014-06-27T00:37:43.123 に答える