Improving depth-resolution, in-plane contrast, and reducing non-uniformity artifacts for wide-angle DBT

Adrian A. Sanchez; Dan Xia; Emil Y. Sidky; Ingrid S. Reiser; John Paul Phillips; Xiaochuan Pan; Zheng Zhang

arxiv: 2605.29133 · v1 · pith:3B4NQ7DWnew · submitted 2026-05-27 · 🧮 math.OC · physics.med-ph

Improving depth-resolution, in-plane contrast, and reducing non-uniformity artifacts for wide-angle DBT

Emil Y. Sidky , Adrian A. Sanchez , John Paul Phillips , Zheng Zhang , Dan Xia , Ingrid S. Reiser , Xiaochuan Pan This is my paper

Pith reviewed 2026-06-29 10:14 UTC · model grok-4.3

classification 🧮 math.OC physics.med-ph

keywords digital breast tomosynthesiswide-angle DBTimage reconstructionsparsity regularizationbackground estimationnon-uniformity artifactsdepth resolution

0 comments

The pith

Augmenting sparsity-regularized reconstruction with background estimation improves depth resolution and reduces non-uniformity in wide-angle DBT.

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

This paper extends a prior two-stage sparsity-regularized iterative reconstruction method for wide-angle digital breast tomosynthesis. The low-resolution first stage now jointly estimates the image and an additional background component whose role is to absorb low-frequency non-uniformity artifacts. Results on one patient data set acquired with a wide-angle system show gains in depth resolution, in-plane contrast, and uniformity. A sympathetic reader would care because clearer separation of structures at different depths and more uniform images could aid interpretation of breast tomosynthesis exams, provided the gains hold beyond the single case.

Core claim

The image reconstruction algorithm is an extension of prior work on sparsity-regularized iterative image reconstruction performed in two stages. The first stage consists of a low-resolution reconstruction that exploits sparsity for quantitative accuracy; this stage is augmented with a formulation that includes the estimation of a background image which absorbs low-frequency artifacts that cause image non-uniformity. The new algorithm is demonstrated on a patient case for which the data are acquired on a wide-angle DBT system, and the results indicate that the algorithm design goals of improved depth resolution, in-plane contrast, and reduced non-uniformity artifacts have been met.

What carries the argument

The two-stage iterative reconstruction in which the low-resolution stage jointly estimates the image and a background image under sparsity regularization so that the background term absorbs low-frequency non-uniformity artifacts.

If this is right

The two-stage method preserves quantitative accuracy from the sparsity-regularized stage while addressing non-uniformity through the background term.
The approach yields measurable gains in depth resolution and in-plane contrast on wide-angle DBT patient data.
Non-uniformity artifacts are reduced by allowing the background component to capture low-frequency variations.
Further empirical results on multiple cases and task-based assessment are required to confirm broader applicability.

Where Pith is reading between the lines

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

The background-estimation idea might transfer to other limited-angle tomography settings that suffer from similar low-frequency shading.
Phantom studies with calibrated depth and contrast targets could quantify any accuracy trade-off introduced by the background term.
Reader studies using detection or characterization tasks would test whether the reported image-quality gains translate to clinical decision performance.

Load-bearing premise

The added background-image estimation absorbs low-frequency non-uniformity artifacts without degrading the quantitative accuracy from the sparsity-regularized low-resolution stage, and the improvement seen on one patient case will hold for other cases.

What would settle it

Reconstruction of additional wide-angle DBT patient cases or a phantom with known ground-truth uniformity and depth-separated features, followed by quantitative metrics showing no gain in depth resolution or contrast or an increase in non-uniformity, would falsify the claim that the goals are met.

read the original abstract

Purpose: This work aims to develop an image reconstruction algorithm for wide-angle digital breast tomosynthesis (DBT) that has improved depth resolution and in-plane contrast while reducing non-uniformity artifacts. Approach: The image reconstruction algorithm is an extension of our prior work on sparsity-regularized iterative image reconstruction. The algorithm is performed in two stages as explained in a prior work. The first stage consists of a low-resolution reconstruction that exploits sparsity for quantitative accuracy. In this work, this first stage is augmented with a formulation that includes the estimation of a "background" image, which absorbs low-frequency artifacts that cause image non-uniformity. Results: The new algorithm is demonstrated on a patient case for which the data are acquired on a wide-angle DBT system. Conclusion: The results on the shown case indicate that the algorithm design goals have been met, but additional empirical results and task-based assessment are needed to strengthen this conclusion.

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.

