Fully 3D PET Image Reconstruction Using A Fourier Preconditioned Conjugate-Gradient Algorithm

Since the data sixes in fully 3D PET imaging are very large, iterative image reconstruction algorithms must converge in very few iterations to be useful. One can improve the convergence rate of the conjugate-gradient (CG) algorithm by incorporating preconditioning operators that approximate the inverse of the Hessian of the objective function. If the 3D cylindrical PET geometry were not truncated at the ends, then the Hessian of the penalized least-squares objective function would be approximately shift-invariant, i.e. G'G would be nearly block-circulant, where G is the system matrix. The authors propose a Fourier preconditioner based on this shift-invariant approximation to the Hessian. Results show that this preconditioner significantly accelerates the convergence of the CG algorithm with only a small increase in computation.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86015/1/Fessler139.pd

Maximum-Likelihood Transmission Image Reconstruction for Overlapping Transmission Beams

In many transmission imaging geometries, the transmitted "beams" of photons overlap on the detector, such that a detector element may record photons that originated in different sources or source locations and thus traversed different paths through the object. Examples include systems based on scanning line sources or on multiple parallel rod sources. The overlap of these beams has been disregarded by both conventional analytical reconstruction methods as well as by previous statistical reconstruction methods. The authors propose a new algorithm for statistical image reconstruction of attenuation maps that explicitly accounts for overlapping beams in transmission scans. The algorithm is guaranteed to monotonically increase the objective function at each iteration. The availability of this algorithm enables the possibility of deliberately increasing the beam overlap so as to increase count rates. Simulated single photon emission tomography transmission scans based on a multiple line source array demonstrate that the proposed method yields improved resolution/noise tradeoffs relative to "conventional" reconstruction algorithms, both statistical and nonstatistical.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85818/1/Fessler78.pd

Maximum Likelihood Transmission Image Reconstruction for Over lapping Transmission Beams

In many transmission imaging geometries, the transmitted “beams” of photons overlap on the detector, such that a detector element may record photons that originated in different sources or source locations and thus traversed different paths through the object, Examples include systems based on scanning line sources or on multiple parallel rod sources. The overlap of these beams has been disregarded by both conventional analytical reconstruction methods as well as by previous statistical reconstruction methods. We propose a new algorithm for statistical image reconstruction of attenuation maps that explicitly accounts for overlapping beams in transmission scans. The algorithm is guaranteed to monotonically increase the objective function at each iteration. The availability of this algorithm enables the possibility of deliberately increasing the beam overlap so as to increase count rates. Simulated SPECT transmission scans based on a multiple line source array demonstrate that the proposed method yields improved resolution/noise tradeoffs relative to “conventional” reconstruction algorithms, both statistical and nonstatistical.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85817/1/Fessler156.pd

Grouped-Coordinate Ascent Algorithms for Penalized-Likelihood Transmission Image Reconstruction

Presents a new class of algorithms for penalized-likelihood reconstruction of attenuation maps from low-count transmission scans. We derive the algorithms by applying to the transmission log-likelihood a version of the convexity technique developed by De Pierro for emission tomography. The new class includes the single-coordinate ascent (SCA) algorithm and Lange's convex algorithm for transmission tomography as special cases. The new grouped-coordinate ascent (GCA) algorithms in the class overcome several limitations associated with previous algorithms. (1) Fewer exponentiations are required than in the transmission maximum likelihood-expectation maximization (ML-EM) algorithm or in the SCA algorithm. (2) The algorithms intrinsically accommodate nonnegativity constraints, unlike many gradient-based methods. (3) The algorithms are easily parallelizable, unlike the SCA algorithm and perhaps line-search algorithms. We show that the GCA algorithms converge faster than the SCA algorithm, even on conventional workstations. An example from a low-count positron emission tomography (PET) transmission scan illustrates the method.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86021/1/Fessler93.pd

Deadtime correction for two multihead Anger cameras in 131I dual‐energy‐window‐acquisition mode

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134914/1/mp8162.pd

Fast Parallelizable Algorithms for Transmission Image Reconstruction

Presents a new class of algorithm for penalized-likelihood reconstruction of attenuation maps from low-count transmission scans. The authors derive the algorithms by applying to the transmission log-likelihood a variation of the convexity technique developed by De Pierro for the emission case. The new algorithms overcome several limitations associated with previous algorithms. (1) Fewer exponentiations are required than in the transmission EM algorithm or in coordinate-ascent algorithms. (2) The algorithms intrinsically accommodate nonnegativity constraints, unlike many gradient-based methods. (3) The algorithms are easily parallelizable, unlike coordinate-ascent algorithms and perhaps line-search algorithms. The authors show that the algorithms converge faster than several alternatives, even on conventional workstations. They give examples from low-count PET transmission scans and from truncated fan-beam SPECT transmission scans.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86006/1/Fessler136.pd

The consequences of a new software package for the quantification of gated-SPECT myocardial perfusion studies

