Bayesian Joint Reconstruction and Uncertainty Quantification for Limited-Data Multi-Contrast X-Ray CT

Abstract

Multi-contrast X-ray computed tomography (CT) combines attenuation contrast with phase and dark-field contrast, which are sensitive to refraction and ultra-small-angle scattering, respectively. Common implementations record several intensity measurements per projection angle by phase or mask stepping, which lengthens acquisition and, at a fixed exposure per measurement, raises the dose. Limited-data acquisition via sparse angular sampling or reduced stepping lowers these costs but makes the data less informative, calling for prior information and uncertainty quantification. We study joint Bayesian reconstruction of attenuation, phase, and dark-field images from separately retrieved sinograms. A joint total variation term in the prior couples edge locations across the images but not their values. Posterior sampling with the No-U-Turn Sampler complements maximum a posteriori estimates with pixel-wise uncertainty. We compare this joint prior with channel-wise total variation alone on synthetic and experimental Talbot-Lau interferometry data under full, sparse-angle, sparse-step, and combined acquisition (up to 24 times fewer exposures). Both maximum a posteriori reconstructions suppressed undersampling artifacts, with synthetic relative errors 1.6-11 times lower than those of filtered backprojection. The joint prior gave modest synthetic accuracy gains, a cleaner experimental dark-field background, 23-30 % lower mean dark-field posterior standard deviation, and 4-38% higher median effective sample size in experimental limited-data runs. In the synthetic study, posterior standard deviations matched or exceeded the errors in flat regions but underestimated them at boundaries, especially in phase. Experimental phase chains mixed slowly under angular undersampling. The formulation also applies to edge illumination with adapted forward and noise models.

Bayesian bilevel operator learning with low-rank adaptation for efficient uncertainty quantification of PDE inverse problems.
Ray Zirui Zhang, Christopher E Miles, Xiaohui Xie +1
Nature communications19 Sep 2026Peer-reviewedCC BYFull textStatistical & nonlinear physicsEurope PMC
Reconstruction of Unobserved Core Electron Density in Tokamak Plasmas with a Neural Operator Transport Surrogate and Uncertainty Quantification
D. Sorokin, A. Granovskiy, M. Stokolesov +3
arXiv17 Oct 2026PreprintCC BY-NC-SA 4.0Full textPlasma physicsarXiv
A differentiable PCA-GPR surrogate framework for full-transient Bayesian inverse uncertainty quantification of black-box nuclear thermal-hydraulic systems
Zikang Yaoa, Eugene Yeea, Fue-Sang Lien
arXiv14 Oct 2026PreprintarXiv licenceFull textComputational physics Fluid dynamicsarXiv
Rapid Uncertainty Quantification on a Latent Field using Fisher Information
Karl Daningburg, A. E. Lovell, Arvind T. Mohan +1
arXiv29 Sep 2026PreprintarXiv licenceFull textComputational physics Nuclear physicsarXiv
OpenTURNS: A Python package for uncertainty quantification
Michaël Baudin, Anne Dutfoy, Régis Lebrun +6
Journal of Open Source Software2 Sep 2026Peer-reviewedCC BY 4.0Full textGeneral mathematics General physicsHAL

Cite

BibTeX RIS

BibTeX

@misc{yazdanian2026bayesian,
  title = {{Bayesian Joint Reconstruction and Uncertainty Quantification for Limited-Data Multi-Contrast X-Ray CT}},
  author = {Hassan Yazdanian and Henrik Mäkinen},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.21406},
  eprint = {2610.21406},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.21406},
  note = {CC BY 4.0},
}

RIS

TY  - UNPB
TI  - Bayesian Joint Reconstruction and Uncertainty Quantification for Limited-Data Multi-Contrast X-Ray CT
AU  - Hassan Yazdanian
AU  - Henrik Mäkinen
PY  - 2026
DA  - 2026/10/15
PB  - arXiv
DO  - 10.48550/arxiv.2610.21406
UR  - https://arxiv.org/abs/2610.21406
AB  - Multi-contrast X-ray computed tomography (CT) combines attenuation contrast with phase and dark-field contrast, which are sensitive to refraction and ultra-small-angle scattering, respectively. Common implementations record several intensity measurements per projection angle by phase or mask stepping, which lengthens acquisition and, at a fixed exposure per measurement, raises the dose. Limited-data acquisition via sparse angular sampling or reduced stepping lowers these costs but makes the data less informative, calling for prior information and uncertainty quantification. We study joint Bayesian reconstruction of attenuation, phase, and dark-field images from separately retrieved sinograms. A joint total variation term in the prior couples edge locations across the images but not their values. Posterior sampling with the No-U-Turn Sampler complements maximum a posteriori estimates with pixel-wise uncertainty. We compare this joint prior with channel-wise total variation alone on synthetic and experimental Talbot-Lau interferometry data under full, sparse-angle, sparse-step, and combined acquisition (up to 24 times fewer exposures). Both maximum a posteriori reconstructions suppressed undersampling artifacts, with synthetic relative errors 1.6-11 times lower than those of filtered backprojection. The joint prior gave modest synthetic accuracy gains, a cleaner experimental dark-field background, 23-30 % lower mean dark-field posterior standard deviation, and 4-38% higher median effective sample size in experimental limited-data runs. In the synthetic study, posterior standard deviations matched or exceeded the errors in flat regions but underestimated them at boundaries, especially in phase. Experimental phase chains mixed slowly under angular undersampling. The formulation also applies to edge illumination with adapted forward and noise models.
ER  -

Abstract

Related

Cite