When using clisops.core.average_shape, I'm getting GEOSException: TopologyException: unable to assign free hole to a shell with a multipolygon. This error is raised by shape_bbox_indexerunary_union operation.
The logic in shape_bbox_indexer is only meant to make a rough subset of the area of interest to reduce data transfers. We can probably compute a simple convex hull to get the rough shape of the bounding box.
huard
commented
1 year ago
Description
When using clisops.core.average_shape, I'm getting
GEOSException: TopologyException: unable to assign free hole to a shellwith a multipolygon. This error is raised byshape_bbox_indexerunary_unionoperation.The logic in
shape_bbox_indexeris only meant to make a rough subset of the area of interest to reduce data transfers. We can probably compute a simple convex hull to get the rough shape of the bounding box.