LP(航空レーザ)のグリッドデータから、砂防フロンティアの区域設定アプリで変換できるTINテキストを作るプログラム
砂防フロンティアの区域設定支援システムに地形データを入れるため、航空レーザ測量の成果を調べた。手元にはLASのほかに、 org.txt 、 grd.txt 、 0.5g.txt 、 lem など、似た名前のファイルが並んでいる。今回はその違いを確認し、 0.5m格子のTXTから三角形のTIN用テキストを作る Pythonツールを用意した。 元データには何が入っていたか 大まかな流れは、レーザで取得した点群から地物などを取り除いて地表面の点を作り、その点をもとに一定間隔の格子標高を作る、というもの。 DEMは地面の高さを数値で表したデータの呼び名 で、ここでは0.5m格子のTXTやLEMがそれに当たる。ファイルごとの違いを、データの並びも含めて整理した。 ファイル例 どんなデータか 記録の並び・形式 主な使い道 09ME6444.las 航空レーザの点群。地面のほか、樹木や建物などに当たった点も含み得る。 LASという バイナリ形式 。点ごとの座標X・Y・Z、反射強度、リターン番号などを記録する。テキストの1行ずつではない。 点群の保管・交換、分類、地表面点の抽出。 09me6444_org.txt オリジナルのレーザ計測点。地表以外も含む不規則な点の集まり。 1点1行のカンマ区切り: id,x,y,z,p 。末尾の p はパルス識別番号。 フィルタリングして地表面点を作る元データ。樹木や建物の高さの解析にも使う。 09me6444_grd.txt グラウンドデータ。地表面として抽出された 不規則な点群 。 grd は格子の意味ではない。 1点1行のカンマ区切り: id,x,y,z 。各点の位置と標高を記録。 地面の点を直接調べる、格子標高やTINを作る材料にする。 09me6444_0.5g.txt 0.5m間隔に整えた 格子標高(DEM) 。 1格子点1行: id,x,y,z,A 。例: 1,19000.25,-80250.25,162.00,0 。 A は地表面点の有無などを示す属性値。 地形表示・断面・傾斜の計算。 今回のTIN作成元 。 09me6444_0.5g.csv LEMに対応する格子データのヘッダー。 格子の行列数、間隔、位置、測量年、座標系番号などの 項目と値 を記録。標高を1点ずつ並べたTXTとは役割が違う。 LEMの標...