cardio 2025.10.1__py3-none-any.whl → 2026.1.0__py3-none-any.whl

cardio/volume.py

@@ -30,6 +30,9 @@ class Volume(Object):
     _mpr_actors: dict[str, list[vtk.vtkImageActor]] = pc.PrivateAttr(
         default_factory=dict
     )
+    _crosshair_actors: dict[str, dict[str, vtk.vtkActor]] = pc.PrivateAttr(
+        default_factory=dict
+    )
 
     @pc.model_validator(mode="after")
     def initialize_volume(self):
@@ -134,6 +137,104 @@ class Volume(Object):
             return self.create_mpr_actors(frame)
         return self._mpr_actors[frame]
 
+    def create_crosshair_actors(self, colors: dict, line_width: float = 1.5) -> dict:
+        """Create 2D crosshair overlay actors for each MPR view.
+
+        Uses screen-space 2D actors that always appear centered in the view.
+
+        Args:
+            colors: Dict mapping view names to RGB tuples
+            line_width: Width of the crosshair lines
+
+        Returns:
+            Dict mapping view names to dicts with line actors
+        """
+        crosshairs = {}
+
+        for view_name in ["axial", "sagittal", "coronal"]:
+            view_crosshairs = {}
+
+            for line_name in ["line1", "line2"]:
+                # Create 2D line using normalized viewport coordinates
+                points = vtk.vtkPoints()
+                lines = vtk.vtkCellArray()
+
+                # Placeholder points - will set based on line orientation
+                if line_name == "line1":
+                    # Vertical line (one of the other planes)
+                    points.InsertNextPoint(0.5, 0.0, 0.0)
+                    points.InsertNextPoint(0.5, 1.0, 0.0)
+                else:
+                    # Horizontal line (other plane)
+                    points.InsertNextPoint(0.0, 0.5, 0.0)
+                    points.InsertNextPoint(1.0, 0.5, 0.0)
+
+                line = vtk.vtkLine()
+                line.GetPointIds().SetId(0, 0)
+                line.GetPointIds().SetId(1, 1)
+                lines.InsertNextCell(line)
+
+                polydata = vtk.vtkPolyData()
+                polydata.SetPoints(points)
+                polydata.SetLines(lines)
+
+                # Use coordinate transform to map normalized coords to viewport
+                coord = vtk.vtkCoordinate()
+                coord.SetCoordinateSystemToNormalizedViewport()
+
+                mapper = vtk.vtkPolyDataMapper2D()
+                mapper.SetInputData(polydata)
+                mapper.SetTransformCoordinate(coord)
+
+                actor = vtk.vtkActor2D()
+                actor.SetMapper(mapper)
+                actor.GetProperty().SetLineWidth(line_width)
+                actor.SetVisibility(False)
+
+                view_crosshairs[line_name] = {
+                    "polydata": polydata,
+                    "actor": actor,
+                }
+
+            # Set colors based on which planes the lines represent
+            if view_name == "axial":
+                view_crosshairs["line1"]["actor"].GetProperty().SetColor(
+                    *colors.get("sagittal", (1, 0, 0))
+                )
+                view_crosshairs["line2"]["actor"].GetProperty().SetColor(
+                    *colors.get("coronal", (0, 1, 0))
+                )
+            elif view_name == "sagittal":
+                view_crosshairs["line1"]["actor"].GetProperty().SetColor(
+                    *colors.get("coronal", (0, 1, 0))
+                )
+                view_crosshairs["line2"]["actor"].GetProperty().SetColor(
+                    *colors.get("axial", (0, 0, 1))
+                )
+            else:  # coronal
+                view_crosshairs["line1"]["actor"].GetProperty().SetColor(
+                    *colors.get("sagittal", (1, 0, 0))
+                )
+                view_crosshairs["line2"]["actor"].GetProperty().SetColor(
+                    *colors.get("axial", (0, 0, 1))
+                )
+
+            crosshairs[view_name] = view_crosshairs
+
+        self._crosshair_actors = crosshairs
+        return crosshairs
+
+    @property
+    def crosshair_actors(self) -> dict:
+        """Get crosshair actors."""
+        return self._crosshair_actors
+
+    def set_crosshairs_visible(self, visible: bool):
+        """Set visibility of all crosshair actors."""
+        for view_crosshairs in self._crosshair_actors.values():
+            for line_data in view_crosshairs.values():
+                line_data["actor"].SetVisibility(visible)
+
     def _get_mpr_coordinate_systems(self):
         """Get coordinate system transformation matrices for MPR views."""
         view_axcodes = {
@@ -197,75 +298,104 @@ class Volume(Object):
 
         return bounds
 
+    def _build_cumulative_rotation(
+        self, rotation_sequence: list, rotation_angles: dict, angle_units=None
+    ) -> np.ndarray:
+        """Build cumulative rotation matrix from sequence of rotations."""
+        from .orientation import AngleUnits
+
+        if angle_units is None:
+            angle_units = AngleUnits.DEGREES
+
+        cumulative_rotation = np.eye(3)
+        if rotation_sequence and rotation_angles:
+            for i, rotation in enumerate(rotation_sequence):
+                angle = rotation_angles.get(i, 0)
+                rotation_matrix = euler_angle_to_rotation_matrix(
+                    EulerAxis(rotation["axis"]), angle, angle_units
+                )
+                cumulative_rotation = cumulative_rotation @ rotation_matrix
+        return cumulative_rotation
+
+    def get_scroll_vector(
+        self,
+        view_name: str,
+        rotation_sequence: list = None,
+        rotation_angles: dict = None,
+        angle_units=None,
+    ) -> np.ndarray:
+        """Get the current normal vector for a view after rotation.
+
+        Args:
+            view_name: One of "axial", "sagittal", "coronal"
+            rotation_sequence: List of rotation definitions
+            rotation_angles: Dict mapping rotation index to angle
+            angle_units: AngleUnits enum (degrees or radians)
+
+        Returns:
+            3D unit vector representing the scroll direction for this view
+        """
+        base_normals = {
+            "axial": np.array([0.0, 0.0, 1.0]),
+            "sagittal": np.array([1.0, 0.0, 0.0]),
+            "coronal": np.array([0.0, 1.0, 0.0]),
+        }
+
+        if view_name not in base_normals:
+            return np.array([0.0, 0.0, 1.0])
+
+        cumulative_rotation = self._build_cumulative_rotation(
+            rotation_sequence, rotation_angles, angle_units
+        )
+        return cumulative_rotation @ base_normals[view_name]
+
     def update_slice_positions(
         self,
         frame: int,
-        axial_pos: float,
-        sagittal_pos: float,
-        coronal_pos: float,
+        origin: list,
         rotation_sequence: list = None,
         rotation_angles: dict = None,
+        angle_units=None,
     ):
         """Update slice positions for MPR views with optional rotation.
 
         Args:
             frame: Frame index
-            axial_pos: Physical position along Z axis (LAS Superior)
-            sagittal_pos: Physical position along X axis (LAS Left)
-            coronal_pos: Physical position along Y axis (LAS Anterior)
+            origin: [x, y, z] position in LPS coordinates (shared by all views)
+            rotation_sequence: List of rotation definitions
+            rotation_angles: Dict mapping rotation index to angle
+            angle_units: AngleUnits enum (degrees or radians), defaults to DEGREES
         """
+        from .orientation import AngleUnits
+
+        if angle_units is None:
+            angle_units = AngleUnits.DEGREES
         if frame not in self._mpr_actors:
             return
 
-        volume_actor = self._actors[frame]
-        image_data = volume_actor.GetMapper().GetInput()
-        bounds = image_data.GetBounds()
-
         actors = self._mpr_actors[frame]
 
-        # Clamp positions to volume bounds
-        axial_pos = max(bounds[4], min(bounds[5], axial_pos))  # Z bounds
-        sagittal_pos = max(bounds[0], min(bounds[1], sagittal_pos))  # X bounds
-        coronal_pos = max(bounds[2], min(bounds[3], coronal_pos))  # Y bounds
-
         # Get coordinate system transformations for each MPR view
         transforms = self._get_mpr_coordinate_systems()
 
-        center = image_data.GetCenter()
-
-        # Step 1: Apply translation to determine slice origins in physical space
-        axial_origin = [center[0], center[1], axial_pos]
-        sagittal_origin = [sagittal_pos, center[1], center[2]]
-        coronal_origin = [center[0], coronal_pos, center[2]]
+        # Build cumulative rotation matrix
+        cumulative_rotation = self._build_cumulative_rotation(
+            rotation_sequence, rotation_angles, angle_units
+        )
 
-        # Step 2: Apply cumulative rotation around the translated origins
-        if rotation_sequence and rotation_angles:
-            cumulative_rotation = np.eye(3)
-            for i, rotation in enumerate(rotation_sequence):
-                angle = rotation_angles.get(i, 0)
-                rotation_matrix = euler_angle_to_rotation_matrix(
-                    EulerAxis(rotation["axis"]), angle
-                )
-                cumulative_rotation = cumulative_rotation @ rotation_matrix
+        # Apply rotation to base transforms
+        axial_transform = cumulative_rotation @ transforms["axial"]
+        sagittal_transform = cumulative_rotation @ transforms["sagittal"]
+        coronal_transform = cumulative_rotation @ transforms["coronal"]
 
-            # Apply rotation to base transforms
-            axial_transform = cumulative_rotation @ transforms["axial"]
-            sagittal_transform = cumulative_rotation @ transforms["sagittal"]
-            coronal_transform = cumulative_rotation @ transforms["coronal"]
-        else:
-            # Use base transforms without rotation
-            axial_transform = transforms["axial"]
-            sagittal_transform = transforms["sagittal"]
-            coronal_transform = transforms["coronal"]
-
-        # Update slices with translated origins and rotated transforms
-        axial_matrix = create_vtk_reslice_matrix(axial_transform, axial_origin)
+        # All views share the same origin
+        axial_matrix = create_vtk_reslice_matrix(axial_transform, origin)
         actors["axial"]["reslice"].SetResliceAxes(axial_matrix)
 
-        sagittal_matrix = create_vtk_reslice_matrix(sagittal_transform, sagittal_origin)
+        sagittal_matrix = create_vtk_reslice_matrix(sagittal_transform, origin)
         actors["sagittal"]["reslice"].SetResliceAxes(sagittal_matrix)
 
-        coronal_matrix = create_vtk_reslice_matrix(coronal_transform, coronal_origin)
+        coronal_matrix = create_vtk_reslice_matrix(coronal_transform, origin)
         actors["coronal"]["reslice"].SetResliceAxes(coronal_matrix)
 
     def update_mpr_window_level(self, frame: int, window: float, level: float):

cardio/volume_property.py

@@ -1,8 +1,6 @@
-# Third Party
 import pydantic as pc
 import vtk
 
-# Internal
 from .blend_transfer_functions import blend_transfer_functions
 from .transfer_function_pair import TransferFunctionPairConfig
 from .types import ScalarComponent

cardio/volume_property_presets.py

@@ -1,11 +1,8 @@
-# System
 import pathlib as pl
 
-# Third Party
 import pydantic as pc
 import tomlkit as tk
 
-# Internal
 from .volume_property import VolumePropertyConfig
 
 

cardio/window_level.py

@@ -1,7 +1,5 @@
-# System
 import functools as ft
 
-# Third Party
 import pydantic as pc
 
 

{cardio-2025.10.1.dist-info → cardio-2026.1.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: cardio
-Version: 2025.10.1
+Version: 2026.1.0
 Summary: A simple web-based viewer for 3D and 4D ('cine') medical imaging data.
 Keywords: Medical,Imaging,3D,4D,Visualization
 Author: Davis Marc Vigneault
@@ -59,12 +59,11 @@ Description-Content-Type: text/markdown
 `cardio` is a simple web-based viewer for 3D and 4D ('cine') medical imaging data,
 built primarily on [trame](https://github.com/kitware/trame),
 [vtk](https://github.com/kitware/vtk), and
-[itk](https://github.com/insightsoftwareconsortium/itk).  `cardio` is able to
-render sequences of mesh files (e.g., `\*.obj` files), segmentation files (e.g.
-`\*nii.gz` files with discrete labels) and volume renderings of grayscale images
-(e.g. \*.nii.gz files with continuous values).  `cardio` is launched from the
-commandline and may be configured either directly from the commandline, via a static
-TOML configuration file, or a combination of the two.
+[itk](https://github.com/insightsoftwareconsortium/itk).  `cardio` can render sequences
+of mesh files (e.g., `\*.obj` files), segmentation files (e.g., `\*nii.gz` files with
+discrete labels) and volume renderings of grayscale images (e.g., `\*.nii.gz` files with
+continuous values).  `cardio` is launched from the commandline and may be configured via
+commandline arguments, a static TOML configuration file, or a combination of the two.
 
 ## Quickstart
 
@@ -76,7 +75,7 @@ $ uv init
 $ uv add cardio
 $ . ./.venv/bin/activate
 (project) cardio --version
-cardio 2025.10.1
+cardio 2026.1.0
 ```
 
 ### Developing

cardio-2026.1.0.dist-info/RECORD

@@ -0,0 +1,30 @@
+cardio/__init__.py,sha256=WZ7MfzISF8ryIC9NdIyub4mP6DQBWsMX766IPcWihXg,601
+cardio/app.py,sha256=h-Xh57Txw3Hom-fDv08YxjXzovuIxY8_5qPy6MJThZs,1354
+cardio/assets/bone.toml,sha256=vv8uVYSHIoKuHkNCoBOkGe2_qoEbXMvQO6ypm3mMOtA,675
+cardio/assets/vascular_closed.toml,sha256=XtaZS_Zd6NSAtY3ZlUfiog3T86u9Ii0oSutU2wBQy78,1267
+cardio/assets/vascular_open.toml,sha256=1M3sV1IGt3zh_3vviysKEk9quKfjF9xUBcIq3kxVHFM,879
+cardio/assets/xray.toml,sha256=siPem0OZ2OkWH0e5pizftpItJKGJgxKJ_S2K0316ubQ,693
+cardio/blend_transfer_functions.py,sha256=fkLDYGMj_QXYs0vmXYT_B1tgTfl3YICLYimtWlxrmbQ,2958
+cardio/color_transfer_function.py,sha256=uTyPdwxi0HjR4wm418cQN9-q9SspkkeqvLNqXaF3zzg,792
+cardio/logic.py,sha256=pjAqiJc_3pfUpB9hXWiSYJTmbG4c1omuUvuT2gshPNo,40431
+cardio/mesh.py,sha256=wYdU0BU84eXrrHp0U0VLwYW7ZpRJ6GbT5Kl_-K6CBzY,9356
+cardio/object.py,sha256=98A32VpFR4UtVqW8dZsRJR13VVyUoJJ20uoOZBgN4js,6168
+cardio/orientation.py,sha256=zOw3R284izAcKXnCCYYFOzOPBqgiE8om8q2JDWTQFoY,6256
+cardio/piecewise_function.py,sha256=X-_C-BVStufmFjfF8IbkqKk9Xw_jh00JUmjC22ZqeqQ,764
+cardio/property_config.py,sha256=YrWIyCoSAfJPigkhriIQybQpJantGnXjT4nJfrtIJco,1689
+cardio/rotation.py,sha256=eUUk74kBpPr-YCZdIAcU5pIX6FKeVGT8dRg2YnaS2yE,7204
+cardio/scene.py,sha256=2oq-zVMSHj9PhNToSMbk8RzJ3jKzo7oaAvDn_U4vcPM,15583
+cardio/screenshot.py,sha256=l8bLgxnU5O0FlnmsyVAzKwM9Y0401IzcdnDP0WqFSTY,640
+cardio/segmentation.py,sha256=KT1ClgultXyGpZDdpYxor38GflY7uCgRzfA8JEGQVaU,6642
+cardio/transfer_function_pair.py,sha256=_J0qXA0InUPpvfWPcW492KGSYAqeb12htCzKBSpWHwo,780
+cardio/types.py,sha256=jxSZjvxEJ03OThfupT2CG9UHsFklwbWeFrUozNXro2I,333
+cardio/ui.py,sha256=uYyO1-Lysi5guFDv4k3HSGJ12Ro27B1yxkmS_PFbd-M,53189
+cardio/utils.py,sha256=ao4a7_vMjGBxTOMhZ7r0D0W4ujiwKPS0i8Xfmn3Gv9k,1497
+cardio/volume.py,sha256=sdQ7gtX4jQDD9U8U95xD9LIV_1hpykj8NTCo1_aIKVM,15035
+cardio/volume_property.py,sha256=-EUMV9sWCaetgGjnqIWxPp39qrxEZKTNDJ5GHUgLMlk,1619
+cardio/volume_property_presets.py,sha256=4-hjo-dukm5sMMmWidbWnVXq0IN4sWpBnDITY9MqUFg,1625
+cardio/window_level.py,sha256=XMkwLAHcmsEYcI0SoHySQZvptq4VwX2gj--ps3hV8AQ,784
+cardio-2026.1.0.dist-info/WHEEL,sha256=4n27za1eEkOnA7dNjN6C5-O2rUiw6iapszm14Uj-Qmk,79
+cardio-2026.1.0.dist-info/entry_points.txt,sha256=xRd1otqKtW9xmidJ4CFRX1V9KWmsBbtxrgMtDColq3w,40
+cardio-2026.1.0.dist-info/METADATA,sha256=xulm4Dgn_80ZFMHSKEBtxnoEfYKFFf3Pfb3aNM1fgzw,3520
+cardio-2026.1.0.dist-info/RECORD,,

cardio-2025.10.1.dist-info/RECORD

@@ -1,29 +0,0 @@
-cardio/__init__.py,sha256=cVW5yofdKIbjqc0sf_VVnIMah9BvAgELy7IzZR50IZQ,602
-cardio/app.py,sha256=TEzgA03EAgI7HSCHhYwYb8tsnAptsrpcysV5CytfAq4,1379
-cardio/assets/bone.toml,sha256=vv8uVYSHIoKuHkNCoBOkGe2_qoEbXMvQO6ypm3mMOtA,675
-cardio/assets/vascular_closed.toml,sha256=XtaZS_Zd6NSAtY3ZlUfiog3T86u9Ii0oSutU2wBQy78,1267
-cardio/assets/vascular_open.toml,sha256=1M3sV1IGt3zh_3vviysKEk9quKfjF9xUBcIq3kxVHFM,879
-cardio/assets/xray.toml,sha256=siPem0OZ2OkWH0e5pizftpItJKGJgxKJ_S2K0316ubQ,693
-cardio/blend_transfer_functions.py,sha256=s5U4hO810oE434wIkPmAP2mrAfqFb4xxxi3hHf_k8og,2982
-cardio/color_transfer_function.py,sha256=KV4j11AXYeaYGeJWBc9I-WZf7Shrm5xjQVq-0bq9Qc8,817
-cardio/logic.py,sha256=a73JJ2b8t-TYggYMGI1IGpn0Q6AFD8I1JNz9fn6QeDk,37557
-cardio/mesh.py,sha256=XIJX2Ldv1Bw82yykE0N4IlQ0xE2r6xitjxsbEvEKrcI,9390
-cardio/object.py,sha256=FvMzhO9gyHhXgME8B6fcnWH74HmsYNoUjpVKUWEpWhM,6191
-cardio/orientation.py,sha256=J3bqZbv8vfl4loGl7ksmuyqWb3zFAz-TVSIahKcg0pc,6145
-cardio/piecewise_function.py,sha256=bwtwgrAMGkgu1krnvsOF9gRMaZb6smsS9jLrgBecSbo,789
-cardio/property_config.py,sha256=XJYcKeRcq8s9W9jqxzVer75r5jBLuvebv780FYdPV8U,1723
-cardio/scene.py,sha256=9jskdEARyJjk7QBIcMt5X2HsnikTAJdoRoLSI0LxcJE,14301
-cardio/screenshot.py,sha256=l8bLgxnU5O0FlnmsyVAzKwM9Y0401IzcdnDP0WqFSTY,640
-cardio/segmentation.py,sha256=SK9f834lld08qYH9l2Ligml212u7MzfcOEgTz8wnpzs,6676
-cardio/transfer_function_pair.py,sha256=90PQXByCL6mMaODo7Yfd-lmdFtCKhcBZbNaiH-PTds0,805
-cardio/types.py,sha256=DYDgA5QmYdU3QQrEgZMouEbMEIf40DJCeXo4V7cDXtg,356
-cardio/ui.py,sha256=wsRujwiNtMeq65R4bx-myxSPDhJT0RskwvBWlGenmHk,42901
-cardio/utils.py,sha256=tFUQ4FxfidTH6GjEIKQwguqhO9T_wJ2Vk0IhbEfxRGA,1616
-cardio/volume.py,sha256=k7NYZCwsmoAbs8mNWDC51gS6nlzPqD9JhDr1iyUQ0cs,10377
-cardio/volume_property.py,sha256=6T2r67SSIDl8F6ZlQvgMCZESLxuXVVAUjOC50lgQEmk,1644
-cardio/volume_property_presets.py,sha256=U2a2MnyCjryzOLEADs3OLSMMmAUnXq82mYK7OVXbQV0,1659
-cardio/window_level.py,sha256=gjk39Iv6jMJ52y6jydOjxBZBsI1nZQMs2CdWWTshQoE,807
-cardio-2025.10.1.dist-info/WHEEL,sha256=4n27za1eEkOnA7dNjN6C5-O2rUiw6iapszm14Uj-Qmk,79
-cardio-2025.10.1.dist-info/entry_points.txt,sha256=xRd1otqKtW9xmidJ4CFRX1V9KWmsBbtxrgMtDColq3w,40
-cardio-2025.10.1.dist-info/METADATA,sha256=0GkZB7yixZKd4nFaAHgWtgCA6VxlV_pbUYdHokFROD4,3540
-cardio-2025.10.1.dist-info/RECORD,,