Pith. sign in

REVIEW 4 major objections 4 minor 32 references

Ordered-subsets Multi-diffusion Model for Sparse-view CT Reconstruction

T0 review · 4 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read The paper claims that splitting CT projection data into ordered subsets, training a diffusion model on each subset, and refining with a whole-sinogram diffusion model reconstructs sparse-view CT images with higher PSNR/SSIM than…

desk verdict OSMM is a plausible new combination of per-subset diffusion models plus a whole-sinogram diffusion model for sparse-view CT; the empirical claim is likely true, but the paper under-reports uncertainty and contains a concrete table error. read the letter →

arxiv 2505.09985 v1 pith:EHHFUHIP submitted 2025-05-15 eess.IV cs.CV

classification eess.IVcs.CV
keywords sparse-viewCTscore-baseddiffusionmodelorderedsubsetssinogramreconstructionunsupervisedlearningimagegenerativepriordataconsistency
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper proposes OSMM, a two-stage diffusion-model pipeline for sparse-view CT reconstruction. Instead of feeding the full sinogram to one score-based diffusion model, it splits the projection views into ordered subsets, trains a separate diffusion model on each subset (MSDM), and then uses a whole-sinogram diffusion model (OWDM) as a global consistency constraint. The authors claim this reduces the difficulty of learning the projection data distribution, preserves fine details, and corrects errors, leading to better reconstruction quality. Experiments on AAPM, CIRS phantom, and preclinical mouse data report that OSMM beats FBP, U-Net, FBPConvNet, patch-based DDPM, and GMSD in PSNR/SSIM across 50–180 views. If correct, the method offers an unsupervised, sparsity-adaptive alternative to supervised and single-model generative approaches.

What carries the argument

The operative mechanism is the ordered-subsets division of the sinogram: projection views $1,\dots,M$ are partitioned into $N$ disjoint subsets in round-robin order, with view $n$ assigned to subset $n$. Each subset is then handled by an independent variance-exploding SDE score model. The MSDM stage applies the $N$ subset predictors with per-step data consistency; the OWDM stage merges the refined subsets and runs one whole-sinogram predictor, again with data consistency. This alternating scheme carries the argument by decomposing the learning task.

What would settle it

A controlled experiment that matches total compute: train a single whole-sinogram diffusion model with the same architecture, parameter count, and total training steps or GPU time as the N subset models together. If the single model matches or exceeds OSMM's PSNR/SSIM on the same sparse-view test data, the claimed benefit of subset division is not real. Also, training OSMM with N=1 inside the same code should reproduce OWDM-only performance; any gap must then be attributed to subset learning.

Watch

Extended reading notes

Core claim

OSMM's core discovery is that ordered-subsets partitioning of the projection data turns one hard diffusion-learning problem into several easier ones: each subset diffusion model (MSDM) reconstructs its own sparse sinogram, recovering fine details, while the one-whole diffusion model (OWDM) keeps the merged sinogram globally consistent. Iterating the two stages with data-consistency updates produces sinograms that, after filtered back projection, have higher PSNR and SSIM than those from FBP, U-Net, FBPConvNet, patch-based DDPM, and GMSD. The authors further claim that the unsupervised training on full-view sinograms transfers across datasets and sparsity levels.

Load-bearing premise

The load-bearing premise is that splitting the projection data into ordered subsets genuinely lowers each diffusion model's learning difficulty while preserving the information needed for fine details; the paper gives this as a heuristic with no statistical or convergence analysis.

Editorial extensions

