HOME EN
Welcome to the SGLI_L1 wiki!
This is a program that converts GCOM-C (SHIKISAI) SGLI Level 1 HDF5 files to GeoTIFF format and performs map projection to ESPG4326. It should work in an environment with Python 3 and GDAL installed.
python L1.py sgli_l1.h5 output_color.tif
If you want to extract other bands, you can specify them as follows:
# Extracting only RGB bands
sgli_bands=['Lt_VN03','Lt_VN05','Lt_VN08']
You can choose bands from Lt_VN01 to Lt_VN11.
We will use OpenStreetMap for position verification. The L1 HDF5 file contains location information, which can be used to minimize the displacement of the coastline.
The program takes GCP (Ground Control Points) every 200 pixels, but you can also choose to place them every 100 pixels. This allows for handling more GCPs compared to using the gdalwarp command. However, it will increase processing time.
Here is an example of the coastline of Madagascar's western region:
The filled area represents the coastline of Madagascar according to OpenStreetMap.
I have added the multithread option to gdalwarp, which improves processing speed. However, please note that the effectiveness may be limited if you have fewer cores available.
By utilizing the multithread option, you can take advantage of parallel processing to speed up the conversion process. This can be particularly beneficial when working with large datasets or running the program on a system with multiple CPU cores.
Please ensure that your system and GDAL version support the multithread option and adjust the settings accordingly to optimize the performance.