remove infinite values from llh map

icecube · Dec 10, 2024 · ae6b1ed · ae6b1ed
1 parent ef5a1bf
commit ae6b1ed
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/skyreader/plot/plot.py b/skyreader/plot/plot.py
@@ -451,14 +451,14 @@ def bounding_box(ra, dec, theta, phi):
             if llh_map:
                 # get rid of nan values only for the contours and
                 # avoiding crashes during plotting
-                #max_map = np.nanmax(equatorial_map)
+                max_map = np.nanmax(equatorial_map)
                 grid_values_for_contours = copy.copy(grid_value)
                 #grid_values_for_contours[
                 #    np.isnan(grid_values_for_contours)
                 #] = max_map
-                #grid_values_for_contours = grid_values_for_contours.clip(
-                #    None, max_map
-                #)
+                grid_values_for_contours = grid_values_for_contours.clip(
+                    None, max_map
+                )
                 contour_levels_for_contours = contour_levels
             else:
                 grid_values_for_contours = np.log(grid_value)