--- a/trafficintelligence/indicators.py	Wed May 22 17:11:40 2024 -0400
+++ b/trafficintelligence/indicators.py	Fri May 24 16:15:38 2024 -0400
@@ -185,12 +185,19 @@
         else:
             return None
 
-    def getInstantOfMostSevereValue(self):
-        '''Returns the instant at which the indicator reaches its most severe value'''
-        if self.mostSevereIsMax:
-            return max(self.values, key=self.values.get)
+    def getInstantOfMostSevereValue(self, minSevereValue = None):
+        '''Returns the instant at which the indicator reaches its most severe value
+        or the instants when value is above minSevereValue (it not None)'''
+        if minSevereValue is None:
+            if self.mostSevereIsMax:
+                return max(self.values, key=self.values.get)
+            else:
+                return min(self.values, key=self.values.get)
         else:
-            return min(self.values, key=self.values.get)
+            if self.mostSevereIsMax:
+                return [t for t in self.values if self.values[t] >= minSevereValue]
+            else:
+                return [t for t in self.values if self.values[t] <= minSevereValue]
 
 # functions to aggregate discretized maps of indicators
 # TODO add values in the cells between the positions (similar to discretizing vector graphics to bitmap)