衛星画像の輝度温度

Tokio Kikuchi

衛星が観測するのは

輝度温度
これを表示するには
データである画像と
較正値 (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 にしておくこと)
緯度経度決定のために、参考画像