GCOM-C SGLI
地図投影・GeoTIFF 出力ツール

取扱説明書
Version 1.2

Oct. 2023
JAXA
GCOM-C project

取扱説明書
改訂履歴 (1/1)
版
1.0
1.01.1

改訂
A-

日付
2021.72021.72022.3頁
24-58911

1.2

A

2022.112023.1052

内容
初版
0.5 画素位置ずれに係る記述を追加
オプション-a, -n, -z, -c を追加
オプション-c のファイル名定義を追加
オプション-a, -n, -z, -c に伴う TIFF タグ設定を変更
XML に Global_attributs の Scene/Image
start/end time を追加
オプション-c の説明を変更
4 項に L2 タイルのリサンプリング指定に係る注意事項を追加

目次
1． はじめに ............................................................................................. 12． 動作環境 ............................................................................................. 1
3． 入力プロダクト/対象データセット .................................................... 1
3.1 L1B シーン/半パスプロダクト ...............................................................................1
3.2 L2 シーンプロダクト ..............................................................................................2
3.3 L2 タイルプロダクト ...............................................................................................2

4． 操作説明 ............................................................................................. 2
5． 出力ファイル ...................................................................................... 75.1 ファイル名規約 ........................................................................................................75.2 GeoTIFF タグ/キー ................................................................................................95.3 XML 付加情報ファイル ......................................................................................... 11

6． ライセンス情報 ................................................................................. 12
7． 付録 .................................................................................................. 13
7.1 SGLI/HDF プロダクト画像データ地図投影方式 .................................................13

取扱説明書1.2 版, Oct. 2023

1． はじめに
本書は「GCOM-C/SGLI HDF プロダクトのための地図投影・GeoTIFF 出力ツール」のユーザ向け取扱説明書です。本ツールは HDF5 ファイルに格納されたラスタ形式の対象データについて、測地緯度経度あるいはポーラステレオ図法による地図投影を行い汎用の GeoTIFF フォーマットファイルとして出力します。ツールは CUI(キャラクターユーザインタフェース)形式であり、ユーザ自身のシステムやスクリプトへの組み込みによる大量処理に対応しています。

2． 動作環境
計算機環境
・CPUIntel Xeon/Core-i7 以上

・メモリ

1GByte 以上

対応 OS
・LinuxKernel 2.6 以上 (64bit)
動作試験は RHEL7/CentOS7 で実施

・Windows

Windows 10 (64bit)

3． 入力プロダクト/対象データセット
本ツールの処理対象となる SGLI HDF5 プロダクトは、L1B シーン/半パス、L2 シーン、L2 タイルです。

3.1

L1B シーン/半パスプロダクト

L1B の VNR/IRS シーンおよび POL 半パスのそれぞれ HDF5 プロダクトにおいて入力対象となるのは以下のデータセットです。
Geometry_data/
Obs_time_*, Sensor_*, Solar_*Image_data/
Lt_*, Land_water_flag, QA_flag
VNR, IRS ではシーン内のすべてのデータを投影しますが、POL では半パスにおいてオプション指定された緯度範囲のみを投影します。POL ではチルト駆動（チルト角切替）中の画像も投影可能ですが、投影されるのは駆動前後の画像であり、駆動中の画像は投影されません。

1取扱説明書1.2 版, Oct. 2023

3.2

L2 シーンプロダクト

L2 シーン(NWLR, IWPR, SST, OKID)のそれぞれ HDF5 プロダクトにおいて入力対象となるのは以下のデータセットです。
Geometry_data/
Sensor_*, Solar_*,Image_data/
NWLR_*, PAR, QA_flag, TAUA_*, CDOM, CHLA, TSM, SST,
Cloud_probability, OKID

3.3

L2 タイルプロダクト

L2 タイルプロダクトでは、HDF グループの Image_data か Geometry_data における、サイズが 4800 x 4800 或いは 1200 x 1200、形式が 8bit 符号なし或いは 16bit
符号あり/なし、の任意のデータセットが指定可能です。出力画像のフレームは、入力データセットの有効データ範囲に基づいて定義されます。指定されたデータセットに有効なデータがない場合、ファイルは出力されません。
L2 タイルプロダクト格納データセットの詳細については以下を参照ください。
https://suzaku.eorc.jaxa.jp/GCOM_C/data/product_std_j.html
なお、本ツールは L2 タイルプロダクトのバージョン 2003 までに含まれる 0.5 画素の位置ずれの補正は行いません。0.5 画素の位置ずれについては以下を参照ください。
https://shikisai.jaxa.jp/faq/faq0086_j.html

