SGLI

5. SGLI daily images

※本ページのデータ作成に使用した原初データは、宇宙航空研究開発機構より提供を受けました。

A) Daily global NDVI images

解像度は約4km（0.0333333333 degrees=1.6′）に落として作成しています。

2025年10月 2025年10月1日

B) Daily RGB images around Japan

解像度は約250m（0.00208333333333 degrees=7.5′′）で作成しています。

2025年12月 2025年12月1日

C) Monthly composited RGB images around Japan

解像度は約250m（0.00208333333333 degrees=7.5′′）で作成しています。
ex. 2023年1月20日の画像=2022年12月21日〜2023年1月20日のコンポジット画像

2022年1月〜12月 2022年8月31日

D) Monthly composited NDVI images

解像度は約4km（0.0333333333 degrees=1.6′）に落として作成しています。
ex. 2023年1月20日の画像=2022年12月21日〜2023年1月20日のコンポジット画像

2022年1月〜12月 2022年8月31日

6. GCOM-C/SGLIの基本情報

・ 2017年12月23日にH-ⅡAロケット37号機で打ち上げられた。
・ GCOM-Cは衛星の名前。
・ SGLIはセンサの名前。
・回帰日数は34日。極域は1日に何度も観測されるが、赤道付近は2.5日に1回くらい観測される。
・ VNR-NPは3鏡筒（11 channel）、VNR-PLは2鏡筒（2 channel）、IRSは走査式（6 channel）。
・解像度はだいたい250m。海域はだいたい1km。偏光鏡筒（VNR-PL）はだいたい1km。
・偏光鏡筒のチルト角は45°。つまり、地球の曲率を考えると衛星天頂角は52.0-60.0°くらいになる。
・偏光鏡筒は北半球では前方視、南半球では後方視。黄道付近でチルトが変わる。
・ VNR-NPとVNR-PLの観測時間のずれは大体2分くらい。
・ GCOM-C軌道（SGLIがどこを観測した or 観測するか）を知りたい場合は、https://suzaku.eorc.jaxa.jp/GCOM_C/orbit/orbit_map_j.html https://shikisai.jaxa.jp/obsplan/obsplan_j.html
・明日観測する場所と衛星幾何を簡単に知りたい場合は、34日前の画像を見れば分かる（ただし軌道修正する時があるので注意）。
・標準プロダクト（L2SG）が作成されるのは、観測から1日後くらい。
・ GCOM-Cの処理状況や制約事項を知りたい場合は、https://shikisai.jaxa.jp/faq/faqlist_j.html
・特定の期間において、熱赤外バンド（T01）の特定画素にノイズが発生したことに伴い、RSRFQではTI01を使う雲識別に異常が発生しているので注意。砂漠とかでcloudやpol_cloudが立っている。
・ RSRF:台風一過で快晴の日はshadowフラグが立っていたりするので注意。
・ RSRF:POLのrecoveryを知るには >>> QA_flagのLTOAK画素相当（16画素）の一つでもPOL_recoveryが立っていればrecoveryの情報が含まれている。
・ RSRFのKチャンネル生成法:4x4に同じ値を格納→画素毎に雲識別、大気補正、BRDF推定を行う→4x4平均してKチャンネルを出力
・ POLのチルト駆動が停止するとLTOAが作成されない時があるらしい。
・ Australiaの湖（Eyre, Gairdner, Torrensなど）、ChileのUyuni塩湖でwater_flagが立っているので注意（元にしている水域マップで水域としているため。変更するのは難しいらしい。）。陸域標準プロダクトは物理量を計算しない。陸域が水域になっていると大気補正時に雲等と間違える事があるのでRSRFQを使うときは注意。→MODISプロダクト等では物理量を計算している。主なLCマップではpermanent water以外はbareやfloddedとしている。