If this is right

  • Sparse-view CT reconstructions from as few as 50–60 views can reach PSNR values above 37 dB on AAPM data, roughly 3 dB higher than the next best generative baseline.
  • The method transfers across scanners and anatomies: models trained only on AAPM data improve CIRS phantom and preclinical mouse reconstructions without retraining.
  • At higher view counts (120, 180, 100) OSMM still gains, implying the subset prior remains useful even when the inverse problem is less ill-posed.
  • Ablation shows both stages matter: removing either MSDM or OWDM lowers PSNR, so the two-stage alternation is a genuine part of the improvement.
  • Increasing the number of subsets from 2 to 3 does not consistently improve results, suggesting there is an optimal split count.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If subset division is the real driver, the same ordered-subsets trick could be applied to other domains with large, redundant measurement spaces, such as MRI k-space or PET sinograms, where per-partition score models might ease training.
  • The paper's 'reduced learning difficulty' rationale predicts that the gain over a single model should grow with the redundancy of the full sinogram; a testable extension is to vary the amount of angular overlap in the projection data and measure the PSNR gain.
  • Because the alternating MSDM/OWDM scheme lacks a convergence guarantee, a practical extension would be to analyze it as a fixed-point iteration or replace the Predictor with a predictor-corrector sampler to see whether the data-consistency steps actually contract.
  • The subset-count effect (N=2 vs N=3) hints at a bias–variance trade-off: more subsets lower per-model difficulty but increase the chance of inconsistent merges; finding the automatic rule for N would turn the heuristic into an algorithm.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper proposes OSMM, an unsupervised ordered-subsets multi-diffusion model for sparse-view CT reconstruction. The method divides full-view sinogram data into multiple view subsets, trains a separate score-based diffusion model (MSDM) on each subset, and then applies a whole-sinogram diffusion model (OWDM) as a global consistency constraint, with alternating data-consistency updates during iterative reconstruction. Experiments on AAPM challenge data, a CIRS phantom, and preclinical mouse data compare OSMM against FBP, U-Net, FBPConvNet, patch-based DDPM, and GMSD, reporting that OSMM achieves the best PSNR/SSIM/MSE at all tested view counts. Ablation studies examine the contributions of MSDM and OWDM and the effect of the number of subsets.

Significance. If the reported empirical results hold, OSMM would provide a generally applicable unsupervised sinogram-domain prior for sparse-view CT that improves on single-model diffusion baselines, with generalization demonstrated across a human-CT test set, a phantom, and a preclinical mouse dataset. The paper has several strengths: the evaluation is on held-out data not used in training, the method is compared against established baselines rather than only against the authors' prior work, ablations isolate the contribution of each stage, and the source code is publicly linked. These strengths make the central claim worthy of serious consideration, and the paper could be a useful contribution to the diffusion-based CT reconstruction literature. However, the significance is currently limited by the lack of statistical reproducibility evidence and by an internal inconsistency in the reporting of one of the headline quantitative results.

major comments (4)
  1. [§IV-C, Tables I–III] The central claim that OSMM 'outperforms' all baselines at every view count rests on single-run point estimates with no standard deviations, no number of repeated experiments, and no significance tests. At the highest view counts the margins over GMSD are small (e.g., Table I: 43.69 vs 41.44; Table II: 46.87 vs 42.86; Table III: 40.20 vs 39.07), so without variance information the superiority claim is not statistically supported. Please report means and standard deviations over at least three independent runs, or otherwise provide a significance analysis.
  2. [§IV-C, Table II, CIRS 180-view row] The CIRS 180-view entry reports PSNR=46.87 dB and MSE=0.21e-3. For normalized images these two quantities are mutually inconsistent: PSNR=46.87 dB corresponds to MSE≈0.021e-3, while MSE=0.21e-3 corresponds to PSNR≈36.8 dB. This internal inconsistency indicates a reporting or computation error in a central results table and must be corrected and verified before the generalization claim can be accepted.
  3. [Abstract and §V] The abstract claims that OSMM outperforms traditional diffusion models in 'noise resilience,' but no experiment with noisy projections is reported anywhere in the paper. The only mention of noisy data is a qualitative statement in §V that 'when projection data is noisy, using too many subsets can hinder convergence.' Either add explicit experiments with noisy sinograms to support the noise-resilience claim, or remove that claim from the abstract.
  4. [§III-A and Algorithm 1] The method's rationale rests on the heuristic that 'dividing the projection data into multi-subsets can more effectively represent prior information' and that per-subset diffusion reduces learning difficulty. This premise is not justified theoretically or empirically, and Algorithm 1's alternating MSDM/OWDM scheme is presented without a convergence analysis or a clear objective that the iterations are guaranteed to minimize. Since the paper's contribution is precisely this decomposition, please add either a formal convergence argument for the alternating scheme or an empirical study (e.g., reconstruction error vs. iteration number, and a comparison of per-subset versus whole-data learning curves) that supports the stated motivation.
