A Monte Carlo and physical phantom evaluation of quantitative In-111 SPECT

Three-dimensional radiobiological dosimetry of kidneys for treatment planning in peptide receptor radionuclide therapy

PURPOSE

Peptide receptor radionuclide therapy (PRRT) delivers high absorbed doses to kidneys and may lead to permanent nephropathy. Reliable dosimetry of kidneys is thus critical for safe and effective PRRT. The aim of this work was to assess the feasibility of planning PRRT based on 3D radiobiological dosimetry (3D-RD) in order to optimize both the amount of activity to administer and the fractionation scheme, while limiting the absorbed dose and the biological effective dose (BED) to the renal cortex.

METHODS

Planar and SPECT data were available for a patient examined with (111)In-DTPA-octreotide at 0.5 (planar only), 4, 24, and 48 h post-injection. Absorbed dose and BED distributions were calculated for common therapeutic radionuclides, i.e., (111)In, (90)Y and (177)Lu, using the 3D-RD methodology. Dose-volume histograms were computed and mean absorbed doses to kidneys, renal cortices, and medullae were compared with results obtained using the MIRD schema (S-values) with the multiregion kidney dosimetry model. Two different treatment planning approaches based on (1) the fixed absorbed dose to the cortex and (2) the fixed BED to the cortex were then considered to optimize the activity to administer by varying the number of fractions.

RESULTS

Mean absorbed doses calculated with 3D-RD were in good agreement with those obtained with S-value-based SPECT dosimetry for (90)Y and (177)Lu. Nevertheless, for (111)In, differences of 14% and 22% were found for the whole kidneys and the cortex, respectively. Moreover, the authors found that planar-based dosimetry systematically underestimates the absorbed dose in comparison with SPECT-based methods, up to 32%. Regarding the 3D-RD-based treatment planning using a fixed BED constraint to the renal cortex, the optimal number of fractions was found to be 3 or 4, depending on the radionuclide administered and the value of the fixed BED. Cumulative activities obtained using the proposed simulated treatment planning are compatible with real activities administered to patients in PRRT.

CONCLUSIONS

The 3D-RD treatment planning approach based on the fixed BED was found to be the method of choice for clinical implementation in PRRT by providing realistic activity to administer and number of cycles. While dividing the activity in several cycles is important to reduce renal toxicity, the clinical outcome of fractionated PRRT should be investigated in the future.

3D calculation of absorbed dose for 131I-targeted radiotherapy: a Monte Carlo study

Various methods, such as those developed by the Medical Internal Radiation Dosimetry (MIRD) Committee of the Society of Nuclear Medicine or employing dose point kernels, have been applied to the radiation dosimetry of (131)I radionuclide therapy. However, studies have not shown a strong relationship between tumour absorbed dose and its overall therapeutic response, probably due in part to inaccuracies in activity and dose estimation. In the current study, the GATE Monte Carlo computer code was used to facilitate voxel-level radiation dosimetry for organ activities measured in an (131)I-treated thyroid cancer patient. This approach allows incorporation of the size, shape and composition of organs (in the current study, in the Zubal anthropomorphic phantom) and intra-organ and intra-tumour inhomogeneities in the activity distributions. The total activities of the tumours and their heterogeneous distributions were measured from the SPECT images to calculate the dose maps. For investigating the effect of activity distribution on dose distribution, a hypothetical homogeneous distribution of the same total activity was considered in the tumours. It was observed that the tumour mean absorbed dose rates per unit cumulated activity were 0.65E-5 and 0.61E-5 mGY MBq(-1) s(-1) for the uniform and non-uniform distributions in the tumour, respectively, which do not differ considerably. However, the dose-volume histograms (DVH) show that the tumour non-uniform activity distribution decreases the absorbed dose to portions of the tumour volume. In such a case, it can be misleading to quote the mean or maximum absorbed dose, because overall response is likely limited by the tumour volume that receives low (i.e. non-cytocidal) doses. Three-dimensional radiation dosimetry, and calculation of tumour DVHs, may lead to the derivation of clinically reliable dose-response relationships and therefore may ultimately improve treatment planning as well as response assessment for radionuclide therapy.

MIRD pamphlet No. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy

In internal radionuclide therapy, a growing interest in voxel-level estimates of tissue-absorbed dose has been driven by the desire to report radiobiologic quantities that account for the biologic consequences of both spatial and temporal nonuniformities in these dose estimates. This report presents an overview of 3-dimensional SPECT methods and requirements for internal dosimetry at both regional and voxel levels. Combined SPECT/CT image-based methods are emphasized, because the CT-derived anatomic information allows one to address multiple technical factors that affect SPECT quantification while facilitating the patient-specific voxel-level dosimetry calculation itself. SPECT imaging and reconstruction techniques for quantification in radionuclide therapy are not necessarily the same as those designed to optimize diagnostic imaging quality. The current overview is intended as an introduction to an upcoming series of MIRD pamphlets with detailed radionuclide-specific recommendations intended to provide best-practice SPECT quantification-based guidance for radionuclide dosimetry.

A method for energy window optimization for quantitative tasks that includes the effects of model-mismatch on bias: application to Y-90 bremsstrahlung SPECT imaging