4． 操作説明ダウンロードしたツールの実行形式ファイル(.exe)を任意のディレクトリ/フォルダに置いてください。インストールは必要ありません。コマンドプロンプト(Linux ではターミナルウィンドウ)にて同実行形式を以下の引数で起動します。
> SGLI_geo_map_linux(win).exe [HDF ファイル名(*.h5)] [オプション]
注意事項：
L2 タイルでは、リサンプリングオプションをデフォルトの BL で処理すると画像端が無効画素となり、隣接タイルをモザイクする際は境界にギャップが発生することがあります。その場合は同オプションに NN を指定してください。

2

取扱説明書1.2 版, Oct. 2023
オプション:
オプション

説明
HDF5 内の処理対象データセットを DATASET に指定します-d DATASET(例: Image_data/Lt_VN01)。 L2 シーン、タイルでは本オプションは必須です。L1B で指定なしの場合はすべての TOA画像データ(Image_data/Lt_*)を処理します。

-o OUTPUTDIR

出力ファイルのディレクトリを OUTPUTDIR に指定します。
指定なしの場合カレントに出力します。
地図投影法をポーラステレオ(PS)に設定します。但し、シーン-pで赤道を跨ぐ場合はエラーとなります。デフォルトは測地緯度経度です。
ピクセルスペーシングを SPACING に設定します。単位は測地
緯度経度の場合は秒(arc-sec)、PS 図法の場合はメートルです。デフォルトは以下の通りです。

-s SPACING

解像度モード:

緯度経度

PS250m モード:

7.5 秒

250m

500m モード:

15 秒

500m

1km モード:

30 秒

1km

指定可能な範囲は、測地緯度経度では 7.5~180 秒、PS では250m ~ 6km です。IRS/L1B で全画像処理時(-d を省略した場合)は 1km 解像度モードデータに対する値を指定します（250m, 500m 解像度画像はそれぞれ指定値の 1/4, 1/2 で投影されます）。
リサンプリング法（NN, BL, CC）を RESAMPLING に指定した番号で選択します。番号の定義は以下の通りです。
0 = NN, 1 = BL, and 2 = CC.
-r RESAMPLING

デフォルトはフラグデータの場合は NN、以外は BL です。L2
タイルプロダクトの処理画像で隣接タイルをモザイクする場
合などは、画像端が無効画素となることによるタイル境界のギャップを防ぐために NN を指定してください。

-m

TOA 輝度データの DN 値 MSB２bits に含まれる迷光補正フラグをマスクします。

3

取扱説明書1.2 版, Oct. 2023
POL/L1B 入力でオプション-d が指定されなかった時に出力するデータ種別を TYPEPOL に指定した番号で選択します。番-b TYPEPOL

号の定義は以下の通りです。
-1 = ストークスデータのみ(Lt_PI*, Lt_PQ*, Lt_PU*)
1 = 偏光角データのみ (Lt_P1*, Lt_P2*)
0 = すべてPOL/L1B 入力時に投影緯度範囲の北端を指定します。POL-u UPPERLAT/L1B 入力時に本オプションは必須です。入力単位は 0.1 度です。
POL/L1B 入力時に投影緯度範囲の南端を指定します。POL-l LOWERLAT/L1B 入力時に本オプションは必須です。入力単位は 0.1 度です。北端との間で指定可能な緯度範囲幅は 5 度～15 度です。
POL/L1B でチルト駆動を含む画像範囲を指定した際に、投影する画像のチルト角を TILT に指定された番号で選択します。
番号の定義は以下の通りです。

-t TILT0 = チルト駆動前のチルト角の画像を投影
1 = チルト駆動後のチルト角の画像を投影
本オプションが指定されない場合（デフォルト）、チルト駆動区間の中心ラインの地上投影位置を境界としてチルト駆動前後のチルト角の画像が投影されます。
処理対象データセットの HDF アトリビュートに設定されたオフセット/スケール(スロープ)を METADATA に指定された種別に基づき GeoTIFF 拡張タグの GDAL_METADATA (42112)に設定します。本オプションは、gdal_transrate のオプショ-a METADATA