minor comments (4)
  1. [§IV-C, Fig. 10 caption] The caption refers to 'GSDM' while the text and tables use 'GMSD'; please harmonize the naming.
  2. [§IV-A] The AAPM data description states that 'the FBP algorithm was employed to obtain artifact-free images at 720 projection angles, serving as the standard reference,' but it is unclear whether the 720-view FBP reconstruction is truly artifact-free; please clarify the reference standard and its relation to the original high-dose CT images.
  3. [§II-C] The subset definitions in the text (e.g., x1 = {1, N+1, 2×N+1, ..., (N−1)×N+1}) appear to contain off-by-one inconsistencies; for N subsets and M views, the last element in x1 should be M−N+1 rather than (N−1)×N+1 when M is a multiple of N. Please fix the notation.
  4. [§IV-C, Table I] The statement in the text that 'Analysis of 90 views reconstructions reveals a similar trend to that observed with 90 views' appears to contain a typo; the second '90 views' should likely be '60 views' or another view count.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: OSMM's improvement claim is empirical and tested on held-out data; the subset-diffusion construction is a new architecture, not a renamed fit.

full rationale

I found no circular step. The derivation chain is self-contained: Eq. (3) is a standard regularized inverse-problem formulation, Eqs. (10)-(14) are the standard VE-SDE score-matching objectives applied to each subset, and the iterative scheme in Eqs. (15)-(21) and Algorithm 1 alternates diffusion-prior refinement with data consistency. The central claim that OSMM outperforms single-model diffusion baselines is supported by quantitative comparisons on held-out AAPM test slices and on CIRS phantom and preclinical mouse data not used in training; no fitted parameter is later renamed as a prediction, and no quantity reported as a prediction is equivalent by construction to a training target. The motivating statement in Section III-A that subset division 'can more effectively represent prior information' is a heuristic rather than a proven theorem, but an unsupported heuristic is not circularity. The only self-referential element is the authors' earlier GMSD model [14] being used as a comparison baseline, but the paper's conclusion does not depend on GMSD being correct; OSMM is evaluated against it as an external competitor. The Table II PSNR/MSE inconsistency and the unsupported 'noise resilience' wording are correctness and evidence concerns, not examples of circular derivation.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

No new physical entities are introduced; the MSDM and OWDM are algorithmic components, not entities with independent falsifiable handles. The main free parameters are architectural and training hyperparameters, several of which are not reported.

free parameters (4)
  • number of subsets N = 2 (3 in comparison)
    Views are divided into N subsets (round-robin); N=2 is used for main experiments, N=3 tested in Table V; no principled selection rule given.
  • regularization weights r1, r2
    Weights balancing data fidelity and MSDM/OWDM regularization in Eq. (15); values not specified in paper.
  • subset regularization factors lambda_n and lambda
    Weights in Eqs. (16) and (20) for per-subset and whole-sinogram regularization; not specified.
  • diffusion hyperparameters (T, sigma_max, architecture, training length)
    Training details inherited from Song et al. VE-SDE but not reported for this implementation; affect the quality of learned priors.
