2023
Rahbaran, Maryam; Kalinowski, Jonathan; Tsui, James Man Git; DeCunha, Joseph; Croce, Kevin; Bergmark, Brian; Devlin, Philip; Enger, Shirin A.
Development Of a Novel Dosimetry Software for Patient-specific Intravascular Brachytherapy Treatment Planning on Optical Coherence Tomography Images Presentation
22.06.2023, (2023 American Brachytherapy Society (ABS) Annual Meeting, Vancouver, Canada).
BibTeX | Tags:
@misc{nokey,
title = {Development Of a Novel Dosimetry Software for Patient-specific Intravascular Brachytherapy Treatment Planning on Optical Coherence Tomography Images},
author = {Maryam Rahbaran and Jonathan Kalinowski and James Man Git Tsui and Joseph DeCunha and Kevin Croce and Brian Bergmark and Philip Devlin and Shirin A. Enger},
year = {2023},
date = {2023-06-22},
note = {2023 American Brachytherapy Society (ABS) Annual Meeting, Vancouver, Canada},
keywords = {},
pubstate = {published},
tppubtype = {presentation}
}
Mégrourèche, Julien; Bekerat, Hamed; Bian, Jingyi; Bui, Alaina; Sankey, Jack; Childress, Lilian; Enger, Shirin A.
Development of a hydrated electron dosimeter for radiotherapy applications: A proof of concept Journal Article
In: Medical Physics, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey,
title = {Development of a hydrated electron dosimeter for radiotherapy applications: A proof of concept},
author = {Julien Mégrourèche and Hamed Bekerat and Jingyi Bian and Alaina Bui and Jack Sankey and Lilian Childress and Shirin A. Enger},
doi = {10.1002/mp.16555},
year = {2023},
date = {2023-06-19},
urldate = {2023-06-19},
journal = {Medical Physics},
abstract = {Background: Hydrated electrons, which are short-lived products of radiolysis in water, increase the optical absorption of water, providing a pathway toward near-tissue-equivalent clinical radiation dosimeters. This has been demonstrated in high-dose-per-pulse radiochemistry research, but, owing to the weak absorption signal, its application in existing low-dose-per-pulse radiotherapy provided by clinical linear accelerators (linacs) has yet to be investigated.
Purpose: The aims of this study were to measure the optical absorption associated with hydrated electrons produced by clinical linacs and to assess the suitability of the technique for radiotherapy (⩽ 1 cGy per pulse) applications.
Methods: 40 mW of 660-nm laser light was sent five passes through deionized water contained in a 10
×
4
×
2 cm3 glass-walled cavity by using four broadband dielectric mirrors, two on each side of the cavity. The light was collected with a biased silicon photodetector. The water cavity was then irradiated by a Varian TrueBeam linac with both photon (10 MV FFF, 6 MV FFF, 6 MV) and electron beams (6 MeV) while monitoring the transmitted laser power for absorption transients. Radiochromic EBT3 film measurements were also performed for comparison.
Results: Examination of the absorbance profiles showed clear absorption changes in the water when radiation pulses were delivered. Both the amplitude and the decay time of the signal appeared consistent with the absorbed dose and the characteristics of the hydrated electrons. By using literature value for the hydrated electron radiation chemical yield (3.0±0.3), we inferred doses of 2.1±0.2 mGy (10 MV FFF), 1.3±0.1 mGy (6 MV FFF), 0.45±0.06 mGy (6 MV) for photons, and 0.47±0.05 mGy (6 MeV) for electrons, which differed from EBT3 film measurements by 0.6%, 0.8%, 10%, and 15.7%, respectively. The half-life of the hydrated electrons in the solution was ∼ 24
μs.
Conclusions: By measuring 660-nm laser light transmitted through a cm-scale, multi-pass water cavity, we observed absorption transients consistent with hydrated electrons generated by clinical linac radiation. The agreement between our inferred dose and EBT3 film measurements suggests this proof-of-concept system represents a viable pathway toward tissue-equivalent dosimeters for clinical radiotherapy applications.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Background: Hydrated electrons, which are short-lived products of radiolysis in water, increase the optical absorption of water, providing a pathway toward near-tissue-equivalent clinical radiation dosimeters. This has been demonstrated in high-dose-per-pulse radiochemistry research, but, owing to the weak absorption signal, its application in existing low-dose-per-pulse radiotherapy provided by clinical linear accelerators (linacs) has yet to be investigated.
Purpose: The aims of this study were to measure the optical absorption associated with hydrated electrons produced by clinical linacs and to assess the suitability of the technique for radiotherapy (⩽ 1 cGy per pulse) applications.
Methods: 40 mW of 660-nm laser light was sent five passes through deionized water contained in a 10
×
4
×
2 cm3 glass-walled cavity by using four broadband dielectric mirrors, two on each side of the cavity. The light was collected with a biased silicon photodetector. The water cavity was then irradiated by a Varian TrueBeam linac with both photon (10 MV FFF, 6 MV FFF, 6 MV) and electron beams (6 MeV) while monitoring the transmitted laser power for absorption transients. Radiochromic EBT3 film measurements were also performed for comparison.
Results: Examination of the absorbance profiles showed clear absorption changes in the water when radiation pulses were delivered. Both the amplitude and the decay time of the signal appeared consistent with the absorbed dose and the characteristics of the hydrated electrons. By using literature value for the hydrated electron radiation chemical yield (3.0±0.3), we inferred doses of 2.1±0.2 mGy (10 MV FFF), 1.3±0.1 mGy (6 MV FFF), 0.45±0.06 mGy (6 MV) for photons, and 0.47±0.05 mGy (6 MeV) for electrons, which differed from EBT3 film measurements by 0.6%, 0.8%, 10%, and 15.7%, respectively. The half-life of the hydrated electrons in the solution was ∼ 24
μs.
Conclusions: By measuring 660-nm laser light transmitted through a cm-scale, multi-pass water cavity, we observed absorption transients consistent with hydrated electrons generated by clinical linac radiation. The agreement between our inferred dose and EBT3 film measurements suggests this proof-of-concept system represents a viable pathway toward tissue-equivalent dosimeters for clinical radiotherapy applications.
Purpose: The aims of this study were to measure the optical absorption associated with hydrated electrons produced by clinical linacs and to assess the suitability of the technique for radiotherapy (⩽ 1 cGy per pulse) applications.
Methods: 40 mW of 660-nm laser light was sent five passes through deionized water contained in a 10
×
4
×
2 cm3 glass-walled cavity by using four broadband dielectric mirrors, two on each side of the cavity. The light was collected with a biased silicon photodetector. The water cavity was then irradiated by a Varian TrueBeam linac with both photon (10 MV FFF, 6 MV FFF, 6 MV) and electron beams (6 MeV) while monitoring the transmitted laser power for absorption transients. Radiochromic EBT3 film measurements were also performed for comparison.
Results: Examination of the absorbance profiles showed clear absorption changes in the water when radiation pulses were delivered. Both the amplitude and the decay time of the signal appeared consistent with the absorbed dose and the characteristics of the hydrated electrons. By using literature value for the hydrated electron radiation chemical yield (3.0±0.3), we inferred doses of 2.1±0.2 mGy (10 MV FFF), 1.3±0.1 mGy (6 MV FFF), 0.45±0.06 mGy (6 MV) for photons, and 0.47±0.05 mGy (6 MeV) for electrons, which differed from EBT3 film measurements by 0.6%, 0.8%, 10%, and 15.7%, respectively. The half-life of the hydrated electrons in the solution was ∼ 24
μs.
Conclusions: By measuring 660-nm laser light transmitted through a cm-scale, multi-pass water cavity, we observed absorption transients consistent with hydrated electrons generated by clinical linac radiation. The agreement between our inferred dose and EBT3 film measurements suggests this proof-of-concept system represents a viable pathway toward tissue-equivalent dosimeters for clinical radiotherapy applications.
Bian, Jingyi; Duran, Juan; Shin, Wook-Geun; Ramos-Méndez, Jose; Sankey, Jack C; Childress, Lilian; Seuntjens, Jan; Enger, Shirin A
GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species Journal Article
In: Physics in Medicine & Biology, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey_33,
title = {GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species},
author = {Jingyi Bian and Juan Duran and Wook-Geun Shin and Jose Ramos-Méndez and Jack C Sankey and Lilian Childress and Jan Seuntjens and Shirin A Enger
},
doi = {10.1088/1361-6560/acd90d},
year = {2023},
date = {2023-06-15},
journal = {Physics in Medicine & Biology},
abstract = {Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.
Berumen, Francisco; Enger, Shirin A.; Beaulieu, Luc
Fast DM,M calculation in LDR brachytherapy using deep learning methods Journal Article
In: Physics in Medicine & Biology, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey_29,
title = {Fast DM,M calculation in LDR brachytherapy using deep learning methods},
author = {Francisco Berumen and Shirin A. Enger and Luc Beaulieu},
doi = {10.1088/1361-6560/accd42},
year = {2023},
date = {2023-05-23},
urldate = {2023-05-23},
journal = {Physics in Medicine & Biology},
abstract = {Objective.The Monte Carlo (MC) method provides a complete solution to the tissue heterogeneity effects in low-energy low-dose rate (LDR) brachytherapy. However, long computation times limit the clinical implementation of MC-based treatment planning solutions. This work aims to apply deep learning (DL) methods, specifically a model trained with MC simulations, to predict accurate dose to medium in medium (DM,M) distributions in LDR prostate brachytherapy.Approach.To train the DL model, 2369 single-seed configurations, corresponding to 44 prostate patient plans, were used. These patients underwent LDR brachytherapy treatments in which125I SelectSeed sources were implanted. For each seed configuration, the patient geometry, the MC dose volume and the single-seed plan volume were used to train a 3D Unet convolutional neural network. Previous knowledge was included in the network as anr2kernel related to the first-order dose dependency in brachytherapy. MC and DL dose distributions were compared through the dose maps, isodose lines, and dose-volume histograms. Features enclosed in the model were visualized.Main results.Model features started from the symmetrical kernel and finalized with an anisotropic representation that considered the patient organs and their interfaces, the source position, and the low- and high-dose regions. For a full prostate patient, small differences were seen below the 20% isodose line. When comparing DL-based and MC-based calculations, the predicted CTVD90metric had an average difference of -0.1%. Average differences for OARs were -1.3%, 0.07%, and 4.9% for the rectumD2cc, the bladderD2cc, and the urethraD0.1cc. The model took 1.8 ms to predict a complete 3DDM,Mvolume (1.18 M voxels).Significance.The proposed DL model stands for a simple and fast engine which includes prior physics knowledge of the problem. Such an engine considers the anisotropy of a brachytherapy source and the patient tissue composition.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Objective.The Monte Carlo (MC) method provides a complete solution to the tissue heterogeneity effects in low-energy low-dose rate (LDR) brachytherapy. However, long computation times limit the clinical implementation of MC-based treatment planning solutions. This work aims to apply deep learning (DL) methods, specifically a model trained with MC simulations, to predict accurate dose to medium in medium (DM,M) distributions in LDR prostate brachytherapy.Approach.To train the DL model, 2369 single-seed configurations, corresponding to 44 prostate patient plans, were used. These patients underwent LDR brachytherapy treatments in which125I SelectSeed sources were implanted. For each seed configuration, the patient geometry, the MC dose volume and the single-seed plan volume were used to train a 3D Unet convolutional neural network. Previous knowledge was included in the network as anr2kernel related to the first-order dose dependency in brachytherapy. MC and DL dose distributions were compared through the dose maps, isodose lines, and dose-volume histograms. Features enclosed in the model were visualized.Main results.Model features started from the symmetrical kernel and finalized with an anisotropic representation that considered the patient organs and their interfaces, the source position, and the low- and high-dose regions. For a full prostate patient, small differences were seen below the 20% isodose line. When comparing DL-based and MC-based calculations, the predicted CTVD90metric had an average difference of -0.1%. Average differences for OARs were -1.3%, 0.07%, and 4.9% for the rectumD2cc, the bladderD2cc, and the urethraD0.1cc. The model took 1.8 ms to predict a complete 3DDM,Mvolume (1.18 M voxels).Significance.The proposed DL model stands for a simple and fast engine which includes prior physics knowledge of the problem. Such an engine considers the anisotropy of a brachytherapy source and the patient tissue composition.
Peppa, Vasiliki; Thomson, Rowan M.; Enger, Shirin A.; Fonseca, Gabriel P.; Lee, Choonik; Lucero, Joseph N. E.; Mourtada, Firas; Siebert, Frank-André; Vijande, Javier; Papagiannis, Panagiotis
A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy Journal Article
In: Medical Physics, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey_34,
title = {A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy},
author = {Vasiliki Peppa and Rowan M. Thomson and Shirin A. Enger and Gabriel P. Fonseca and Choonik Lee and Joseph N. E. Lucero and Firas Mourtada and Frank-André Siebert and Javier Vijande and Panagiotis Papagiannis},
doi = {10.1002/mp.16455},
year = {2023},
date = {2023-05-17},
urldate = {2023-05-17},
journal = {Medical Physics},
abstract = {Purpose: To provide the first clinical test case for commissioning of 192 Ir brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow.
Acquisition and validation methods: A computational patient phantom model was generated from a clinical multi-catheter 192 Ir HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic 192 Ir HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution.
Data format and usage notes: The dataset is available online at http://irochouston.mdanderson.org/rpc/BrachySeeds/BrachySeeds/index.html,https://doi.org/10.52519/00005. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations.
Potential applications: The dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Purpose: To provide the first clinical test case for commissioning of 192 Ir brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow.
Acquisition and validation methods: A computational patient phantom model was generated from a clinical multi-catheter 192 Ir HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic 192 Ir HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution.
Data format and usage notes: The dataset is available online at http://irochouston.mdanderson.org/rpc/BrachySeeds/BrachySeeds/index.html,https://doi.org/10.52519/00005. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations.
Potential applications: The dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.
Acquisition and validation methods: A computational patient phantom model was generated from a clinical multi-catheter 192 Ir HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic 192 Ir HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution.
Data format and usage notes: The dataset is available online at http://irochouston.mdanderson.org/rpc/BrachySeeds/BrachySeeds/index.html,https://doi.org/10.52519/00005. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations.
Potential applications: The dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.
Antaki, Majd; Renaud, Marc-André; Morcos, Marc; Seuntjens, Jan; Enger, Shirin A.
Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem Journal Article
In: Physics in Medicine & Biology, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey_28,
title = {Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem},
author = {Majd Antaki and Marc-André Renaud and Marc Morcos and Jan Seuntjens and Shirin A. Enger },
doi = {10.1088/1361-6560/acbc63},
year = {2023},
date = {2023-03-13},
urldate = {2023-03-13},
journal = {Physics in Medicine & Biology},
abstract = {Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.
Bui, Alaina; Bekerat, Hamed; Childress, Lilian; Sankey, Jack; Seuntjens, Jan; Enger, Shirin A.
Effects of incoming particle energy and cluster size on the G-value of hydrated electrons Journal Article
In: Physics in Medicine & Biology, 2023.
Abstract | Links | BibTeX | Tags:
@article{nokey_27,
title = {Effects of incoming particle energy and cluster size on the G-value of hydrated electrons},
author = {Alaina Bui and Hamed Bekerat and Lilian Childress and Jack Sankey and Jan Seuntjens and Shirin A. Enger},
url = {https://doi.org/10.1016/j.ejmp.2023.102540},
doi = {10.1016/j.ejmp.2023.102540},
year = {2023},
date = {2023-02-16},
urldate = {2023-02-16},
journal = {Physics in Medicine & Biology},
abstract = {In hydrated electron (e-aq) dosimetry, absorbed radiation dose to water is measured by monitoring the concentration of radiation-induced e-aq. However, to obtain accurate dose, the radiation chemical yield of e-aq, G(e-aq), is needed for the radiation quality/setup under investigation. The aim of this study was to investigate the time-evolution of the G-values for the main generated reactive species during water radiolysis using GEANT4-DNA. The effects of cluster size and linear energy transfer (LET) on G(e-aq) were examined. Validity of GEANT4-DNA for calculation of G(e-aq) for clinically relevant energies was studied. Three scenarios were investigated with different phantom sizes and incoming electron energies (1 keV to 1 MeV). The time evolution of G(e-aq) was in good agreement with published data and did not change with decreasing phantom size. The time-evolution of the G-values increases with increasing LET for all radiolytic species. The particle tracks formed with high-energy electrons are separated and the resulting reactive species develop independently in time. With decreasing energy, the mean separation distance between reactive species decreases. The particle tracks might not initially overlap but will overlap shortly thereafter due to diffusion of reactive species, increasing the probability of e-aq recombination with other species. This also explains the decrease of G(e-aq) with cluster size and LET. Finally, if all factors are kept constant, as the incoming electron energy increases to clinically relevant energies, G(e-aq) remains similar to its value at 1 MeV, hence GEANT4-DNA can be used for clinically relevant energies.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
In hydrated electron (e-aq) dosimetry, absorbed radiation dose to water is measured by monitoring the concentration of radiation-induced e-aq. However, to obtain accurate dose, the radiation chemical yield of e-aq, G(e-aq), is needed for the radiation quality/setup under investigation. The aim of this study was to investigate the time-evolution of the G-values for the main generated reactive species during water radiolysis using GEANT4-DNA. The effects of cluster size and linear energy transfer (LET) on G(e-aq) were examined. Validity of GEANT4-DNA for calculation of G(e-aq) for clinically relevant energies was studied. Three scenarios were investigated with different phantom sizes and incoming electron energies (1 keV to 1 MeV). The time evolution of G(e-aq) was in good agreement with published data and did not change with decreasing phantom size. The time-evolution of the G-values increases with increasing LET for all radiolytic species. The particle tracks formed with high-energy electrons are separated and the resulting reactive species develop independently in time. With decreasing energy, the mean separation distance between reactive species decreases. The particle tracks might not initially overlap but will overlap shortly thereafter due to diffusion of reactive species, increasing the probability of e-aq recombination with other species. This also explains the decrease of G(e-aq) with cluster size and LET. Finally, if all factors are kept constant, as the incoming electron energy increases to clinically relevant energies, G(e-aq) remains similar to its value at 1 MeV, hence GEANT4-DNA can be used for clinically relevant energies.
2022
Morén, Björn; Antaki, Majd; Famulari, Gabriel; Morcos, Marc; Larsson, Torbjörn; Enger, Shirin A; Tedgren, Åsa Carlsson
In: Medical Physics, 2022.
Abstract | Links | BibTeX | Tags:
@article{nokey_30,
title = {Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy},
author = {Björn Morén and Majd Antaki and Gabriel Famulari and Marc Morcos and Torbjörn Larsson and Shirin A Enger and Åsa Carlsson Tedgren},
doi = {10.1002/mp.16134},
year = {2022},
date = {2022-12-22},
journal = {Medical Physics},
abstract = {Background: Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT).
Purpose: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.
Methods: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.
Results: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra (D0.1cc) and rectum (D0.1cc and D1cc) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.
Conclusions: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Background: Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT).
Purpose: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.
Methods: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.
Results: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra (D0.1cc) and rectum (D0.1cc and D1cc) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.
Conclusions: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.
Purpose: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.
Methods: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.
Results: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra (D0.1cc) and rectum (D0.1cc and D1cc) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.
Conclusions: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.
Carroll, Liam; Enger, Shirin A.
Compact radiation detector could expedite the use of dynamic PET Online
2022, visited: 16.12.2022.
@online{nokey_31,
title = {Compact radiation detector could expedite the use of dynamic PET},
author = {Liam Carroll and Shirin A. Enger},
url = {https://physicsworld.com/a/compact-radiation-detector-could-expedite-the-use-of-dynamic-pet/},
doi = {10.1002/mp.16055},
year = {2022},
date = {2022-12-16},
urldate = {2022-12-16},
journal = {Medical Physics},
keywords = {},
pubstate = {published},
tppubtype = {online}
}
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.
Abstract | Links | BibTeX | Tags:
@article{nokey,
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},
url = {https://doi-org.proxy3.library.mcgill.ca/10.1002/mp.16055
},
doi = {10.1002/mp.16055},
year = {2022},
date = {2022-10-17},
urldate = {2022-10-17},
journal = {Medical Physics },
abstract = {Background:
Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.
Purpose:
The aim of this study was to A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.
Methods:
The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: 18F, 11C, 15O and 68Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half max (FWHM) of the signal. In addition, contribution by the annihilation photons vs positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.
Results:
The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for 11C simulation. All simulations had a significant annihilation photon contribution. The signal for 18F and 11C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.
Conclusions:
The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Background:
Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.
Purpose:
The aim of this study was to A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.
Methods:
The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: 18F, 11C, 15O and 68Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half max (FWHM) of the signal. In addition, contribution by the annihilation photons vs positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.
Results:
The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for 11C simulation. All simulations had a significant annihilation photon contribution. The signal for 18F and 11C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.
Conclusions:
The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.
Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.
Purpose:
The aim of this study was to A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.
Methods:
The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: 18F, 11C, 15O and 68Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half max (FWHM) of the signal. In addition, contribution by the annihilation photons vs positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.
Results:
The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for 11C simulation. All simulations had a significant annihilation photon contribution. The signal for 18F and 11C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.
Conclusions:
The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.
Sebastien, Quetin
Artificial Intelligence-based Brachytherapy Presentation
From New avenues in the non-operative management of patients with rectal cancer Conference, 14.10.2022.
BibTeX | Tags:
@misc{nokey,
title = {Artificial Intelligence-based Brachytherapy},
author = {Quetin Sebastien},
editor = {New avenues in the non-operative management of patients with rectal cancer: Time for discussion},
year = {2022},
date = {2022-10-14},
urldate = {2022-10-14},
howpublished = {From New avenues in the non-operative management of patients with rectal cancer Conference},
keywords = {},
pubstate = {published},
tppubtype = {presentation}
}
Zou, Yujing
TransMedTech Excellence Scholarship (Doctoral Award) Journal Article
In: 2022.
@article{nokey,
title = {TransMedTech Excellence Scholarship (Doctoral Award)},
author = {Yujing Zou },
url = {https://transmedtech.org/en/training/transmedtech-institute-excellence-scholarships/},
year = {2022},
date = {2022-09-01},
urldate = {2022-09-01},
organization = {The Institut TransMedTech},
key = {award},
type = {award},
keywords = {award},
pubstate = {published},
tppubtype = {article}
}
Kalinowski, Jonathan
McGill Faculty of Medicine and Health Sciences Internal Studentship award
2022.
Abstract | Links | BibTeX | Tags:
@award{nokey,
title = {McGill Faculty of Medicine and Health Sciences Internal Studentship},
author = {Jonathan Kalinowski},
url = {https://www.mcgill.ca/medhealthsci-gradstudies/funding-opportunities/graduate-students/internal-studentships},
year = {2022},
date = {2022-08-15},
urldate = {2022-08-15},
organization = {McGill University Faculty of Medicine and Health Sciences},
abstract = {Internal Studentships are open to highly qualified Faculty of Medicine graduate students who are registered full-time in a research training program (Thesis) leading to an M.Sc or PhD degree.
},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Internal Studentships are open to highly qualified Faculty of Medicine graduate students who are registered full-time in a research training program (Thesis) leading to an M.Sc or PhD degree.
Rahbaran, Maryam
Graduate Excellence Award award
2022.
@award{nokey,
title = {Graduate Excellence Award },
author = {Maryam Rahbaran},
year = {2022},
date = {2022-08-08},
abstract = {Merit-based recruitment award for first year MSc students. },
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Merit-based recruitment award for first year MSc students.
Sebastien, Quetin; Zou, Yujing
Deep Learning Framework : Tensorboard and Pytorch Lightning Workshop
2022.
@workshop{nokey,
title = {Deep Learning Framework : Tensorboard and Pytorch Lightning},
author = {Quetin Sebastien and Yujing Zou},
url = {https://www.youtube.com/watch?v=8q09b-Yqly4&list=PLVH7T2_su-vkHLGQXJ0gHijbhjLJOCbaq&index=18
https://mcmedhacks.com/},
year = {2022},
date = {2022-07-20},
urldate = {2022-07-20},
howpublished = {from McMedHacks 2022},
keywords = {},
pubstate = {published},
tppubtype = {workshop}
}
Zou, Yujing; Alvarez, David-Santiago Ayala
ESTRO 2022: innovations in brachytherapy Journal Article
In: ESTRO Newsletter, pp. 754–767, 2022.
@article{zou2022brachyestro,
title = {ESTRO 2022: innovations in brachytherapy},
author = {Yujing Zou and David-Santiago Ayala Alvarez },
url = {https://www.estro.org/About/Newsroom/Newsletter/Brachytheraphy/ESTRO-2022-innovations-in-brachytherapy-Brachyther},
year = {2022},
date = {2022-06-29},
urldate = {2022-06-29},
journal = {ESTRO Newsletter},
pages = {754--767},
publisher = {The European SocieTy for Radiotherapy and Oncology},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Sebastien, Quetin
Deep Learning Framework : Pytorch tensors and Autograd Workshop
2022.
@workshop{nokey,
title = {Deep Learning Framework : Pytorch tensors and Autograd},
author = {Quetin Sebastien },
url = {https://www.youtube.com/watch?v=3X0ZEfY-nuc&list=PLVH7T2_su-vkHLGQXJ0gHijbhjLJOCbaq&index=12
https://mcmedhacks.com/},
year = {2022},
date = {2022-06-29},
urldate = {2022-06-29},
howpublished = {from McMedHacks 2022},
keywords = {},
pubstate = {published},
tppubtype = {workshop}
}
Sebastien, Quetin; Zou, Yujing
Introduction to medical image processing with Python : DICOM and histopathology images Workshop
2022.
@workshop{nokey,
title = {Introduction to medical image processing with Python : DICOM and histopathology images},
author = {Quetin Sebastien and Yujing Zou},
url = {https://www.youtube.com/watch?v=oazONk9JpFg&list=PLVH7T2_su-vkHLGQXJ0gHijbhjLJOCbaq&index=7
https://mcmedhacks.com/},
year = {2022},
date = {2022-06-17},
urldate = {2022-06-17},
howpublished = {from McMedHacks 2022},
keywords = {},
pubstate = {published},
tppubtype = {workshop}
}
Zou, Yujing
Fonds de recherche du Québec - Santé (FRQS) PhD doctoral training scholarship award
2022, ($84,000 for 2022 - 2026).
@award{nokey,
title = {Fonds de recherche du Québec - Santé (FRQS) PhD doctoral training scholarship},
author = {Yujing Zou},
url = {https://repertoire.frq.gouv.qc.ca/offres/rechercheOffres.do?methode=afficher},
year = {2022},
date = {2022-06-02},
urldate = {2022-06-02},
organization = {Fonds de recherche du Québec},
note = {$84,000 for 2022 - 2026},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Zou, Yujing
Fonds de recherche du Québec - Nature et technologies (FRQNT) PhD doctoral training scholarship award
2022, ($84,000 for 2022 - 2026; Declined due to acceptance of FRQS).
@award{nokey,
title = {Fonds de recherche du Québec - Nature et technologies (FRQNT) PhD doctoral training scholarship},
author = {Yujing Zou},
url = {https://repertoire.frq.gouv.qc.ca/offres/rechercheOffres.do?methode=afficher},
year = {2022},
date = {2022-06-02},
urldate = {2022-06-02},
organization = {Fonds de recherche du Québec},
note = {$84,000 for 2022 - 2026; Declined due to acceptance of FRQS},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Zou, Yujing; Lecavalier-Barsoum, Magali; Pelmus, Manuela; Maleki, Farhad; Enger, Shirin A.
Curietherapies 2022.
Abstract | Links | BibTeX | Tags:
@conference{nokey,
title = {Predictive modeling of post radiation-therapy recurrence for gynecological cancer patients using clinical and histopathology imaging features},
author = {Yujing Zou and Magali Lecavalier-Barsoum and Manuela Pelmus and Farhad Maleki and Shirin A. Enger},
url = {https://www.researchgate.net/publication/361436138_Predictive_modeling_of_post_radiation-therapy_recurrence_for_gynecological_cancer_patients_using_clinical_and_histopathology_imaging_features},
year = {2022},
date = {2022-05-23},
urldate = {2022-05-23},
organization = {Curietherapies},
abstract = {Purpose: To build a machine-learning (ML) classifier to predict the clinical endpoint of post-Radiation-Therapy (RT) recurrence of gynecological cancer patients, while exploring the outcome predictability of cell spacing and nuclei size pre-treatment histopathology image features and clinical variables. Materials and Methods: Thirty-six gynecological (i.e., cervix, vaginal, and vulva) cancer patients (median age at diagnosis = 59.5 years) with a median follow-up time of 25.7 months, nine of which (event rate of 25%) experienced post-RT recurrence, were included in this analysis. Patient-specific nuclei size and cell spacing distributions from cancerous and non-tumoral regions of pre-treatment hematoxylin and eosin (H&E) stained digital histopathology Whole-Slide-Images (WSI) were extracted. The mean and standard deviation of these distributions were computed as imaging features for each WSI. Clinical features of clinical and radiological stage at the time of radiation, p16 status, age at diagnosis, and cancer type were also obtained. Uniquely, a Tree-based Pipeline Optimization Tool (TPOT) AutoML approach, including hyperparameter tuning, was implemented to find the best performing pipeline for this class-imbalanced and small dataset. A Radial Basis Function Kernel (RBF) sampler (gamma = 0.25) was applied to combined imaging and clinical input variables for training. The resulting features were fed into an XGBoost (ie., eXtreme gradient-boosting) classifier (learning rate = 0.1). Its outputs were propagated as “synthetic features” followed by polynomial feature transforms. All raw and transformed features were trained with a decision tree classification algorithm. Results of model evaluation metrics from a 10-fold stratified shuffle split cross-validation were averaged. A permutation test (n=1000) was performed to validate the significance of the classification scores. Results: Our model achieved a 10-fold stratified shuffle split cross-validation scores of 0.87 for mean accuracy, 0.92 for mean balanced accuracy, 0.78 for precision, 1 for recall, 0.85 for F1 score, and 0.92 for Area Under the Curve of Receiver Operating Characteristics Curve, to predict our patient cohort’s post-RT recurrence binary outcome. A p-value of 0.036 was obtained from the permutation test. This implies real dependencies between our combined imaging and clinical features and outcomes which were learned by the classifier, and the primising model performance was not by chance. Conclusions: Despite the small dataset and low event rate, as a proof of concept, we showed that a decision-tree-based ML classification algorithm using an XGBoost algorithm is able to utilize combined (cell spacing & nuclei size) imaging and clinical features to predict post-RT outcomes for gynecological cancer patients.},
keywords = {},
pubstate = {published},
tppubtype = {conference}
}
Purpose: To build a machine-learning (ML) classifier to predict the clinical endpoint of post-Radiation-Therapy (RT) recurrence of gynecological cancer patients, while exploring the outcome predictability of cell spacing and nuclei size pre-treatment histopathology image features and clinical variables. Materials and Methods: Thirty-six gynecological (i.e., cervix, vaginal, and vulva) cancer patients (median age at diagnosis = 59.5 years) with a median follow-up time of 25.7 months, nine of which (event rate of 25%) experienced post-RT recurrence, were included in this analysis. Patient-specific nuclei size and cell spacing distributions from cancerous and non-tumoral regions of pre-treatment hematoxylin and eosin (H&E) stained digital histopathology Whole-Slide-Images (WSI) were extracted. The mean and standard deviation of these distributions were computed as imaging features for each WSI. Clinical features of clinical and radiological stage at the time of radiation, p16 status, age at diagnosis, and cancer type were also obtained. Uniquely, a Tree-based Pipeline Optimization Tool (TPOT) AutoML approach, including hyperparameter tuning, was implemented to find the best performing pipeline for this class-imbalanced and small dataset. A Radial Basis Function Kernel (RBF) sampler (gamma = 0.25) was applied to combined imaging and clinical input variables for training. The resulting features were fed into an XGBoost (ie., eXtreme gradient-boosting) classifier (learning rate = 0.1). Its outputs were propagated as “synthetic features” followed by polynomial feature transforms. All raw and transformed features were trained with a decision tree classification algorithm. Results of model evaluation metrics from a 10-fold stratified shuffle split cross-validation were averaged. A permutation test (n=1000) was performed to validate the significance of the classification scores. Results: Our model achieved a 10-fold stratified shuffle split cross-validation scores of 0.87 for mean accuracy, 0.92 for mean balanced accuracy, 0.78 for precision, 1 for recall, 0.85 for F1 score, and 0.92 for Area Under the Curve of Receiver Operating Characteristics Curve, to predict our patient cohort’s post-RT recurrence binary outcome. A p-value of 0.036 was obtained from the permutation test. This implies real dependencies between our combined imaging and clinical features and outcomes which were learned by the classifier, and the primising model performance was not by chance. Conclusions: Despite the small dataset and low event rate, as a proof of concept, we showed that a decision-tree-based ML classification algorithm using an XGBoost algorithm is able to utilize combined (cell spacing & nuclei size) imaging and clinical features to predict post-RT outcomes for gynecological cancer patients.
Zou, Yujing; Lecavalier-Barsoum, Magali; Pelmus, Manuela; Maleki, Farhad; Enger, Shirin A.
Young Investigator Competition Winner at the Curietherapies Conference award
2022.
Abstract | Links | BibTeX | Tags:
@award{nokey,
title = {Young Investigator Competition Winner at the Curietherapies Conference },
author = {Yujing Zou and Magali Lecavalier-Barsoum and Manuela Pelmus and Farhad Maleki and Shirin A. Enger },
url = {https://www.researchgate.net/publication/360979157_SP-0014_McMedHacks_Deep_learning_for_medical_image_analysis_workshops_and_Hackathon_in_radiation_oncology},
year = {2022},
date = {2022-05-23},
urldate = {2022-05-23},
organization = {Curietherapies},
abstract = {Purpose/Objective: The McMedHacks workshop and presentation series was created to teach individuals from various backgrounds about deep learning (DL) for medical image analysis in May, 2021. Material/Methods: McMedHacks is a free and student-led 8-week summer program. Registration for the event was open to everyone, including a form to survey participants’ area of expertise, country of origin, level of study, and level of programming skills. The weekly workshops were instructed by 8 students and experts assisted by 20 mentors who provided weekly tutorials. Recent developments in DL and medical physics were highlighted by 21 leaders from industry and academia. A virtual grand challenge Hackathon took place at the end of the workshop series. All events were held virtually and recorded on Zoom to accommodate all time zones and locations. The workshops were designed as interactive coding demos and shared through Google Colab notebooks. Results: McMedHacks gained 356 registrations from participants of 38 different countries (Fig. 1) from undergraduates, to PhDs and MDs. A vast number of disciplines and professions were represented, dominated by medical physics students, academic, and clinical medical physicists (Fig. 2). Sixty-nine participants earned a certificate of completion by having engaged with at least 12 of all 14 events. The program received participant feedback average scores of 4.768, 4.478, 4.579, 4.292, 4.84 out of five for the qualities of presentation, workshop session, tutorial and mentor, assignments, and course delivery, respectively. The eight-week long workshop’s duration allowed participants to digest the taught materials in a continuous manner as opposed to bootcamp-style conference workshops. Conclusion: The overwhelming interest and engagement for the McMedHacks workshop series from the Radiation Oncology (RadOnc) community illustrates a demand for Artificial Intelligence (AI) education in RadOnc. The future of RadOnc clinics will inevitably integrate AI. Therefore, current RadOnc professionals, and student and resident trainees should be prepared to understand basic AI principles and its applications to troubleshoot, innovate, and collaborate. McMedHacks set an excellent example of promoting open and multidisciplinary education, scientific communication, and leadership for integrating AI education into the RadOnc community on an international level. Therefore, we advocate for implementation of AI curriculums in professional education programs such as Commission on Accreditation of Medical Physics Education Programs (CAMPEP). Furthermore, we encourage experts from around the world in the field of AI, or RadOnc, or both, to take initiatives like McMedHacks to collaborate and push forward AI education in their departments and lead practical workshops, regardless of their levels of education.},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Purpose/Objective: The McMedHacks workshop and presentation series was created to teach individuals from various backgrounds about deep learning (DL) for medical image analysis in May, 2021. Material/Methods: McMedHacks is a free and student-led 8-week summer program. Registration for the event was open to everyone, including a form to survey participants’ area of expertise, country of origin, level of study, and level of programming skills. The weekly workshops were instructed by 8 students and experts assisted by 20 mentors who provided weekly tutorials. Recent developments in DL and medical physics were highlighted by 21 leaders from industry and academia. A virtual grand challenge Hackathon took place at the end of the workshop series. All events were held virtually and recorded on Zoom to accommodate all time zones and locations. The workshops were designed as interactive coding demos and shared through Google Colab notebooks. Results: McMedHacks gained 356 registrations from participants of 38 different countries (Fig. 1) from undergraduates, to PhDs and MDs. A vast number of disciplines and professions were represented, dominated by medical physics students, academic, and clinical medical physicists (Fig. 2). Sixty-nine participants earned a certificate of completion by having engaged with at least 12 of all 14 events. The program received participant feedback average scores of 4.768, 4.478, 4.579, 4.292, 4.84 out of five for the qualities of presentation, workshop session, tutorial and mentor, assignments, and course delivery, respectively. The eight-week long workshop’s duration allowed participants to digest the taught materials in a continuous manner as opposed to bootcamp-style conference workshops. Conclusion: The overwhelming interest and engagement for the McMedHacks workshop series from the Radiation Oncology (RadOnc) community illustrates a demand for Artificial Intelligence (AI) education in RadOnc. The future of RadOnc clinics will inevitably integrate AI. Therefore, current RadOnc professionals, and student and resident trainees should be prepared to understand basic AI principles and its applications to troubleshoot, innovate, and collaborate. McMedHacks set an excellent example of promoting open and multidisciplinary education, scientific communication, and leadership for integrating AI education into the RadOnc community on an international level. Therefore, we advocate for implementation of AI curriculums in professional education programs such as Commission on Accreditation of Medical Physics Education Programs (CAMPEP). Furthermore, we encourage experts from around the world in the field of AI, or RadOnc, or both, to take initiatives like McMedHacks to collaborate and push forward AI education in their departments and lead practical workshops, regardless of their levels of education.
Zou, Yujing
Biological & Biomedical Engineering PhD Recruitment Award award
2022, (These awards are designed to help recruit top students to our program and are offered to applicants wishing to start in Fall or Winter. The standard recruitment awards are $10,000/year for three years for Doctoral students.).
@award{nokey,
title = {Biological & Biomedical Engineering PhD Recruitment Award },
author = {Yujing Zou},
url = {https://www.mcgill.ca/bbme/programs/funding#BME-Recruitment-Award},
year = {2022},
date = {2022-05-10},
urldate = {2022-05-10},
organization = {McGill Biological & Biomedical Engineering department },
note = {These awards are designed to help recruit top students to our program and are offered to applicants wishing to start in Fall or Winter. The standard recruitment awards are $10,000/year for three years for Doctoral students.},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Jafarzadeh, Hossein
Biological & Biomedical Engineering PhD Recruitment Award award
2022.
@award{nokey,
title = {Biological & Biomedical Engineering PhD Recruitment Award },
author = {Hossein Jafarzadeh },
url = {https://www.mcgill.ca/bbme/programs/funding#BME-Recruitment-Award},
year = {2022},
date = {2022-05-10},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
Enger, Shirin A.; Famulari, Gabriel
Delivery system for intensity modulated high dose rate brachytherapy with intermediate energy brachytherapy isotopes Patent
2022, (US Patent 11,324,966).
BibTeX | Tags:
@patent{enger2022delivery,
title = {Delivery system for intensity modulated high dose rate brachytherapy with intermediate energy brachytherapy isotopes},
author = {Shirin A. Enger and Gabriel Famulari},
year = {2022},
date = {2022-05-01},
urldate = {2022-05-01},
publisher = {Google Patents},
note = {US Patent 11,324,966},
keywords = {},
pubstate = {published},
tppubtype = {patent}
}
Zou, Yujing
Graduate Research Enhancement and Travel Awards (GREAT Awards) award
2022, (In 2009, Graduate and Postdoctoral Studies introduced the Graduate Research Enhancement and Travel awards (GREAT awards) program in consultation with the Faculty Deans and Associate Deans. These awards cover dissemination of research through graduate student presentations at conferences, and other graduate student research-enhancement activities, such as travel for fieldwork, archival inquiry and extra-mural collaborative research. GREAT budgets are allocated to Faculties each year, and as such, are managed directly by the Associate Dean's office.).
@award{nokey,
title = {Graduate Research Enhancement and Travel Awards (GREAT Awards)},
author = {Yujing Zou},
url = {https://www.mcgill.ca/gps/funding/fac-staff/awards/great},
year = {2022},
date = {2022-04-14},
note = {In 2009, Graduate and Postdoctoral Studies introduced the Graduate Research Enhancement and Travel awards (GREAT awards) program in consultation with the Faculty Deans and Associate Deans. These awards cover dissemination of research through graduate student presentations at conferences, and other graduate student research-enhancement activities, such as travel for fieldwork, archival inquiry and extra-mural collaborative research. GREAT budgets are allocated to Faculties each year, and as such, are managed directly by the Associate Dean's office.},
keywords = {},
pubstate = {published},
tppubtype = {award}
}
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.
Rahbaran, Maryam; Kalinowski, Jonathan; Tsui, James Man Git; DeCunha, Joseph; Enger, Shirin A.
Monte-Carlo Based Simulations of the Uncertainties in Clinical Water-Based Intravascular Brachytherapy Dosimetry Presentation
11.04.2022.
@misc{nokey,
title = {Monte-Carlo Based Simulations of the Uncertainties in Clinical Water-Based Intravascular Brachytherapy Dosimetry},
author = {Maryam Rahbaran and Jonathan Kalinowski and James Man Git Tsui and Joseph DeCunha and Shirin A. Enger},
year = {2022},
date = {2022-04-11},
urldate = {2022-04-11},
journal = {MCMA},
abstract = {"Introduction
Coronary artery disease (CAD) is the most common form of cardiovascular disease and is caused by excess plaque along the arterial wall, blocking blood flow to the heart (stenosis). Percutaneous transluminal coronary angioplasty widens a narrowed artery, leaving behind metal stents (1). However, in-stent restenosis (ISR) may occur due to damage to the arterial wall tissue, triggering neointimal hyperplasia which produces fibrotic and calcified plaques, narrowing the artery again. Drug-eluting stents (DES) slowly release medication to inhibit neointimal hyperplasia to prevent ISR but they fail in 3% to 20% of cases (2). Intravascular brachytherapy (IVBT), which uses b-emitting radionuclides to prevent ISR, is used in these failed cases. However, current dosimetry for IVBT is water based and does not consider attenuation of the radiation by heterogeneities such as the IVBT device guidewire, non-uniform distribution of calcified plaques, and stent material, or the angular dependence of dose distribution (3, 4, 5). The aim of this study was to investigate the uncertainties in clinical water based IVBT dosimetry, considering the effect of heterogeneities on dose distribution.
Materials & Methods
An inhouse Monte-Carlo based dosimetry package for IVBT applications based on Geant4 10.04 (patch 2) was developed. Patient’s artery was modelled as a 32 mm long, 8.4 mm diameter cylinder comprised of three layers: tunica media, represented with muscle, tunica intima, represented with fibrotic plaque, and tunica adventitia, represented with collagen. These layers had mass densities 1.06 g/cm3, 1.22 g/cm3 and 1.07 g/cm3 respectively. The innermost layer consisted of calcified plaque of density 1.45 g/cm3 with varying thicknesses between 0.9 and 1.9 mm with an eccentric shape and a rough surface. The stents had similar composition to Boston Scientific Synergy stents and were modelled to not overlap. The Novoste Beta-Cath 3.5F IVBT device model was used, which has a 90Sr90Y source. The geometry is shown in Figure 1a. A cylindrical scoring geometry was implemented. Two set of simulations were performed. In the first simulation called water phantom, the entire system consisted of water with unit density, and dose to water was calculated similar to the clinical water based dosimetry. In the second simulation called the artery model proper material and mass densities were assigned to each component. To ensure uncertainties below 0.8% within a 1 mm radial distance to the source and 2% within 4.2 mm from the source, 100 million decay events were simulated. The Penelope physics list was used to simulate the electromagnetic interactions between particles. Average, minimum, and maximum dose was calculated at 2.0 mm from the source center and directly and 1 mm behind the outermost stents and guidewire. Absorbed dose was normalized to 23 Gy at 2.0 mm from the source center.
Results
International Conference on Monte Carlo Techniques for Medical Applications, 2022
Compared to the water phantom (Figure 1b), average dose in the artery model (Figure 1c) was attenuated by 50.9% at 2 mm from the source centre and directly behind the guidewire and outermost stent by 66.2%, and by 69.5% 1 mm behind this region. There was significant variation in dose around the source due to the guidewire attenuating dose the most, and heterogeneous distribution of calcification.
Discussion & Conclusions
Dosimetry for IVBT based on dose rate in water is not accurate. Heterogeneities need to be considered to deliver adequate dose to the lesion area. Stent material, heterogenous distribution of calcification and the off cantered placement of the guidewire affects the uniformity of dose distribution around the source. Patients may benefit from personalized treatment planning taking dose-attenuating by different tissue/material heterogeneities into account.
References
[1] Virani, Salim, S., et al. ""Heart Disease and Stroke Statistics—2020 Update"". Circulation, vol. 141, no. 9, March 03, 2020, pp. e336. doi: 10.1161/CIR.0000000000000757.
[2] Lee M, Banka G. In-stent restenosis. Interv Cardiol Clin 2016;5: 211e220.
[3] Chiu-Tsao ST, Schaart DR, Soares CG, et al. Dose calculation formalisms and consensus dosimetry parameters for intravascular brachytherapy dosimetry: Recommendations of the AAPM Therapy Physics Committee Task Group No. 149. Med Phys 2007;34: 4126e4157.
[4] Rivard MJ, Coursey BM, DeWerd LA, et al. Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys 2004;31:633e674.
[5] Nath R, Amols H, Coffey C, et al. Intravascular brachytherapy physics: Report of the AAPM Radiation Therapy Committee Task group No. 60. Med Phys 1999;26:119e152.’
[6] Agostinelli S, Allison J, Amako K, et al. Geant4da simulation toolkit. Nucl Instrum Methods Phys Res 2003;506:230e303."},
keywords = {},
pubstate = {published},
tppubtype = {presentation}
}
"Introduction
Coronary artery disease (CAD) is the most common form of cardiovascular disease and is caused by excess plaque along the arterial wall, blocking blood flow to the heart (stenosis). Percutaneous transluminal coronary angioplasty widens a narrowed artery, leaving behind metal stents (1). However, in-stent restenosis (ISR) may occur due to damage to the arterial wall tissue, triggering neointimal hyperplasia which produces fibrotic and calcified plaques, narrowing the artery again. Drug-eluting stents (DES) slowly release medication to inhibit neointimal hyperplasia to prevent ISR but they fail in 3% to 20% of cases (2). Intravascular brachytherapy (IVBT), which uses b-emitting radionuclides to prevent ISR, is used in these failed cases. However, current dosimetry for IVBT is water based and does not consider attenuation of the radiation by heterogeneities such as the IVBT device guidewire, non-uniform distribution of calcified plaques, and stent material, or the angular dependence of dose distribution (3, 4, 5). The aim of this study was to investigate the uncertainties in clinical water based IVBT dosimetry, considering the effect of heterogeneities on dose distribution.
Materials & Methods
An inhouse Monte-Carlo based dosimetry package for IVBT applications based on Geant4 10.04 (patch 2) was developed. Patient’s artery was modelled as a 32 mm long, 8.4 mm diameter cylinder comprised of three layers: tunica media, represented with muscle, tunica intima, represented with fibrotic plaque, and tunica adventitia, represented with collagen. These layers had mass densities 1.06 g/cm3, 1.22 g/cm3 and 1.07 g/cm3 respectively. The innermost layer consisted of calcified plaque of density 1.45 g/cm3 with varying thicknesses between 0.9 and 1.9 mm with an eccentric shape and a rough surface. The stents had similar composition to Boston Scientific Synergy stents and were modelled to not overlap. The Novoste Beta-Cath 3.5F IVBT device model was used, which has a 90Sr90Y source. The geometry is shown in Figure 1a. A cylindrical scoring geometry was implemented. Two set of simulations were performed. In the first simulation called water phantom, the entire system consisted of water with unit density, and dose to water was calculated similar to the clinical water based dosimetry. In the second simulation called the artery model proper material and mass densities were assigned to each component. To ensure uncertainties below 0.8% within a 1 mm radial distance to the source and 2% within 4.2 mm from the source, 100 million decay events were simulated. The Penelope physics list was used to simulate the electromagnetic interactions between particles. Average, minimum, and maximum dose was calculated at 2.0 mm from the source center and directly and 1 mm behind the outermost stents and guidewire. Absorbed dose was normalized to 23 Gy at 2.0 mm from the source center.
Results
International Conference on Monte Carlo Techniques for Medical Applications, 2022
Compared to the water phantom (Figure 1b), average dose in the artery model (Figure 1c) was attenuated by 50.9% at 2 mm from the source centre and directly behind the guidewire and outermost stent by 66.2%, and by 69.5% 1 mm behind this region. There was significant variation in dose around the source due to the guidewire attenuating dose the most, and heterogeneous distribution of calcification.
Discussion & Conclusions
Dosimetry for IVBT based on dose rate in water is not accurate. Heterogeneities need to be considered to deliver adequate dose to the lesion area. Stent material, heterogenous distribution of calcification and the off cantered placement of the guidewire affects the uniformity of dose distribution around the source. Patients may benefit from personalized treatment planning taking dose-attenuating by different tissue/material heterogeneities into account.
References
[1] Virani, Salim, S., et al. ""Heart Disease and Stroke Statistics—2020 Update"". Circulation, vol. 141, no. 9, March 03, 2020, pp. e336. doi: 10.1161/CIR.0000000000000757.
[2] Lee M, Banka G. In-stent restenosis. Interv Cardiol Clin 2016;5: 211e220.
[3] Chiu-Tsao ST, Schaart DR, Soares CG, et al. Dose calculation formalisms and consensus dosimetry parameters for intravascular brachytherapy dosimetry: Recommendations of the AAPM Therapy Physics Committee Task Group No. 149. Med Phys 2007;34: 4126e4157.
[4] Rivard MJ, Coursey BM, DeWerd LA, et al. Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys 2004;31:633e674.
[5] Nath R, Amols H, Coffey C, et al. Intravascular brachytherapy physics: Report of the AAPM Radiation Therapy Committee Task group No. 60. Med Phys 1999;26:119e152.’
[6] Agostinelli S, Allison J, Amako K, et al. Geant4da simulation toolkit. Nucl Instrum Methods Phys Res 2003;506:230e303."
Coronary artery disease (CAD) is the most common form of cardiovascular disease and is caused by excess plaque along the arterial wall, blocking blood flow to the heart (stenosis). Percutaneous transluminal coronary angioplasty widens a narrowed artery, leaving behind metal stents (1). However, in-stent restenosis (ISR) may occur due to damage to the arterial wall tissue, triggering neointimal hyperplasia which produces fibrotic and calcified plaques, narrowing the artery again. Drug-eluting stents (DES) slowly release medication to inhibit neointimal hyperplasia to prevent ISR but they fail in 3% to 20% of cases (2). Intravascular brachytherapy (IVBT), which uses b-emitting radionuclides to prevent ISR, is used in these failed cases. However, current dosimetry for IVBT is water based and does not consider attenuation of the radiation by heterogeneities such as the IVBT device guidewire, non-uniform distribution of calcified plaques, and stent material, or the angular dependence of dose distribution (3, 4, 5). The aim of this study was to investigate the uncertainties in clinical water based IVBT dosimetry, considering the effect of heterogeneities on dose distribution.
Materials & Methods
An inhouse Monte-Carlo based dosimetry package for IVBT applications based on Geant4 10.04 (patch 2) was developed. Patient’s artery was modelled as a 32 mm long, 8.4 mm diameter cylinder comprised of three layers: tunica media, represented with muscle, tunica intima, represented with fibrotic plaque, and tunica adventitia, represented with collagen. These layers had mass densities 1.06 g/cm3, 1.22 g/cm3 and 1.07 g/cm3 respectively. The innermost layer consisted of calcified plaque of density 1.45 g/cm3 with varying thicknesses between 0.9 and 1.9 mm with an eccentric shape and a rough surface. The stents had similar composition to Boston Scientific Synergy stents and were modelled to not overlap. The Novoste Beta-Cath 3.5F IVBT device model was used, which has a 90Sr90Y source. The geometry is shown in Figure 1a. A cylindrical scoring geometry was implemented. Two set of simulations were performed. In the first simulation called water phantom, the entire system consisted of water with unit density, and dose to water was calculated similar to the clinical water based dosimetry. In the second simulation called the artery model proper material and mass densities were assigned to each component. To ensure uncertainties below 0.8% within a 1 mm radial distance to the source and 2% within 4.2 mm from the source, 100 million decay events were simulated. The Penelope physics list was used to simulate the electromagnetic interactions between particles. Average, minimum, and maximum dose was calculated at 2.0 mm from the source center and directly and 1 mm behind the outermost stents and guidewire. Absorbed dose was normalized to 23 Gy at 2.0 mm from the source center.
Results
International Conference on Monte Carlo Techniques for Medical Applications, 2022
Compared to the water phantom (Figure 1b), average dose in the artery model (Figure 1c) was attenuated by 50.9% at 2 mm from the source centre and directly behind the guidewire and outermost stent by 66.2%, and by 69.5% 1 mm behind this region. There was significant variation in dose around the source due to the guidewire attenuating dose the most, and heterogeneous distribution of calcification.
Discussion & Conclusions
Dosimetry for IVBT based on dose rate in water is not accurate. Heterogeneities need to be considered to deliver adequate dose to the lesion area. Stent material, heterogenous distribution of calcification and the off cantered placement of the guidewire affects the uniformity of dose distribution around the source. Patients may benefit from personalized treatment planning taking dose-attenuating by different tissue/material heterogeneities into account.
References
[1] Virani, Salim, S., et al. ""Heart Disease and Stroke Statistics—2020 Update"". Circulation, vol. 141, no. 9, March 03, 2020, pp. e336. doi: 10.1161/CIR.0000000000000757.
[2] Lee M, Banka G. In-stent restenosis. Interv Cardiol Clin 2016;5: 211e220.
[3] Chiu-Tsao ST, Schaart DR, Soares CG, et al. Dose calculation formalisms and consensus dosimetry parameters for intravascular brachytherapy dosimetry: Recommendations of the AAPM Therapy Physics Committee Task Group No. 149. Med Phys 2007;34: 4126e4157.
[4] Rivard MJ, Coursey BM, DeWerd LA, et al. Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys 2004;31:633e674.
[5] Nath R, Amols H, Coffey C, et al. Intravascular brachytherapy physics: Report of the AAPM Radiation Therapy Committee Task group No. 60. Med Phys 1999;26:119e152.’
[6] Agostinelli S, Allison J, Amako K, et al. Geant4da simulation toolkit. Nucl Instrum Methods Phys Res 2003;506:230e303."
Daoud, Youstina; Carroll, Liam; Enger, Shirin A.
A Radiation detector simulation toolkit for calculating the Arterial Input Function during Dynamic Positron Emission Tomography Proceedings Article
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).
"
Enger, Shirin A.; Sankey, Jack; Childress, Lilian; Megroureche, Julien
Radiation dosimeter Patent
2022, (US Patent App. 17/298,743).
BibTeX | Tags:
@patent{enger2022radiation,
title = {Radiation dosimeter},
author = {Shirin A. Enger and Jack Sankey and Lilian Childress and Julien Megroureche},
year = {2022},
date = {2022-02-01},
urldate = {2022-02-01},
publisher = {Google Patents},
note = {US Patent App. 17/298,743},
keywords = {},
pubstate = {published},
tppubtype = {patent}
}
Jafarzadeh, Hossein; Mao, Ximeng; Enger, Shirin A.
Bayesian Optimization in Treatment Planning of High Dose Rate Brachytherapy Proceedings Article
In: MEDICAL PHYSICS, pp. E200–E200, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{jafarzadeh2022bayesian,
title = {Bayesian Optimization in Treatment Planning of High Dose Rate Brachytherapy},
author = { Hossein Jafarzadeh and Ximeng Mao and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E200--E200},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Thibodeau-Antonacci, Alana; Enger, Shirin A.; Bekerat, Hamed; Vuong, Te
Gafchromic film and scintillator detector measurements in phantom with a novel intensity-modulated brachytherapy endorectal shield Proceedings Article
In: MEDICAL PHYSICS, pp. 5688–5689, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{thibodeau2022gafchromic,
title = {Gafchromic film and scintillator detector measurements in phantom with a novel intensity-modulated brachytherapy endorectal shield},
author = {Alana Thibodeau-Antonacci and Shirin A. Enger and Hamed Bekerat and Te Vuong},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {8},
pages = {5688--5689},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Martinez, Victor Daniel Diaz; Cyr, Melodie; Slobodan, Devic; Tomic, Nada; Lewis, David F; Enger, Shirin A.
Use of the Monte Carlo Method to Relate GAFCHROMIC (R) EBT3 Film Response to Absorbed Dose for Alpha Particle Dosimetry Proceedings Article
In: MEDICAL PHYSICS, pp. 5653–5653, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{martinez2022use,
title = {Use of the Monte Carlo Method to Relate GAFCHROMIC (R) EBT3 Film Response to Absorbed Dose for Alpha Particle Dosimetry},
author = {Victor Daniel Diaz Martinez and Melodie Cyr and Devic Slobodan and Nada Tomic and David F Lewis and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {8},
pages = {5653--5653},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Bian, Jingyi; Duran, Juan; Shin, Wook-Geun; Ramos-Mendez, Jose; Childress, Lilian; Sankey, Jack; Seuntjens, Jan; Enger, Shirin A.
Update of the GEANT4-DNA Software for Simulation of Radiation Chemical Yield for Reactive Water Radiolysis Species at Different Temperature and PH Proceedings Article
In: MEDICAL PHYSICS, pp. E911–E912, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{bian2022update,
title = {Update of the GEANT4-DNA Software for Simulation of Radiation Chemical Yield for Reactive Water Radiolysis Species at Different Temperature and PH},
author = {Jingyi Bian and Juan Duran and Wook-Geun Shin and Jose Ramos-Mendez and Lilian Childress and Jack Sankey and Jan Seuntjens and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E911--E912},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
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 Proceedings Article
In: MEDICAL PHYSICS, pp. 5643–5644, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@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}
}
Dev, Sachin; Behmand, Behnaz; Pater, Piotr; Enger, Shirin A.
To Establish Correlation Between Physical Microdosimetric Quantities Such as Dose Mean Lineal Energy, Target Size and Biological Endpoints Such as DNA Double Strand Breaks in HeLa Cells Irradiated with Iridium-192 High Dose Rate Brachytherapy Source and 225kV X-Rays Proceedings Article
In: MEDICAL PHYSICS, pp. E345–E345, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{dev2022establish,
title = {To Establish Correlation Between Physical Microdosimetric Quantities Such as Dose Mean Lineal Energy, Target Size and Biological Endpoints Such as DNA Double Strand Breaks in HeLa Cells Irradiated with Iridium-192 High Dose Rate Brachytherapy Source and 225kV X-Rays},
author = {Sachin Dev and Behnaz Behmand and Piotr Pater and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E345--E345},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Babik, Joud; Chabaytah, Naim; Behmand, Behnaz; Connell, Tanner; Evans, Michael; Ruo, Russell; Poirier, Y; Enger, Shirin A.
Characterization of the Relative Biological Effectiveness of a Range of Photon Energies for Irradiation of HeLa and PC-3 Cell Lines Proceedings Article
In: MEDICAL PHYSICS, pp. E980–E980, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{babik2022characterization,
title = {Characterization of the Relative Biological Effectiveness of a Range of Photon Energies for Irradiation of HeLa and PC-3 Cell Lines},
author = {Joud Babik and Naim Chabaytah and Behnaz Behmand and Tanner Connell and Michael Evans and Russell Ruo and Y Poirier and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E980--E980},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
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 Proceedings Article
In: MEDICAL PHYSICS, pp. E828–E828, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@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}
}
Berumen-Murillo, Francisco; Enger, Shirin A.; Beaulieu, Luc
Sub-Second D (M, M) Calculation for LDR Prostate Brachytherapy Using Deep Learning Methods Proceedings Article
In: MEDICAL PHYSICS, pp. E163–E163, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{berumen2022sub,
title = {Sub-Second D (M, M) Calculation for LDR Prostate Brachytherapy Using Deep Learning Methods},
author = {Francisco Berumen-Murillo and Shirin A. Enger and Luc Beaulieu},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E163--E163},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Han, Dae Yup; Damato, Antonio; Enger, Shirin A.; Roberts, K; Prisciandaro, J; Han, D; Nunez, D
Emerging and Re-Emerging Brachytherapy Treatments Proceedings Article
In: MEDICAL PHYSICS, pp. E356–E356, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{han2022emerging,
title = {Emerging and Re-Emerging Brachytherapy Treatments},
author = {Dae Yup Han and Antonio Damato and Shirin A. Enger and K Roberts and J Prisciandaro and D Han and D Nunez},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E356--E356},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
Zou, Yujing; Lecavalier-barsoum, Magali; Pelmus, Manuela; Enger, Shirin A.
Patient-Specific Nuclei Size and Cell Spacing Distribution Extraction From Histopathology Whole Slide Images for Treatment Outcome Prediction Modelling Proceedings Article
In: MEDICAL PHYSICS, pp. E266–E266, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
Abstract | Links | BibTeX | Tags:
@inproceedings{zou2022patient,
title = {Patient-Specific Nuclei Size and Cell Spacing Distribution Extraction From Histopathology Whole Slide Images for Treatment Outcome Prediction Modelling},
author = {Yujing Zou and Magali Lecavalier-barsoum and Manuela Pelmus and Shirin A. Enger},
url = {https://w4.aapm.org/meetings/2022AM/programInfo/programAbs.php?sid=10686&aid=66642},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E266--E266},
organization = {WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA},
abstract = {MO-C930-IePD-F5-1 (Monday, 7/11/2022) 9:30 AM - 10:00 AM [Eastern Time (GMT-4)]
Exhibit Hall | Forum 5
Purpose: To deliver a fully automated and generalizable approach extracting patient-specific nuclei size (ns) and cell spacing (cs) distributions from cancerous and non-tumoral regions of hematoxylin and eosin (H&E) stained digital histopathology Whole-Slide-Images (WSI) for gynecological cancer multiscale treatment outcome modelling.
Methods: Each pre-treatment gigapixel H&E WSI digitized at 40 x magnification (0.2482 microns/pixel) were divided into 5000 x 5000-pixel patches. Within each patch, the nucleus centers were identified by a difference of gaussians blob detection algorithm obtaining Delaunay triangulations and Voronoi diagrams providing cs radius. The ns radius was computed from stained pixels dominated by hematoxylin content with an automatic thresholding algorithm. With multiprocessing CPUs on a PC for each WSI, eight feature types were calculated preserving biopsies tissue heterogeneity: the mean and standard deviation of cs and ns distributions concatenated from all patches for cancerous and non-tumoral regions. This method was applied to 40 patients (1 WSI per patient) with treatment outcomes of post radiation-therapy (RT) recurrence (n = 9),and death (n = 8)).
Results: The WSI cancerous region cs distribution mean among patients without post-RT recurrence has a median of 6.64 microns, and those with post-RT recurrence with a median of 7.11 microns. This indicates the potential of utilizing such distribution features in treatment outcome prognosis modelling. Furthermore, at the third quartile, the WSI non-tumoral region ns distribution standard deviation among patients without post-RT recurrence has a value of 2.16 microns, and 1.46 microns for those with post-RT recurrence.
Conclusion: Our approach derives patient-specific microscopic data distributions from histopathology WSI that can be directly associated with retrospective patient outcomes. They are complementary and spatially orthogonal to information served by other medical imaging modalities such as CT, MR, and Ultrasound. Therefore, it has the unique potential to augment treatment outcome model inference when properly fused with radiological scans.
Funding Support, Disclosures, and Conflict of Interest: CIHR grant number 103548 and Canada Research Chairs Program (grant #252135)
Keywords
Image Analysis, Feature Extraction, Radiation Therapy
Taxonomy
IM/TH- Image Analysis (Single Modality or Multi-Modality): Imaging biomarkers and radiomics
Contact Email
yujing.zou@mail.mcgill.ca},
keywords = {},
pubstate = {published},
tppubtype = {inproceedings}
}
MO-C930-IePD-F5-1 (Monday, 7/11/2022) 9:30 AM - 10:00 AM [Eastern Time (GMT-4)]
Exhibit Hall | Forum 5
Purpose: To deliver a fully automated and generalizable approach extracting patient-specific nuclei size (ns) and cell spacing (cs) distributions from cancerous and non-tumoral regions of hematoxylin and eosin (H&E) stained digital histopathology Whole-Slide-Images (WSI) for gynecological cancer multiscale treatment outcome modelling.
Methods: Each pre-treatment gigapixel H&E WSI digitized at 40 x magnification (0.2482 microns/pixel) were divided into 5000 x 5000-pixel patches. Within each patch, the nucleus centers were identified by a difference of gaussians blob detection algorithm obtaining Delaunay triangulations and Voronoi diagrams providing cs radius. The ns radius was computed from stained pixels dominated by hematoxylin content with an automatic thresholding algorithm. With multiprocessing CPUs on a PC for each WSI, eight feature types were calculated preserving biopsies tissue heterogeneity: the mean and standard deviation of cs and ns distributions concatenated from all patches for cancerous and non-tumoral regions. This method was applied to 40 patients (1 WSI per patient) with treatment outcomes of post radiation-therapy (RT) recurrence (n = 9),and death (n = 8)).
Results: The WSI cancerous region cs distribution mean among patients without post-RT recurrence has a median of 6.64 microns, and those with post-RT recurrence with a median of 7.11 microns. This indicates the potential of utilizing such distribution features in treatment outcome prognosis modelling. Furthermore, at the third quartile, the WSI non-tumoral region ns distribution standard deviation among patients without post-RT recurrence has a value of 2.16 microns, and 1.46 microns for those with post-RT recurrence.
Conclusion: Our approach derives patient-specific microscopic data distributions from histopathology WSI that can be directly associated with retrospective patient outcomes. They are complementary and spatially orthogonal to information served by other medical imaging modalities such as CT, MR, and Ultrasound. Therefore, it has the unique potential to augment treatment outcome model inference when properly fused with radiological scans.
Funding Support, Disclosures, and Conflict of Interest: CIHR grant number 103548 and Canada Research Chairs Program (grant #252135)
Keywords
Image Analysis, Feature Extraction, Radiation Therapy
Taxonomy
IM/TH- Image Analysis (Single Modality or Multi-Modality): Imaging biomarkers and radiomics
Contact Email
yujing.zou@mail.mcgill.ca
Exhibit Hall | Forum 5
Purpose: To deliver a fully automated and generalizable approach extracting patient-specific nuclei size (ns) and cell spacing (cs) distributions from cancerous and non-tumoral regions of hematoxylin and eosin (H&E) stained digital histopathology Whole-Slide-Images (WSI) for gynecological cancer multiscale treatment outcome modelling.
Methods: Each pre-treatment gigapixel H&E WSI digitized at 40 x magnification (0.2482 microns/pixel) were divided into 5000 x 5000-pixel patches. Within each patch, the nucleus centers were identified by a difference of gaussians blob detection algorithm obtaining Delaunay triangulations and Voronoi diagrams providing cs radius. The ns radius was computed from stained pixels dominated by hematoxylin content with an automatic thresholding algorithm. With multiprocessing CPUs on a PC for each WSI, eight feature types were calculated preserving biopsies tissue heterogeneity: the mean and standard deviation of cs and ns distributions concatenated from all patches for cancerous and non-tumoral regions. This method was applied to 40 patients (1 WSI per patient) with treatment outcomes of post radiation-therapy (RT) recurrence (n = 9),and death (n = 8)).
Results: The WSI cancerous region cs distribution mean among patients without post-RT recurrence has a median of 6.64 microns, and those with post-RT recurrence with a median of 7.11 microns. This indicates the potential of utilizing such distribution features in treatment outcome prognosis modelling. Furthermore, at the third quartile, the WSI non-tumoral region ns distribution standard deviation among patients without post-RT recurrence has a value of 2.16 microns, and 1.46 microns for those with post-RT recurrence.
Conclusion: Our approach derives patient-specific microscopic data distributions from histopathology WSI that can be directly associated with retrospective patient outcomes. They are complementary and spatially orthogonal to information served by other medical imaging modalities such as CT, MR, and Ultrasound. Therefore, it has the unique potential to augment treatment outcome model inference when properly fused with radiological scans.
Funding Support, Disclosures, and Conflict of Interest: CIHR grant number 103548 and Canada Research Chairs Program (grant #252135)
Keywords
Image Analysis, Feature Extraction, Radiation Therapy
Taxonomy
IM/TH- Image Analysis (Single Modality or Multi-Modality): Imaging biomarkers and radiomics
Contact Email
yujing.zou@mail.mcgill.ca
Thibodeau-Antonacci, Alana; Vuong, Te; Liontis, B; Rayes, F; Pande, S; Enger, Shirin A.
Development of a Novel MRI-Compatible Applicator for Intensity Modulated Rectal Brachytherapy Proceedings Article
In: MEDICAL PHYSICS, pp. E240–E240, WILEY 111 RIVER ST, HOBOKEN 07030-5774, NJ USA 2022.
BibTeX | Tags:
@inproceedings{thibodeau2022development,
title = {Development of a Novel MRI-Compatible Applicator for Intensity Modulated Rectal Brachytherapy},
author = {Alana Thibodeau-Antonacci and Te Vuong and B Liontis and F Rayes and S Pande and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
booktitle = {MEDICAL PHYSICS},
volume = {49},
number = {6},
pages = {E240--E240},
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.
BibTeX | Tags:
@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}
}
Weishaupt, Luca L.; Vuong, Te; Thibodeau-Antonacci, Alana; Garant, A; Singh, K; Miller, C; Martin, A; Schmitt-Ulms, F; Enger, Shirin A.
PO-1325 Automated rectal tumor segmentation with inter-observer variability-based uncertainty estimates Journal Article
In: Radiotherapy and Oncology, vol. 170, pp. S1120–S1121, 2022.
BibTeX | Tags:
@article{weishaupt2022po,
title = {PO-1325 Automated rectal tumor segmentation with inter-observer variability-based uncertainty estimates},
author = {Luca L. Weishaupt and Te Vuong and Alana Thibodeau-Antonacci and A Garant and K Singh and C Miller and A Martin and F Schmitt-Ulms and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Radiotherapy and Oncology},
volume = {170},
pages = {S1120--S1121},
publisher = {Elsevier},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Zou, Yujing; Weishaupt, Luca; Enger, Shirin A.
SP-0014 McMedHacks: Deep learning for medical image analysis workshops and Hackathon in radiation oncology Journal Article
In: Radiotherapy and Oncology, vol. 170, pp. S4–S5, 2022.
Abstract | Links | BibTeX | Tags:
@article{zou2022sp,
title = {SP-0014 McMedHacks: Deep learning for medical image analysis workshops and Hackathon in radiation oncology},
author = {Yujing Zou and Luca Weishaupt and Shirin A. Enger},
url = {https://www-sciencedirect-com.proxy3.library.mcgill.ca/science/article/pii/S0167814022038695?via%3Dihub},
doi = {10.1016/S0167-8140(22)03869-5},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Radiotherapy and Oncology},
volume = {170},
pages = {S4--S5},
publisher = {Elsevier},
abstract = {Purpose/Objective: The McMedHacks workshop and presentation series was created to teach individuals from various backgrounds about deep learning (DL) for medical image analysis in May, 2021. Material/Methods: McMedHacks is a free and student-led 8-week summer program. Registration for the event was open to everyone, including a form to survey participants’ area of expertise, country of origin, level of study, and level of programming skills. The weekly workshops were instructed by 8 students and experts assisted by 20 mentors who provided weekly tutorials. Recent developments in DL and medical physics were highlighted by 21 leaders from industry and academia. A virtual grand challenge Hackathon took place at the end of the workshop series. All events were held virtually and recorded on Zoom to accommodate all time zones and locations. The workshops were designed as interactive coding demos and shared through Google Colab notebooks. Results: McMedHacks gained 356 registrations from participants of 38 different countries (Fig. 1) from undergraduates, to PhDs and MDs. A vast number of disciplines and professions were represented, dominated by medical physics students, academic, and clinical medical physicists (Fig. 2). Sixty-nine participants earned a certificate of completion by having engaged with at least 12 of all 14 events. The program received participant feedback average scores of 4.768, 4.478, 4.579, 4.292, 4.84 out of five for the qualities of presentation, workshop session, tutorial and mentor, assignments, and course delivery, respectively. The eight-week long workshop’s duration allowed participants to digest the taught materials in a continuous manner as opposed to bootcamp-style conference workshops. Conclusion: The overwhelming interest and engagement for the McMedHacks workshop series from the Radiation Oncology (RadOnc) community illustrates a demand for Artificial Intelligence (AI) education in RadOnc. The future of RadOnc clinics will inevitably integrate AI. Therefore, current RadOnc professionals, and student and resident trainees should be prepared to understand basic AI principles and its applications to troubleshoot, innovate, and collaborate. McMedHacks set an excellent example of promoting open and multidisciplinary education, scientific communication, and leadership for integrating AI education into the RadOnc community on an international level. Therefore, we advocate for implementation of AI curriculums in professional education programs such as Commission on Accreditation of Medical Physics Education Programs (CAMPEP). Furthermore, we encourage experts from around the world in the field of AI, or RadOnc, or both, to take initiatives like McMedHacks to collaborate and push forward AI education in their departments and lead practical workshops, regardless of their levels of education.},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Purpose/Objective: The McMedHacks workshop and presentation series was created to teach individuals from various backgrounds about deep learning (DL) for medical image analysis in May, 2021. Material/Methods: McMedHacks is a free and student-led 8-week summer program. Registration for the event was open to everyone, including a form to survey participants’ area of expertise, country of origin, level of study, and level of programming skills. The weekly workshops were instructed by 8 students and experts assisted by 20 mentors who provided weekly tutorials. Recent developments in DL and medical physics were highlighted by 21 leaders from industry and academia. A virtual grand challenge Hackathon took place at the end of the workshop series. All events were held virtually and recorded on Zoom to accommodate all time zones and locations. The workshops were designed as interactive coding demos and shared through Google Colab notebooks. Results: McMedHacks gained 356 registrations from participants of 38 different countries (Fig. 1) from undergraduates, to PhDs and MDs. A vast number of disciplines and professions were represented, dominated by medical physics students, academic, and clinical medical physicists (Fig. 2). Sixty-nine participants earned a certificate of completion by having engaged with at least 12 of all 14 events. The program received participant feedback average scores of 4.768, 4.478, 4.579, 4.292, 4.84 out of five for the qualities of presentation, workshop session, tutorial and mentor, assignments, and course delivery, respectively. The eight-week long workshop’s duration allowed participants to digest the taught materials in a continuous manner as opposed to bootcamp-style conference workshops. Conclusion: The overwhelming interest and engagement for the McMedHacks workshop series from the Radiation Oncology (RadOnc) community illustrates a demand for Artificial Intelligence (AI) education in RadOnc. The future of RadOnc clinics will inevitably integrate AI. Therefore, current RadOnc professionals, and student and resident trainees should be prepared to understand basic AI principles and its applications to troubleshoot, innovate, and collaborate. McMedHacks set an excellent example of promoting open and multidisciplinary education, scientific communication, and leadership for integrating AI education into the RadOnc community on an international level. Therefore, we advocate for implementation of AI curriculums in professional education programs such as Commission on Accreditation of Medical Physics Education Programs (CAMPEP). Furthermore, we encourage experts from around the world in the field of AI, or RadOnc, or both, to take initiatives like McMedHacks to collaborate and push forward AI education in their departments and lead practical workshops, regardless of their levels of education.
Enger, Shirin A.
SP-0706 Creating a'successful'work environment Journal Article
In: Radiotherapy and Oncology, vol. 170, pp. S620, 2022.
BibTeX | Tags:
@article{enger2022sp,
title = {SP-0706 Creating a'successful'work environment},
author = {Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Radiotherapy and Oncology},
volume = {170},
pages = {S620},
publisher = {Elsevier},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Weishaupt, Luca L.; Sayed, Hisham Kamal; Mao, Ximeng; Choo, Richard; Stish, Bradley J; Enger, Shirin A; Deufel, Christopher
Approaching automated applicator digitization from a new angle: Using sagittal images to improve deep learning accuracy and robustness in high-dose-rate prostate brachytherapy Journal Article
In: Brachytherapy, 2022.
BibTeX | Tags:
@article{weishaupt2022approaching,
title = {Approaching automated applicator digitization from a new angle: Using sagittal images to improve deep learning accuracy and robustness in high-dose-rate prostate brachytherapy},
author = {Luca L. Weishaupt and Hisham Kamal Sayed and Ximeng Mao and Richard Choo and Bradley J Stish and Shirin A Enger and Christopher Deufel},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Brachytherapy},
publisher = {Elsevier},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Vuong, Te; Garant, A; Khosrow-Khavar, F; Devic, S; Enger, Shirin A.; Boutros, M; Cohen, A; Miller, CS; Friedman, G; Galiatsatos, P; others,
A141 IS SURGERY STILL THE ONLY TREATMENT OPTION FOR CURABLE RECTAL CANCER? Journal Article
In: Journal of the Canadian Association of Gastroenterology, vol. 5, no. Suppl 1, pp. 13, 2022.
BibTeX | Tags:
@article{vuong2022a141,
title = {A141 IS SURGERY STILL THE ONLY TREATMENT OPTION FOR CURABLE RECTAL CANCER?},
author = {Te Vuong and A Garant and F Khosrow-Khavar and S Devic and Shirin A. Enger and M Boutros and A Cohen and CS Miller and G Friedman and P Galiatsatos and others },
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Journal of the Canadian Association of Gastroenterology},
volume = {5},
number = {Suppl 1},
pages = {13},
publisher = {Oxford University Press},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Weishaupt, Luca L; Vuong, Te; Thibodeau-Antonacci, Alana; Garant, A; Singh, KS; Miller, C; Martin, A; Enger, Shirin A.
A121 QUANTIFYING INTER-OBSERVER VARIABILITY IN THE SEGMENTATION OF RECTAL TUMORS IN ENDOSCOPY IMAGES AND ITS EFFECTS ON DEEP LEARNING Journal Article
In: Journal of the Canadian Association of Gastroenterology, vol. 5, no. Supplement_1, pp. 140–142, 2022.
BibTeX | Tags:
@article{weishaupt2022a121,
title = {A121 QUANTIFYING INTER-OBSERVER VARIABILITY IN THE SEGMENTATION OF RECTAL TUMORS IN ENDOSCOPY IMAGES AND ITS EFFECTS ON DEEP LEARNING},
author = {Luca L Weishaupt and Te Vuong and Alana Thibodeau-Antonacci and A Garant and KS Singh and C Miller and A Martin and Shirin A. Enger},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Journal of the Canadian Association of Gastroenterology},
volume = {5},
number = {Supplement_1},
pages = {140--142},
publisher = {Oxford University Press US},
keywords = {},
pubstate = {published},
tppubtype = {article}
}
Rodrigues-Machado, Fernanda C; Pestre, Pauline; Dumont, Vincent; Bernard, Simon; Janitz, Erika; Scanlon, Liam; Enger, Shirin A.; Childress, Lilian; Sankey, Jack
Sideband cavity absorption readout (SideCAR) with a robust frequency lock Journal Article
In: Optics Express, vol. 30, no. 2, pp. 754–767, 2022.
BibTeX | Tags:
@article{rodrigues2022sideband,
title = {Sideband cavity absorption readout (SideCAR) with a robust frequency lock},
author = {Fernanda C Rodrigues-Machado and Pauline Pestre and Vincent Dumont and Simon Bernard and Erika Janitz and Liam Scanlon and Shirin A. Enger and Lilian Childress and Jack Sankey},
year = {2022},
date = {2022-01-01},
urldate = {2022-01-01},
journal = {Optics Express},
volume = {30},
number = {2},
pages = {754--767},
publisher = {Optical Society of America},
keywords = {},
pubstate = {published},
tppubtype = {article}
}