Uncertainty Quantification of Inclusion Boundaries in the Context of X-Ray Tomography

DOI10.1137/21M1433782zbMath1524.65022arXiv2107.06607OpenAlexW3177689450MaRDI QIDQ6109154

Per Christian Hansen, Yiqiu Dong, Babak Maboudi Afkham

Publication date: 30 June 2023

Published in: SIAM/ASA Journal on Uncertainty Quantification (Search for Journal in Brave)

Abstract: In this work, we describe a Bayesian framework for reconstructing the boundaries of piecewise smooth regions in the X-ray computed tomography (CT) problem in an infinite-dimensional setting. In addition to the reconstruction, we are also able to quantify the uncertainty of the predicted boundaries. Our approach is goal oriented, meaning that we directly detect the discontinuities from the data, instead of reconstructing the entire image. This drastically reduces the dimension of the problem, which makes the application of Markov Chain Monte Carlo (MCMC) methods feasible. We show that our method provides an excellent platform for challenging X-ray CT scenarios (e.g., in case of noisy data, limited angle, or sparse angle imaging). We investigate the performance and accuracy of our method on synthetic data as well as on real-world data. The numerical results indicate that our method provides an accurate method in detecting boundaries of piecewise smooth regions and quantifies the uncertainty in the prediction.

Full work available at URL: https://arxiv.org/abs/2107.06607

zbMATH Keywords

inverse problemsBayesian frameworkX-ray CTgoal-oriented UQWhittle-Matérn field

Mathematics Subject Classification ID

Probabilistic models, generic numerical methods in probability and statistics (65C20) Monte Carlo methods (65C05) Biomedical imaging and signal processing (92C55) Applications of Markov chains and discrete-time Markov processes on general state spaces (social mobility, learning theory, industrial processes, etc.) (60J20) Radon transform (44A12)

Cites Work

• Unnamed Item
• Unnamed Item
• Unnamed Item
• Unnamed Item
• An Explicit Link between Gaussian Fields and Gaussian Markov Random Fields: The Stochastic Partial Differential Equation Approach
• The Bayesian formulation of EIT: analysis and algorithms
• Besov priors for Bayesian inverse problems
• Whittle-Matérn priors for Bayesian statistical inversion with applications in electrical impedance tomography
• A Mumford-Shah level-set approach for the inversion and segmentation of X-ray tomography data
• Hierarchical Bayesian level set inversion
• Inference from iterative simulation using multiple sequences
• Rank-normalization, folding, and localization: an improved \(\widehat{R}\) for assessing convergence of MCMC (with Discussion)
• Randomized approaches to accelerate MCMC algorithms for Bayesian inverse problems
• Hyperparameter estimation in Bayesian MAP estimation: parameterizations and consistency
• Discretization-invariant Bayesian inversion and Besov space priors
• Cauchy difference priors for edge-preserving Bayesian inversion
• The Mathematics of Computerized Tomography
• Inverse problems: A Bayesian perspective
• An Analysis of Infinite Dimensional Bayesian Inverse Shape Acoustic Scattering and Its Numerical Approximation
• Centrosymmetric (Cross-Symmetric) Matrices, Their Basic Properties, Eigenvalues, and Eigenvectors
• Local Tomography
• Local Tomography II
• Active contours without edges
• X-ray Computed Tomography using a Sparsity Enforcing Prior Model Based on Haar Transformation in a Bayesian Framework
• A joint reconstruction and lambda tomography regularization technique for energy-resolved x-ray imaging
• Advanced Lectures on Machine Learning
• Inconsistent Estimation and Asymptotically Equal Interpolations in Model-Based Geostatistics
• Lower-bound gradient estimates for first-order Hamilton-Jacobi equations and applications to the regularity of propagating fronts
• MCMC methods for functions: modifying old algorithms to make them faster

Related Items (5)

A Bayesian approach for consistent reconstruction of inclusions ⋮ CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python ⋮ Object based Bayesian full-waveform inversion for shear elastography ⋮ Tomographic inverse problems: mathematical challenges and novel applications. Abstracts from the workshop held April 30 -- May 5, 2023 ⋮ Uncertainty quantification for goal-oriented inverse problems via variational encoder-decoder networks

This page was built for publication: Uncertainty Quantification of Inclusion Boundaries in the Context of X-Ray Tomography