火山DEM10mメッシュのインポート

●HEADERの確認

●インポート

A）LEMファイルをインポート

LEMファイル LEMファイルはテキスト形式のデータ（拡張子はlem）です。このファイルでは、DEMデータが次のように並べられています。 レコード番号 標高値 標高値 ・・・ 標高値 レコード番号 標高値 標高値 ・・・ 標高値 レコード番号 標高値 標高値 ・・・ 標高値 ・ ・ ・ ・・・ ・ ・ ・ ・ ・・・ ・ ・ ・ ・ ・・・ ・ レコード番号 標高値 標高値 ・・・ 標高値 このデータは、海部及び空白部を「−９９９９」で記述してします。また、標高値は１０ｃｍ単位で表現しています。 ▼「レコード番号」は必要ないので、ExcelでDemデータ（.LEM）を開き、「レコード番号」を削除します。A列目が「レコード番号」となりますので、A列目を削除します。 削除しましたら、テキスト（タブ区切り）として保存します。 ▼メインメニューから、[メイン]＞[入力（インポート）．．.]を選択します。 ▼インポートするデータを選択します。＜インポート（入力）＞ウィンドウの[入力形式]の欄にある[フィルタ]で“ラスタ”を選択し、 さらに“TEXT-RASTER”を選択し、[ファイルの選択…]を押します。 ▼＜ファイル（複数）を選択＞ウィンドウが現れます。Excelで編集したテキストファイルを選択します。 ▼選択しましたら、次へ進みます。 ▼＜インポートパラメータ＞ウィンドウが現れますが、まずはインポートするためのパラメータをHEADERフォルダにある“HEADERファイル”から 読み取ります。ここでは「岩手山」をインポートするので、「岩手山」のHEADERファイルを表示しています。 HEADERファイルの“南北方向の点数”がライン数で“東西方向の点数”がカラム数になります。 パラメータは次のようになります。 ライン（Y方向)数およびカラム（X方向)数 ：６００×８００ ライン方向およびカラム方向のセルサイズ ：１０m ピクセルのビット数 ：１６ビット符号付き ライン数およびカラム数はデータによって異なりますのでご注意ください。 設定しましたらインポートを始めます。 ▼インポートしたデータを表示するとこのようになります。見やすいように、陰影表示しています。標高値は１０cm単位で表示しています。

B）DATファイルをインポート

DATファイル DATファイルはバイナリ形式で保存されたデータ（拡張子はdat）です。このdatファイルは“ビッグエンディアン”（２バイトデータは １バイトずつ分けて記録されますが、その際上位のバイトから記録されています）で記録されていますので、この点に注意してインポートします。 LEMファイルと同じように、海部及び空白部を「−９９９９」で記述してします。また、標高値は１０ｃｍ単位で表現しています。 ▼メインメニューから、[メイン]＞[入力（インポート）．．.]を選択します。 ▼インポートするデータを選択します。＜インポート（入力）＞ウィンドウの[入力形式]の欄にある[フィルタ]で“ラスタ”を選択し、 続けて“SIMPLE-ARRAY”を選択し、[ファイルの選択…]を押します。 ▼＜ファイル（複数）を選択＞ウィンドウが現れます。「岩手山」をインポートしますので、DATフォルダにある“IWATE.DAT”を選択します。 ▼ファイルを選択しましたら次へ進みます。 ▼＜インポートパラメータ＞ウィンドウが現れます。インポートするためのパラメータをREADMEフォルダにある“READMEファイル”から 読み取ります。「岩手山」をインポートするために必要となるパラメータは、以下のとおりです。 ライン（Y方向)数およびカラム（X方向)数 ：６００×８００ ライン方向およびカラム方向のセルサイズ ：１０ｍ ピクセルのビット数 ：１６ビット符号付き バイトオーダー ：高−低 なお、ライン数およびカラム数はデータによって異なるので、インポートするデータに置き換えてください。 以上、設定しましたらインポートを始めます。 ▼インポートしたデータを表示するとこのようになります。見やすいように、陰影表示しています。標高値は１０cm単位で表示しています。

●ジオリファレンス

▼HEADERファイルに書かれている座標情報を確認します。ここに書かれている座標情報は新世界測地系（JGD２０００）に対応しています。 ▼メインメニューから[メイン]＞[ジオリファレンス（座標付与）...]を選択します。 ▼＜オブジェクトのジオリファレンス処理＞ウィンドウの[ファイル]＞[開く...]から、先ほどインポートしたDemデータを開きます。 ▼＜ジオリファレンス オプション＞ウィンドウで座標参照系を指定します。「岩手山」の座標参照系は“JGD２０００/Japan-１９ �]”です。 ▼HEADERファイルに書かれているラスタの左下の座標と右上の座標を参考に、４隅の座標を入力します。 注意！ ・HEADERファイルに書かれている（X,Y)座標はcm単位で書かれていますが、この作業ではm単位で座標値を入力します。 ・この座標系は平面直角座標系です。そのため、「X座標が南北方向」を示し、「Y座標が東西方向」を示します。 参考までに、「岩手山」の4隅の座標は次のようになります。座標を入力しましたら、保存します。 位置 ライン カラム Northing Easting 左下 600 0 -19692.43 10950.78 左上 0 0 -13692.46 10950.78 右下 600 800 -19692.43 18950.87 右上 0 800 -13692.46 18950.87

●標高値の調整

▼メインメニューから[ラスタ]＞[組み合わせ（演算）]＞[一般演算式．．．]を選択します。 ▼＜ラスタの組み合わせ演算＞ウィンドウでの[ラスタ（複数）．．．]から、インポートしたDEMデータを選択します。 選択しましたら、パラメータを設定します。 ラスタに対し、どのような演算を行うか設定します。 ラスタ全体に同じ数をかける場合、タイプで“代数”を選択し、操作で“スケール／オフセット”を選択します。 続けてパラメータを設定します。スケールファクタにラスタ全体へかける数を入力します。今回は“０.１”です。 ラスタ全体に０.１をかけるとラスタ値は小数点をとるので、ラスタタイプは“３２ビット浮動小数点”を選択します。 まとめると、次のようになります。 タイプ ：代数 操作 ：スケール／オフセット スケールファクター ：０.1 ラスタタイプ ：３２ビット浮動小数点 設定しましたら、実行します。 ５０ｍDEMデータと重ねて表示すると、このようになります。 ５０ｍDEMデータと１０ｍDEMデータを比較するため、同じ場所を拡大表示しました。上のウィンドウが１０ｍDEMデータで、下のウィンドウが５０ｍDEMデータです。１０ｍDEMデータのほうがより細かく地形を表現していることがよくわかります。