in flattened probability map, fill nans with the minimum between map …

…minimum or 1e-16
icecube · Dec 10, 2024 · cbec55e · cbec55e
1 parent c42e683
commit cbec55e
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/skyreader/plot/plot.py b/skyreader/plot/plot.py
@@ -832,8 +832,9 @@ def bounding_box(ra, dec, theta, phi):
             column_names = ['2DLLH']
         else:
             # avoid excessively heavy data format for the flattened map
-            equatorial_map[np.isnan(equatorial_map)] = np.nanmin(
-                equatorial_map
+            min_prob = np.nanmin(equatorial_map)
+            equatorial_map[np.isnan(equatorial_map)] = np.min(
+                [min_prob, 1e-16]
             )
             equatorial_map[equatorial_map < 1e-16] = np.mean(
                 equatorial_map[equatorial_map < 1e-16]