From 242dad8a50ec5839a3ee7c82ec55e7e464a210dc Mon Sep 17 00:00:00 2001
From: Igor Tatarnikov <61896994+IgorTatarnikov@users.noreply.github.com>
Date: Fri, 16 Aug 2024 09:42:30 +0100
Subject: [PATCH] Update hemisphere slicing (#375)

* Changed the way the midpoint is calculated from center of mass to center of bounding box

* Added check for symmetrical atlas
---
 brainrender/scene.py | 19 ++++++++++---------
 1 file changed, 10 insertions(+), 9 deletions(-)

diff --git a/brainrender/scene.py b/brainrender/scene.py
index dcf4a0a2..b6deb655 100644
--- a/brainrender/scene.py
+++ b/brainrender/scene.py
@@ -287,16 +287,17 @@ def add_brain_region(
         actors = self.add(*regions)
 
         # slice to keep only one hemisphere
-        if hemisphere == "right":
-            plane = self.atlas.get_plane(
-                pos=self.root._mesh.center_of_mass(), norm=(0, 0, 1)
-            )
-        elif hemisphere == "left":
-            plane = self.atlas.get_plane(
-                pos=self.root._mesh.center_of_mass(), norm=(0, 0, -1)
-            )
-
         if hemisphere in ("left", "right"):
+            if self.atlas.metadata["symmetric"]:
+                mesh_center = (
+                    self.root._mesh.bounds().reshape((3, 2)).mean(axis=1)
+                )
+            else:
+                mesh_center = self.root._mesh.center_of_mass()
+
+            normal = (0, 0, 1) if hemisphere == "right" else (0, 0, -1)
+            plane = self.atlas.get_plane(pos=mesh_center, norm=normal)
+
             if not isinstance(actors, list):
                 actors._mesh.cut_with_plane(
                     origin=plane.center,