The changes made in natcap/pygeoprocessing#428 caused an HRA test to fail, because HRA's geometry simplification function was always setting the vector geometry type to Unknown. pygeoprocessing.zonal_statistics will now require a type of Polygon or MultiPolygon, so this raised an error.
I think we can copy over the geometry type from the input vector, which should always be Polygon or MultiPolygon for the AOI. It may still be Unknown in the case of habitat or stressor vectors that contain multiple geometry types.
The changes made in natcap/pygeoprocessing#428 caused an HRA test to fail, because HRA's geometry simplification function was always setting the vector geometry type to
Unknown.pygeoprocessing.zonal_statisticswill now require a type ofPolygonorMultiPolygon, so this raised an error.I think we can copy over the geometry type from the input vector, which should always be
PolygonorMultiPolygonfor the AOI. It may still beUnknownin the case of habitat or stressor vectors that contain multiple geometry types.