Learning Whole Heart Mesh Generation From Patient Images for Computational Simulations

Probabilistic learning of the Purkinje network from the electrocardiogram

The identification of the Purkinje conduction system in the heart is a challenging task, yet essential for a correct definition of cardiac digital twins for precision cardiology. Here, we propose a probabilistic approach for identifying the Purkinje network from non-invasive clinical data such as the standard electrocardiogram (ECG). We use cardiac imaging to build an anatomically accurate model of the ventricles; we algorithmically generate a rule-based Purkinje network tailored to the anatomy; we simulate physiological electrocardiograms with a fast model; we identify the geometrical and electrical parameters of the Purkinje-ECG model with Bayesian optimization and approximate Bayesian computation. The proposed approach is inherently probabilistic and generates a population of plausible Purkinje networks, all fitting the ECG within a given tolerance. In this way, we can estimate the uncertainty of the parameters, thus providing reliable predictions. We test our methodology in physiological and pathological scenarios, showing that we are able to accurately recover the ECG with our model. We propagate the uncertainty in the Purkinje network parameters in a simulation of conduction system pacing therapy. Our methodology is a step forward in creation of digital twins from non-invasive data in precision medicine. An open source implementation can be found at http://github.com/fsahli/purkinje-learning.

Dynamic Virtual Simulation with Real-Time Haptic Feedback for Robotic Internal Mammary Artery Harvesting

Coronary heart disease, a leading global cause of mortality, has witnessed significant advancement through robotic coronary artery bypass grafting (CABG), with the internal mammary artery (IMA) emerging as the preferred "golden conduit" for its exceptional long-term patency. Despite these advances, robotic-assisted IMA harvesting remains challenging due to the absence of force feedback, complex surgical maneuvers, and proximity to the beating heart. This study introduces a novel virtual simulation platform for robotic IMA harvesting that integrates dynamic anatomical modeling and real-time haptic feedback. By incorporating a dynamic cardiac model into the surgical scene, our system precisely simulates the impact of cardiac pulsation on thoracic cavity operations. The platform features high-fidelity representations of thoracic anatomy and soft tissue deformation, underpinned by a comprehensive biomechanical framework encompassing fascia, adipose tissue, and vascular structures. Our key innovations include a topology-preserving cutting algorithm, a bidirectional tissue coupling mechanism, and dual-channel haptic feedback for electrocautery simulation. Quantitative assessment using our newly proposed Spatial Asymmetry Index (SAI) demonstrated significant behavioral adaptations to cardiac motion, with dynamic scenarios yielding superior SAI values compared to static conditions. These results validate the platform's potential as an anatomically accurate, interactive, and computationally efficient solution for enhancing surgical skill acquisition in complex cardiac procedures.

AI-Powered Multimodal Modeling of Personalized Hemodynamics in Aortic Stenosis

Aortic stenosis (AS) is the most common valvular heart disease in developed countries. High-fidelity preclinical models can improve AS management by enabling therapeutic innovation, early diagnosis, and tailored treatment planning. However, their use is currently limited by complex workflows necessitating lengthy expert-driven manual operations. Here, we propose an AI-powered computational framework for accelerated and democratized patient-specific modeling of AS hemodynamics from computed tomography (CT). First, we demonstrate that the automated meshing algorithms can generate task-ready geometries for both computational and benchtop simulations with higher accuracy and 100 times faster than existing approaches. Then, we show that the approach can be integrated with fluid-structure interaction and soft robotics models to accurately recapitulate a broad spectrum of clinical hemodynamic measurements of diverse AS patients. The efficiency and reliability of these algorithms make them an ideal complementary tool for personalized high-fidelity modeling of AS biomechanics, hemodynamics, and treatment planning.

SDF4CHD: Generative modeling of cardiac anatomies with congenital heart defects

