A BoundingBox requested through zampy must contain at least 2x2 pixels

[retsios]

SpatialBounds(52.3, 7.1, 52.2, 6.9) is too small for zampy to handle. One will eventually get the error "ValueError: arange: cannot compute length". This is caused by the way infer_resolution()

zampy/src/zampy/utils/regrid.py

Line 21 in 7ad91da

def infer_resolution(dataset: xr.Dataset) -> tuple[float, float]:

in regrid.py works. It is needed to provide a boundingbox that includes at least 2 pixels, both horizontally and vertically. E.g. SpatialBounds(53, 7, 52, 6) would work, since dataset CAMS has 0.75 degree pixels.

[BSchilperoort]

Hi Bas,

Thanks for raising this issue! I assume the resulting dataset only has a single point of data (in space), in which case we should fill the new grid with the single available datapoint.
The 2x1 pixels case will be a little bit less trivial to properly resolve 🤔

I still have to move over the regridding routines in Zampy to https://github.com/EXCITED-CO2/xarray-regrid
Before I do so I will add these edgecases to xarray-regrid's unit tests and write a solution.