PowerPoint プレゼンテーション

(1)

しきさい画像のGeoTIFF変換について

2020年10月16日 A

1 

しきさい画像は、画像の種類に応じたGDALコマンド、パラメータを用いて、

GeoTIFF変換できます。

以下の画像１～画像４のGeoTIFF変換例を紹介します。

【画像１】L1B、L2シーン画像（センサ観測座標系）

① L2 SST (海面水温) シーン画像の変換例

・・頁2

【画像２】HDFViewで出力したPNG画像（センサ観測座標系）

➁ HDFViewで出力したL2 SST(海面水温) PNG画像の変換例・・頁4

【画像３】L3画像（等緯度・経度座標系(EQR)）

➂ L3 NDVI (正規化植生指数) 画像の変換例

・・頁6

【画像４】L2タイル画像（等面積座標系(EQA)）

➃ L2 EVI (拡張植生指数)画像の変換例

・・頁7

A

付録：「しきさい」の投影座標系

(2)

① L2 SST (海面水温) シーン画像の変換例（センサ観測座標系）

１）GDALINFOでデータセット名を取得

Driver: HDF5/Hierarchical Data Format Release 5

Files: GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5 Size is 512, 512 ・・ Subdatasets: SUBDATASET_1_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Latitude SUBDATASET_1_DESC=[783x501] //Geometry_data/Latitude (32-bit floating-point)

SUBDATASET_2_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Longitude SUBDATASET_2_DESC=[783x501] //Geometry_data/Longitude (32-bit floating-point)

SUBDATASET_3_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Obs_time SUBDATASET_3_DESC=[783x501] //Geometry_data/Obs_time (16-bit integer)

SUBDATASET_4_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Sensor_zenith SUBDATASET_4_DESC=[783x501] //Geometry_data/Sensor_zenith (16-bit integer)

SUBDATASET_5_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Solar_zenith SUBDATASET_5_DESC=[783x501] //Geometry_data/Solar_zenith (16-bit integer)

SUBDATASET_6_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Image_data/QA_flag SUBDATASET_6_DESC=[7820x5000] //Image_data/QA_flag (16-bit unsigned integer)

SUBDATASET_7_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Image_data/SST SUBDATASET_7_DESC=[7820x5000] //Image_data/SST (16-bit unsigned integer)

gdalinfo GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5

２）GDAL_TRANSLATEでVRT (GDAL Virtual Format)へ変換

緯度 (Latitude)、経度 (longitude)及び海面水温データ(SST)それぞれのVRTファイルを

作成します。

gdal_translate -of VRT HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Latitude out_latitude.vrt

a.

b.

c.

a.

_{任意のファイル名.vrt} 画像ファイル名

ここでは、L2シーン画像のGeoTIFF変換例を紹介します。

A

(3)

3 ① L2 SST (海面水温) シーン画像の変換例（センサ観測座標系）

３）VRTファイルの編集

４）GDALWARPでGeoTIFF変換及び再投影

GDALWARPコマンドで、３）で編集したVRTファイルを

GeoTIFF変換及びEPSG:4326へ再投影します。

３）で保存したファイル名

２）で変換した海面水温データ(SST) のVRTファイルに緯度経度ファイル等の情報をメモ帳等

で追加後、上書き保存します。

<VRTDataset rasterXSize="5000" rasterYSize="7820">

<Metadata>

<MDI key="Geometry_data_Grid_interval">2500 </MDI>

<MDI key="Geometry_data_Grid_interval_unit">meter</MDI>

：

<VRTDataset rasterXSize="5000" rasterYSize="7820">

<Metadata domain="GEOLOCATION">

<MDI key="X_DATASET">

out_longitude.vrt

</MDI>

<MDI key="X_BAND">1</MDI>

<MDI key="Y_DATASET">

out_latitude.vrt

</MDI>

<MDI key="Y_BAND">1</MDI>

<MDI key="PIXEL_OFFSET">0.5</MDI>

