1. Uncertainty quantification of CT regularized reconstruction within the Bayesian framework
- Author
-
Negin Khoeiniha, Patricio Guerrero, and Wim Dewulf
- Subjects
Technology - Abstract
Computed Tomography (CT) reconstruction is an important inverse problem in industrial imaging, requiring robust methods to address different sources of error in the data and model. Among the various reconstruction approaches that tackle different challenges in CT modeling [1], such as limitations in the data, statistical methods are known for their ability to model various properties [2, 3] of the input data and assess their impact on the reconstruction outcome [4]. These methods can also incorporate regularization techniques, which, for instance, force smoothness in the solution by penalizing heavy oscillations in pixel values. Given our interest in the uncertainty Quantification (UQ) of CT reconstruction, we employ statistical methods within the Bayesian framework for solving this inverse problem. First, we define a model for CT reconstruction and formulate a statistical framework for the problem. We then reconstruct our data using appropriate methods and quantify the uncertainty in the results. Our primary focus is on the effect of noise on the reconstruction and the corresponding uncertainty it induces. Additionally, we aim to utilize computationally feasible sampling techniques to analyze the distribution of the solution, enabling an in-depth evaluation of the results. To achieve these goals, we apply a rapid regularized Markov Chain Monte Carlo (MCMC) reconstruction method [4, 5], employing the Metropolis-Adjusted Langevin Algorithm (MALA) [6] and its Lipschitz-adaptive variant (LipMALA) [7]. This approach produces a volumetric model where each voxel is represented by a probability distribution, which can be transformed into a triplet of gray-value models: one for the central value and two for the bounds of the confidence interval. Bi-directional and uni-directional length measurements derived from these gray-value models applied to real CT data yield task-specific measurement uncertainties. This method significantly reduces computational and storage demands compared to classic Monte Carlo simulations while incorporating regularization techniques. Experimental results using aluminum cylindrical step gauge [8] data acquired with a Nikon 225 CT system validate our approach. This work provides a scalable, statistically grounded methodology for UQ in CT reconstruction, offering enhanced reliability for industrial applications.
- Published
- 2025
- Full Text
- View/download PDF