Review of convergent beam tomography in single photon emission computed tomography.

Investigation of convergent-beam single photon emission computed tomography (SPECT) is actively being pursued to evaluate its clinical potentials. Fan-beam, cone-beam, pin-hole and astigmatic collimators are being used with rotating gamma cameras having large crystal areas, to increase the sensitivity for emission and transmission computed tomography of small organs such as the thyroid, brain or heart. With new multi-detector SPECT systems, convergent-beam geometry offers the ability to simultaneously obtain emission and transmission data necessary to quantify uptake of radiopharmaceutical distributions in the heart. The development of convergent-beam geometry in SPECT requires the integration of hardware and software. In considering hardware, the optimum detector system for cone-beam tomography is a system that satisfies the data sufficiency condition for which the scanning trajectory intersects any plane passing through the reconstructed region of interest. However, the major development of algorithms has been for the data insufficient case of single planar orbit acquisitions. The development of these algorithms have made possible the preliminary evaluation of this technology and the imaging of brain and heart are showing significant potential for the clinical application of cone-beam tomography. Presently, significant research activity is pursuing the development of algorithms for data acquisitions that satisfy the data sufficiency condition and that can be implemented easily and inexpensively on clinical SPECT systems.     

Early exercise thallium-201 single photon emission computed tomography in unstable angina: a prospective study.

To evaluate the safety and diagnostic value of early symptom-limited exercise electrocardiography (ECG) and exercise thallium-201 single photon emission computed tomography (SPECT) in unstable angina (UA), 39 patients were studied prospectively soon after stabilization on medical treatment. No patient had a history of myocardial infarction (MI) or revascularization and patients with left bundle branch block were excluded. Exercise ECG and exercise thallium-201 SPECT were performed 8 +/- 4 days and 11 +/- 6 days respectively after admission to hospital. Seventeen out of 39 patients (44%) had positive exercise ECGs and 22 (56%) negative or inconclusive ones. Exercise thallium-201 SPECT was positive in 26 patients (67%) and negative in 13 patients (33%). Thirty-one patients underwent coronary arteriography and 24 of them proved to have significant coronary artery disease (CAD). The sensitivity, specificity and positive predictive value of exercise ECG in detecting CAD are 62%, 86%, and 94% respectively while the corresponding results are 96%, 100%, and 100% for exercise thallium-201 SPECT. Therefore, it is concluded that the early symptom-limited exercise test is safe in medically stabilized patients with UA. Early exercise thallium-201 SPECT is highly sensitive and predictive of the presence of significant CAD among patients in the early recovery phase of UA and can be used in selecting this group of patients for coronary angiography and other therapeutic strategies.     

Rapid calculation of detectability in Bayesian single photon emission computed tomography

We consider the calculation of lesion detectability using a mathematical model observer, the channelized Hotelling observer (CHO), in a signal-known-exactly/background-known-exactly detection task for single photon emission computed tomography (SPECT). We focus on SPECT images reconstructed with Bayesian maximum a posteriori methods. While model observers are designed to replace time-consuming studies using human observers, the calculation of CHO detectability is usually accomplished using a large number of sample images, which is still time consuming. We develop theoretical expressions for a measure of detectability, the signal-to-noise-ratio (SNR) of a CHO observer, that can be very rapidly evaluated. Key to our expressions are approximations to the reconstructed image covariance. In these approximations, we use methods developed in the PET literature, but modify them to reflect the different nature of attenuation and distance-dependent blur in SPECT. We validate our expressions with Monte Carlo methods. We show that reasonably accurate estimates of the SNR can be obtained at a computational expense equivalent to approximately two projection operations, and that evaluating SNR for subsequent lesion locations requires negligible additional computation.     

Convolutional image reconstruction for quantitative single photon emission computed tomography

Two image reconstruction algorithms have been investigated. They are based on filtered backprojection, and are useful when the tissue attenuation is considered to be uniform in the object. The first method uses a weighted backprojection, the weighting factor being determined in such a way that the photon attenuation is compensated with low noise propagation. The parameters involved in the convolution kernel and the correction function were determined by a computer iteration program. The second method, which is a simplified version of the first, uses conventional backprojection, and takes a shorter computation time than the first method. The statistical noise of an image can be minimised by suitable positioning of the coordinate origin for the reconstruction. The theory of the two methods, their performance on statistical noise and some results of mathematical and experimental phantom studies are described.     

Brain single photon emission computed tomography findings in depressive pseudodementia patients