<MDI key="LINE_OFFSET">0.5</MDI>

<MDI key="PIXEL_STEP">10</MDI>

<MDI key="LINE_STEP">10</MDI>

</Metadata>

<Metadata>

<MDI key="Geometry_data_Grid_interval">2500 </MDI>

<MDI key="Geometry_data_Grid_interval_unit">meter</MDI>

：

＜追加前＞

＜追加後＞

太字は、2)で作成した緯度、経度のファイル名

gdalwarp –of Gtiff -t_srs EPSG:4326 out_SST.vrt SST_output.tif

任意の出力ファイル名

＜QGISでの出力ファイル表示例＞

(4)

１）GDALINFOでデータセット名を取得

Files: GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5 Size is 512, 512 ・・ Subdatasets: SUBDATASET_1_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Latitude SUBDATASET_1_DESC=[783x501] //Geometry_data/Latitude (32-bit floating-point)

SUBDATASET_2_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Longitude SUBDATASET_2_DESC=[783x501] //Geometry_data/Longitude (32-bit floating-point)

・・

SUBDATASET_7_NAME=HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Image_data/SST SUBDATASET_7_DESC=[7820x5000] //Image_data/SST (16-bit unsigned integer)

gdalinfo GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5

２）GDAL_TRANSLATEでVRT (GDAL Virtual Format)へ変換

緯度 (Latitude)、経度 (longitude)及びHDFViewで出力したPNG (海面水温データ(SST))そ

れぞれのVRTファイルを作成します。

gdal_translate -of VRT HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Latitude out_latitude.vrt

a.

b.

a.

gdal_translate -of VRT HDF5:"GC1SG1_202006080125H04910_L2SG_SSTDQ_1002.h5"://Geometry_data/Longitude out_longitude.vrt

b.

元画像ファイル名

➁ HDFViewで出力したL2 SST(海面水温) PNG画像（センサ観測座標系）

ここでは、HDFViewで出力したPNG画像のGeoTIFF変換例を紹介します。

PNG画像の元画像ファイルの緯度、経度データセット名を取得します。

任意のファイル名.vrt

A

(5)

5 ３）VRTファイルの編集

４）GDALWARPでGeoTIFF変換及び再投影

GDALWARPコマンドで、VRTファイルを

GeoTIFF変換及びEPSG:4326へ再投影します。

３）で保存したファイル名

２）で変換したHDFViewで出力した海面水温データ(SST)のVRTファイルに緯度経度ファイル

等の情報をメモ帳等で追加後、上書き保存します。

<VRTDataset rasterXSize="5000" rasterYSize="7820">

<VRTRasterBand dataType="Byte" band="1">

<ColorInterp>Red</ColorInterp>

<SimpleSource>

<SourceFilename relativeToVRT="1">SST.png</SourceFilename>

：

<VRTDataset rasterXSize="5000" rasterYSize="7820">

<Metadata domain="GEOLOCATION">

<MDI key="X_DATASET">

out_longitude.vrt

</MDI>

<MDI key="X_BAND">1</MDI>

<MDI key="Y_DATASET">

out_latitude.vrt

</MDI>

<MDI key="Y_BAND">1</MDI>

<MDI key="PIXEL_OFFSET">0.5</MDI>

<MDI key="LINE_OFFSET">0.5</MDI>

<MDI key="PIXEL_STEP">10</MDI>

<MDI key="LINE_STEP">10</MDI>

</Metadata>

<VRTRasterBand dataType="Byte" band="1">

<ColorInterp>Red</ColorInterp>

<SimpleSource>

<SourceFilename

relativeToVRT="1">SST.png</SourceFilename>

：

＜追加前＞

＜追加後＞

太字は、2)で作成したファイル名

gdalwarp –of Gtiff -t_srs EPSG:4326 SST_png.vrt SST_png_output.tif

任意の出力ファイル名

＜QGISでの出力ファイル表示例＞

(6)

➂ L3 NDVI (正規化植生指数) 画像の変換例（等緯度・経度座標系(EQR)）