Congenital heart disease (CHD) encompasses a spectrum of cardiovascular structural abnormalities, often requiring customized treatment plans for individual patients. Computational modeling and analysis of these unique cardiac anatomies can improve diagnosis and treatment planning and may ultimately lead to improved outcomes. Deep learning (DL) methods have demonstrated the potential to enable efficient treatment planning by automating cardiac segmentation and mesh construction for patients with normal cardiac anatomies. However, CHDs are often rare, making it challenging to acquire sufficiently large patient cohorts for training such DL models. Generative modeling of cardiac anatomies has the potential to fill this gap via the generation of virtual cohorts; however, prior approaches were largely designed for normal anatomies and cannot readily capture the significant topological variations seen in CHD patients. Therefore, we propose a type- and shape-disentangled generative approach suitable to capture the wide spectrum of cardiac anatomies observed in different CHD types and synthesize differently shaped cardiac anatomies that preserve the unique topology for specific CHD types. Our DL approach represents generic whole heart anatomies with CHD type-specific abnormalities implicitly using signed distance fields (SDF) based on CHD type diagnosis. To capture the shape-specific variations, we then learn invertible deformations to morph the learned CHD type-specific anatomies and reconstruct patient-specific shapes. After training with a dataset containing the cardiac anatomies of 67 patients spanning 6 CHD types and 14 combinations of CHD types, our method successfully captures divergent anatomical variations across different types and the meaningful intermediate CHD states across the spectrum of related CHD diagnoses. Additionally, our method demonstrates superior performance in CHD anatomy generation in terms of CHD-type correctness and shape plausibility. It also exhibits comparable generalization performance when reconstructing unseen cardiac anatomies. Moreover, our approach shows potential in augmenting image-segmentation pairs for rarer CHD types to significantly enhance cardiac segmentation accuracy for CHDs. Furthermore, it enables the generation of CHD cardiac meshes for computational simulation, facilitating a systematic examination of the impact of CHDs on cardiac functions.

Robust automated calcification meshing for personalized cardiovascular biomechanics

Calcification has significant influence over cardiovascular diseases and interventions. Detailed characterization of calcification is thus desired for predictive modeling, but calcium deposits on cardiovascular structures are still often manually reconstructed for physics-driven simulations. This poses a major bottleneck for large-scale adoption of computational simulations for research or clinical use. To address this, we propose an end-to-end automated image-to-mesh algorithm that enables robust incorporation of patient-specific calcification onto a given cardiovascular tissue mesh. The algorithm provides a substantial speed-up from several hours of manual meshing to ~1 min of automated computation, and it solves an important problem that cannot be addressed with recent template-based meshing techniques. We validated our final calcified tissue meshes with extensive simulations, demonstrating our ability to accurately model patient-specific aortic stenosis and Transcatheter Aortic Valve Replacement. Our method may serve as an important tool for accelerating the development and usage of personalized cardiovascular biomechanics.

LinFlo-Net: A Two-Stage Deep Learning Method to Generate Simulation Ready Meshes of the Heart

We present a deep learning model to automatically generate computer models of the human heart from patient imaging data with an emphasis on its capability to generate thin-walled cardiac structures. Our method works by deforming a template mesh to fit the cardiac structures to the given image. Compared with prior deep learning methods that adopted this approach, our framework is designed to minimize mesh self-penetration, which typically arises when deforming surface meshes separated by small distances. We achieve this by using a two-stage diffeomorphic deformation process along with a novel loss function derived from the kinematics of motion that penalizes surface contact and interpenetration. Our model demonstrates comparable accuracy with state-of-the-art methods while additionally producing meshes free of self-intersections. The resultant meshes are readily usable in physics based simulation, minimizing the need for postprocessing and cleanup.

Image2Flow: A proof-of-concept hybrid image and graph convolutional neural network for rapid patient-specific pulmonary artery segmentation and CFD flow field calculation from 3D cardiac MRI data

Computational fluid dynamics (CFD) can be used for non-invasive evaluation of hemodynamics. However, its routine use is limited by labor-intensive manual segmentation, CFD mesh creation, and time-consuming simulation. This study aims to train a deep learning model to both generate patient-specific volume-meshes of the pulmonary artery from 3D cardiac MRI data and directly estimate CFD flow fields. This proof-of-concept study used 135 3D cardiac MRIs from both a public and private dataset. The pulmonary arteries in the MRIs were manually segmented and converted into volume-meshes. CFD simulations were performed on ground truth meshes and interpolated onto point-point correspondent meshes to create the ground truth dataset. The dataset was split 110/10/15 for training, validation, and testing. Image2Flow, a hybrid image and graph convolutional neural network, was trained to transform a pulmonary artery template to patient-specific anatomy and CFD values, taking a specific inlet velocity as an additional input. Image2Flow was evaluated in terms of segmentation, and the accuracy of predicted CFD was assessed using node-wise comparisons. In addition, the ability of Image2Flow to respond to increasing inlet velocities was also evaluated. Image2Flow achieved excellent segmentation accuracy with a median Dice score of 0.91 (IQR: 0.86-0.92). The median node-wise normalized absolute error for pressure and velocity magnitude was 11.75% (IQR: 9.60-15.30%) and 9.90% (IQR: 8.47-11.90), respectively. Image2Flow also showed an expected response to increased inlet velocities with increasing pressure and velocity values. This proof-of-concept study has shown that it is possible to simultaneously perform patient-specific volume-mesh based segmentation and pressure and flow field estimation using Image2Flow. Image2Flow completes segmentation and CFD in ~330ms, which is ~5000 times faster than manual methods, making it more feasible in a clinical environment.

Computational approaches for mechanobiology in cardiovascular development and diseases

The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.