BACKGROUND: Recently, there have been studies suggesting that depressive pseudodementia would include early-stage dementing disorder. Through the comparison of the 99mTc-HMPAO single photon emission computed tomography (SPECT) image of depressive pseudodementia subjects, healthy comparison subjects, depressed subjects free of cognitive impairment, and dementia of Alzheimer's type (DAT) subjects, we aimed to see part of pathophysiology of the depressive pseudodementia of elderly patients. METHODS: Study subjects consisted of seven patients with depressive pseudodementia, seven healthy comparison subjects, seven patients with depression free of cognitive impairment, and eleven patients with DAT. Depression patients were diagnosed according to DSM-III-R. DAT patients were diagnosed by DSM III-R and NINCDS-ADRDA criteria of DAT. Other measures for assessment include Hamilton Rating Scale for Depression and Mini Mental State Exam. All underwent 99mTc-HMPAO SPECT scan. The images of each group were analyzed using statistical parametric mapping of Friston, which compares the images on voxel-by-voxel basis. RESULTS: The results were as follows (1) The DAT group showed significant decreases of cerebral blood flow (CBF) in the right frontal, right temporal region, and both parietal regions as compared with control group (P < 0.05). (2) The depression group showed a significant decrease of CBF in the left frontal region as compared with control group (P < 0.05). (3) The depressive pseudodementia group showed significant decreases of CBF in both parietal regions as compared with control group (P < 0.05). (4) The depressive pseudodementia group showed significant decreases of CBF in the right temporal region and both parietal regions as compared with depression group (P < 0.05). (5) The DAT group showed significant decreases of CBF in the right temporal region, both frontal regions, and both parietal regions as compared with depressive pseudodementia group (P < 0.05). LIMITATIONS: The small number of subjects may make it difficult to generalize from our results. Because decreased CBF in depressive pseudodementia is found while the subjects were depressed, we cannot tell whether it is a state marker or a trait marker. CONCLUSIONS: The depressive pseudodementia group showed decreased CBF in the temporo-parietal region, similar to that of the DAT group and different from that of the depression group.     

Quantitative pulmonary single photon emission computed tomography for radiotherapy applications.

Pulmonary imaging using single photon emission computed tomography (SPECT) is the focus of current radiotherapy research, including dose-response analysis and three-dimensional (3D) radiation treatment planning. Improvement in the quantitative capability of SPECT may help establish its potential role in this application as well as others requiring accurate knowledge of pulmonary blood flow. The purposes of this study were to quantitatively evaluate SPECT filtered backprojection (FBP) and ordered subset-expectation maximization (OS-EM) reconstruction implementations for measuring absolute activity concentration in lung phantom experiments, and to incorporate quantitative SPECT techniques in 3D-RTP for lung cancer. Quantitative FBP (nonuniform iterative Chang attenuation compensation, scatter correction, and 3D postreconstruction Metz filtering) and OS-EM implementations were compared with a "clinical" implementation of FBP (uniform multiplicative Chang attenuation compensation and post-reconstruction von Hann filtering), for their ability to improve quantification of inactive and active spherical defects in the lungs of an anthropomorphic torso phantom. Activity concentration estimates were found to depend on many factors, such as region of interest size, scatter subtraction constant (k), postreconstruction deconvolution filtering and, in the case of OS-EM, total number of iterations. In general, reconstruction implementations incorporating compensation for nonuniform attenuation and scatter provided reduced bias relative to the clinical implementation. Potential applications to lung radiotherapy, including dose-functional histograms and treatment planning are also discussed. SPECT has the potential to provide accurate estimates of lung activity distributions that, together with improved image quality, may be useful for the study and prediction of therapeutic response.     

NREM sleep EEG frequency spectral correlates of sleep complaints in primary insomnia subtypes

STUDY OBJECTIVES: To determine whether the frequency spectrum of the sleep EEG is a physiologic correlate of 1) the degree to which individuals with persistent primary insomnia (PPI) underestimate their sleep time compared with the traditionally scored polysomnogram (PSG) and 2) the sleep complaints in PPI subjects who have relatively long traditionally scored PSG sleep times and relatively greater underestimation of sleep time. DESIGN: We compared EEG frequency spectra from REM and NREM sleep in PPI subjects subtyped as subjective insomnia sufferers (those with relatively long total sleep time and relative underestimation of sleep time compared with PSG), and objective insomnia sufferers (those with relatively short PSG total sleep time) with EEG frequency spectra in normals. We also studied the correlation between these indices and the degree of underestimation of sleep. Further, we determined the degree to which sleep EEG indexes related to sleep complaints. SETTING: Duke University Medical Center Sleep Laboratory. PARTICIPANTS: Normal (N=20), subjective insomnia (N=12), and objective insomnia (N=18) subjects. INTERVENTIONS: N/A MEASUREMENTS AND RESULTS: Lower delta and greater alpha, sigma, and beta NREM EEG activity were found in the patients with subjective insomnia but not those with objective insomnia, compared with the normal subjects. These results were robust to changes in the subtyping criteria. No effects were found for REM spectral indexes. Less delta non- REM EEG activity predicted greater deviation between subjective and PSG estimates of sleep time across all subjects. For the subjective insomnia subjects, diminished low-frequency and elevated higher frequency non- REM EEG activity was associated with their sleep complaints. CONCLUSIONS: NREM EEG frequency spectral indexes appear to be physiologic correlates of sleep complaints in patients with subjective insomnia and may reflect heightened arousal during sleep.     

