- TOP >
- SHIKISAI GeoTIFF >
- SHIKISAI GeoTIFF-Land
SHIKISAI GeoTIFF-Land | Ocean | Atmosphere | Cryosphere | Level1B
Level-2 / Level-3 Land products
Click on the product name to see the details of the product.
If there are SD arrays (data sets) for multiple observation wavelengths, a representative one is shown.
VGI (Vegetation index) [NDVI/ EVI/ SDI]
Level-2 products are provided as a single HDF5 file including EVI (Enhanced Vegetation Index), SDI (Shadow Index), and NDVI (Normalized Difference Vegetation Index) by tile units.
Level-2 statistical products are provided as separate geophysical products, each in tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
Level-3 products are provided as individual physical quantity products.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
JASMES: NetCDF4
Physical quantities
NDVI (Normalized Difference Vegetation Index)
Indices indicating vegetation cover
Resolution
Level-2: 250m
Level-2 statistics: 250m
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
Physical quantities
EVI (Enhanced Vegetation Index)
Indices indicating vegetation activity
Resolution
Level-2: 250m
Level-2 statistics: 1km
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
Physical quantities
SDI (ShaDow Index)
An index indicating shadow fraction of vegetation area inferred from spectral reflectance.
Resolution
Level-2: 250m
Level-2 statistics: 250m
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
AGB (Above Ground Biomass) [AGB/ VRI]）
Level-2 products are provided as a single HDF5 file including VRI (Vegetation Roughness Index) and AGB (Above Ground Biomass) by tile units.
Level-2 statistical products are provided as separate geophysical products, by tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
Level-3 products are provided as individual physical quantity products.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
JASMES: NetCDF4
Physical quantities
AGB (Above Ground Biomass)
Dry weight of above-ground vegetation
Resolution
Level-2: 1km
Level-2 statistics: 1km
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
Physical quantities
VRI (Vegetation Roughness Index)
An index indicating vertical plant structure observed from multiple angles
Resolution
Level-2: 1km
Level-2 statistics: 1km
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
LAI (Leaf Area Index) [FPAR/ LAI]
Level-2 products are provided as a single HDF5 file including FPAR (Fraction of absorbed PAR) and LAI (Leaf Area Index) by tile units.
Level2 statistical products are provided as separate geophysical products, by tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
Level-3 products are provided as individual physical quantity products.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
JASMES: NetCDF4
Physical quantities
FPAR (Fraction of absorbed PAR)
Fraction of photosynthetically active radiation absorbed by vegetation
Resolution
Level-2: 250m
Level-2 statistics: 1km
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
Physical quantities
LAI (Leaf Area Index)
The sum of the one-sided green leaf area per unit ground area
Resolution
Level-2: 250m
Level-2 statistics: 250m
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
LST (Land Surface Temperature)
Level-2 and Level-2 statistical products are provided as HDF5 files by tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
Level-3 products are provided as HDF5 files of global products.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
JASMES: NetCDF4
Physical quantities
LST (Land Surface Temperature)
The temperature of the terrestrial land surface.
Resolution
Level-2: 250m
Level-2 statistics: 250m
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
RSRF (Land Surface Reflectance)
Level-2 products and Level2 statistical products are provided as a single HDF5 file including visible near-infrared, short-wavelength infrared, etc. by tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
Level-3 products are provided as a global area product for each geophysical quantity.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
Physical quantities
RSRF (Atmospheric corrected reflectance)
Land surface reflectance corrected for the effects of atmospheric scattering and absorption.
Resolution
Level-2: 250m
Level-2 statistics: 250m (RT01-02, RV01-11, GEOI, GEOV),
1km (RS01-04, RP01-02, RN08, RN11, GEOP)
Level-3: 1/24deg (4.6km)
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
Level-3
(daily statistics)
Level-3
(8-day statistics)
Level-3
(1-month statistics)
LTOA (Top of atmosphere radiance)
Level-2 are provided as a single HDF5 file including visible near-infrared, short-wavelength infrared, etc.by tile units.
Level2 statistical products are provided are provided as separate geophysical products, by tile units.
For more information on tile products, see FAQ https://shikisai.jaxa.jp/faq/faq0080.html.
G-Portal: HDF5
G-Portal Processing Request: GeoTIFF(Individual files)
JASMES: NetCDF4
Physical quantities
LTOA (Top of atmosphere radiance)
Top of atmosphere (TOA) radiance
Resolution
Level-2: 250m, 1km
Level-2 statistics: 250m
Level-2
Level-2 statistics
(8-day)
Level-2 statistics
(1-month)
The following methods are available to use GeoTIFF files of “SHIKISAI” Level-2/Level-3 data.
* How to convert standard data obtained from G-Portal
- -CUI: SGLI Map projection & GeoTIFF conversion Tool (The link to the download site is provided from "Tools " below.)
- -GUI: Earth Observation Data Conversion Tool (Click on the tool name to go to the download page.)
* How to convert netCDF data obtained from the JASMES site to GeoTIFF
- JASMES netCDF data can be read by the Add Layer function in QGIS.
- GeoTIFF can be converted by QGIS's format conversion function, etc. Global data is centered at 180deg longitude, therefore a reference coordinate system must be set if necessary.
- Please see the following FAQ for information on converting on the command line.
- https://shikisai.jaxa.jp/faq/faq0092.html
* How to obtain GeoTIFF from G-Portal processing requests
- Request processing of GeoTIFF data from the "Processing" button in the data search results.
- Please refer to the following FAQs for the products covered.
- https://shikisai.jaxa.jp/faq/faq0044.html
Tools
This section provides examples of GeoTIFF conversions using the following tools.
Click on the tool name to go to the download page.
| Target product | Tool name | |
|---|---|---|
Level-2 Level-2 statistics |
SGLI Map projection & GeoTIFF conversion Tool |
Link to the "Tools & Documents" page in G-Portal. Platforms: Windows, Linux |
Level-3 |
QGIS(Quantum GIS)/gdal | Link to an external site. Platforms: Windows, Linux |
This section explains converting a GeoTIFF file using the "Tools" above and the HDF5 sample data below.
Data used in the explanation
Please obtain the following data from G-Portal.
Level-2 LST:
Observed 08/07/2022 Hokkaido, Japan (Tile No.0428) Nighttime
GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000.h5 (4MB)
Level-2 LTOA:
Observed 05/03/2022 Around the Sea of Japan (Tile No.0528) Daytime
GC1SG1_20220503D01D_T0528_L2SG_LTOAK_2011.h5 (42MB)
Level-3 LST:
Monthly statistics for August 2022 Global Nighttime
GC1SG1_20220801A01D_D0000_3MSG_LST_F_3000.h5 (11MB)
1) Command for Level-2/Level-2 statistic products (single SD array) conversion
Using the Windows command prompt or the equivalent, go to the directory where "SGLI_geo_map_win.exe" is located.
(The Windows command prompt can be accessed from Windows Start > Windows System > Command Prompt)
Below is an example where SGLI_geo_map_win.exe and the data folder containing the data are in the sgli_data folder directly on the C drive.
The command is an example of map projection of an SD array LST in Image_data and output as an LZW compressed GeoTIFF file.
The pixel value of 65535 is set for an invalid value, and the slope and offset in the attribute (text information) are stored in the GDAL_METADATA tag.
The values can be found in the <Slope> and <Offset> tags in the XML file output with the GeoTIFF conversion.
>cd C:\sgli_data
>SGLI_geo_map_win.exe C:\sgli_data\data\GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000.h5 -d Image_data/LST -r 0 -a default -n 65535 -z -o C:\sgli_data
For details on optional commands, please refer to the User's manual.
- SGLI_geo_map_win.exe: SGLI Map projection & GeoTIFF conversion Tool for Windows
- GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000.h5: Product name obtained from G-Portal
- -d Image_data/LST: SD array name (dataset name) in Image_data to be converted
- -r 0: Resampling method (NN)
- -a default: Set conversion coefficients (slope/offset) for physical quantities. If the software supports it, like QGIS, the converted values will be displayed (the data will retain their original values).
- -n 65535: Set an invalid value
- -z: Apply LZW compression
- -o C:\sgli_data: Set any output destination
See "GeoTIFF Tag Information" for the tag information stored in the converted GeoTIFF file.
Converted image displayed in QGIS
2) Combining multiple SD arrays of Level-2/Level-2 statistic products into one file
Using the Windows command prompt or the equivalent, go to the directory where "SGLI_geo_map_win.exe" is located.
(The Windows command prompt can be accessed from Windows Start > Windows System > Command Prompt)
Below is an example where SGLI_geo_map_win.exe and the data folder containing the data are in the sgli_data folder on the C drive.
The command is an example of map projection of SD arrays Lt_VN03, Lt_VN05, and Lt_VN08 in Image_data and output as a single LZW compressed GeoTIFF file.
The pixel value of 65535 is set for an invalid value, and the slope and offset values in the attribute (text information) are stored in the GDAL_METADATA tag.
The values are in the tags in the XML file output with the GeoTIFF conversion.
Multiple arrays can be output to a single GeoTIFF file using the QA_flag on LTOA/RSRF products or single files.
* To output two files, one geophysical and one QA_flag, in a single GeoTIFF and display them in QGIS, set Layer Properties > Symbology > Rendering Type to single band gray.
>cd C:\sgli_data
>SGLI_geo_map_win.exe C:\sgli_data\data\GC1SG1_20220503D01D_T0528_L2SG_LTOAK_2011.h5 -c VN08,VN05,VN03 -r 0 -a reflectance -n 65535 -z -o C:\sgli_data
For details on optional commands, please refer to the User's manual.
- SGLI_geo_map_win.exe: SGLI Map projection & GeoTIFF conversion Tool for Windows
- GC1SG1_20220503D01D_T0528_L2SG_LTOAK_2011.h5: Product name obtained from G-Portal
- -c VN08,VN05,VN03: SD array names (dataset names) in Image_data to be converted
- -r 0: Resampling method (NN)
- -a reflectance: Set the reflectance conversion coefficients (slope_reflectance, offset_reflectance). If the software supports it, like QGIS, the converted values will be displayed (the data will retain their original values).
- -n 65535: Set an invalid value
- -z: Apply LZW compression
- -o C:\sgli_data: Set any output destination
See "GeoTIFF Tag Information" for the tag information stored in the converted GeoTIFF file.
Converted image displayed in QGIS
3) The gdal command during Level-3 product conversion
To use the gdal command, start OSgeo4W Shell.
(The OSgeo4W Shell can be accessed from Windows Start > QGIS 3.22.8 > OSgeo4W Shell for a typical installation. The above GeoTIFF output tool is also available from the OSgeo4W Shell.)
Navigate to the directory where the data to be converted to GeoTIFF is stored.
Below is an example where the data to be converted is in the sgli_data folder on the C drive.
The command is an example of map projection of the SD array LST_AVE in Image_data and output as an LZW compressed GeoTIFF file.
The pixel value of 65535 is set for an invalid value, and the slope and offset values in the attribute (text information) are stored in the GDAL_METADATA tag.
The slope and offset can be checked with the gdalinfo command, HDFView (tool download required), etc.
For example, please refer to "How to Confirm SD array name of Level-2/Level-3 land products."
>cd C:\sgli_data
>gdal_translate -of GTiff -a_srs EPSG:4326 -a_ullr -180 90 180 -90 -a_nodata 65535 -a_scale 0.02 -a_offset 0 -co COMPRESS=LZW HDF5:"GC1SG1_20220801A01D_D0000_3MSG_LST_F_3000.h5"://Image_data/LST_AVE L3LST_A01M.tif
- gdal_translate: Format conversion commands
- -of GTiff: Set GeoTIFF output
- -a_srs EPSG:4326: Set the reference coordinate system for the input file.
- -a_ullr -180 90 180 -90: Set the upper left and lower right coordinate values (longitude and latitude).
- -a_nodata 65535: Set an invalid value
- -a_scale, -a_offset: Set the slope and offset values.
- -co COMPRESS=LZW: Set LZW compression.
- HDF5:"GC1SG1_20220801A01D_D0000_3MSG_LST_F_3000.h5"://Image_data/LST_AVE: Information on SUBDATASET_X_NAME to be converted when displayed in gdalinfo (X: number)
- L3LST_A01M.tif: Sets the output file name.
Converted image displayed in QGIS
4) The gdal command for applying QA_flag to Level-2/Level-3 products
Here is an example of applying QA_flag information to LST data created with the "SGLI Map projection & GeoTIFF conversion Tool."
Convert the LST and QA_flag to GeoTIFF using option "-c" in the "SGLI Map projection & GeoTIFF conversion Tool."
The following command extracts data for which the 10-bit (residual 2K or more) and 15th-bit (no input data) flags in one pixel set in QA_flag of LST are "0." The flags are "0" or "1."
Information on LST's QA_flag can be found in the "Products and Algorithms" section of the EORC GCOM-C website (https://suzaku.eorc.jaxa.jp/GCOM_C/data/update/Algorithm_LST_en.html).
>cd C:\sgli_data
>gdal_calc -A GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000_LST_QA_flag.tif --A_band=2 --outfile=LST0428A01D_20220807QA.tif --calc="(((A >> 10) & 1) == 0) * ((A >> 15) == 0) * 1" --NoDataValue=65535
- gdal_calc: gdal command to perform calculations in the data
- -A: Set the file name used in the formula to "A."
- GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000_LST_QA_flag.tif: File name to be calculated
- --A_band=2: Set the order in GeoTIFF where QA_flag is stored
- --outfile: Set any output file name for calculation results
- --calc: Set up calculation formulas
- (A >> 10) & 1: Shift the bit information by 10 bits to the right. Leave only the upper bits. Also, do an AND calculation of the bits and extract the 10th bit.
- A >> 15: Shift the bit information by 15 bits to the right.
- (((A >> 10) & 1) == 0) * ((A >> 15) == 0) * 1: Returns 1 if the parentheses are satisfied. The last "*1" is set to "1" if the condition is met.
- --NoDataValue: Set an invalid value
Example of bit calculation
* A >> 10 or right_shift(A, 10) ------ If A has the pixel value (decimal: 18431) "01000111111111111111", it is shifted right by 10 bits. The data will be for "010001".
* A & 1 or bitwise_and(A, 1) ------ Performs a bitwise operation of "1" (decimal: 1) on the pixel value of A. Extract the 10th bit of the original data.
| 010001 | (17) | ||
| & | 000001 | (1) | |
| --------------------- | |||
| 000001 | (1) | ||
(0 and 0 --> 0, 1 and 1 --> 1, 0 and 1 --> 0)
Create LST data using the created QA_flag data.
>gdal_calc -A GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000_LST_QA_flag.tif --A_band=1 -B LST0428A01D_20220807QA.tif --outfile=LST0428A01D_20220807addQA.tif --calc="(A <= 65534) * (B == 1) * A + (A <= 65534) * (B == 0) * 65535 + (A > 65534) * (B == 1) * 65535 + (A > 65534) * (B == 0) * 65535" --NoDataValue=65535
- gdal_calc: gdal command to perform calculations in the data
- -A, -B: Set "A" and "B" as file names to be used in formulas
- GC1SG1_20220807A01D_T0428_L2SG_LST_Q_3000_LST_QA_flag.tif: Set to "A" the name of the file to be calculated
- LST0428A01D_20220807QA.tif: Set the file name from which flags are extracted to "B."
- --A_band=1: Set the order in GeoTIFF where LST is stored
- --outfile: Set any output file name for calculation results
- --calc: Set up calculation formulas
- (A <= 65534) * (B == 1) * A: If the A data is within the valid range and the B data is "1", the A value is left unchanged. Calculations using slope/offset are also possible.
- (A <= 65534) * (B == 0) * 65535: Data, where A data are within the valid range and B data are "0," are assigned to the invalid value.
- (A > 65534) * (B == 1) * 65535: The data where A data are outside the valid value and B data are "1" are assigned to the invalid value.
- (A > 65534) * (B == 0) * 65535: Data where A data are outside the valid values and B data are "0" are assigned to invalid values.
Set slope/offset with gdal_edit.
>gdal_edit -scale 0.002 -offset 0 LST0428A01D_20220807addQA.tif
QGIS display of LST image with QA_flag applied
The application of QA_flag for Level-3 products is the same as for Level-2. Apply QA_flag as necessary.
The following is an example of excluding pixels with no input data and pixels with error flags.
>gdal_translate -of GTiff -a_srs EPSG:4326 -a_ullr -180 90 180 -90 -co COMPRESS=LZW HDF5:"GC1SG1_20220801A01M_D0000_3MSG_LST_F_3000.h5"://Image_data/LST_QA_flag L3LST_A01M_QA_flag.tif
>gdal_calc -A L3LST_A01M.tif -B L3LST_A01M_QA_flag_cal.tif --outfile=L3LST_A01M_addQA.tif --calc="(A <= 65534) * (B == 1) * A + (A <= 65534) * (B == 0) * 65535 + (A > 65534) * (B == 1) * 65535 + (A > 65534) * (B == 0) * 65535" --NoDataValue=65535
- gdal_translate: Format conversion commands.
For a description of the commands, see "3) gdal command example for Level-3 product conversion." - gdal_calc: gdal command to perform calculations in the data
- -A, -B: Set "A" and "B" as file names to be used in formulas
- L3LST_A01M.tif: Set to "A" the name of the file to be calculated
- L3LST_A01M_QA_flag_cal.tif: Set the file name from which flags are extracted to "B."
- --outfile: Set any output file name for calculation results
- --calc: Set up calculation formulas
- (A <= 65534) * (B == 1) * A: If the A data are within the valid range and the B data are "1," the A value is left unchanged. Calculations using slope/offset are also possible.
- (A <= 65534) * (B == 0) * 65535: Data, where A data are within the valid range and B data are "0," are assigned to the invalid value.
- (A > 65534) * (B == 1) * 65535: The data where A data is outside the valid value and B data is "1" is assigned to the invalid value.
- (A > 65534) * (B == 0) * 65535: Data where A data are outside the valid values and B data are "0" are assigned to invalid values.
Set slope/offset with gdal_edit.
>gdal_edit -scale 0.02 -offset 0 L3LST_A01M_addQA.tif
QGIS display of L3 LST image with QA_flag applied
5) How to Confirm SD array names of Level-2/Level-3 Land products
For SD array names, please refer to "SD Array Names of Level-2/Level-3 Land Products" below.
For more information, please visit the EORC GCOM-C "Standard Products and Algorithms" website.
The SD array can be checked in the following ways.
For GeoTIFF conversion of Level 3 products, please use gdalinfo of gdal.
| gdalinfo | HDFView |
|---|---|
If you have a typical QGIS installation, please use the OSGeo4W Shell. (Start Menu > QGIS 3.22.8 > OSGeo4W Shell) for version 3.22 Go to the directory containing the data to be converted and enter gdalinfo followed by the file name. Image_data_xxx_Offset, Slope corresponds to offset, scale in gdal. After "Image_data/" in SUBDATASET is the SD array name; for conversion in gdal, specify the underlined part. |
Example of Level-3 LST display. (1) SD array name in Image_data (2) SD array name in Geometry_data (only folder if not included) (3) Additional information on LST_AVE, slope, offset information, etc. |
>gdalinfo GC1SG1_20220801A01M_D0000_3MSG_LST_F_3000.h5
Image_data_LST_AVE_Error_DN=65535 |
|
6）SD array names of Level-2/Level-3 Land products
| Product name | Physical quantity | ID | SD array name(dataset name) | |
|---|---|---|---|---|
| Level-2 | LTOA (Top of atmosphere radiance) |
LTOA (Top of atmosphere radiance) |
LTOA | Land_water_flag Lt_P(1,2)_(0,m60,p60) Lt_P(I,Q,U)(01,02) Lt_SW(01-04) Lt_TI(01,02) Lt_VN(01-11) Lt_VN(08,11)P Statistic_data_(SWI,TIR,VNI) |
RSRF (Atmospheric corrected reflectance) |
RSRF (Atmospheric corrected reflectance) |
RSRF | Land_water_flag Angstrom PAR Rp_PL(01,02) Rs_PI(01,02) Rs_SW(01-04) Rs_VN(01-11) Rs_VN(08,11)P SWR Tau_500 Tb_TI(01,02) |
|
VGI (Vegetation index) |
EVI (Enhanced Vegetation Index) |
VGI_ | EVI |
|
NDVI (Normalized Difference Vegetation Index) |
NDVI |
|||
SDI (Shadow Index) |
SDI |
|||
AGB (Above Ground Biomass) |
AGB (Above Ground Biomass) |
AGB_ | AGB |
|
VRI (Vegetation Roughness Index) |
VRI |
|||
LAI (Leaf Area Index) |
FPAR (Fraction of absorbed PAR) |
LAI_ | FAPAR |
|
LAI (Leaf Area Index) |
LAI Overstory_LAI |
|||
| LST (Land Surface Temperature) | LST (Land Surface Temperature) |
LST_ | E01 E02 LST |
|
| Level-2 statistics | LTOA (Top of atmosphere radiance) |
LTOA (Top of atmosphere radiance) |
LTOA | Lt_P(I,Q,U)(01,02) Lt_SW(01-04) Lt_TI(01,02) Lt_VN(01-11) Lt_VN(08,11)P |
RSRF (Land Surface Reflectance) |
RSRF (Land Surface Reflectance) |
RV(01-11) | Rs_VN(01-11)_AVE |
|
| RS(01-04) | Rs_SW(01-04)_AVE |
|||
| RT(01,02) | Tb_TI(01,02)_AVE |
|||
| GEO(V,I,P) | Relative_azimuth_AVE Sensor_zenith_AVE Solar_zenith_AVE |
|||
| RN(08,11) | Rs_VN(08,11)P_AVE |
|||
| RP(01,02) | Rs_PI(01,02)_AVE |
|||
| SWR_ | SWR_AVE |
|||
| EVI (Enhanced Vegetation Index) | EVI (Enhanced Vegetation Index) |
EVI_ | EVI_AVE |
|
| NDVI (Normalized Difference Vegetation Index) | NDVI (Normalized Difference Vegetation Index) |
NDVI | NDVI_AVE |
|
| SDI (Shadow Index) | SDI (Shadow Index) |
SDI_ | SDI_AVE |
|
| AGB (Above Ground Biomass) | AGB (Above Ground Biomass) |
AGB_ | AGB_AVE |
|
| VRI (Vegetation Roughness Index) | VRI (Vegetation Roughness Index) |
VRI_ | VRI_AVE |
|
| FPAR (Fraction of absorbed PAR) | FPAR (Fraction of absorbed PAR) |
FPAR | FAPAR_AVE |
|
| LAI (Leaf Area Index) | LAI (Leaf Area Index) |
LAI_ | LAI_AVE |
|
| LST (Land Surface Temperature) | LST (Land Surface Temperature) |
LST_ | LST_AVE |
|
| Level-3 | RSRF (Atmospheric corrected reflectance) | RSRF (Atmospheric corrected reflectance) |
RP(01,02) | Rs_PI(01,02)_AVE |
| RS(01-04) | Rs_SW(01-04)_AVE |
|||
| RV(01-11) | Rs_VN(01-11)_AVE | |||
| RN(08,11) | Rs_VN(08,11)P_AVE | |||
| RT(01,02) | Tb_TI(01,02)_AVE | |||
| SNZV | Sensor_zenith_AVE | |||
| SNZ(I,P) | Sensor_zenith_(IR,PL)_AVE | |||
| SLZ(I,V) | Solar_zenith_AVE | |||
| SLZP | Solar_zenith_PL_AVE | |||
| RLAV | Absolute_relative_azimuth_AVE | |||
| RLA(I,P) | Absolute_relative_azimuth_(IR,PL)_AVE | |||
| EVI (Enhanced Vegetation Index) | EVI (Enhanced Vegetation Index) |
EVI | EVI_AVE | |
| NDVI (Normalized Difference Vegetation Index) | NDVI (Normalized Difference Vegetation Index) |
NDVI | NDVI_AVE | |
| SDI (Shadow Index) | SDI (Shadow Index) |
SDI | SDI_AVE | |
| AGB (Above Ground Biomass) | AGB (Above Ground Biomass) |
AGB | AGB_AVE | |
| VRI (Vegetation Roughness Index) | VRI (Vegetation Roughness Index) |
VRI | VRI_AVE | |
| FPAR (Fraction of absorbed PAR) | FPAR (Fraction of absorbed PAR) |
FPAR | FAPAR_AVE | |
| LAI (Leaf Area Index) | LAI (Leaf Area Index) |
LAI | LAI_AVE | |
| LST (Land Surface Temperature) | LST (Land Surface Temperature) |
LST | LST_AVE |
7) Converting to a physical quantity
When converted to GeoTIFF using the SGLI Map projection & GeoTIFF conversion Tool, the coefficients (slope/offset) of the geophysical conversion are output to an XML file created at the same time.
Display in a text editor or browser.
Below is an example of LST XML.
The values of <Slope> and <Offset> are substituted for "Pixel Value * Slope + Offset" to convert the pixel value into a physical quantity. The unit is Kelvin [K].
<Data_description>Sea Surface Temperature[SST]: SST[degree]=DN*Slope+Offset
<Slope>2.000000e-002</Slope>
<Offset>0.000000e+000</Offset>
<Unit>Kelvin</Unit>
When the gdal command is used, the value converted to Kelvin [K] with the following command is stored.
The slope can be left as "2.000000e-002".
>gdal_calc -A input.tif --outfile=output.tif --type=Float32 --calc="A * 0.02" --NoDataValue=65535
8) Checking the data quality flag (QA Flag)
Product quality flags (QA_Flag) are stored as 16- or 8-bit ON (1) and OFF (0) information.
Product quality information can be found in the "QA_flag information" section of each product page on the "Standard Products & Algorithms" (https://suzaku.eorc.jaxa.jp/GCOM_C/data/product_std.html) website.
| Bit No. | 15 | 14 | 13 | 12 | 11 | 10 | 9 | 8 | 7 | 6 | 5 | 4 | 3 | 2 | 1 | 0 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Example | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
9）GeoTIFF Tag Information
| Classification | Tag ID | Name | Definition/Value | |
|---|---|---|---|---|
| TIFF Field | 258 |
BitPerSample |
Number of bits per component. |
8 or 16 32（when converting physical quantity） |
259 |
Compression |
Compression method |
1: No compression 5: LZW（when compression is specified） |
|
257 |
ImageLength |
Image length |
Height (vertical length) of an image expressed in line units. |
|
256 |
ImageWidth |
Image width |
Width (horizontal length) of an image expressed in pixel units. |
|
274 |
Orientation |
Scanning orientation |
1: TopLeft（The image is stored as it is shown） |
|
262 |
PhotometricInterpretation |
Photometric type |
1: BlackIsZero（Display 0 in grayscale image as black） |
|
284 |
PlanarConfiguration |
Order of storing image data |
1: Pixel priority mode (pixel-by-pixel) e.g.RGBRGBRGB…… 2: Plane priority mode (plane-by-plane) e.g.RRR……GGG……BBB…… |
|
339 |
SampleFormat |
Type of data |
1: Unsigned integer data 3: single precision floating point number (when converting physical quantity) |
|
277 |
SamplesPerPixel |
Number of samples per pixel |
Number of stored channels/band |
|
278 |
RowsPerStrip |
The number of lines per strip |
The number of lines of image data per strip. |
|
42112 |
GDAL_METADATA |
GDAL_METADATA |
The slope and offset values are stored according to the settings. |
|
42113 |
GDAL_NODATA |
GDAL_NODATA |
Used in the GDAL library and treated as transparent. |
|
| GeoKey | 34735 |
GeoKeyDirectoryTag |
GeoKey Directory |
Store GeoKey directory |
34736 |
GeoDoubleParamsTag |
GeoDouble parameter |
Stores double-valued GeoKeys |
|
34737 |
GeoAsciiParamsTag |
GeoAscii Parameters |
Stores ASCII-valued GeoKeys |
|
33922 |
ModelTiepointTag |
Tie point |
Tie point information(I,J,K,X,Y,Z). Image position (I,J,K) and geographic coordinates (X,Y,Z) |
|
33550 |
ModelPixelScaleTag |
Pixel scale in map linear units per pixel |
Horizontal and vertical pixel spacing (ScaleX, ScaleY, ScaleZ) |
|
1024 |
GTModelTypeGeoKey |
Type of geographic coordinate information |
1: ModelTypeProjected（Projected coordinate system） 2: ModelTypeGeographic（Latitude-longitude coordinate system） |
|
1025 |
GTRasterTypeGeoKey |
Type of raster |
1: RasterPixelIsArea (One pixel represents an area in the real world) |
|
1026 |
GTCitationGeoKey |
CRS Citation |
Geographic (Lat/Lon) / WGS84 PS / WGS84 (When PS is selected) |
|
2054 |
GeogAngularUnitsGeoKey |
Coordinate unit(angle) identification code. |
9102: Degree |
|
2048 |
GeographicTypeGeoKey |
Type of coordinate system |
4326: WGS84 |
|
2049 |
GeogCitationGeoKey |
Geographic coordinate system citation |
WGS 84 (When converted by gdal) |
|
3072 |
ProjectedCSTypeGeoKey |
Projected coordinate reference system |
32767: User-defined (When PS is selected) |
|
3073 |
PCSCitationGeoKey |
Projection citation |
Geographic (Lat/Lon) / WGS84 PS / WGS84 (When PS is selected) |
|
3074 |
ProjectionGeoKey |
Projection method |
32767: User-defined (When PS is selected) |
|
3075 |
ProjCoordTransGeoKey |
Coordinate transformation method |
15: CT_PolarStereographic (When PS is selected) |
|
3076 |
ProjLinearUnitsGeoKey |
Linear unit |
9001: Linear_Meter (When PS is selected) |
|
3081 |
ProjNatOriginLatGeoKey |
Latitude of map-projection Natural origin |
71(Northern hemisphere), -71(Southern hemisphere) (When PS is selected) |
|
3082 |
ProjFalseEastingGeoKey |
Gives the easting coordinate of the map projection Natural origin |
0 (When PS is selected) |
|
3083 |
ProjFalseNorthingGeoKey |
Gives the northing coordinate of the map projection Natural origin |
0 (When PS is selected) |
|
3095 |
ProjStraightVertPoleLongGeoKey |
Longitude at Straight Vertical Pole. For polar stereographic |
0 (When PS is selected) |
|