Quantitative Yttrium-90 ((90)Y) bremsstrahlung single photon emission computed tomography (SPECT) imaging has shown great potential to provide reliable estimates of (90)Y activity distribution for targeted radionuclide therapy dosimetry applications. One factor that potentially affects the reliability of the activity estimates is the choice of the acquisition energy window. In contrast to imaging conventional gamma photon emitters where the acquisition energy windows are usually placed around photopeaks, there has been great variation in the choice of the acquisition energy window for (90)Y imaging due to the continuous and broad energy distribution of the bremsstrahlung photons. In quantitative imaging of conventional gamma photon emitters, previous methods for optimizing the acquisition energy window assumed unbiased estimators and used the variance in the estimates as a figure of merit (FOM). However, for situations, such as (90)Y imaging, where there are errors in the modeling of the image formation process used in the reconstruction there will be bias in the activity estimates. In (90)Y bremsstrahlung imaging this will be especially important due to the high levels of scatter, multiple scatter, and collimator septal penetration and scatter. Thus variance will not be a complete measure of reliability of the estimates and thus is not a complete FOM. To address this, we first aimed to develop a new method to optimize the energy window that accounts for both the bias due to model-mismatch and the variance of the activity estimates. We applied this method to optimize the acquisition energy window for quantitative (90)Y bremsstrahlung SPECT imaging in microsphere brachytherapy. Since absorbed dose is defined as the absorbed energy from the radiation per unit mass of tissues in this new method we proposed a mass-weighted root mean squared error of the volume of interest (VOI) activity estimates as the FOM. To calculate this FOM, two analytical expressions were derived for calculating the bias due to model-mismatch and the variance of the VOI activity estimates, respectively. To obtain the optimal acquisition energy window for general situations of interest in clinical (90)Y microsphere imaging, we generated phantoms with multiple tumors of various sizes and various tumor-to-normal activity concentration ratios using a digital phantom that realistically simulates human anatomy, simulated (90)Y microsphere imaging with a clinical SPECT system and typical imaging parameters using a previously validated Monte Carlo simulation code, and used a previously proposed method for modeling the image degrading effects in quantitative SPECT reconstruction. The obtained optimal acquisition energy window was 100-160 keV. The values of the proposed FOM were much larger than the FOM taking into account only the variance of the activity estimates, thus demonstrating in our experiment that the bias of the activity estimates due to model-mismatch was a more important factor than the variance in terms of limiting the reliability of activity estimates.

Accuracy and precision of radioactivity quantification in nuclear medicine images

The ability to reliably quantify activity in nuclear medicine has a number of increasingly important applications. Dosimetry for targeted therapy treatment planning or for approval of new imaging agents requires accurate estimation of the activity in organs, tumors, or voxels at several imaging time points. Another important application is the use of quantitative metrics derived from images, such as the standard uptake value commonly used in positron emission tomography (PET), to diagnose and follow treatment of tumors. These measures require quantification of organ or tumor activities in nuclear medicine images. However, there are a number of physical, patient, and technical factors that limit the quantitative reliability of nuclear medicine images. There have been a large number of improvements in instrumentation, including the development of hybrid single-photon emission computed tomography/computed tomography and PET/computed tomography systems, and reconstruction methods, including the use of statistical iterative reconstruction methods, which have substantially improved the ability to obtain reliable quantitative information from planar, single-photon emission computed tomography, and PET images.

Half-time myocardial perfusion SPECT imaging with attenuation and Monte Carlo-based scatter correction

PURPOSE

To test the potential of a new reconstruction algorithm with Monte Carlo-based scatter correction in half-time myocardial perfusion single-photon emission computed tomography (SPECT).

MATERIALS AND METHODS

The mathematical four-dimensional NURBS-based Cardiac-Torso phantom and the SIMIND Monte Carlo simulation package were used to simulate full-time and half-time SPECT projection data. The data were reconstructed using the standard ordered subset expectation maximization-based algorithm and the new Monte Carlo-based algorithm. Defect contrast, myocardium versus ventricle contrast and resolution were calculated. In addition to the simulation studies, full-time and half-time SPECT projection data of 30 patients were reconstructed with the standard and the new method. The patient data were qualitatively evaluated by four nuclear medicine experts on a scale from 1 (poor quality) to 5 (high quality).

RESULTS

The new reconstruction method with half-time data produced higher contrast and better resolution in the simulations and also achieved higher qualitative scores in the patient study than the standard reconstruction with full-time data.

CONCLUSION

Half-time myocardial perfusion imaging using the new reconstruction algorithm with Monte Carlo-based scatter correction produced images with superior quality when compared with full-time imaging with standard reconstruction.

Tumor dosimetry and response for 153Sm-ethylenediamine tetramethylene phosphonic acid therapy of high-risk osteosarcoma

(153)Sm-ethylenediamine tetramethylene phosphonic acid ((153)Sm-EDTMP) therapy for osteosarcoma is being investigated. In this study, we analyzed the influence of (153)Sm-EDTMP administered activity (AA), osteosarcoma tumor density, mass, and the shape of the tumor on absorbed dose (AD). We also studied the biologic implication of the nonuniform tumor AD distribution using radiobiologic modeling and examined the relationship between tumor AD and response.

METHODS

Nineteen tumors in 6 patients with recurrent, refractory osteosarcoma enrolled in a phase I or II clinical trial of (153)Sm-EDTMP were analyzed using the 3-dimensional radiobiologic dosimetry (3D-RD) software package. Patients received a low dose of (153)Sm-EDTMP (37.0-51.8 MBq/kg), followed on hematologic recovery by a second, high dose (222 MBq/kg). Treatment response was evaluated using either CT or MRI after each therapy. SPECT/CT of the tumor regions were obtained at 4 and 48 h or 72 h after (153)Sm-EDTMP therapy for 3D-RD analysis. Mean tumor AD was also calculated using the OLINDA/EXM unit-density sphere model and was compared with the 3D-RD estimates.

RESULTS

On average, a 5-fold increase in the AA led to a 4-fold increase in the mean tumor AD over the high- versus low-dose-treated patients. The range of mean tumor AD and equivalent uniform dose (EUD) for low-dose therapy were 1.48-14.6 and 0.98-3.90 Gy, respectively. Corresponding values for high-dose therapy were 2.93-59.3 and 1.89-12.3 Gy, respectively. Mean tumor AD estimates obtained from OLINDA/EXM were within 5% of the mean AD values obtained using 3D-RD. On an individual tumor basis, both mean AD and EUD were positively related to percentage tumor volume reduction (P = 0.031 and 0.023, respectively).

CONCLUSION

The variations in tumor density, mass, and shape seen in these tumors did not affect the mean tumor AD estimation significantly. The tumor EUD was approximately 2- and 3-fold lower than the mean AD for low- and high-dose therapy, respectively. A dose-response relationship was observed for transient tumor volume shrinkage.

EQPlanar: a maximum-likelihood method for accurate organ activity estimation from whole body planar projections