ン-a_scale/-a_offset と等価です。METADATA に指定する種別は以下の通りです。
default = HDF の Offset/Slope をセット
reflectance = HDF の Offset_reflectance/Slope_reflectance をセット
画像内の無効値として指定する NODATA を GeoTIFF 拡張タ
グの GDAL_NODATA (42113)に設定します(例: -n 65535)。

-n NO_DATA

本オプションは gdal_transrate のオプション-a_nodata と等価です。設定可能な無効値は出力ファイル毎に 1 つのみです。

-z

画像データの LZW による内部圧縮を適用します。

4

取扱説明書1.2 版, Oct. 2023
DATASET に指定された HDF5 内複数チャネルのデータをチャネル結合 GeoTIFF ファイルとして出力します。本オプションの対象は Image_data グループの以下です。
L1B or L2/LTOA = Lt_*と QA_flag,
L2/RSRF = Rs_*, Tb_*, Rp_*と QA_flag,
その他 = QA_flag とそれに該当する一つの uint16 画像
チャネル名は空白を開けずにコンマ(,)で区切って指定します。L1B および L2 タイルの LTOA/RSRF で指定可能なチャネ
ル名は以下の通りです。
VN01,VN02,VN03,VN04,VN05,VN06,VN07,
VN08,VN09,VN10,VN11,VN08P,VN11P,
SW01,SW02,SW03,SW04,TI01,TI02,
P1_p60,P1_0,P1_m60,P2_p60,P2_0,P2_m60,
PI01,PQ01,PU01,PI02,PQ02,PU02,PL01,PL02,-c DATASETQA_flag(例: -c VN03,VN05,VN08,QA_flag).
番号が昇順で連続するチャネル名の場合はハイフン(-)を使用可能です。
(例: -c VN01-11,SW01-03).
チャネル名先頭文字(V, S, T, P)以外の文字についてはワイルドカード(*,?)を使用可能です。
(例: -c "SW*,TI*", -c "P?01").
ピクセルスペーシングが異なる画像が指定された場合はピクセルスペーシングが最小のデータの投影画像に対して他の投影画像をオーバーサンプリングして格納します。
L2 シーンおよび LTOA/RSRF 以外の L2 タイルの場合は、対象の画像データと該当する QA_flag を指定します。
(例: -c SST,QA_flag)
オプション-c と-d が同時に指定された場合は-c が有効となり-d は無視されます。

-h

使用法を表示して終了します。

-v

バージョンおよびライセンス情報を表示して終了します。

5

取扱説明書1.2 版, Oct. 2023
起動例:
1) L1B/VNR プロダクト入力、迷光フラグマスク、リサンプリング法 NN、ピクセ
ルスペーシング 10 秒、出力ディレクトリ”./output”> SGLI_geo_map_linux.exe [L1B/VNR].h5 -m -r 0 -s 10.0 -o ./output
2) L1B/VNR プロダクト入力、Lt_VN11 指定、迷光フラグマスク> SGLI_geo_map_linux.exe [L1B/VNR].h5 -d Image_data/Lt_VN11 -m3) L1B/VNR プロダクト入力、Sensor_zenith 指定> SGLI_geo_map_linux.exe [L1B/VNR].h5 -d Geometry_data/Sensor_zenith
4) L1B/IRS プロダクト入力、Lt_TI01 指定、PS 図法> SGLI_geo_map_linux.exe [L1B/IRS].h5 -d Image_data/Lt_TI01 -p
5) L1B/POL プロダクト入力、投影緯度範囲 N45.0°~N30.0°> SGLI_geo_map_linux.exe [L1B/POL].h5 -u 45.0 -l 30.06) L2/LTOA プロダクト入力、Lt_VN11 指定、迷光フラグマスク> SGLI_geo_map_linux.exe [L2/LTOA].h5 -d Image_data/Lt_VN11 -m

6

取扱説明書1.2 版, Oct. 2023

5． 出力ファイル
処理完了後、入力 HDF5 内の指定データセット毎に地図投影済み画像データを格納した GeoTIFF 形式ファイル (.tif)、 および付加情報を格納した XML 形式ファイル(.xml)を出力します。画像データが測地緯度経度投影で経度 180 度を跨ぐ場合は、経度 180 度で東西 2 つのファイルに分割してそれぞれ出力します。

