Accurate dosimetry in 131I radionuclide therapy using patient-specific, 3-dimensional methods for SPECT reconstruction and absorbed dose calculation

Abstract

(131)I radionuclide therapy studies have not shown a strong relationship between tumor absorbed dose and response, possibly due to inaccuracies in activity quantification and dose estimation. The goal of this work was to establish the accuracy of (131)I activity quantification and absorbed dose estimation when patient-specific, 3-dimensional (3D) methods are used for SPECT reconstruction and for absorbed dose calculation.

Methods: Clinically realistic voxel-phantom simulations were used in the evaluation of activity quantification and dosimetry. SPECT reconstruction was performed using an ordered-subsets expectation maximization (OSEM) algorithm with compensation for scatter, attenuation, and 3D detector response. Based on the SPECT image and a patient-specific density map derived from CT, 3D dosimetry was performed using a newly implemented Monte Carlo code. Dosimetry was evaluated by comparing mean absorbed dose estimates calculated directly from the defined phantom activity map with those calculated from the SPECT image of the phantom. Finally, the 3D methods were applied to a radioimmunotherapy patient, and the mean tumor absorbed dose from the new calculation was compared with that from conventional dosimetry obtained from conjugate-view imaging.

Results: Overall, the accuracy of the SPECT-based absorbed dose estimates in the phantom was >12% for targets down to 16 mL and up to 35% for the smallest 7-mL tumor. To improve accuracy in the smallest tumor, more OSEM iterations may be needed. The relative SD from multiple realizations was <3% for all targets except for the smallest tumor. For the patient, the mean tumor absorbed dose estimate from the new Monte Carlo calculation was 7% higher than that from conventional dosimetry.

Conclusion: For target sizes down to 16 mL, highly accurate and precise dosimetry can be obtained with 3D methods for SPECT reconstruction and absorbed dose estimation. In the future, these methods can be applied to patients to potentially establish correlations between tumor regression and the absorbed dose statistics from 3D dosimetry.

FIGURE 1

Two typical slices of voxel phantom show superimposed tumors. (A) Defined activity map. (B) Defined density map. (C) SPECT image. (D) Absorbed dose-rate map.

FIGURE 2

Evaluation procedure.

FIGURE 3

Counts in target VOI plotted as function of iteration number for HE collimator (A) and UHE collimator (B).

FIGURE 4

One slice of RIT patient CT image, SPECT image, and corresponding dose-rate map.

FIGURE 5

Differential dose-volume histogram corresponding to tumor outlined in patient image of Figure 4.