Abstract:Purpose: We present a probabilistic approach to medical image analysis that requires, and makes use of, explicit prior information provided by a medical expert. Depending on the choice of prior model the method can be used for image enhancement, analysis, and segmentation. Methods: The methodology is based on a probabilistic approach to medical image analysis, that allows integration of 1) arbitrarily complex prior information (for which realizations can be generated), 2) information about a convolution operator of the imaging system, and 3) information about the noise in the reconstructed image into a posterior probability density. The method was demonstrated on positron emission tomography (PET) images obtained from a phantom and a patient with lung cancer. The likelihood model (multivariate log-normal) and the convolution operator were derived from phantom data. Two examples of prior information were used to show the potential of the method. The extended Metropolis-Hastings algorithm, a Markov chain Monte Carlo method, was used to generate realizations of the posterior distribution of the tracer activity concentration. Results: A set of realizations from the posterior was used as the base of a quantitative PET image analysis. The mean and variance of activity concentrations were computed, as well as the probability of high tracer uptake and statistics on the size and activity concentration of high uptake regions. For both phantom and in vivo images, the estimated images of mean activity concentrations appeared to have reduced noise levels, and a sharper outline of high activity regions, as compared to the original PET. The estimated variance of activity concentrations was high at the edges of high activity regions. Conclusions: The methodology provides a probabilistic approach for medical image analysis that explicitly takes into account medical expert knowledge as prior information. The presented first results indicate the potential of the method to improve the detection of small lesions. The methodology allows for a probabilistic measure of the size and activity level of high uptake regions, with possible long-term perspectives for early detection of cancer, as well as treatment, planning, and follow-up.

PET Image Reconstruction: Mean, Variance, and Optimal Minimax Criterion

Robust Framework for PET Image Reconstruction Incorporating System and Measurement Uncertainties.

Uncertainty Penalized Weighted Least Squares Framework For Pet Reconstruction Under Uncertain System Models

PET image reconstruction: a robust state space approach.

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.

Nonlocal Total Variation Based Dynamic PET Image Reconstruction with Low-Rank Constraints

A Robust Least Squares Framework for PET Image Reconstruction

Low dose PET reconstruction with total variation regularization.

Nonlocal Low-Rank and Total Variation Constrained PET Image Reconstruction

Rethinking PET Image Reconstruction: Ultra-Low-Dose, Sinogram and Deep Learning

Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging.

PEt image reconstruction based on particle filter framework

EMISTA-Based Quantitative PET Reconstruction.

3D Tensor Based Nonlocal Low Rank Approximation in Dynamic PET Reconstruction.

Simultaneous Reconstruction and Segmentation of Dynamic Pet Via Low-Rank and Sparse Matrix Decomposition

State-Space Reconstruction of Pet Parametric Maps

Probabilistic deconvolution of PET images using informed priors

Estimation of the Local Statistical Noise in Positron Emission Tomography Revisited: Practical Implementation.

Ordered Subset Expectation Maximization Algorithm for Positron Emission Tomographic Image Reconstruction Using Belief Kernels

Joint Reconstruction of Dynamic PET Activity and Kinetic Parametric Images Using Total Variation Constrained Dictionary Sparse Coding

PET Denoising and Uncertainty Estimation Based on NVAE Model Using Quantile Regression Loss