5.1

ファイル名規約

ファイル命名規約は以下の通りです。
L1B POL オプション–t 無し:
GID_ULAT_LLAT_SDS.tif/.xml
GID = 入力 HDF ファイルのグラニュール ID
ULAT = オプション-u で指定された投影北端緯度 (S90.0~N90.0)
LLAT = オプション-l で指定された投影南端緯度 (S90.0~N90.0)
SDS = 指定データセット名例:
GC1SG1_201806250054L04800_1BSG_POLDK_1002_N45.0_N30.0_Lt_P1_0.tif
GC1SG1_201806250054L04800_1BSG_POLDK_1002_N45.0_N30.0_Lt_P1_0.xml
L1B POL オプション–t 有り:
GID_ULAT_LLAT_TILT_SDS.tif/.xml
TILT = チルト駆動前或いは後のチルト角 (T-45~T+45)
例:
GC1SG1_~(上記サンプルと同じ)~_N45.0_N30.0_T+45_Lt_P1_0.tif
GC1SG1_~(上記サンプルと同じ)~_N45.0_N30.0_T+45_Lt_P1_0.xml
上記以外:
GID_SDS.tif/.xml例:
GC1SG1_201806250119R04810_1BSG_VNRDQ_1002_Lt_VN11.tif
GC1SG1_201806250119R04810_1BSG_VNRDQ_1002_Lt_VN11.xml

7

取扱説明書1.2 版, Oct. 2023

経度 180 度を跨いで 2 つに分割されたファイル:
経度 180 度を境に東側(~180 度)および西側(-180 度~)に分割されたファイルについて、それぞれ拡張子を除くファイル名末尾に“_E”および”_W”を付加します。
例:
GC1SG1_202107030044H04008_1BSG_IRSDQ_2003_Lt_TI01_E.tif
GC1SG1_202107030044H04008_1BSG_IRSDQ_2003_Lt_TI01_E.xml
GC1SG1_202107030044H04008_1BSG_IRSDQ_2003_Lt_TI01_W.tif
GC1SG1_202107030044H04008_1BSG_IRSDQ_2003_Lt_TI01_W.xml
チャネル結合オプション”-c”指定時の SDS 名:
オプションの DATASET に指定した文字列について、カンマおよびワイルドカード文字(*,?)がそれぞれʼ_ʼおよびʼxʼに置き換えられた後に、そのまま SDS 名に使用されます。L1B および L2 タイルの LTOA/RSRF ではデータ種別 ID(Lt_, Rs_, Tb_,Rp_)が種別毎チャネル名各先頭に挿入されます。QA_flag が指定された場合は、その文字列が SDS 名の最後に追加されます。
例:
オプション指定 -c VN03,VN05,VN08
GC1SG1_202103100138K05510_1BSG_VNRDQ_2002_Lt_VN03_VN05_VN08.tif
GC1SG1_202103100138K05510_1BSG_VNRDQ_2002_Lt_VN03_VN05_VN08.xml

オプション指定: -c "SW*,TI*,QA_flag"
GC1SG1_202103100138K05510_1BSG_IRSDQ_2002_Lt_SWx_Tix_QA_flag.tif
GC1SG1_202103100138K05510_1BSG_IRSDQ_2002_Lt_SWx_Tix_QA_flag.xml

オプション指定: -c "V*,S*,T*,P*"GC1SG1_20210310D01D_T0529_L2SG_RSRFQ_2000_Rs_Vx_Sx_Tb_Tx_Rs_Px.tif
GC1SG1_20210310D01D_T0529_L2SG_RSRFQ_2000_Rs_Vx_Sx_Tb_Tx_Rs_Px.xml

8

取扱説明書1.2 版, Oct. 2023

5.2

GeoTIFF タグ/キー

出力 GeoTIFF ファイルに設定されるタグ/キーについて表 5.2-1 および表 5.2-2 に示します。
表 5.2-1 TIFF タグ設定
タグ値

TIFFTAG_IMAGEWIDTHTIFFTAG_IMAGELENGTHTIFFTAG_BITSPERSAMPLE

image width
image height
8 or 16