assumptions (4)
  • domain assumption CT imaging is linear: x = AI (Eq. 1).
    Standard CT model; the system matrix A and Radon transform are assumed known and exact.
  • domain assumption Score-based SDE framework with VE-SDE learns the true score function.
    The paper relies on the standard score-matching result that s_theta approximates grad log p_t; network approximation error is not quantified.
  • ad hoc to paper Subset division reduces learning difficulty and improves fine-detail reconstruction.
    Stated in Section III-A as motivation; no formal guarantee that per-subset models are easier to train or that their combination is superior.
  • ad hoc to paper Alternating optimization with data consistency reaches a good solution.
    Algorithm 1 is a heuristic alternating scheme; convergence to an optimum of Eq. (15) is not established.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Ordered-subsets Multi-diffusion Model for Sparse-view CT Reconstruction." pith.science (2026). https://pith.science/paper/EHHFUHIP

@misc{pith2026250509985,
  author       = {Pith},
  title        = {Pith review of: Ordered-subsets Multi-diffusion Model for Sparse-view CT Reconstruction},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EHHFUHIP}},
  note         = {Machine review of arXiv:2505.09985}
}
read the original abstract

Score-based diffusion models have shown significant promise in the field of sparse-view CT reconstruction. However, the projection dataset is large and riddled with redundancy. Consequently, applying the diffusion model to unprocessed data results in lower learning effectiveness and higher learning difficulty, frequently leading to reconstructed images that lack fine details. To address these issues, we propose the ordered-subsets multi-diffusion model (OSMM) for sparse-view CT reconstruction. The OSMM innovatively divides the CT projection data into equal subsets and employs multi-subsets diffusion model (MSDM) to learn from each subset independently. This targeted learning approach reduces complexity and enhances the reconstruction of fine details. Furthermore, the integration of one-whole diffusion model (OWDM) with complete sinogram data acts as a global information constraint, which can reduce the possibility of generating erroneous or inconsistent sinogram information. Moreover, the OSMM's unsupervised learning framework provides strong robustness and generalizability, adapting seamlessly to varying sparsity levels of CT sinograms. This ensures consistent and reliable performance across different clinical scenarios. Experimental results demonstrate that OSMM outperforms traditional diffusion models in terms of image quality and noise resilience, offering a powerful and versatile solution for advanced CT imaging in sparse-view scenarios.

Figures

Figures reproduced from arXiv: 2505.09985 by the authors.

