Abstract:The diagnosis of medical conditions and subsequent treatment often involves radionuclide imaging techniques. To refine localisation accuracy and improve diagnostic confidence, compared with the use of a single scanning technique, a combination of two (or more) techniques can be used but with a higher risk of misalignment. For this to be reliable and accurate, recorded data undergo processing to suppress noise and enhance resolution. A step in image processing techniques for such inverse problems is the inclusion of smoothing. Standard approaches, however, are usually limited to applying identical models globally. In this study, we propose a novel Laplace and Gaussian mixture prior distribution that incorporates different smoothing strategies with the automatic model-based estimation of mixture component weightings creating a locally adaptive model. A fully Bayesian approach is presented using multi-level hierarchical modelling and Markov chain Monte Carlo (MCMC) estimation methods to sample from the posterior distribution and hence perform estimation. The proposed methods are assessed using simulated γ -eye TM camera images and demonstrate greater noise reduction than existing methods but without compromising resolution. As well as image estimates, the MCMC methods also provide posterior variance estimates and hence uncertainty quantification takes into consideration any potential sources of variability. The use of mixture prior models, part Laplace random field and part Gaussian random field, within a Bayesian modelling approach is not limited to medical imaging applications but provides a more general framework for analysing other spatial inverse problems. Locally adaptive prior distributions provide a more realistic model, which leads to robust results and hence more reliable decision-making, especially in nuclear medicine. They can become a standard part of the toolkit of everyone working in image processing applications.

Momentum-Space Renormalization Group Transformation in Bayesian Image Modeling by Gaussian Graphical Model

Learning Gaussian Networks

Statistical analysis of the expectation-maximization algorithm with loopy belief propagation in Bayesian image modeling

Inverse Renormalization Group Transformation in Bayesian Image Segmentations

Sublinear Computational Time Modeling by Momentum-Space Renormalization Group Theory in Statistical Machine Learning Procedures

Efficient Bayesian Computation by Proximal Markov Chain Monte Carlo: When Langevin Meets Moreau

Multivariate, Heteroscedastic Empirical Bayes via Nonparametric Maximum Likelihood

Mixtures of conditional Gaussian scale mixtures applied to multiscale image representations

A Non-Gaussian Bayesian Filter Using Power and Generalized Logarithmic Moments

Mixture prior distributions and Bayesian models for robust radionuclide image processing

Learning Gaussian mixture model with a maximization-maximization algorithm for image classification

Bayesian Generalized Kernel Mixed Models

Partial Gaussian Graphical Model Estimation

A Multivariate Non-Gaussian Bayesian Filter Using Power Moments

Bayesian Metric Learning for Uncertainty Quantification in Image Retrieval

The Bridged Posterior: Optimization, Profile Likelihood and a New Approach to Generalized Bayes

Bayesian Non-linear Latent Variable Modeling via Random Fourier Features

Accelerated Bayesian Imaging by Relaxed Proximal-Point Langevin Sampling

Markov Chain Monte Carlo for Automated Face Image Analysis

Parameter Estimation For Multivariate Generalized Gaussian Distributions