Accurate MRI-free PET Quantification and Reliability Assessment via CT-to-MRI Synthesis and Conformal Prediction
Organizations: Seoul National University, Seoul, Republic of Korea · Yonsei University, Seoul, Republic of Korea
Abstract
Standardized uptake value ratio (SUVR) quantification of amyloid PET commonly relies on a subject-specific T1-weighted MRI, yet a usable MRI is not always available. We therefore combined CT-to-MRI synthesis with conformal prediction for MRI-free PET quantification with reliability assessment. A slice-consistent Brownian Bridge Diffusion Model was strengthened with multi-dataset training and data-source conditioning to enrich the anatomical prior learned from MRI, and with expanded multi-slice CT context to exploit the through-plane continuity of subject-specific anatomy in CT. Synthetic MRI was used for FreeSurfer-based anatomical parcellation and PET quantification, while deep-ensemble uncertainty was calibrated with conformal prediction to estimate a plausible range for the corresponding real MRI-based SUVR. On C-PiB PET from OASIS-3, the proposed method achieved the highest synthesis and parcellation fidelity among the compared methods. Mean cortical SUVR from the proposed synthetic MRI agreed closely with the real MRI-based reference (MAE , Pearson ) and yielded the same amyloid-positivity classification as the real MRI in all 87 test subjects at the prespecified cut-off. Conformal prediction achieved empirical coverage close to the prespecified levels, and at 90% target coverage, 85 of 87 subjects had prediction intervals that supported an unambiguous amyloid-positivity decision. Zero-shot evaluation on an independent external cohort showed that the proposed CT-to-MRI model could support PET quantification without retraining, although conformal calibration required target-cohort recalibration to recover the intended coverage. These findings support CT-to-MRI-based MRI-free PET quantification as a promising approach for providing accurate SUVR measurements and an explicit assessment of their reliability.
Figures & tables
| Severance ( ) | OASIS-3 ( ) | |||
| Method | PSNR | PSNR | ||
| Kim2025 | ||||
| LSGAN | ||||
| ALDM | ||||
| Choo2024 | ||||
| Ours | ||||
| Region | Metric | Kim2025 | LSGAN | ALDM | Choo2024 | Ours |
| Cerebral cortex | DSC | |||||
| HD95 (mm) | ||||||
| Cerebral WM | DSC | |||||
| HD95 (mm) | ||||||
| Cerebellar GM (ref.) | DSC | |||||
| HD95 (mm) |
| Pipeline | MAE | Mean bias | 95% LoA | Conc. | FN/FP | |
| PET-only | ||||||
| PET-template (GAAIN-derived) | 0.952 | 89.7% | 0/9 | |||
| PET-template (OASIS-derived) | 0.973 | 96.6% | 0/3 | |||
| CT-based | ||||||
| Presotto2018 | 0.989 | 97.7% | 1/1 | |||
| Matsuda2021 | 0.982 | 98.9% | 0/1 | |||
| Model | Frontal | Parietal | Lateral Temporal | Precuneus | Anterior Cingulate | Posterior Cingulate |
| Kim2025 | ||||||
| LSGAN | ||||||
| ALDM | ||||||
| Choo2024 | ||||||
| Ours |
Appendix figures & tables12 assets
Supplementary material from the paper’s appendix.
Appendix
| Item | Setting |
| Framework | PyTorch 2.6 |
| Optimizer | Adam, lr , bfloat16 autocast; EMA decay 0.995 |
| Effective batch size | 16 (a mini-batch of 2 with 8-step gradient accumulation on a single A100 GPU) |
| Synthesis input/output | 15 axial CT slices central 3 MRI slices (H3) |
| Diffusion backbone | BBDM, U-Net score network |
| Diffusion process | Brownian bridge with steps, linear , maximum variance 1, the network is trained to predict with an loss |
| Setting | Choo2024 | ALDM | Kim2025 | LSGAN |
| Implementation source | Released CT2MRI code [ S1 ] | Released network modules [ S2 ] | Re-implemented from the publication [ S3 ] | 3D pix2pix with least-squares loss [ S4 , S5 ] |
| Input output | 3 CT slices 3 MRI slices | 3D CT 3D MRI | 3 CT slices central MRI slice | 3D CT 3D MRI |
| Network | 2D BBDM U-Net; MRI histogram condition | 3D first-stage encoder/decoder with MS-SPADE; latent diffusion U-Net | 2.5D residual U-Net (16/32/64 encoder channels; six residual blocks) | 3D U-Net generator and 3D PatchGAN discriminator |
| Training objective | Brownian bridge, ; | First stage ; second-stage -prediction MSE, | Least-squares adversarial | |
| Optimizer | Adam, , ; EMA 0.995 | AdamW, , weight decay ; EMA 0.9999 | Adam, , | Adam, , |
| Effective batch | 16 (2 8 accumulation) | 4 (1 4 accumulation) | 16 | 4 (1 4 accumulation) |
| Implementation | Template or normalization settings | Target atlas / reference region |
| PET-only template-based normalization | ||
| PET-template (GAAIN-derived) | PiB template constructed from the public GAAIN AD-100 and YC-0 calibration PET/MRI datasets | FreeSurfer DK parcellation of the MNI152 template / GAAIN cerebellar GM |
| PET-template (OASIS-derived) | PiB template from 197 OASIS-3 subjects. Normalization used 8-mm source smoothing and default regularization. | Population atlas from 200 OASIS-3 FreeSurfer parcellations / cerebellar cortex |
| CT-based anatomical normalization | ||
| Presotto2018 | HU set to . Bias correction disabled. Gaussians per tissue class were 1/1/2/2/2/2. | DK atlas from atlasreader (cvs_avg35 in MNI152) / cerebellar cortex |
| Matsuda2021 | Skull-based initialization and CT coregistration to avg152T1. Default segmentation settings without HU thresholding. | DK atlas from atlasreader (cvs_avg35 in MNI152) / cerebellar cortex |
| Severance ( ) | OASIS-3 ( ) | |||
| Model | PSNR | PSNR | ||
| BBDM [ S1 ] | ||||
| H1 | ||||
| H1 H2 | ||||
| H1 H2 H3 (proposed) | ||||
| H4 (explored) | ||||
| Pipeline | Frontal | Parietal | Lat. temporal | Precuneus | Ant. cingulate | Post. cingulate |
| PET-template (GAAIN-derived) | ||||||
| PET-template (OASIS-derived) | ||||||
| Presotto2018 | ||||||
| Matsuda2021 | ||||||
| Ours |
| 80% target | 90% target | 95% target | ||||
| Region | Coverage | Width | Coverage | Width | Coverage | Width |
| Frontal | 77.0% (67/87) [67.1, 84.6] | 88.5% (77/87) [80.1, 93.6] | 95.4% (83/87) [88.8, 98.2] | |||
| Parietal | 83.9% (73/87) [74.8, 90.2] | 88.5% (77/87) [80.1, 93.6] | 89.7% (78/87) [81.5, 94.5] | |||
| Lateral temporal | 77.0% (67/87) [67.1, 84.6] | 90.8% (79/87) [82.9, 95.3] | 93.1% (81/87) [85.8, 96.8] | |||
| Precuneus | 85.1% (74/87) [76.1, 91.1] | 92.0% (80/87) [84.3, 96.0] | 94.3% (82/87) [87.2, 97.5] | |||
| Anterior cingulate | 82.8% (72/87) [73.5, 89.3] | 96.6% (84/87) [90.3, 98.8] | 96.6% (84/87) [90.3, 98.8] | |||
| Cut-off | Ref. +/– | FP/FN | Confident | Indet. | Captured | Conf. disag. |
| 1.30 | 24/63 | 0/0 | 84/87 | 3/87 | – | 0/84 |
| 1.35 | 22/65 | 1/1 | 83/87 | 4/87 | 2/2 | 0/83 |
| 1.40 | 21/66 | 1/0 | 86/87 | 1/87 | 1/1 | 0/86 |
| 1.42 | 21/66 | 0/0 | 85/87 | 2/87 | – | 0/85 |
| 1.45 | 20/67 | 0/0 | 85/87 | 2/87 | – | 0/85 |
| 1.50 | 20/67 | 0/0 | 86/87 | 1/87 | – | 0/86 |
| Measure | Target | OASIS-3 cal. GARD | GARD-specific calibration | ||
| Coverage | Width | Coverage | Width | ||
| Composite MCSUVR | |||||
| MCSUVR | 80% | 42.7% | 0.075 | % | 0.160 |
| 90% | 62.2% | 0.119 | % | 0.203 | |
| 95% | 78.0% | 0.151 | % | 0.252 | |
| Regional SUVR at 90% target | |||||
Medical Image Synthesis
Medical image synthesis creates, translates, reconstructs, and edits clinical images to fill missing data and support analysis.
See all 312 papers