Help needed! I am trying to convert the full disk HSD (Himawari Standard Data) file from Himawari-8 to a WGS84 TIFF. When using Satpy Scene, it only outputs a disk-shaped TIFF with an unknown geographic coordinate system. I need assistance from expe

[ll ll]

mport sys
from datetime import datetime
from datetime import timedelta
import glob
from satpy.scene import Scene
import numpy as np


scen = Scene(filenames=[r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0110.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0210.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0310.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0410.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0510.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0610.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0710.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0810.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S0910.DAT',
                        r'**\HS_H08_20210309_0500_B01_FLDK_R10_S1010.DAT',
                        ],
                        reader='ahi_hsd',
                        reader_kwargs={'mask_space': False})#reader_kwargs={'calib_mode': 'update'}
scen.load(['B01'])
scen.load(scen.available_dataset_ids())


# band_data = scen['B01'].compute()
# band_data_array = band_data.values
scen.save_datasets(writer='geotiff', dtype= np.float32, enhance= False, base_dir=r"**")