Desk Editor's Note private letter to a colleague

Incremental extension of the authors' own sparsity work, adding a background term and showing it on one patient case without numbers or comparisons.

read the letter

The core addition here is a background-image estimation step folded into the first stage of their existing two-stage sparsity-regularized reconstruction. That step is meant to soak up low-frequency non-uniformity while leaving the quantitative low-resolution output intact. The paper shows the result on a single wide-angle DBT patient case and states that depth resolution, in-plane contrast, and uniformity all look better.

What the work actually supplies is a straightforward engineering tweak on top of the group's prior papers. The two-stage structure and the sparsity term are already published; the new piece is the joint estimation of the background image. The authors are explicit that this is an extension and that more cases plus task-based reader studies are still needed.

The obvious limitation is the validation. There are no ROI measurements, no contrast-to-noise numbers, no error bars, and no side-by-side comparison of the algorithm with and without the background term. Because the background is estimated jointly, any signal leakage would quietly affect the very quantitative values the first stage is supposed to protect, yet nothing in the presented material rules that out. A single qualitative case cannot carry that claim.

This paper is for people already working on DBT reconstruction algorithms who want to see the next small step in this particular line of sparsity methods. A reader outside that narrow group will not get much from it. The manuscript is coherent on its own terms and the authors do not overstate what they have shown, so it is worth sending to referees who can ask for the missing quantitative checks and additional cases.

Referee Report

2 major / 1 minor

Summary. The manuscript extends the authors' prior sparsity-regularized iterative reconstruction for wide-angle DBT by adding a background-image estimation step in the low-resolution stage to absorb low-frequency non-uniformity artifacts. The two-stage algorithm is demonstrated on data from a single patient case acquired on a wide-angle DBT system. Visual results are presented as indicating that the design goals of improved depth resolution, in-plane contrast, and reduced non-uniformity artifacts have been met, although the conclusion notes that additional empirical results and task-based assessment are required.

Significance. If the background estimation step can be shown to remove only low-frequency artifacts while preserving the quantitative accuracy delivered by the sparsity-regularized stage, the method would address a practical limitation in wide-angle DBT reconstruction. The work is a direct continuation of the group's earlier publications and supplies a concrete algorithmic modification, but the single-case qualitative demonstration provides limited external validation.

major comments (2)

[Results] Results section: the demonstration consists of qualitative images from a single patient case with no quantitative metrics (contrast, resolution, or artifact measures), no error bars, and no direct comparison of reconstructions with versus without the background term.
[Abstract and algorithm description] Abstract and § on algorithm formulation: the central claim that the background image 'absorbs low-frequency artifacts' without degrading quantitative accuracy from the sparsity stage is not tested; no ROI statistics, ground-truth reference, or ablation study is supplied to rule out mid-frequency signal leakage into the background image.

minor comments (1)

[Introduction and Methods] The manuscript repeatedly refers to 'a prior work' for the two-stage framework; add a concise self-contained summary of the baseline algorithm so that the novel background-estimation term can be evaluated independently.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed review and constructive comments on our manuscript. We address each major comment below in a point-by-point manner. The work is presented as a preliminary algorithmic extension demonstrated on a single patient case, consistent with the stated conclusion that further empirical validation is required.

read point-by-point responses

Referee: [Results] Results section: the demonstration consists of qualitative images from a single patient case with no quantitative metrics (contrast, resolution, or artifact measures), no error bars, and no direct comparison of reconstructions with versus without the background term.

Authors: We agree that the results are limited to qualitative visual assessment on a single patient case without quantitative metrics, error bars, or explicit with/without-background comparisons. This scope is stated in the manuscript conclusion, which notes the need for additional empirical results and task-based assessment. The figures provide side-by-side visual evidence of the effect of the background term on non-uniformity while preserving structures, but we acknowledge that quantitative metrics would require a larger cohort and are outside the current demonstration. revision: no
Referee: [Abstract and algorithm description] Abstract and § on algorithm formulation: the central claim that the background image 'absorbs low-frequency artifacts' without degrading quantitative accuracy from the sparsity stage is not tested; no ROI statistics, ground-truth reference, or ablation study is supplied to rule out mid-frequency signal leakage into the background image.

