無限に尽きるボロノイ多角形?
私はvoronoi.pyスクリプト(Bill SimonによるSteve Fortuneのコードの翻訳)の実装に取り組んでいます。QGISリンクによれば、それは基本的に私が作業しているコードと同じですが、無限に拡張できるポリゴンを生成する方法は2つで異なります。QGISの場合、サイトポイント範囲を超えた範囲で境界矩形が生成されると想定しています。私の場合、興味のあるサイトポイントを囲む円として偽のサイトポイントを生成します。サークルポイントの半径は、サイトポイントの範囲の直径の2倍です。2つの間の結果を添付の図に示します。左側にQGIS実装、右側に私の実装です。 。 ご覧のとおり、QGIS実装はサイトポイントのない奇妙なポリゴンを生成します(両方の実装の例としてサイト473を参照)。これは明らかに間違っていますが、実装が正しいかどうかもわかりません。私はArcMapのArcInfoライセンスにアクセスできないため、ArcMapがネイティブに生成するものと比較することはできず、テストするためにサイトポイントシェープファイルをアタッチする方法がわかりません。したがって、エッジポリゴンと「無限の長方形」を強制するのに役立つ「無限の円」の使用についてのコメントがあれば、感謝します。さらに良いことに、他のGISプラットフォームで出力をテストすることに興味がある人がいれば、使用したサイトポイントのシェープファイルを提供できます。 問題の一部のように見えるのは、サイトポイントを超えるポリゴンがオーバーラップしているため、塗りつぶしなしでシンボル化すると、異なるポリゴンとして表示されることです。次の図は、サイト473および415についてのポリゴンが選択されたときに何が起こるかを示しています。したがって、元の質問は残っていますが、新しいしわがあります。 これはバグ(または「不完全性」)であり、0%のバッファーが使用されている場合、QGISで正しい結果になります。まだ解決する必要があるのは、円形の「無限大ポイント」の密集したパターンが、長方形の「無限大エッジ」と同じ結果をもたらすかどうかです。