pyflwdir.regions.region_slices#

pyflwdir.regions.region_slices(regions: ndarray) → tuple[ndarray, list[tuple[slice, ...]]][source]#

Return array slices bounding each unique label in regions.

NOTE: a region must be a connected area with the same ID, where IDs are integer values greater than zero.

Parameters:

regions (2D array of int) – Raster of positive region IDs.

Returns:

  • lbs (1D array) – array of the unique region IDs

  • slices (list of tuples) – Each tuple contains slices, one for each dimension