Optimizing targeted radionuclide therapy requires patient-specific estimation of organ doses. The organ doses are estimated from quantitative nuclear medicine imaging studies, many of which involve planar whole body scans. We have previously developed the quantitative planar (QPlanar) processing method and demonstrated its ability to provide more accurate activity estimates than conventional geometric-mean-based planar (CPlanar) processing methods using physical phantom and simulation studies. The QPlanar method uses the maximum likelihood-expectation maximization algorithm, 3D organ volume of interests (VOIs), and rigorous models of physical image degrading factors to estimate organ activities. However, the QPlanar method requires alignment between the 3D organ VOIs and the 2D planar projections and assumes uniform activity distribution in each VOI. This makes application to patients challenging. As a result, in this paper we propose an extended QPlanar (EQPlanar) method that provides independent-organ rigid registration and includes multiple background regions. We have validated this method using both Monte Carlo simulation and patient data. In the simulation study, we evaluated the precision and accuracy of the method in comparison to the original QPlanar method. For the patient studies, we compared organ activity estimates at 24 h after injection with those from conventional geometric mean-based planar quantification using a 24 h post-injection quantitative SPECT reconstruction as the gold standard. We also compared the goodness of fit of the measured and estimated projections obtained from the EQPlanar method to those from the original method at four other time points where gold standard data were not available. In the simulation study, more accurate activity estimates were provided by the EQPlanar method for all the organs at all the time points compared with the QPlanar method. Based on the patient data, we concluded that the EQPlanar method provided a substantial increase in accuracy of organ activity estimates from 24 h planar images compared to the CPlanar using 24 h SPECT as the golden standard. For other time points, where no golden standard is available, better agreement between estimated and measured projections was observed by using the EQPlanar method compared to the QPlanar method. This phenomenon is consistent with the improvement in goodness of fit seen in both simulation data and 24 h patient data. Therefore, this indicates the improved reliability of organ activity estimates obtained though the EQPlanar method.

Development and evaluation of a model-based downscatter compensation method for quantitative I-131 SPECT

PURPOSE

The radionuclide 131I has found widespread use in targeted radionuclide therapy (TRT), partly due to the fact that it emits photons that can be imaged to perform treatment planning or posttherapy dose verification as well as beta rays that are suitable for therapy. In both the treatment planning and dose verification applications, it is necessary to estimate the activity distribution in organs or tumors at several time points. In vivo estimates of the 131I activity distribution at each time point can be obtained from quantitative single-photon emission computed tomography (QSPECT) images and organ activity estimates can be obtained either from QSPECT images or quantification of planar projection data. However, in addition to the photon used for imaging, 131I decay results in emission of a number of other higher-energy photons with significant abundances. These higher-energy photons can scatter in the body, collimator, or detector and be counted in the 364 keV photopeak energy window, resulting in reduced image contrast and degraded quantitative accuracy; these photons are referred to as downscatter. The goal of this study was to develop and evaluate a model-based downscatter compensation method specifically designed for the compensation of high-energy photons emitted by 131I and detected in the imaging energy window.

METHODS

In the evaluation study, we used a Monte Carlo simulation (MCS) code that had previously been validated for other radionuclides. Thus, in preparation for the evaluation study, we first validated the code for 131I imaging simulation by comparison with experimental data. Next, we assessed the accuracy of the downscatter model by comparing downscatter estimates with MCS results. Finally, we combined the downscatter model with iterative reconstruction-based compensation for attenuation (A) and scatter (S) and the full (D) collimator-detector response of the 364 keV photons to form a comprehensive compensation method. We evaluated this combined method in terms of quantitative accuracy using the realistic 3D NCAT phantom and an activity distribution obtained from patient studies. We compared the accuracy of organ activity estimates in images reconstructed with and without addition of downscatter compensation from projections with and without downscatter contamination.

RESULTS

We observed that the proposed method provided substantial improvements in accuracy compared to no downscatter compensation and had accuracies comparable to reconstructions from projections without downscatter contamination.

CONCLUSIONS

The results demonstrate that the proposed model-based downscatter compensation method is effective and may have a role in quantitative 131I imaging.

Reduction of collimator correction artefacts with bayesian reconstruction in spect

Poor resolution of single photon emission computed tomography (SPECT) has degraded its use in clinical practice. Collimator correction has been shown to improve the reconstructed resolution, but the correction can generate ringing artefacts, which lower image quality. This paper investigates whether Bayesian reconstruction methods could reduce these artefacts. We have applied and tested three Bayesian reconstruction methods: smoothing prior, median root prior, and anatomical prior. To demonstrate the efficacy of these methods, we compared their physical and visual performance both in phantom and patient studies. All the three Bayesian reconstruction methods reduced the collimator correction artefacts. Images reconstructed using the smoothing prior and the median root prior had slightly lower contrast than the standard reconstruction with collimator correction, whereas the anatomical prior produced images with good resolution and contrast.

Image quantification for radiation dose calculations--limitations and uncertainties

Radiation dose calculations in nuclear medicine depend on quantification of activity via planar and/or tomographic imaging methods. However, both methods have inherent limitations, and the accuracy of activity estimates varies with object size, background levels, and other variables. The goal of this study was to evaluate the limitations of quantitative imaging with planar and single photon emission computed tomography (SPECT) approaches, with a focus on activity quantification for use in calculating absorbed dose estimates for normal organs and tumors. To do this we studied a series of phantoms of varying complexity of geometry, with three radionuclides whose decay schemes varied from simple to complex. Four aqueous concentrations of ⁹⁹mTc, ¹³¹I, and ¹¹¹In (74, 185, 370, and 740 kBq mL⁻¹) were placed in spheres of four different sizes in a water-filled phantom, with three different levels of activity in the surrounding water. Planar and SPECT images of the phantoms were obtained on a modern SPECT/computed tomography (CT) system. These radionuclides and concentration/background studies were repeated using a cardiac phantom and a modified torso phantom with liver and "tumor" regions containing the radionuclide concentrations and with the same varying background levels. Planar quantification was performed using the geometric mean approach, with attenuation correction (AC), and with and without scatter corrections (SC and NSC). SPECT images were reconstructed using attenuation maps (AM) for AC; scatter windows were used to perform SC during image reconstruction. For spherical sources with corrected data, good accuracy was observed (generally within ±10% of known values) for the largest sphere (11.5 mL) and for both planar and SPECT methods with ⁹⁹mTc and ¹³¹I, but were poorest and deviated from known values for smaller objects, most notably for ¹¹¹In. SPECT quantification was affected by the partial volume effect in smaller objects and generally showed larger errors than the planar results in these cases for all radionuclides. For the cardiac phantom, results were the most accurate of all of the experiments for all radionuclides. Background subtraction was an important factor influencing these results. The contribution of scattered photons was important in quantification with ¹³¹I; if scatter was not accounted for, activity tended to be overestimated using planar quantification methods. For the torso phantom experiments, results show a clear underestimation of activity when compared to previous experiment with spherical sources for all radionuclides. Despite some variations that were observed as the level of background increased, the SPECT results were more consistent across different activity concentrations. Planar or SPECT quantification on state-of-the-art gamma cameras with appropriate quantitative processing can provide accuracies of better than 10% for large objects and modest target-to-background concentrations; however when smaller objects are used, in the presence of higher background, and for nuclides with more complex decay schemes, SPECT quantification methods generally produce better results.

