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.

Non-negative Matrix Factorization using Partial Prior Knowledge for Radiation Dosimetry

Nonnegative-Matrix-Factorization-Based Hyperspectral Unmixing with Partially Known Endmembers

Nonnegative Matrix Factorization with Data-Guided Constraints for Hyperspectral Unmixing.

Unmixing methods based on nonnegativity and weakly mixed pixels for astronomical hyperspectral datasets

L1/2 Sparsity Constrained Nonnegative Matrix Factorization for Hyperspectral Unmixing

Unsupervised Hyperspectral Unmixing Based on Constrained Nonnegative Matrix Factorization and Particle Swarm Optimization

Probabilistic deconvolution of PET images using informed priors

A Complexity Constrained Nonnegative Matrix Factorization for Hyperspectral Unmixing

Nonlinear Hyperspectral Unmixing Based on Sparse Non-Negative Matrix Factorization

Inertia-Constrained Pixel-by-Pixel Nonnegative Matrix Factorisation: A Hyperspectral Unmixing Method Dealing with Intra-Class Variability

Adaptive $textit{L}_{bf 1/2}$ Sparsity-Constrained NMF With Half-Thresholding Algorithm for Hyperspectral Unmixing

Blind Spectral Unmixing in Terahertz Domain Using Nonnegative Matrix Factorization

Mixture prior distributions and Bayesian models for robust radionuclide image processing

Incremental Kernel Non-Negative Matrix Factorization for Hyperspectral Unmixing

A Non-Negative Matrix Factorization Approach For Hyperspectral Unmixing With Partial Known Endmembers

Nonlinear Hyperspectral Unmixing With Robust Nonnegative Matrix Factorization

A Spectral Library and Method for Sparse Unmixing of Hyperspectral Images in Fluorescence Guided Resection of Brain Tumors

Multiple Clustering Guided Nonnegative Matrix Factorization for Hyperspectral Unmixing

Hyperspectral Signal Unmixing Based on Constrained Non-Negative Matrix Factorization Approach

Hyperspectral Unmixing via Nonnegative Matrix Factorization with Handcrafted and Learnt Priors

Minimum distance constrained nonnegative matrix factorization for hyperspectral data unmixing