pith. machine review for the scientific record. sign in

arxiv: 2604.24437 · v1 · submitted 2026-04-27 · ⚛️ physics.med-ph

Recognition: unknown

Weighted Directional Total Nuclear Variation for Joint Yttrium-90 PET/SPECT Reconstruction with CTAC-derived Guidance

D Deidda, D R McGowan, J Anton-Rodriguez, K Thielemans, S Arridge, S Porter

Authors on Pith no claims yet

Pith reviewed 2026-05-07 17:09 UTC · model grok-4.3

classification ⚛️ physics.med-ph
keywords Yttrium-90PET reconstructionSPECT reconstructionjoint reconstructiontotal nuclear variationCT guidanceimage regularisationpersonalised dosimetry
0
0 comments X

The pith

A weighted directional total nuclear variation regularizer improves recovery of small lesions in joint Y-90 PET/SPECT reconstruction by aligning edges with CT-derived directions.

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

The paper sets out to show that a new joint regularizer can exploit the complementary strengths of low-count Y-90 PET and higher-count but blurred bremsstrahlung SPECT by enforcing geometry-consistent edges across the two modalities. CT attenuation maps supply the directional information that tells the algorithm where activity changes should occur together, while a nuclear-norm penalty on the joint gradient matrix avoids forcing the intensities themselves to match. Data-driven normalisation handles the large difference in scale between the modalities. In phantom tests the method raises recovery coefficients relative to simpler directional and nuclear-norm penalties and beats sequential kernel reconstruction for the smallest spheres; in nine patients it raises tumour-to-background ratio at matched background activity. These gains matter because accurate post-treatment activity maps are required for personalised dosimetry after selective internal radiation therapy.

Core claim

The central claim is that the weighted directional total nuclear variation (w-dTNV) regulariser, formed by taking the nuclear norm of the dual-modality Jacobian whose columns are the CT-directed and modality-normalised gradients of the PET and SPECT images, produces higher recovery coefficients for small spheres than directional total variation or weighted total nuclear variation and higher tumour-to-background ratios than sequential hybrid kernel expectation maximisation at comparable background levels.

What carries the argument

The weighted directional total nuclear variation (w-dTNV), which penalises the nuclear norm of a 2-by-3 Jacobian matrix assembled from the CT-weighted gradients of the PET and SPECT images after data-driven per-voxel scaling to equalise modality ranges.

If this is right

  • In the NEMA IEC phantom with 20 bootstrapped noise realisations, w-dTNV raises recovery coefficients above those of dTV and w-TNV and exceeds SHKEM for the smallest spheres.
  • Across 45 lesions in nine post-SIRT patients, w-dTNV produces higher tumour-to-background ratios than SHKEM at comparable background activity.
  • The CT-guided coupling supplies a practical route to more stable quantitative Y-90 activity maps for personalised dosimetry.

Where Pith is reading between the lines

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

  • The same CT-derived directional weighting could be tested on other paired modalities where one channel has high spatial resolution but low counts and the other has high counts but spatial blur.
  • If the normalisation step proves robust across scanner vendors, the method could be deployed without scanner-specific tuning.
  • Motion mismatch between the CT and the emission scans would be expected to degrade the directional guidance and could be checked by introducing controlled misalignment in phantom data.

Load-bearing premise

That the directions taken from the CT attenuation map correctly indicate where the true activity edges should align in both PET and SPECT, and that normalising the two modalities from separate preliminary reconstructions does not bias the final joint solution.

What would settle it

A larger patient study in which tumour-to-background ratio no longer increases under w-dTNV while background activity remains comparable, or a phantom experiment in which recovery coefficients for the smallest spheres fall below those of SHKEM, would falsify the claimed benefit.

Figures

Figures reproduced from arXiv: 2604.24437 by D Deidda, D R McGowan, J Anton-Rodriguez, K Thielemans, S Arridge, S Porter.