A treatment planning method for sequentially combining radiopharmaceutical therapy and external radiation therapy

PURPOSE

Effective cancer treatment generally requires combination therapy. The combination of external beam therapy (XRT) with radiopharmaceutical therapy (RPT) requires accurate three-dimensional dose calculations to avoid toxicity and evaluate efficacy. We have developed and tested a treatment planning method, using the patient-specific three-dimensional dosimetry package 3D-RD, for sequentially combined RPT/XRT therapy designed to limit toxicity to organs at risk.

METHODS AND MATERIALS

The biologic effective dose (BED) was used to translate voxelized RPT absorbed dose (D(RPT)) values into a normalized total dose (or equivalent 2-Gy-fraction XRT absorbed dose), NTD(RPT) map. The BED was calculated numerically using an algorithmic approach, which enabled a more accurate calculation of BED and NTD(RPT). A treatment plan from the combined Samarium-153 and external beam was designed that would deliver a tumoricidal dose while delivering no more than 50 Gy of NTD(sum) to the spinal cord of a patient with a paraspinal tumor.

RESULTS

The average voxel NTD(RPT) to tumor from RPT was 22.6 Gy (range, 1-85 Gy); the maximum spinal cord voxel NTD(RPT) from RPT was 6.8 Gy. The combined therapy NTD(sum) to tumor was 71.5 Gy (range, 40-135 Gy) for a maximum voxel spinal cord NTD(sum) equal to the maximum tolerated dose of 50 Gy.

CONCLUSIONS

A method that enables real-time treatment planning of combined RPT-XRT has been developed. By implementing a more generalized conversion between the dose values from the two modalities and an activity-based treatment of partial volume effects, the reliability of combination therapy treatment planning has been expanded.

The effect of volume-of-interest misregistration on quantitative planar activity and dose estimation

In targeted radionuclide therapy (TRT), dose estimation is essential for treatment planning and tumor dose response studies. Dose estimates are typically based on a time series of whole-body conjugate view planar or SPECT scans of the patient acquired after administration of a planning dose. Quantifying the activity in the organs from these studies is an essential part of dose estimation. The quantitative planar (QPlanar) processing method involves accurate compensation for image degrading factors and correction for organ and background overlap via the combination of computational models of the image formation process and 3D volumes of interest defining the organs to be quantified. When the organ VOIs are accurately defined, the method intrinsically compensates for attenuation, scatter and partial volume effects, as well as overlap with other organs and the background. However, alignment between the 3D organ volume of interest (VOIs) used in QPlanar processing and the true organ projections in the planar images is required. The aim of this research was to study the effects of VOI misregistration on the accuracy and precision of organ activity estimates obtained using the QPlanar method. In this work, we modeled the degree of residual misregistration that would be expected after an automated registration procedure by randomly misaligning 3D SPECT/CT images, from which the VOI information was derived, and planar images. Mutual information-based image registration was used to align the realistic simulated 3D SPECT images with the 2D planar images. The residual image misregistration was used to simulate realistic levels of misregistration and allow investigation of the effects of misregistration on the accuracy and precision of the QPlanar method. We observed that accurate registration is especially important for small organs or ones with low activity concentrations compared to neighboring organs. In addition, residual misregistration gave rise to a loss of precision in the activity estimates that was on the order of the loss of precision due to Poisson noise in the projection data. These results serve as a lower bound on the effects of misregistration on the accuracy and precision of QPlanar activity estimate and demonstrate that misregistration errors must be taken into account when assessing the overall precision of organ dose estimates.

Effects of shortened acquisition time on accuracy and precision of quantitative estimates of organ activity

PURPOSE

Quantitative estimation of in vivo organ uptake is an essential part of treatment planning for targeted radionuclide therapy. This usually involves the use of planar or SPECT scans with acquisition times chosen based more on image quality considerations rather than the minimum needed for precise quantification. In previous simulation studies at clinical count levels (185 MBq 111In), the authors observed larger variations in accuracy of organ activity estimates resulting from anatomical and uptake differences than statistical noise. This suggests that it is possible to reduce the acquisition time without substantially increasing the variation in accuracy.

METHODS

To test this hypothesis, the authors compared the accuracy and variation in accuracy of organ activity estimates obtained from planar and SPECT scans at various count levels. A simulated phantom population with realistic variations in anatomy and biodistribution was used to model variability in a patient population. Planar and SPECT projections were simulated using previously validated Monte Carlo simulation tools. The authors simulated the projections at count levels approximately corresponding to 1.5-30 min of total acquisition time. The projections were processed using previously described quantitative SPECT (QSPECT) and planar (QPlanar) methods. The QSPECT method was based on the OS-EM algorithm with compensations for attenuation, scatter, and collimator-detector response. The QPlanar method is based on the ML-EM algorithm using the same model-based compensation for all the image degrading effects as the QSPECT method. The volumes of interests (VOIs) were defined based on the true organ configuration in the phantoms. The errors in organ activity estimates from different count levels and processing methods were compared in terms of mean and standard deviation over the simulated phantom population.

RESULTS

There was little degradation in quantitative reliability when the acquisition time was reduced by half for the QSPECT method (the mean error changed by < 1%, e.g., 0.9%-0.3% = 0.6% for the spleen). The magnitude of the errors and variations in errors for large organ with high uptake were still acceptable for 1.5 min scans, even though the errors were slightly larger than those for the 30 min scans (i.e., < 2% for liver, < 3% for heart). The errors over the ranges of scan times studied for the QPlanar method were all within 0.3% for all organs.

CONCLUSIONS

