feat: Add right atrium appendage landmarks for Rodero 2021 cases (#1202)

Co-authored-by: pyansys-ci-bot <[email protected]>

Author: mhoeijm

doc/source/changelog/1202.added.md

@@ -0,0 +1 @@
+Add right atrium appendage landmarks

examples/preprocessor/compute-atrial-fibers_pr.py

@@ -104,9 +104,14 @@
 # compute left atrium fiber
 la = simulator.compute_left_atrial_fiber()
 
+# Import the appendage landmarks.
+from ansys.health.heart.pre.database_utils import right_atrium_appendage_landmarks
+
+# Get the right atrium appendage landmark of the first case of Rodero2021.
+right_atrium_appendage_coordinates = right_atrium_appendage_landmarks.get("Rodero2021").get(1)
+
 # Appendage apex point should be manually given to compute right atrium fiber
-appendage_apex = [39, 29, 98]
-ra = simulator.compute_right_atrial_fiber(appendage_apex)
+ra = simulator.compute_right_atrial_fiber(appendage=right_atrium_appendage_coordinates)
 
 ###############################################################################
 # .. note::
@@ -150,8 +155,16 @@
 mesh = ra.ctp()
 streamlines = mesh.streamlines(vectors="e_l", source_radius=50, n_points=5000)
 tubes = streamlines.tube()
-plotter.add_mesh(mesh, opacity=0.5, color="white")
-plotter.add_mesh(tubes, color="red")
+plotter.add_mesh(mesh, opacity=0.5, color="white", label="myocardium")
+plotter.add_mesh(tubes, color="red", label="fibers")
+plotter.add_mesh(
+    pv.PolyData(right_atrium_appendage_coordinates),
+    color="blue",
+    point_size=20,
+    render_points_as_spheres=True,
+    label="right atrium appendage",
+)
+plotter.add_legend()
 plotter.show()
 
 ###############################################################################

examples/simulator/ep-mechanics-simulator-fullheart.py

@@ -101,11 +101,18 @@
 ###############################################################################
 # Compute the fiber orientation
 # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# Compute the fiber orientation on the entire model.
+
+# Import the appendage landmarks.
+from ansys.health.heart.pre.database_utils import right_atrium_appendage_landmarks
+
+# Get the right atrium appendage landmark of the first case of Rodero2021.
+right_atrium_appendage_coordinates = right_atrium_appendage_landmarks.get("Rodero2021").get(1)
 
 # Compute fiber orientation in the ventricles and atria.
 simulator.compute_fibers()
 simulator.compute_left_atrial_fiber()
-simulator.compute_right_atrial_fiber(appendage=[39, 29, 98])
+simulator.compute_right_atrial_fiber(appendage=right_atrium_appendage_coordinates)
 
 # Switch the atria to active.
 simulator.model.left_atrium.fiber = True

examples/simulator/ep-simulator-fullheart_rodero.py

@@ -107,6 +107,12 @@
 # ~~~~~~~~~~~~~~~~~~~~~~~~~
 # Compute fiber orientation and plot the fibers on the entire model.
 
+# Import the appendage landmarks.
+from ansys.health.heart.pre.database_utils import right_atrium_appendage_landmarks
+
+# Get the right atrium appendage landmark of the first case of Rodero2021.
+right_atrium_appendage_coordinates = right_atrium_appendage_landmarks.get("Rodero2021").get(1)
+
 # Compute ventricular fibers.
 simulator.compute_fibers(method="D-RBM")
 
@@ -116,7 +122,7 @@
 simulator.model.right_atrium.fiber = True
 simulator.model.left_atrium.fiber = True
 simulator.compute_left_atrial_fiber()
-simulator.compute_right_atrial_fiber(appendage=[39, 29, 98])
+simulator.compute_right_atrial_fiber(appendage=right_atrium_appendage_coordinates)
 simulator.model.plot_fibers(n_seed_points=1000)
 
 ###############################################################################

src/ansys/health/heart/pre/database_utils.py

@@ -627,3 +627,30 @@ def get_compatible_input(
     "Superior vena cava inlet": 16,
     "Inferior vena cava inlet": 17,
 }
+
+# Rodero2021 left atrial appendage landmarks. These landmarks are used as input
+# to compute the fibers on the left atrium.
+right_atrium_appendage_landmarks = {
+    "Rodero2021": {
+        1: [39, 29, 98],
+        2: [37, 41, 103],
+        3: [39, 60, 99],
+        4: [33, 34, 102],
+        5: [25, 44, 91],
+        6: [26, 36, 86],
+        7: [33, 28, 91],
+        8: [35, 39, 94],
+        9: [39, 27, 100],
+        10: [33, 36, 96],
+        11: [30, 50, 91],
+        12: [28, 42, 93],
+        13: [27, 46, 89],
+        14: [44, 23, 101],
+        15: [37, 32, 88],
+        16: [36, 39, 98],
+        17: [28, 30, 94],
+        18: [26, 30, 83],
+        19: [39, 33, 88],
+        20: [32, 38, 95],
+    }
+}