Add max DSM metric

city_metrix/layers/layer.py

@@ -65,6 +65,9 @@ def mean(self):
 
     def count(self):
         return self._zonal_stats("count")
+
+    def max(self):
+        return self._zonal_stats("max")
 
     def _zonal_stats(self, stats_func):
         if box(*self.zones.total_bounds).area <= MAX_TILE_SIZE**2:

city_metrix/metrics/__init__.py

@@ -3,4 +3,5 @@
 from .built_land_with_high_land_surface_temperature import built_land_with_high_land_surface_temperature
 from .mean_tree_cover import mean_tree_cover
 from .urban_open_space import urban_open_space
-from .natural_areas import natural_areas
+from .natural_areas import natural_areas
+from .max_dsm import max_dsm

city_metrix/metrics/max_dsm.py

@@ -0,0 +1,14 @@
+from geopandas import GeoDataFrame, GeoSeries
+
+from city_metrix.layers import AlosDSM
+from exactextract import exact_extract
+
+def max_dsm(zones: GeoDataFrame) -> GeoSeries:
+    """
+    Get max AlosDSM.
+    :param zones: GeoDataFrame with geometries to collect zonal stats on
+    :return: Pandas Series of hights
+    """
+    max_height = exact_extract(AlosDSM(), zones, ["max"])
+
+    return max_height

environment.yml

@@ -21,3 +21,4 @@ dependencies:
   - pip=23.3.1
   - pip:
     - cartoframes==1.2.5
+    - git+https://github.com/isciences/exactextract