REVIEW 3 major objections 5 minor 25 references
Quantum-Assisted Tomographic Image Refinement with Limited Qubits for High-Resolution Imaging
T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read Region-by-region quantum refinement reconstructs a 100×100 CT image from 2,500 qubits.
desk verdict The region-wise QUBO refinement genuinely cuts qubit count for 100x100 reconstructions, but the advertised 90% dose reduction does not follow from the method—the refinement stage uses the full-resolution sinogram. 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 load-bearing mechanism is the region-wise QUBO built from a zero-masked sinogram difference. For a selected region $S_l$, the pixels are encoded with $m$ binary variables each, the region's contribution is isolated as $D(\theta,s)=P(\theta,s)-P_z(\theta,s)$, and the objective $\sum_{\theta,s}(I_P(\theta,s)-D(\theta,s))^2$ is expanded into $\mathbf{x}^\top\mathbf{Q}\mathbf{x}+\mathbf{c}^\top\mathbf{x}$. Because only $n^2 m$ variables enter each solve, qubit count depends on region size rather than final image size; repeated passes over all regions, with interpolation and filtering between passes, drive convergence.
What would settle it
Measure the radiation exposure of the actual acquisition that produces the full-resolution sinogram $P$ used in Equations (6)–(9), then check whether the downscaling step changes delivered dose; if the dose is fixed by acquiring $P$, the claimed 90% reduction is absent. A direct alternative is to acquire a true 50×50 sinogram on a phantom and test whether the refinement pipeline still matches the 100×100 reconstruction.
Extended reading notes
Core claim
The central discovery is that the qubit bottleneck in quantum tomography can be bypassed by treating reconstruction as coarse-to-fine refinement over spatial regions rather than as one large binary optimization. Each refinement step builds a local QUBO whose objective is the squared difference between the projection of a qubit-encoded region and the target contribution $D(\theta,s)=P(\theta,s)-P_z(\theta,s)$ obtained by subtracting a zero-masked sinogram from the original sinogram. With 100 projection angles, a 100×100 binary phantom is reconstructed from a 50×50 start after two passes over four regions, and a sparse-view 50-angle integer phantom is recovered after smoother interpolation and Gaussian filtering. These results are offered as evidence that 2,500 and 7,500 qubits respectively suffice for 100×100 binary and integer-valued images, and that the pipeline extends naturally to sparse-view acquisition.
Load-bearing premise
The dose-reduction claim rests on the assumption that reducing the sinogram that enters the reconstruction also reduces what is physically acquired; the paper's experiments downscale an already-acquired full-resolution sinogram, so the dose saving is assumed rather than measured.
Editorial extensions
If this is right
- A 100×100 binary tomographic image can be reconstructed with 2,500 qubits, and a 100×100 integer-valued image with 7,500 qubits, under the four-region scheme demonstrated.
- Sparse-view sinograms with half the projection angles pass through the same region-wise refinement pipeline without algorithmic changes.
- If the low-resolution sinogram is acquired directly, the paper's $(1-n/N)$ formula implies starting a 500×500 reconstruction at 50×50 would require roughly 90% less projection data and radiation dose.
- Overlapping or adaptively chosen regions are supported, so regions of clinical interest could be refined selectively rather than the whole image.
- Initialization smoothness matters: nearest-neighbor upscaling stalls refinement, while smoother interpolation plus Gaussian filtering ($\sigma=1$) lets the second pass converge.
Reading between the lines
- The dose-reduction claim is not established by the reported experiments: those experiments acquire a full-resolution sinogram and then downscale it, so the radiation exposure is set by the full acquisition, not by the downscaled reconstruction input.
- Because the region-wise loop is a block-coordinate descent over spatially disjoint blocks, a classical QUBO solver could run the same decomposition; comparing classical and quantum solvers on the same region subproblems would isolate what quantum annealing actually adds.
- A natural testable extension is to acquire a true low-resolution sinogram at the detector level and check whether the same refinement pipeline recovers the 100×100 target; if it does, the dose-reduction claim would move from ratio arithmetic to measurement.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a hybrid quantum-classical tomographic reconstruction framework intended to reduce both qubit counts and radiation exposure. The method first downscales a sinogram, solves a low-resolution QUBO using the D-Wave hybrid solver, upscales the resulting image, and then iteratively refines spatial subregions by solving local QUBO problems derived from the original full-resolution sinogram. Two variants are presented: full-view sinogram downscaling and sparse-view acquisition with detector-axis downscaling. Experiments on 100x100 binary and integer-valued Shepp-Logan phantoms show visually accurate reconstructions using 2,500 and 7,500 qubits per region, respectively. The abstract and discussion further claim that reconstructing a 500x500 image from a 50x50 initialization could reduce projection data, and hence radiation dose, by up to 90%.
Significance. If the central qubit-reduction result holds, the contribution is a useful incremental step for quantum-annealing-based tomography: replacing a full 10,000-variable QUBO with four 2,500-variable regional QUBOs is practically meaningful for current D-Wave hardware, and the authors provide detailed D-Wave energy logs and a transparent iterative procedure for two phantom cases. However, the advertised clinical benefit, up to 90% radiation dose reduction, is not supported by the algorithm as written or by the experiments. The 500x500 case is never run, and the refinement stage fits each region to the original full-resolution sinogram rather than to a physically reduced-dose acquisition. The paper therefore contributes a qubit-reduction technique and a data-reduction heuristic, but not a demonstrated low-dose imaging protocol.
major comments (3)
- [Section 2.2.1, Eqs. (6)-(9); Abstract; Section IV] The claimed radiation-dose reduction is unsupported because the refinement objective in Eqs. (6)-(9) minimizes the difference between the superposed region sinogram and the original full-resolution sinogram P. Patient dose is determined by how P is physically acquired, so downscaling P after acquisition does not reduce exposure. The sparse-view experiment reduces the number of projection angles from 100 to 50, which is a 50% reduction in views, and detector-axis pooling does not by itself lower dose; no experiment approaches the claimed 90% reduction. Please either revise the abstract and discussion to describe data-volume reduction in post-processing, or redesign the pipeline so refinement is performed against a sinogram physically acquired at reduced resolution and demonstrate the dose saving with an experiment.
- [Section I and Section IV, 90% formula] The 90% reduction figure is arithmetically inconsistent with the paper's own downscaling model. Section 2.2.1 downscales both the angular and detector dimensions, with N_S = n_S d1 and M_S = m_S d2, so the number of sinogram elements scales as (n/N)^2, not as n/N. For N=500 and n=50, two-dimensional downscaling gives a 99% reduction in sinogram elements, while the claimed 90% corresponds only to one-dimensional detector downscaling. The paper should define precisely which acquisition parameters are assumed to scale with dose and reconcile the formula (1 - n/N) with the two-dimensional downscaling used in the method.
- [Section III, Table I and Figs. 3-4] The accuracy claims rely on visual inspection and on the gap between the target minimum energy and the D-Wave minimum energy. The energy gap is a self-consistency check of the QUBO solve, not a measure of tomographic fidelity, because the target minimum is computed from the same objective being solved. No quantitative image error (e.g., MAE, MSE, SSIM) is reported for the final reconstructions. Please add such metrics, and state explicitly whether the binary reconstruction is exactly equal to the phantom or merely visually identical.
minor comments (5)
- [Section 2.2.1, Eq. (10)] The phrase "moving the constant term to the left-hand side" is confusing; the constant term is subtracted from the expanded expression to obtain the QUBO form. Please rephrase for clarity.
- [Table I] Target and D-Wave energies are reported with up to ten decimal places, which likely exceeds the precision justified by the sinogram data and the solver. Please report energies with an appropriate number of significant digits.
- [Figure 2(b)] The caption says "1 × d vertical patch," but the text describes pooling along the detector/position axis. Please clarify the orientation and terminology so the figure matches the text.
- [Section IV, simple decomposition] The estimates for maximum reconstructible image size under a 10,000-qubit constraint and the O(10^6) annealing computations are stated without derivation. Please either provide a short derivation or label these as rough heuristics with clear assumptions.
- [End matter] The section titled "APPENDIX AND THE USE OF SUPPLEMENTAL FILES" is empty and should be removed, and the header date "NOVEMBER 2020" does not match the submission year.
Circularity Check
The advertised 90% dose reduction is the chosen resolution ratio (1 − n/N) restated as a result, while the refinement objective fits the original full-resolution sinogram P, making the dose claim circular by construction.
-
self definitional
[Abstract; Section 2.2.1, Eq. (6); Section 4 (Discussion)]
"By initiating the reconstruction from a lower-resolution sinogram and refining it progressively, the total radiation exposure can theoretically be reduced by up to (1 − 𝑛/𝑁) × 100% where 𝑛 and 𝑁 denote the number of pixels along one spatial dimension of the initial and final images, respectively. ... The QUBO objective for region 𝑆𝑙 is then formulated using the difference between the original and zero-masked sinograms as follows: ∑ [{𝐼𝑃 − (𝑃 − 𝑃𝑧)}(𝜃, 𝑠)]2"
The 90% reduction is exactly 1 − n/N for the chosen n = 50 and N = 500, so the claimed 'reduction in projection data' is the definitional ratio of initial to final image dimensions, not a measured consequence of the reconstruction process. Moreover, Eq. (6) minimizes against the original full-resolution sinogram P(θ,s), so the radiation dose is fixed by the acquisition of P; downscaling P after acquisition cannot lower patient exposure. The headline dose-reduction prediction is therefore a restatement of the user-chosen input sizes, equivalent to its input by construction.
full rationale
The core reconstruction pipeline is not circular: it defines a QUBO from sinogram differences, solves it region-wise with a D-Wave hybrid solver, and tests against known phantoms. The qubit counts (2,500 for the binary case and 7,500 for the integer-valued case) follow directly from the stated pixel encodings and region sizes, and the convergence of the solver toward the reported target minima is an empirical check on the optimization, though the target minima are computed from the same QUBO formulation rather than from an independent source. The circularity is confined to the second headline benefit: the 'up to 90% reduction in radiation dose' is the arithmetic quantity (1 − n/N), asserted for a 500×500 case that is not actually experimented on, while the refinement objective in Eq. (6) uses the original full-resolution sinogram P. Self-citations to the authors' earlier QUBO tomography papers are present but are not load-bearing for this defect; those prior results are externally published and are not invoked as uniqueness theorems. Overall score 6: one central advertised prediction reduces by construction, while the qubit-efficiency contribution remains self-contained and non-circular.
Assumptions & free parameters
free parameters (4)
- sinogram downscaling ratio (d1, d2 or n/N) =
d1=d2=2, n=50, N=100
- downscaled sinogram normalization factor =
1/2
- Gaussian filter sigma =
sigma=1
- number of refinement iterations =
2 per experiment
assumptions (6)
- domain assumption Parallel-beam Radon transform with line integrals is an accurate model of the imaging process.
- ad hoc to paper The 1/2 scaling of the downscaled sinogram correctly represents the low-resolution image's line integrals.
- domain assumption The QUBO formulation and binary encoding from refs [10, 11] are correct and applicable here.
- domain assumption D-Wave hybrid solver returns a solution close enough to the global minimum for the refinement to converge.
- ad hoc to paper Region-wise refinement with fixed surrounding pixels converges to the correct image.
- domain assumption The method applies to fan-beam, cone-beam, synchrotron, and electron tomography without modification.
Cite this review
Pith. "Pith review of Quantum-Assisted Tomographic Image Refinement with Limited Qubits for High-Resolution Imaging." pith.science (2026). https://pith.science/paper/UFRZVOP5
@misc{pith2026250420654,
author = {Pith},
title = {Pith review of: Quantum-Assisted Tomographic Image Refinement with Limited Qubits for High-Resolution Imaging},
year = {2026},
howpublished = {\url{https://pith.science/paper/UFRZVOP5}},
note = {Machine review of arXiv:2504.20654}
}
read the original abstract
We propose a quantum-assisted reconstruction framework for high-resolution tomographic imaging that significantly reduces both qubit requirements and radiation exposure. Conventional quantum reconstruction methods require solving QUBO (Quadratic Unconstrained Binary Optimization) problems over full-resolution image grids, which limits scalability under current hardware constraints. Our method addresses this by combining sinogram downscaling with region-wise iterative refinement, allowing reconstruction to begin from a reduced-resolution sinogram and image, then progressively upscaled and optimized region by region. Experimental validation on binary and integer-valued Shepp-Logan phantoms demonstrates accurate reconstructions under both dense and sparsely sampled projection conditions using significantly fewer qubits. We observed that nearest-neighbor interpolation may cause edge artifacts that hinder convergence, which can be mitigated by smoother interpolation and Gaussian filtering. Notably, reconstructing a 500 by 500 image from a 50 by 50 initialization demonstrates the potential for up to 90% reduction in projection data, corresponding to a similar reduction in radiation dose. These findings highlight the practicality and scalability of the proposed method for quantum-enhanced tomographic reconstruction, offering a promising direction for low-dose, high-fidelity imaging with current-generation quantum devices.
Reference graph
Works this paper leans on
-
[1]
Mitarai, K., Negoro, M., Kitagawa, M., & Fujii, K. (2018). Quantum circuit learning. Physical Review A, 98(3), 032309
2018
-
[2]
Wurtz, J., & Love, P. J. (2021). Classically optimal variational quantum algorithms. IEEE Transactions on Quantum Engineering, 2, 1-7
work page 2021
-
[3]
C., Barends, R.,
Arute, F., Arya, K., Babbush, R., Bacon, D., Bardin, J. C., Barends, R., ... & Martinis, J. M. (2019). Quantum supremacy using a programmable superconducting processor. Nature, 574(7779), 505-510
2019
-
[4]
Pelofske, E., Hahn, G., & Djidjev, H. N. (2022). P arallel quantum annealing. Scientific Reports, 12(1), 4499
work page 2022
-
[5]
Lee, H., & Jun, K. (2025). Range dependent Hamiltonian algorithms for numerical QUBO formulation. Scientific Reports, 15(1), 8819
work page 2025
-
[6]
Lee, H., & Jun, K. (2024). QUBO Refinement: Achieving Superior Precision through Iterative Quantum Formulation with Limited Qubits. arXiv preprint arXiv:2411.16138
arXiv 2024
-
[7]
O’Malley, D. & Vesselinov, V. V. ToQ.jl: A high -level programming language for D -Wave machines based on Julia, in 2016 IEEE High Performance Extreme Computi ng Conference (HPEC) (2016) pp. 1 –7 https://doi.org/10.1109/HPEC.2016.7761616
arXiv 2016
-
[8]
& Bernoudy, W
McGeoch, C., Farre, P. & Bernoudy, W. D-Wave hybrid solver service+ advantage: Technology update (Tech Rep, 2020)
2020
Show all 25 references
-
[9]
D-wave ocean software documentation,
“ D-wave ocean software documentation,” https://docs.ocean.dwavesys.com, accessed: 2024-12-03
2024
-
[10]
Jun, K. (2023). A highly accurate quantum optimization algorithm for CT image reconstruction based on sinogram patterns. Scientific Reports, 13(1), 14407
2023
-
[11]
Jun, K. (2023). Quantum optimization algorithms for CT image segmentation from X-ray data. arXiv preprint arXiv:2306.05522
2023 arXiv
-
[12]
& Quirynen, M
Jacops, R., Mraiwa, N., Steenberghe, D., Gijbels, F. & Quirynen, M. Appearance, locationi course and mor- fology of the mandibular incisive canal on spiral CT scan. Dentomaxillofac. Radiol. 31, 322–327 (2002)
2002
-
[13]
Hoff, J. A. et al. Age and gender distributions of coronary artery calcium detected by electron beam tomography in 35,246 adults. Am. J. Cardiol. 87, 1335–1339 (2001)
2001
-
[14]
Böhm, T. et al. Quantitative synchrotron X -ray tomography of the material-tissue interface in rat cortex implanted with neural probes. Sci. Rep. 9, 7646 (2019)
2019
-
[15]
Wang, X. et al. High performance model based image reconstruction. ACM SIGPLAN Not. 51(8), 1–12 (2016)
2016
-
[16]
& Timmer, J
Schomberg, H. & Timmer, J. The gridding method for image reconstruction by Fourier transformation. IEEE Trans. Med. Imaging 14(3), 596–607 (1995)
1995
-
[17]
Singh, R. et al. Artificial intelligence in image reconstruction: The change is here. Phys. Med. 79, 113–125 (2020)
2020
-
[18]
Alignment Solution for CT Image Reconstr uction using Fixed Point and Virtual Rotation Axis
Jun, K. and Yoon, S. “Alignment Solution for CT Image Reconstr uction using Fixed Point and Virtual Rotation Axis”, Sci. Rep. 7, 41218 (2017)
2017
-
[19]
Alignment theory of parallel -beam computed tomography image reconstruction for elastic -type objects using virtual focusing method
Jun, K. and Kim, D. “Alignment theory of parallel -beam computed tomography image reconstruction for elastic -type objects using virtual focusing method”, Plos One, 13(6), e0198259 (2018)
2018
-
[20]
Virtual multi-alignment theory of parallel-beam CT image reconstruction for elastic objects
Jun, K. and Jung, J. “Virtual multi-alignment theory of parallel-beam CT image reconstruction for elastic objects”, Sci. Rep. 9 (2019) 6847
2019
-
[21]
Virtual multi -alignment theory of parallel -beam CT image reconstruction for rigid objects
Jun, K. “Virtual multi -alignment theory of parallel -beam CT image reconstruction for rigid objects”, Sci. Rep. 9 (2019) 13518
2019
-
[22]
Jun, K. (2024). QUBO formulations for a system of linear equations. Results in Control and Optimization, 14, 100380. 10 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. xx, NO. x, 2020
2024
-
[23]
& Fuchs, T
Dremel, K., Prjamkov, D., Firsching, M., Weule, M., Lang, T., Papadaki, A., ... & Fuchs, T. O. (2025). Utilizing Quantum Annealing in Computed Tomography Image Reconstruction. IEEE Transactions on Quantum Engineering
2025
-
[24]
Lee, H., & Jun, K. (2025). Quantum Supremacy in Tomographic Imaging: Advances in Quantum Tomography Algorithms. arXiv preprint arXiv:2502.04830
2025 arXiv
-
[25]
Jun, K ., & Lee, H. (2023). HUBO and QUBO models for prime factorization. Scientific Reports, 13(1), 10080
2023
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.