AVNIR2のCEOSフォーマット(BSQ方式)を読み込むプログラムをCとかで書く。gdalを使って、なんとか頑張る。北極の画像を使ったときの工程)。バイナリBSQ方式)に、ゲインaとオフセットbが格納されている。http://www.eorc.jaxa.jp/ALOS/doc/jformat.htm放射輝度は、ゲインとオフセットを使って、以下のように計算される。太陽天頂角AVNIR-2のF0(太陽の放射照度スペクトルを、AVNIR2のバンド応答関数で重み付き積分したもの。1076.5太陽天頂角は、summary.txtに書いてある。太陽地球間距離dは、観測時刻から計算する。天文単位での太陽地球間距離を求めるには、平均値D0(=1[AU])を観測日の太陽地球間距離Dで割ればよい。以下のような近似式がある(たぶん、時系列のdをフーリエ級数展開したもの)。Aを365日で1周期としているけど、閏年(DOYが366のとき)はどうするんだろう?(D0/D)^2=1.00011+0.034221*cosA+0.00128*sinA+0.000719*cos2A+0.000077*sin2AA=(2π/365)*DOY大気補正(地表面反射率への変換)検討中。AVNIR2の応答関数が必要。JAXA村上さんに聞くのが一番早いかも。→聞いてみたら、一般公開されてないので、田殿さんに聞いてみたらとのこと。その他UTMの画像なのに、画像の4隅の座標が緯度経度でしか書かれていない。4隅の緯度経度が小数点以下第3位までしか示されていない。