Digital anthropomorphic phantoms of non-rigid human respiratory and voluntary body motion for investigating motion correction in emission imaging.

The development of methods for correcting patient motion in emission tomography has been receiving increased attention. Often the performance of these methods is evaluated through simulations using digital anthropomorphic phantoms, such as the commonly used extended cardiac torso (XCAT) phantom, which models both respiratory and cardiac motion based on human studies. However, non-rigid body motion, which is frequently seen in clinical studies, is not present in the standard XCAT phantom. In addition, respiratory motion in the standard phantom is limited to a single generic trend. In this work, to obtain a more realistic representation of motion, we developed a series of individual-specific XCAT phantoms, modeling non-rigid respiratory and non-rigid body motions derived from the magnetic resonance imaging (MRI) acquisitions of volunteers. Acquisitions were performed in the sagittal orientation using the Navigator methodology. Baseline (no motion) acquisitions at end-expiration were obtained at the beginning of each imaging session for each volunteer. For the body motion studies, MRI was again acquired only at end-expiration for five body motion poses (shoulder stretch, shoulder twist, lateral bend, side roll, and axial slide). For the respiratory motion studies, an MRI was acquired during free/regular breathing. The magnetic resonance slices were then retrospectively sorted into 14 amplitude-binned respiratory states, end-expiration, end-inspiration, six intermediary states during inspiration, and six during expiration using the recorded Navigator signal. XCAT phantoms were then generated based on these MRI data by interactive alignment of the organ contours of the XCAT with the MRI slices using a graphical user interface. Thus far we have created five body motion and five respiratory motion XCAT phantoms from the MRI acquisitions of six healthy volunteers (three males and three females). Non-rigid motion exhibited by the volunteers was reflected in both respiratory and body motion phantoms with a varying extent and character for each individual. In addition to these phantoms, we recorded the position of markers placed on the chest of the volunteers for the body motion studies, which could be used as external motion measurement. Using these phantoms and external motion data, investigators will be able to test their motion correction approaches for realistic motion obtained from different individuals. The non-uniform rational B-spline data and the parameter files for these phantoms are freely available for downloading and can be used with the XCAT license.

Tumor pretargeting in mice using MORF conjugated CC49 antibody and radiolabeled complimentary cMORF effector.

AIM: Using the antiCEA antibody MN14, a LS174T mouse tumor model has been successfully targeted with (99m)Tc for imaging and ¹88Re for radiotherapy by phosphorodiamidate morpholino oligomers (MORF)/complementary MORF (cMORF) pretargeting strategy. This investigation evaluated the antiTAG-72 antibody CC49 as an alternative to MN14 for this application. METHODS: Both CC49 and MN14 were labeled with ¹¹¹In via SCN-benzyl-DTPA and their biodistributions were compared to that of MN14 labeled via DTPA anhydride. Since the accessibility of the antibody to the effector is required for optimization of pretargeting, the internalization of both MORF-CC49 and MORF-MN14 antibodies in LS174T cells were evaluated in culture. In addition, the accessible concentration of MORF-CC49 antibody in tumor was determined in a series of pretargeting studies with escalating dosages of the [(99m)Tc]cMORF effector. Finally, using these results and our semi-empirical model, an imaging study was performed under optimal pretargeting conditions. RESULTS: The biodistribution of ¹¹¹In to trace the MN14 antibody depended significantly on the labeling method. Furthermore, both MORF-CC49 and MORF-MN14 antibodies showed rapid internalization in culture. Fortunately, the accessibility in tumor was found to be less seriously reduced in vivo. In a pretargeting study under optimal conditions, both by imaging and by necropsy, the [(99m)Tc]cMORF effector accumulated predominantly in the tumor of pretargeted mice. Normal tissue accumulations were minimal except in kidneys, liver, and a segment of intestines. CONCLUSION: MORF pretargeting with CC49 was equally successful in the LS174T tumor model to the MORF pretargeting with MN14. The MORF-CC49 antibody may therefore be considered for future investigations toward early clinical trials.

