Skip to content

map_raster yield incorrect ancillary projection/interpolation with dask=2025.9.1 #261

Description

@Skealz

For another wind estimation processor that use similar code as xsar, I noticed an issue with ancillary data in the resulting product. I investigated and found that map_raster would interpolate ancillary data incorrectly.
The problematic section of the map_raster is :
https://github.com/umr-lops/xsar/blob/develop/src/xsar/base_dataset.py#L922C1-L933C18
I could confirm that this section yield incorrect results by forcing to go through the first section of the "if" block, and observed radically different results.

Going through the first part of the "if" block (xarray interp) - OK :
Image

Going through the second part (RectBivariateSpline) - NOT OK :
Image

I investigated further and finally tested using an older version of dask package.
Using dask=2025.3.0 fix the issue without modifying the code.

Going through the second part (RectBivariateSpline) with dask=2025.3.0 - OK :
Image

Short term fix :
Pin to dask 2025.3.0
EDIT : dask 2025.7.0 works fine too.

Long term fix :
For each xsar releases, pin important dependencies, at carefully tested versions

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions