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

Xing Rong, Yong Du, Michael Ljungberg, Erwann Rault, Stefaan Vandenberghe, Eric Frey

Research output: Contribution to journalArticle

Abstract

Purpose: Yttrium-90 (90Y) 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., 111In) or imaging agents (e.g., 99mTc MAA) are typically used for treatment planning. It would, however, be useful to image 90Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the 90Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of 90Y 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 90Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. Methods: 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 90Y 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 90Y 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. Results: 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 131I using the same collimator. Conclusions: The proposed 90Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for 131I. The method may be useful in verifying organ doses for targeted radionuclide therapy using 90Y.

Original languageEnglish (US)
Pages (from-to)2346-2358
Number of pages13
JournalMedical Physics
Volume39
Issue number5
DOIs
StatePublished - May 2012

Fingerprint

Single-Photon Emission-Computed Tomography
Photons
Radioisotopes
Yttrium
Monte Carlo Method
Therapeutics

Keywords

  • model-based compensation
  • Monte Carlo simulation
  • quantitative SPECT
  • targeted radionuclide therapy
  • Yttrium-90 bremsstrahlung

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

Development and evaluation of an improved quantitative 90Y bremsstrahlung SPECT method. / Rong, Xing; Du, Yong; Ljungberg, Michael; Rault, Erwann; Vandenberghe, Stefaan; Frey, Eric.

In: Medical Physics, Vol. 39, No. 5, 05.2012, p. 2346-2358.

Research output: Contribution to journalArticle

Rong, Xing ; Du, Yong ; Ljungberg, Michael ; Rault, Erwann ; Vandenberghe, Stefaan ; Frey, Eric. / Development and evaluation of an improved quantitative 90Y bremsstrahlung SPECT method. In: Medical Physics. 2012 ; Vol. 39, No. 5. pp. 2346-2358.
@article{897e7bcb84404c70a5f5bb3ce4f830f7,
title = "Development and evaluation of an improved quantitative 90Y bremsstrahlung SPECT method",
abstract = "Purpose: Yttrium-90 (90Y) 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., 111In) or imaging agents (e.g., 99mTc MAA) are typically used for treatment planning. It would, however, be useful to image 90Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the 90Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of 90Y 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 90Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. Methods: 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 90Y 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 90Y 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. Results: 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 131I using the same collimator. Conclusions: The proposed 90Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for 131I. The method may be useful in verifying organ doses for targeted radionuclide therapy using 90Y.",
keywords = "model-based compensation, Monte Carlo simulation, quantitative SPECT, targeted radionuclide therapy, Yttrium-90 bremsstrahlung",
author = "Xing Rong and Yong Du and Michael Ljungberg and Erwann Rault and Stefaan Vandenberghe and Eric Frey",
year = "2012",
month = "5",
doi = "10.1118/1.3700174",
language = "English (US)",
volume = "39",
pages = "2346--2358",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "5",

}

TY - JOUR

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

AU - Rong, Xing

AU - Du, Yong

AU - Ljungberg, Michael

AU - Rault, Erwann

AU - Vandenberghe, Stefaan

AU - Frey, Eric

PY - 2012/5

Y1 - 2012/5

N2 - Purpose: Yttrium-90 (90Y) 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., 111In) or imaging agents (e.g., 99mTc MAA) are typically used for treatment planning. It would, however, be useful to image 90Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the 90Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of 90Y 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 90Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. Methods: 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 90Y 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 90Y 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. Results: 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 131I using the same collimator. Conclusions: The proposed 90Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for 131I. The method may be useful in verifying organ doses for targeted radionuclide therapy using 90Y.

AB - Purpose: Yttrium-90 (90Y) 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., 111In) or imaging agents (e.g., 99mTc MAA) are typically used for treatment planning. It would, however, be useful to image 90Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the 90Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of 90Y 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 90Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. Methods: 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 90Y 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 90Y 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. Results: 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 131I using the same collimator. Conclusions: The proposed 90Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for 131I. The method may be useful in verifying organ doses for targeted radionuclide therapy using 90Y.

KW - model-based compensation

KW - Monte Carlo simulation

KW - quantitative SPECT

KW - targeted radionuclide therapy

KW - Yttrium-90 bremsstrahlung

UR - http://www.scopus.com/inward/record.url?scp=84861638196&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84861638196&partnerID=8YFLogxK

U2 - 10.1118/1.3700174

DO - 10.1118/1.3700174

M3 - Article

VL - 39

SP - 2346

EP - 2358

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 5

ER -