An Evaluation of Iterative Reconstruction Strategies on Mediastinal Lesion Detection Using Hybrid Ga-67 SPECT Images.

Hybrid LROC studies can be used to more realistically assess the impact of reconstruction strategies, compared to those constructed with digital phantoms. This is because hybrid data provides the background variability that is present in clinical imaging, as well as, control over critical imaging parameters, required to conduct meaningful tests. Hybrid data is obtained by adding Monte Carlo simulated lesions to disease free clinical projection data. Due to Ga-67 being a particularly challenging radionuclide for imaging, we use Ga-67 hybrid SPECT data to study the effectiveness of the various correction strategies developed to account for degradations in SPECT imaging. Our data was obtained using GE-VG dual detector SPECT-CT camera. After determining a target lesion contrast we conduct pilot LROC studies to obtain a near-optimal set of reconstruction parameters for the different strategies individually. These near-optimal parameters are then used to reconstruct the final evaluation study sets. All LROC study results reported here were obtained employing human observers only. We use final LROC study results to assess the impact of attenuation compensation, scatter compensation and detector resolution compensation on data reconstructed with the RBI-EM algorithm. We also compare these with FBP reconstructions of the same dataset. Our experiment indicates an improvement in detection accuracy, as various degradations inherent in the image acquisition process are compensated for in the reconstruction process.

Estimation of the Rigid-Body Motion from Three-Dimensional Images Using a Generalized Center-of-Mass Points Approach.

We present an analytical method for the estimation of rigid-body motion in sets of three-dimensional SPECT and PET slices. This method utilizes mathematically defined generalized center-of-mass points in images, requiring no segmentation. It can be applied to compensation of the rigid-body motion in both SPECT and PET, once a series of 3D tomographic images are available. We generalized the formula for the center-of-mass to obtain a family of points co-moving with the object's rigid-body motion. From the family of possible points we chose the best three points which resulted in the minimum root-mean-square difference between images as the generalized center-of-mass points for use in estimating motion. The estimated motion was used to sum the sets of tomographic images, or incorporated in the iterative reconstruction to correct for motion during reconstruction of the combined projection data. For comparison, the principle-axes method was also applied to estimate the rigid-body motion from the same tomographic images. To evaluate our method for different noise levels, we performed simulations with the MCAT phantom. We observed that though noise degraded the motion-detection accuracy, our method helped in reducing the motion artifact both visually and quantitatively. We also acquired four sets of the emission and transmission data of the Data Spectrum Anthropomorphic Phantom positioned at four different locations and/or orientations. From these we generated a composite acquisition simulating periodic phantom movements during acquisition. The simulated motion was calculated from the generalized center-of-mass points calculated from the tomographic images reconstructed from individual acquisitions. We determined that motion-compensation greatly reduced the motion artifact. Finally, in a simulation with the gated MCAT phantom, an exaggerated rigid-body motion was applied to the end-systolic frame. The motion was estimated from the end-diastolic and end-systolic images, and used to sum them into a summed image without obvious artifact. Compared to the principle-axes method, in two of the three comparisons with anthropomorphic phantom data our method estimated the motion in closer agreement to than of the Polaris system than the principal-axes method, while the principle-axes method gave a more accurate estimation of motion in most cases for the MCAT simulations. As an image-driven approach, our method assumes angularly complete data sets for each state of motion. We expect this method to be applied in correction of respiratory motion in respiratory gated SPECT, and respiratory or other rigid-body motion in PET.

Factors Influencing Lesion Detection in SPECT Lung Images.

