QGISを使ってGCOM-C(しきさい)の地表面温度データを画像で表示する方法
この資料では、QGIS（Quantum GIS）を使ってGCOM-C(しきさい)のレベル２の地表面温度プロダクトを読み込み、画像化する方法を示します。
QGISについては、以下のG-Portalページの「ツール」→「Quantum GIS」に記載されているリンク先をご参照ください。
https://gportal.jaxa.jp/gpr/information/tool

また、必要なデータをG-Portalで検索し、ダウンロードする操作については細かくは示しておりませんので、必要に応じて下記の「G-Portalユーザ向け取扱説明書」をご参照下さい。
https://gportal.jaxa.jp/gpr/assets/mng_upload/COMMON/upload/GPortalUserManual_ja.pdf

JAXA/G-Portalサポートデスク
z-gportal-support@ml.jaxa.jp

1

① QGISに読み込ませる地表面温度プロダクトを用意する
まず、QGISに読み込ませるための「GeoTIFF」形式の地表面温度プロダクトを用意します。
入手方法については、以下の2通りの手順があります。
A. G-Portalより、地表面温度プロダクトをGeoTIFF形式にする加工要求を実行し、ダウンロードする。
B. G-Portalより、地表面温度プロダクトをダウンロードし、 SGLI 地図投影・GeoTIFF出力ツールを使ってHDF5形式のプロダクトをGeoTIFF形式に変換する。

2

①-A(1/2) G-Portalより加工要求を実行しGeoTIFF形式のデータを入手する(1)

QGISで画像として表示させたいGCOM-C(しきさい)の地表面温度データをG-Portalで検索する。

(2)

検索結果が表示されたら、所望データの行にある「加工」ボタンをクリックする。

「加工」ボタンを押す※「ダウンロード」ボタンをクリックした場合、HDF5形式でデータがダウンロードされます。
3

①-A(2/2) G-Portalより加工要求を実行しGeoTIFF形式のデータを入手する(3) 以下のように加工要求を実行し、生産完了メールを受信後、GeoTIFFファイルをダウンロードする。

「GeoTIFF」を選択し、「加工要求」ボタンを押す。
4①-B(1/2) G-Portalから取得したHDF5プロダクトをGeoTIFFに変換する
コマンドプロンプトを立ち上げ、SGLI 地図投影・GeoTIFF出力ツール(Ver.1.2)を使ってHDF5のプロダクトをGeoTIFFに変換する。
ツールの詳細はG-Portal Webページ「ツール・ドキュメント⇒GCOM-Cタブ」にある取扱説明書を参照。
https://gportal.jaxa.jp/gpr/information/tool#GCOM-C
C:¥xxxxx¥Desktop¥GCOM-C_L2-LST>SGLI_geo_map_win.exe GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000.h5 -d Image_data/LST========================== Input parameters ==========================
Input file= GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000.h5
Projection coordinates = Geodetic Latitude/Longitude
Pixel spacing= Default (7.5sec/15sec/30sec)
例として2022年8月10日の日本域
Resampling method= BL
Data selection= Image_data/LSTのタイルデータを用います。
Output directory=.
======================================================================
product_type = L2/tile
tile = 0529
pid = GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000
hemi = 0
p=0, Image_data/LST
data rank = 2
data size = 4800/4800
data type = unsigned short
INFO: type of attr is int -> cast to unsigned short! function (H5_common_read2.c) line = 272
max_valid_DN = 65534
x range = 129.4125000000, 141.6000000000
y range = 30.0000000000, 38.8083333333
outsize = 5850, 4228
output = [.¥GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000_LST.tif]
output = [.¥GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000_LST.xml]
done...

5

①-B(2/2) G-Portalから取得したHDF5プロダクトをGeoTIFFに変換する
SGLI 地図投影・GeoTIFF出力ツールのコマンドオプションについて処理対象データセットの HDF アトリビュートに設定されたオフセット/スケール(スロープ)をMETADATAに指定された種別に基づき GeoTIFF 拡張タグの GDAL_METADATA (42112)に設定します。
上記オプションを指定した場合、本資料18ページ目に記載している計算式において、slope/offset 値の計算を省略することができます。
ただし、温度の単位がケルビン(K)となっていることにご注意ください。
【コマンド実行例】>SGLI_geo_map_win.exe GC1SG1_20220810D01D_T0529_L2SG_LST_Q_3000.h5 -d Image_data/LST –a default

オプションの詳細は前項に記載している取扱説明書の4ページを参照。
https://gportal.jaxa.jp/gpr/information/tool#GCOM-C

6