Figure 5
Figure 5. Visual representation of the forward diffusion procedures of the proposed [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 29 canonical work pages

  1. [1]

    Estimated risks of radiation-induced fatal cancer from pediatric CT,

    D. J. Brenner, C. D. Elliston, E. J. Hall, and W. E. Berdon, “Estimated risks of radiation-induced fatal cancer from pediatric CT,” Am. J. Roent- genol., vol. 176, no. 2, pp. 289-296, 2001

  2. [2]

    Accurate image reconstruction from few-views and limited -angle data in divergent -beam CT,

    E. Y . Sidky, C. M. Kao, and X. Pan, “Accurate image reconstruction from few-views and limited -angle data in divergent -beam CT,” J. X-ray Sci Technol., vol. 14, no. 2, pp. 119-139, 2006

  3. [3]

    Adaptive statistical iterative reconstruction technique for radiation dose reduction in chest CT: a pilot study,

    S. Singh, M. K. Kalra, M. D. Gilman, J. Hsieh, H. H. Pien, S. R. Digu- marthy, and J. A. Shepard, “Adaptive statistical iterative reconstruction technique for radiation dose reduction in chest CT: a pilot study,” Radi- ology, vol. 259, no. 2, pp. 565-573, 2011

  4. [4]

    A. P. Dempster, ‘‘Maximum likelihood from incomplete data via the EM algorithm,’’ J. Roy. Stat. Soc. B. Methodol., vol. 39, no. 1, pp. 1-22,1977

  5. [5]

    Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,

    E. Candes, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Trans. Inf. Theory, vol. 52, no. 2, pp. 489-509, 2006

  6. [6]

    A. P. Dempster, ‘‘Maximum likelihood from incomplete data via the EM algorithm,’’ J. Roy. Stat. Soc. B, Methodol., vol. 39, no. 1, pp. 1-22, 1977

  7. [7]

    Deep con- volutional neural network for inverse problems in imaging,

    J. Kyong Hwan, M. T. McCann, E. Froustey, and M. Unser, “Deep con- volutional neural network for inverse problems in imaging,” IEEE Trans. Image Process., vol. 26, no. 9, pp. 4509-4522, 2017

  8. [8]

    A sparse-view CT reconstruction method based on combination of DenseNet and deconvo- lution,

    Z. Zhang, X. Liang, X. Dong, Y . Xie, and G. Cao, “A sparse-view CT reconstruction method based on combination of DenseNet and deconvo- lution,” IEEE Trans. Med. Imaging, vol. 37, no. 6, pp. 1407-1417, 2018

Show all 32 references
  1. [9]

    Low-dose CT with a residual encoder-decoder convolutional neural net- work

    H. Chen, Y. Zhang, M.K. Kalra, F. Lin, Y. Chen, P Liao, J Zhou, G Wang. “Low-dose CT with a residual encoder-decoder convolutional neural net- work.” IEEE Trans. Med. Imaging, vol.36, no.12, pp. 2524-2535, 2017

  2. [10]

    Multi-domain integrative Swin transformer network for sparse-view tomographic reconstruction,

    J. Pan, H. Zhang, W. Wu, Z. Gao, W. Wu. "Multi-domain integrative Swin transformer network for sparse-view tomographic reconstruction," Patterns, vol. 3, no.6, pp. 100498, 2022

  3. [11]

    Deep unsupervised learning using nonequilibrium thermodynamics,

    J. Sohl -Dickstein, E. Weiss, N. Maheswaranathan , and S. Ganguli, “Deep unsupervised learning using nonequilibrium thermodynamics,” ICML., pp 2256–2265, 2015

  4. [12]

    Patch-based denoising diffusion prob- abilistic model for sparse-view CT reconstruction

    W. Xia, W. Cong, and G. Wang. "Patch-based denoising diffusion prob- abilistic model for sparse-view CT reconstruction." arXiv preprint arXiv:2211.10388, 2022

  5. [13]

    Sub -volume-based denoising diffusion probabilistic model for cone-beam CT reconstruction from in- complete data,

    W. Xia, C. Niu, W. Cong, and G. Wang, “Sub -volume-based denoising diffusion probabilistic model for cone-beam CT reconstruction from in- complete data,” CoRR, 2023

  6. [14]

    Generative modeling in sinogram domain for sparse-view CT reconstruction,

    B. Guan, C. Yang, L. Zhang, S. Niu, M. Zhang, and Q. Liu. “Generative modeling in sinogram domain for sparse-view CT reconstruction,” IEEE Trans. Radiat. Plasma Med. Sci., vol. 8, no.2, 195 – 207, 2023

  7. [15]

    Stage-by-stage wavelet opti- mization refinement diffusion model for sparse-view CT reconstruction

    K. Xu, S. Lu, B. Huang, W. Wu, Q. Liu, “Stage-by-stage wavelet opti- mization refinement diffusion model for sparse-view CT reconstruction.” IEEE Trans. Med. Imaging, doi:10.1109, 2024

  8. [16]

    Sparse-view spectral CT reconstruction and material decomposition based on multi-channel SGM

    Y. Liu, X. Zhou, C. Wei, Q. Xu, “Sparse-view spectral CT reconstruction and material decomposition based on multi-channel SGM.” IEEE Trans. Med. Imaging, doi:10.1109, 2024

  9. [17]

    DPER: Diffusion prior driven neural representation for lim- ited angle and sparse view CT reconstruction

    C. Du, X. Lin, Q. Wu, X. Tian, Y. Su, Z. Luo, H. Wei, S.K. Zhou, J Y u, Y Zhang, “DPER: Diffusion prior driven neural representation for lim- ited angle and sparse view CT reconstruction.” arXiv preprint arXiv:2404.17890, 2024

  10. [18]

    Score-based generative modeling through stochastic differential equations,

    Y . Song, J. Sohl-Dickstein, D. P. Kingma, A. Kumar, S. Ermon, and B. Poole, “Score-based generative modeling through stochastic differential equations,” ICLR, 2021

  11. [19]

    Deep generative learn- ing via Schrödinger bridge,

    G. Wang, Y . Jiao, Q. Xu, Y . Wang, and C. Yang, “Deep generative learn- ing via Schrödinger bridge,” ICML, pp. 10794–10804, 2021

  12. [20]

    Heavy -tailed den oising score matching,

    J. Deasy, N. Simidjievski, and P. Liò, “Heavy -tailed den oising score matching,” arXiv preprint arXiv:2112.09788, 2021

  13. [21]

    Maxi- mum likelihood training of implicit nonlinear diffusion models,

    D. Kim, B. Na, S. J. Kwon, D. Lee, W. Kang, and I. -C. Moon, “Maxi- mum likelihood training of implicit nonlinear diffusion models,” arXiv preprint arXiv:2205.13699, 2022

  14. [22]

    Stochastic Differential Equations. Hoboken

    L. Arnold, “Stochastic Differential Equations. Hoboken .” NJ, USA: Wiley, 1974

  15. [23]

    Stochastic Differential Equations: An introduction with applications

    B. Oksendal, “Stochastic Differential Equations: An introduction with applications.” Berlin, Germany: Springer Science & Business Media , 2013

  16. [24]

    Reverse-time diffusion equation models,

    B. D. Anderson, “Reverse-time diffusion equation models,” Stochastic Processes Appl., vol. 12, no. 3, pp. 313-326, 1982

  17. [25]

    Accelerated image reconstruction using or- dered subsets of projection data

    H.M. Hudson, R.S. Larkin, “Accelerated image reconstruction using or- dered subsets of projection data.” IEEE Trans. Med. Imag, vol.13, no.4, pp. 601-609, 1994

  18. [26]

    Accessed: Apr

    Low Dose CT Grand Challenge. Accessed: Apr. 6, 2017. [Online]. Avail- able: http://www.aapm.org/GrandChallenge/LowDoseCT/

  19. [27]

    Fast calculation of the exact radiological path fora three- dimensional CT array,

    R. L. Siddon, “Fast calculation of the exact radiological path fora three- dimensional CT array,” Med. Phys., vol. 12, no. 2, pp. 252-255, 1985

  20. [28]

    A fast algorithm to calculate the exact radiological path through a pixel or voxel space,

    F. Jacobs, E. Sundermann, B. de Sutter, M. Christiaens, and I. Lemahieu, “A fast algorithm to calculate the exact radiological path through a pixel or voxel space,” J. Comput. Inf. Technol., vol. 6, no. 1, pp. 89-94, 1998

  21. [29]

    Operator discretization library (ODL),

    J. Adler, H. Kohr, and O. Oktem, “Operator discretization library (ODL),”

  22. [30]

    Deep-neural-network-based sinogram synthesis for sparse -view CT image reconstruction,

    H. Lee, J. Lee, H. Kim, B. Cho, and S. Cho, “Deep-neural-network-based sinogram synthesis for sparse -view CT image reconstruction, ” IEEE Trans. Radiat. Plasma Med. Sci., vol. 3, no. 2, pp. 109-119, 2019

  23. [31]

    Statistically regulated and adaptive EM reconstruction for emission computed tomography

    D. J. Kadrmas, “Statistically regulated and adaptive EM reconstruction for emission computed tomography.” IEEE Trans. Nucl. Sci. , vol.48, no.3, pp. 790-798, 2001

  24. [2017]

    Available: https://github.com/odlgroup/odl

    [Online]. Available: https://github.com/odlgroup/odl

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.