Figure 1
Figure 1. Figure 1: Left: VOIs (warm background (BG), lung insert (LI), view at source ↗
Figure 2
Figure 2. Figure 2: Visualisation of the five identified lesions and the view at source ↗
Figure 3
Figure 3. Figure 3: Visual comparison of NEMA phantom reconstructions view at source ↗
Figure 4
Figure 4. Figure 4: Recovery coefficient (RC) for NEMA phantom spheres for w-dTNV, dTV, w-TNV and SHKEM (shown at subiteration 27). Also shown are the equivalent RCs for a high-count, long-OSEM reconstruction (24ss, 240it). The grey line demonstrates perfect recovery. Statistical significance is reported for paired comparisons versus w-dTNV (Benjamini–Hochberg adjusted): ∗ p < 0.05, ∗∗ p < 0.01, ∗∗∗ p < 0.001. In contrast, SH… view at source ↗
Figure 5
Figure 5. Figure 5: Mean voxel-wise coefficient of variation (CoV) view at source ↗
Figure 6
Figure 6. Figure 6: Impact of varying PET weighting (ωPET) while holding SPECT regularisation constant (top row) and varying SPECT weighting (ωSPECT) while holding PET regularisation constant (bottom row) on TBR and CoVbkg for lesion 1 (first and second columns) and lesion 2 (third and fourth columns). The first and third columns show w-dTNV and the second and fourth show w-TNV. Axial views. Coronal views view at source ↗
Figure 7
Figure 7. Figure 7: Visual comparison of reconstruction methods for patient SIRT3 lesion 1 (top row in each quadrant) and lesion 2 view at source ↗
Figure 8
Figure 8. Figure 8: TBR mean (left) for all lesions and patients, and view at source ↗
Figure 9
Figure 9. Figure 9: Patient study. CTAC-derived attenuation map view at source ↗
Figure 10
Figure 10. Figure 10: NEMA phantom (central axial slice). CTAC-derived view at source ↗
read the original abstract

Quantitative post-treatment activity imaging is essential for personalised dosimetry after Yttrium-90 selective internal radiation therapy (SIRT). Yttrium-90 PET offers high spatial resolution but is extremely low-count, whereas bremsstrahlung SPECT has higher count statistics but is degraded by blur, scatter, and septal penetration. Since both modalities image the same microsphere distribution, joint synergistic reconstruction can exploit their physical coupling, with CT attenuation correction (CTAC) providing additional anatomical guidance. We propose weighted directional total nuclear variation (w-dTNV), a joint variational regulariser for coupled PET/SPECT reconstruction with CTAC-guided anisotropy. w-dTNV penalises the nuclear norm of a dual-modality Jacobian, promoting co-located, geometry-consistent edges without forcing intensity correlation. Directionality is derived from the CTAC attenuation map $\mu$ and applied to PET/SPECT gradients, allowing efficient per-voxel spectral computations. PET/SPECT scale disparity is mitigated using data-driven modality normalisation from preliminary reconstructions. We evaluate w-dTNV on a NEMA IEC phantom with 20 bootstrapped PET noise realisations and on 9 post-SIRT patients with 45 lesions, against dTV, w-TNV, and sequential hybrid kernel expectation maximisation (SHKEM). In the phantom, w-dTNV improves recovery coefficients over dTV and w-TNV, and improves recovery over SHKEM for the smallest spheres. In patients, w-dTNV gives higher tumour-to-background ratios than SHKEM at comparable background activity. These results suggest that CTAC-guided synergistic variational coupling improves lesion recovery and clinical lesion contrast, offering a practical route towards more stable post-SIRT Yttrium-90 activity estimates for personalised dosimetry.

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

1 major / 1 minor

Summary. The manuscript proposes weighted directional total nuclear variation (w-dTNV), a joint variational regularizer that applies a weighted nuclear norm to the dual-modality Jacobian of PET and SPECT gradients, with anisotropy directions derived from the CT attenuation map and scale disparity addressed via data-driven normalization from preliminary reconstructions. It reports improved recovery coefficients over dTV, w-TNV, and SHKEM on a NEMA IEC phantom (20 noise realizations) and higher tumour-to-background ratios than SHKEM at comparable background activity on 9 post-SIRT patients (45 lesions).

Significance. If the quantitative gains hold after addressing methodological details, the approach could improve quantitative accuracy for post-SIRT Y-90 dosimetry by exploiting complementary PET resolution and SPECT count statistics under CTAC anatomical guidance, without forcing intensity correlation between modalities. The dual-modality testing on phantom and patient data with multiple baselines is a positive aspect.

major comments (1)
  1. Abstract and Methods (normalization paragraph): The data-driven modality normalisation from preliminary reconstructions is load-bearing for the central claim, as the nuclear norm couples modalities through joint singular values; any residual scale mismatch or artifact bias from the unspecified preliminary reconstructions can directly influence which edges are promoted, making it unclear whether reported gains in recovery coefficients and TBR are attributable to w-dTNV/CTAC guidance or to the normalisation step itself. An explicit description of the preliminary reconstruction method, normalisation factor computation, and a sensitivity test would be required to support the attribution.
minor comments (1)
  1. Abstract: The quantitative improvements are stated qualitatively (e.g., 'improves recovery coefficients'); adding specific values or effect sizes with statistical details would aid assessment of clinical relevance.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and positive assessment of the potential impact of w-dTNV. We address the single major comment below.

read point-by-point responses

  1. Referee: Abstract and Methods (normalization paragraph): The data-driven modality normalisation from preliminary reconstructions is load-bearing for the central claim, as the nuclear norm couples modalities through joint singular values; any residual scale mismatch or artifact bias from the unspecified preliminary reconstructions can directly influence which edges are promoted, making it unclear whether reported gains in recovery coefficients and TBR are attributable to w-dTNV/CTAC guidance or to the normalisation step itself. An explicit description of the preliminary reconstruction method, normalisation factor computation, and a sensitivity test would be required to support the attribution.

    Authors: We agree that the data-driven normalization requires fuller documentation to support attribution of the reported gains. In the revised manuscript we will expand the Methods section with: (i) the exact preliminary reconstruction method (sequential hybrid-kernel EM applied independently to each modality), (ii) the normalization-factor computation (ratio of mean activity concentrations inside a liver VOI drawn on the preliminary images), and (iii) a sensitivity study showing that recovery coefficients and TBR remain stable for normalization factors varied by ±20 % around the data-driven value. These additions will make clear that the quantitative improvements arise from the CTAC-guided joint nuclear-norm coupling rather than from the scaling step alone. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation and evaluation remain self-contained

full rationale

The paper defines w-dTNV as a new joint variational regulariser that penalises the nuclear norm of the dual-modality Jacobian with CTAC-derived directional weights, then mitigates scale disparity via a separate data-driven normalisation step from preliminary reconstructions. The central claims consist of measured empirical gains in recovery coefficients (phantom) and tumour-to-background ratios (patients) against independent baselines (dTV, w-TNV, SHKEM). No equation or step equates the reported improvements to the normalisation factors or to any fitted parameter by construction; the regulariser is motivated from first principles and tested on external datasets. No self-citation chain, ansatz smuggling, or renaming of known results is load-bearing for the claimed performance gains.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Full text unavailable; specific free parameters, axioms, and invented entities cannot be audited in detail. The abstract implies data-driven normalization factors and CT-derived directions but provides no explicit list or values.

pith-pipeline@v0.9.0 · 5645 in / 1166 out tokens · 69826 ms · 2026-05-07T17:09:24.863998+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

57 extracted references · 29 canonical work pages

  1. [1]

    International recommendations for personalised selective internal radiation therapy of primary and metastatic liver diseases with yttrium-90 resin microspheres,

    H. Levillain, O. Bagni, C. M. Deroose, A. Dieudonn ´e, S. Gnesin, O. S. Grosser, S. C. Kappadath, A. Kennedy, N. Kokabi, D. M. Liu, D. C. Madoff, A. Mahvash, A. Martinez de la Cuesta, D. C. E. Ng, P. M. Paprottka, C. Pettinato, M. Rodr ´ıguez-Fraile, R. Salem, B. Sangro, L. Strigari, D. Y . Sze, B. J. de Wit van der veen, and P. Flamen, “International rec...

    work page doi:10.1007/s00259-020-05163-5 2021

  2. [2]

    Haste, M

    P. Haste, M. Tann, S. Persohn, T. LaRoche, V . Aaron, T. Maux- ion, N. Chauhan, M. R. Dreher, and M. S. Johnson, “Correlation of Technetium-99m Macroaggregated Albumin and Yttrium-90 Glass Mi- crosphere Biodistribution in Hepatocellular Carcinoma: A Retrospective Review of Pretreatment Single Photon Emission CT and Posttreatment Positron Emission Tomograp...

    2017

  3. [3]

    A new internal pair production branching ratio of 90Y: The development of a non-destructive assay for 90Y and 90Sr,

    R. G. Selwyn, R. J. Nickles, B. R. Thomadsen, L. A. DeWerd, and J. A. Micka, “A new internal pair production branching ratio of 90Y: The development of a non-destructive assay for 90Y and 90Sr,”Applied Radiation and Isotopes, vol. 65, no. 3, pp. 318–327, Mar. 2007

    2007

  4. [4]

    Development and evaluation of an improved quantitative (90)Y bremsstrahlung SPECT method,

    X. Rong, Y . Du, M. Ljungberg, E. Rault, S. Vandenberghe, and E. C. Frey, “Development and evaluation of an improved quantitative (90)Y bremsstrahlung SPECT method,”Medical Physics, vol. 39, no. 5, pp. 2346–2358, May 2012. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES, VOL. XX, NO. X, FEBRUARY 2026 13 Axial slice. Coronal slice. Fig. 9: Patie...

    2012

  5. [5]

    Deidda, A

    D. Deidda, A. M. Denis-Bacelar, A. J. Fenwick, K. M. Ferreira, W. Heetun, B. F. Hutton, D. R. McGowan, A. P. Robinson, J. Scuffham, K. Thielemans, and R. Twyman, “Triple modality image reconstruction of PET data using SPECT, PET, CT information increases lesion uptake in images of patients treated with radioembolization with 90Y micro-spheres,”EJNMMI Phys...

    work page doi:10.1186/s40658-023-00549-4 2023

  6. [6]

    Hybrid PET-MR list-mode kernelized expectation maximization reconstruction,

    D. Deidda, N. A. Karakatsanis, P. M. Robson, Y .-J. Tsai, N. Efthimiou, K. Thielemans, Z. A. Fayad, R. G. Aykroyd, and C. Tsoumpas, “Hybrid PET-MR list-mode kernelized expectation maximization reconstruction,” Inverse Problems, vol. 35, no. 4, p. 044001, Apr. 2019. [Online]. Available: https://iopscience.iop.org/article/10.1088/1361-6420/ab013f

    work page doi:10.1088/1361-6420/ab013f 2019

  7. [7]

    Joint MR-PET Reconstruction Using a Multi-Channel Image Regularizer,

    F. Knoll, M. Holler, T. Koesters, R. Otazo, K. Bredies, and D. K. Sod- ickson, “Joint MR-PET Reconstruction Using a Multi-Channel Image Regularizer,”IEEE Transactions on Medical Imaging, vol. 36, no. 1, pp. 1–16, Jan. 2017, conference Name: IEEE Transactions on Medical Imaging

    2017

  8. [8]

    (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, Jun. 2021. [Online]. Available: https://royalsocietypublishing .org/doi/abs/10.1098/rsta.2020.0205

    work page doi:10.1098/rsta.2020.0205 2021

  9. [9]

    Evaluation of quantitative 90Y SPECT based on experimental phantom studies,

    D. Minarik, K. S. Gleisner, and M. Ljungberg, “Evaluation of quantitative 90Y SPECT based on experimental phantom studies,” Physics in Medicine & Biology, vol. 53, no. 20, p. 5689, Sep. 2008. [Online]. Available: https://dx.doi.org/10.1088/0031-9155/53/20/008

    work page doi:10.1088/0031-9155/53/20/008 2008

  10. [10]

    P., LAIRD, N

    A. P. Dempster, N. M. Laird, and D. B. Rubin, “Maximum Likelihood from Incomplete Data Via the EM Algorithm,”Journal of the Royal Statistical Society: Series B (Methodological), vol. 39, no. 1, pp. 1–22, Sep. 1977. [Online]. Available: https://onlinelibrary.wiley.com/doi/10.1 111/j.2517-6161.1977.tb01600.x

    work page arXiv 1977

  11. [11]

    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

    1982

  12. [12]

    Accelerated image reconstruction using ordered subsets of projection data,

    H. Hudson and R. Larkin, “Accelerated image reconstruction using ordered subsets of projection data,”IEEE Transactions on Medical Imaging, vol. 13, no. 4, pp. 601–609, Dec. 1994, conference Name: IEEE Transactions on Medical Imaging

    1994

  13. [13]

    Stopping Rule for the MLE Algorithm Based on Statistical Hypothesis Testing,

    E. Veklerov and J. Llacer, “Stopping Rule for the MLE Algorithm Based on Statistical Hypothesis Testing,”IEEE Transactions on Medical Imaging, vol. 6, no. 4, pp. 313–319, Dec. 1987. [Online]. Available: https://ieeexplore.ieee.org/abstract/document/4307849

    work page arXiv 1987

  14. [14]

    Influence of Reconstruction Iterations on 18F-FDG PET/CT Standardized Uptake Values,

    C. J. Jaskowiak, J. A. Bianco, S. B. Perlman, and J. P. Fine, “Influence of Reconstruction Iterations on 18F-FDG PET/CT Standardized Uptake Values,”Journal of Nuclear Medicine, vol. 46, no. 3, pp. 424–428, Mar

  15. [15]

    Available: https://jnm.snmjournals.org/content/46/3/424

    [Online]. Available: https://jnm.snmjournals.org/content/46/3/424

  16. [16]

    Nonlinear total variation based noise removal algorithms,

    L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,”Physica D: Nonlinear Phenomena, vol. 60, no. 1, pp. 259–268, Nov. 1992. [Online]. Available: https://www.sciencedirect.com/science/article/pii/016727899290242F

    work page arXiv 1992

  17. [17]

    Total-Variation Regularization in Positron Emission Tomography,

    E. Jonsson, S.-c. Huang, and T. Chan, “Total-Variation Regularization in Positron Emission Tomography,” 1998

    1998

  18. [18]

    Multi-modality imaging with structure-promoting regularisers,

    M. J. Ehrhardt, “Multi-modality imaging with structure-promoting regularisers,”arXiv:2007.11689 [cs, eess, math], 2020. [Online]. Available: http://arxiv.org/abs/2007.11689

    work page arXiv 2007

  19. [19]

    Utilizing MRI information to estimate F18-FDG distributions in rat flank tumors,

    J. Bowsher, H. Yuan, L. Hedlund, T. Turkington, G. Akabani, A. Badea, W. Kurylo, C. Wheeler, G. Cofer, M. Dewhirst, and G. Johnson, “Utilizing MRI information to estimate F18-FDG distributions in rat flank tumors,” inIEEE Symposium Conference Record Nuclear Science 2004., vol. 4, Oct. 2004, pp. 2488–2492 V ol. 4

    2004

  20. [20]

    PET Reconstruction With an Anatomical MRI Prior Using Parallel Level Sets,

    M. J. Ehrhardt, P. Markiewicz, M. Liljeroth, A. Barnes, V . Kolehmainen, J. S. Duncan, L. Pizarro, D. Atkinson, B. F. Hutton, S. Ourselin, K. Thielemans, and S. R. Arridge, “PET Reconstruction With an Anatomical MRI Prior Using Parallel Level Sets,”IEEE Transactions on Medical Imaging, vol. 35, no. 9, pp. 2189–2199, Sep. 2016, conference Name: IEEE Transa...

    2016

  21. [21]

    PET Image Reconstruction Using Kernel Method,

    G. Wang and J. Qi, “PET Image Reconstruction Using Kernel Method,” IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES, VOL. XX, NO. X, FEBRUARY 2026 14 IEEE Transactions on Medical Imaging, vol. 34, no. 1, pp. 61–71, Jan. 2015, conference Name: IEEE Transactions on Medical Imaging

    2026

  22. [22]

    Anatomically-aided PET reconstruction using the kernel method,

    W. Hutchcroft, G. Wang, K. T. Chen, C. Catana, and J. Qi, “Anatomically-aided PET reconstruction using the kernel method,” Physics in Medicine and Biology, vol. 61, no. 18, pp. 6668–6683, Sep

  23. [23]

    Available: https://iopscience.iop.org/article/10.1088/003 1-9155/61/18/6668

    [Online]. Available: https://iopscience.iop.org/article/10.1088/003 1-9155/61/18/6668

    work page doi:10.1088/003

  24. [24]

    Hybrid kernelised expectation maximisation for Bremsstrahlung SPECT reconstruction in SIRT with 90Y micro- spheres,

    D. Deidda, A. M. Denis-Bacelar, A. J. Fenwick, K. M. Ferreira, W. Heetun, B. F. Hutton, A. P. Robinson, J. Scuffham, and K. Thielemans, “Hybrid kernelised expectation maximisation for Bremsstrahlung SPECT reconstruction in SIRT with 90Y micro- spheres,”EJNMMI Physics, vol. 9, no. 1, p. 25, Apr. 2022. [Online]. Available: https://doi.org/10.1186/s40658-022-00452-4

    work page doi:10.1186/s40658-022-00452-4 2022

  25. [25]

    Theranostic SPECT reconstruction for improved resolution: application to radionuclide therapy dosimetry,

    H. Marquis, D. Deidda, A. Gillman, K. P. Willowson, Y . Gholami, T. Hioki, E. Eslick, K. Thielemans, and D. L. Bailey, “Theranostic SPECT reconstruction for improved resolution: application to radionuclide therapy dosimetry,”EJNMMI Physics, vol. 8, no. 1, p. 16, Dec. 2021. [Online]. Available: https: //ejnmmiphys.springeropen.com/articles/10.1186/s40658-0...

    work page doi:10.1186/s40658-021-00362-x 2021

  26. [26]

    Joint reconstruction of PET-MRI by exploiting structural similarity,

    M. Ehrhardt, K. Thielemans, L. Pizarro, D. Atkinson, S. Ourselin, B. Hutton, and S. Arridge, “Joint reconstruction of PET-MRI by exploiting structural similarity,”Inverse Problems, vol. 31, p. 015001, Jan. 2015

    2015

  27. [27]

    Non-convex joint-sparsity regularization for synergistic PET and SENSE MRI reconstruction,

    A. Mehranian and A. Reader, “Non-convex joint-sparsity regularization for synergistic PET and SENSE MRI reconstruction,”Journal of Nuclear Medicine, vol. 57, no. supplement 2, pp. 639–639, May 2016. [Online]. Available: https://jnm.snmjournals.org/content/57/supplement 2/639

    2016

  28. [28]

    Joint Reconstruction of Multi-channel, Spectral CT Data via Constrained Total Nuclear Variation Minimization,

    D. S. Rigie and P. J. La Rivi `ere, “Joint Reconstruction of Multi-channel, Spectral CT Data via Constrained Total Nuclear Variation Minimization,”Physics in medicine and biology, vol. 60, no. 5, pp. 1741–1762, Mar. 2015. [Online]. Available: https: //www.ncbi.nlm.nih.gov/pmc/articles/PMC4669200/

    2015

  29. [29]

    Total Nuclear Variation and Jacobian Extensions of Total Variation for Vector Fields,

    K. M. Holt, “Total Nuclear Variation and Jacobian Extensions of Total Variation for Vector Fields,”IEEE Transactions on Image Processing, vol. 23, no. 9, pp. 3975–3989, Sep. 2014, conference Name: IEEE Transactions on Image Processing

    2014

  30. [30]

    Derivatives of Spectral Functions,

    A. S. Lewis, “Derivatives of Spectral Functions,”Mathematics of Operations Research, vol. 21, no. 3, pp. 576–588, Aug. 1996. [Online]. Available: https://pubsonline.informs.org/doi/10.1287/moor.21.3.576

    work page doi:10.1287/moor.21.3.576 1996

  31. [31]

    Radioembolization and the Dynamic Role of 90Y PET/CT,

    A. S. Pasciak, A. C. Bourgeois, J. M. McKinney, T. T. Chang, D. R. Osborne, S. N. Acuff, and Y . C. Bradley, “Radioembolization and the Dynamic Role of 90Y PET/CT,”Frontiers in Oncology, vol. 4, Feb

  32. [32]

    Available: https://www.frontiersin.org/journals/oncolog y/articles/10.3389/fonc.2014.00038/full

    [Online]. Available: https://www.frontiersin.org/journals/oncolog y/articles/10.3389/fonc.2014.00038/full

    work page doi:10.3389/fonc.2014.00038/full 2014

  33. [33]

    Accelerating Stochastic Gradient Descent using Predictive Variance Reduction,

    R. Johnson and T. Zhang, “Accelerating Stochastic Gradient Descent using Predictive Variance Reduction,” inAdvances in Neural Information Processing Systems, vol. 26. Curran Associates, Inc., 2013

    2013

  34. [34]

    Optimising Subset Selection in Synergistic Emission Tomography Reconstruction,

    S. D. Porter, D. Deidda, S. Arridge, and K. Thielemans, “Optimising Subset Selection in Synergistic Emission Tomography Reconstruction,” in2024 IEEE Nuclear Science Symposium (NSS), Medical Imaging Conference (MIC) and Room Temperature Semiconductor Detector Conference (RTSD), Oct. 2024, pp. 1–2, iSSN: 2577-0829. [Online]. Available: https://ieeexplore.ie...

    work page arXiv 2024

  35. [35]

    Fast EM-like methods for maximum

    A. De Pierro and M. Yamagishi, “Fast EM-like methods for maximum ”a posteriori” estimates in emission tomography,”IEEE Transactions on Medical Imaging, vol. 20, no. 4, pp. 280–288, Apr. 2001. [Online]. Available: http://ieeexplore.ieee.org/document/921477/

    2001

  36. [36]

    Globally convergent ordered subsets algorithms: application to tomography,

    S. Ahn and J. A. Fessler, “Globally convergent ordered subsets algorithms: application to tomography,” in2001 IEEE Nuclear Science Symposium Conference Record (Cat. No.01CH37310). IEEE, 2002, pp. 1064–1068. [Online]. Available: http://ieeexplore.ieee.org/docume nt/1009736/

    work page arXiv 2002

  37. [37]

    Fast PET Reconstruction with Variance Reduction and Prior-Aware Preconditioning,

    M. J. Ehrhardt, Z. Kereta, and G. Schramm, “Fast PET Reconstruction with Variance Reduction and Prior-Aware Preconditioning,” Jun. 2025, number: arXiv:2506.04976 arXiv:2506.04976 [physics]. [Online]. Available: http://arxiv.org/abs/2506.04976

    work page arXiv 2025

  38. [38]

    Assessment of bootstrap resampling performance for PET data,

    P. J. Markiewicz, A. J. Reader, and J. C. Matthews, “Assessment of bootstrap resampling performance for PET data,”Physics in Medicine & Biology, vol. 60, no. 1, p. 279, Dec. 2014. [Online]. Available: https://iopscience.iop.org/article/10.1088/0031-9155/60/1/279

    work page doi:10.1088/0031-9155/60/1/279 2014

  39. [39]

    Optimization of image reconstruction for yttrium-90 SIRT on a LYSO PET/CT system using a Bayesian penalized likelihood reconstruction algorithm,

    L. M. Rowley, K. M. Bradley, P. Boardman, A. Hallam, and D. R. McGowan, “Optimization of image reconstruction for yttrium-90 SIRT on a LYSO PET/CT system using a Bayesian penalized likelihood reconstruction algorithm,”Journal of Nuclear Medicine, Sep. 2016. [Online]. Available: https://jnm.snmjournals.org/content/early/2016/09/ 28/jnumed.116.176552

    2016

  40. [40]

    Random coincidence estimation from single event rates on the Discovery ST PET/CT scanner,

    C. Stearns, D. McDaniel, S. Kohlmyer, P. Arul, B. Geiser, and V . Shanmugam, “Random coincidence estimation from single event rates on the Discovery ST PET/CT scanner,” in2003 IEEE Nuclear Science Symposium. Conference Record, vol. 5, Oct. 2003, pp. 3067–3069 V ol.5. [Online]. Available: https://ieeexplore.ieee.org/document/1352545

    work page arXiv 2003

  41. [41]

    3D implementation of Scatter Estimation in 3D PET,

    M. Iatrou, R. M. Manjeshwar, S. G. Ross, K. Thielemans, and C. W. Stearns, “3D implementation of Scatter Estimation in 3D PET,” in2006 IEEE Nuclear Science Symposium Conference Record, vol. 4, Oct. 2006, pp. 2142–2145. [Online]. Available: https://ieeexplore.ieee.org/document/4179452

    work page arXiv 2006

  42. [42]

    Simultaneous Multi-Bed MAP Reconstruction with CT-Guided Directional TV Prior for Y-90 PET SIRT,

    S. Porter, D. Deidda, D. McGowan, S. Arridge, and K. Thielemans, “Simultaneous Multi-Bed MAP Reconstruction with CT-Guided Directional TV Prior for Y-90 PET SIRT,”The PET is Wonderful Journal, vol. 2, no. 1, Oct. 2025. [Online]. Available: https: //journals.ed.ac.uk/PiWJournal/article/view/10841

    2025

  43. [43]

    Sensitivity, resolution, and linearity of the scintillation camera,

    H. O. Anger, “Sensitivity, resolution, and linearity of the scintillation camera,”IEEE Transactions on Nuclear Science, vol. NS-13, no. 3, pp. 380–392, 1966

    1966

  44. [44]

    Fast simulation of yttrium-90 bremsstrahlung photons with GATE,

    E. Rault, S. Staelens, R. Van Holen, J. De Beenhouwer, and S. Vandenberghe, “Fast simulation of yttrium-90 bremsstrahlung photons with GATE,”Medical Physics, vol. 37, no. 6Part1, pp. 2943–2950, 2010, eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1118/1.3431998. [On- line]. Available: https://onlinelibrary.wiley.com/doi/abs/10.1118/1.3431 998

    work page doi:10.1118/1.3431998 2010

  45. [45]

    SIMIND Monte Carlo based image reconstruction,

    M. Ljungberg, “SIMIND Monte Carlo based image reconstruction,” Journal of Nuclear Medicine, vol. 56, no. supplement 3, pp. 43–43, May 2015. [Online]. Available: https://jnm.snmjournals.org/content/56 /supplement 3/43

    2015

  46. [46]

    A deep neural network for fast and accurate scatter estimation in quantitative SPECT/CT under challenging scatter conditions,

    H. Xiang, H. Lim, J. A. Fessler, and Y . K. Dewaraja, “A deep neural network for fast and accurate scatter estimation in quantitative SPECT/CT under challenging scatter conditions,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 47, no. 13, pp. 2956–2967, Dec. 2020. [Online]. Available: https://doi.org/10.1007/s00259-020-04840-9

    work page doi:10.1007/s00259-020-04840-9 2020

  47. [47]

    TotalSegmentator: Robust Segmentation of 104 Anatomic Structures in CT Images,

    J. Wasserthal, H.-C. Breit, M. T. Meyer, M. Pradella, D. Hinck, A. W. Sauter, T. Heye, D. T. Boll, J. Cyriac, S. Yang, M. Bach, and M. Segeroth, “TotalSegmentator: Robust Segmentation of 104 Anatomic Structures in CT Images,”Radiology: Artificial Intelligence, vol. 5, no. 5, p. e230024, Jul. 2023. [Online]. Available: https: //pmc.ncbi.nlm.nih.gov/article...

    2023

  48. [48]

    Controlling the false discovery rate: a practical and powerful approach to multiple testing.Journal of the Royal Statistical Society: Series B (Methodological), 57(1):289–300, 1995

    Y . Benjamini and Y . Hochberg, “Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing,” Journal of the Royal Statistical Society. Series B (Methodological), vol. 57, no. 1, pp. 289–300, 1995. [Online]. Available: https: //www.jstor.org/stable/2346101

    work page arXiv 1995

  49. [49]

    Random-effects models for longitudinal data,

    N. M. Laird and J. H. Ware, “Random-effects models for longitudinal data,”Biometrics, vol. 38, no. 4, pp. 963–974, Dec. 1982

    1982

  50. [50]

    Fitting Linear Mixed-Effects Models,

    “Fitting Linear Mixed-Effects Models,” inMixed-Effects Models in S and S-PLUS, J. C. Pinheiro and D. M. Bates, Eds. New York, NY: Springer, 2000, pp. 133–199. [Online]. Available: https://doi.org/10.1007/0-387-22747-4 4

    work page doi:10.1007/0-387-22747-4 2000

  51. [51]

    Impact of 90Y PET gradient-based tumor segmentation on voxel-level dosimetry in liver radioembolization,

    J. K. Mikell, R. K. Kaza, P. L. Roberson, K. C. Younge, R. N. Srinivasa, B. S. Majdalany, K. C. Cuneo, D. Owen, T. Devasia, M. J. Schipper, and Y . K. Dewaraja, “Impact of 90Y PET gradient-based tumor segmentation on voxel-level dosimetry in liver radioembolization,” EJNMMI Physics, vol. 5, no. 1, p. 31, Nov. 2018. [Online]. Available: https://doi.org/10....

    work page doi:10.1186/s40658-018-0230-y 2018

  52. [52]

    SIRF: Synergistic Image Reconstruction Framework,

    E. Ovtchinnikov, D. Atkinson, C. Kolbitsch, B. A. Thomas, O. Bertolli, C. O. da Costa-Luis, N. Efthimiou, R. Fowler, E. Pasca, P. Wadhwa, E. Emond, J. Matthews, C. Prieto, A. J. Reader, C. Tsoumpas, M. Turner, and K. Thielemans, “SIRF: Synergistic Image Reconstruction Framework,” in2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MI...

    work page arXiv 2017

  53. [53]

    STIR: Software for tomographic image reconstruction release 2,

    K. Thielemans, C. Tsoumpas, S. Mustafovic, T. Beisel, P. Aguiar, N. Dikaios, and M. Jacobson, “STIR: Software for tomographic image reconstruction release 2,”Physics in medicine and biology, vol. 57, pp. 867–83, Feb. 2012

    2012

  54. [54]

    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, Jan. 2024. [Online]. Available: https://www.frontiersin.org/journals/nuclear-medicine/articles/10.3389/f nume.2023.1324562/full

    work page doi:10.3389/f 2024

  55. [55]

    Integration of advanced 3D SPECT modeling into the open-source STIR framework,

    B. M. Fuster, C. Falcon, C. Tsoumpas, L. Livieratos, P. Aguiar, A. Cot, D. Ros, and K. Thielemans, “Integration of advanced 3D SPECT modeling into the open-source STIR framework,”Medical Physics, vol. 40, no. 9, p. 092502, Sep. 2013. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES, VOL. XX, NO. X, FEBRUARY 2026 15

    2013

  56. [56]

    Core Imaging Library - Part I: a versatile Python framework for tomographic imaging,

    J. S. Jørgensen, E. Ametova, G. Burca, G. Fardell, E. Papoutsellis, E. Pasca, K. Thielemans, M. Turner, R. Warr, W. R. B. Lionheart, and P. J. Withers, “Core Imaging Library - Part I: a versatile Python framework for tomographic imaging,”Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 379, no. 2204,...

    work page doi:10.1098/rsta.2020.0192 2021

  57. [57]

    The New Theories of Economics,

    V . Pareto, “The New Theories of Economics,”Journal of Political Economy, Sep. 1897. [Online]. Available: https://www.journals.uchicag o.edu/doi/10.1086/250454

    work page doi:10.1086/250454