An earlier localization ROC (LROC) study that found attenuation correction (AC) degraded the detection of solitary pulmonary nodules (SPN) in hybrid SPECT lung images had several potential shortcomings related to the simulation methods. We sought to address these issues with a revised LROC study. Clinical Tc-99m NeoTect scans acquired with a simultaneous transmission-emission protocol defined the normal cases in a single-slice LROC study. Abnormal cases contained a simulated 1-cm lung lesion. Four rescaled-block-iterative EM (RBI) reconstruction strategies applied: 1) AC, scatter correction (SC), and resolution compensation (RC); 2) AC only; 3) RC only; and 4) no corrections (NC). Images from these strategies underwent 3D Gaussian post-smoothing. Performances were defined by the average area under the LROC curve obtained from three human observers. The strategy ranking in order of decreasing performance was: 1) RBI with RC; 2) RBI with all corrections; 3) RBI with AC; and 4) RBI with no corrections. A multireader-multicase (MRMC) analysis only found significant patient and patient-strategy effects. The conflicting results concerning AC from this study and the previous one may revolve around lesion masking effects, which, by design, were not a factor in the current study.

The influence of attenuation and scatter compensation on the apparent distribution of Tc-99m sestamibi in cardiac slices.

BACKGROUND: Our objective was to study the differences in relative count distributions in the left ventricular walls with attenuation compensation (AC) versus AC and triple-energy-window scatter compensation (SC), compared with standard filtered backprojection (FBP). METHODS AND RESULTS: Two hundred patients identified as having normal cardiac perfusion with FBP after undergoing either pharmacologically or physiologically induced stress were included in this study. Projection data were reconstructed with FBP, 10 iterations of ordered-subset expectation-maximization (OSEM) with AC, and OSEM with AC+SC. A comparison was made of average percentage of maximum counts within each of 9 regions of CEqual (Marconi Medical Systems, Inc, Cleveland, Ohio) polar maps (ie, the apex, 4 midventricular regions, and 4 basal regions). Compared with OSEM(AC), a slight decrease at the apex exists when SC is included. The elevated inferior-to-anterior count ratio in the midventricular and basal regions noted with OSEM(AC) decreased to close to 1.0 with OSEM(AC+SC). The anterior-to-lateral ratio for both regions was closest to 1.0 for OSEM(AC+SC). In the midventricular region, the lateral-to-septal ratio decreased further below 1.0 with OSEM(AC+SC) than it did with OSEM(AC). This was the only basal ratio not to improve to close to 1.0 with OSEM(AC+SC). In a subset of patients identified at the time of clinical reading as having a possible attenuation-caused decrease in the inferior region, AC elevated the inferior-to-anterior ratio to above 1.0 for the midventricular region. AC+SC resulted in a ratio of near 1.0 for this region. In another subset of patients identified as having anterior attenuation artifacts, compensation methods (either AC or AC+SC) failed to show an improvement compared with FBP. CONCLUSIONS: AC and SC improve the uniformity of the polar map, especially by bringing the inferior-to-anterior ratio closer to 1.0. Further investigation is necessary to determine the cause of the increased midventricular septal polar map count. In addition, the subset of patients identified as having breast-like attenuation artifacts causing a decreased polar map count in the anterior wall (relative to the inferior wall) also needs further attention.

Improved image quality and computation reduction in 4-D reconstruction of cardiac-gated SPECT images.

Spatiotemporal reconstruction of cardiac-gated SPECT images permits us to obtain valuable information related to cardiac function. However, the task of reconstructing this four-dimensional (4-D) data set is computation intensive. Typically, these studies are reconstructed frame-by-frame: a nonoptimal approach because temporal correlations in the signal are not accounted for. In this work, we show that the compression and signal decorrelation properties of the Karhunen-Loève (KL) transform may be used to greatly simplify the spatiotemporal reconstruction problem. The gated projections are first KL transformed in the temporal direction. This results in a sequence of KL-transformed projection images for which the signal components are uncorrelated along the time axis. As a result, the 4-D reconstruction task is simplified to a series of three-dimensional (3-D) reconstructions in the KL domain. The reconstructed KL components are subsequently inverse KL transformed to obtain the entire spatiotemporal reconstruction set. Our simulation and clinical results indicate that KL processing provides image sequences that are less noisy than are conventional frame-by-frame reconstructions. Additionally, by discarding high-order KL components that are dominated by noise, we can achieve savings in computation time because fewer reconstructions are needed in comparison to conventional frame-by-frame reconstructions.