Monte Carlo evaluation of Compton scatter subtraction in single photon emission computed tomography

In single photon emission computed tomography (SPECT), Compton scattering produces a background that degrades the image quality and contributes erroneously to quantitative measurements. A clinically implementable compensation algorithm has previously been reported that subtracts a Compton scatter image, acquired in an energy window set below the energy of the photopeak, from the primary image acquired in the photopeak window. We present an evaluation and justification of the assumptions made in the previous empirical development of the subtraction algorithm. A Monte Carlo model of the SPECT system in which the Compton scattered vents may be followed independently of the nonscattered events was used to evaluate this subtraction technique. Simulation shows that the assumptions made in the experimental application of this algorithm were valid. Specifically (1) the "scatter" energy window used in the experiment (91-125 keV for imaging Tc-99m) contains only scattered events, (2) the shape of the line spread function (LSF) for the events in the scatter window is a reasonable approximation to the shape of the scatter in the photopeak window, and (3) the ratio of the number of scattered events in the photopeak window to the number of events in the scatter window is 0.57, close to the value of 0.5 derived heuristically. Thus, Monte Carlo simulation validates the basic assumptions underlying the empirical implementation of the scatter subtraction algorithm.     

Area weighted convolutional interpolation for data reprojection in single photon emission computed tomography

A data reprojection algorithm has been developed for use in single photon emission computed tomography on an array processor equipped computer system. The algorithm makes use of an accurate representation of pixel activity (uniform square pixel model of intensity distribution), and is rapidly performed due to the efficient handling of an array-based algorithm and the fast Fourier transform on parallel processing hardware. The algorithm consists of using a pixel driven nearest-neighbor projection operation to an array of subdivided projection bins. The subdivided project bin array is then convolved with the angle-dependent projection of the area of a uniform square pixel and compressed to original bin size. The new algorithm has thus been named the area weighted convolution (AWC) method of interpolation. When compared to nearest-neighbor and linear interpolation algorithms, the new AWC algorithm was found to be more accurate, having an accuracy approaching that of the line length algorithm. It also yielded an easier and more efficient implementation on parallel hardware than line length or linear interpolation, with faster execution times than either.     

A comparison of attenuation correction methods for quantitative single photon emission computed tomography

Transverse section tomograms of experimental phantoms and patients have been obtained using a GE 400T camera and a filtered back-projection reconstruction technique. These tomograms have been compared with the corresponding sections reconstructed from the same tomographic projection data, but using iterative algorithms with correction for photon attenuation. The comparison assesses the importance of including a correction for attenuation as well as demonstrating how closely a simple geometric attenuation correction, applied to the filtered back-projection reconstruction method, approximates to a more accurate correction incorporated in the computation of line integrals during iterative reconstruction. A comparison is also made between the behaviour of reconstruction algorithms with simulated projection data and real data in terms of convergence properties, and some shortcomings arising from simulation are noted.     

Effects of acute hypobaric hypoxia on regional cerebral blood flow distribution: a single photon emission computed tomography study in humans

