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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [§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.
- [§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.
- [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.
- [§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)
- [§IV-C, Fig. 10 caption] The caption refers to 'GSDM' while the text and tables use 'GMSD'; please harmonize the naming.
- [§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.
- [§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.
- [§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
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
free parameters (4)
- number of subsets N =
2 (3 in comparison)
- regularization weights r1, r2
- subset regularization factors lambda_n and lambda
- diffusion hyperparameters (T, sigma_max, architecture, training length)
assumptions (4)
- domain assumption CT imaging is linear: x = AI (Eq. 1).
- domain assumption Score-based SDE framework with VE-SDE learns the true score function.
- ad hoc to paper Subset division reduces learning difficulty and improves fine-detail reconstruction.
- ad hoc to paper Alternating optimization with data consistency reaches a good solution.
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
Reference graph
Works this paper leans on
-
[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
work page 2001
-
[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
work page 2006
-
[3]
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
work page 2011
-
[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
work page 1977
-
[5]
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
work page 2006
-
[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
work page 1977
-
[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
work page 2017
-
[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
work page 2018
Show all 32 references
-
[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
2017
-
[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
2022
-
[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
2015
-
[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
2022 arXiv
-
[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
2023
-
[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
2023
-
[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
2024
-
[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
2024
-
[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
2024 arXiv
-
[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
2021
-
[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
2021
-
[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
2021 arXiv
-
[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
2022 arXiv
-
[22]
Stochastic Differential Equations. Hoboken
L. Arnold, “Stochastic Differential Equations. Hoboken .” NJ, USA: Wiley, 1974
1974
-
[23]
Stochastic Differential Equations: An introduction with applications
B. Oksendal, “Stochastic Differential Equations: An introduction with applications.” Berlin, Germany: Springer Science & Business Media , 2013
2013
-
[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
1982
-
[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
1994
-
[26]
Accessed: Apr
Low Dose CT Grand Challenge. Accessed: Apr. 6, 2017. [Online]. Avail- able: http://www.aapm.org/GrandChallenge/LowDoseCT/
2017
-
[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
1985
-
[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
1998
-
[29]
Operator discretization library (ODL),
J. Adler, H. Kohr, and O. Oktem, “Operator discretization library (ODL),”
-
[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
2019
-
[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
2001
-
[2017]
Available: https://github.com/odlgroup/odl
[Online]. Available: https://github.com/odlgroup/odl
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.