LROC analysis of detector-response compensation in SPECT.

Localization ROC (LROC) observer studies examined whether detector response compensation (DRC) in ordered-subset, expectation-maximization (OSEM) reconstructions helps in the detection and localization of hot tumors. Simulated gallium (Ga-67) images of the thoracic region were used in the study. The projection data modeled the acquisition of attenuated 93- and 185-keV photons with a medium-energy parallel-hole collimator, but scatter was not modeled. Images were reconstructed with five strategies: 1) OSEM with no DRC; 2) OSEM preceded by restoration filtering; 3) OSEM with iterative DRC; 4) OSEM with an ideal DRC; and 5) filtered backprojection (FBP) with no DRC. All strategies included attenuation correction. There were four LROC studies conducted. In a study using a single tumor activity, the ideal DRC offered the best performance, followed by iterative DRC, restoration filtering, OSEM with no DRC, and FBP. Statistical significance at the 5% level was found between all pairs of strategies except for restoration filtering and OSEM with no DRC. A similar ranking was found for a more realistic study using multiple tumor activities. Additional studies considered the effects of OSEM iteration number and tumor activity on the detection improvement that iterative DRC offered with respect to OSEM with no DRC.

Comparing filtered backprojection and ordered-subsets expectation maximization for small-lesion detection and localization in 67Ga SPECT.

UNLABELLED: Iterative reconstruction of SPECT images has recently become clinically available as an alternative to filtered backprojection (FBP). However, there is conflicting evidence on whether iterative reconstruction, such as with the ordered-subsets expectation maximization (OSEM) algorithm, improves diagnostic performance over FBP. The study objective was to determine if the detection and localization of small lesions in simulated thoracic gallium SPECT images are better with OSEM reconstruction than with FBP, both with and without attenuation correction (AC). METHODS: Images were simulated using an analytic projector acting on the mathematic cardiac torso computer phantom. Perfect scatter rejection was assumed. Lesion detection accuracy was assessed using localization receiver operating characteristic methodology. The images were read by 5 nuclear medicine physicians. For each reconstruction strategy and for each observer, data were collected in 2 viewing sessions of 100 images. Two-way ANOVA and, when indicated, the Scheffé multiple comparisons test were applied to check for significant differences. RESULTS: Little difference in the accuracy of detection or localization was seen between FBP with and without AC. OSEM with AC extended the contrast range for accurate lesion detection and localization over that of the other methods investigated. Without AC, no significant difference between OSEM and FBP reconstruction was detected. CONCLUSION: OSEM with AC may improve the detection and localization of thoracic gallium-labeled lesions over FBP reconstruction.

A mathematical model of motion of the heart for use in generating source and attenuation maps for simulating emission imaging.

