2.4.2

• Fixed an issue where MFD flow direction was producing many nodata holes given
  a large-enough DEM. These nodata holes would then propagate to flow
  accumulation and stream extraction, producing very disjointed stream
  networks. #350
• Improved progress logging in MFD flow direction, MFD flow accumulation, MFD
  stream extraction to prevent integer overflows in percentages and improve
  the readability of log messages. #246
• Fixed an issue introduced in 2.4.1 where zonal_statistics results could
  vary depending on the bounding boxes of the disjoint polygon sets calculated
  from the input vector.
• Added a new function, pygeoprocessing.align_bbox, which pads a bounding
  box to align with the grid of a given geotransform.
• Added support for int64 and uint64 (GDAL 3.5+)
  #352
• Added support for signed bytes (GDAL 3.7+)
  #329
• Fixed a bug where calculate_slope did not work as expected when the DEM
  raster's nodata value is NaN (#352)