Dan Baston

I'm looking into how we might accomplish this using OverlayNG, using a simple case of a square split with a line. It is straightfoward to mark a set of OverlayEdges...

I had been thinking of having `GDALNoDataMaskBand` read the requested window into a memory dataset at full resolution, reclassify the values to 0 or 255, and then read from the...

I'm having trouble reproducing this. What platform are you on? I notice that there is a typo in the warning message.

> the code before https://github.com/libgeos/geos/commit/60edf0e5f040e065193bd5beb79f7fa8febb1cd5 did not initialize the coordinates We didn't initialize them explicitly, but they were still initialized. If I create a `std::vector(5)` or `std::array`, the five elements...

It works for me on the netCDF files I have with multiple variables. Do you have a file you can share so I can reproduce the error? In the not-too-distant...

This works for me: ``` exactextract -r "d2m:NETCDF:/home/dan/Downloads/era5_d2m_modified.nc:d2m[1]" -p /home/dan/data/gadm36_level_0.gpkg -f GID_0 -s "mean(d2m)" -o /tmp/test.csv ``` ``` GID_0,d2m_mean ABW,-nan AFG,4685.02294921875 AGO,26304.724609375 AIA,-nan ALA,3962.07885742188 ... ``` This is the latest...

> But I am getting a warning about a pybind11 path inconsistency upon running the cmake commmand when installing exact extract. It doesn't sound related, but please feel free to...

With https://github.com/isciences/exactextract/pull/59 you can do, e.g. ``` from exactextract import exact_extract import fiona import rioxarray import xarray xds = xarray.open_dataset('/tmp/era5.nc') polys = fiona.open('/home/dan/data/ne_50m_admin_0_countries.shp') results = exact_extract(xds['d2m'], polys, 'mean') ``` this...

Any thoughts on an API that would work well for shapely?

I didn't add a progress callback because > 90% of the time was inside GEOS. I've thought about making the GDAL GEOS conversions more efficient (WKB is a slow way...