From 622efc6414b1e9efb45a230c5f76d10bcce25a9f Mon Sep 17 00:00:00 2001
From: Giacomo Sommani <79135283+G-Sommani@users.noreply.github.com>
Date: Wed, 3 Dec 2025 15:10:07 +0100
Subject: [PATCH 1/3] Sort grid values and related arrays for second time

Sort grid values and corresponding arrays in descending order after application of Gaussian smoothing and interpolation from the equatorial map.
---
 skyreader/utils/handle_map_data.py | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/skyreader/utils/handle_map_data.py b/skyreader/utils/handle_map_data.py
index e1375ddf..da3bd44e 100644
--- a/skyreader/utils/handle_map_data.py
+++ b/skyreader/utils/handle_map_data.py
@@ -138,6 +138,12 @@ def extract_map(
         )
         grid_value = grid_value.clip(min_map, None)
 
+        sorting_indices = np.argsort(-grid_value)
+        grid_value = grid_value[sorting_indices]
+        grid_dec = grid_dec[sorting_indices]
+        grid_ra = grid_ra[sorting_indices]
+        uniq_array = uniq_array[sorting_indices]
+
     return grid_value, grid_ra, grid_dec, equatorial_map, uniq_array
 
 

From e7eedc1247cad4ea056e0649d61d31035b0ef6b3 Mon Sep 17 00:00:00 2001
From: Giacomo Sommani <79135283+G-Sommani@users.noreply.github.com>
Date: Wed, 3 Dec 2025 15:14:20 +0100
Subject: [PATCH 2/3] Avoid repetitions in sorting of grid values and
 coordinates

---
 skyreader/utils/handle_map_data.py | 19 +++++++------------
 1 file changed, 7 insertions(+), 12 deletions(-)

diff --git a/skyreader/utils/handle_map_data.py b/skyreader/utils/handle_map_data.py
index da3bd44e..1992ee25 100644
--- a/skyreader/utils/handle_map_data.py
+++ b/skyreader/utils/handle_map_data.py
@@ -90,13 +90,7 @@ def extract_map(
     grid_value: np.ndarray = np.asarray(grid_value_list)
     uniq_array: np.ndarray = np.asarray(uniq_list)
 
-    sorting_indices = np.argsort(grid_value)
-    grid_value = grid_value[sorting_indices]
-    grid_dec = grid_dec[sorting_indices]
-    grid_ra = grid_ra[sorting_indices]
-    uniq_array = uniq_array[sorting_indices]
-
-    min_value = grid_value[0]
+    min_value = np.min(grid_value)
 
     if remove_min_val or (not llh_map):
         # renormalize
@@ -111,6 +105,7 @@ def extract_map(
         # show 2 * delta_LLH
         grid_value = grid_value * 2.
         equatorial_map *= 2.
+        sorting_indices = np.argsort(grid_value)
     else:
         # Convert to probability
         equatorial_map = np.exp(-1. * equatorial_map)
@@ -137,12 +132,12 @@ def extract_map(
             equatorial_map, np.pi/2 - grid_dec, grid_ra
         )
         grid_value = grid_value.clip(min_map, None)
-
         sorting_indices = np.argsort(-grid_value)
-        grid_value = grid_value[sorting_indices]
-        grid_dec = grid_dec[sorting_indices]
-        grid_ra = grid_ra[sorting_indices]
-        uniq_array = uniq_array[sorting_indices]
+
+    grid_value = grid_value[sorting_indices]
+    grid_dec = grid_dec[sorting_indices]
+    grid_ra = grid_ra[sorting_indices]
+    uniq_array = uniq_array[sorting_indices]
 
     return grid_value, grid_ra, grid_dec, equatorial_map, uniq_array
 

From 4f6bb19d0af261f22e64dc3799bfc3aaeafb803d Mon Sep 17 00:00:00 2001
From: Giacomo Sommani <79135283+G-Sommani@users.noreply.github.com>
Date: Wed, 3 Dec 2025 15:42:07 +0100
Subject: [PATCH 3/3] Use np.nanmin for minimum value calculation

Change min function to nanmin to handle NaN values.
---
 skyreader/utils/handle_map_data.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/skyreader/utils/handle_map_data.py b/skyreader/utils/handle_map_data.py
index 1992ee25..f6cee0a4 100644
--- a/skyreader/utils/handle_map_data.py
+++ b/skyreader/utils/handle_map_data.py
@@ -90,7 +90,7 @@ def extract_map(
     grid_value: np.ndarray = np.asarray(grid_value_list)
     uniq_array: np.ndarray = np.asarray(uniq_list)
 
-    min_value = np.min(grid_value)
+    min_value = np.nanmin(grid_value)
 
     if remove_min_val or (not llh_map):
         # renormalize