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.

Uncertainty Analysis in 3D SPECT Reconstruction based on Probabilistic Programming

[Analytical Algorithms Improving the Reconstruction of SPECT Image].

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

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

Uncertainty Propagation from Projections to Region Counts in Tomographic Imaging: Application to Radiopharmaceutical Dosimetry

Probabilistic deconvolution of PET images using informed priors

SVD Reconstruction Algorithm in 3D SPECT Imaging

Simultaneous reconstruction of tissue attenuation and radioactivity maps in SPECT.

Estimating Uncertainty in PET Image Reconstruction via Deep Posterior Sampling

Efficient determination of the uncertainty for the optimization of SPECT system design: a subsampled fisher information matrix

PET image reconstruction: a robust state space approach.

Nonlinear Dual Reconstruction of SPECT Activity and Attenuation Images

An efficient analytical calculation of probability matrix in 2D SPECT

Efficient fully 3-D iterative SPECT reconstruction with Monte Carlo-based scatter compensation

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

Application of Monte Carlo Algorithms to Cardiac Imaging Reconstruction

Spatial Auto-Regressive Analysis of Correlation in 3-D PET With Application to Model-Based Simulation of Data

A Higher-Order Polynomial Method for SPECT Reconstruction.

Probabilistic 3D segmentation for aleatoric uncertainty quantification in full 3D medical data

Effective Emission Tomography Image Reconstruction Algorithms for SPECT Data

A projection-domain low-count quantitative SPECT method for alpha-particle emitting radiopharmaceutical therapy