人工衛星画像の入力

　GRASSには，LANDSATのMSSやTM，およびSPOTなどのデータをインポート（変換入力）するコマンドが用意されています．これらは，i.tape.で始まるコマンドです．その他，あらかじめ設定されていない衛星画像のフォーマットに対応するために，i.tape.otherというコマンドがあります．

　ここでは，i.tape.otherを用いて日本の人工衛星であるJERS-1のOPS（光学センサー）データをCD-ROMからGRASSへインポートする方法の例を示します（CD-ROMはRESTEC等から購入できます）．コマンド名の中にtapeという単語があるように，本当は磁気テープ装置から読み取るようになっていますが，CD-ROMを含む普通のファイルでも読めます．通常CCTテープは先頭から順に読み込むようになっていますが，CD-ROMや通常のファイルでは，順番に読むかわり，ファイル名がついています．このため，ファイルのスキップなど通常の操作とは若干感じが違うので注意してください．なお，画像データの取り込みは，xy座標系で行ないます．地図などと合わせて利用する場合には，その後にUTM座標系に位置合わせをするなどして変換し，使用します．

　grass 4.2.1 > i.tape.other　↓

　Enter tape device name: /cd-rom/ops/opsimage_b2_p66_r241

　　　（ファイル名をフルパスで入力します．）

　Please mount and load tape, then hit RETURN-->　↓

---------------------改ページ------------------

　　　（必要に応じて，記載してください．必要なければESCキーを押します．）

　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

　

　0___ number of tape files to be skipped

　　　（CD-ROM等のファイルの場合，スキップする必要はありません．）

　1___ number of records in the remaining files to be skipped

　　　（JERS-1のデータでは，１レコード目がヘッダーなので，スキップします．）

　32__ number of bytes in each record to be skipped

　　　（各レコードの最初の部分は，画像ではないのでスキップします．

　　　　そのbyte数を指定します）

　band files

　

　1___ number of bands on the tape

　　　（ファイルにはいくつのバンドが入っているかを指定します．）

　

　data format

　_ band sequential (BSQ) | mark one with an x

　x band interleaved (BIL) |

　　　（データのフォーマット形式を指定します．xを入力します．）

　

　0_____ if you select BSQ format and all the bands are in a single file,

　enter the total number of records in the file. Otherwise enter 0

　　　（必要に応じてデータの総レコード数を入力します．ここでは0を入力しています．）

　4680__ length (in bytes) of the longest record on the tape

　　　（レコードの長さをbyte単位で入力します．ここでは4680を入力しています．）

　1___ blocking factor of data in the file

　AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE

　(OR <Ctrl-C> TO CANCEL)

　ESCキーを押す．

---------------------改ページ------------------

　Please mark an x by the bands you want extracted

　

　x 1

　　　　（1種類のデータしかないので，1にxを入力します．）

　AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE

　(OR <Ctrl-C> TO CANCEL)

　ESCキーを押す．

---------------------改ページ------------------

　select a prefix/group for the band cell files to be created

　Enter 'list' for a list of existing imagery groups

　Enter 'list -f' for a verbose listing

　Hit RETURN to cancel request

　> band2

　　　　（GRASSで利用するラスターファイル名を入力します．）

---------------------改ページ------------------

　EXTRACT

　

　please select region of the image to extract

　

　

　start row: 1______

　end row: 4200___

　

　start col: 1______

　end col: 4500___

　　　　　（取り出す画像の列(row)行と列(col)の最初と最後の値を指定します．

　　　　　　ここの例では，JERS-1の全データを取りこんでいます．）

　

　AFTER COMPLETING ALL ANSWERS, HIT <ESC> TO CONTINUE

　(OR <Ctrl-C> TO CANCEL)

　ESCキーを押す．

---------------------改ページ------------------

以上です．この方法で取り出すと，何故か画像の領域が負に変わってしまいます．r.supportコマンドで，ヘッダー部分の変更を行ない領域を正しく変えてください．

　

第11回実習へ戻る