TIFFTAG_COMPRESSIONCOMPRESSION_NONE orCOMPRESSION_LZW
PHOTOMETRIC_MINISBLACKORIENTATION_TOPLEFT
1 to n (if option -c specified)image height
SAMPLEFORMAT_INT orSAMPLEFORMAT_UINT
1 or 2 (if option -c specified)
xml data
text data
6 parameters of model tie point tag
3 parameters of model pixel scaletag
text data

TIFFTAG_PHOTOMETRICTIFFTAG_ORIENTATIONTIFFTAG_SAMPLESPERPIXELTIFFTAG_ROWSPERSTRIPTIFFTAG_SAMPLEFORMATTIFFTAG_PLANARCONFIGTIFFTAG_GDAL_METADATA*TIFFTAG_GDAL_NODATA*GTIFF_TIEPOINTSGTIFF_PIXELSCALEGTIFF_ASCIIPARAMS*) 拡張タグ

9

取扱説明書1.2 版, Oct. 2023
表 5.2-2 GeoTIFF キー設定
キー値

GTModelTypeGeoKey
GTRasterTypeGeoKey
GeographicTypeGeoKey
GeogAngularUnitsGeoKey
GTCitationGeoKey
PCSCitationGeoKey
GeogLinearUnitsGeoKey (PS only)
ProjectedCSTypeGeoKey (PS only)
ProjectionGeoKey (PS only)
ProjCoordTransGeoKey (PS only)
ProjLinearUnitsGeoKey (PS only)

ModelTypeProjected
RasterPixelIsArea
GCS_WGS_84Angular_Degree
text datatext data
Linear_Meter
KvUserDefined
KvUserDefined
CT_PolarStereographic
Linear_Meter
latitude of the image projection
plane(= 71 or -71)
longitude as the vertical axis of theimage frame (= 0)00

ProjNatOriginLatGeoKey (PS only)
ProjStraightVertPoleLongGeoKey (PSonly)
ProjFalseEastingGeoKey (PS only)
ProjFalseNorthingGeoKey (PS only)

GeoTIFF フォーマットの参照情報は以下の通りです。
GeoTIFF Format Specification
GeoTIFF Revision 1.0
Specification Version: 1.8.2
Last Modified: 28 December, 2000.
http://geotiff.maptools.org/spec/geotiffhome.html (accessed Jul. 1, 2021)

10

取扱説明書1.2 版, Oct. 2023

5.3XML 付加情報ファイル

出力 XML ファイルに設定されるタグ/内容について表 5.3-1 に示します。
表 5.3-1 出力 XML ファイル設定のタグ/内容
タグ
GeoTIFFAncillary
Process_information
Software_name
Software_version
Processing_time
Map_projection
Pixel_spacing
Resampling_method
MSB2bitsMask
Upper_latitude
Lower_latitude
Fixed_tilt_angle
Data_information
Granule_ID
Dataset_name /Dataset_nn_name
HDF_attributes
Global_attributes

内容

ツール名
ツールバージョン
処理日付/時刻（CPU の地方時）
地図投影法 (“Geodetic Latitude/Longitude” or“Polar Stereographic”)
ピクセルスペーシング (in degrees for Lat/Lonor in meters for PS)
リサンプリング法 (“NN”, “BL”, or “CC”)
迷光フラグマスク有無(“Applied” or “None”)
オプション-u による指定北端緯度(L1B/POL)
オプション-l による指定南端緯度(L1B/POL)
投影チルト角(“Pre-angle of tlit-driving” or“Post-angle of tilt-driving”) (L1B/POL のオプション–t 指定時)
入力 HDF プロダクトのグラニュール ID
指定 HDF データセット名
チャネル結合ファイルはチャネル番号(nn)を挿入

HDF グループにセットされた Scene/Image のstart time アトリビュートを格納、タグにはHDF のアトリビュート名をコピー
HDF グループにセットされた Scene/Image のScene_end_time/end time アトリビュートを格納、タグには HDF
Image_end_timeのアトリビュート名をコピー
Image_data or Geometry_data (指定データセットの上位グループ名)
attributes
HDF グループにセットされたすべてのアトリビ(HDF グループのア ュートを格納、タグには HDF のアトリビュート
トリビュート名)
名をコピー
SDS_name (指定データセット名)
attributes
HDF データセットにセットされたすべてのアト(HDF データセット
リビュートを格納、タグには HDF のアトリビュのアトリビュート
ート名をコピー名)
Scene_start_time/
Image_start_time

