人工衛星画像の入力

 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
 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のデータでは,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単位で入力します.この値は,異なる場合があります.)

 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コマンドで,ヘッダー部分の変更を行ない領域を正しく変えてください).最後に,正しく入力されているか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:
1______
 end row: 5965___

 start col: 1______
 end col: 6920___

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

なお,ここで示した値は,場合により異なることがありますので注意してください.


戻る