１）GDALINFOでデータセット名を取得

Files: GC1SG1_20200401D01M_D0000_3MSG_NDVIF_1001.h5 Size is 512, 512 ・・ Subdatasets: SUBDATASET_1_NAME=HDF5:"GC1SG1_20200401D01M_D0000_3MSG_NDVIF_1001.h5"://Image_data/NDVI_AVE SUBDATASET_1_DESC=[4320x8640] //Image_data/NDVI_AVE (16-bit unsigned integer)

SUBDATASET_2_NAME=HDF5:"GC1SG1_20200401D01M_D0000_3MSG_NDVIF_1001.h5"://Image_data/NDVI_QA_flag SUBDATASET_2_DESC=[4320x8640] //Image_data/NDVI_QA_flag (8-bit unsigned character)

gdalinfo GC1SG1_20200401D01M_D0000_3MSG_NDVIF_1001.h5

２）GDAL_TRANSLATEでGeoTIFFへ変換

以下のコマンドでデータセット名を指定し、GeoTIFF変換します。

gdal_translate -of Gtiff –a_srs EPSG:4326 –a_ullr -180 90 180 -90

HDF5:"GC1SG1_20200401D01M_D0000_3MSG_NDVIF_1001.h5"://Image_data/NDVI_AVE NDVI_output.tif

入力ファイルの参照座標系入力ファイルのデータセット名画像ファイル名入力ファイルの左上 (X, Y) 、右下 (X, Y) 出力ファイル名

＜QGISでの出力ファイル表示例＞

ここでは、L3画像のGeoTIFF変換例を紹介します。

A

(7)

7 ➃ L2 EVI (拡張植生指数) 画像の変換例（等面積座標系(EQA)/タイル）

２）GDALINFOでデータセット名を取得

Subdatasets:

SUBDATASET_1_NAME=HDF5:"GC1SG1_20200701D01M_T0428_L2SG_EVI_Q_2000.h5"://Image_data/EVI_AVE SUBDATASET_1_DESC=[4800x4800] //Image_data/EVI_AVE (16-bit unsigned integer)

gdalinfo GC1SG1_20200701D01M_T0428_L2SG_EVI_Q_2000.h5

３）GDAL_TRANSLATEでGeoTIFFへ変換

以下のコマンドでGCP点、データセット名を指定し、GeoTIFF変換します。

gdal_translate -of Gtiff -gcp 0 0 155.5707804 49.99895833 -gcp 4799 0 171.1244553 49.99895833 -gcp

4799 4799 143.5961321 40.00104167 -gcp 0 4799 130.5445343 40.00104167 -gcp 2400 0 163.3492384

49.99895833 -gcp 4799 2400 155.5595384 44.99895833 -gcp 2400 4799 137.071693 40.00104167 -gcp 0

2400 141.4205745 44.99895833 -gcp 2400 2400 148.4915296 44.99895833

HDF5:"GC1SG1_20200701D01M_T0428_L2SG_EVI_Q_2000.h5"://Image_data/EVI_AVE output.tif

入力ファイルのデータセット名画像ファイル名出力ファイル名

＜QGISでの出力ファイル表示例＞

ここでは、L2タイル画像のGeoTIFF変換例を紹介します。

A

タイル番号

１）タイル画像の緯度・経度計算

気候変動観測衛星「しきさい」 (GCOM-C) データ利用ハンドブックの「4.1.4.1 レベル 2 プロダクト作成単位」

から、対象タイル画像の分解能に応じて緯度・経度を計算します。

GCP計算箇所

← 9点の緯度・経度を GCPとして設定

EPSG:4326へ再投影します。

gdalwarp -of Gtiff –t_srs EPSG:4326 output.tif output2.tif

出力ファイル名入力ファイル名 (0, 0) (2400, 0) (4799, 0) (4799, 4799) (2400, 4799) (0, 4799) (0, 2400) (2400, 2400) (4799, 2400)

分解能Q(250m)例

(8)