※ SGLIの解像度（赤道直下。合っているかどうかは保証しない。）。
・ VNR(NP): FOV中央
[along-track] resolution = 249.4m, resampling interval ≒ resolution
[cross-track] resolution = 249.4m, resampling interval ≒ resolution
・ VNR(NP): FOV端
[along-track] resolution = 314.9m, resampling interval ≒ 250m
[cross-track] resolution = 412.0m, resampling interval ≒ resolution
※along-track方向とcross-track方向でresolutionが異なる。along-track方向ではresampling領域が隣接画素で重なる。
・ IRS: FOV中央
[along-track] resolution = 249.4m, resampling interval ≒ resolution
[cross-track] resolution = 249.4m, resampling interval ≒ resolution
・ IRS: FOV端
[along-track] resolution = 340.3m, resampling interval ≒ 250m
[cross-track] resolution = 492.7m, resampling interval ≒ resolution
※mechanical scanning なのでVNR(NP)より観測幅をやや広くとっている。
・ VNR(PL): FOV中央
[along-track] resolution = 2496.7m, resampling interval ≒ 1000m
[cross-track] resolution = 1512.6m, resampling interval ≒ resolution
※along-track方向ではresampling領域が隣接画素で重なる。
・ VNR(PL): FOV端
[along-track] resolution = 2982.5m, resampling interval ≒ 1000m
[cross-track] resolution = 2195.7m, resampling interval ≒ 1765m
※along-track方向とcross-track方向でresolutionが異なる。along-track方向ではresampling領域が隣接画素で重なる。

※ SGLIデータをダウンロードする方法。
・例えば、2023/04/01-05/31のVGI_Q（T0901〜T0917のタイル）をDLしたい時は、次のようなシェルを書いて実行すれば明日の朝までにはDLできているはず（G-Portalへのユーザー登録や鍵交換はユーザーマニュアルに従ってあらかじめやっておく）。 shファイル

※ Tile[VV,HH]内の画素[i,j]について緯度経度を計算する方法。
・ VV:南北方向のタイル番号（0始まり）、HH:東西方向のタイル番号（0始まり）
・ i:南北方向の画素番号（0始まり）、j:東西方向の画素番号（0始まり）
・ dims[0]:南北方向の画素数、dims[1]:東西方向の画素数（Qタイルの場合は4800、Kタイルの場合は1200）
・ ver_resolution:南北方向の解像度（ = 180.0/(18.0*dims[0]): Qタイルの場合は0.00208333333333、Kタイルの場合は0.00833333333333）
・ hor_resolution:東西方向の解像度（ = 360.0/(36.0*dims[1]): Qタイルの場合は0.00208333333333、Kタイルの場合は0.00833333333333）
int i;/* i */
int j;/* j */
int iv;/* VV */
int ih;/* HH */
float lat;/* 緯度 */
float lon;/* 経度 */
hsize_t dims[2];
float ver_resolution;
float r;
float M_PI=3.14159...;/* 円周率 */
lat = 90.0 - ( iv * dims[0] + i + 0.5 ) * ver_resolution;
r = cos( lat * M_PI / 180.0 );
lon = ( ( j + 0.5 ) * hor_resolution + ih * 10.0 - 180.0 ) / r;

※ GCOM-C標準プロダクトのデータやタイル抜け等の問題

・ 2020/04/30 --- 南米北西部のタイル(T0710,0711,0810,0910,1010)の陸地に水域フラグが立っている。なので陸域標準プロダクトはエラー値。
・ 2023/09/30-,2022/04/13-05/21,2022/01/18-03/05他 --- 全タイル。TIRバンドの異常によるRSRFの雲識別異常（stripe状）。なので陸域標準プロダクトはstripe状のエラー値。
・ 2022/03 --- GC1SG1_20220301D01M_Tvvhh_L2SG_RVxxQ/RTxxQ/RSxxQ...のデータ入力不足による欠損多発：全タイル。
・ 2023/04/16 --- GC1SG1_20230416D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトの抜け：T0516,T0615,T0616,T0617,T0715,T0716。
・ 2023/04/26 --- GC1SG1_20230426D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトの抜け：T0117,0118,0216,0217,0317,0318。
・ 2023/04/26 --- GC1SG1_20230426D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトのデータ欠損（すべてエラー値）：T0218。
・ 2023/04/27 --- GC1SG1_20230427D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトの抜け：T0717,0718,0817,0817。
・ 2023/04/29 --- GC1SG1_20230429D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトの抜け：T0911。
・ 2023/05/26 --- GC1SG1_20230526D01D_Tvvhh_L2SG_LTOAQ...と後続プロダクトの抜け：アジア東部のタイル抜け。
・ 2023/09/16 --- GC1SG1_20230916D01D_Tvvhh_L2SG_RSRFQ...と後続プロダクト（Rs_VN11が異常値なので後続も異常値）：T0508。07/28以降のタイルに時々発生。
・ 2023/10/30 --- GC1SG1_20231030D01D_Tvvhh_L2SG_RSRFQ...と後続プロダクトのデータ異常：T1011。
・ 2023/11/22 --- GC1SG1_20231122D01D_Tvvhh_L2SG_RSRFQ...と後続プロダクト（Rs_VN11が異常値なので後続も異常値）：T0619。