These data indicate that, for the purposes of organ activity estimation, acquisition times could be reduced at least by a factor of 2 for the QSPECT and QPlanar methods with little effect on the errors in organ activity estimates. The acquisition time can be further reduced for the QPlanar method, assuming well-registered VOIs are available and the activity distribution in organs can be treated as uniform. Although the differences in accuracy and precision were statistically significant for all the acquisition times shorter than 30 min, the magnitude of the changes in accuracy and precision were small and likely not clinically important. The reduction in SPECT acquisition time justified by this study makes the use of SPECT a more clinically practical alternative to conventional planar scanning for targeted radiotherapy treatment planning.

The impact of 3D volume of interest definition on accuracy and precision of activity estimation in quantitative SPECT and planar processing methods

Accurate and precise estimation of organ activities is essential for treatment planning in targeted radionuclide therapy. We have previously evaluated the impact of processing methodology, statistical noise and variability in activity distribution and anatomy on the accuracy and precision of organ activity estimates obtained with quantitative SPECT (QSPECT) and planar (QPlanar) processing. Another important factor impacting the accuracy and precision of organ activity estimates is accuracy of and variability in the definition of organ regions of interest (ROI) or volumes of interest (VOI). The goal of this work was thus to systematically study the effects of VOI definition on the reliability of activity estimates. To this end, we performed Monte Carlo simulation studies using randomly perturbed and shifted VOIs to assess the impact on organ activity estimates. The 3D NCAT phantom was used with activities that modeled clinically observed (111)In ibritumomab tiuxetan distributions. In order to study the errors resulting from misdefinitions due to manual segmentation errors, VOIs of the liver and left kidney were first manually defined. Each control point was then randomly perturbed to one of the nearest or next-nearest voxels in three ways: with no, inward or outward directional bias, resulting in random perturbation, erosion or dilation, respectively, of the VOIs. In order to study the errors resulting from the misregistration of VOIs, as would happen, e.g. in the case where the VOIs were defined using a misregistered anatomical image, the reconstructed SPECT images or projections were shifted by amounts ranging from -1 to 1 voxels in increments of with 0.1 voxels in both the transaxial and axial directions. The activity estimates from the shifted reconstructions or projections were compared to those from the originals, and average errors were computed for the QSPECT and QPlanar methods, respectively. For misregistration, errors in organ activity estimations were linear in the shift for both the QSPECT and QPlanar methods. QPlanar was less sensitive to object definition perturbations than QSPECT, especially for dilation and erosion cases. Up to 1 voxel misregistration or misdefinition resulted in up to 8% error in organ activity estimates, with the largest errors for small or low uptake organs. Both types of VOI definition errors produced larger errors in activity estimates for a small and low uptake organs (i.e. -7.5% to 5.3% for the left kidney) than for a large and high uptake organ (i.e. -2.9% to 2.1% for the liver). We observed that misregistration generally had larger effects than misdefinition, with errors ranging from -7.2% to 8.4%. The different imaging methods evaluated responded differently to the errors from misregistration and misdefinition. We found that QSPECT was more sensitive to misdefinition errors, but less sensitive to misregistration errors, as compared to the QPlanar method. Thus, sensitivity to VOI definition errors should be an important criterion in evaluating quantitative imaging methods.

Arterial wall dosimetry for non-Hodgkin lymphoma patients treated with radioimmunotherapy

Tumors in non-Hodgkin lymphoma (NHL) patients are often proximal to the major blood vessels in the abdomen or neck. In external-beam radiotherapy, these tumors present a challenge because imaging resolution prevents the beam from being targeted to the tumor lesion without also irradiating the artery wall. This problem has led to potentially life-threatening delayed toxicity. Because radioimmunotherapy has resulted in long-term survival of NHL patients, we investigated whether the absorbed dose (AD) to the artery wall in radioimmunotherapy of NHL is of potential concern for delayed toxicity. SPECT resolution is not sufficient to enable dosimetric analysis of anatomic features of the thickness of the aortic wall. Therefore, we present a model of aortic wall toxicity based on data from 4 patients treated with (131)I-tositumomab.

METHODS

Four NHL patients with periaortic tumors were administered pretherapeutic (131)I-tositumomab. Abdominal SPECT and whole-body planar images were obtained at 48, 72, and 144 h after tracer administration. Blood-pool activity concentrations were obtained from regions of interest drawn on the heart on the planar images. Tumor and blood activity concentrations, scaled to therapeutic administered activities-both standard and myeloablative-were input into a geometry and tracking model (GEANT, version 4) of the aorta. The simulated energy deposited in the arterial walls was collected and fitted, and the AD and biologic effective dose values to the aortic wall and tumors were obtained for standard therapeutic and hypothetical myeloablative administered activities.

RESULTS

Arterial wall ADs from standard therapy were lower (0.6-3.7 Gy) than those typical from external-beam therapy, as were the tumor ADs (1.4-10.5 Gy). The ratios of tumor AD to arterial wall AD were greater for radioimmunotherapy by a factor of 1.9-4.0. For myeloablative therapy, artery wall ADs were in general less than those typical for external-beam therapy (9.4-11.4 Gy for 3 of 4 patients) but comparable for 1 patient (32.6 Gy).

CONCLUSION

Blood vessel radiation dose can be estimated using the software package 3D-RD combined with GEANT modeling. The dosimetry analysis suggested that arterial wall toxicity is highly unlikely in standard dose radioimmunotherapy but should be considered a potential concern and limiting factor in myeloablative therapy.

Evaluation of a method for activity estimation in Sm-153 EDTMP imaging