11

取扱説明書1.2 版, Oct. 2023

6． ライセンス情報
本ツールには HDF5 および GeoTIFF のそれぞれファイルインタフェースのためのオープンソースライブラリを使用しています。これらのライセンス情報は以下の通りです。
Copyright Notice for HDF5 (Hierarchical Data Format 5) Software Library========================================================
HDF5 (Hierarchical Data Format 5) Software Library and Utilities
Copyright (c) 2006-2018, The HDF Group.
NCSA HDF5 (Hierarchical Data Format 5) Software Library and Utilities
Copyright (c) 1998-2006, The Board of Trustees of the University of Illinois.
All rights reserved.

https://www.hdfgroup.org/licenses (accessed Jul. 1, 2021)

libgeotiff Licensing==============
Code by Frank Warmerdam has this copyright notice:* Copyright (c) 1999, Frank Warmerdam
Code by Niles Ritter is under this licence:*Written By: Niles D. Ritter.

* Copyright (c) 1995, Niles D. Ritter

https://github.com/OSGeo/libgeotiff/blob/master/libgeotiff/LICENSE(accessed Jul. 1, 2021)

12

取扱説明書1.2 版, Oct. 2023

7． 付録
7.1

SGLI/HDF プロダクト画像データ地図投影方式

本章では、ツール内で実施されている地図投影アルゴリズムの概要について記述します。

L1B シーン/半パスプロダクト
L1B シ ー ン /半 パ ス プ ロ ダ ク ト の 地 図 投 影 は 、 Rational Polynomial Coefficients(RPC)を内部的に作成することで行われます。RPC は地上空間の三次元座標位置(X, Y,Z)と画像の 2 次元座標位置(L, P)とを関連づける 3 次の有理多項式であり、最新のL2/LTOA 処理アルゴリズムにおけるオルソ投影でも使用されています。RPC の変換式について以下に示します。

𝑓𝑓 �𝑋𝑋 ′ ,𝑌𝑌 ′ ,𝑍𝑍 ′ �

𝑓𝑓 �𝑋𝑋 ′ ,𝑌𝑌 ′ ,𝑍𝑍 ′ �

𝐿𝐿′ = 𝑓𝑓1 (𝑋𝑋 ′ ,𝑌𝑌 ′ ,𝑍𝑍 ′ ), 𝑃𝑃′ = 𝑓𝑓3 (𝑋𝑋 ′ ,𝑌𝑌 ′ ,𝑍𝑍 ′)24𝑓𝑓(𝑋𝑋, 𝑌𝑌, 𝑍𝑍) = 𝑐𝑐1 + 𝑐𝑐2 𝑌𝑌 + 𝑐𝑐3 𝑋𝑋 + 𝑐𝑐4 𝑍𝑍 + 𝑐𝑐5 𝑋𝑋𝑋𝑋 + 𝑐𝑐6 𝑌𝑌𝑌𝑌+𝑐𝑐7 𝑍𝑍𝑍𝑍 + 𝑐𝑐8 𝑌𝑌 2 + 𝑐𝑐9 𝑋𝑋 2 +𝑐𝑐10 𝑍𝑍 2+ 𝑐𝑐11 𝑋𝑋𝑋𝑋𝑋𝑋 + 𝑐𝑐12 𝑌𝑌 3 + 𝑐𝑐13 𝑋𝑋 2 𝑌𝑌 + 𝑐𝑐14 𝑍𝑍 2 𝑌𝑌 + 𝑐𝑐15 𝑌𝑌 2 𝑋𝑋 + 𝑐𝑐16 𝑋𝑋 3 +𝑐𝑐17 𝑍𝑍 2 𝑋𝑋+ 𝑐𝑐18 𝑌𝑌 2 𝑍𝑍 + 𝑐𝑐19 𝑋𝑋 2 𝑍𝑍 + 𝑐𝑐20 𝑍𝑍 3𝐿𝐿′ =𝐿𝐿−𝐿𝐿𝑜𝑜
𝐿𝐿𝑠𝑠, 𝑃𝑃′ =𝑃𝑃−𝑃𝑃𝑜𝑜𝑃𝑃𝑠𝑠, 𝑋𝑋′ =𝑋𝑋−𝑋𝑋𝑜𝑜𝑋𝑋𝑠𝑠, 𝑌𝑌′ =𝑌𝑌−𝑌𝑌𝑜𝑜
𝑌𝑌𝑠𝑠, 𝑍𝑍′ =𝑍𝑍−𝑍𝑍𝑜𝑜𝑍𝑍.

