Mathias Hauser

And thanks for the feedback b.t.w.!

Points I'd be interested in for the comparison. - polygons - one-at-a-time - several-at-once - grid - regular grid - irregular grid - unstructured grids - overlapping regions - fractional...

```python import pandas as pd attr = dict() # regionmask attr["regionmask"] = { "polygons": "many", "regular grid": "yes", "irregular grid": "yes", "unstructured grids": "yes", "overlapping regions": "yes", "point-on-line": "rasterio", "fractional...

**python-rasterstats** - https://github.com/perrygeo/python-rasterstats - wrapper around `rasterio.features.rasterize` - tailored to working directly with tif and shp files

With the release of v0.10.0 I am closing this issue. Technically this was already fixed in v0.9.0 but this has seen some internal improvements. I am still considering to always...

As reference: ```python if not as_3D: mask = list() for num in numbers: mask_coarse = mask_sampled.values == num mask.append(mask_coarse[np.newaxis, ...]) mask_sampled = np.concatenate(mask, axis=0) else: mask_sampled = mask_sampled.values ```

Simpler repro (no datetime, no resample) - but I don't know why it happens, either. ```python import xarray as xr data1 = xr.DataArray([np.nan,1], dims='x', coords={'x': [0, 6]}) data2 = data1.astype('float32')...

Ok this is a scipy problem - do you want to raise a issue in scipy? ```python import scipy as sp import numpy as np xi = np.array([0, 6]) yi...

Or #7963?

I would like to merge this but after checking #7963 and #8301 I am a bit confused. - Should I update the error message to accommodate for the "second chunk...