SPARC: Slice-to-volume Pipeline for Automated Reconstruction of gated 3D+time fetal Cardiac MRI
Authors: Arnaud Boutillon, Naomi Clarke, Tomas Woodgate, Daniel West, Alina Schneider, Rachael Franklin, Anthony Price, Jo Hajnal, +3 more
Organizations: Biomedical Computing Research Department, School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom · Imaging Physics and Engineering Research Department, School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom · Early Life Imaging Research Department, School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom · Medical Physics and Clinical Engineering, Guys and St Thomas’ NHS Foundation Trust, London, United Kingdom
Fetal cardiac MRI (fCMR) provides valuable diagnostic information complementary to echocardiography, particularly for complex congenital heart disease (CHD). Dynamic cine imaging captures cardiac motion essential for assessment of cardiac function; however, the reconstruction of 3D+time cine volumes from 2D+time acquired slices remains challenging due to unpredictable fetal motion and the absence of automated and robust processing tools suitable for clinical deployment. We present the SPARC pipeline (Slice-to-volume Pipeline for Automated Reconstruction of gated 3D+time fetal Cardiac MRI) which combines physics-informed slice-to-volume reconstruction (SVR) of Doppler ultrasound (DUS) gated stacks of slices, assisted by deep learning (DL) models for thoracic segmentation and anatomical reorientation. The proposed SVR algorithm achieves a tenfold reduction in reconstruction time relative to existing frame-wise approaches (4.8±1.0 vs 49.0±14.1 min, p<0.0001) while improving the reconstruction quality. Thoracic segmentation performance using ensemble aggregation exceeded inter-rater agreement (Dice 84.7±3.9% vs 81.4±7.7%, p<0.05), while anatomical reorientation achieved a success rate of 90.1%. End-to-end evaluation on a large held-out clinical cohort (n=121) demonstrated fully automatic processing in 82.6% of cases with a mean runtime of 7.1±1.3 min, compatible with clinical deployment. The complete SPARC pipeline is publicly available as a Docker container https://hub.docker.com/r/aboutill/sparc and is currently deployed at our institution as a clinical research tool.
Figures & tables
Figure 3: Overview of the proposed SPARC pipeline.
Domain
Evaluation
Subjects n
5-fold cross- validation set
Held-out test set
Unique subjects
Source
Thoracic segmentation
116 [483]
23 (92)
141
Cardiac segmentation
98
22
Reorientation
91
21
Target
Thoracic segmentation
25 [222]
5 (36)
60
Reconstruction
–
10
Table 1: Dataset splits for each evaluated task in the source (free-running; n=141 subjects) and target (Doppler-gated; n=199 subjects) domains. Numbers of stacks are shown in brackets.
Domain
Ensemble
DSC [%] ↑
HD95% [mm] ↓
ASSD [mm] ↓
Test
Training
Source
Source
✓
84.9±7.0
4.1±3.0
1.1±0.7
Joint
✓
85.7±6.0
4.1±3.2
1.1±0.7
Target
Source
✓
79.2±7.3
10.8±10.6
2.1±2.1
Target
✓
81.8±4.5
6.5±2.8
1.4±0.4
Joint
✓
84.3±3.9
5.9±2.8
1.2±0.4
Table 2: Quantitative evaluation of thoracic segmentation models on held-out test sets from the source (free-running; m=92 stacks across n=23 subjects) and target (Doppler-gated; m=36 stacks across n=5 subjects) domains shown as mean ± standard deviation. The best performance within each domain is shown in bold (linear mixed-effect model with Bonferroni correction, p<0.05 ).
Method
van Amerom
Proposed
Model assumptions
Motion states
Frame-wise
Slice-wise
tPSF
✓
✗
Step complexity
Slice/frame-to-volume registration
O(KT)
O(K)
Super-resolution reconstruction
O(NT2)
O(NT)
Leave-stack -out
NCC ↑
0.89±0.04
0.91±0.03
NRMSE ↓
0.13±0.03
0.11±0.02
Table 3: Quantitative evaluation of the reconstruction algorithms on the target domain test set (Doppler-gated; n=10 subjects) shown as mean ± standard deviation. The best performance within each domain is shown in bold (paired t -test, p<0.05 )
Domain
Ensemble
DSC [%] ↑
HD95% [mm] ↓
ASSD [mm] ↓
CD [mm] ↓
Test
Training
Source
Source
✓
83.4±4.9
4.8±1.9
1.7±0.4
4.2±2.3
Joint
✓
83.8±5.0
4.7±1.5
1.6±0.4
4.0±2.3
Target
Source
✓
20.0±30.8
27.8±12.0
14.4±7.7
39.5±22.3
Target
✓
73.8±9.0
8.2±3.6
3.2±1.4
12.3±5.2
Joint
✓
71.0±7.1
9.2±1.4
3.5±0.9
11.1±7.6
Table 4: Quantitative evaluation of cardiac segmentation models on held-out test sets from the source (free-running; n=22 subjects) and target (Doppler-gated; n=5 subjects) domains shown as mean ± standard deviation.
Domain
Ensemble
GD [ ° ] ↓
CD [mm] ↓
Test
Training
Source
Source
Chordal
13.1±6.3
1.6±1.5
Joint
Chordal
10.5±5.8
1.1±1.0
Target
Source
Chordal
74.9±52.8
5.2±4.5
Target
Chordal
33.4±36.4
4.4±4.5
Joint
Chordal
16.0±11.8
2.9±2.5
Table 5: Quantitative evaluation of reorientation models on held-out test sets from the source (free-running; n=21 subjects) and target (Doppler-gated; n=7 subjects) domains shown as mean ± standard deviation. Each image was sampled 250 times using independently drawn random rigid transformations. The best performance within each domain is shown in bold (linear mixed-effect model comparison with Bonferroni correction, p<0.05 ). † Target domain GD inter-rater agreement was computed on a single reorientation per subject ( n=7 ) and is not directly comparable to network performance assessed over whole range of rotations.
SPARC pipeline evaluation on clinical cohort ( n=121)
1) Manual intervention summary
Step
n ( % )
Time per case
Manual stack exclusion
6 (5.0%)
∼ 30 s
Manual template selection
3 (2.5%)
∼ 60 s
Manual thoracic segmentation
3 (2.5%)
∼ 120 s
Manual anatomical reorientation
12 (9.9%)
∼ 60 s
Table 6: End-to-end SPARC pipeline evaluation on a held-out clinical cohort ( n=121 subjects). 1) Manual intervention summary: the number and estimated duration of manual interventions required at each pipeline stage. 2) Template stack selection evaluation. 3) Thoracic segmentation evaluation. 4) Slice-to-volume reconstruction evaluation with mean processing time. 5) Anatomical reorientation evaluation stratified by cardiac centre localisation ablation. 6) End-to-end pipeline summary: overall counts and mean processing time for fully automatic and semi-automatic runs.
Figure S2: Geodesic distance (GD) of the Transfer approach as a function of simulated rotation magnitude ∣θ∣ , evaluated on the held-out target domain test set (Doppler-gated; n=7 subjects, 250 randomly sampled rigid transformations per subject). Each point represents one simulated transformation, with subjects colour-coded.
Reconstructing cardiac motion from CMR sequences is critical for diagnosis, prognosis, and intervention. Existing methods rely on complete CMR stacks to infer full heart motion, limiting their applicability during intervention when only sparse observations are available. We present TetHeart, the first end-to-end framework for unified 4D heart mesh recovery from both offline full-stack and intra-procedural sparse-slice observations. Our method leverages deformable tetrahedra to capture shape and motion in a coherent space shared across cardiac structures. Before a procedure, it initializes detailed, patient-specific heart meshes from high-quality full stacks, which can then be updated using whatever slices can be obtained in real-time, down to a single one during the procedure. TetHeart incorporates several key innovations: (i) an attentive slice-adaptive 2D-3D feature assembly mechanism that integrates information from arbitrary numbers of slices at any position; (ii) a distillation strategy to ensure accurate reconstruction under extreme sparsity; and (iii) a weakly supervised motion learning scheme requiring annotations only at keyframes, such as the end-diastolic and end-systolic phases. Trained and validated on three large public datasets and evaluated zero-shot on additional private interventional and public datasets without retraining, TetHeart achieves state-of-the-art accuracy and strong generalization in both pre- and intra-procedural settings. Code and dataset is available at https://github.com/Scalsol/TetHeart.
Yihong Chen, Jiancheng Yang, Deniz Sayin Mercadier +3
School of Computer and Communication Science, EPFL · ELLIS Institute Finland, Aalto University, Finland · University of North Carolina at Charlotte +1
Reconstructing a three-dimensional left-ventricular (LV) endocardial surface from cardiac magnetic resonance (CMR) data is challenging when supervision is available on only a small number of axial slices. Through-plane geometry is weakly constrained, and automatically generated two-dimensional masks can propagate segmentation errors into the recovered shape. We present MR-RS-SDFR, a per-case implicit signed distance field (SDF) framework that reconstructs a continuous LV surface from a CMR volume and sparse axial weak masks. The method first builds a cross-slice SDF initialization from axial and longitudinal geometric cues and then refines the field using two complementary signals: MRI edge-field normal alignment, which provides an image-derived boundary cue independent of the weak masks, and differentiable reslice Dice and contour consistency, which preserve agreement with the observed planes. We evaluate three weak-mask generators -- LOO TransUNet, LOO nnU-Net, and an off-the-shelf Medical SAM3 model used without MM-WHS-specific training or fine-tuning -- and five sparsity levels from 4 to 64 axial planes. In the sparse-16 setting, final MR-RS-SDFR reconstruction reaches 0.928 Dice and 3.80mm HD95 with Medical SAM3 masks. The upstream generators do not exhibit a single common ranking across 2D and dense 3D segmentation, and nnU-Net- and Medical-SAM3-driven sparse reconstruction achieve the same mean final Dice despite different upstream error profiles. Across all three sparse-16 mask sources, MR-RS-SDFR is numerically better than protocol-matched full GHD+DVS in both Dice and HD95. Final Dice improves markedly from sparse-4 to sparse-16 and then saturates at the reported precision through sparse-64. These results support MRI-guided per-case SDF refinement as a reconstruction strategy that remains effective across weak-mask generators and supervision densities.
Quanxin Zheng, Shuai Zhao
1China Electronics (Beijing) Information Technology Research Institute Co., Ltd.
Accurate 3D+t whole-heart mesh reconstruction from cine MRI is a clinically crucial yet technically challenging task. The difficulty of this task arises from two coupled factors: inherently sparse sampling of 3D cardiac anatomy by 2D image slices and the tight coupling between cardiac shape and motion. Current cardiac image-to-mesh approaches typically reconstruct only a subset of cardiac chambers or a single phase of the cardiac cycle. In this work, we propose CineMesh4D, a novel end-to-end 4D (3D+t) pipeline that directly reconstructs patient-specific whole-heart mesh from multi-view 2D cine MRI via cross-domain mapping. Specifically, we introduce a differentiable rendering loss that enables supervision of 3D+t whole-heart mesh from multi-view sparse contours of cine MRI. Furthermore, we develop a dual-context temporal block that fuses global and local cardiac temporal information to capture high-dimensional sequential patterns. In quantitative and qualitative evaluations, CineMesh4D outperforms existing approaches in terms of reconstruction quality and motion consistency, providing a practical pathway for personalized real-time cardiac assessment. The code will be publicly released once the manuscript is accepted.
Xiaoyue Liu, Xiaohan Yuan, Mark Y Chan +2
Department of Biomedical Engineering, National University of Singapore, Singapore · School of Automation, Southeast University, Nanjing, China · Department of Medicine, National University of Singapore, Singapore +1