natcap/pygeoprocessing

Logging for `zonal_statistics` slightly incorrect

Closed this issue · 0 comments

I have a case where I'm running pygeoprocessing.zonal_statistics with:

  • A stack of pre-aligned rasters
  • Non-overlapping geometries, indicated to zonal_stats through polygons_might_overlap=False

I can confirm that there's only one rasterization happening, which is great! But there's also some logging that makes it seem like there's going to be more processing involved:

02/01/2024 19:55:41  (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3929) INFO rasterizing disjoint polygon set 1 of 40 set 1.0% complete ()

It turns out that there's never a set 2 of 40, so I think the logging is just using the wrong variable name.