From a482ffd6948730dabc7eed1d3110bc3a43fb350c Mon Sep 17 00:00:00 2001 From: G-Sommani Date: Tue, 10 Dec 2024 15:40:05 +0100 Subject: [PATCH] remove nans in llh map only when calculating contours. Otherwise, nans are kept --- skyreader/plot/plot.py | 10 +++------- skyreader/utils/handle_map_data.py | 5 ----- 2 files changed, 3 insertions(+), 12 deletions(-) diff --git a/skyreader/plot/plot.py b/skyreader/plot/plot.py index d0b5520f..5fbbe227 100644 --- a/skyreader/plot/plot.py +++ b/skyreader/plot/plot.py @@ -449,16 +449,12 @@ def bounding_box(ra, dec, theta, phi): if not circular: if llh_map: - # get rid of nan values only for the contours and - # avoiding crashes during plotting - max_map = np.nanmax(equatorial_map) + # get rid of nan values only for the contours + # this avoids crashes during plotting 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 - #) + ] = np.nanmax(equatorial_map) contour_levels_for_contours = contour_levels else: grid_values_for_contours = np.log(grid_value) diff --git a/skyreader/utils/handle_map_data.py b/skyreader/utils/handle_map_data.py index 4cef11cc..661c4f14 100644 --- a/skyreader/utils/handle_map_data.py +++ b/skyreader/utils/handle_map_data.py @@ -111,11 +111,6 @@ def extract_map( # show 2 * delta_LLH grid_value = grid_value * 2. equatorial_map *= 2. - - #max_map = np.nanmax(equatorial_map) - #equatorial_map[np.isnan(equatorial_map)] = max_map - #grid_value[np.isnan(grid_value)] = max_map - #grid_value = grid_value.clip(None, max_map) else: # Convert to probability equatorial_map = np.exp(-1. * equatorial_map)