複数の二次元多角形領域を調整修正（互いに素で、平面を余さず被覆）

こんにちは。
「複数の二次元多角形領域(polygon)が互いに素で、平面を余さず被覆」の集合となるよう調整修正を施す方法を探してみました。

目的は polygon同士のギャップや重なりを無くすことです。これらは検出が容易ではないですが、結合させてみるとギャップの存在が分かります¹。

polygon の集合は GeoJSON 形式を仮定しました（入出力ファイル）。

• ogr2ogr (GDAL) で、ST_MakeValidを施し、各 polygon の self-intersection を修正する。

$ ogr2ogr -f GeoJSON -explodecollections output.json input.json -dialect sqlite -sql "select ST_MakeValid(geometry) as geometry from OGRGeoJSON"

• pprepair (planar partition repair) ：入力形式は自動的に各種に対応しますが、出力形式は今のところ、 cmake の引数で指定した一種類です。入力ファイル内で "properties": {} のように空になっていると segment fault になります。

$ cmake -DDRIVER="GeoJSON" .
$ make
$ ./pprepair -h
$ ./pprepair -i input.json -v
$ ./pprepair -i input.json -o output.json -rrlb

---

1. 実例は、「結合でできた「隙間」を排除する」。 ↩