GCOM-Cデータのダウンロード
GCOM-Cとは
JAXAが2017年12月23日に打ち上げた気候変動観測衛星「しきさい」です。正式名称を地球環境変動観測ミッション:Global Change Observation Mission Climate (GCOM-C) といいます。可視から熱赤外光を観測する多波長光学放射計(SGLI:Second generation GLobal Imager)を搭載し,雲,エアロゾル,海色,植生,雪氷等を全地球規模で長期間継続的に観測することで,将来の気温上昇量の正確な予測に必要となる放射収支,および炭素循環の変動メカニズムの解明に貢献することを目的としています。
GCOM-Cデータのダウンロード
ftpでデータサーバにアクセスする
まずは,アカウントをG-Portal (https://gportal.jaxa.jp/gpr/) で作成する。このHPからダウンロードすることが一番簡単。
ここでは,ftpによるダウンロード方法について (自分のためのメモ程度のざっくりとした) 説明をします。
$ sftp -oPort=2051 UserName@ftp.gportal.jaxa.jp
password: **********
ftp接続できたら,ダウンロードしたいファイルを探します。いくつかのファイルだけを入手したいときは,そのままgetしましょう。複数ファイルが入手したいときは,ひとつひとつダウンロードしていくと面倒くさいので,ここでディレクトリやファイル名のルールを確認して,wgetしちゃおう。
Example:
sftp> ls standard/GCOM-C/GCOM-C.SGLI/L2.LAND.LAI_/1/2018/11/10/
sftp> get standard/GCOM-C/GCOM-C.SGLI/L2.LAND.LAI_/1/2018/11/10/GC1SG1_20181110D01D_T1413_L2SG_LAI_Q_1000.h5
wgetで一気にデータをダウンロードする
# 1 month composites of the GCOM-C RGB band data
$ cd ${DATA_DIR}/L3
$ for band in RV03 RV05 RV08; do
# Ascending (A) とdescending (D) のそれぞれのファイルがあるので,ダウンロードするときには,気をつけること!
# (Descending (D) を使おう! Ascending (A) だと,北半球の北のほうのデータしかない。これは私が可視光の反射率を使っているためなので,ご自分の状況に応じて判断してください。)
# GC1SG1_20190801[A or D]01M_D0000_3MSG_${band}F_1001.h5
$ wget -c --user=UserName --password=anonymous ftp://ftp.gportal.jaxa.jp/standard/GCOM-C/GCOM-C.SGLI/L3.LAND.${band}/1/2019/08/GC1SG1_20190801D01M_D0000_3MSG_${band}F_1001.h5 -P ${DATA_DIR}
$ done
# 8-day statistics of the GCOM-C RGB band data
$ cd ${DATA_DIR}/L3
$ for band in RV03 RV05 RV08; do
$ wget -c --user=UserName --password=anonymous ftp://ftp.gportal.jaxa.jp/standard/GCOM-C/GCOM-C.SGLI/L3.LAND.${band}/1/2019/08/GC1SG1_20190805D08D_D0000_3MSG_${band}F_1001.h5 -P ${DATA_DIR}
$ done
# LAI data
$ cd ${DATA_DIR}
$ for i in 08 09 10; do # month
$ for TILE in 0428 0529; do
$ wget -c --user=UserName --password=anonymous ftp://ftp.gportal.jaxa.jp/standard/GCOM-C/GCOM-C.SGLI/L2.LAND.LAI_.Statistics/1/2019/${i}/GC1SG1_2019????D08D_T${TILE}_L2SG_LAI_Q_1000.h5 -P ${DATA_DIR}
$ done
$ done
よく使うGCOM-Cのタイル (Tile ID)
北海道: 0428,本州+四国: 0529,中国地方: 0528,沖縄: 0629, モンゴル: 0425, タイ: 0727
H5ファイルの中身を確認する
$ h5dump -n GC1SG1_20190801A01M_D0000_3MSG_RV03F_1001.h5
h5dump -n [h5ファイル名]
#### ファイルの中身が表示される
#HDF5 "GC1SG1_20190801A01M_D0000_3MSG_RV03F_1001.h5" {
#FILE_CONTENTS {
# group /
# group /Geometry_data
# group /Global_attributes
# group /Image_data
# dataset /Image_data/Rs_VN03_AVE
# dataset /Image_data/Rs_VN03_QA_flag
# group /Processing_attributes
# }
#}
使うdataset (/Image_data/Rs_VN03_AVE) のパラメーターを確認する
$ h5dump -A --dataset=/Image_data/Rs_VN03_AVE GC1SG1_20190801A01M_D0000_3MSG_RV03F_1001.h5
h5dump -A --dataset=[dataset名] [h5ファイル名]
#### datasetのパラメーターが表示される
#### "Offset","Slope","Unit","Error_DN"等の重要な情報をここで入手する。
GCOM-Cデータ(h5ファイル)を使うためのpythonツールは下のURLからダウンロードできます。
The python tool to use the h5 file of SGLI can be downloaded from the followings:
https://github.com/K0gata/SGLI_Python_Open_Tool
GCOM-CのLevel2のHDF5ファイルをgeotiffに変換するプログラムはG-portal (https://gportal.jaxa.jp/gpr/information/tool?lang=ja#GCOM-C) からダウンロードできます。
The h5 file can be converted into a geotiff file using the converter. It can be downloaded from G-portal on the JAXA homepage: https://gportal.jaxa.jp/gpr/information/tool?lang=en#GCOM-C.
使い方の例:
./SGLI_geo_map_linux.exe GC1SG1_20240328D01D_T0529_L2SG_RSRFQ_3002.h5 -r NN -d /Image_data/Rs_VN08
※-r はリサンプリングのオプション 「nearest neighbor」、-d は、geotifファイルにしたい「HDFファイルのデータセット」の名前。データセット名は、上で書いた「h5dump -n」コマンドや、「gdalinfo」コマンドで確認。gdalinfoの場合は、SubdatasetsのNameを確認。
※-rオプションは、-dオプションよりも前に書く。