fix issue 44

mnp/evaluation/subselection_output.py

@@ -429,10 +429,18 @@ class HotSpotsMap(SubselectionMap):
             Map with key population count per pixel
         """
 
-        self.geomap.array = sparse.csr_array(
+        hotspot_array = (
             self.key_population_count_map.geomap.array
             / self.assigned_species_map.geomap.array
         )
+
+        # Note (assigned_species > 0) / (o keypopulations) = np.nan.
+        # Replace with float32.max, so parsed to NoData later on.
+        hotspot_array[np.isnan(hotspot_array)] = np.finfo(
+            rio.dtypes.get_minimum_dtype([0.0, 1.0])
+        ).max
+
+        self.geomap.array = sparse.csr_array(hotspot_array)
         return self
 
     def save_qml(self):