Semiquantitative analysis of myocardial perfusion scintigraphy (MPS) has reduced inter- and intraobserver variability, and enables researchers to compare parameters in the same patient over time, or between groups of patients. There are several software packages available that are designed to process MPS data and quantify parameters. In this study the performances of two systems, quantitative gated SPECT (QGS) and 4D-MSPECT, in the processing of clinical patient data and phantom data were compared. The clinical MPS data of 148 consecutive patients were analysed using QGS and 4D-MSPECT to determine the end-diastolic volume, end-systolic volume and left ventricular ejection fraction. Patients were divided into groups based on gender, body mass index, heart size, stressor type and defect type. The AGATE dynamic heart phantom was used to provide reference values for the left ventricular ejection fraction. Although the correlations were excellent (correlation coefficients 0.886 to 0.980) for all parameters, significant differences (p < 0.001) were found between the systems. Bland-Altman plots indicated that 4D-MSPECT provided overall higher values of all parameters than QGS. These differences between the systems were not significant in patients with a small heart (end-diastolic volume < 70 ml). Other clinical factors had no direct influence on the relationship. Additionally, the phantom data indicated good linear responses of both systems. The discrepancies between these software packages were clinically relevant, and influenced by heart size. The possibility of such discrepancies should be taken into account when a new quantitative software system is introduced, or when multiple software systems are used in the same institution.Vascular Biology and Interventio

In vivo mapping of cholinergic terminals in normal aging, Alzheimer's disease, and Parkinson's disease

To map presynaptic cholinergic terminal densities in normal aging (n = 36), Alzheimer's disease (AD) (n = 22), and Parkinson's disease (PD) (n = 15), we performed single-photon emission computed tomography using [ 123 I]iodoben-zovesamicol (IBVM), an in vivo marker of the vesicular acetylcholine transporter. We used coregistered positron emission tomography with [ 18 F]fluorodexyglucose for metabolic assessment and coregistered magnetic resonance imaging for atrophy assessment. In controls (age, 22–91 years), cortical IBVM binding declined only 3.7% per decade. In AD, cortical binding correlated inversely with dementia severity. In mild dementia, binding differed according to age of onset, but metabolism did not. With an onset age of less than 65 years, binding was reduced severely throughout the entire cerebral cortex and hippocapus (about 30%), but with an onset age of 65 years or more, binding reductions were restricted to temporal cortex and hippocampus. In PD without dementia, binding was reduced only in parietal and occipital cortex, but demented PD subjects had extensive cortical binding decreases similar to early-onset AD. We conclude that cholinergic neuron integrity can be monitored in living AD and PD patients, and that it is not so devastated in vivo as suggested by postmortem choline acetylransferase activity (50–80%).Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/50361/1/410400309_ftp.pd

A study of SPECT/CT camera stability for quantitative imaging

BACKGROUND: The purpose of this study was twofold: to evaluate the quantitative stability of a SPECT/CT gamma camera over time and to determine if daily flood acquisitions can reliably serve as calibration factors for quantitative SPECT. Using a cylindrical water phantom filled with measured amounts of (99m)Tc, factors were calculated to convert counts/cc to activity/cps. Measurements were made over an 18-month period. System sensitivity data calculated from (57)Co daily quality assurance (DQA) flood acquisitions were then compared to the (99m)Tc calibration factors to determine the relationship of the factors. RESULTS: The coefficient of variation is 2.7 % for the (99m)Tc cylinder conversion factors and 2.6 % for the (57)Co DQA flood data. The greatest difference between the cylinder conversion factors and the flood data is less than 3 %. CONCLUSIONS: Based on the results, the camera was stable within 3 % over an 18-month time period. The daily flood source acquisitions can be a reliable source for tracking camera stability and may provide information on updating the calibration factor for quantitative imaging

Small average differences in attenuation corrected images between men and women in myocardial perfusion scintigraphy: a novel normal stress database

<p>Abstract</p> <p>Background</p> <p>The American Society of Nuclear Cardiology and the Society of Nuclear Medicine state that incorporation of attenuation-corrected (AC) images in myocardial perfusion scintigraphy (MPS) will improve image quality, interpretive certainty, and diagnostic accuracy. However, commonly used software packages for MPS usually include normal stress databases for non-attenuation corrected (NC) images but not for attenuation-corrected (AC) images. The aim of the study was to develop and compare different normal stress databases for MPS in relation to NC vs. AC images, male vs. female gender, and presence vs. absence of obesity. The principal hypothesis was that differences in mean count values between men and women would be smaller with AC than NC images, thereby allowing for construction and use of gender-independent AC stress database.</p> <p>Methods</p> <p>Normal stress perfusion databases were developed with data from 126 male and 205 female patients with normal MPS. The following comparisons were performed for all patients and separately for normal weight vs. obese patients: men vs. women for AC; men vs. women for NC; AC vs. NC for men; and AC vs. NC for women.</p> <p>Results</p> <p>When comparing AC for men vs. women, only minor differences in mean count values were observed, and there were no differences for normal weight vs. obese patients. For all other analyses major differences were found, particularly for the inferior wall.</p> <p>Conclusions</p> <p>The results support the hypothesis that it is possible to use not only gender independent but also weight independent AC stress databases.</p