RPC の作成は、L1B の HDF5 データに格納されている SGLI センサモデルパラメー
タ、GCOM-C 衛星軌道姿勢データ、各機器アライメントモデルパラメータ、によって構築される厳密な物理投影モデルに基づいて行われます。RPC 作成のための投影モデ
ルの概念図について図 7.1-1 に示します。SGLI シーン/半パスの比較的広い画像範囲（スワスおよびライン数）をカバーするために、RPC は分割した画像毎に複数作成されてそれぞれ投影に使用されます。

13

取扱説明書1.2 版, Oct. 2023
Pimage space
(L1B)

L

focal point(satellite) orbitZ
objectspace

6000m
elevation
layer0m

YX
Geodetic / Map Coordinates

図 7.1-1

RPC 作成のための投影モデル概念図

地図投影では、まず、選択された地図座標上の出力画像のフレームについて、シーン四隅の緯度経度を用いて定義します。なお、POL の半パスデータでは、四隅緯度経度の代わりにオプションで指定された緯度範囲に基づき同フレームを定義します。次いで、出力画像上に等間隔のグリッドを設定し、それぞれのグリッドに対応する入力画像上座標値について RPC を用いて算出します。この際、出力画像上のグリッドに該当する地図座標値を、RPC 入力として定義された地図座標値に変換します。
数値地形モデル(Digital Elevation Model: DEM)は入力されないため、RPC に入力する高さ値(Z)は 0(ゼロ)に固定され、画像データは準拠楕円体モデル(GRS80)上に投影されます。すなわち、DEM 入力によるオルソ(地形)補正が行われている LTOA 由来のL2 タイルプロダクトとは異なり、本ツールで出力される画像では同補正は行われません。
最後に、出力画像全画素に対応する入力画像座標値について、RPC で求めた等間隔グリッドの対応座標値からのバイリニア補間により決定し、これらのデジタル値を入力画像のデジタル値からリサンプリングします。
L2 シーンプロダクト
L2 シーンプロダクトには、RPC 作成に必要な物理投影モデルを構築するためのセンサモデルや衛星軌道姿勢等のデータが含まれていません。よって、代わりに HDF プロ
ダクトの画像データ上 10 画素間隔グリッドの値として Geometry_data グループに14

取扱説明書1.2 版, Oct. 2023
格納された緯度経度データを使用して投影を行います。シーン内の任意の緯度経度に該当する画像上のアドレスは、その緯度経度を内包する 4 つの緯度経度グリッドデータからバイリニア補間を用いて計算されます。
投影手順は、出力画像の等間隔グリッドに対応する入力画像座標値の計算以外は、同計算に RPC を用いた L1B シーン/半パスプロダクトと同じです。但し、同計算に緯度経度 10 画素間隔グリッドデータ上でのサーチ処理が入る L2 シーンプロダクトの処理時間は、同サーチ処理が不要な RPC を用いた L1B シーン/半パスプロダクトの処理時間に対して遅くなります。なお、L1B シーン/半パスプロダクト同様にオルソ(地形)補正は行われません。
L2 タイルプロダクト
L2 タイルプロダクトでは、入力タイルのサンソン図法と測地緯度経度座標との間の以下変換式を用いて直接的に地図投影を行います。

x=m(λ cos ϕ '−10h + 180)
+ 0.5 ,
10y=n(90 − 10v − ϕ )+ 0 .5 ,10ϕ ′ = ϕ ⋅ π / 180
m=n=4800 (250m mode)
1200 (1km mode).

ここで、(x, y)は画像座標値、(ϕ, λ)は測地緯度経度、(v, h)はタイルの垂直水平番号です。画像座標の原点は、タイル左上隅画素の中心=(1, 1)としています。
ポーラステレオ図法への投影時は、測地緯度経度との間の変換を追加して処理を行います。出力画像のフレームは、入力データセットの有効データ範囲に基づいて定義されます。

15