人工衛星画像の入力
GRASSには,LANDSATのMSSやTM,およびSPOTなどのデータを入力するコマンドが用意されています.これらは,
i.tape.で始まるコマンドです.その他,あらかじめ設定されていない衛星画像のフォーマットに対応するために,i.tape.otherというコマンドがあります.ここでは,
i.tape.otherを用いて日本の人工衛星であるJERS-1のOPS(光学センサー)データ,およびLANDSATのTMデータをCD-ROMからGRASSへ入力する方法の例を示します(CD-ROMはRESTEC等から購入できます).コマンド名の中にtapeという単語があるように,本当は磁気テープ装置から読み取るようになっていますが,CD-ROMを含む普通のファイルでも読めます.通常CCTテープは先頭から順に読み込むようになっていますが,CD-ROMや通常のファイルでは,順番に読むかわりにファイル名がついています.このため,ファイルのスキップなど通常の操作とは若干感じが違うので注意してください.なお,画像データの取り込みは,xy座標系で行ないます.地図などと合わせて利用する場合には,その後にUTM座標系に位置合わせをするなどして変換し使用します.(1)i.tape.otherコマンドを用いたJERS-1画像入力
GRASS
を起動します.つぎにi.tape.otherコマンドを実行します.GRASS 5.0beta6 > i.tape.other
(ファイル名をフルパスで入力します.)
Please mount and load tape, then hit RETURN-->
Please enter the following information
TAPE IDENTIFICATION
_______________________________________________________________________
IMAGE DESCRIPTION
_______________________________________________________________________
_______________________________________________________________________
_______________________________________________________________________
TITLE FOR THE EXTRACTED CELL FILES
______________________________________________________________________
AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE
(OR <Ctrl-C> TO CANCEL)
ESC
キーを押す. GENERIC TAPE EXTRACTION
tape layout
(CD-ROM等のファイルの場合,スキップする必要はありません.)
1___ number of records in the remaining files to be skipped
(JERS-1のデータでは,1レコード目がヘッダーなので,スキップします.)
32__ number of bytes in each record to be skipped
(各レコードの最初の部分は,画像ではないのでスキップします.ここには,そのbyte数を指定します)
band files
(ファイルにはいくつのバンドが入っているかを指定します.)
(必要に応じてデータの総レコード数を入力します.ここでは
0を入力しています.)4680__ length (in bytes) of the longest record on the tape
(レコードの長さをbyte単位で入力します.この値は,異なる場合があります.)
1___ blocking factor of data in the file
ESC
キーを押す.(1種類のデータしかないので,1にxを入力します.)
AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE
(OR <Ctrl-C> TO CANCEL)
ESC
キーを押す.(GRASSで利用するラスターファイル名を入力します.)
(取り出す画像の列(row)行と列(col)の最初と最後の値を指定します.ここの例では,JERS-1の全データを取りこんでいます.)
AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE
(OR <Ctrl-C> TO CANCEL)
ESC
キーを押す.以上で入力作業は終了です(現時点でのプログラムのバグではないかと思いますが,この方法で取り出すと,なぜか画像の領域が負に変わってしまいます.r.supportコマンドで,ヘッダー部分の変更を行ない領域を正しく変えてください).最後に,正しく入力されているか
d.rastコマンドを使って確認しましょう(2)i.tape.otherコマンドを用いたTM画像入力
基本的に,JERS-1と同じですが,パラメータが下記のように異なります.以下には,JERSと異なる場所のみ示し,青色で数字等が示してあります.
★ファイル名が異なります.
Enter tape device name:
/cdrom/scene001/imgy_01.dat(ファイル名が異なります.)
★テープのレイアウトが一部異なります.
GENERIC TAPE EXTRACTION
tape layout
7020__ length (in bytes) of the longest record on the tape
(レコードの長さをbyte単位で入力します.この値は,異なる場合があります.)
★画像の大きさが異なります.
please select region of the image to extract
start row:
(取り出す画像の列(row)行と列(col)の最初と最後の値を指定します.)
なお,ここで示した値は,場合により異なることがありますので注意してください.