Authors: The background term is introduced specifically in the low-resolution stage as a smooth component intended to capture low-frequency non-uniformities, leaving the sparsity-regularized stage to recover higher-frequency quantitative content. This separation follows from the two-stage formulation in our prior work. We acknowledge that no ROI statistics, ground-truth comparisons, or ablation study are provided to quantify potential mid-frequency leakage, as the demonstration uses clinical patient data without reference standards. The visual results are offered as supporting evidence that structures remain intact. revision: no

Circularity Check

0 steps flagged

No circularity; algorithmic extension with single-case empirical demo is self-contained

full rationale

The paper presents a two-stage reconstruction algorithm as an explicit extension of the authors' prior sparsity-regularized method, with the new element being an added background-image estimation term to absorb low-frequency artifacts. No mathematical derivation, uniqueness theorem, or first-principles prediction is offered that reduces to its own inputs by construction; the performance claim is supported solely by qualitative inspection of one patient dataset. Self-citations to the base algorithm are present but do not bear the load of the new claim, which is assessed empirically rather than through any fitted parameter or renamed result. The manuscript itself notes the need for additional validation, confirming the absence of a closed self-referential loop.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

Review performed on abstract only; therefore the ledger is necessarily incomplete and based on explicit statements in the abstract. The background image is introduced as a new modeling component whose independent validation is not shown.

free parameters (1)

regularization weights for sparsity and background estimation
These control the trade-off between data fidelity, sparsity, and artifact absorption; their specific values are not reported and must be chosen or fitted.

axioms (1)

domain assumption Sparsity in a suitable transform domain yields quantitative accuracy in the low-resolution stage
Invoked via reference to prior work; assumed to remain valid when the background component is added.

invented entities (1)

background image no independent evidence
purpose: To absorb low-frequency artifacts that cause image non-uniformity
New modeling variable introduced in the first-stage optimization; no independent evidence (e.g., phantom validation or theoretical guarantee) is provided in the abstract.

pith-pipeline@v0.9.1-grok · 5714 in / 1510 out tokens · 30814 ms · 2026-06-29T10:14:09.759326+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

14 extracted references · 2 canonical work pages · 1 internal anchor

[1]

H. H. Barrett and K. J. Myers, Foundations of Image Science , Wiley series in pure and applied optics, Wiley-Interscience, Hoboken, N.J. (2003)

2003
[2]

