Abstract:In this work, we present a method for approximating constrained maximum entropy (ME) reconstructions of SPECT data with modifications to a block-iterative maximum a posteriori (MAP) algorithm. Maximum likelihood (ML)-based reconstruction algorithms require some form of noise smoothing. Constrained ME provides a more formal method of noise smoothing without requiring the user to select parameters. In the context of SPECT, constrained ME seeks the minimum-information image estimate among those whose projections are a given distance from the noisy measured data, with that distance determined by the magnitude of the Poisson noise. Images that meet the distance criterion are referred to as feasible images. We find that modeling of all principal degrading factors (attenuation, detector response, and scatter) in the reconstruction is critical because feasibility is not meaningful unless the projection model is as accurate as possible. Because the constrained ME solution is the same as a MAP solution for a particular value of the MAP weighting parameter, beta, the constrained ME solution can be found with a MAP algorithm if the correct value of beta is found. We show that the RBI-MAP algorithm, if used with a dynamic scheme for estimating beta, can approximate constrained ME solutions in 20 or fewer iterations. We compare results for various methods of achieving feasible images on a simulation of Tl-201 cardiac SPECT data. Results show that the RBI-MAP ME approximation provides images and quantitative estimates close to those from a slower algorithm that gives the true ME solution. Also, we find that the ME results have higher spatial resolution and greater high-frequency noise content than a feasibility-based stopping rule, feasibility-based low-pass filtering, and a quadratic Gibbs prior with beta selected according to the feasibility criterion. We conclude that fast ME approximation is possible using either RBI-MAP with the dynamic procedure or a feasibility-based stopping rule, and that such reconstructions may be particularly useful in applications where resolution is critical.

Fast Estimation of Diffusion Tensors under Rician noise by the EM algorithm

An improved EM algorithm for solving MLE in constrained diffusion kurtosis imaging of human brain

Learning Diffusion Model from Noisy Measurement using Principled Expectation-Maximization Method

DiffusionEIT: Diffusion Model for Electrical Impedance Tomography

Diffusion Information Theoretic Learning for Distributed Estimation over Network

$ρ$-Diffusion: A diffusion-based density estimation framework for computational physics

ML-ESTIMATION OF THE T-DISTRIBUTION USING EM AND ITS EXTENSIONS, ECM AND ECME

Fast Diffusion Model For Seismic Data Noise Attenuation

ML ESTIMATION OF THE t DISTRIBUTION USING EM AND ITS EXTENSIONS, ECM AND ECME

Maximum likelihood estimates of diffusion coefficients from single-particle tracking experiments

Fast maximum entropy approximation in SPECT using the RBI-MAP algorithm

Rician Denoising Diffusion Probabilistic Models For Sodium Breast MRI Enhancement

An Expectation-Maximization Algorithm for Training Clean Diffusion Models from Corrupted Observations

Expectation maximization based framework for joint localization and parameter estimation in single particle tracking from segmented images

CDiffMR: Can We Replace the Gaussian Noise with K-Space Undersampling for Fast MRI?

Nonlinear Gradient Field Estimation in Diffusion MRI Tensor Simulation

Improved Patch Denoising Diffusion Probabilistic Models for Magnetic Resonance Fingerprinting

Application of Rician Correction Non-local Mean Algorithm in Diffusion-weighted Image

Regularized EM algorithm for sparse parameter estimation in nonlinear dynamic systems with application to gene regulatory network inference

Highly Accelerated MRI via Implicit Neural Representation Guided Posterior Sampling of Diffusion Models

Statistical Spatially Inhomogeneous Diffusion Inference