Single Photon Emission Computed Tomography (SPECT) and radiopharmaceutical stabilizing agents allowed us to investigate regional cerebral blood flow (CBF) distribution in six resting healthy subjects during acute laboratory hypobaric hypoxic conditions. In the hypobaric experiment stabilized 99mTc-D, L-hexamethyl-propylene amine oxime was injected 40 min after reaching hypoxic conditions corresponding to an altitude of 5500 m above sea level. Arterial blood sample was taken after five additional minutes. Mean arterial oxygen pressure and haemoglobin saturation were 28 mmHg and 56%, respectively. The control experiment was performed similarly, apart from barometric pressure and blood gas analysis. We analysed CBF distribution in 12 regions of functional interest bilaterally in frontal, parietal, temporal, occipital cortex, in the hippocampus, in the basal ganglia and other central structures of brain. No overall effect of hypoxia on normalized regional CBF distribution in the considered regions was found. Motor cortex (Brodmann 4) and basal ganglia were the only regions in which hypobaric hypoxia significantly increased relative distribution of the radiopharmaceutical [F(1,5)=18.30; P < 0.008 and F(1,5)=10.85; P < 0.022, respectively]. Despite severe hypoxia, we did not observe any major regional CBF redistribution. We found a small relative increase in blood flow to the motor cortex and the basal ganglia, at rest after 40 min of hypobaric hypoxia, suggesting a preferential compensatory mechanism of these functional regions of brain.     

Increased EEG sigma and beta power during NREM sleep in primary insomnia

The hyperarousal model of primary insomnia suggests that a deficit of attenuating arousal during sleep might cause the experience of non-restorative sleep. In the current study, we examined EEG spectral power values for standard frequency bands as indices of cortical arousal and sleep protecting mechanisms during sleep in 25 patients with primary insomnia and 29 good sleeper controls. Patients with primary insomnia demonstrated significantly elevated spectral power values in the EEG beta and sigma frequency band during NREM stage 2 sleep. No differences were observed in other frequency bands or during REM sleep. Based on prior studies suggesting that EEG beta activity represents a marker of cortical arousal and EEG sleep spindle (sigma) activity is an index of sleep protective mechanisms, our findings may provide further evidence for the concept that a simultaneous activation of wake-promoting and sleep-protecting neural activity patterns contributes to the experience of non-restorative sleep in primary insomnia.     

Increased EEG sigma and beta power during NREM sleep in primary insomnia

The hyperarousal model of primary insomnia suggests that a deficit of attenuating arousal during sleep might cause the experience of non-restorative sleep. In the current study, we examined EEG spectral power values for standard frequency bands as indices of cortical arousal and sleep protecting mechanisms during sleep in 25 patients with primary insomnia and 29 good sleeper controls. Patients with primary insomnia demonstrated significantly elevated spectral power values in the EEG beta and sigma frequency band during NREM stage 2 sleep. No differences were observed in other frequency bands or during REM sleep. Based on prior studies suggesting that EEG beta activity represents a marker of cortical arousal and EEG sleep spindle (sigma) activity is an index of sleep protective mechanisms, our findings may provide further evidence for the concept that a simultaneous activation of wake-promoting and sleep-protecting neural activity patterns contributes to the experience of non-restorative sleep in primary insomnia.     

Hexa-methyl-propylene-amine-oxime (HMPAO) single photon emission computed tomography (SPECT) in epilepsy

Summary Twenty-eight normal volunteers and 74 seizure patients were investigated with hexa-methyl-propylene-amine-oxime (HMPAO) brain single photon emission computed tomography (SPECT). Fifty-four patients suffered from partial seizures and 20 patients had generalized seizures. Indices describing regional tracer distribution (RIs) were calculated in all investigated subjects. Regions whose RI exceeded the mean normal RI +/-3 SD were defined as abnormally perfused. In normals a significant interhemispheric asymmetry of HMPAO deposition was found, with higher RI values in the right frontal, parietal, temporal and occipital regions. In 59.3% of partial seizure patients abnormal RIs were found. Low RIs were detected predominantly in the frontal and temporal cortex, while elevated RIs were observed predominantly in the anterior basal ganglia. Only in 20% of the cases with generalized seizures, abnormal RIs were found. In one patient an ictal and 3 follow-up SPECT studies were obtained. Here SPECT results indicated transient rCBF changes between the ictal and seizure free state. EEG and SPECT foci were ipsilaterally located in 69.2% of the partial seizure cases. The results indicate that HMPAO brain SPECT is valuable for the detection of rCBF abnormalities in seizure patients and that patients with partial seizures have mostly several abnormally perfused areas in their brains.     

Iterative three-dimensional expectation maximization restoration of single photon emission computed tomography images: application in striatal imaging