E. J. Cand \`e s, J. Romberg, and T. Tao, ``Robust uncertainty principles: E xact signal reconstruction from highly incomplete frequency information,'' IEEE Trans. Inf. Theory 52 (2), 489--509 (2006)

2006
[3]

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. Tech. 14 (2), 119--139 (2006)

2006
[4]

E. Y. Sidky and X. Pan, ``Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization,'' Phys. Med. Biol. 53 (17), 4777--4807 (2008)

2008
[5]

X. Pan, E. Y. Sidky, and M. Vannier, ``Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction?,'' Inv. Prob. 25 (12), art. no. 123009 (2009)

2009
[6]

Zhang, B

Z. Zhang, B. Chen, D. Xia, et al. , ``Directional- TV algorithm for image reconstruction from limited-angular-range data,'' Med. Image Anal. 70 , art. no. 102030 (2021)

2021
[7]

E. Y. Sidky, X. Wu, X. Duan, et al. , ``Accurate volume image reconstruction for digital breast tomosynthesis with directional-gradient and pixel sparsity regularization,'' J Med. Imaging 12 , art. no. S13013 (2025)

2025
[8]

Chambolle and T

A. Chambolle and T. Pock, ``A first-order primal-dual algorithm for convex problems with applications to imaging,'' J. Math. Imaging Vis. 40 , 120--145 (2011)

2011
[9]

E. Y. Sidky, J. H. J rgensen, and X. Pan, ``Convex optimization problem prototyping for image reconstruction in computed tomography with the C hambolle-- P ock algorithm,'' Phys. Med. Biol. 57 , 3065--3091 (2012)

2012
[10]

E. Y. Sidky, J. P. Phillips, Z. Zhang, et al. , ``Efficient primal-dual algorithm for imaging applications with matrix stacking, applied to DBT image reconstruction,'' arXiv preprint arXiv:2604.23063 (2026)

work page internal anchor Pith review Pith/arXiv arXiv 2026
[11]

E. Y. Sidky, J. P. Phillips, W. Zhou, et al. , ``A signal detection model for quantifying overregularization in nonlinear image reconstruction,'' Med. Phys. 48 , 6312--6323 (2021)

2021
[12]

Iyadomi, R

T. Iyadomi, R. Parada, A. Kim, et al. , ``Accelerating low-frequency convergence for limited-angle DBT via two-channel fidelity in PDHG ,'' arXiv preprint arXiv:2603.23955 (2026). Accepted paper to the 9th International Conference on Image Formation in X-Ray Computed Tomography

work page arXiv 2026
[13]

S. D. Rose, A. A. Sanchez, E. Y. Sidky, et al. , ``Investigating simulation-based metrics for characterizing linear iterative reconstruction in digital breast tomosynthesis,'' Med. Phys. 44 (9), e279--e296 (2017)

2017
[14]

A. A. Sanchez, E. Y. Sidky, and X. Pan, ``Use of the hotelling observer to optimize image reconstruction in digital breast tomosynthesis,'' J. Med. Imaging 3 , art. no. 011008 (2016)

2016

[1] [1]

H. H. Barrett and K. J. Myers, Foundations of Image Science , Wiley series in pure and applied optics, Wiley-Interscience, Hoboken, N.J. (2003)

2003

[2] [2]

E. J. Cand \`e s, J. Romberg, and T. Tao, ``Robust uncertainty principles: E xact signal reconstruction from highly incomplete frequency information,'' IEEE Trans. Inf. Theory 52 (2), 489--509 (2006)

2006

[3] [3]

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. Tech. 14 (2), 119--139 (2006)

2006

[4] [4]

E. Y. Sidky and X. Pan, ``Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization,'' Phys. Med. Biol. 53 (17), 4777--4807 (2008)

2008

[5] [5]

X. Pan, E. Y. Sidky, and M. Vannier, ``Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction?,'' Inv. Prob. 25 (12), art. no. 123009 (2009)

2009

[6] [6]

Zhang, B

Z. Zhang, B. Chen, D. Xia, et al. , ``Directional- TV algorithm for image reconstruction from limited-angular-range data,'' Med. Image Anal. 70 , art. no. 102030 (2021)

2021

[7] [7]

E. Y. Sidky, X. Wu, X. Duan, et al. , ``Accurate volume image reconstruction for digital breast tomosynthesis with directional-gradient and pixel sparsity regularization,'' J Med. Imaging 12 , art. no. S13013 (2025)

2025

[8] [8]

Chambolle and T

A. Chambolle and T. Pock, ``A first-order primal-dual algorithm for convex problems with applications to imaging,'' J. Math. Imaging Vis. 40 , 120--145 (2011)

2011

[9] [9]

E. Y. Sidky, J. H. J rgensen, and X. Pan, ``Convex optimization problem prototyping for image reconstruction in computed tomography with the C hambolle-- P ock algorithm,'' Phys. Med. Biol. 57 , 3065--3091 (2012)

2012

[10] [10]

E. Y. Sidky, J. P. Phillips, Z. Zhang, et al. , ``Efficient primal-dual algorithm for imaging applications with matrix stacking, applied to DBT image reconstruction,'' arXiv preprint arXiv:2604.23063 (2026)

work page internal anchor Pith review Pith/arXiv arXiv 2026

[11] [11]

E. Y. Sidky, J. P. Phillips, W. Zhou, et al. , ``A signal detection model for quantifying overregularization in nonlinear image reconstruction,'' Med. Phys. 48 , 6312--6323 (2021)

2021

[12] [12]

Iyadomi, R

T. Iyadomi, R. Parada, A. Kim, et al. , ``Accelerating low-frequency convergence for limited-angle DBT via two-channel fidelity in PDHG ,'' arXiv preprint arXiv:2603.23955 (2026). Accepted paper to the 9th International Conference on Image Formation in X-Ray Computed Tomography

work page arXiv 2026

[13] [13]

S. D. Rose, A. A. Sanchez, E. Y. Sidky, et al. , ``Investigating simulation-based metrics for characterizing linear iterative reconstruction in digital breast tomosynthesis,'' Med. Phys. 44 (9), e279--e296 (2017)

2017

[14] [14]

A. A. Sanchez, E. Y. Sidky, and X. Pan, ``Use of the hotelling observer to optimize image reconstruction in digital breast tomosynthesis,'' J. Med. Imaging 3 , art. no. 011008 (2016)

2016