PostGIS補間/三角測量オプション

DEM構築を自動化するオープンソースワークフローを探しています。LIDARデータセットに含まれる一連のサイトがあります。各サイトにサイト固有のDEMを作成し、プロセスを自動化することを検討しています。

これまでのところ、自動化されています：

• PostGISへのデータの読み込み（LIDARおよびサイトポイント）
• サイト境界の作成（st_bufferとst_envelopeの組み合わせ）

現在、各サイトでデータサブセットを補間し、Surfer7グリッドファイルにエクスポートするオプションを探しています。

現在、QuantumGISを使用してPostGISデータベースのデータをクエリし、csvファイルにエクスポートしてから、Surfer7でデータを手動でインポートおよびグリッド化します。うまくいけば、これも自動化できます。

そのため、このユースケースでは、LIDARデータをDEMに三角化する必要があります。考えられる他の潜在的なケースには、さまざまな補間方法が関係します。したがって、逆距離とクリギングのオプションがある場合、興味があります。

これは私たちにとって非常に学習の練習です-私たちは赤ちゃんのステップで働いています！

この自動化は、WindowsまたはLinuxのSAGAで実行できます。SAGAは、ポイントクラウド形式（spc）としてSAGA形式でインポートされたビッグポイントデータセットを処理できます。次に、補間または三角測量のための関数があり、Surferにエクスポートします。自動化は、SAGAコマンドツールを使用して実行されます。

.las / lazに元のLIDARデータがある場合、lastools、http： //www.cs.unc.edu/~isenburg/lastools/を使用して、大きなLiDARデータを処理できます。直接動作するように見えるlas2tinとlas2demがありますが、私はそれらを使用していません。以下を使用して、一連の.lazファイルをGRASS内のラスターDEMにアセンブルしました。

cat *.laz |~/Downloads/SRC/lastools/bin/las2txt -stdin  -parse xyztc -header pound -otxt |invproj -v -f %.9f +init=epsg:2284 |sed -e 's// /' | r.in.xyz input=- output=lidarJunk method=min type=FCELL x=1 y=2 z=3 fs=\    zscale=1.0 percent=100 --overwrite  
# the sed and fs=\ bits are for dealing with reformatting spaces and tabs between invproj and r.in.xyz

Lastoolsは適切に自動化され、DEMまたはTINのデータを選択した範囲にクリップできるように見えます。