②GeoTIFFをフリーツールであるQGIS(https://qgis.org/ja/site/)で読み込み、画像化する<前準備 ～QGISで読み込む“陸地のデータ”を用意する～>
以下のサイトから陸地のデータ（QGISでデータを画像化する際、陸に色を塗るために使う）をダウンロードする。
https://tapiquen-sig.jimdofree.com/english-version/free-downloads/world/

7

ここをクリックし、右のような画面が出たら、「ファイルを保存する」を選び、「OK」を押す。（お使いのブラウザによって保存場所が異なるのでご注意ください。）

8

ファイルが保存された場所のフォルダを開き、「World_Countries.rar」というファイルを見つける。
次に、そのファイルをダブルクリックして展開する。すると、下記のようなファイルが出てくる。
（このうち「 World_Countries.shp 」を後程使用する。）
・World_Countries.dbf
・World_Countries.prj
・World_Countries.sbn
・World_Countries.sbx
・World_Countries.shp
・World_Countries.shp.xml
・World_Countries.shx
・World_Countries_README.txt

9③取得したデータをQGISで読み込ませる
１、QGISを起動します。（ここでは例としてバージョン3.10のQGISを使って説明しますが、お使いのバージョンにより画面構成が多少異なる可能性があります。）

10

「新規プロジェクト」をダブルクリック

11

手順①で作ったGeoTIFFファイル、および手順②で前準備として取得したファイルのうち「 World_Countries.shp 」をここにドラッグ＆ドロップする。

12

日本部分を拡大すると、このように表示される。

13

「ラスタ」から「ラスタ計算機」を選択する。

14

このような画面が開きます。
ここで行うことは、データに入っている値を地表面温度の値に変換するための“式を作る操作”です。

15

①①「バンド」の枠内にあるファイル名をダブルクリックする。（すると、下の「式」の枠内にそれが表示される。）

16

②「*」ボタン（×の意味）をクリックする。（すると、下の「式」の枠内にそれが表示される。）②

17

③「*」の右横にカーソルを当てて、1回クリックする。すると入力可能な状態になるので、「0.02」と入力する。
※データに格納された値にslope値（LSTの場合は0.02）を掛け、offset値（LSTの場合0）を足すことで温度の値に変換できます。ただし、そのままでは単位がKですので、273を引くことで℃に変換できます。

③

https://suzaku.eorc.jaxa.jp/GCOM_C/data/update/Algorithm_LST_ja.html

18

④「-」ボタンをクリックする。
⑤「式」の欄に「-」が表示されたら、その右側を1回クリックし、入力可能な状態とする。
⑥「273」と入力する。

④

⑤
⑥19

⑦

⑦「ラスタレイヤ」欄の「出力レイヤ」の右端にある「…」を1回クリックする。

20

⑧適当な場所（デスクトップ上でも可）を指定し、適当なファイル名（例えば「GCOM-C_LST」など）を付けて保存ボ
タンを押す。

21

⑨「OK」ボタンを押す。

⑨22

しばらく計算処理が動いたのち、このような画面が表示される。

23

ここからは、画像に色をつける方法の説明です。

先程保存したファイル名の部分をダブルクリックする。
24

このようなウインドウが開く。

25

「レンダリングタイプ」から「単バンド疑似カラー」を選択する。

26

細かく見たい地表面温度の値をそれぞれ入力する。この例では、最小値（最低温度）を20℃、最大値（最高温度）を60℃に設定している。

27

この状態だと、地表面温度を画像化した際、温度が低いほど赤、高いほど青に近い色になり、イメージと逆になってしまう。
この赤～青を逆転させるためには、カラーランプの右側にある「▼」を押し、表示された項目の中にある「カラーランプを反転」を押す。

28

低い温度が青、高い温度が赤になった。

ここにチェックを入れる。

29

ここにある「輝度」「彩度」「コントラスト」を調整することで、カラーランプの濃淡などを細かく調整することもできます。

30

「適用」をクリックし、「OK」をクリックする。

31

結果として、このような色が着いた画像が表示される。
※色合いは多少変更しています。

変換前の値が表示されないように元データのチェックを外します。
32

※この先の説明は、色の濃淡の分かりやすさを考慮して、濃い色合いにしたものをサンプルとして使います。
この「＋」と「ー」でズームできます。
また、その2つ左にある「手のマーク」をクリックすることで、画像の表示範囲を動かすことができます。

33

このようにして作成した地表面温度のデータを保存するには、メニューにある「プロジェクト」から「名前をつけて保存」を選択し、適当な名前を付けて適当なフォルダに保存する。

※次に同じファイル開きたいときは、保存したファイルをダブルクリックすれば、QGISが自動的に起動して開きます。

34

表示している範囲を画像ファイルとして出力したい場合には、メニューにある「プロジェクト」から「インポートとエクスポート」を選択し、「地図を画像にエクスポート」をクリックする。

35

適当な名前を付けて「保存」ボタンを押す。

「保存」ボタンをクリックする。

36

最後に、メニューにある「プロジェクト」から「QGISを終了」を選択する。

37