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.

Reconstructing images from projections using the maximum-entropy method. Numerical simulations of low-aspect astrotomography

Minimum Cross-Entropy Algorithm for Image Reconstruction from Incomplete Projection

New-generation maximum entropy method: a Lagrangian-based algorithm for dynamic reconstruction of interferometric data

Maximum-entropy image reconstruction using wavelets

New-generation Maximum Entropy Method (ngMEM): a Lagrangian-based algorithm for dynamic reconstruction of interferometric data

On the Image Reconstruction of Capacitively Coupled Electrical Resistance Tomography (CCERT) with Entropy Priors

Application of the maximum entropy method to spectral-domain optical coherence tomography for enhancing axial resolution.

High resolution VLBI polarisation imaging of AGN with the Maximum Entropy Method

N-dimensional maximum-entropy tomography via particle sampling

A Comparison of the MEM and the MP-CLEAN Methods for the Image Reconstruction of Extended Sources with Simulated Visibility Data

Shape-enhanced Maximum Intensity Projection.

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

High-dimensional maximum-entropy phase space tomography using normalizing flows

Maximum a Posteriori Reconstruction Algorithms in Emission CT

MEM and CLEAN Imaging of VLBA Polarisation Observations of Compact Active Galactic Nuclei

Maximum-entropy reconstruction of gravitational lenses using shear and/or magnification data

A simulation study on image reconstruction in magnetic particle imaging with field-free-line encoding

Tomographic reconstruction from Poisson distributed data: a fast and convergent EM-TV dual approach

Closed form solution of the maximum entropy equations with application to fast radio astronomical image formation

Image Reconstruction Based on Sequential Monte Carlo Principle for EMT

Image Reconstruction Using a Mixture Score Function (MSF)