Conversation
|
Thanks @jacobbieker for this. Have you tested this against EUMETSAT? |
|
@peterdudfield Yeah, it works. Here is the output from the first two timesteps of running the consumer on the 1000m high res MTG data: |
| time = pd.to_datetime(rounded_time) | ||
| da = da.assign_coords({"time": time}).expand_dims("time") | ||
|
|
||
| da.attrs.pop("_FillValue", None) |
There was a problem hiding this comment.
im not sure this is right? I saw some errors like this with the upgrade from zarr2 to zarr3. What error do you get if you dont include this?
There was a problem hiding this comment.
If this isn't here, it gives an error of
16:58:40.876 ERROR [main.py:119] | Error: failed to prevent overwriting existing key _FillValue in attrs. This is probably an encoding field used by xarray to describe how a variable is serialized. To proceed, remove this key from the variable's attributes manually.
Removing the FillValue here fixes that and allows for the dataarrays to be added to the Icechunk store successfully
There was a problem hiding this comment.
a funny, does that happen for you on the other satellites?
What python version you running? What zarr, xarray and icechunk version you running?
There was a problem hiding this comment.
Its with Zarr 3.1.1, Xarray 2025.7.1, icechunk 1.1.1
It doesn't with the other ones, but the MTG format is quite a bit different than the MSG ones. In #44, this also does show up for other geostationary satellites as well.
|
I wonder if this can actually be achieved with the existing odegree-12 config, and it's just a matter of changing the RESOLUTION option to 500 when starting the download. I appreciate the coordinates would then be wrong, but I think with icechunk those aren't used anymore anyway... |
Pull Request
Description
Adds the MTG High resolution config, with the correct product ID to get the 4 high res channels.
How Has This Been Tested?
Running the configuration.
Checklist: