constraints on values of grid_value

icecube · Sep 13, 2024 · a9a5dfb · a9a5dfb
1 parent 2ffef7c
commit a9a5dfb
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/skyreader/plot/plot.py b/skyreader/plot/plot.py
@@ -552,6 +552,8 @@ def bounding_box(ra, dec, theta, phi):
         grid_value = healpy.get_interp_val(
             equatorial_map, np.pi/2 - grid_dec, grid_ra
         )
+        grid_value[np.isnan(grid_value)]=0.0
+        grid_value = grid_value.clip(0.0, None)
         sorted_values = np.sort(equatorial_map)[::-1]
 
         # Calculate the contours