Aligning Fetal Anatomy with Kinematic Tree Log-Euclidean PolyRigid Transforms
Authors: Yingcheng Liu, Athena Taymourtash, Yang Liu, Esra Abaci Turk, William M. Wells, Leo Joskowicz, P. Ellen Grant, Polina Golland
Organizations: Computer Science and Artificial Intelligence Lab, MIT, Cambridge, USA · California Institute of Technology, Pasadena, USA · Boston Children’s Hospital and Harvard Medical School, Boston, USA · The Hebrew University of Jerusalem, Jerusalem, Israel
Automated analysis of articulated bodies is crucial in medical imaging. Existing surface-based models often ignore internal volumetric structures and rely on deformation methods that lack anatomical consistency guarantees. To address this problem, we introduce a differentiable volumetric body model based on the Skinned Multi-Person Linear (SMPL) formulation, driven by a new Kinematic Tree-based Log-Euclidean PolyRigid (KTPolyRigid) transform. KTPolyRigid resolves Lie algebra ambiguities associated with large, non-local articulated motions, and encourages smooth, bijective volumetric mappings. Evaluated on 53 fetal MRI volumes, KTPolyRigid yields deformation fields with significantly fewer folding artifacts. Furthermore, our framework enables robust groupwise image registration and a label-efficient, template-based segmentation of fetal organs. It provides a robust foundation for standardized volumetric analysis of articulated bodies in medical imaging.
Figures & tables
Figure 1: Method overview. KTPolyRigid ΦT flattens the fetal body into a canonical T-pose. The flow-based mapping ΦP further standardizes the body shape to the population mean Vˉ . Averaging the population-level canonical images Js produces a population mean image Jˉ . Groupwise registration optimizes the body pose θ to minimize image variance across subjects. Segmentation is achieved by total transformation Φ=ΦT∘ΦP .
Figure 2: Example canonical image I~s of a subject s . Left three panels: Three coronal cross-sections (full and zoomed-in views). Right panel: mid-sagittal cross-section.
Figure 3: Local distortions. Spatial log determinant Jacobian field log2det[∂xΦT(x)] of the deformation ΦT for subject in Fig. 2 .
Method
% folds ( ↓ )
STDlog2∣detJ∣ ( ↓ )
GPU Time (ms, ↓ )
GPU Mem. (MB, ↓ )
LBS
1.88±0.57
0.89±0.08
0.33±0.03
24.27±0.95
PolyRigid [ 2 ]
5.37±3.85
0.90±0.18
93.50±6.13
3332±129
KTPolyRigid (ours)
1.38±0.40
0.74±0.06
31.19±2.69
1551±66
Table 1: Deformation field statistics. Mean ± standard deviation are reported. ( ↓ ) indicates that lower values represent better performance.
MIT Computer Science and Artificial Intelligence Laboratory · Massachusetts General Hospital, Harvard Medical School · Boston Children’s Hospital, Harvard Medical School
MIT Computer Science and Artificial Intelligence Laboratory, Cambridge, USA · Boston Children’s Hospital and Harvard Medical School, Boston, USA · Inria, Université Côte d’Azur, Sophia Antipolis, France