search

pytroll / satpy

Python package for earth-observing satellite data processing
http://satpy.readthedocs.org/en/latest/
GNU General Public License v3.0
1.06k stars 292 forks source link

Read the EUMETSAT Test Data for MTG Using Satpy #2809

Open FredFiro88 opened 3 months ago

FredFiro88 commented 3 months ago

Hello everyone, I was able to read the old data normally using Satpy, but I'm having issues reading the latest available data (FDHSI & HRFI ) with the fci_l1c_nc reader. Should I use different methods or what should I do ? https://sftp.eumetsat.int/public/folder/UsCVknVOOkSyCdgpMimJNQ/User-Materials/Test-Data/MTG/MTG_FCI_L1C_RC72-20240113_TD-505_Feb2024/ If anyone has managed to read this data and extract some products , I would appreciate it if they could share the script with us . Thank you in advance !

gerritholl commented 3 months ago

Can you please add the code you are running and the error message you are getting? What issues are you having?

FredFiro88 commented 3 months ago

Can you please add the code you are running and the error message you are getting? What issues are you having?

from satpy.scene import Scene from satpy.writers import to_image from satpy.resample import get_area_def from satpy import find_files_and_readers from IPython.display import display import glob import warnings import matplotlib.pyplot as plt import warnings warnings.filterwarnings('ignore')

fnames = glob.glob('./MTG/*.nc') from satpy import available_readers scn = Scene(reader='fci_l1c_nc', filenames=fnames)

this is the error when i use FDHSI data : WARNING:satpy.readers:Don't know how to open the following files: {'./MTG/W_XX-EUMETSAT-Darmstadt,IMG+SAT,MTI1+FCI-1C-RRAD-FDHSI-FD--CHK-TRAIL--DIS-NC4E_C_EUMT_20240113115255_IDPFI_OPE_20240113115007_20240113115924_N_JLS_C_0072_0041.nc'} this is the error when i use HRFI data : WARNING:satpy.readers.yaml_reader:No filenames found for reader: fci_l1c_nc

gerritholl commented 3 months ago

Your code is correct (but missing the loading, resampling, saving part). The warning message is fine; satpy doesn't use or need the trail (and indeed cannot read it), just the body files. HRFI data should work; do you have the latest satpy?

ameraner commented 3 months ago

You may also need to add import hdf5plugin in order to be able to open compressed files. See also the example here: https://satpy.readthedocs.io/en/stable/examples/fci_l1c_natural_color.html