This manuscript documents the alteration of the heart model of the three-dimensional (3D) mathematical cardiac torso (MCAT) phantom to represent cardiac motion. The objective of the inclusion of motion was to develop a digital simulation of the heart such that the impact of cardiac motion on single-photon emission computed tomography (SPECT) imaging could be assessed and methods of quantitating cardiac function could be investigated. The motion of the gated 3D MCAT's (gMCAT) heart is modeled using 128 separate and evenly spaced time samples from a blood volume curve approximating an average heart cycle. Sets of adjacent time samples can be grouped together to represent a single time interval within the heart cycle. Maximum and minimum chamber volumes were selected to be similar to those of a normal healthy person while the total heart volume stayed constant during the cardiac cycle. Myocardial mass was conserved during the cardiac cycle and the bases of the ventricles were modeled as moving towards the static apex. The orientation of the 3D MCAT heart was changed during contraction to rotate back and forth around the long axis through the center of the left ventricle (LV) using the end systolic time interval as the time point at which to reverse direction. Simple respiratory motion was also introduced by changing the orientation of the long axis of the heart to represent its variation with respiration. Heart models for 24 such orientations spanning the range of motion during the respiratory cycle were averaged together for each time sample to represent the blurring of the heart during the acquisition of multiple cardiac cycles. Finally, an option to model apical thinning of the myocardium was included. As an illustration of the application of the gMCAT phantom, the gated heart model was evaluated by measuring myocardial wall thickening. A linear relationship was obtained between maximum myocardial counts and myocardial thickness, similar to published results. Similar results were obtained for full width at half maximum (FWHM) measurements. With the presence of apical thinning, an apparent increase in counts in the apical region compared to the other heart walls in the absence of attenuation compensation turns into an apparent decrease in counts with attenuation compensation. The apical decrease was more prominent in end systole (ES) than end diastole (ED) due to the change in the partial volume effect. These observations agree with clinical trends. It is concluded that the gMCAT phantom can be used to study the influence of various physical parameters on radionuclide perfusion imaging.

An investigation of the estimation of ejection fractions and cardiac volumes by a quantitative gated SPECT software package in simulated gated SPECT images.

BACKGROUND: The purpose of this investigation was to determine the accuracy of the estimation of ejection fractions (EFs) and left ventricular volumes from a commercially available software package (Quantitative Gated SPECT [QGS]) as a function of different true EFs, count level in the acquisitions, severity and location of perfusion defects, increasing hepatic activity, and modified wall motion. METHODS AND RESULTS: The dynamic mathematic cardiac-torso digital phantom was used to create three-dimensional source and attenuation maps representing the distribution of a technetium-99m-labeled cardiac perfusion agent in the chest. Three hearts with varying end-systolic volumes were used to investigate different EFs. Perfusion defects were created as localized uptake within selected portions of the cardiac walls, scaled to the desired fraction of the normal wall uptake, and subtracted from the normal distribution. The hepatic uptake was increased up to five times of the normal heart uptake to investigate the influence of a "hot" liver. Alteration of lateral wall motion was also investigated. A three-dimensional projector that included the influence of distance-dependent spatial resolution and nonuniform attenuation was then used to create projection images. The projections were scaled to the desired acquisition count level, and Poisson noise was added. Automatic determination of EF slightly overestimated the true EF for normal count levels by 3% to 7% of the true EF and underestimated the true EF by up to 9% for very low count levels for 180-degree reconstructions. The accuracy for determining the volumes was not as high as for the EFs (an average error of 12% was observed). The calculated EFs were relatively accurate for perfusion defects of 50% or less. When perfusion defects exceeded 50%, extracardiac counts were included in the heart contours, causing larger underestimations of EF. With removal of the extracardiac counts, the EFs increased. With a hepatic uptake of two or more times the heart uptake, no meaningful EF could be obtained. Either drawing a single region of interest for every slice or use of the manual mode with constrain option could remarkably improve the estimation. The accuracy of the calculation of EF and volumes for the heart with stationary wall was fairly high but decreased significantly when coupled with perfusion defects. CONCLUSION: It is concluded that the QGS program evaluates the functional parameter of EF accurately. The biggest limitations occurred in determining the appropriate cardiac contour if areas with very high extracardiac counts were present in the heart slices, and when a greater than 50% decrease occurred in uptake for perfusion defects.

Reducing the influence of the partial volume effect on SPECT activity quantitation with 3D modelling of spatial resolution in iterative reconstruction.

