1.数値地図50mメッシュ(標高)の入力
国土地理院刊行の数値地図50mメッシュのような緯度経度分割によるデータをGRASSのUTM座標系へ入力する場合はr.in.llを用います.r.in.llはバイナリー形式のラスター型データを入力するコマンドであるため,それに対応したデータを作成しなければなりません.それらの手順を以下に示します.
(1)変換プログラムを用いてGRASSのバイナリー形式に変換
つぎのプログラムをダウンロードし,解凍します.プログラムは2つあり,アスキー形式のラスターデータに変換するdmj2gras.fとアスキー形式のラスターデータをGRASSのバイナリ形式に変換するasc2gbin.fです.
プログラムのダウンロード(ここをクリックしてください;1.4KByte)
tar -zxvf giscvprog.tar.gz ファイルの解凍
両プログラムとも,フォートラン77の形式で書かれていますので,コンパイルが必要です.コンパイルは,それぞれつぎのように行ないます(文頭のf77はシステムにより異なり,g77の場合もあります).
f77 dmj2gras.f -o dmj2gras.out ↓
f77 asc2gbin.f -o asc2gbin.out ↓
(f77 プログラムソース.f -o 実行形.outです.
赤色部分を入力する,↓はリターン)コンパイル後,dmj2gras.out,asc2gbin.outの順に実行します.なお,GRASSコマンドr.in.llに必要なパラメータはdmj2gras.outの実行時に表示されます.各実行例を以下に示します.まず,CD-ROMからアスキー形式でデータを抽出します.
dmj2gras.out
↓ (プログラムを起動します)Input file name =
/cdrom/data/5235/523500.mem ↓(変換するcdromのファイル名を入力します)
Output file name =
dem.asc ↓(出力するアスキー形式のファイル名を入力します)
sea area code =
0 ↓ (海コード(-9999)の標高を指定します)convert to meter (divide by 10)? [y]/n =
↓(もとの標高データは10倍されているので,10で割るか
どうかの指示を与えます.デフォルトはyです.)
(r.in.llに必要な情報(以下の
緑色の文字)が次のように表示されます)number of byte per cell =
4south-west corner of lati. & long. =
sw,34:40:00.75N,135:15:01.125Enumber of rows and columns of data =
200,200resolution of data latres,lonres =
1.5,2.25spheroid =
besselつぎに,バイナリー形式にデータを変更します.
asc2gbin.out
↓File name for input (ascii raster format) =
dem.asc ↓File name for output(Grass binary format) =
dem.bin ↓ここで作成したバイナリー形式のファイル(例:dem.bin)をGRASSコマンドr.in.llを用いてUTM座標形に変換・入力します.
(2)r.in.llコマンドを用いてGRASSに入力
r.in.llコマンドを実行します.
grass4.2.1>
r.in.ll ↓OPTION: input file
key: input
required: YES
enter option >
dem.bin ↓(バイナリー形式のラスター型データファイル名を入力します.)
----- 省 略(以降,入力の確認は省略します) -----
OPTION: output raster file
key: output
required: YES
Enter a new raster file name
Enter 'list' for a list of existing raster files
Hit RETURN to cancel request
>
elevation ↓(GRASSで利用するラスターファイル名を入力します.)
OPTION: number of bytes per cell
key: bpc
required: YES
enter option >
4 ↓(1つの格子(セル)データが用いるbyte数を入力します.
前述のプログラムを用いた場合は4byteです.)
OPTION: one corner latitude and longitude of the input
format: {nw|ne|sw|se},dd:mm:ss{N|S},ddd:mm:ss{E|W}
key: corner
format: corner,lat,lon
required: YES
enter option >
sw,34:40:00.75N,135:15:01.125E ↓(データ領域の4隅のいずれかの緯度経度の値を書式に従って入力します.
例では,sw(南西)の値が示されています.なお,この時の値は隅の格
子セルの中心の座標値を意味します(r.in.asciiとは異なります).
従って,swの場合にはデータ領域の南西の位置座標に1格子の大きさの
半分を足した値を指定します.)
OPTION: number of rows and columns in the input file
key: dimension
format: rows,cols
required: YES
enter option >
200,200 ↓(データの大きさ(格子数)の行・列数を入力します.)
OPTION: resolution of the input (in arc seconds)
key: res
format: latres,lonres
required: YES
enter option >
1.5,2.25 ↓(データの緯度,経度の分解能(格子の大きさ)を秒単位で入力します.
国土地理院の50mメッシュでは,それぞれ1.5,2.25です.)
OPTION: spheroid
key: spheroid
required: YES
enter option >
bessel ↓(変換時に参照する楕円体名を入力します.besselを用います.)
FLAG: Set the following flag?
Signed data (high bit means negative value)?(y/n) [n]
y ↓(データに負の値がある場合はyを入力します.)
以上でラスター型データファイルおよびデータのサポートファイルが作成され,変換および登録が終了します.