pith. sign in

arxiv: 2505.18782 · v6 · submitted 2025-05-24 · ⚛️ physics.med-ph

Joint Reconstruction of Activity and Attenuation in PET by Diffusion Posterior Sampling in Wavelet Coefficient Space

Pith reviewed 2026-05-19 14:14 UTC · model grok-4.3

classification ⚛️ physics.med-ph
keywords joint reconstructionPETattenuation correctiondiffusion posterior samplingwavelet diffusion modeltime-of-flightemission tomography
0
0 comments X

The pith

Joint PET activity and attenuation maps can be reconstructed from emission data alone via diffusion posterior sampling in wavelet space.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a technique to recover both the radioactive activity distribution and tissue attenuation properties in PET scans directly from the emission measurements, eliminating the need for a separate CT or MR anatomical scan. It achieves this by pairing a pre-trained diffusion model that operates in wavelet coefficient space with posterior sampling to address the joint inverse problem in three dimensions. A reader would care because this simplifies the imaging workflow, avoids misalignment artifacts, and lowers extra radiation exposure from auxiliary scans. On simulated data the approach produces higher-quality, noise-free results than standard maximum-likelihood joint reconstruction or U-Net post-processing, particularly when time-of-flight information is present.

Core claim

We propose a joint reconstruction of activity and attenuation approach that relies solely on emission data by combining a wavelet diffusion model and diffusion posterior sampling to produce fully three-dimensional reconstructions. Experimental results on simulated data show the method outperforms maximum likelihood activity and attenuation and MLAA-UNet, yielding high-quality noise-free images across count settings with time-of-flight; it also reconstructs non-TOF data with noticeable degradation in low-count conditions and demonstrates feasibility on real Biograph mMR data with joint scatter estimation.

What carries the argument

Wavelet diffusion model combined with diffusion posterior sampling applied to the joint activity-attenuation estimation problem from PET emission data.

If this is right

  • Stand-alone PET imaging becomes feasible without auxiliary anatomical scans while maintaining quantification accuracy when time-of-flight data is available.
  • High-quality noise-free reconstructions are obtained even in low-count regimes provided time-of-flight information is used.
  • The framework can process non-TOF data, although reconstruction quality drops markedly under low-count conditions.
  • Real-data results with simultaneous scatter estimation indicate readiness for clinical workflow integration.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If generalization holds, hybrid PET/CT systems could reduce patient dose by omitting the CT component in selected protocols.
  • The same posterior-sampling strategy in wavelet space might transfer to related tomographic inverse problems such as SPECT attenuation correction.
  • Large-scale testing across scanner vendors and patient populations would be the next practical step to establish robustness.

Load-bearing premise

The pre-trained wavelet diffusion model generalizes to unseen patient anatomies, scanner geometries, and count levels without significant domain shift or retraining.

What would settle it

Direct comparison of the estimated attenuation maps against ground-truth maps derived from co-registered CT scans on a collection of real clinical PET patient studies, checking whether the voxel-wise errors remain within clinical tolerance for accurate quantification.

Figures

Figures reproduced from arXiv: 2505.18782 by Alexandre Bousse, Antoine De Paepe, Baptiste Laurent, Catherine Cheze-Le-Rest, Cl\'ementine Phung-Ngoc, Dimitris Visvikis, Hong-Phuong Dang, Olivier Saut, Thibaut Merlin.