Quantitative parameters such as the maximum and total counts in a volume are influenced by the partial volume effect. The magnitude of this effect varies with the non-stationary and anisotropic spatial resolution in SPECT slices. The objective of this investigation was to determine whether iterative reconstruction which includes modelling of the three-dimensional (3D) spatial resolution of SPECT imaging can reduce the impact of the partial volume effect on the quantitation of activity compared with filtered backprojection (FBP) techniques which include low-pass, and linear restoration filtering using the frequency distance relationship (FDR). The iterative reconstruction algorithms investigated were maximum-likelihood expectation-maximization (MLEM), MLEM with ordered subset acceleration (ML-OS), and MLEM with acceleration by the rescaled-block-iterative technique (ML-RBI). The SIMIND Monte Carlo code was used to simulate small hot spherical objects in an elliptical cylinder with and without uniform background activity as imaged by a low-energy ultra-high-resolution (LEUHR) collimator. Centre count ratios (CCRs) and total count ratios (TCRs) were determined as the observed counts over true counts. CCRs were unstable while TCRs had a bias of approximately 10% for all iterative techniques. The variance in the TCRs for ML-OS and ML-RBI was clearly elevated over that of MLEM, with ML-RBI having the smaller elevation. TCRs obtained with FDR-Wiener filtering had a larger bias (approximately 30%) than any of the iterative reconstruction methods but near stationarity is also reached. Butterworth filtered results varied by 9.7% from the centre to the edge. The addition of background has an influence on the convergence rate and noise properties of iterative techniques.

Evaluation of right and left ventricular volume and ejection fraction using a mathematical cardiac torso phantom.

UNLABELLED: The availability of gated SPECT has increased the interest in the determination of volume and ejection fraction of the left ventricle (LV) for clinical diagnosis. However, the same indices for the right ventricle (RV) have been neglected. The objective of this investigation was to use a mathematical model of the anatomical distribution of activity in gated blood-pool imaging to evaluate the accuracy of two ventricular volume and ejection fraction determination methods. In this investigation, measurements from the RV were emphasized. METHODS: The mathematical cardiac torso phantom, developed to study LV myocardium perfusion, was modified to simulate the radioactivity distribution of a 99mTc-gated blood-pool study. Twenty mathematical cardiac torso phantom models of the normal heart with different LV volumes (122.3 +/- 11.0 ml), RV volumes (174.6 +/- 22.3 ml) and stroke volumes (75.7 +/- 3.3 ml) were randomly generated to simulate variations among patients. An analytical three-dimensional projector with attenuation and system response was used to generate SPECT projection sets, after which noise was added. The projections were simulated for 128 equidistant views in a 360 degrees rotation mode. RESULTS: The radius of rotation was varied between 24 and 28 cm to mimic such variation in patient acquisitions. The 180 degrees and 360 degrees projection sets were reconstructed using the filtered backprojection reconstruction algorithm with Butter-worth filtering. Comparison was made with and without application of the iterative Chang attenuation correction algorithm. Volumes were calculated using a modified threshold and edge detection method (hybrid threshold), as well as a count-based method. A simple background correction procedure was used with both methods. CONCLUSION: Results indicate that cardiac functional parameters can be measured with reasonable accuracy using both methods. However, the count-based method had a larger bias than the hybrid threshold method when RV parameters were determined for 180 degrees reconstruction without attenuation correction. This bias improved after attenuation correction. The count-based method also tended to overestimate the end systolic volume slightly. An improved background correction could possibly alleviate this bias.

A Monte Carlo evaluation of the channel ratio scatter correction method.

Several methods exist to eliminate the contribution of scattered photons during imaging. One of these, the channel ratio (CR) scatter correction technique, uses the change in the ratio of counts from two symmetrical adjacent energy windows straddling the energy photopeak. The accuracy of the results depends upon the assumption that the ratio of the scatter components in the two windows (H value) is constant and is independent of the relative size of the scatter contribution. In this study a Monte Carlo simulation was used to investigate the behaviour of the scatter component for different source sizes at different depths. Four disc sources containing a 99Tcm solution were simulated at different depths as imaged with a scintillation camera. Two 10% energy windows with 5% offsets to either side of the 140 keV photopeak of 99Tcm were used. The ratio of the scattered counts in the lower energy window to the scattered counts in the upper window (true H value) was determined from the simulation for each source at every depth. Since it is impossible to measure the true H value at different organ depths during a clinical study, the use of an average H value was considered. Scatter correction was applied to the images simulated at the various depths in water. The geometric mean was calculated and attenuation correction performed assuming mono-exponential attenuation. For quantitation purposes the corrected counts were expressed in terms of a references source. The choice of the reference source yielding the best quantitative results was also investigated. Results of this Monte Carlo simulation study show that although the true H value depends on both source size and depth of the source in the scattering medium, the CR scatter correction technique can be applied successfully when an average H value is used.

