pybob.ddem_tools

pybob.ddem_tools.area_alt_dist(DEM, glacier_shapes, glacier_inds=None, bin_width=None)[source]

Calculate an Area-Altitude Distribution for a glacier outline(s), given an input DEM.

Parameters:
  • DEM (pybob.GeoImg) – input DEM.
  • glacier_shapes (array-like) – mask representing glacier outlines. Can be boolean or integer depending on whether one or several AADs should be calculated.
  • glacier_inds (array-like) – array representing glacier indices in glacier_shapes. If unspecified, only one AAD is returned.
  • bin_width (float) – width of elevation bands to calculate area distribution in. If unspecified, result will be the minimum of 50m or 10% of the elevation range.
Returns bins, aads:
 

array, or list of arrays, representing elevation bands and glacier area (in DEM horizontal units) per elevation band.