Single photon emission computed tomography imaging suffers from poor spatial resolution and high statistical noise. Consequently, the contrast of small structures is reduced, the visual detection of defects is limited and precise quantification is difficult. To improve the contrast, it is possible to include the spatially variant point spread function of the detection system into the iterative reconstruction algorithm. This kind of method is well known to be effective, but time consuming. We have developed a faster method to account for the spatial resolution loss in three dimensions, based on a postreconstruction restoration method. The method uses two steps. First, a noncorrected iterative ordered subsets expectation maximization (OSEM) reconstruction is performed and, in the second step, a three-dimensional (3D) iterative maximum likelihood expectation maximization (ML-EM) a posteriori spatial restoration of the reconstructed volume is done. In this paper, we compare to the standard OSEM-3D method, in three studies (two in simulation and one from experimental data). In the two first studies, contrast, noise, and visual detection of defects are studied. In the third study, a quantitative analysis is performed from data obtained with an anthropomorphic striatal phantom filled with 123-I. From the simulations, we demonstrate that contrast as a function of noise and lesion detectability are very similar for both OSEM-3D and OSEM-R methods. In the experimental study, we obtained very similar values of activity-quantification ratios for different regions in the brain. The advantage of OSEM-R compared to OSEM-3D is a substantial gain of processing time. This gain depends on several factors. In a typical situation, for a 128 x 128 acquisition of 120 projections, OSEM-R is 13 or 25 times faster than OSEM-3D, depending on the calculation method used in the iterative restoration. In this paper, the OSEM-R method is tested with the approximation of depth independent resolution. For the striatum this approximation is appropriate, but for other clinical situations we will need to include a spatially varying response. Such a response is already included in OSEM-3D.     

Modifying constrained least-squares restoration for application to single photon emission computed tomography projection images

Image restoration methods have been shown to increase the contrast of nuclear medicine images by decreasing the effects of scatter and septal penetration. Image restoration can also reduce the high-frequency noise in the image. This study applies constrained least-squares (CLS) restoration to the projection images of single photon emission computed tomography (SPECT). In a previous study, it was noted that CLS restoration has the potential advantage of automatically adapting to the blurred object. This potential is confirmed using planar images. CLS restoration is then modified to improve its performance when applied to SPECT projection image sets. The modification was necessary because the Poisson noise in low count SPECT images causes considerable variation in the CLS filter. On phantom studies, count-dependent Metz restoration was slightly better than the modified CLS restoration method, according to measures of contrast and noise. However, CLS restoration was generally judged as yielding the best results when applied to clinical studies, apparently because of its ability to adapt to the image being restored.     

Psychological stress is associated with heightened physiological arousal during NREM sleep in primary insomnia

The objective of this study was to evaluate cross-sectional relationships among symptoms of psychological stress, sleep, and physiological arousal during non-rapid eye movement (NREM) sleep in a sample of 30 patients with chronic, primary insomnia (mean age, 30.2 years, 60% female). Study measures included indexes of subjective stress, visually scored sleep, and physiological arousal during NREM sleep: quantitative electroencephalogram (QEEG) and quantitative electrocardiogram (QEKG) measures. Psychological stress was more strongly related to indexes of physiological arousal during NREM sleep than to visually scored measures of sleep. Higher levels of perceived stress were associated with decreased EEG delta power (rho = -0.50, p < .01) and increased EEG beta power (rho = 0.38, p < .05). Increased frequency of stress-related avoidance behaviors was associated with decreased EKG high-frequency power (rho = -0.46, p < .05). Although QEEG measures were significantly correlated with sleep maintenance (QEEG delta power rho = 0.45, p < .01; QEEG beta power rho = -0.54, p < .01) and time spent in delta sleep (QEEG delta power rho = 0.65, p < .001; QEEG beta power rho = -0.65, p < .001), QEKG measures were unrelated to visually scored measures of sleep. Perceived stress and stress-related avoidance behaviors were associated with multiple indexes of physiological arousal during NREM sleep in patients with chronic, primary insomnia.     

Correlation between different memory systems and striatal dopamine D2/D3 receptor density: a single photon emission computed tomography study

OBJECTIVE: Although studies have indicated that the human basal ganglia have a specific role in different memory systems, the functional significance of the striatal dopamine activities for the basal ganglia remains less clear. This study assessed the relationship between measures of striatal dopamine activities and indices of different memory systems in healthy individuals. METHOD: Single photon emission computed tomography and [123I]IBZM (iodobenzamide) were used to assess the striatal dopamine D2/D3 receptor density in 62 healthy volunteers aged between 19 and 61 years. All subjects underwent a Wechsler Memory Scale-Revised test. RESULTS: Dopamine D2/D3 receptor densities in the striatum decline with age. Stepwise regression analysis showed that verbal delayed recall and working memory account for most of the variance in dopamine D2/D3 measurements. These relationships remain significantly after controlling for age effects. CONCLUSIONS: Brain striatal dopamine activities are also significantly associated with various memory systems, in addition to motor functions. This may explain why patients with neuropsychiatric disorders may have both memory and motor impairments.