Help about Netcdf Data

[markjinz]

Hello I am encountering issues with transforming HRIT segment data into NetCDF format. Specifically, I need to create NetCDF data that adheres to WMO standards or Grib2 format using the Satpy library. Although I have attempted to generate this data using the EPTC library, I am still facing problems with the time variable and other variables.

Could someone provide a script or guidance on how to successfully create NetCDF data from HRIT segment data using Satpy or another reliable method? Any tips or advice on resolving issues related to time and other variables would be greatly appreciated.

[mraspaud]

Hi @markjinz

Something like this can work to generate a netcdf file:

    scn = Scene(filenames=filenames,
                reader='seviri_l1b_hrit',
                )

    composites = ["VIS006", "VIS008", ...]
    scn.load(composites)
    rscn = scn.resample("eurol")
    rscn.save_datasets(writer="cf")

Or have you tried that already?

[markjinz]

Hi @markjinz

Something like this can work to generate a netcdf file:

    scn = Scene(filenames=filenames,
                reader='seviri_l1b_hrit',
                )

    composites = ["VIS006", "VIS008", ...]
    scn.load(composites)
    rscn = scn.resample("eurol")
    rscn.save_datasets(writer="cf")

Or have you tried that already?

hello @mraspaud this is my code : from satpy import Scene import glob filenames = glob.glob('./raw02/H202405281015') scn = Scene(filenames=filenames, reader='seviri_l1b_hrit') scn.load(['IR_120', 'IR_108']) scn.save_datasets(writer='cf', datasets=['IR_120', 'IR_108'], filename='seviri_test.nc', exclude_attrs=['raw_metadata'],groups={'ir': ['IR_120', 'IR108']}) but i have this problem : **/home/mark/anaconda3/envs/finland/lib/python3.10/site-packages/satpy/writers/cfwriter.py:571: FutureWarning: The default behaviour of the CF writer will soon change to not compress data by default. warnings.warn("The default behaviour of the CF writer will soon change to not compress data by default.", No time dimension in datasets of group ir, skipping time bounds creation. [None, None]** And i have a heavy file as output , the main idea that i want convert my HRIT segments (8 files each channel ) to 1 file Netcdf or Grib but at the result i have a file with many times in seconds