Mathias Hauser

Take a look at https://github.com/mathause/filefinder - however, I am not sure it matches your specs & under the hood it also uses `glob` but it also parses the patterns such...

Indeed, there is a a bit of ocean not part of any region - but it's the same with both versions of the naturalearthdata. v4.1 is orange and v5.0 is...

> @mathause have you interacted with that community at all? No I didn't so far, but you have ;-) nvkelso/natural-earth-vector#319

I was able to simplify the problem a bit (no transform, simpler coordinates): ``` python from shapely import geometry from rasterio import features import numpy as np outline = np.array([...

I was not able to construct an inconsistent example with `all_touched=True`, also when shifting the `Polygon` by 0.5. I can try to reproduce directly in `gdal`.

xarray (I think) and regionmask would likely welcome this change. A proper deprecation could be a bit of a pain. You could check if `get_logger().level != default` and issue a...

Just to support the point - [xarray sets the pooch log-level](https://github.com/pydata/xarray/blob/d994273b5d61bbdc2fde5050af5922deb66db274/xarray/tutorial.py#L121) when tutorial data is downloaded: ```python logger = pooch.get_logger() logger.setLevel("WARNING") ``` Setting the log-level for other users was certainly...

> Would be easy if we could just force named arguments? Yes that would be a good idea. A deprecation would be nice but for the backends refactor (`open_dataset`) this...

Not sure if this is too naïve from my side but can you just pass N+1 levels to the function?

Thanks for your report. I think the issue is that dask cannot correctly infer the dtype of the result - or at least not it's length (maybe because it does...