Absolute activity evaluation is fundamental for internal radionuclide dosimetry when patient-specific therapy optimization is wanted. Often, quantification is attempted with 3D SPECT image based (IB) methods, but the true concentration values can be underestimated due to the partial volume effect (PVE). This is especially true when small diffuse lesions are present. In this paper, we describe a 3D region of interest (ROI) based quantification method (LS-ROI), which estimates the ROI concentration values directly from the projection data acquired in the tomographic scan once ROIs have been segmented on a CT and/or a SPECT image. The method, which has inherent PVE correction capabilities, was applied both on simulated and on real phantom data. Simulations reflected the case of a patient with bone metastases treated with 153Sm-EDTMP: Both the activity in the metastases and the total retention in the skeleton were evaluated. Thirty noisy data sets were produced in order to evaluate the accuracy and precision of the method. The effect of region segmentation errors on estimated concentrations was thoroughly investigated. Real data were acquired on a NEMA phantom, where a cylindrical central region (283 cm3) simulated the bone and two spheres (10.3 and 25.5 cm3) simulated the metastases. The results obtained with the LS-ROI method were compared with those of a conventional 3D IB method and those of a quantitative conjugate view approach derived from LS-ROI and applied to the anterior and posterior views acquired in the tomographic scan (LS-ROI anterior-posterior: LS-ROI-AP). Simulations showed that when the geometry of regions is known, the LS-ROI method recovered the simulated concentration values within 20%, while the IB method underestimated the concentration in high activity small lesions by as much as 49%. Segmentation errors, up to 44% of the true region volume, produced a higher variation in LS-ROI estimates than in IB ones; however, the overall bias of the LS-ROI estimates (< or = 25%) remained lower than that of IB estimates. In the case of the evaluation of the total retention in the skeleton, the LS-ROI method recovered the simulated value within 2%, while IB underestimated it up to 13%. In all the cases, the LS-ROI-AP method showed an accuracy comparable with that of the LS-ROI one, and a worse precision just because of the lower number of counts used in the analysis. However, a worsening of LS-ROI-AP performances was demonstrated in the case of strong overlap of regions: In this case, a bias of up to 40% was observed. The results obtained on real phantom data confirmed the simulation results: The IB method underestimated activity up to 47% in the smallest sphere, while the bias was reduced to 13% with LS-ROI and LS-ROI-AP estimates. The good quantification capabilities of the LS-ROI method can be useful for absolute activity quantification in the case of small active diffused lesions and constitute the basis for the development of an accurate patient-specific planning strategy in internal radionuclide treatments, provided there is a reliable segmentation of lesions.

Comparison of organ residence time estimation methods for radioimmunotherapy dosimetry and treatment planning--patient studies

The estimation of organ residence time is essential for high-dose myeloablative regimens in radioimmunotherapy (RIT). Frequently, this estimation is based on a series of simple planar scans and planar processing. The authors previously performed a simulation study which demonstrated that the accuracy of this methodology is limited compared to a hybrid planar/SPECT residence time estimation method. In this work the authors applied this hybrid method to data from a clinical trial of high-dose myeloablative yttrium-90 ibritumomab tiuxetan therapy. Image data acquired from 18 patients were comprised of planar scans at five time points ranging from 1 to 144 h postinjection and abdominal and thoracic SPECT/CT scans obtained at 24 h postinjection. The simple planar processing method used in this work was based on the geometric mean method with energy window based scatter compensation. No explicit background subtraction nor object or source thickness corrections were performed. The SPECT projections were reconstructed using iterative reconstruction with compensations for attenuation, scatter, and full collimator-detector response. Large differences were observed when residence times were estimated using the simple planar method compared to the hybrid method. The differences were not constant but varied in magnitude and sign. For the dose-limiting organ (liver), the average difference was -18% and variation in the difference was 19%, similar to the differences observed in a previously reported simulation study. The authors also looked at the relationship between the weight of the patient and the liver residence time and found that there was no meaningful correlation for either method. This indicates that weight would not be an adequate proxy for an experimental estimate of residence time when choosing the activity to administer for therapy. The authors conclude that methods such as the simple planar method used here are inadequate for RIT treatment planning. More sophisticated methods, such as the hybrid SPECT/planar method investigated here, are likely to be better predictors of organ dose and, as a result, organ toxicities.

Striatal dopamine D2 binding is related to frontal BOLD response during updating of long-term memory representations

Multi-modal brain imaging was used to examine the relation between individual differences in resting-state striatal dopamine D2 binding and the magnitude of prefrontal BOLD activation during updating of long-term memory (LTM) representations. Increased activity in the left prefrontal cortex was observed when LTM updating was required, and there was a positive correlation between striatal D2 activity and the magnitude of left prefrontal activity during updating. These findings support predictions from neurocomputational models of a relation of dopaminergic neurotransmission to transient cognitive operations and related brain activity.

Evaluation of quantitative imaging methods for organ activity and residence time estimation using a population of phantoms having realistic variations in anatomy and uptake

Estimating organ residence times is an essential part of patient-specific dosimetry for radioimmunotherapy (RIT). Quantitative imaging methods for RIT are often evaluated using a single physical or simulated phantom but are intended to be applied clinically where there is variability in patient anatomy, biodistribution, and biokinetics. To provide a more relevant evaluation, the authors have thus developed a population of phantoms with realistic variations in these factors and applied it to the evaluation of quantitative imaging methods both to find the best method and to demonstrate the effects of these variations. Using whole body scans and SPECT/CT images, organ shapes and time-activity curves of 111In ibritumomab tiuxetan were measured in dosimetrically important organs in seven patients undergoing a high dose therapy regimen. Based on these measurements, we created a 3D NURBS-based cardiac-torso (NCAT)-based phantom population. SPECT and planar data at realistic count levels were then simulated using previously validated Monte Carlo simulation tools. The projections from the population were used to evaluate the accuracy and variation in accuracy of residence time estimation methods that used a time series of SPECT and planar scans, Quantitative SPECT (QSPECT) reconstruction methods were used that compensated for attenuation, scatter, and the collimator-detector response. Planar images were processed with a conventional (CPlanar) method that used geometric mean attenuation and triple-energy window scatter compensation and a quantitative planar (QPlanar) processing method that used model-based compensation for image degrading effects. Residence times were estimated from activity estimates made at each of five time points. The authors also evaluated hybrid methods that used CPlanar or QPlanar time-activity curves rescaled to the activity estimated from a single QSPECT image. The methods were evaluated in terms of mean relative error and standard deviation of the relative error in the residence time estimates taken over the phantom population. The mean errors in the residence time estimates over all the organs were < 9.9% (pure QSPECT), < 13.2% (pure QPLanar), < 7.2% (hybrid QPlanar/QSPECT), < 19.2% (hybrid CPlanar/QSPECT), and 7%-159% (pure CPlanar). The standard deviations of the errors for all the organs over all the phantoms were < 9.9%, < 11.9%, < 10.8%, < 22.0%, and < 107.9% for the same methods, respectively. The processing methods differed both in terms of their average accuracy and the variation of the accuracy over the population of phantoms, thus demonstrating the importance of using a phantom population in evaluating quantitative imaging methods. Hybrid CPlanar/QSPECT provided improved accuracy compared to pure CPlanar and required the addition of only a single SPECT acquisition. The QPlanar or hybrid QPlanar/QSPECT methods had mean errors and standard deviations of errors that approached those of pure QSPECT while providing simplified image acquisition protocols, and thus may be more clinically practical.

