Closed danielFlemstrom closed 4 months ago
data.xvec.zonal_stats assumes that geometries and data has same CRS.
As a consequence we need to update aggregate.py::aggregate_spatial with
... gdf = gdf.to_crs(data.rio.crs) <------ This one needs to be added geometries = gdf.geometry.values positional_parameters = {"data": 0} vec_cube = data.xvec.zonal_stats( ...
Example when we get all null out of aggregate spatial
IN: data with SWEREF geometries in WGS84 reducer median OUT: polygon misses the data and we get all null
data.xvec.zonal_stats assumes that geometries and data has same CRS.
As a consequence we need to update aggregate.py::aggregate_spatial with
Example when we get all null out of aggregate spatial
IN: data with SWEREF geometries in WGS84 reducer median OUT: polygon misses the data and we get all null