Figure 1
Figure 1. Figure 1: Proposed reconstruction framework for JRAA-DPS. A WDM is trained to generate activity-attenuation image pairs and [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Experiment 1—Reconstructions in HC setting with relative difference with respect to the reference images. are shown in [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Experiment 1—Reconstructions in LC setting with relative difference with respect to the reference images. (same here ) 0 0.2 0.4 0.6 0.8 1 SSIM (Activity) HC MLAA (TOF) MLAA-UNet (TOF) JRAA-DPS (TOF) JRAA-DPS (non-TOF) LC 20 30 40 0 0.2 0.4 0.6 0.8 1 PSNR SSIM (Attenuation) 20 30 40 PSNR [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Experiment 1—Normalized absolute bias against normalized STD for different values of the hyperparameters (cf. Section IV-D2). All reconstructions used TOF data [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Experiment 1—CRC against RC for a synthetic 15-mm￾diameter hot lesion added into the liver for different values of the hyperparameters (cf. Section IV-D2). The background was defined as a 70-mm-diameter spherical region of interest (ROI) surrounding the lesion. Mean reconstructions with JRAA-DPS are displayed against the GT reference. All reconstructions used TOF data. Method Reconstruction duration (min) … view at source ↗
Figure 7
Figure 7. Figure 7: Experiment 1—Mean, bias and standard deviation maps of 30 JRAA-DPS reconstructions of a single measurement. Activity λ PSNR=32.85 SSIM=0.906 PSNR=32.85 SSIM=0.906 PSNR=37.23 SSIM=0.953 PSNR=37.23 SSIM=0.953 PSNR=34.23 SSIM=0.896 PSNR=34.23 SSIM=0.896 PSNR=34.56 SSIM=0.911 PSNR=34.56 SSIM=0.911 0 2500 5000 7500 10000 12500 15000 17500 20000 Activity (Bq/mL) 40 20 0 20 40 Signed relative error (%) Attenuatio… view at source ↗
Figure 8
Figure 8. Figure 8: Experiment 2—Reconstructions from Biograph mMR (non-TOF). The bed attenuation was added in the forward model following (24). MLAA attenuation support was derived from the reference MR-derived attenuation map [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Experiment 2–Activity and attenuation profiles along a tumor (cf. green line in [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
read the original abstract

Attenuation correction (AC) is necessary for accurate activity quantification in positron emission tomography (PET). Conventional reconstruction methods typically rely on attenuation maps derived from a co-registered computed tomography (CT) or magnetic resonance (MR) scan. However, this additional scan may complicate the imaging workflow, introduce misalignment artifacts and increase radiation exposure. In this paper, we propose a joint reconstruction of activity and attenuation (JRAA) approach that eliminates the need for auxiliary anatomical imaging by relying solely on emission data. This framework combines wavelet diffusion model (WDM) and diffusion posterior sampling (DPS) to reconstruct fully three-dimensional (3-D) data. Experimental results on simulated data show our method outperforms maximum likelihood activity and attenuation (MLAA) and MLAA-UNet with U-Net-based post processing, and yields high-quality noise-free reconstructions across various count settings with time-of-flight (TOF). It is also able to reconstruct non-TOF data, although the reconstruction quality significantly degrades in low-count (LC) conditions, limiting its practical effectiveness in such settings. Nonetheless, a non-TOF Biograph mMR real data reconstruction with joint scatter estimation highlights the potential of the method for clinical applications. This approach represents a step towards stand-alone PET imaging by reducing the dependence on anatomical modalities while maintaining quantification accuracy, even in LC scenarios when TOF information is available. Our code is available on GitHub at https://github.com/clemphg/jraa-dps.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a joint reconstruction of activity and attenuation (JRAA) method for PET that relies solely on emission data by combining a pre-trained wavelet diffusion model (WDM) with diffusion posterior sampling (DPS). It reports that the approach outperforms MLAA and MLAA-UNet on simulated data, produces high-quality noise-free reconstructions across count levels when TOF information is available, and demonstrates feasibility on one real non-TOF Biograph mMR scan with joint scatter estimation, positioning it as a step toward stand-alone PET imaging without auxiliary anatomical scans.

Significance. If the performance claims are substantiated, the work would address a clinically relevant problem by reducing dependence on co-registered CT/MR scans, thereby simplifying workflows and lowering radiation exposure. The technical choice of operating diffusion models in wavelet coefficient space for 3-D joint reconstruction is novel for this application, and the public release of code supports reproducibility. However, the current evidence base is preliminary, relying primarily on qualitative comparisons and limited validation data.

major comments (2)
  1. [Abstract] Abstract: the central claim of outperformance over MLAA and MLAA-UNet is stated only in qualitative terms without reporting specific quantitative metrics (e.g., RMSE, SSIM, or bias values), statistical significance tests, error bars, or details on the number of simulated realizations, training data composition, or hyperparameter selection. This absence directly weakens the ability to assess the magnitude and reliability of the reported improvements.
  2. [Results] Results section (simulated and real data experiments): validation is confined to a small set of simulated phantoms plus a single real non-TOF scan, with no quantitative cross-anatomy, cross-scanner, or cross-count-level ablation studies. Because the WDM is pre-trained, any domain shift in wavelet statistics for unseen patient anatomies or scanner geometries (including TOF vs. non-TOF response) would propagate into the DPS-conditioned samples and could invalidate the generalization claim for stand-alone PET viability.
minor comments (2)
  1. [Abstract / Introduction] The abstract and introduction would benefit from a brief statement of the precise loss or objective used to train the WDM and how the DPS conditioning is formulated for the joint activity-attenuation problem.
  2. [Figures] Figure captions and axis labels should explicitly indicate whether reconstructions are shown for TOF or non-TOF data and the corresponding count level to improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which help clarify the presentation of our results. We agree that quantitative details strengthen the abstract and have revised it accordingly. For the validation concerns, we have expanded the results and discussion sections to include additional quantitative analyses and explicit limitations while maintaining the scope of this feasibility study.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim of outperformance over MLAA and MLAA-UNet is stated only in qualitative terms without reporting specific quantitative metrics (e.g., RMSE, SSIM, or bias values), statistical significance tests, error bars, or details on the number of simulated realizations, training data composition, or hyperparameter selection. This absence directly weakens the ability to assess the magnitude and reliability of the reported improvements.

    Authors: We agree that the abstract would benefit from quantitative support. In the revised manuscript we have added representative RMSE, SSIM and bias values for activity and attenuation maps (with standard deviations across 10 noise realizations), noted that paired t-tests showed statistically significant improvements (p < 0.05), and included brief statements on training data (50 simulated 3-D phantoms) and hyperparameter choices (guidance scale and number of DPS steps). revision: yes

  2. Referee: [Results] Results section (simulated and real data experiments): validation is confined to a small set of simulated phantoms plus a single real non-TOF scan, with no quantitative cross-anatomy, cross-scanner, or cross-count-level ablation studies. Because the WDM is pre-trained, any domain shift in wavelet statistics for unseen patient anatomies or scanner geometries (including TOF vs. non-TOF response) would propagate into the DPS-conditioned samples and could invalidate the generalization claim for stand-alone PET viability.

    Authors: We acknowledge the limited validation set. The revised results section now reports quantitative metrics across three count levels (high, medium, low) with error bars and an explicit ablation on count rate. We have added a dedicated limitations paragraph discussing potential domain shift in wavelet statistics for unseen anatomies or scanner geometries and have tempered the generalization language to emphasize that the work demonstrates feasibility rather than broad clinical readiness. Comprehensive cross-anatomy and cross-scanner experiments would require additional multi-center data beyond the present study. revision: partial

Circularity Check

0 steps flagged

No circularity in derivation; method uses independent pre-trained prior and separate validation data

full rationale

The paper presents a joint activity-attenuation reconstruction method that combines a pre-trained wavelet diffusion model with diffusion posterior sampling conditioned on emission data. Performance is evaluated via comparisons to MLAA and MLAA-UNet on simulated phantoms and one real non-TOF dataset, with no equations or claims reducing the reported outperformance to quantities fitted directly from the evaluation data. The central premise relies on the generalization of the pre-trained WDM, which is treated as an external component rather than derived from the test cases; validation phantoms are independent of any fitting steps described. No self-citation chains, self-definitional loops, or fitted-input-as-prediction patterns appear in the abstract or described workflow, rendering the derivation self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the assumption that a diffusion model pre-trained on wavelet coefficients of PET emission data can serve as a sufficiently accurate prior for posterior sampling of both activity and attenuation; no explicit free parameters or new physical entities are declared in the abstract.

axioms (1)
  • domain assumption A pre-trained wavelet diffusion model captures the statistical distribution of realistic PET emission data sufficiently well to act as a prior for joint reconstruction.
    Invoked when diffusion posterior sampling is used to recover activity and attenuation from emission measurements alone.

pith-pipeline@v0.9.0 · 5833 in / 1183 out tokens · 43423 ms · 2026-05-19T14:14:38.077877+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · 2 internal anchors

  1. [1]

    History and future technical innovation in positron emission tomography,

    T. Jones and D. Townsend, “History and future technical innovation in positron emission tomography,” Journal of Medical Imaging , vol. 4, no. 1, pp. 011 013– 011 013, 2017

  2. [2]

    Time-of-flight PET data determine the attenuation sinogram up to a constant,

    M. Defrise, A. Rezaei, and J. Nuyts, “Time-of-flight PET data determine the attenuation sinogram up to a constant,” Physics in Medicine & Biology , vol. 57, no. 4, p. 885, 2012

  3. [3]

    Simultaneous reconstruction of activity and attenuation in time-of-flight PET,

    A. Rezaei, M. Defrise, G. Bal, C. Michel, M. Conti, C. Watson, and J. Nuyts, “Simultaneous reconstruction of activity and attenuation in time-of-flight PET,” IEEE transactions on medical imaging , vol. 31, no. 12, pp. 2224–2233, 2012

  4. [4]

    ML-reconstruction for TOF-PET with simultaneous estimation of the attenuation factors,

    A. Rezaei, M. Defrise, and J. Nuyts, “ML-reconstruction for TOF-PET with simultaneous estimation of the attenuation factors,” IEEE transactions on medical imaging, vol. 33, no. 7, pp. 1563–1572, 2014

  5. [5]

    Attenuation correction in emission tomography using the emission data—a review,

    Y . Berker and Y . Li, “Attenuation correction in emission tomography using the emission data—a review,” Medical Physics , vol. 43, no. 2, pp. 807–832, 2016

  6. [6]

    Joint estimation of activity, attenuation and motion in respiratory-self-gated time-of-flight PET,

    M. Elhamiasl, F. Jolivet, A. Rezaei, M. Fieseler, K. Sch ¨afers, J. Nuyts, G. Schramm, and F. Boada, “Joint estimation of activity, attenuation and motion in respiratory-self-gated time-of-flight PET,” Physics in Medicine & Biology , vol. 70, no. 7, p. 075 003, 2025

  7. [7]

    Deep learning for PET image reconstruction,

    A. J. Reader, G. Corda, A. Mehranian, C. da Costa-Luis, S. Ellis, and J. A. Schnabel, “Deep learning for PET image reconstruction,” IEEE Transactions on Radiation and Plasma Medical Sciences , vol. 5, no. 1, pp. 1–25, 2020

  8. [8]

    A review on low-dose emission tomography post-reconstruction denoising with neural network approaches,

    A. Bousse, V . S. S. Kandarpa, K. Shi, K. Gong, J. S. Lee, C. Liu, and D. Visvikis, “A review on low-dose emission tomography post-reconstruction denoising with neural network approaches,” IEEE Transactions on Radiation and Plasma Medical Sciences , vol. 8, no. 4, pp. 333–347, 2024

  9. [9]

    Deep-learning-based methods of attenuation correction for SPECT and PET,

    X. Chen and C. Liu, “Deep-learning-based methods of attenuation correction for SPECT and PET,” Journal of Nuclear Cardiology , vol. 30, no. 5, pp. 1859–1878, 2023

  10. [10]

    Direct attenuation correction of brain PET images using only emission data via a deep convolutional encoder-decoder (Deep-DAC),

    I. Shiri, P. Ghafarian, P. Geramifar, K. H. -Y . Leung, M. Ghelichoghli, M. Oveisi, A. Rahmim, and M. R. Ay, “Direct attenuation correction of brain PET images using only emission data via a deep convolutional encoder-decoder (Deep-DAC),” European radiology, vol. 29, pp. 6867–6879, 2019

  11. [11]

    Deep learning-based attenuation correction in the absence of structural information for whole-body positron emission tomography imaging,

    X. Dong, Y . Lei, T. Wang, K. Higgins, T. Liu, W. J. Curran, H. Mao, J. A. Nye, and X. Yang, “Deep learning-based attenuation correction in the absence of structural information for whole-body positron emission tomography imaging,” Physics in Medicine & Biology , vol. 65, no. 5, p. 055 011, 2020

  12. [12]

    A deep learning approach for 18F-FDG PET attenuation correction,

    F. Liu, H. Jang, R. Kijowski, G. Zhao, T. Bradshaw, and A. B. McMillan, “A deep learning approach for 18F-FDG PET attenuation correction,” EJNMMI physics, vol. 5, pp. 1–15, 2018

  13. [13]

    Synthetic CT generation from non-attenuation corrected PET images for whole-body PET imaging,

    X. Dong, T. Wang, Y . Lei, K. Higgins, T. Liu, W. J. Curran, H. Mao, J. A. Nye, and X. Yang, “Synthetic CT generation from non-attenuation corrected PET images for whole-body PET imaging,” Physics in Medicine & Biology , vol. 64, no. 21, p. 215 016, 2019

  14. [14]

    Improving the accuracy of simultaneously reconstructed activity and attenuation maps using deep learning,

    D. Hwang, K. Y . Kim, S. K. Kang, S. Seo, J. C. Paeng, D. S. Lee, and J. S. Lee, “Improving the accuracy of simultaneously reconstructed activity and attenuation maps using deep learning,” Journal of Nuclear Medicine , vol. 59, no. 10, pp. 1624–1629, 2018

  15. [15]

    Diffusion models beat GANs on image synthesis,

    P. Dhariwal and A. Nichol, “Diffusion models beat GANs on image synthesis,” Advances in neural information processing systems , vol. 34, pp. 8780–8794, 2021

  16. [16]

    Diffusion models in medical imaging: A comprehensive survey,

    A. Kazerouni, E. K. Aghdam, M. Heidari, R. Azad, M. Fayyaz, I. Hacihaliloglu, and D. Merhof, “Diffusion models in medical imaging: A comprehensive survey,” Medical Image Analysis , p. 102 846, 2023

  17. [17]

    Score-based generative models for PET image reconstruction,

    I. R. Singh, A. Denker, R. Barbano, ˇZ. Kereta, B. Jin, K. Thielemans, P. Maass, and S. Arridge, “Score-based generative models for PET image reconstruction,” Machine Learning for Biomedical Imaging , vol. 2, pp. 547–585, Special Issue for Generative Models 2024, ISSN : 2766-905X

  18. [18]

    Likelihood-scheduled score-based generative modeling for fully 3D PET image reconstruction,

    G. Webber, Y . Mizuno, O. D. Howes, A. Hammers, A. P. King, and A. J. Reader, “Likelihood-scheduled score-based generative modeling for fully 3D PET image reconstruction,” IEEE transactions on medical imaging , 2025

  19. [19]

    CT reconstruction using diffusion posterior sampling conditioned on a nonlinear measurement model,

    S. Li, X. Jiang, M. Tivnan, G. J. Gang, Y . Shen, and J. W. Stayman, “CT reconstruction using diffusion posterior sampling conditioned on a nonlinear measurement model,” Journal of Medical Imaging , vol. 11, no. 4, pp. 043 504– 043 504, 2024

  20. [20]

    Material decomposition in photon-counting computed tomography with diffusion models: Comparative study and hybridization with variational regularizers,

    C. Vazia, T. Dassow, A. Bousse, J. Froment, B. Vedel, F. Vermet, A. Perelli, J. -P. Tasu, and D. Visvikis, “Material decomposition in photon-counting computed tomography with diffusion models: Comparative study and hybridization with variational regularizers,” arXiv preprint arXiv:2503.15383 , 2025

  21. [21]

    Joint reconstruction of the activity and the attenuation in PET by diffusion posterior sampling: A feasibility study,

    C. Phung-Ngoc, A. Bousse, A. De Paepe, H. -P. Dang, O. Saut, and D. Visvikis, “Joint reconstruction of the activity and the attenuation in PET by diffusion posterior sampling: A feasibility study,” arXiv preprint arXiv:2412.11776 , 2024

  22. [22]

    Joint reconstruction of activity and attenuation for PET imaging with diffusion prior,

    S. Bae, J. S. Lee, and K. Gong, “Joint reconstruction of activity and attenuation for PET imaging with diffusion prior,” in 2025 IEEE 22nd International Symposium on Biomedical Imaging (ISBI) , IEEE, 2025, pp. 1–4

  23. [23]

    WDM: 3D wavelet diffusion models for high-resolution medical image synthesis,

    P. Friedrich, J. Wolleb, F. Bieder, A. Durrer, and P. C. Cattin, “WDM: 3D wavelet diffusion models for high-resolution medical image synthesis,” in MICCAI Workshop on Deep Generative Models , Springer, 2024, pp. 11–21

  24. [24]

    A concave prior penal- izing relative differences for maximum-a-posteriori reconstruction in emission tomography,

    J. Nuyts, D. Bequ ´e, P. Dupont, and L. Mortelmans, “A concave prior penal- izing relative differences for maximum-a-posteriori reconstruction in emission tomography,” IEEE Transactions on nuclear science , vol. 49, no. 1, pp. 56–60, 2002

  25. [25]

    Maximum likelihood reconstruction for emission tomography,

    L. A. Shepp and Y . Vardi, “Maximum likelihood reconstruction for emission tomography,” IEEE transactions on medical imaging , vol. 1, no. 2, pp. 113–122, 1982

  26. [26]

    A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography,

    A. De Pierro, “A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography,” IEEE Transactions on Medical Imaging, vol. 14, no. 1, pp. 132–137, 1995

  27. [27]

    (An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods,

    S. R. Arridge, M. J. Ehrhardt, and K. Thielemans, “(An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods,” Philosophical Transactions of the Royal Society A , vol. 379, no. 2200, p. 20 200 205, 2021

  28. [28]

    Practical joint reconstruction of activity and attenuation with autonomous scaling for time-of-flight PET,

    Y . Li, S. Matej, and J. S. Karp, “Practical joint reconstruction of activity and attenuation with autonomous scaling for time-of-flight PET,” Physics in Medicine & Biology , vol. 65, no. 23, p. 235 037, 2020. 15 Activity λ Attenuation µ Reference No joint scatter refinement Scanner- derived PSNR=35.21 SSIM=0.90 PSNR=35.21 SSIM=0.90 PSNR=18.78 SSIM=0.66 P...

  29. [29]

    Denoising diffusion probabilistic models,

    J. Ho, A. Jain, and P. Abbeel, “Denoising diffusion probabilistic models,” Advances in neural information processing systems , vol. 33, pp. 6840–6851, 2020

  30. [30]

    Generative modeling by estimating gradients of the data distribution,

    Y . Song and S. Ermon, “Generative modeling by estimating gradients of the data distribution,” Advances in neural information processing systems , vol. 32, 2019

  31. [31]

    Denoising Diffusion Implicit Models

    J. Song, C. Meng, and S. Ermon, “Denoising diffusion implicit models,” arXiv preprint arXiv:2010.02502, 2020

  32. [32]

    Diffusion posterior sampling for synergistic recon- struction in spectral computed tomography,

    C. Vazia, A. Bousse, B. Vedel, F. Vermet, Z. Wang, T. Dassow, J. -P. Tasu, D. Visvikis, and J. Froment, “Diffusion posterior sampling for synergistic recon- struction in spectral computed tomography,” in 2024 IEEE 21st international symposium on biomedical imaging (ISBI 2024). IEEE , 2024

  33. [33]

    Diffusion posterior sampling for general noisy inverse problems,

    H. Chung, J. Kim, M. T. Mccann, M. L. Klasky, and J. C. Ye, “Diffusion posterior sampling for general noisy inverse problems,” in In Proceedings of the Eleventh International Conference on Learning Representations , ICLR, 2023

  34. [34]

    Decomposed diffusion sampler for accelerating large-scale inverse problems,

    H. Chung, S. Lee, and J. C. Ye, “Decomposed diffusion sampler for accelerating large-scale inverse problems,” in The Twelfth International Conference on Learning Representations, ICLR 2024, Vienna, Austria, May 7-11, 2024 , OpenReview.net, 2024

  35. [35]

    A Survey on Diffusion Models for Inverse Problems

    G. Daras, H. Chung, C.-H. Lai, Y . Mitsufuji, J. C. Ye, P. Milanfar, A. G. Dimakis, and M. Delbracio, “A survey on diffusion models for inverse problems,” arXiv preprint arXiv:2410.00083, 2024

  36. [36]

    arXiv preprint arXiv:2111.08005 , year=

    Y . Song, L. Shen, L. Xing, and S. Ermon, “Solving inverse problems in medical imaging with score-based generative models,” arXiv preprint arXiv:2111.08005 , 2021

  37. [37]

    PET image denoising based on 3D denoising diffusion probabilistic model: Evaluations on total-body datasets,

    B. Yu, S. Ozdemir, Y . Dong, W. Shao, K. Shi, and K. Gong, “PET image denoising based on 3D denoising diffusion probabilistic model: Evaluations on total-body datasets,” in International Conference on Medical Image Computing and Computer-Assisted Intervention , Springer, 2024, pp. 541–550

  38. [38]

    Solving inverse problems using pre-trained 2D diffusion models,

    H. Chung, D. Ryu, M. T. McCann, M. L. Klasky, and J. C. Ye, “Solving inverse problems using pre-trained 2D diffusion models,” in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition , 2023, pp. 22 542–22 551

  39. [39]

    Improving 3D imaging with pre-trained perpendicular 2D diffusion models,

    S. Lee, H. Chung, M. Park, J. Park, W. -S. Ryu, and J. C. Ye, “Improving 3D imaging with pre-trained perpendicular 2D diffusion models,” in Proceedings of the IEEE/CVF International Conference on Computer Vision , 2023, pp. 10 710– 10 720

  40. [40]

    PET image denoising based on denoising diffusion probabilistic model,

    K. Gong, K. Johnson, G. El Fakhri, Q. Li, and T. Pan, “PET image denoising based on denoising diffusion probabilistic model,” European Journal of Nuclear Medicine and Molecular Imaging , vol. 51, no. 2, pp. 358–368, 2024

  41. [41]

    DiffusionBlend: Learning 3D image prior through position-aware diffusion score blending for 3D computed tomography reconstruction,

    B. Song, J. Hu, Z. Luo, J. Fessler, and L. Shen, “DiffusionBlend: Learning 3D image prior through position-aware diffusion score blending for 3D computed tomography reconstruction,” Advances in Neural Information Processing Systems , vol. 37, pp. 89 584–89 611, 2024

  42. [42]

    Brain imaging generation with latent diffusion models,

    W. H. Pinaya, P.-D. Tudosiu, J. Dafflon, P. F. Da Costa, V . Fernandez, P. Nachev, S. Ourselin, and M. J. Cardoso, “Brain imaging generation with latent diffusion models,” in MICCAI Workshop on Deep Generative Models , Springer, 2022, pp. 117–126

  43. [43]

    Solving linear inverse problems provably via posterior sampling with latent diffusion models,

    L. Rout, N. Raoof, G. Daras, C. Caramanis, A. Dimakis, and S. Shakkottai, “Solving linear inverse problems provably via posterior sampling with latent diffusion models,” Advances in Neural Information Processing Systems , vol. 36, pp. 49 960–49 990, 2023

  44. [44]

    Solving blind inverse problems: Adaptive diffusion models for motion-corrected sparse-view 4DCT,

    A. De Paepe, A. Bousse, C. Phung-Ngoc, and D. Visvikis, “Solving blind inverse problems: Adaptive diffusion models for motion-corrected sparse-view 4DCT,” arXiv preprint arXiv:2501.12249 , 2025. 16

  45. [45]

    Adaptive diffusion models for sparse-view motion-corrected head cone-beam CT,

    A. De Paepe, A. Bousse, C. Phung-Ngoc, Y . Mellak, and D. Visvikis, “Adaptive diffusion models for sparse-view motion-corrected head cone-beam CT,” arXiv e-prints, arXiv–2504, 2025

  46. [46]

    Correction methods for random coincidences in fully 3D whole-body PET: Impact on data and image quality,

    D. Brasse, P. E. Kinahan, C. Lartizien, C. Comtat, M. Casey, and C. Michel, “Correction methods for random coincidences in fully 3D whole-body PET: Impact on data and image quality,” Journal of nuclear medicine , vol. 46, no. 5, pp. 859–867, 2005

  47. [47]

    New, faster, image-based scatter correction for 3D PET,

    C. C. Watson, “New, faster, image-based scatter correction for 3D PET,” IEEE Transactions on Nuclear Science , vol. 47, no. 4, pp. 1587–1594, 2000

  48. [48]

    Perfor- mance analysis of an improved 3-D PET Monte Carlo simulation and scatter correction,

    C. Holdsworth, C. Levin, M. Janecek, M. Dahlbom, and E. Hoffman, “Perfor- mance analysis of an improved 3-D PET Monte Carlo simulation and scatter correction,” IEEE transactions on nuclear science , vol. 49, no. 1, pp. 83–89, 2002

  49. [49]

    PET scatter estimation using deep learning U-Net architecture,

    B. Laurent, A. Bousse, T. Merlin, S. Nekolla, and D. Visvikis, “PET scatter estimation using deep learning U-Net architecture,” Physics in Medicine & Biology, vol. 68, no. 6, p. 065 004, 2023

  50. [50]

    Evaluation of deep learning-based scatter correction on a long-axial field-of- view PET scanner,

    B. Laurent, A. Bousse, T. Merlin, A. Rominger, K. Shi, and D. Visvikis, “Evaluation of deep learning-based scatter correction on a long-axial field-of- view PET scanner,” European Journal of Nuclear Medicine and Molecular Imaging, pp. 1–14, 2025

  51. [51]

    Energy-based scatter estimation in clinical PET,

    J. J. Hamill, J. Cabello, S. Surti, and J. S. Karp, “Energy-based scatter estimation in clinical PET,” Medical Physics , vol. 51, no. 1, pp. 54–69, 2024

  52. [52]

    Joint material decomposition and scatter estimation for spectral CT,

    A. Lorenzon, S. Z. Liu, X. Jiang, G. J. Gang, and J. W. Stayman, “Joint material decomposition and scatter estimation for spectral CT,” in Conference proceedings. International Conference on Image F ormation in X-Ray Computed Tomography, vol. 2024, 2024, p. 186

  53. [53]

    PARALLELPROJ—an open-source framework for fast calculation of projections in tomography,

    G. Schramm and K. Thielemans, “PARALLELPROJ—an open-source framework for fast calculation of projections in tomography,” Frontiers in Nuclear Medicine , vol. 3, p. 1 324 562, 2024

  54. [54]

    PET attenuation coefficients from CT images: Experimental evaluation of the transformation of CT into PET 511-kev attenuation coefficients,

    C. Burger, G. Goerres, S. Schoenes, A. Buck, A. Lonn, and G. V on Schulthess, “PET attenuation coefficients from CT images: Experimental evaluation of the transformation of CT into PET 511-kev attenuation coefficients,” European journal of nuclear medicine and molecular imaging , vol. 29, no. 7, pp. 922–927, 2002

  55. [55]

    CT-free attenuation correction of 13N-ammonia cardiac PET images using conditional denoising diffusion implicit model with logarithmic linear normalization,

    H. Sun, X. Hong, Y . Lu, F. Wang, A. Sanaat, W. Xu, S. Wang, H. Zaidi, and L. Lu, “CT-free attenuation correction of 13N-ammonia cardiac PET images using conditional denoising diffusion implicit model with logarithmic linear normalization,” Computer Methods and Programs in Biomedicine , p. 109 188, 2025

  56. [56]

    Deep learning-based attenuation correction for whole-body PET—a multi-tracer study with 18F-FDG, 68Ga- DOTATATE, and 18F-Fluciclovine,

    T. Toyonaga, D. Shao, L. Shi, J. Zhang, E. M. Revilla, D. Menard, J. Ankrah, K. Hirata, M. -K. Chen, J. A. Onofrey, et al. , “Deep learning-based attenuation correction for whole-body PET—a multi-tracer study with 18F-FDG, 68Ga- DOTATATE, and 18F-Fluciclovine,”European journal of nuclear medicine and molecular imaging , vol. 49, no. 9, pp. 3086–3097, 2022

  57. [57]

    Transmission imaging for integrated PET-MR systems,

    S. L. Bowen, N. Fuin, M. A. Levine, and C. Catana, “Transmission imaging for integrated PET-MR systems,” Physics in Medicine & Biology , vol. 61, no. 15, p. 5547, 2016

  58. [58]

    Improved PET/MRI accuracy by use of static transmission source in empirically derived hardware attenuation correction,

    A. Farag, R. T. Thompson, J. D. Thiessen, F. S. Prato, and J. Th´eberge, “Improved PET/MRI accuracy by use of static transmission source in empirically derived hardware attenuation correction,” EJNMMI physics , vol. 8, no. 1, p. 24, 2021

  59. [59]

    Attenuation correction for human PET/MRI studies,

    C. Catana, “Attenuation correction for human PET/MRI studies,” Physics in Medicine & Biology , vol. 65, no. 23, 23TR02, 2020

  60. [60]

    Impact of time-of-flight PET on quantification errors in MR imaging–based attenuation correction,

    A. Mehranian and H. Zaidi, “Impact of time-of-flight PET on quantification errors in MR imaging–based attenuation correction,” Journal of Nuclear Medicine , vol. 56, no. 4, pp. 635–641, 2015

  61. [61]

    Maximum-likelihood joint image reconstruction/motion estimation in attenuation-corrected respiratory gated PET/CT using a single attenuation map,

    A. Bousse, O. Bertolli, D. Atkinson, S. Arridge, S. Ourselin, B. F. Hutton, and K. Thielemans, “Maximum-likelihood joint image reconstruction/motion estimation in attenuation-corrected respiratory gated PET/CT using a single attenuation map,” IEEE Transactions on Medical Imaging , vol. 35, no. 1, pp. 217–228, 2016

  62. [62]

    Maximum-likelihood joint image reconstruction and motion estimation with misaligned attenuation in TOF-PET/CT,

    A. Bousse, O. Bertolli, D. Atkinson, S. Arridge, S. Ourselin, B. F. Hutton, and K. Thielemans, “Maximum-likelihood joint image reconstruction and motion estimation with misaligned attenuation in TOF-PET/CT,” Physics in Medicine & Biology , vol. 61, no. 3, pp. L11–19, 2016