コマンドラインから、初期オフセットと最終オフセットを指定してサブラインを検索します

Aug 19 2020

基本的に次のようなCSVファイルがあります。

id,id_wkt,start,end,wkt
11,11,0,43,"LINESTRING(3.3623676 51.3132714,3.3626663 51.3130068,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.3643705 51.3117777)"
12,11,43,217,"LINESTRING(3.3623676 51.3132714,3.3626663 51.3130068,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.3643705 51.3117777)"

これは、WKT行文字列を持つレコードで構成されます。一部のWKTは複数回表示され、実際のジオメトリは、それぞれの列の開始オフセットと終了オフセット(ラインに沿ったメートル単位の長さ)を使用してラインストリングから取得する必要があります。

これらの長さのオフセットに基づいてサブラインストリングを抽出したいと思います。これは、OGRなどのコマンドラインツールを介して可能ですか?

回答

Socci Aug 18 2020 at 22:08

いくつかの調査の後、私は実際に以下を使用して解決策を見つけましたogr2ogr

ogr2ogr output.csv -lco GEOMETRY=AS_WKT input.csv -dialect SQLITE -sql "SELECT id, id_wkt, ST_Line_Substring(GEOMETRY, start / GeodesicLength(GEOMETRY), end / GeodesicLength(GEOMETRY)) FROM input"

それは与えます

WKT,id,id_wkt
"LINESTRING (3.3623676 51.3132714,3.3626663 51.3130068,3.36274075608485 51.3129464307024)","11","11"
"LINESTRING (3.36274075608485 51.3129464307024,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.36436609040842 51.3117804449745)","12","11"