REVIEW 3 major objections 5 minor 10 references
Signal-to-noise and spatial resolution in in-line imaging. 2. Phase-contrast tomography
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Phase-contrast CT lifts 3D image SNR by up to about 11 times in breast tissue.
desk verdict A substantive PB-CT paper with a new dose-normalized quality metric and honest experimental data, but the strong gain claims rest on an unvalidated monomorphous assumption and one partially circular detector calibration. 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 engine of the argument is the homogeneous transport-of-intensity equation (TIE-Hom), also called Paganin's method: for a sample with a constant ratio $\gamma = \delta/\beta$, free-space propagation maps the object-plane intensity to $I_R = (1-a^2\nabla_\perp^2)I_0$, and inversion is a convolution with a 2D or 3D filter (a modified Bessel $K_0$ function in 2D, a Yukawa potential in 3D). Because the TIE-Hom operator commutes with the X-ray projection operator, the 3D reconstruction can be written as $\beta = \mathcal{P}^{-1}(1-a^2\nabla^2)^{-1} C_R$: the noise-suppressing Paganin filter and the noise-amplifying CT ramp filter act on the same frequencies, with the net effect that $\mathrm{SNR}^2$ in PB-CT becomes nearly resolution-independent. The paper's quantitative predictions are organized by the dimensionless ratio $\gamma/N_F$, where $\gamma \approx 869$ for breast tissue at 32 keV and $N_F$ is the minimal Fresnel number; when $\gamma/N_F$ is 10–100, each detected photon can carry many bits of phase information, which is the source of the 'beneficial violation' of the noise-resolution uncertainty relation.
What would settle it
Scan a calibration phantom with known inserts having two or more different phase-to-absorption ratios under the same 32 keV, 5 m, 4 mGy protocol. If single-ratio TIE-Hom retrieval reproduces the measured 3D SNR gains of about 4–11 but the reconstructed absorption values deviate from the known insert values by more than the stated resolution, the monomorphous assumption is the cause; a multi-material retrieval that removes the deviations would then confirm that part of the gain was an artifact.
Extended reading notes
Core claim
The paper's central claim is that in three-dimensional propagation-based phase-contrast CT, applying the homogeneous transport-of-intensity (TIE-Hom, Paganin) retrieval before filtered back-projection converts the phase signal into a genuine SNR advantage at fixed dose and resolution, while also taming the noise-amplifying ramp filter of CT. Quantitatively, the authors define the biomedical imaging quality $Q_C = Q_S \, C_m \, (R_{\mathrm{ab,air}}/R_{\mathrm{ab,tissue}})^2$ (eq. 14), with $Q_S$ the intrinsic quality equal to $\mathrm{SNR}^2/(I_{\mathrm{in}}\Delta^n)$, $C_m$ the sample contrast, and the dose terms normalizing by absorbed dose or mean glandular dose. They measure 3D SNR gain factors $G_3 \approx 4.3$ (Xineos flat-panel detector) and $G_3 \approx 10.9$ (Eiger photon-counting detector) in mastectomy samples at 32 keV, 4 mGy MGD, and 5 m propagation, values that lie between the two theoretical estimates in eqs. (22) and (23). $Q_S$ rises from 0.93 to 1.77 and 3.91, above the Epanechnikov bound (a constant just above unity) that limits absorption-only imaging. The authors also show $Q_C$ for a CT volume is typically much smaller than unity because of CT's ill-posedness, but the measured gain factor transfers directly to $Q_C$, so the improvement is a real imaging-quality gain rather than a noise-reduction artifact.
Load-bearing premise
The whole calculation treats each breast sample as one homogeneous material with a single fixed phase-to-absorption ratio (about 869 at 32 keV), and if real breast tissue violates that, Paganin retrieval can introduce artifacts and the measured gain no longer reflects true signal improvement.
Editorial extensions
If this is right
- At equal dose and spatial resolution, PB-CT with TIE-Hom retrieval can raise 3D SNR by roughly an order of magnitude over attenuation-based CT in breast tissue, corresponding to a potential ~100-fold dose reduction for matched image quality because dose scales as $\mathrm{SNR}^2$.
- The biomedical imaging quality $Q_C$ gives a single dimensionless number for comparing imaging chains on a given sample class; for PB-CT it is maximized when the mean transmission through the sample is $e^{-2}\approx 0.135$ ($\mu L=2$).
- For reconstructed CT volumes, $Q_C$ is inherently small (about $10^{-3}$ without retrieval in this experiment) because filtered back-projection amplifies high-frequency noise; projected or slab-averaged reconstructions have much higher $Q_C$, so the dose-fractionation theorem does not hold for CT.
- Detector choice matters through spatial resolution and PSF shape: the photon-counting Eiger detector with 75 $\mu$m pixels gave roughly 2.5 times the 3D gain of the flat-panel Xineos detector, consistent with Fresnel-number scaling.
- The framework provides quantitative targets for choosing propagation distance, detector resolution, and the TIE-Hom parameter $a$ when moving PB-CT from synchrotron beamlines toward clinical X-ray sources.
Reading between the lines
- The single-$\gamma$ monomorphous assumption is the load-bearing approximation: if real breast tissue mixes compositions within a voxel, part of the measured SNR gain may come from the retrieval smoothing over mismatched $\delta/\beta$ ratios, and a multi-material retrieval would separate the true phase signal from that artifact.
- The same $Q_C$ metric could be extended to polychromatic clinical X-ray sources, where beam hardening changes the effective $\delta/\beta$ ratio; a spectral or calibration correction would be needed, and $Q_C$ would then allow a fair comparison between synchrotron and clinical PB-CT.
- Because the gain scales with $\gamma/N_F$, propagation distance and detector pixel size are not independent knobs: an optimal Fresnel number should match the Paganin filter width to the detector PSF, a prediction testable by scanning a fixed phantom at several distances.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This second part of the authors' study on in-line phase-contrast imaging applies the noise-resolution uncertainty (NRU) framework to propagation-based phase-contrast computed tomography (PB-CT). It introduces a dimensionless 'biomedical imaging quality characteristic' Q_C that combines contrast, SNR, spatial resolution, and radiation dose, and derives analytical expressions for Q_C in CT and PB-CT (eqs. 14-24). The theory is tested on synchrotron PB-CT data from two human mastectomy samples scanned at 32 keV with 4 mGy mean glandular dose using a flat-panel (Xineos) and a photon-counting (Eiger) detector. The measured 3D SNR gain factors are 4.3 (Xineos) and 10.9 (Eiger), and these lie between the theoretical estimates from eqs. (22) and (23). The paper concludes that PB-CT can deliver large, dose-relevant SNR gains over conventional CT while preserving spatial resolution.
Significance. The paper makes a useful contribution by proposing a sample- and dose-aware quality metric for biomedical phase-contrast imaging and by providing an explicit theoretical derivation of its value in CT, linking the smallness of Q_C to the ill-posedness of CT reconstruction (Appendix C). The experimental data are valuable: they are from full intact mastectomy samples at a clinically relevant dose, and the Eiger comparison is a genuinely non-trivial test, with the measured gain of 10.9 falling between the independent theoretical estimates of 9.5 and 11.6. The theoretical predictions for Q_C (9.5e-3) and the measured value (7.2e-3) agree to within about 30%, which is reasonable given the modelling approximations. The main weaknesses concern the support for the claim that the SNR gains are achieved at unchanged spatial resolution: the Xineos resolution is partly estimated from the same SNR data being compared, and no resolution check is performed after phase retrieval in the reconstructed volumes.
major comments (3)
- [Section 6, after Table 1] The Xineos spatial resolution is re-estimated as 212 μm from the measured SNR ratio via the NRU, and this value is then used to compute the Fresnel number and the theoretical gain factors for Xineos. This makes the Xineos gain comparison partially self-fulfilling: the measured SNR ratio is consistent with the NRU by construction, so the subsequent agreement between the measured 3D gain (4.3) and the theoretical range (3.2-4.9) is not an independent validation for that detector. An independent measurement of the Xineos resolution in the object plane (e.g. from an edge or line-pair phantom) is needed to break the circularity.
- [Section 5, eq. (20) and following text; Section 6] The reduction of the gain factor G_3 to the ratio of measured SNRs relies on the assertion that both the signal and the spatial resolution are identical before and after TIE-Hom retrieval. That assertion requires the sample to be strictly monomorphous, with a single ratio gamma = delta/beta valid at every point. The mastectomy samples are mixtures of adipose and glandular tissue, and the value gamma = 869.4 used here is the ratio of the differences in delta and beta between the two tissue types, not a locally valid material constant. For a non-monomorphous object the inverse TIE-Hom operator does not recover the true high-frequency signal; the measured SNR gain may therefore partly be a resolution-for-noise trade-off. No spatial resolution or artifact quantification is reported in the reconstructed volumes after retrieval, so the central claim that the gains of 4.3 and 10.9 are achieved 'without loss of spatial resolution' is currently not supported for these samples. The authors should measure the resolution in the reconstructed volumes (e.g. with a resolution phantom or via Fourier ring correlation) or otherwise bound the resolution loss.
- [Section 6, Tables 1 and 2] The SNR values and gain factors are reported without error bars or a noise model. The intra-scan spread in Table 2 is small (e.g. Eiger gains 10.7-11.0, Xineos gains 4.0-4.5), but the central quantitative claims ('approximately 4.3' and '10.9') and the comparison with theory (including the 30% discrepancy between theoretical and measured Q_C for Eiger, 9.5e-3 vs 7.2e-3) cannot be statistically evaluated without uncertainties. A per-pixel noise model or at least the standard deviation of the measured SNR values over independent slices should be reported.
minor comments (5)
- [Figure 1 caption] The caption states 'MGD 4 μGy', while the text and Figure 2 caption state 4 mGy; this is presumably a unit typo and should be corrected.
- [Section 6] The displayed computation of the incident fluence I_in is typeset in a way that is hard to follow ('8 10^3 Gy / (600 7 10^17 Gym)'); please rewrite the arithmetic with explicit units and intermediate steps so that the result 0.19 μm^-2 can be verified.
- [Section 5] The statement that the gain factor is independent of the parameter a in TIE-Hom retrieval is qualified by the conditions for validity of eq. (8), but this qualification is easy to miss; please state explicitly that the invariance holds only when the retrieval parameter is within the regime where the noise-resolution duality applies.
- [Appendix C vs Section 5] Equation (24) in the main text and equation (C3) in Appendix C present essentially the same result; the duplication is not a problem, but the authors should cross-reference them explicitly to avoid confusion.
- [Section 6] The value of the breast tissue gamma = 869.4 is cited from an online calculator (TS-Imaging, 2025); since this parameter is central to the quantitative comparison, the underlying elemental composition and density assumptions should be stated, or at least a range of values relevant to the 0-100% glandularity interval should be reported.
Circularity Check
Xineos theoretical gain comparison is partly calibrated: its spatial resolution is inferred from measured SNR via NRU and then fed back into the gain formulas.
-
fitted input called prediction
[Section 6, paragraph after Table 1; use of calibrated Δ[Xineos] in eqs. (22)-(23)]
"Using the NRU, eq.(8), with the measured value Δ[Eiger] ≈ 83.2 μm and the ratio of average measured SNRs for Xineos and Eiger from Table 1, the value of Δ in the object plane for the Xineos detector can be estimated as Δ[Xineos] ≈ 83.2 μm × (85.6 / 33.6) ≈ 212.0 μm. ... We will use the value of Δ[Xineos] ≈ 212.0 μm in the calculations below, which corresponds to N_{F,MN}[Xineos] ≈ 240.1 and G2[Xineos] ≈ (869.4 / 240.1)^{1/2} ≈ 1.9."
The Xineos detector resolution is not independently measured; it is inferred from the measured 2D SNR ratio using the NRU. Substituting that calibrated Δ into the gain expression G2 = (γ/N_F)^{1/2} with N_F = Δ²/(λR) yields G2[Xineos] = G2[Eiger]/(SNR_X/SNR_E) ≈ 4.8/2.55 ≈ 1.9, which is exactly the average 2D gain already measured in Table 1. The same fitted Δ is then inserted into eqs. (22) and (23) to report 'theoretical' 3D gain factors of 3.2 and 4.9 for Xineos. The Xineos theory-experiment comparison is therefore partly self-fulfilling: a parameter required for the prediction was calibrated from the very SNR data to which the theoretical gain is compared. The Eiger comparison uses an independently assumed pixel/PSF-based resolution and remains unaffected.
full rationale
The central experimental result, the measured 3D SNR gains of about 4.3 (Xineos) and 10.9 (Eiger) in Table 2, is a direct measurement and is not circular. The Eiger theoretical comparison is based on an independently assumed detector resolution (75 μm pixel / rectangular PSF) and gives a genuine prediction (measured 10.9 vs theoretical 9.5-11.6). The new biomedical quality metric QC is introduced as an explicit definition, not derived from itself. The TIE-Hom/Paganin retrieval and NRU framework are externally established and mathematically derived in prior work; citing them is not a circularity. The one genuine circular step is confined to the Xineos column of the theoretical comparison: Δ[Xineos] is estimated from the measured SNR ratio via NRU and then used to compute the 'theoretical' gains, making that specific comparison partially self-fulfilling. This affects a secondary consistency check rather than the measured gain claim itself, so the paper is only mildly circular overall. The monomorphous-sample assumption is a correctness/validity concern, not a circularity.
Assumptions & free parameters
free parameters (1)
- Effective spatial resolution of Xineos detector in object plane (Delta[Xineos]) =
212.0 micrometers
assumptions (6)
- domain assumption The imaged object is monomorphous/homogeneous, i.e., the delta-to-beta ratio is constant within the sample, so the TIE-Hom equation (4) applies with a single parameter gamma = delta/beta.
- domain assumption The noise-resolution duality (NRU), eq. (8), holds: SNR squared divided by Delta to the power n is invariant under linear photon-number-conserving filters and is bounded by unity for absorption imaging.
- domain assumption Photon detection statistics are Poissonian, and spatial ergodicity allows SNR and variance to be estimated from flat-field areas.
- domain assumption Paraxial scalar wave propagation with thin-object approximation and near-Fresnel region applies.
- standard math CT reconstruction is linear with known ramp filter; a specific FBP implementation with nearest-neighbor interpolation gives the constant 12/pi squared in eq. (C1)/(24).
- standard math Fourier analysis, Parseval's theorem, and kernel properties (Bessel and Yukawa functions) are valid.
invented entities (1)
-
Biomedical imaging quality characteristic Q_C
Cite this review
Pith. "Pith review of Signal-to-noise and spatial resolution in in-line imaging. 2. Phase-contrast tomography." pith.science (2026). https://pith.science/paper/2MKWE5RN
@misc{pith2026250620277,
author = {Pith},
title = {Pith review of: Signal-to-noise and spatial resolution in in-line imaging. 2. Phase-contrast tomography},
year = {2026},
howpublished = {\url{https://pith.science/paper/2MKWE5RN}},
note = {Machine review of arXiv:2506.20277}
}
read the original abstract
In the first part of this paper, quantitative aspects of propagation-based phase-contrast imaging (PBI) were investigated using theoretical and numerical approaches, as well as experimental two-dimensional PBI images collected with plane monochromatic X-rays at a synchrotron beamline. In this second part, signal-to-noise ratio, spatial resolution and contrast are studied in connection with the radiation dose in three-dimensional PBI images of breast tissue samples obtained using propagation-based phase-contrast computed tomography (PB-CT) with energy-integrating and photon-counting detectors. The analysis is based on the theory of PBI and PB-CT using the homogeneous Transport of Intensity equation (Paganin's method). A biomedical image quality characteristic, suitable for quantitative assessment of X-ray images of biological samples, is introduced and applied. The key factors leading to high values of the biomedical imaging quality in PBI and to relatively low values of the same quality metric in CT imaging are identified and discussed in detail. This study is aimed primarily at developing tools for quantitative assessment and optimization of medical PB-CT imaging, initially at synchrotron facilities, with the prospect of subsequent transfer of the technology to medical clinics.
Figures
Reference graph
Works this paper leans on
-
[1]
Arhatari, B.D., Stevenson, A.W., Abbey, B., Nesterets, Y.I., Maksimenko, A., Hall, C.J., Thompson, D., Mayo, S.C., Fiala, T., Quiney, H.M., Taba, S.T., Lewis, S.J., Brennan, P.C., Dimmock, M., Hausermann, D. & Gureyev, T.E. (2021). Appl. Sci. 11,
work page 2021
-
[6]
This theoretical estimation confirms that ,3 0[]CDQ in the present configuration should indeed be expected to be much smaller than unity. It was already noted in Section 4 that ,3 ~1CDQ requires 3~ 10SNR , while in the present configuration we had 0[ ] 0.52SNR . The key reason for the relatively small value of ,3CDQ in CT reconstruction is the ill-po...
work page 2009
-
[8]
A related question is about the measured values of ,3 0[]CDQ being significantly smaller than unity. It can be shown from eq.(B3), or from the results in (Gureyev et al., 2016), that 2 2 2 4 2 0[ ] (12 / ) / a pixSNR n , (C1) where 02k= . Note that the mean number of photons collected in a single detector pixel with the size pix during the ...
work page 2016
-
[40]
Paganin D., Gureyev T.E., Mayo S.C., Stevenson A.W., Nesterets Ya.I. & Wilkins, S.W. (2004). J. Micros. 214, 315-327. Pakzad, A., Turnbull, R., Mutch, S.J., Leatham, T.A., Lockie, D., Fox, J., Kumar, B., Häusermann, D., Hall, C.J., Maksimenko, A., Arhatari, B.D., Nesterets, Y.I., Entezam, A., Taba, S.T., Brennan, P.C., Gureyev, T.E. & Quiney, H.M. (2025)....
work page Pith review arXiv 2004
-
[88]
Gureyev, T.E., Paganin, D.M., Myers, G.R., Nesterets, Ya.I. & Wilkins, S.W. (2006). Appl. Phys. Lett. 89, 034102. Gureyev, T.E., Nesterets, Y.I., Stevenson, A.W., Miller, P.R., Pogany, A. & Wilkins, S.W. (2008). Opt. Expr. 16, 3223-3241. Gureyev, T.E., Nesterets, Y.I., de Hoog, F., Schmalz, G., Mayo, S.C., Mohammadi, S. & Tromba, G. (2014). Opt. Expr. 22,...
work page 2006
-
[126]
DOI: https://dx.doi.org/10.18434/T4D01F; Johns, P.C. & Yaffe, M.J. (1985). Med. Phys. 12, 289-296. Kato, H. (2014). Japan. J. Radiolog. Technol. 70, 684-691. Kitchen, M. J., Buckley, G. A., Gureyev, T. E., Wallace, M. J., Andres-Thio, N., Uesugi, K., Yagi, N. & Hooper, S. B. (2017). Sci. Rep. 7, 15953. Liu, Q., Suleiman, M. E., McEntee, M. F. & Soh, B. P....
doi:10.18434/t4d01f 1985
-
[1995]
does not hold for the types of CT imaging considered in the present paper, which is again explained by the mathematical ill-posedness of CT reconstruction; see also a related discussion in (Gureyev et al., 2018)
work page 2018
-
[2000]
that the noise variance is equal to the integral of noise power spectrum (NPS). In particular, 2 ,, ()RR Wd = qq , where , ()RW q is the NPS of () r 38 reconstructed from a PB-CT scan collected at z = R. An explicit expression for , ()RW q can be obtained in conjunction with eq.(6), in the case of a PB-CT scan with aM projections uniformly distri...
work page 2014
Show all 10 references
-
[4120]
& Gureyev, T.E
Baran, P., Pacile, S., Nesterets, Y.I., Mayo, S.C., Dullin, C., Dreossi, D., Arfelli, F., Thompson, D., Lockie, D., McCormack, M., Taba, S.T., Brun, F., Pinamonti, M., Nickson, C., Hall, C., Dimmock, M., Zanconati, F., Cholewa, M., Quiney, H., Brennan, P.C., Tromba, G. & Gurey...
2017
-
[7890]
G., Quiney, H.M
Gureyev, T.E., Brown, H. G., Quiney, H.M. & Allen, L.J. (2022). J. Opt. Soc. Am. A 39, C143-C155. 35 Gureyev, T.E., Paganin, D.M. & Quiney, H.M. (2024). J. Synchrotron Rad., 31, 896-909. Higgins, P.D., Attix, F.H., Hubbell, J.H and Seltzer, S.W., Berger, M.J. & Sibata, C.H. (1...
2022
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.