Development and evaluation of an improved quantitative 90Y bremsstrahlung SPECT method

Abstract

Yttrium-90 ((90)Y) is one of the most commonly used radionuclides in targeted radionuclide therapy (TRT). Since it decays with essentially no gamma photon emissions, surrogate radionuclides (e.g., (111)In) or imaging agents (e.g., (99m)Tc MAA) are typically used for treatment planning. It would, however, be useful to image (90)Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the (90)Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of (90)Y bremsstrahlung photons using single photon emission computed tomography (SPECT) imaging. The continuous and broad energy distribution of bremsstrahlung photons, however, imposes substantial challenges on accurate quantification of the activity distribution. The aim of this work was to develop and evaluate an improved quantitative (90)Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. Accurate modeling of image degrading factors such as object attenuation and scatter and the collimator-detector response is essential to obtain quantitatively accurate images. All of the image degrading factors are energy dependent. Thus, the authors separated the modeling of the bremsstrahlung photons into multiple categories and energy ranges. To improve the accuracy, the authors used a bremsstrahlung energy spectrum previously estimated from experimental measurements and incorporated a model of the distance between (90)Y decay location and bremsstrahlung emission location into the SIMIND code used to generate the response functions and kernels used in the model. This improved Monte Carlo bremsstrahlung simulation was validated by comparison to experimentally measured projection data of a (90)Y line source. The authors validated the accuracy of the forward projection model for photons in the various categories and energy ranges using the validated Monte Carlo (MC) simulation method. The forward projection model was incorporated into an iterative ordered subsets-expectation maximization (OS-EM) reconstruction code to allow for quantitative SPECT reconstruction. The resulting code was validated using both a physical phantom experiment with spherical objects in a warm background and a realistic anatomical phantom simulation. In the physical phantom study, the authors evaluated the method in terms of quantitative accuracy of activity estimates in the spheres; in the simulation study, the authors evaluated the accuracy and precision of activity estimates from various organs and compared them to results from a previously proposed method. The authors demonstrated excellent agreement between the experimental measurement and Monte Carlo simulation. In the XCAT phantom simulation, the proposed method achieved much better accuracy in the modeling (error in photon counts was -1.1 %) compared to a previously proposed method (errors were more than 20 %); the quantitative accuracy of activity estimates was excellent for all organs (errors were from -1.6 % to 11.9 %) and comparable to previously published results for (131)I using the same collimator. The proposed (90)Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for (131)I. The method may be useful in verifying organ doses for targeted radionuclide therapy using (90)Y.