The effect of exercise on normal splenic volume measured with SPECT.

In a study group of 20 healthy young men, splenic volume was determined with SPECT before and after exercise. A randomly chosen control group of 10 comparable men was studied similarly, but without exercise intervention. The mean splenic volume did not change significantly in the control group (i.e., from 292.9 ml to 282.1 ml [P = 0.75]). The mean splenic volume decreased 60.1 ml from 279.4 to 219.3 ml (21.5%) in the study group and this was highly significant (P = 0.01). Although exercise induced splenic autotransfusion is generally considered to be unimportant in humans, significant splenic contractility was observed with this technique. In the normal individual with a large spleen or with functional splenomegaly, the contractility response may become more important and can now be measured in a quick and easy manner.

Determination of glomerular filtration rate with radionuclide renography and direct urinary activity quantitation.

OBJECTIVE: The direct urinary activity quantitation method is quick (approximately 40 minutes), requires only a single blood sample, is performed as part of standard renal scanning and shows high accuracy compared with 24-hour creatinine clearance. The purpose was to evaluate the practical application and accuracy of this technique at our clinic. DESIGN: Direct urinary activity quantitation was done in patients scheduled for routine radionuclide renography and compared to standard multiple-blood-sample techniques by means of Cr-51-EDTA and Tc-99m-DTPA. SETTING: Academic Medical Complex, Department of Nuclear Medicine, Universitas Hospital, Bloemfontein. PARTICIPANTS: Fifteen patients scheduled for routine radionuclide renography (glomerular filtration rate (GFR) determination) were voluntarily enrolled in the study. The GFRs of selected patients varied over a wide range. Possible obstructive uropathy was excluded. MAIN OUTCOME MEASURES: GFRs obtained by the direct urinary method were compared with GFRs determined by multisample Cr-51-EDTA and Tc-99m DTPA. RESULTS: GFRs from the direct urinary method compared with multisample Tc-99m-DTPA showed differences from -19.85 to 22.95 ml/min with a mean of 0.2 (+/- 12.25) ml/min (r = 0.93). When compared with multisample Cr-51 EDTA, differences ranged from -34.35 to 21.05 ml/min with a mean of -4.25 (+/- 16.08) ml/min (r = 0.90). CONCLUSION: The direct urinary activity technique is CONCLUSION: The direct urinary activity technique is easily applied and highly accurate compared with previous standardised multisample radionuclide techniques for determination of GFR.

Difference between splenic volume measured at necropsy and that measured in vivo by radionuclide tomography.

Reference values for splenic volume used in this study are based on necropsy measurements made in 1970. Because the volumes (measured by radionuclide tomography) seemed to be consistently greater than the necropsy values, the splenic volume in 35 healthy male volunteers ranging in age from 18-30 years (median age, 21 years) was studied. Their mean (SD) splenic volume was 281 (73) ml compared with 138 (34) ml in the reference group. The mean splenic volume in healthy volunteers seems to be significantly greater than in cadavers, and these results suggest that reference range for splenic volume in vivo should be revised.

Verification of a varying threshold edge detection SPECT technique for spleen volume: a comparison with computed tomography volumes.

