GLIDE-Reg: Global-to-Local Deformable Registration Using Co-Optimized Foundation and Handcrafted Features
Authors: Yunzheng Zhu, Aichi Chien, Kimaya kulkarni, Luoting Zhuang, Stephen Park, Ricky Savjani, Daniel Low, William Hsu
Organizations: Medical & Imaging Informatics, Department of Radiological Sciences, David Geffen School of Medicine at UCLA, USA · Magnetic Resonance Research Labs, Department of Radiological Sciences, David Geffen School of Medicine at UCLA, USA · Department of Radiation Oncology, David Geffen School of Medicine at UCLA, Los Angeles, USA
Deformable registration is crucial in medical imaging. Several existing applications include lesion tracking, probabilistic atlas generation, and treatment response evaluation. However, current methods often lack robustness and generalizability across two key factors: spatial resolution and differences in anatomical coverage. We jointly optimize a registration field and a learnable dimensionality reduction module so that compressed VFM embeddings remain registration-relevant, and fuse these global semantic cues with MIND local descriptors. GLIDE-Reg achieves average dice similarity coefficients (DSC) across 6 anatomical structures of 0.859, 0.862, and 0.901 in two public cohorts (Lung250M and NLST) and one institution cohort (UCLA5DCT), and outperforms the state-of-the-art DEEDS (0.834, 0.858, 0.900) with relative improvements of 3.0%, 0.5%, and 0.1%. For target registration errors, GLIDE-Reg achieves 1.58 mm on Lung250M landmarks (compared to 1.25 mm on corrField and 1.91 mm on DEEDS) and 1.11 mm on NLST nodule centers (compared to 1.11 mm on DEEDS). The substantiated performance on the nodule centers also demonstrates its robustness across challenging downstream tasks, such as nodule tracking, which is an essential prior step for early-stage lung cancer diagnosis.
Figures & tables
Dataset
NLST
Lung250M
UCLA5DCT
Study Disease
Indeterminate Pulmonary Nodule
COPD
N/A
Study Age
55 to 74
48 to 80
N/A
Collection Date
1999 to 2001
After 2007 (Phase I or II)
N/A
Task (Time Diff.)
Longitudinal ( 1 to 2 years)
Insp. Exp. ( ∼30 mins)
Free-Breathing ( <5 mins)
Data Type
Low-Dose CT
Breath-Hold CT
Free-Breath CT
Resolution
512×512×[67,196]
512×512×[112,135]
500×500×[342,411]
Table 1: Specifications of Datasets used for this study.
Figure 1: GLIDE-Reg Pipeline. A pair of 3D images Ifix and Imov are passed to a global feature extractor and a local feature extractor (MIND) for extracting the optimal global feature representations GFfix∗ , GFmov∗ and local feature representations LFfix , LFmov , respectively. Then the features are registered with 1) feature-independent coupled convex optimization and 2) fused Adam optimization in Local-Global Registration and output the Iwarp . DR stands for dimensionality reduction.
Lung250M
NLST
[blinded]5DCT
Method
TRE ↓
DSC ↑
%∣J∣<0↓
Runtime
TRE ↓
DSC ↑
%∣J∣<0↓
Runtime
DSC ↑
%∣J∣<0↓
Landmarks
Mean6
Nodule
Mean6
Mean6
No Reg.
16.25±3.84 *
0.517±0.080 *
−
−
18.93±10.27 *
0.472±0.100 *
−
−
0.767±0.049 *
−
Affine
30.36±12.12 *
0.488±0.106 *
−
<5 sec
30.36±23.88 *
0.726±0.045 *
−
<10 sec
0.751±0.050 *
−
SyN [ 1 ]
13.95±3.63 *
0.634±0.085 *
0.0
<1 min
17.65±10.24 *
0.534±0.125 *
<0.001
<2 min
0.838±0.036 *
<0.001
NiftyReg [ 21 ]
14.15±3.93 *
0.563±0.086 *
0.565±0.006
<2.5 min
26.93±14.25 *
0.396±0.088 *
1.846±1.992
<4 min
0.807±0.032 *
0.010±0.013
Table 2: Comparison of registration performance on three datasets evaluated by DSC( ↑ ) and TRE( ↓ ) (in mm ), and nonpositive Jacobian determinant %∣J∣<0 ( ↓ ). No Reg. represents the baseline performance for the image pairs. Methods with top two performance for TRE and DSC in each column are bolded and underlined , respectively. Statistical testings are computed per pair on the TRE and average DSC using Wilcoxon signed-rank test and adjusted by Benjamini-Hochberg correction with false discovery rate of 0.05 . Methods that are outperformed by or outperform GLIDE-Reg with statistical significance ( p<0.05 ) are marked with ∗ and †, respectively.
Figure 2: Comparison of the best performing 9 algorithms (best four in the learning-based and all five feature-based instance-optimization algorithms) on the NLST dataset.
Figure 3: Comparison of the feature-based instance-optimization methods on one challenging case from Lung250M, where all learning-based methods failed.
Intelligent Modelling & Analysis Group (IMA), School of Computer Science, University of Nottingham, UK · Lab for Uncertainty in Data and Decision Making (LUCID), Nottingham Biomedical Research Centre (BRC), School of Medicine, University of Nottingham, UK
Department of Information Technology, Uppsala University, Uppsala, Sweden · Department of Medical Physics, Uppsala University Hospital, Uppsala, Sweden · Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden