Home

Welcome to the SGLI_L1 wiki!

GCOM-C（しきさい）SGLIのL2 hdf5のファイルをgeotiffに変換し、ESPG4326に地図投影するプログラムです。
python3系とgdalがインストールされている環境であれば動くと思います。

使い方

python L1.py sgli_l1.h5 output_color.tif

ほかのバンドのを抜き出したい場合、
#RGBの3バンドだけ抽出する
sgli_bands=['Lt_VN03','Lt_VN05','Lt_VN08']

で別のバンドを指定します。
Lt_VN01~Lt_VN11まで選べます。

位置の検証

Open Street Mapを使って検証します。
L1のhdf5のファイルには位置情報も入っており、それを使えば海岸線の位置ずれはあまりないと思います。
プログラムではgcpを200ピクセルおきにとっていますが、100ピクセルおきに入れることもできます。
おそらく、gdalwarpコマンドを叩くより多くのGCPを扱える。ただし、処理に時間がかかるようになります。

べた塗のものがopen street mapのマダガスカル島です。