SPECT enables quantitation of organ volume with radionuclide techniques using threshold edge detection methods. Previous phantom studies showed that a negative correlation exists between volume and threshold value. In those studies, the use of calibration curves were believed to correct for volume dependence on threshold values. The aim of this study was to evaluate the accuracy of spleen volume determination in 20 patients with SPECT by employing a varying threshold edge detection technique with volumes derived from CT. All patients had both radionuclide and CT examinations that were reconstructed with a filtered backprojection algorithm. During SPECT reconstruction, transverse slices were obtained with attenuation correction (Method A) and without attenuation correction (Method B). CT volumes were calculated from manually drawn regions of interest, whereas SPECT volumes were calculated with an automated algorithm using previously determined calibration curves. A confidence interval for calculated SPECT volumes also was calculated because of possible errors in the threshold value. The spleen volumes studied ranged from 91.2 ml to 1660.1 ml. Regression analysis yielded equations of CT = 0.97 SPECT + 7.07 (r = 0.996) and CT = 1.05 SPECT - 19.25 (r = 0.990) between CT and SPECT spleen volumes with a standard error of the y estimates of 31.10 ml and 54.47 ml, respectively. A mean percentage difference of 10.5% +/- 7.6% and 11.4% +/- 6.6% in spleen volume was obtained for Methods A and B in comparison with CT spleen volumes. The threshold value varied between 40.9% and 32.4% for Method A and between 41.2% and 28.5% for Method B because the spleen volume is increased. The varying threshold edge detection technique described in this paper can be implemented successfully in the clinical setting.

The channel ratio method of scatter correction for radionuclide image quantitation.

The accuracy of quantitation of radionuclide distributions in human tissue with the scintillation camera is decreased by attenuation and scatter of photons. If scatter correction is applied satisfactorily, narrow beam attenuation can be applied. In this article, a scatter correction technique, the channel ratio (CR) method, is introduced. The CR scatter correction method is proposed for quantitation of the radionuclide distribution in organs. The improvement in the geometrical resolution was measured and examples of clinical images are presented. In this method, the change in the ratio of counts from two symmetrical adjacent energy windows straddling the energy photopeak was used to eliminate the contribution of scattered photons during imaging with 99mTc. The theory and methods for the empirical affirmation are described. To apply the CR scatter correction method, two constants, the ratio of primary photons G and the ratio of scattered photons H in the same windows, were determined. Different sized sources in varying depths of water were imaged. When the source activities were quantified after scatter correction with the CR method, the measurements ranged from 96%-108% in comparison to the reference value in 100 mm water. The scatter fraction increased from 0.20 in 10 mm water to 1.44 in 200 mm water. The geometrical resolution expressed as full width at tenth maximum in 150 mm water improved by 30.4% and was restored to the value of the geometrical resolution in air. The CR scatter correction method is a simple method to correct for scatter in order to facilitate accurate quantitation of the radionuclide distribution during imaging with a scintillation camera.

The effects of different correction techniques on absolute volume determination with SPECT using a threshold edge detection method.

Quantitation of planar radionuclide images is hampered by structures containing radioactivity which overlie or underlie the organ of interest. The introduction of single photon emission computerized tomography (SPECT) overcame this problem to a large extent and enhanced the contrast of the images. Attenuation of photons, however, degrades the resultant SPECT images and correction methods for photon absorption and scatter were subsequently proposed. The different correction methods have variable effects on the reconstructed images. If threshold techniques are used to quantitate organ volume, i.e., combining pixels with the same percentage of the maximum pixel count in the volume, the selected threshold values which give the most accurate volume determination, will be affected by the specific correction method used. In this study, the effect of various SPECT image correction methods on threshold was investigated. A thorax phantom containing volumes ranging from 30 to 1200 ml was used. Threshold values varying from 45.6% (210 ml without any correction) to 23.7% (1200 ml with a combination of scatter subtraction and attenuation correction) were used to produce correct quantitation when different methods were investigated. A negative correlation was found between threshold and volume. This reduction in threshold was most prominent when scatter and attenuation correction were combined. This study shows that correction methods for attenuation of photons influence the threshold value for volume quantitation and the use of a constant threshold value could lead to underestimation of larger volumes.