Three-dimensional imaging-based radiobiological dosimetry

Targeted radionuclide therapy holds promise as a new treatment for cancer. Advances in imaging are making it possible for researchers to evaluate the spatial distribution of radioactivity in tumors and normal organs over time. Matched anatomical imaging, such as combined single-photon emission computed tomography/computed tomography and positron emission tomography/computed tomography, has also made it possible to obtain tissue density information in conjunction with the radioactivity distribution. Coupled with sophisticated iterative reconstruction algorithms, these advances have made it possible to perform highly patient-specific dosimetry that also incorporates radiobiological modeling. Such sophisticated dosimetry techniques are still in the research investigation phase. Given the attendant logistical and financial costs, a demonstrated improvement in patient care will be a prerequisite for the adoption of such highly-patient specific internal dosimetry methods.

Accuracy of quantitative reconstructions in SPECT/CT imaging

The goal of this study was to determine the quantitative accuracy of our OSEM-APDI reconstruction method based on SPECT/CT imaging for Tc-99m, In-111, I-123, and I-131 isotopes. Phantom studies were performed on a SPECT/low-dose multislice CT system (Infinia-Hawkeye-4 slice, GE Healthcare) using clinical acquisition protocols. Two radioactive sources were centrally and peripherally placed inside an anthropometric Thorax phantom filled with non-radioactive water. Corrections for attenuation, scatter, collimator blurring and collimator septal penetration were applied and their contribution to the overall accuracy of the reconstruction was evaluated. Reconstruction with the most comprehensive set of corrections resulted in activity estimation with error levels of 3-5% for all the isotopes.

Improved activity estimation with MC-JOSEM versus TEW-JOSEM in 111In SPECT

We have previously developed a fast Monte Carlo (MC)-based joint ordered-subset expectation maximization (JOSEM) iterative reconstruction algorithm, MC-JOSEM. A phantom study was performed to compare quantitative imaging performance of MC-JOSEM with that of a triple-energy-window approach (TEW) in which estimated scatter was also included additively within JOSEM, TEW-JOSEM. We acquired high-count projections of a 5.5 cm3 sphere of 111In at different locations in the water-filled torso phantom; high-count projections were then obtained with 111In only in the liver or only in the soft-tissue background compartment, so that we could generate synthetic projections for spheres surrounded by various activity distributions. MC scatter estimates used by MC-JOSEM were computed once after five iterations of TEW-JOSEM. Images of different combinations of liver/background and sphere/background activity concentration ratios were reconstructed by both TEW-JOSEM and MC-JOSEM for 40 iterations. For activity estimation in the sphere, MC-JOSEM always produced better relative bias and relative standard deviation than TEW-JOSEM for each sphere location, iteration number, and activity combination. The average relative bias of activity estimates in the sphere for MC-JOSEM after 40 iterations was -6.9%, versus -15.8% for TEW-JOSEM, while the average relative standard deviation of the sphere activity estimates was 16.1% for MC-JOSEM, versus 27.4% for TEW-JOSEM. Additionally, the average relative bias of activity concentration estimates in the liver and the background for MC-JOSEM after 40 iterations was -3.9%, versus -12.2% for TEW-JOSEM, while the average relative standard deviation of these estimates was 2.5% for MC-JOSEM, versus 3.4% for TEW-JOSEM. MC-JOSEM is a promising approach for quantitative activity estimation in 111In SPECT.

Comparison of residence time estimation methods for radioimmunotherapy dosimetry and treatment planning--Monte Carlo simulation studies

Estimating the residence times in tumor and normal organs is an essential part of treatment planning for radioimmunotherapy (RIT). This estimation is usually done using a conjugate view whole body scan time series and planar processing. This method has logistical and cost advantages compared to 3-D imaging methods such as Single photon emission computed tomography (SPECT), but, because it does not provide information about the 3-D distribution of activity, it is difficult to fully compensate for effects such as attenuation and background and overlapping activity. Incomplete compensation for these effects reduces the accuracy of the residence time estimates. In this work we compare residence times estimates obtained using planar methods to those from methods based on quantitative SPECT (QSPECT) reconstructions. We have previously developed QSPECT methods that provide compensation for attenuation, scatter, collimator-detector response, and partial volume effects. In this study we compared the use of residence time estimation methods using QSPECT to planar methods. The evaluation was done using the realistic NCAT phantom with organ time activities that model (111)In ibritumomab tiuxetan. Projection data were obtained using Monte Carlo simulations (MCS) that realistically model the image formation process including penetration and scatter in the collimator-detector system. These projection data were used to evaluate the accuracy of residence time estimation using a time series of QSPECT studies, a single QSPECT study combined with planar scans and the planar scans alone. The errors in the residence time estimates were 3.8%, 15%, and 2%-107% for the QSPECT, hybrid planar/QSPECT, and planar methods, respectively. The quantitative accuracy was worst for pure planar processing and best for pure QSPECT processing. Hybrid planar/QSPECT methods, where a single QSPECT study was combined with a series of planar scans, provided a large and statistically significant improvement in quantitative accuracy for most organs compared to the planar scans alone, even without sophisticated attention to background subtraction or thickness corrections in planar processing. These results indicate that hybrid planar/QSPECT methods are generally superior to pure planar methods and may be an acceptable alternative to performing a time series of QSPECT studies.

Tumor dosimetry on SPECT (186)Re-HEDP scans: variations in the results from the reconstruction methods used

The aim of this work was to estimate tumor-absorbed doses delivered from the administration of fixed activities of (186)Re-HEDP for the treatment of bone metastases from prostate cancer. The variations and reproducibility in the estimated absorbed dose owing to the reconstruction algorithm used (OSEM vs. FBP) were also analysed. For this aim, correction methods for scatter and attenuation were kept identical, whereas the same calibration data and thresholding techniques were used to obtain quantification. This study was carried out in spinal and pelvic lesions of 7 patients. For comparison, the absorbed doses, based upon the maximum and the mean voxel count, were calculated, resulting in the absorbed dose (maximum): 60 Gy (sigma = 30 Gy) and 33 Gy (sigma = 15 Gy) for OSEM and FBP, respectively, and absorbed dose (mean): 26 Gy (sigma = 12 Gy) and 17 Gy (sigma = 7 Gy) with OSEM and FBP, respectively. We concluded that the administration of fixed activity resulted in a range of absorbed doses, and we showed that, despite using the same approach, the choice of the reconstruction algorithm can result in differences higher than 50% in the estimated tumor-absorbed doses. In conclusion, the need for a standardization of the methodology used for the calculations is emphasized by this work, especially when comparisons between patients and different centers are of interest.

