slc.par ファイルから国土基盤情報(標高10m)データを使って
必要範囲の DEM データ (GAMMA 形式) を作成するスクリプト
(GAMMA)
シェルスクリプト名
- get_DEM_par_UTF8.sh
- ※本スクリプトは、GAMMA ソフトウェアでのみ動作します。
使用前の準備
- パラメータファイル:「Mosaic_FG-GML-DEM10B-20200409_extents_vec-20201019_a-ATT_dd_UTF8.csv」
- 基盤地図情報 10mメッシュ(標高)の各タイル (GeoTIFF) の図郭名、TIFF ファイル名、四隅の緯度経度情報を含むリストファイル。
- 処理にあたってスクリプトから参照します(必須)。
- 通常は、スクリプトと同じ場所(カレントディレクトリ)にあれば問題はありません。
- 特別なディレクトリ (TOOLSなど) に置く場合は、スクリプト内のパスを修正してください。
################################################################# # DEMファイル名(GeoTIFF)と範囲のリストファイル(.csv)のパス: kiban10m_param="./Mosaic_FG-GML-DEM10B-20200409_extents_vec-20201019_a-ATT_dd_UTF8.csv" #################################################################
- 使用する環境に応じて、スクリプト内の、国土基盤情報(標高10m)データフォルダー(ディレクトリ)パスを修正してください。
################################################################# # 国土基盤情報(標高10m)データのパス: dir_DEM_TIFF="/G/FG-GML-DEM10B-20200409_GeoTIFF" #################################################################
実行手順
$ sh ./get_DEM_TIFFlist_UTF8.sh
と入力、Enter キーを押します。
($ はシェルスクリプトのプロンプト)- 上記パラメータファイルが正しく設定されていないと、以下のメッセージが表示されます。
####################################################### DEMと範囲の対応表(.csv)が見つかりません! 「Mosaic_FG-GML-DEM10B-20200409_extents_vec-20201019_a-ATT_dd_UTF8.csv」 パスを正しく設定してください! #######################################################
- 国土基盤情報(標高10m)データのパスが正しく設定されていないと、以下のメッセージが表示されます。
####################################################### 国土基盤情報(標高10m)のディレクトリが見つかりません! データ(GeoTIFF)のあるパスを正しく設定してください! #######################################################
- 上記エラーが出た場合は、スクリプト内の変数を自分の環境に合わせて修正します。
- 上記がすべてクリアされると、以下のような使用法の案内 (Usage) が表示されます。

- 実際に .par ファイルを引数に置いて実行すると、処理が進みます。

- 処理が終わると、作成された DEM ファイルが表示されます。
作業ファイル
カレントディレクトリには以下の 5個のファイルが作られます。
- slc_corners_out.txt
- target_meshcode.lst
- DEM.vrt
- DEM.EQA.dem(成果品)
- DEM.EQA.dem_par(成果品)
最後の 2つが、使用した slc.par ファイルの範囲で取り出した GAMMA-SAR 形式の DEM データです。
座標系は、緯度経度座標 (EQA) です。
作成
株式会社オープン・ジー・アイ・エス