Authors: Galib Ahmed, Istiak Ahmed, Aritra Islam Saswato, Asib Mostakim Fony, Kazi Shahriar Sanjid, Md. Tanzim Hossain, Md. Anwarul Islam, Md. Nishan Khan, Md. Misbah Khan, Labiba Faiza Karim, Jobaer Rahman, S M Hasibul Hoque, Rahnuma Shahrin Rista, Kamruzzaman Rumman, Md Arifur Rahman, Syed Md. Akram Hussain, Mohammad Ashrafuzzaman Khan, M. Monir Uddin
Organizations: Department of Mathematics and Physics, North South University, Dhaka, 1229, Bangladesh · Big-Matrix Lab, Department of Mathematics and Physics, North South University, Dhaka, 1229, Bangladesh · Department of Electrical and Computer Engineering, North South University, Dhaka, 1229, Bangladesh · Department of Radiation Oncology, Dhaka, National Institute of Cancer Research & Hospital, 1212, Bangladesh · Department of Data Science, Erlangen, Friedrich-Alexander University, 91054, Germany · Square Cancer Centre, Square Hospitals Limited,, Dhaka 1205, Bangladesh · Department of Computer Science & Engineering,, Louisiana Tech University, 201 Mayfield Ave, Ruston, LA 71272, USA · Department of Radiation Oncology, Labaid Cancer Hospital & Super Speciality Center,, Dhaka 1205, Bangladesh · Department of Oncology & Radiotherapy, Bangladesh Specialized Hospital Limited,, Dhaka 1207, Bangladesh · Department of Clinical Oncology,, Bangladesh Medical University, Dhaka 1000, Bangladesh
We developed DAMM-Net++, a 2.5D architecture for thoracic OAR and target volume segmentation that addresses three persistent challenges in radiotherapy auto-contouring: inter-slice surface incoherence, systematic failure on small low-contrast targets, and the absence of per-case reliability signals. The central component is an anatomy-change-aware bidirectional selective state-space memory that models through-plane anatomical change and selectively propagates context along the axial slice sequence. A boundary-aware decoder sharpens near-surface predictions, and an uncertainty head provides calibrated per-voxel confidence for clinical triage. We evaluated 2,146 patients across four centers, an independent external cohort of 112 patients, and a multicenter reader study involving 17 radiation oncologists on 305 cases. The model achieves a mean Dice of 0.955 and HD95 of 3.78 mm, with the largest gains on low-contrast organs-at-risk (OARs) and target volumes where through-plane context is most critical. The uncertainty head is well-calibrated and supports case-level triage. In the reader study, AI assistance reduced contouring time by 75-80% across experience levels and raised junior-reader IoU from 0.861 to 0.925, matching the unedited model. External validation showed a modest internal-to-external drop (< 5%) with calibrated uncertainty transferring without recalibration. The complete deployment pipeline from DICOM ingestion to TPS-compatible RTSTRUCT export has been integrated into the clinical workflow at a partner hospital, where it is used to assist with contouring. These results suggest that anatomically motivated inter-slice memory, paired with uncertainty-guided review, offers a clinically viable path for thoracic auto-contouring.
Figures & tables
Institution
Patients
Proportion (%)
Acquisition Years
Role
Internal Cohort
Bangladesh Medical University (BMU)
1,447
67.4
2020–2026
Internal
Square Hospital Limited (SHL)
309
14.4
2018–2026
Internal
Labaid Hospital (LH)
278
13.0
2019–2025
Internal
External Validation Cohort
United Hospital (UHL)
112
N/A
2026
External Validation
Table 1: Institutional distribution of the multicenter thoracic radiotherapy planning cohort. BMU, SHL, and LH comprised the internal cohort split 70/15/15; UHL was reserved exclusively for external validation.
Fig. 1 : Overview of the proposed DAMM-Net++ architecture. The input is a sequence of T consecutive axial CT slices, each passed through a weight-shared ConvNeXt encoder. The Multi-Scale Inter-Slice Transition branch computes inter-slice feature differences; the Bidirectional Selective State-Space Memory accumulates anatomy-change-informed context in both through-plane directions; and the Memory-Guided Boundary-Aware Decoder reconstructs full-resolution predictions. The network produces segmentation, boundary, signed distance map, presence, and uncertainty outputs.
Fig. 2 : Detailed architecture of three core modules. (a) Multi-Scale Inter-Slice Transition Branch: multi-scale difference processing with channel and spatial attention. (b) Bidirectional Selective State-Space Memory: forward and backward through-plane sweeps with gated fusion. (c) Memory-Guided Cross-Attention: querying current features against memory context.
Model
Transition signal
Selective SSM
Dice ↑
Δ vs SSM only
Base
N/A
×
0.9047
N/A
SSM only
×
✓
0.9251
N/A
Difference + SSM
Real
✓
0.9350
↑0.0099∗∗
Randomized + SSM
Random
✓
0.9254
↓0.0003ns
Shuffled + SSM
Corrupted
✓
0.9282
↑0.0031∗
Full model
Real
✓
0.9545
↑0.0294∗∗
Table 3: Controlled ablation isolating the inter-slice transition signal. All variants include the ConvNeXt encoder and MGBA decoder. The selective SSM is disabled in the Base condition. Blue and rose mark the best and second-best values, respectively; green arrows indicate improvement relative to the SSM-only baseline; red arrows indicate regression. Critical comparison: real transition vs random/corrupted controls.
Structure
Base
+Mem
Full
Δ (Full − Base)
Spinal cord
0.8901
0.9204
0.9350
↑0.0449
Esophagus
0.8210
0.8788
0.9058
↑0.0848
Heart
0.9601
0.9668
0.9699
↑0.0098
Lung (R)
0.9704
0.9741
0.9764
↑0.0060
Lung (L)
0.9758
0.9788
0.9810
↑0.0052
Trachea
0.8642
0.9187
0.9425
↑0.0783
Table 4: Per-structure Dice. Base removes the inter-slice transition and memory modules; +Mem restores memory only; Full is DAMM-Net++. Gains concentrate on structures with high inter-slice variability. Blue and rose mark the best and second-best configuration per structure; green values show the Full-minus-Base improvement.
Variant
Dice ↑
ACE ↓
NLL ↓
AUROC ↑
Softmax-entropy proxy
0.9502
0.081
0.214
0.712
Variance head, std. loss
0.9518
0.052
0.176
0.784
Unc. head + UA loss (ours)
0.9545
0.031
0.142
0.861
Table 5: Uncertainty quality. Top : confidence calibration and error-detection performance. ACE is the absolute calibration error between predicted confidence and observed IoU. Bottom : Dice on the retained region as low-confidence voxels are routed to review. Blue and rose mark the best and second-best values; arrows quantify the gain over accepting all voxels.
Fig. 3 : Per-patient DSC distributions ( ↑ higher is better) for DAMM-Net++ and four baselines across nine thoracic structures. Paired violin plots show individual trajectories, box-and-whisker summaries, and paired Wilcoxon brackets (red: p<0.05 ; gray: ns) with Hedges’ g annotations. DAMM-Net++ achieves significantly higher DSC on eight of nine structures; Heart is the sole exception.
Struc- ture
Metric
Ours
nnUNet
nnMamba
UNETR++
Swin- UNETR
SC
DSC
0.935±0.014
0.912±0.015∥
0.917±0.012∥
0.918±0.010∥
0.914±0.013∥
IoU
0.878±0.024
0.841±0.022∥
0.847±0.021∥
0.846±0.018∥
0.843±0.021∥
HD95
1.585±0.504
7.523±8.50∥
7.721±9.00∥
7.234±8.20∥
7.698±8.80∥
ASD
0.671±0.217
6.321±8.00∥
6.497±8.50∥
5.912±7.80∥
6.654±8.40∥
NSD
0.995±0.007
0.977±0.036∥
0.978±0.035∥
0.979±0.035∥
0.976±0.037∥
ESO
DSC
0.906±0.020
0.908±0.016
0.903±0.017‡
0.902±0.018‡
0.907±0.016
Table 6: Per-structure, per-metric comparison of five models on the multi-institutional thoracic dataset (background excluded). Values are mean ± standard deviation over the test set; bold indicates the best result. Superscripts report paired Wilcoxon tests against DAMM-Net++, Holm-adjusted within each 36-comparison metric family: ∥padj<0.001 , ‡padj<0.01 , †padj<0.05 , and ns not significant. Effect sizes use Hedges’ g with BCa 95% CIs from 1,000 resamples. Structure abbreviations: SC, spinal cord; ESO, esophagus; HRT, heart; RL, right lung; LL, left lung; TRA, trachea; BDY, body; GTV, gross tumor volume; CTV, clinical target volume.
Fig. 4 : Through-plane consistency of the proposed model’s predictions on a representative thoracic case. Rows show axial (top), coronal (middle), and sagittal (bottom) planes reformatted from the same volumetric prediction. Columns show (1) original CT, (2) ground-truth overlay, and (3) DAMM-Net++ prediction overlay.
Fig. 5 : Contour overlay comparison on a representative GTV axial slice. Left: CT with ground truth and predictions; legend values report per-case IoU. Right: magnified boundary view. Color coding: Ground Truth (green), DAMM-Net++ (red), nnUNet (royal Blue), UNETR++ (yellow), SwinUNETR (magenta), nnMamba (cyan).
Fig. 6 : Voxel-level error maps for all five models. Rows: DAMM-Net++, nnUNet, UNETR++, SwinUNETR, nnMamba. Columns: (1) original CT, (2) ground truth, (3) prediction, (4) error overlay, (5) zoomed boundary error. Green, blue, and pink denote body, CTV, and GTV, respectively; red indicates prediction–reference disagreement.
Metric
Expert ( n=5 )
Senior ( n=5 )
Junior ( n=7 )
Model (unedited)
p† Exp vs Jun
p† Sen vs Jun
p† Model vs Jun
IoU: No-AI
Mean ± SD
0.9490 ± 0.0058
0.9047 ± 0.0092
0.8606 ± 0.0133
N/A
< 0.0001
< 0.0001
N/A
Median [IQR]
0.9493 [0.9450–0.9530]
0.9046 [0.8981–0.9107]
0.8606 [0.8515–0.8702]
N/A
IoU: AI-Assisted
Mean ± SD
0.9702 ± 0.0038
0.9487 ± 0.0062
0.9245 ± 0.0091
0.9226 ± 0.0070
0.0021
< 0.0001
< 0.0001
Median [IQR]
0.9702 [0.9676–0.9727]
0.9486 [0.9447–0.9531]
0.9246 [0.9184–0.9306]
N/A
Table 7: Reader study summary. Values are mean ± SD unless stated otherwise. IoU: Intersection over Union. Model (unedited): DAMM-Net++ predictions without reader editing. Expert ≥ 20 years ( n=5 ), Senior 10–15 years ( n=5 ), Junior < 10 years ( n=7 ).
Medical image segmentation delineates organs, lesions, and other structures in scans to support measurement, diagnosis, treatment planning, and research.
We present a clinically deployed end-to-end auto-contouring system for cervical cancer radiotherapy planning, anchored by the Boundary-Aware Transformer with Region-Aware Mamba (BAT-RM), a hybrid architecture that integrates Sobel-gated boundary attention, a linear-time, multi-directional Mamba module for long-range context, and a boundary-skeleton-guided fusion gate. This design achieves linear-time complexity for long-range context modeling, avoiding the quadratic cost of full spatial self-attention. The full pipeline spans multi-institutional data collection, rigorous inter-rater quality assurance, external validation in an independent cohort, and a web-based clinical interface natively compatible with Varian, RayStation, and Monaco. Against four baselines, BAT-RM achieves superior performance across seven anatomical classes, with statistically significant improvements in target volumes, including GTV and CTV, and in organs at risk such as the rectum and bladder. A prospective multi-center reader study involving 13 radiation oncologists demonstrated that AI assistance elevates junior oncologists' IoU from 0.899 to 0.965, approaching senior-level accuracy, while reducing contouring time by more than 80%. The system also reduced expert consultation rates and improved inter-reader consistency, reflecting gains in both efficiency and quality assurance. Following clinical deployment at a partner hospital, the system reduced patient wait times from days to hours without additional staffing, enabling same-day or next-day initiation of treatment for routine cases. BAT-RM demonstrates that a rigorous research pipeline, from data curation to clinical deployment, can translate directly into measurable patient benefit in resource-constrained settings where the demand for radiotherapy far exceeds specialist capacity.
Istiak Ahmed, Kazi Shahriar Sanjid, Galib Ahmed +16
Department of Electrical and Computer Engineering, North South University, Dhaka, 1229, Bangladesh · Big-Matrix Lab, Department of Mathematics and Physics, North South University, Dhaka, 1229, Bangladesh · Department of Data Science, Erlangen, Friedrich-Alexander University, 91054, Germany +8
Four-dimensional computed tomography (4DCT) captures the full respiratory cycle of thoracic anatomy, yet current Internal Target Volume contouring workflows process each phase in isolation, discarding temporal coherence and leaving contours vulnerable to phase-specific artifacts. We present a lightweight framework that applies parameter-efficient fine-tuning to the Segment Anything Model 3 (SAM 3) via low-rank adaptation (LoRA) to align its text-prompted segmentation with the medical domain using only seven annotated 3D CT volumes. Furthermore, the framework incorporates a hard negative mining strategy to improve boundary discrimination in low-contrast thoracic regions. At inference, phase-wise predictions are refined through phase-coherent temporal filtering and spatial connectivity analysis. Since respiratory motion is continuous and periodic, genuine anatomy appears in contiguous blocks of phases, whereas transient artifacts appear sporadically and are thus effectively suppressed. Experiments on pulmonary and cardiac structures yield median Dice scores of 0.968 and 0.910 with 95th-percentile Hausdorff distances of 0.998 mm and 2.931 mm, respectively. The proposed framework effectively eliminates the severe false-positive predictions inherent in the zero-shot inference of the unadapted SAM 3. With only seven annotated volumes, the framework retains over 95% of full-data accuracy, and the entire pipeline is trainable on a single consumer-grade GPU, demonstrating a scalable, data-efficient solution for adaptive radiotherapy.
Changwoo Song
Oncosoft Inc., Seoul, South Korea · Department of Computer Science & Engineering, Chungnam National University, Daejeon, South Korea
Background: Accurate segmentation of the Left Anterior Descending (LAD) artery in 3D free-breathing, non-contrast CT is critical for cardiac dose sparing in thoracic radiotherapy. The LAD is extremely small, has poor soft-tissue contrast, and varies substantially across patients; even manual contours show limited inter-observer agreement, underscoring the ambiguity of the vessel boundaries. Purpose: To develop a transformer-based framework that improves LAD delineation in low-contrast, imbalanced CT through local-global context modeling and uncertainty-guided optimization. Methods: We propose NA-UNETR, a 3D transformer-based segmentation model whose Neighborhood Attention (NA) and Dilated NA (DiNA) blocks jointly capture fine structural detail and long-range context. Given the scarcity of annotated LAD data, the model is pretrained on 1,000 CTA volumes of general coronary anatomy and fine-tuned with LoRA-based parameter-efficient adaptation on 20 free-breathing institutional CT scans. A composite Dice-Focal and Hausdorff loss, dynamically balanced via homoscedastic uncertainty, improves overlap and boundary accuracy. Results: NA-UNETR reached 45.64% Dice, 38.16 mm HD95, and 10.01 mm ASD, improving Dice by 3.10 percentage points over nnU-Net and reducing HD95 by 2.96 mm relative to Swin UNETR, with the strongest boundary accuracy among all models and improved centerline stability. On ImageCAS it achieved 79.49% Dice, 8.89 mm HD95, and 1.02 mm ASD. Ablations confirmed that residual blocks, variable kernels, and uncertainty-weighted loss each contributed. Conclusions: NA-UNETR balances local precision and global context for thin, low-contrast LAD structures, offering a computationally efficient framework for substructure-level cardiac segmentation in radiotherapy planning.
Rafi Ibn Sultan, Chengyin Li, Yiannos Demetriou +6
Department of Computer Science, Wayne State University, Detroit, MI, 48202, USA · Department of Radiation Oncology, Henry Ford Health, Detroit, MI, 48202, USA · Department of Clinical Oncology, Alexandria University, Alexandria, 21526, Egypt +4