From aacb0615f6be6ffc3bfba9d5f3d49e8326c4df0c Mon Sep 17 00:00:00 2001
From: Ahmet Nihat Simsek <ahmet.n.simsek@gmail.com>
Date: Mon, 23 Sep 2024 15:51:25 +0200
Subject: [PATCH] fix: statisitcal map centroid calculation resulting in NaN
 points

---
 siibra/volumes/parcellationmap.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/siibra/volumes/parcellationmap.py b/siibra/volumes/parcellationmap.py
index c73ec7cb..68e49024 100644
--- a/siibra/volumes/parcellationmap.py
+++ b/siibra/volumes/parcellationmap.py
@@ -593,7 +593,7 @@ def compute_centroids(self) -> Dict[str, point.Point]:
             with QUIET:
                 mapimg = self.fetch(index=index)  # returns a mask of the region
             maparr = np.asanyarray(mapimg.dataobj)
-            centroid_vox = np.mean(np.where(maparr == 1), axis=1)
+            centroid_vox = np.mean(np.nonzero(maparr), axis=1)
             assert regionname not in centroids
             centroids[regionname] = point.Point(
                 np.dot(mapimg.affine, np.r_[centroid_vox, 1])[:3], space=self.space