1. 今月のピックアップ画像

今月のSGLIアラート

A) 2025年8月：フランス オード県ナルボンヌ近郊の火災 [42.8°N〜43.4°N, 2.6°E〜3.2°E]

B) 2025年7月：トルコ北西部の火災 [40.0°N〜41.0°N, 29.5°E〜30.5°E]

C) 2025年8月：アメリカ カリフォルニア州の火災 [34.5°N〜35.5°N, 119.8°W〜120.8°W]

D) 2024年〜2025年：ブラジルの植生減少 [6.5°S〜7.5°S, 58.5°W〜59.5°W]

既存の Global Land Cover map （グローバル土地被覆図）からは分からない変化

北海道胆振東部地震[2018年9月6日]後の植生回復の様子

2026年1月3日、4日：関東地方の積雪

2025年12月：人吉盆地の霧（雲海）

2025年12月：亀岡盆地の霧（雲海）

2023年5月5日：霧ヶ峰高原の林野火災前後の様子

2023年4月11日、12日：Kamchatka半島の火山噴火前後の様子

2023年4月から5月にかけての東日本の展葉の様子([R:G:B]=[VN08:VN05:VN03])

2020年〜2023年の4月から5月にかけての東南アジアLST(LSTによる10-days MVC)

2023年2月27日：東日本

2022年4月：アメリカNewMexico州アルバカーキ近郊の森林火災

2. Global land Cover

1) タイル[Tvvhh=T0529(aroud Japan)]

A) Global Land Cover

B) Global Land Cover Change （強調表示しています）

2) 他のタイル

3. SGLIアラート

1) タイル[Tvvhh=T1012(aroud Bolivia and Brazil)]

B) NDVI年間アラート

A) NDVI月間アラート

2) 他のタイル

4. 画像集

1) 2020年5月20日： 北海道滝上町のシバザクラ

2) 2020年10月21日： 東北地方紅葉

3a) 2021年12月22日： 人吉盆地の霧その1

3b) 2022年01月10日： 人吉盆地の霧その2

4) 2022年1月7日： 関東地方積雪

5) 2022年3月： 阿蘇山周辺の野焼き

6) 2022年3月5日： 渡良瀬の野焼き

7) 2023年3月29日： アラビア半島の灌漑農業

8) 2023年5月3日： 西ヨーロッパ

★ 土地被覆の見え方： Around France

5. SGLI daily images

A) Daily global NDVI images

B) Daily RGB images around Japan

C) Monthly composited RGB images around Japan

D) Monthly composited NDVI images

6. GCOM-C/SGLIの基本情報

※ GCOM-C標準プロダクトのデータやタイル抜け等の問題

7. 便利なツール集

A) 2025年8月：フランスオード県ナルボンヌ近郊の火災 [42.8°N〜43.4°N, 2.6°E〜3.2°E]

C) 2025年8月：アメリカカリフォルニア州の火災 [34.5°N〜35.5°N, 119.8°W〜120.8°W]

1) 2020年5月20日：北海道滝上町のシバザクラ

2) 2020年10月21日：東北地方紅葉

3a) 2021年12月22日：人吉盆地の霧その1

3b) 2022年01月10日：人吉盆地の霧その2

4) 2022年1月7日：関東地方積雪

5) 2022年3月：阿蘇山周辺の野焼き

6) 2022年3月5日：渡良瀬の野焼き

7) 2023年3月29日：アラビア半島の灌漑農業

8) 2023年5月3日：西ヨーロッパ