Fast hybrid SPECT simulation including efficient septal penetration modelling (SP-PSF)

Single photon emission computed tomography (SPECT) images are degraded by the detection of scattered photons and photons that penetrate the collimator septa. In this paper, a previously proposed Monte Carlo software that employs fast object scatter simulation using convolution-based forced detection (CFD) is extended towards a wide range of medium and high energy isotopes measured using various collimators. To this end, a fast method was developed for incorporating effects of septal penetrating (SP) photons. The SP contributions are obtained by calculating the object attenuation along the path from primary emission to detection followed by sampling a pre-simulated and scalable septal penetration point spread function (SP-PSF). We found that with only a very slight reduction in accuracy, we could accelerate the SP simulation by four orders of magnitude. To achieve this, we combined: (i) coarse sampling of the activity and attenuation distribution; (ii) simulation of the penetration only for a coarse grid of detector pixels followed by interpolation and (iii) neglection of SP-PSF elements below a certain threshold. By inclusion of this SP-PSF-based simulation it became possible to model both primary and septal penetrated photons while only 10% extra computation time was added to the CFD-based Monte Carlo simulator. As a result, a SPECT simulation of a patient-like distribution including SP now takes less than 5 s per projection angle on a dual processor PC. Therefore, the simulator is well-suited as an efficient projector for fully 3D model-based reconstruction or as a fast data-set generator for applications such as image processing optimization or observer studies.

Implementation of an iterative scatter correction, the influence of attenuation map quality and their effect on absolute quantitation in SPECT

We investigated the accuracy of qSPECT, a quantitative SPECT reconstruction algorithm we have developed which employs corrections for collimator blurring, photon attenuation and scatter, and provides images in units of absolute radiotracer concentrations (kBq cm(-3)). Using simulated and experimental phantom data with characteristics similar to clinical cardiac perfusion data, we studied the implementation of a scatter correction (SC) as part of an iterative reconstruction protocol. Additionally, with experimental phantom studies we examined the influence of CT-based attenuation maps, relative to those obtained from conventional SPECT transmission scans, on SCs and quantitation. Our results indicate that the qSPECT estimated scatter corrections did not change appreciably after the third iteration of the reconstruction. For the simulated data, qSPECT concentrations agreed with images reconstructed using ideal, scatter-free, simulated data to within 6%. For the experimental data, we observed small systematic differences in the scatter fractions for data using different combinations of SCs and attenuation maps. The SCs were found to be significantly influenced by errors in image coregistration. The reconstructed concentrations using CT-based corrections were more quantitatively accurate than those using attenuation maps from conventional SPECT transmission scans. However, segmenting the attenuation maps from SPECT transmission scans could provide sufficient accuracy for most applications.

Correction for scatter and septal penetration using convolution subtraction methods and model-based compensation in 123I brain SPECT imaging-a Monte Carlo study

Scatter and septal penetration deteriorate contrast and quantitative accuracy in single photon emission computed tomography (SPECT). In this study four different correction techniques for scatter and septal penetration are evaluated for 123I brain SPECT. One of the methods is a form of model-based compensation which uses the effective source scatter estimation (ESSE) for modelling scatter, and collimator-detector response (CDR) including both geometric and penetration components. The other methods, which operate on the 2D projection images, are convolution scatter subtraction (CSS) and two versions of transmission dependent convolution subtraction (TDCS), one of them proposed by us. This method uses CSS for correction for septal penetration, with a separate kernel, and TDCS for scatter correction. The corrections are evaluated for a dopamine transporter (DAT) study and a study of the regional cerebral blood flow (rCBF), performed with 123I. The images are produced using a recently developed Monte Carlo collimator routine added to the program SIMIND which can include interactions in the collimator. The results show that the method included in the iterative reconstruction is preferable to the other methods and that the new TDCS version gives better results compared with the other 2D methods.

Comparison of conventional, model-based quantitative planar, and quantitative SPECT image processing methods for organ activity estimation using In-111 agents

Accurate quantification of organ radionuclide uptake is important for patient-specific dosimetry. The quantitative accuracy from conventional conjugate view methods is limited by overlap of projections from different organs and background activity, and attenuation and scatter. In this work, we propose and validate a quantitative planar (QPlanar) processing method based on maximum likelihood (ML) estimation of organ activities using 3D organ VOIs and a projector that models the image degrading effects. Both a physical phantom experiment and Monte Carlo simulation (MCS) studies were used to evaluate the new method. In these studies, the accuracies and precisions of organ activity estimates for the QPlanar method were compared with those from conventional planar (CPlanar) processing methods with various corrections for scatter, attenuation and organ overlap, and a quantitative SPECT (QSPECT) processing method. Experimental planar and SPECT projections and registered CT data from an RSD Torso phantom were obtained using a GE Millenium VH/Hawkeye system. The MCS data were obtained from the 3D NCAT phantom with organ activity distributions that modelled the uptake of (111)In ibritumomab tiuxetan. The simulations were performed using parameters appropriate for the same system used in the RSD torso phantom experiment. The organ activity estimates obtained from the CPlanar, QPlanar and QSPECT methods from both experiments were compared. From the results of the MCS experiment, even with ideal organ overlap correction and background subtraction, CPlanar methods provided limited quantitative accuracy. The QPlanar method with accurate modelling of the physical factors increased the quantitative accuracy at the cost of requiring estimates of the organ VOIs in 3D. The accuracy of QPlanar approached that of QSPECT, but required much less acquisition and computation time. Similar results were obtained from the physical phantom experiment. We conclude that the QPlanar method, based on 3D organ VOIs and accurate models of the projection process, provided a substantial increase in accuracy of organ activity estimates from planar images compared to CPlanar processing and had accuracy approaching that of QSPECT.