衛星画像の輝度温度

衛星が観測するのは

• 輝度温度
• これを表示するには
• データである画像と
• 較正値 (Calibration data) が必要

較正表

• 画像＝データ
• 0 〜 255 の数値 (オリジナルは 0 〜 1023)
• 保存書庫 -> 換算表
• ダウンロードした CAL データを見てみる
• gunzip MT1R......CAL.dat.gz -> gz が取れる (cat/more で読める)
• ./calview MT1R......CAL.dat -> グラフ表示

画像をデータとして扱う

• 画像の元になっているのは、赤外線で見たときの温度
• 画素の値 -> 温度の換算表 -> CAL
• GAME データを取ってきて利用
• GAME と CAL
• 同じ日、同じ時刻の IR1 と CAL を保存 (作業ディレクトリに) 例:

     http://weather.is.kochi-u.ac.jp/sat/CAL/2010/Jul/MT1R10071203CAL.dat.gz
     http://weather.is.kochi-u.ac.jp/sat/GAME/2010/Jul/IR1/MT1R10071203IR1.pgm.gz
  

解凍

• gunzip で解凍:

     gunzip MT1R10071203*.gz
  

グラフ表示プログラム

• gameview を準備してあるので、:

     ./gameview MT1R10071203IR1.pgm
  

• で、表示 (10度K おきの等高線で表示されるはず)

グラフ表示／上級編

• GAME データは 北緯70度から南緯20度、東経70度から東経160度の広範囲なので、切り出して拡大表示したい。
• --top --bottom --left --right で範囲を指定できる。
• 例 北緯40度から20度 東経120度から140度の場合
• ./gameview --top 40 --bottom 20 --left 120 --right 140 MT1R10071203IR1.pgm

さらに

• 等高線（等温線）の間隔を 5K おきにしたいときは
• --interval 5 で指定する
• ./gameview --top 40 --bottom 20 --left 120 --right 140 --interval 5 MT1R10071203IR1.pgm
• IR2, IR3, IR4 の表示も可能 VIS は0-255の数値をそのまま使用
• 出力画像ファイルを --outfile example.png で指定すると画像ファイルで保存 (拡張子は .png にしておくこと)
• 緯度経度決定のために、参考画像