Alternating minimization algorithm with automatic relevance determination for transmission tomography under Poisson noise

DOI10.1137/141000038zbMATH Open1330.94009arXiv1412.8464OpenAlexW2963981262MaRDI QIDQ3454486FDOQ3454486

Authors: Yan Kaganovsky, Shaobo Han, Soysal Degirmenci, David G. Politte, David Brady, Joseph A. O'Sullivan, Lawrence Carin 

Publication date: 25 November 2015

Published in: SIAM Journal on Imaging Sciences (Search for Journal in Brave)

Abstract: We propose a globally convergent alternating minimization (AM) algorithm for image reconstruction in transmission tomography, which extends automatic relevance determination (ARD) to Poisson noise models with Beer's law. The algorithm promotes solutions that are sparse in the pixel/voxel-differences domain by introducing additional latent variables, one for each pixel/voxel, and then learning these variables from the data using a hierarchical Bayesian model. Importantly, the proposed AM algorithm is free of any tuning parameters with image quality comparable to standard penalized likelihood methods. Our algorithm exploits optimization transfer principles which reduce the problem into parallel 1D optimization tasks (one for each pixel/voxel), making the algorithm feasible for large-scale problems. This approach considerably reduces the computational bottleneck of ARD associated with the posterior variances. Positivity constraints inherent in transmission tomography problems are also enforced. We demonstrate the performance of the proposed algorithm for x-ray computed tomography using synthetic and real-world datasets. The algorithm is shown to have much better performance than prior ARD algorithms based on approximate Gaussian noise models, even for high photon flux.

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

Recommendations

• Tomographic reconstruction from Poisson distributed data: a fast and convergent EM-TV dual approach
• Accelerated alternating minimization algorithm for Poisson noisy image recovery
• Preconditioned alternating projection algorithms for maximum \textit{a posteriori} ECT reconstruction
• A fast algorithm for deconvolution and Poisson noise removal
• An algorithm for hybrid regularizers based image restoration with Poisson noise.

zbMATH Keywords

Poisson noisealternating minimizationtransmission tomographyautomatic relevance determinationoptimization transferx-ray CT

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) Convex programming (90C25) Nonconvex programming, global optimization (90C26) Nonlinear programming (90C30) Ill-posedness and regularization problems in numerical linear algebra (65F22) Image processing (compression, reconstruction, etc.) in information and communication theory (94A08) Numerical solution to inverse problems in abstract spaces (65J22)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Elements of Information Theory
• On the convergence properties of the EM algorithm
• Title not available (Why is that?)
• Pattern recognition and machine learning.
• Hierarchical modeling and analysis for spatial data
• 10.1162/15324430152748236
• Linear and nonlinear programming.
• Robust regression using iteratively reweighted least-squares
• Principles of computerized tomography imaging
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• A probing method for computing the diagonal of a matrix inverse.
• Bayesian learning for neural networks
• An estimator for the diagonal of a matrix
• Newton’s Method with a Model Trust Region Modification
• Large scale Bayesian inference and experimental design for sparse linear models
• Latent Variable Bayesian Models for Promoting Sparsity
• Convergence theorems for generalized alternating minimization procedures
• Fast Model-Based X-Ray CT Reconstruction Using Spatially Nonhomogeneous ICD Optimization

Cited In (3)

• Accelerated alternating minimization algorithm for Poisson noisy image recovery
• A combined first and fractional order regularization method for mixed Poisson-White spike noisy image restoration
• A variational method for Abel inversion tomography with mixed Poisson-Laplace-Gaussian noise

Uses Software

• Matlab
• spBayes
• VARD
• PRMLT

This page was built for publication: Alternating minimization algorithm with automatic relevance determination for transmission tomography under Poisson noise

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3454486)