Current Projects
Development of a non-invasive radiation detector to measure the arterial input function for dynamic positron emission tomography
Dynamic Positron Emission Tomography (dPET) is a nuclear medicine technique where a radioactive tracer molecule is injected into a patient that can be used to image the patient’s biological function. Invaluable information in areas such as oncology, cardiology and neurology can be obtained with dPET, allowing for the use of personalized medicine techniques. Unfortunately, dPET requires measurement of the amount of radiation in a patient’s arteries. This is normally done through invasive arterial blood sampling which is a complicated technique that not all clinics can perform. The goal of this project is to develop a non-invasive radiation detector that will measure this information non-invasively.
2022
Carroll, Liam; Enger, Shirin A.
Second prize at the International Conference on Monte Carlo Techniques for Medical Applications award
2022.
@award{nokey,
title = {Second prize at the International Conference on Monte Carlo Techniques for Medical Applications },
author = {Liam Carroll and Shirin A. Enger},
year = {2022},
date = {2022-04-13},
urldate = {2022-04-13},
journal = {International Conference on Monte Carlo Techniques for Medical Applications },
abstract = {Introduction
Geant4[1] is a Monte Carlo toolkit that provides a flexible platform to design radiation transport simulations. This flexibility requires a high level of complexity when writing new user-codes. New users must undergo extensive training to begin writing useful simulations. The aim of this study was to develop a modular radiation simulation software package called MaRSS based on Geant4 user-code to serve as both an educational tool and as a simulation tool for medical radiation detector simulations.
Materials & Methods
MaRSS builds on Geant4 using Penelope electromagnetic physics models and cross-sections[2]. To give the users possibility to change simulation parameters without changing the source code, MaRSS is equipped with a set of messenger classes which are intercom modules provided by Geant4 to configure applications and provide user interactivity with the code. These messenger classes add additional user commands that can be used to add or remove volumes from the simulation geometry and associated sensitive detectors. The sensitive detectors are objects that Geant4 uses to save simulation results. A number of default volumes and detectors are included in MaRSS that can be used for the design of scintillating fiber-based radiation detectors. MaRSS adds to the existing Geant4 sensitive detector code by creating a new base class called RunSD that is inherited by all sensitive detector objects. RunSD is implemented such that the code needed to initialize a sensitive detector is included in two classes, normally, this code is spread out in several classes. A similar approach is taken with geometrical volumes. To validate MaRSS, range in water of positrons emitted from four radioisotopes commonly used for positron emission tomography was calculated and compared with published work: Fluorine-18 (18F), Carbon-11 (11C), Oxygen-15 (15O) and Gallium-68 (68Ga). A sphere with a radius of 1 m was filed with water. For each radioisotope, 100 million decay events were simulated originating at the center of the simulated water sphere. The resulting positrons were allowed to annihilate. Two energy cuts were simulated, 1 keV and 0.1 keV. Two histograms were created, one using the annihilation locations and another with the energy of the emitted positrons. Results were compared with published positron range values calculated with a PENELOPE[2]-based Monte Carlo software called PeneloPET [3], an analytical expression to estimate the range of positrons described by Cal-Gonzales et.al.[3], a simulation by Lehnert et. al. [4] written using GATE[5] and a separate Monte Carlo software written by Champion and Le Loirec[6] that directly simulates the formation of positronium in water to calculate positron range.
Results
Figure 1 shows the mean calculated positron ranges and calculated positron emission energies compared to literature values. Simulated positron energy means were within 1.8% of literature values. Simulated ranges were within 2% of GATE simulation[4].
References
[1] S. Agostinelli et al., “Geant4—a simulation toolkit,” Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip., vol. 506, no. 3, pp. 250–303, Jul. 2003.
[2] J. Baró, J. Sempau, J. M. Fernández-Varea, and F. Salvat, “PENELOPE: An algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter,” Nucl. Inst. Methods Phys. Res. B, vol. 100, no. 1, pp. 31–46, May 1995.
[3] J. Cal-González et al., “Positron range estimations with PeneloPET,” Phys. Med. Biol., vol. 58, no. 15, pp. 5127–5152, 2013.
[4] W. Lehnert, M.-C. Gregoire, A. Reilhac, and S. R. Meikle, “Analytical positron range modelling in heterogeneous media for PET Monte Carlo simulation,” Phys. Med. Biol., vol. 56, no. 11, p. 3313, May 2011.
[5] D. Strul, G. Santin, D. Lazaro, V. Breton, and C. Morel, “GATE (geant4 application for tomographic emission): a PET/SPECT general-purpose simulation platform,” Nucl. Phys. B - Proc. Suppl., vol. 125, pp. 75–79, Sep. 2003.
[6] C. Champion and C. Le Loirec, “Positron follow-up in liquid water: II. Spatial and energetic study for the most important radioisotopes used in PET,” Phys. Med. Biol., vol. 52, no. 22, pp. 6605–6625, Nov. 2007.
},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Introduction
Geant4[1] is a Monte Carlo toolkit that provides a flexible platform to design radiation transport simulations. This flexibility requires a high level of complexity when writing new user-codes. New users must undergo extensive training to begin writing useful simulations. The aim of this study was to develop a modular radiation simulation software package called MaRSS based on Geant4 user-code to serve as both an educational tool and as a simulation tool for medical radiation detector simulations.
Materials & Methods
MaRSS builds on Geant4 using Penelope electromagnetic physics models and cross-sections[2]. To give the users possibility to change simulation parameters without changing the source code, MaRSS is equipped with a set of messenger classes which are intercom modules provided by Geant4 to configure applications and provide user interactivity with the code. These messenger classes add additional user commands that can be used to add or remove volumes from the simulation geometry and associated sensitive detectors. The sensitive detectors are objects that Geant4 uses to save simulation results. A number of default volumes and detectors are included in MaRSS that can be used for the design of scintillating fiber-based radiation detectors. MaRSS adds to the existing Geant4 sensitive detector code by creating a new base class called RunSD that is inherited by all sensitive detector objects. RunSD is implemented such that the code needed to initialize a sensitive detector is included in two classes, normally, this code is spread out in several classes. A similar approach is taken with geometrical volumes. To validate MaRSS, range in water of positrons emitted from four radioisotopes commonly used for positron emission tomography was calculated and compared with published work: Fluorine-18 (18F), Carbon-11 (11C), Oxygen-15 (15O) and Gallium-68 (68Ga). A sphere with a radius of 1 m was filed with water. For each radioisotope, 100 million decay events were simulated originating at the center of the simulated water sphere. The resulting positrons were allowed to annihilate. Two energy cuts were simulated, 1 keV and 0.1 keV. Two histograms were created, one using the annihilation locations and another with the energy of the emitted positrons. Results were compared with published positron range values calculated with a PENELOPE[2]-based Monte Carlo software called PeneloPET [3], an analytical expression to estimate the range of positrons described by Cal-Gonzales et.al.[3], a simulation by Lehnert et. al. [4] written using GATE[5] and a separate Monte Carlo software written by Champion and Le Loirec[6] that directly simulates the formation of positronium in water to calculate positron range.
Results
Figure 1 shows the mean calculated positron ranges and calculated positron emission energies compared to literature values. Simulated positron energy means were within 1.8% of literature values. Simulated ranges were within 2% of GATE simulation[4].
References
[1] S. Agostinelli et al., “Geant4—a simulation toolkit,” Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip., vol. 506, no. 3, pp. 250–303, Jul. 2003.
[2] J. Baró, J. Sempau, J. M. Fernández-Varea, and F. Salvat, “PENELOPE: An algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter,” Nucl. Inst. Methods Phys. Res. B, vol. 100, no. 1, pp. 31–46, May 1995.
[3] J. Cal-González et al., “Positron range estimations with PeneloPET,” Phys. Med. Biol., vol. 58, no. 15, pp. 5127–5152, 2013.
[4] W. Lehnert, M.-C. Gregoire, A. Reilhac, and S. R. Meikle, “Analytical positron range modelling in heterogeneous media for PET Monte Carlo simulation,” Phys. Med. Biol., vol. 56, no. 11, p. 3313, May 2011.
[5] D. Strul, G. Santin, D. Lazaro, V. Breton, and C. Morel, “GATE (geant4 application for tomographic emission): a PET/SPECT general-purpose simulation platform,” Nucl. Phys. B - Proc. Suppl., vol. 125, pp. 75–79, Sep. 2003.
[6] C. Champion and C. Le Loirec, “Positron follow-up in liquid water: II. Spatial and energetic study for the most important radioisotopes used in PET,” Phys. Med. Biol., vol. 52, no. 22, pp. 6605–6625, Nov. 2007.
Geant4[1] is a Monte Carlo toolkit that provides a flexible platform to design radiation transport simulations. This flexibility requires a high level of complexity when writing new user-codes. New users must undergo extensive training to begin writing useful simulations. The aim of this study was to develop a modular radiation simulation software package called MaRSS based on Geant4 user-code to serve as both an educational tool and as a simulation tool for medical radiation detector simulations.
Materials & Methods
MaRSS builds on Geant4 using Penelope electromagnetic physics models and cross-sections[2]. To give the users possibility to change simulation parameters without changing the source code, MaRSS is equipped with a set of messenger classes which are intercom modules provided by Geant4 to configure applications and provide user interactivity with the code. These messenger classes add additional user commands that can be used to add or remove volumes from the simulation geometry and associated sensitive detectors. The sensitive detectors are objects that Geant4 uses to save simulation results. A number of default volumes and detectors are included in MaRSS that can be used for the design of scintillating fiber-based radiation detectors. MaRSS adds to the existing Geant4 sensitive detector code by creating a new base class called RunSD that is inherited by all sensitive detector objects. RunSD is implemented such that the code needed to initialize a sensitive detector is included in two classes, normally, this code is spread out in several classes. A similar approach is taken with geometrical volumes. To validate MaRSS, range in water of positrons emitted from four radioisotopes commonly used for positron emission tomography was calculated and compared with published work: Fluorine-18 (18F), Carbon-11 (11C), Oxygen-15 (15O) and Gallium-68 (68Ga). A sphere with a radius of 1 m was filed with water. For each radioisotope, 100 million decay events were simulated originating at the center of the simulated water sphere. The resulting positrons were allowed to annihilate. Two energy cuts were simulated, 1 keV and 0.1 keV. Two histograms were created, one using the annihilation locations and another with the energy of the emitted positrons. Results were compared with published positron range values calculated with a PENELOPE[2]-based Monte Carlo software called PeneloPET [3], an analytical expression to estimate the range of positrons described by Cal-Gonzales et.al.[3], a simulation by Lehnert et. al. [4] written using GATE[5] and a separate Monte Carlo software written by Champion and Le Loirec[6] that directly simulates the formation of positronium in water to calculate positron range.
Results
Figure 1 shows the mean calculated positron ranges and calculated positron emission energies compared to literature values. Simulated positron energy means were within 1.8% of literature values. Simulated ranges were within 2% of GATE simulation[4].
References
[1] S. Agostinelli et al., “Geant4—a simulation toolkit,” Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip., vol. 506, no. 3, pp. 250–303, Jul. 2003.
[2] J. Baró, J. Sempau, J. M. Fernández-Varea, and F. Salvat, “PENELOPE: An algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter,” Nucl. Inst. Methods Phys. Res. B, vol. 100, no. 1, pp. 31–46, May 1995.
[3] J. Cal-González et al., “Positron range estimations with PeneloPET,” Phys. Med. Biol., vol. 58, no. 15, pp. 5127–5152, 2013.
[4] W. Lehnert, M.-C. Gregoire, A. Reilhac, and S. R. Meikle, “Analytical positron range modelling in heterogeneous media for PET Monte Carlo simulation,” Phys. Med. Biol., vol. 56, no. 11, p. 3313, May 2011.
[5] D. Strul, G. Santin, D. Lazaro, V. Breton, and C. Morel, “GATE (geant4 application for tomographic emission): a PET/SPECT general-purpose simulation platform,” Nucl. Phys. B - Proc. Suppl., vol. 125, pp. 75–79, Sep. 2003.
[6] C. Champion and C. Le Loirec, “Positron follow-up in liquid water: II. Spatial and energetic study for the most important radioisotopes used in PET,” Phys. Med. Biol., vol. 52, no. 22, pp. 6605–6625, Nov. 2007.
Daoud, Youstina; Carroll, Liam; Enger, Shirin A.
A Radiation detector simulation toolkit for calculating the Arterial Input Function during Dynamic Positron Emission Tomography Inproceedings
In: International Conference on Monte Carlo Techniques for Medical Applications, 2022.
@inproceedings{nokey,
title = {A Radiation detector simulation toolkit for calculating the Arterial Input Function during Dynamic Positron Emission Tomography},
author = {Youstina Daoud and Liam Carroll and Shirin A. Enger},
year = {2022},
date = {2022-04-10},
booktitle = {International Conference on Monte Carlo Techniques for Medical Applications},
abstract = {"Introduction
Dynamic Positron Emission Tomography (dPET) is a functional imaging modality that provides an accurate assessment of patients’ physiological activities and response to treatments such as cancer, cardiac diseases and Alzheimer’s disease. It requires the measurement of the time-course activity concentration of the positron emitting PET radioisotopes in the patient’s arterial plasma, called the Arterial Input Function (AIF). The gold standard measurement of the AIF requires blood samples from the patient during the dPET. In our group, we are developing a non-invasive radiation detector that, placed on a patient’s wrist during the dPET scan, measures the number of positrons and photons escaping the radial artery and calculates the AIF. We have also developed a Modular Radiation Simulation Software for detector simulations called MaRSS that allows the user to run a Geant4-based Monte Carlo simulation, to calculate the AIF. Using the Monte Carlo method, MaRSS simulates a radioactive source decay in the radial artery and scores the amount of radiation escaping the radial artery and reaching the detector placed on the simulated patient’s wrist phantom. The wrist phantom is designed as a cylinder containing 2 holes that simulates the radial artery and vein. The shape and the depth of the radial artery vary between patients and proper knowledge of the distance between the radial artery and the skin, as well as its surface area, is important to accurately design the wrist phantom. Therefore, our aim was to develop a graphical user interface (GUI) allowing the user to import 2D ultrasound scans of a patient’s wrist, provide tools to measure the distance between the radial artery and the skin as well as the radial artery’s surface area and to create the necessary input file to MaRSS. The GUI provides MaRSS with a patient specific and more accurate wrist phantom, providing a patient-specific and more accurate calculation of the AIF without knowledge of C++ or Geant4.
Materials & Methods
The GUI elements were implemented using the multi-platform application and widget toolkit Qt 5 [1]. The C++ library, VTK 8.2.0 [2] was integrated in the GUI, which enables the user to import and manipulate the 2D ultrasound images. The toolkit comprises a measurement tool, a visualization window, a detector tab, a radiation source tab and MaRSS which is its simulation tool. To create an accurate wrist phantom, three 2D – cross secctional ultrasound scans of the patient’s wrist at 2 cm, 4 cm and 6 cm from the wrist crease and 1 longitudinal scan along the radial artery may be acquired and saved in DICOM format. In our case the BK3000 ultrasound system is used. These scans are imported into the GUI by selecting the folder that contains the images. Using the measurement functionalities shown in the top left corner of Figure 1, the surface of the radial artery is measured by drawing an ellipse on the artery’s boundary, then the toolkit measures the surface of the drawn ellipse and displays it in the Measurement window. The artery’s depth is also measured and displayed by drawing a straight line between the artery’s boundary and the skin. Using the left and right arrows, the user can navigate through the selected folder and measure the artery’s surface and depth on the other scans. The top right corner of the GUI shown in Figure 1, illustrates a Detector tab and a Source tab. The Detector tab allows the user to import a detector in STL format and place it on the ultrasound scan to simulate different setups of the detector, this functionality is still under development and
International Conference on Monte Carlo Techniques for Medical Applications, 2022
is optional. The Source tab allows the user to add the radioactive source used during the dPET by entering its mass number and its atomic number. After completing the 3 mandatory steps : import of the scan, measurement of different parameters extracted from the scan and choice of the radioactive source, the user can run the simulation by clicking on Run Simulation in the Simulation menu. The toolkit runs the MaRSS and creates the wrist phantom using the artery’s surface and depth measured by the user, then starts the decay of the chosen source placed randomly inside the artery.
Results
This toolkit allows the user to import 2D ultrasound scans and measure the radial artery’s surface and depth along the wrist, choose the radioactive source from the Source drop-down menu and specify the detector position. An input file to the MaRSS is thus created providing the required information to simulate the wrist phantom, the source and the detector’s position in MaRSS. The Run Simulation tab displays the output of the simulation in the GUI making it the only used tool for setting up the simulation and viewing the results.
Discussion & Conclusions
This toolkit enables the user to run a Geant4 Monte Carlo based simulation for detector development applications in 3 easy steps, not requiring any programming knowledge.
References
[1] Blanchette J, Summerfield M. C++ GUI programming with Qt 4: Prentice Hall Professional; 2006.
[2] Schroeder WJ, Avila LS, Hoffman W. Visualizing with VTK: a tutorial. IEEE Computer graphics and applications. 2000;20(5):20-7.
Acknowledgements This research was undertaken,in part, thanks to funding from the Canada Research Chairs Program (grant # 252135) as well as CHRP (NSERC+CIHR grant 170620).
"},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
"Introduction
Dynamic Positron Emission Tomography (dPET) is a functional imaging modality that provides an accurate assessment of patients’ physiological activities and response to treatments such as cancer, cardiac diseases and Alzheimer’s disease. It requires the measurement of the time-course activity concentration of the positron emitting PET radioisotopes in the patient’s arterial plasma, called the Arterial Input Function (AIF). The gold standard measurement of the AIF requires blood samples from the patient during the dPET. In our group, we are developing a non-invasive radiation detector that, placed on a patient’s wrist during the dPET scan, measures the number of positrons and photons escaping the radial artery and calculates the AIF. We have also developed a Modular Radiation Simulation Software for detector simulations called MaRSS that allows the user to run a Geant4-based Monte Carlo simulation, to calculate the AIF. Using the Monte Carlo method, MaRSS simulates a radioactive source decay in the radial artery and scores the amount of radiation escaping the radial artery and reaching the detector placed on the simulated patient’s wrist phantom. The wrist phantom is designed as a cylinder containing 2 holes that simulates the radial artery and vein. The shape and the depth of the radial artery vary between patients and proper knowledge of the distance between the radial artery and the skin, as well as its surface area, is important to accurately design the wrist phantom. Therefore, our aim was to develop a graphical user interface (GUI) allowing the user to import 2D ultrasound scans of a patient’s wrist, provide tools to measure the distance between the radial artery and the skin as well as the radial artery’s surface area and to create the necessary input file to MaRSS. The GUI provides MaRSS with a patient specific and more accurate wrist phantom, providing a patient-specific and more accurate calculation of the AIF without knowledge of C++ or Geant4.
Materials & Methods
The GUI elements were implemented using the multi-platform application and widget toolkit Qt 5 [1]. The C++ library, VTK 8.2.0 [2] was integrated in the GUI, which enables the user to import and manipulate the 2D ultrasound images. The toolkit comprises a measurement tool, a visualization window, a detector tab, a radiation source tab and MaRSS which is its simulation tool. To create an accurate wrist phantom, three 2D – cross secctional ultrasound scans of the patient’s wrist at 2 cm, 4 cm and 6 cm from the wrist crease and 1 longitudinal scan along the radial artery may be acquired and saved in DICOM format. In our case the BK3000 ultrasound system is used. These scans are imported into the GUI by selecting the folder that contains the images. Using the measurement functionalities shown in the top left corner of Figure 1, the surface of the radial artery is measured by drawing an ellipse on the artery’s boundary, then the toolkit measures the surface of the drawn ellipse and displays it in the Measurement window. The artery’s depth is also measured and displayed by drawing a straight line between the artery’s boundary and the skin. Using the left and right arrows, the user can navigate through the selected folder and measure the artery’s surface and depth on the other scans. The top right corner of the GUI shown in Figure 1, illustrates a Detector tab and a Source tab. The Detector tab allows the user to import a detector in STL format and place it on the ultrasound scan to simulate different setups of the detector, this functionality is still under development and
International Conference on Monte Carlo Techniques for Medical Applications, 2022
is optional. The Source tab allows the user to add the radioactive source used during the dPET by entering its mass number and its atomic number. After completing the 3 mandatory steps : import of the scan, measurement of different parameters extracted from the scan and choice of the radioactive source, the user can run the simulation by clicking on Run Simulation in the Simulation menu. The toolkit runs the MaRSS and creates the wrist phantom using the artery’s surface and depth measured by the user, then starts the decay of the chosen source placed randomly inside the artery.
Results
This toolkit allows the user to import 2D ultrasound scans and measure the radial artery’s surface and depth along the wrist, choose the radioactive source from the Source drop-down menu and specify the detector position. An input file to the MaRSS is thus created providing the required information to simulate the wrist phantom, the source and the detector’s position in MaRSS. The Run Simulation tab displays the output of the simulation in the GUI making it the only used tool for setting up the simulation and viewing the results.
Discussion & Conclusions
This toolkit enables the user to run a Geant4 Monte Carlo based simulation for detector development applications in 3 easy steps, not requiring any programming knowledge.
References
[1] Blanchette J, Summerfield M. C++ GUI programming with Qt 4: Prentice Hall Professional; 2006.
[2] Schroeder WJ, Avila LS, Hoffman W. Visualizing with VTK: a tutorial. IEEE Computer graphics and applications. 2000;20(5):20-7.
Acknowledgements This research was undertaken,in part, thanks to funding from the Canada Research Chairs Program (grant # 252135) as well as CHRP (NSERC+CIHR grant 170620).
"
Dynamic Positron Emission Tomography (dPET) is a functional imaging modality that provides an accurate assessment of patients’ physiological activities and response to treatments such as cancer, cardiac diseases and Alzheimer’s disease. It requires the measurement of the time-course activity concentration of the positron emitting PET radioisotopes in the patient’s arterial plasma, called the Arterial Input Function (AIF). The gold standard measurement of the AIF requires blood samples from the patient during the dPET. In our group, we are developing a non-invasive radiation detector that, placed on a patient’s wrist during the dPET scan, measures the number of positrons and photons escaping the radial artery and calculates the AIF. We have also developed a Modular Radiation Simulation Software for detector simulations called MaRSS that allows the user to run a Geant4-based Monte Carlo simulation, to calculate the AIF. Using the Monte Carlo method, MaRSS simulates a radioactive source decay in the radial artery and scores the amount of radiation escaping the radial artery and reaching the detector placed on the simulated patient’s wrist phantom. The wrist phantom is designed as a cylinder containing 2 holes that simulates the radial artery and vein. The shape and the depth of the radial artery vary between patients and proper knowledge of the distance between the radial artery and the skin, as well as its surface area, is important to accurately design the wrist phantom. Therefore, our aim was to develop a graphical user interface (GUI) allowing the user to import 2D ultrasound scans of a patient’s wrist, provide tools to measure the distance between the radial artery and the skin as well as the radial artery’s surface area and to create the necessary input file to MaRSS. The GUI provides MaRSS with a patient specific and more accurate wrist phantom, providing a patient-specific and more accurate calculation of the AIF without knowledge of C++ or Geant4.
Materials & Methods
The GUI elements were implemented using the multi-platform application and widget toolkit Qt 5 [1]. The C++ library, VTK 8.2.0 [2] was integrated in the GUI, which enables the user to import and manipulate the 2D ultrasound images. The toolkit comprises a measurement tool, a visualization window, a detector tab, a radiation source tab and MaRSS which is its simulation tool. To create an accurate wrist phantom, three 2D – cross secctional ultrasound scans of the patient’s wrist at 2 cm, 4 cm and 6 cm from the wrist crease and 1 longitudinal scan along the radial artery may be acquired and saved in DICOM format. In our case the BK3000 ultrasound system is used. These scans are imported into the GUI by selecting the folder that contains the images. Using the measurement functionalities shown in the top left corner of Figure 1, the surface of the radial artery is measured by drawing an ellipse on the artery’s boundary, then the toolkit measures the surface of the drawn ellipse and displays it in the Measurement window. The artery’s depth is also measured and displayed by drawing a straight line between the artery’s boundary and the skin. Using the left and right arrows, the user can navigate through the selected folder and measure the artery’s surface and depth on the other scans. The top right corner of the GUI shown in Figure 1, illustrates a Detector tab and a Source tab. The Detector tab allows the user to import a detector in STL format and place it on the ultrasound scan to simulate different setups of the detector, this functionality is still under development and
International Conference on Monte Carlo Techniques for Medical Applications, 2022
is optional. The Source tab allows the user to add the radioactive source used during the dPET by entering its mass number and its atomic number. After completing the 3 mandatory steps : import of the scan, measurement of different parameters extracted from the scan and choice of the radioactive source, the user can run the simulation by clicking on Run Simulation in the Simulation menu. The toolkit runs the MaRSS and creates the wrist phantom using the artery’s surface and depth measured by the user, then starts the decay of the chosen source placed randomly inside the artery.
Results
This toolkit allows the user to import 2D ultrasound scans and measure the radial artery’s surface and depth along the wrist, choose the radioactive source from the Source drop-down menu and specify the detector position. An input file to the MaRSS is thus created providing the required information to simulate the wrist phantom, the source and the detector’s position in MaRSS. The Run Simulation tab displays the output of the simulation in the GUI making it the only used tool for setting up the simulation and viewing the results.
Discussion & Conclusions
This toolkit enables the user to run a Geant4 Monte Carlo based simulation for detector development applications in 3 easy steps, not requiring any programming knowledge.
References
[1] Blanchette J, Summerfield M. C++ GUI programming with Qt 4: Prentice Hall Professional; 2006.
[2] Schroeder WJ, Avila LS, Hoffman W. Visualizing with VTK: a tutorial. IEEE Computer graphics and applications. 2000;20(5):20-7.
Acknowledgements This research was undertaken,in part, thanks to funding from the Canada Research Chairs Program (grant # 252135) as well as CHRP (NSERC+CIHR grant 170620).
"
Martinez, Victor Daniel Diaz; Carroll, Liam; Enger, Shirin A.
Monte Carlo Simulation of the 224Ra Decay Chain and the Diffusion of 220Rn for Diffusing Alpha-Emitters Radiotherapy Inproceedings
In: MEDICAL PHYSICS, pp. E828–E828, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
@inproceedings{martinez2022monte,
title = {Monte Carlo Simulation of the 224Ra Decay Chain and the Diffusion of 220Rn for Diffusing Alpha-Emitters Radiotherapy},
author = {Victor Daniel Diaz Martinez and Liam Carroll and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E828--E828},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Carroll, Liam; Enger, Shirin A.
Monte Carlo Simulations of a Non-Invasive Positron Detector to Measure the Arterial Input Function for Dynamic PET Inproceedings
In: Journal of Physics: Conference Series, pp. 012005, IOP Publishing 2022.
@inproceedings{carroll2022monte,
title = {Monte Carlo Simulations of a Non-Invasive Positron Detector to Measure the Arterial Input Function for Dynamic PET},
author = {Liam Carroll and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {Journal of Physics: Conference Series},
volume = {2167},
number = {1},
pages = {012005},
organization = {IOP Publishing},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Carroll, Liam; Enger, Shirin A.
Non-invasive measurement of the arterial input function for dynamic positron emission tomography: Simulation of clinical workflow Inproceedings
In: MEDICAL PHYSICS, pp. 5643–5644, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
@inproceedings{carroll2022non,
title = {Non-invasive measurement of the arterial input function for dynamic positron emission tomography: Simulation of clinical workflow},
author = {Liam Carroll and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {8},
pages = {5643--5644},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Daoud, Youstina; Carroll, Liam; Enger, Shrin A.
PO-1617 Mapping of the human wrist to develop a non-invasive radiation detector for Dynamic PET application Journal Article
In: Radiotherapy and Oncology, vol. 170, pp. S1405–S1406, 2022.
@article{daoud2022po,
title = {PO-1617 Mapping of the human wrist to develop a non-invasive radiation detector for Dynamic PET application},
author = {Youstina Daoud and Liam Carroll and Shrin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Radiotherapy and Oncology},
volume = {170},
pages = {S1405--S1406},
publisher = {Elsevier},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Carroll, Liam; Enger, Shirin A
Simulation of a novel, non-invasive radiation detector to measure the arterial input function for dynamic PET Journal Article
In: Medical Physics, 2022.
@article{carroll2022simulation,
title = {Simulation of a novel, non-invasive radiation detector to measure the arterial input function for dynamic PET},
author = {Liam Carroll and Shirin A Enger},
year = {2022},
date = {2022-01-01},
journal = {Medical Physics},
publisher = {Wiley Online Library},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
2021
Thibodeau-Antonacci, Alana; Jafarzadeh, Hossein; Carroll, Liam; Weishaupt, Luca L.
Mitacs Globalink Research Award award
2021.
@award{Thibodeau-Antonacci2021c,
title = {Mitacs Globalink Research Award},
author = {Alana Thibodeau-Antonacci and Hossein Jafarzadeh and Liam Carroll and Luca L. Weishaupt},
url = {https://www.mitacs.ca/en/programs/globalink/globalink-research-award},
year = {2021},
date = {2021-07-01},
urldate = {2021-07-01},
organization = {MITACS},
abstract = {The Mitacs Globalink Research Award (GRA) supports research collaborations between Canada and select partner organizations and eligible countries and regions. It was awarded to Alana Thibodeau-Antonacci, Hossein Jafarzadeh, Liam Carroll and Luca L. Weishaupt.
Under the joint supervision of a home and host professor, successful senior undergraduate students, graduate students, as well as postdoctoral fellows will receive a $6,000 research award to conduct a 12- to 24-week research project in the other country. Awards are offered in partnership with Mitacs’s Canadian academic partners (and, in some cases, with Mitacs’s international partners) and are subject to available funding. },
howpublished = {Mitacs},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
The Mitacs Globalink Research Award (GRA) supports research collaborations between Canada and select partner organizations and eligible countries and regions. It was awarded to Alana Thibodeau-Antonacci, Hossein Jafarzadeh, Liam Carroll and Luca L. Weishaupt.
Under the joint supervision of a home and host professor, successful senior undergraduate students, graduate students, as well as postdoctoral fellows will receive a $6,000 research award to conduct a 12- to 24-week research project in the other country. Awards are offered in partnership with Mitacs’s Canadian academic partners (and, in some cases, with Mitacs’s international partners) and are subject to available funding.
Under the joint supervision of a home and host professor, successful senior undergraduate students, graduate students, as well as postdoctoral fellows will receive a $6,000 research award to conduct a 12- to 24-week research project in the other country. Awards are offered in partnership with Mitacs’s Canadian academic partners (and, in some cases, with Mitacs’s international partners) and are subject to available funding.
2020
Carroll, Liam; Croteau, Etienne; Kertzscher, Gustavo; Sarrhini, Otman; Turgeon, Vincent; Lecomte, Roger; Enger, Shirin A.
In: Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB), vol. 76, pp. 92–99, 2020, ISSN: 1724-191X.
@article{carroll_cross-validation_2020,
title = {Cross-validation of a non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic positron emission tomography},
author = {Liam Carroll and Etienne Croteau and Gustavo Kertzscher and Otman Sarrhini and Vincent Turgeon and Roger Lecomte and Shirin A. Enger},
doi = {10.1016/j.ejmp.2020.06.009},
issn = {1724-191X},
year = {2020},
date = {2020-08-01},
journal = {Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB)},
volume = {76},
pages = {92--99},
abstract = {Kinetic modeling of positron emission tomography (PET) data can assess index rate of uptake, metabolism and predict disease progression more accurately than conventional static PET. However, it requires knowledge of the time-course of the arterial blood radioactivity concentration, called the arterial input function (AIF). The gold standard to acquire the AIF is by invasive means. The purpose of this study was to validate a previously developed dual readout scintillating fiber-based non-invasive positron detector, hereinafter called non-invasive detector (NID), developed to determine the AIF for dynamic PET measured from the human radial artery. The NID consisted of a 3 m long plastic scintillating fiber with each end coupled to a 5 m long transmission fiber followed by a silicon photomultiplier. The scintillating fiber was enclosed inside the grooves of a plastic cylindrical shell. Two sets of experiments were performed to test the NID against a previously validated microfluidic positron detector. A closed-loop microfluidic system combined with a wrist phantom was used. During the first experiment, the three PET radioisotopes 18F, 11C and 68Ga were tested. After optimizing the detector, a second series of tests were performed using only 18F and 11C. The maximum pulse amplitude to electronic noise ratio was 52 obtained with 11C. Linear regressions showed a linear relation between the two detectors. These preliminary results show that the NID can accurately detect positrons from a patient's wrist and has the potential to non-invasively measure the AIF during a dynamic PET scan. The accuracy of these measurements needs to be determined.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Kinetic modeling of positron emission tomography (PET) data can assess index rate of uptake, metabolism and predict disease progression more accurately than conventional static PET. However, it requires knowledge of the time-course of the arterial blood radioactivity concentration, called the arterial input function (AIF). The gold standard to acquire the AIF is by invasive means. The purpose of this study was to validate a previously developed dual readout scintillating fiber-based non-invasive positron detector, hereinafter called non-invasive detector (NID), developed to determine the AIF for dynamic PET measured from the human radial artery. The NID consisted of a 3 m long plastic scintillating fiber with each end coupled to a 5 m long transmission fiber followed by a silicon photomultiplier. The scintillating fiber was enclosed inside the grooves of a plastic cylindrical shell. Two sets of experiments were performed to test the NID against a previously validated microfluidic positron detector. A closed-loop microfluidic system combined with a wrist phantom was used. During the first experiment, the three PET radioisotopes 18F, 11C and 68Ga were tested. After optimizing the detector, a second series of tests were performed using only 18F and 11C. The maximum pulse amplitude to electronic noise ratio was 52 obtained with 11C. Linear regressions showed a linear relation between the two detectors. These preliminary results show that the NID can accurately detect positrons from a patient's wrist and has the potential to non-invasively measure the AIF during a dynamic PET scan. The accuracy of these measurements needs to be determined.
Carroll, Liam
First prize for best presentation at the congrès annuel de l'AQPMC award
2020.
@award{Carroll2020b,
title = {First prize for best presentation at the congrès annuel de l'AQPMC},
author = {Liam Carroll},
year = {2020},
date = {2020-07-01},
urldate = {2020-07-01},
organization = {Congrès Annuel de l'AQPMC},
abstract = {Nous tenons particulièrement à féliciter Liam Carroll, étudiant à McGill Medical Physics qui a remporté le prix de la meilleure présentation.},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Nous tenons particulièrement à féliciter Liam Carroll, étudiant à McGill Medical Physics qui a remporté le prix de la meilleure présentation.
Carroll, Liam
Excellence Award from Department of Biological and Biomedical Engineering award
2020.
@award{Carroll2020,
title = {Excellence Award from Department of Biological and Biomedical Engineering},
author = {Liam Carroll},
year = {2020},
date = {2020-01-01},
urldate = {2020-01-01},
organization = {McGill University},
school = {Bioengineering and Biomedical Engineering Program},
abstract = {The award recognizes excellent performance in the PhD Program of the BBME/BME graduate program. The selection is based on a combination of overall excellence, performance in the program, grades, progress in project and publications. },
howpublished = {Department of Biological and Biomedical Engineering},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
The award recognizes excellent performance in the PhD Program of the BBME/BME graduate program. The selection is based on a combination of overall excellence, performance in the program, grades, progress in project and publications.
2019
Turgeon, Vincent; Kertzscher, Gustavo; Carroll, Liam; Hopewell, Robert; Massarweh, Gassan; Enger, Shirin A.
Characterization of scintillating fibers for use as positron detector in positron emission tomography Journal Article
In: Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB), vol. 65, pp. 114–120, 2019, ISSN: 1724-191X.
@article{turgeon_characterization_2019,
title = {Characterization of scintillating fibers for use as positron detector in positron emission tomography},
author = {Vincent Turgeon and Gustavo Kertzscher and Liam Carroll and Robert Hopewell and Gassan Massarweh and Shirin A. Enger},
doi = {10.1016/j.ejmp.2019.08.009},
issn = {1724-191X},
year = {2019},
date = {2019-09-01},
journal = {Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB)},
volume = {65},
pages = {114--120},
abstract = {PURPOSE: Manual and automatic blood sampling at different time intervals is considered the gold standard to determine the arterial input function (AIF) in dynamic positron emission tomography (PET). However, blood sampling is characterized by poor time resolution and is an invasive procedure. The aim of this study was to characterize the scintillating fibers used to develop a non-invasive positron detector.
METHODS: The detector consists of a scintillating fiber coupled at each end to transmission fiber-optic cables that are connected to photo multiplier tubes in a dual readout setup. The detector is designed to be wrapped around the wrist of the patient undergoing dynamic PET. The attenuation length and bending losses were measured with excitation from gamma radiation (137Cs) and ultraviolet (UV) light. The response to positron-emitting radio-tracers was evaluated with 18F and 11C.
RESULTS: The attenuation length for a 3.0 m and 1.5 m long scintillating fiber both coincides with the attenuation length given by the manufacturer when excited with the 137Cs source, but not with the UV source due to the differences in scintillation mechanisms. The bending losses are smaller than the measurement uncertainty for the 137Cs source irradiation, and increase when the bending radius decrease for the UV source irradiation. The signal-to-noise ratio for 18F and 11C solutions are 1.98 and 22.54 respectively. The measured decay constant of 11C agrees with its characteristic value.
CONCLUSION: The performed measurements in the dual readout configuration suggest that scintillating fibers may be suitable to determine the AIF non-invasively.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
PURPOSE: Manual and automatic blood sampling at different time intervals is considered the gold standard to determine the arterial input function (AIF) in dynamic positron emission tomography (PET). However, blood sampling is characterized by poor time resolution and is an invasive procedure. The aim of this study was to characterize the scintillating fibers used to develop a non-invasive positron detector.
METHODS: The detector consists of a scintillating fiber coupled at each end to transmission fiber-optic cables that are connected to photo multiplier tubes in a dual readout setup. The detector is designed to be wrapped around the wrist of the patient undergoing dynamic PET. The attenuation length and bending losses were measured with excitation from gamma radiation (137Cs) and ultraviolet (UV) light. The response to positron-emitting radio-tracers was evaluated with 18F and 11C.
RESULTS: The attenuation length for a 3.0 m and 1.5 m long scintillating fiber both coincides with the attenuation length given by the manufacturer when excited with the 137Cs source, but not with the UV source due to the differences in scintillation mechanisms. The bending losses are smaller than the measurement uncertainty for the 137Cs source irradiation, and increase when the bending radius decrease for the UV source irradiation. The signal-to-noise ratio for 18F and 11C solutions are 1.98 and 22.54 respectively. The measured decay constant of 11C agrees with its characteristic value.
CONCLUSION: The performed measurements in the dual readout configuration suggest that scintillating fibers may be suitable to determine the AIF non-invasively.
METHODS: The detector consists of a scintillating fiber coupled at each end to transmission fiber-optic cables that are connected to photo multiplier tubes in a dual readout setup. The detector is designed to be wrapped around the wrist of the patient undergoing dynamic PET. The attenuation length and bending losses were measured with excitation from gamma radiation (137Cs) and ultraviolet (UV) light. The response to positron-emitting radio-tracers was evaluated with 18F and 11C.
RESULTS: The attenuation length for a 3.0 m and 1.5 m long scintillating fiber both coincides with the attenuation length given by the manufacturer when excited with the 137Cs source, but not with the UV source due to the differences in scintillation mechanisms. The bending losses are smaller than the measurement uncertainty for the 137Cs source irradiation, and increase when the bending radius decrease for the UV source irradiation. The signal-to-noise ratio for 18F and 11C solutions are 1.98 and 22.54 respectively. The measured decay constant of 11C agrees with its characteristic value.
CONCLUSION: The performed measurements in the dual readout configuration suggest that scintillating fibers may be suitable to determine the AIF non-invasively.
Carroll, Liam; Croteau, Etienne; Kertzscher, Gustavo; Sarrhini, Otman; Lecomte, Roger; Enger, Shirin A.
Validation of a Non-Invasive Positron Detector for Use with Dynamic PET Under Clinical Conditions Presentation
AAPM 61st Annual Meeting, 14.07.2019.
@misc{nokey,
title = {Validation of a Non-Invasive Positron Detector for Use with Dynamic PET Under Clinical Conditions},
author = {Liam Carroll and Etienne Croteau and Gustavo Kertzscher and Otman Sarrhini and Roger Lecomte and Shirin A. Enger},
year = {2019},
date = {2019-07-14},
howpublished = {AAPM 61st Annual Meeting},
keywords = {},
pubstate = {published},
tppubtype = {presentation}
}
Carroll, Liam; Croteau, Etienne; Kertzscher, Gustavo; Sarrhini, Otman; Lecomte, Roger; Enger, Shirin A.
Cross-validation of a non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic PET imaging Presentation
Society of Nuclear Medicine and Molecular Imaging (SNMMI), 24.06.2019.
@misc{Carroll2019c,
title = {Cross-validation of a non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic PET imaging},
author = {Liam Carroll and Etienne Croteau and Gustavo Kertzscher and Otman Sarrhini and Roger Lecomte and Shirin A. Enger},
year = {2019},
date = {2019-06-24},
howpublished = {Society of Nuclear Medicine and Molecular Imaging (SNMMI)},
keywords = {},
pubstate = {published},
tppubtype = {presentation}
}
Carroll, Liam; Kertzscher, Gustavo; Enger, Shirin A.
Winners of the Marika Zelenka Roy Innovation Prize award
2019.
@award{Carroll2019b,
title = {Winners of the Marika Zelenka Roy Innovation Prize},
author = {Liam Carroll and Gustavo Kertzscher and Shirin A. Enger},
year = {2019},
date = {2019-01-01},
urldate = {2019-01-01},
organization = {McGill University},
abstract = {(50 k CAD) at McGill University Clinical Innovation Competition for the detector BetaSense (A non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic positron emission tomography imaging). },
howpublished = {McGill University Clinical Innovation Competition},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
(50 k CAD) at McGill University Clinical Innovation Competition for the detector BetaSense (A non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic positron emission tomography imaging).
Carroll, Liam; Kertzscher, Gustavo; Enger, Shirin A.
Winners of the Centech propulsion program award
2019.
@award{Carroll2019,
title = {Winners of the Centech propulsion program},
author = {Liam Carroll and Gustavo Kertzscher and Shirin A. Enger},
year = {2019},
date = {2019-01-01},
urldate = {2019-01-01},
organization = {Centech},
abstract = {Two years’ business incubator, (50 k CAD) + access to consulting firms to build a commercial version of BetaSense that can be used in a clinical environment for patient measurements. },
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Two years’ business incubator, (50 k CAD) + access to consulting firms to build a commercial version of BetaSense that can be used in a clinical environment for patient measurements.