Abstract:System designs in single photon emission tomography (SPECT) can be evaluated based on the fundamental trade-off between bias and variance that can be achieved in the reconstruction of emission tomograms. This trade off can be derived analytically using the Cramer-Rao type bounds, which imply the calculation and the inversion of the Fisher information matrix (FIM). The inverse of the FIM expresses the uncertainty associated to the tomogram, enabling the comparison of system designs. However, computing, storing and inverting the FIM is not practical with 3-D imaging systems. In order to tackle the problem of the computational load in calculating the inverse of the FIM, a method based on the calculation of the local impulse response and the variance, in a single point, from a single row of the FIM, has been previously proposed for system design. However this approximation (circulant approximation) does not capture the global interdependence between the variables in shift-variant systems such as SPECT, and cannot account e.g., for data truncation or missing data. Our new formulation relies on subsampling the FIM. The FIM is calculated over a subset of voxels arranged in a grid that covers the whole volume. Every element of the FIM at the grid points is calculated exactly, accounting for the acquisition geometry and for the object. This new formulation reduces the computational complexity in estimating the uncertainty, but nevertheless accounts for the global interdependence between the variables, enabling the exploration of design spaces hindered by the circulant approximation. The graphics processing unit accelerated implementation of the algorithm reduces further the computation times, making the algorithm a good candidate for real-time optimization of adaptive imaging systems. This paper describes the subsampled FIM formulation and implementation details. The advantages and limitations of the new approximation are explored, in comparison with the circulant approximation, in the context of design optimization of a parallel-hole collimator SPECT system and of an adaptive imaging system (similar to the commercially available D-SPECT).

Efficient Calculation of Resolution and Variance in 2-D Circular-Orbit SPECT

[Analytical Algorithms Improving the Reconstruction of SPECT Image].

Short-Scan Spect Imaging With Non-Uniform Attenuation And 3d Distance-Dependent Spatial Resolution

Solving Outside-Axial-Field-Of-View Scatter Correction Problem in Pet Via Digital Experimentation

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

An efficient analytical calculation of probability matrix in 2D SPECT

A Class of Analytical Methods That Compensate for Attenuation and Spatially-Variant Resolution in 2D SPECT

Use of Monte Carlo Simulations for the Scatter Events Analysis of PET Scanners

Evaluation of imaging protocol for ECT based on CS image reconstruction algorithm

Image Reconstruction in 3D Short-Scan SPECT

/Spl Pi/-Scheme Short-Scan SPECT and Image Reconstruction with Nonuniform Attenuation

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

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

Simultaneous reconstruction of tissue attenuation and radioactivity maps in SPECT.

Low dose PET reconstruction with total variation regularization.

Effective Compensation for Physical Effects in 3d Single-Photon Emission Computed Tomography

Toward Optimal Noniterative Reconstruction for 3D SPECT with Uniform Attenuation and Distance-Dependent Spatial Resolution

Interior SPECT- Exact and Stable ROI Reconstruction from Uniformly Attenuated Local Projections.

Π-Scheme Short-Scan SPECT and Image Reconstruction

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

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