1
|
Georgi P, Tedgren ÅC, Persson L, Johansen JG. Determination of intrinsic energy dependence of point-like inorganic scintillation detector in brachytherapy. Med Phys 2024; 51:5059-5069. [PMID: 38197459 DOI: 10.1002/mp.16904] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Revised: 12/05/2023] [Accepted: 12/07/2023] [Indexed: 01/11/2024] Open
Abstract
BACKGROUND Inorganic scintillation detectors (ISDs) are promising for in vivo dosimetry in brachytherapy (BT). ISDs have fast response, providing time resolved dose rate information, and high sensitivity, attributed to high atomic numbers. However, the conversion of the detector signal to absorbed dose-to-water is highly dependent on the energy spectrum of the incident radiation. This dependence is comprised of absorbed dose energy dependence, obtainable with Monte Carlo (MC) simulation, and the absorbed dose-to-signal conversion efficiency or intrinsic energy dependence requiring measurements. Studies have indicated negligible intrinsic energy dependence of ZnSe:O-based ISDs in Ir-192 BT. A full characterization has not been performed earlier. PURPOSE This study characterizes the intrinsic energy dependence of ZnSe:O-based ISDs for kV X-ray radiation qualities, with energies relevant for BT. METHODS Three point-like ISDs made from fiber-coupled cuboid ZnSe:O-based scintillators were calibrated at the Swedish National Metrology Laboratory for ionizing radiation. The calibration was done in terms of air kerma free-in-air,K air ${K}_{{\mathrm{air}}}$ , in 13 X-ray radiation qualities, Q $Q$ , from 25 to 300 kVp (CCRI 25-250 kV and ISO 4037 N-series), and in terms of absorbed dose to water,D w ${D}_{\mathrm{w}}$ , in a Co-60 beam,Q 0 ${Q}_0$ . The mean absorbed dose to the ISDs, relative toK air ${K}_{{\mathrm{air}}}$ andD w ${D}_{\mathrm{w}}$ , were obtained with the MC code TOPAS (Geant4) using X-ray spectra obtained with SpekPy software and laboratory filtration data and a generic Co-60 source. The intrinsic energy dependence was determined as a function of effective photon energy,E e f f ${E}_{eff}$ , (relative to Co-60). The angular dependence of the ISD signal was measured in a 25 kVp (0.20 mm Al HVL) and 135 kVp beam (0.48 mm Cu HVL), by rotating the ISDs 180° around the fiber's longitudinal axis (perpendicular to the beam). A full 360° was not performed due to setup limitations. The impact of detector design was quantified with MC simulation. RESULTS Above 30 keVE e f f ${E}_{eff}$ the intrinsic energy dependence varied with less than 5 ± 4% from unity for all detectors (with the uncertainty expressed as the mean of all expanded measurement uncertainties for individualE e f f ${E}_{eff}$ above 30 keV, k = 2). Below 30 keV, it decreased with up to 17% and inter-detector variations of 13% were observed, likely due to differences in detector geometry not captured by the simulations using nominal geometry. In the 25 kVp radiation quality, the ISD signal varied with 24% over a ∼45° rotation. For 135 kVp, the corresponding variation was below 3%. Assuming a 0.05 mm thicker layer of reflective paint around the sensitive volume changed the absorbed dose with 6.3% at the lowestE e f f ${E}_{eff}$ , and with less than 2% at higher energies. CONCLUSION The study suggests that the ISDs have an intrinsic energy dependence relative to Co-60 lower than 5 ± 4% in radiation qualities withE e f f ${E}_{eff}\ $ > 30 keV. Therefore, they could in principle be calibrated in a Co-60 beam quality and transferred to such radiation qualities with correction factors determined only by the absorbed dose energy dependence obtained from MC simulations. This encourages exploration of the ISDs' applications in intensity modulated BT with Yb-169 or other novel intermediate energy isotopes.
Collapse
Affiliation(s)
- Peter Georgi
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Åsa Carlsson Tedgren
- Department of Medicine, Health, and Caring Sciences, Linköping University, Linköping, Sweden
- Department of Oncology-Pathology, the Karolinska Institute, Stockholm, Sweden
- Department of Medical Physics and Nuclear Medicine, Karolinska Hospital, Stockholm, Sweden
| | - Linda Persson
- Department of Oncology-Pathology, the Karolinska Institute, Stockholm, Sweden
- National Regulation Department, Swedish Radiation Safety Authority, Stockholm, Sweden
| | - Jacob Graversen Johansen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| |
Collapse
|
2
|
Georgi PD, Jørgensen EB, Heidotting M, Tanderup K, Kertzscher G, Johansen JG. A simple calibration routine for small inorganic scintillation detectors for in vivo dosimetry during brachytherapy. Brachytherapy 2024:S1538-4721(24)00070-9. [PMID: 38853063 DOI: 10.1016/j.brachy.2024.05.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 04/04/2024] [Accepted: 05/07/2024] [Indexed: 06/11/2024]
Abstract
BACKGROUND In vivo dosimetry (IVD) is rarely performed in brachytherapy (BT), allowing potential dose misadministration to go unnoticed. This study presents a clinical routine-calibration method of detectors for IVD in high (HDR) and pulsed dose rate (PDR) Ir-192 BT. PURPOSE To evaluate the dosimetric precision and feasibility of an in-clinic calibration routine of detectors for IVD in afterloading BT. METHODS Calibrations were performed in a PMMA phantom with two needles inserted 20 mm apart. The source was loaded in one of the needles at 15 dwells for 10 s. The detector was placed in the other needle, and its signal was recorded. The mean signal at each dwell position was fitted to the expected dose rate with the calibration factor and the detector's longitudinal position being free parameters. The method was tested with an inorganic scintillation detector using one Ir-192 FlexiSource HDR and two Ir-192 GammaMedPlus PDR sources and followed by validation measurements in water. RESULTS The standard measurement uncertainty (k = 1) of the calibration factor in absolute terms (Gy/s) was 3.2/3.4% for the HDR/PDR source. The uncertainty was dominated by source strength uncertainty, and the precision of the method was <1%. The mean ± 1SD of the difference in measured and expected dose rate during validation was 1.5 ± 4.7% (HDR) and 0.0 ± 4.1% (PDR) with a positional uncertainty in the setup of 0.33/0.23 mm (HDR/PDR) (k = 1). CONCLUSION A precise and feasible in-clinic calibration method for IVD and source strength consistency tests in BT was presented.
Collapse
Affiliation(s)
- Peter D Georgi
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark.
| | - Erik B Jørgensen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | | | - Kari Tanderup
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Jacob G Johansen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| |
Collapse
|
3
|
Tho D, Bélanger C, Jørgensen EB, Tanguay J, Rosales HML, Beddar S, Johansen JG, Kertzscher G, Lavallée MC, Beaulieu L. Establishing a fingerprinting method for fast catheter identification in HDR brachytherapy in vivo dosimetry. Brachytherapy 2024; 23:165-172. [PMID: 38281894 DOI: 10.1016/j.brachy.2023.10.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Revised: 08/23/2023] [Accepted: 10/06/2023] [Indexed: 01/30/2024]
Abstract
PURPOSE To use quantities measurable during in vivo dosimetry to build unique channel identifiers, that enable detection of brachytherapy errors. MATERIALS AND METHODS Treatment plan of 360 patients with prostate cancer who underwent high-dose-rate brachytherapy (range, 16-25 catheters; mean, 17) were used. A single point virtual dosimeter was placed at multiple positions within the treatment geometry, and the source-dosimeter distance and dwell time were determined for each dwell position in each catheter. These values were compared across all catheters, dwell position by dwell position, simulating a treatment delivery. A catheter was considered uniquely identified if, for a given dwell position, no other catheters had the same measured values. The minimum number of dwell positions needed to identify a specific catheter and the optimal dosimeter location uniquely were determined. The radial (r) and vertical (z) dimensions of the source-dosimeter distance were also examined for their utility in discriminating catheters. RESULTS Using a virtual dosimeter with no uncertainties, all catheters were identified in 359 of the 360 cases with 9 dwell position measurements. When only the dwell time were measured, all catheters were uniquely identified after 1 dwell position. With a 2-mm spatial accuracy (r,z), all catheters were identified in 94% of the plans. Simultaneous measurement of source-dosimeter distance and dwell time ensured full catheter identification in all plans ranging from 2 to 6 dwell positions. The number of dwell positions needed to uniquely identify all catheters was lower when the distance from the implant center was higher. CONCLUSIONS The most efficient fingerprinting approach involved combining source-dosimeter distance (i.e., source tracking) and dwell time. The further the dosimeter is placed from the center of the implant the better it can uniquely identify catheters.
Collapse
Affiliation(s)
- Daline Tho
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada.
| | - Cédric Bélanger
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada
| | - Erik B Jørgensen
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Jérémie Tanguay
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada
| | - Haydee M L Rosales
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada
| | - Sam Beddar
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX
| | - Jacob G Johansen
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Marie-Claude Lavallée
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada
| | - Luc Beaulieu
- Centre Intégré De Cancérologie, CHU De Québec, Université Laval, Centre De Recherche Chu De Québec, Québec, Canada; Département De Physique, De Génie Physique Et D'optique, Centre De Recherche Sur Le Cancer, Québec, Canada
| |
Collapse
|
4
|
Kaveckyte V, Jørgensen EB, Kertzscher G, Johansen JG, Tedgren ÅC. Monte Carlo characterization of high atomic number inorganic scintillators for in vivo dosimetry in 192 Ir brachytherapy. Med Phys 2022; 49:4715-4730. [PMID: 35443079 DOI: 10.1002/mp.15674] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 03/01/2022] [Accepted: 04/06/2022] [Indexed: 11/05/2022] Open
Abstract
BACKGROUND There is increased interest in vivo dosimetry for 192 Ir brachytherapy (BT) treatments using high atomic number (Z) inorganic scintillators. Their high light output enables construction of small detectors with negligible stem effect and simple readout electronics. Experimental determination of absorbed-dose energy dependence of detectors relative to water is prevalent, but it can be prone to high detector positioning uncertainties and does not allow for decoupling of absorbed-dose energy dependence from other factors affecting detector response. PURPOSE To investigate which measurement conditions and detector properties could affect their absorbed-dose energy dependence in BT in vivo dosimetry. METHODS We used a general-purpose MC code penelope for the characterization of high-Z inorganic scintillators with the focus on ZnSe (Z¯=32). Two other promising media CsI (Z¯=54) and Al2 O3 (Z¯=11) were included for comparison in selected scenarios. We determined absorbed-dose energy dependence of crystals relative to water under different scatter conditions (calibration phantom 12 × 12 × 30 cm3 , characterization phantoms 20 × 20 × 20 cm3 , 30 × 30 × 30 cm3 , 40 × 40 × 40 cm3 , and patient-like elliptic phantom 40 × 30 × 25 cm3 ). To mimic irradiation conditions during prostate treatments, we evaluated whether the presence of pelvic bones and calcifications affect ZnSe response. ZnSe detector design influence was also investigated. RESULTS In contrast to low-Z organic and medium-Z inorganic scintillators, ZnSe and CsI media have substantially greater absorbed-dose energy dependence relative to water. The response was phantom-size dependent and changed by 11 % between limited- and full-scatter conditions for ZnSe, but not for Al2 O3 . For a given phantom size, a part of the absorbed-dose energy dependence of ZnSe is caused not due to in-phantom scatter but due to source anisotropy. Thus, the absorbed-dose energy dependence of high-Z scintillators is a function of not only the radial distance but also the polar angle. Pelvic bones did not affect ZnSe response, whereas large and intermediate size calcifications reduced it by 9 % and 5 %, respectively, when placed midway between the source and the detector. CONCLUSIONS Unlike currently prevalent low- and medium-Z scintillators, high-Z crystals are sensitive to characterization and in vivo measurement conditions. However, good agreement between MC data for ZnSe in the present study and experimental data for ZnSe:O by Jørgensen et al (2021) suggest that detector signal is proportional to the average absorbed dose to the detector cavity. This enables an easy correction for non-TG43-like scenarios (e.g., patient sizes and calcifications) through MC simulations. Information that should be provided to the clinic by the detector vendors. This article is protected by copyright. All rights reserved.
Collapse
Affiliation(s)
- Vaiva Kaveckyte
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, SE-581 85, Sweden
| | - Erik B Jørgensen
- Department of Clinical Medicine, Aarhus University, Aarhus, DK-8000, Denmark.,Department of Oncology, Aarhus University Hospital, Aarhus, DK-8000, Denmark
| | - Gustavo Kertzscher
- Department of Oncology, Aarhus University Hospital, Aarhus, DK-8000, Denmark
| | - Jacob G Johansen
- Department of Clinical Medicine, Aarhus University, Aarhus, DK-8000, Denmark.,Department of Oncology, Aarhus University Hospital, Aarhus, DK-8000, Denmark
| | - Åsa Carlsson Tedgren
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, SE-581 85, Sweden.,Department of Medical Radiation Physics and Nuclear Medicine, Karolinska University Hospital, Stockholm, SE-171 76, Sweden.,Department of Oncology-Pathology, Karolinska Institute, Stockholm, SE-171 76, Sweden
| |
Collapse
|
5
|
Lekatou A, Peppa V, Karaiskos P, Pantelis E, Papagiannis P. On the potential of 2D ion chamber arrays for high-dose rate remote afterloading brachytherapy quality assurance. Phys Med Biol 2022; 67. [PMID: 35334474 DOI: 10.1088/1361-6560/ac612d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 03/25/2022] [Indexed: 01/08/2023]
Abstract
Objective. To investigate the potential of 2D ion chamber arrays to serve as a standalone tool for the verification of source strength, positioning and dwell time, within the framework of192Ir high-dose rate brachytherapy device quality assurance (QA).Approach.A commercially available ion chamber array was used. Fitting of a 2D Lorentzian peak function to experimental data from a multiple source dwell position irradiation on a frame-by-frame basis, facilitated tracking of the source center orthogonal projection on the array plane. For source air kerma strength verification, Monte Carlo simulation was employed to obtain a chamber array- and source-specific correction factor of calibration with a 6 MV photon beam. This factor converted the signal measured by each ion chamber element to air kerma in free space. A source positioning correction was also applied to lift potential geometry mismatch between experiment and Monte Carlo simulation.Main results.Spatial and temporal accuracy of source movement was verified within 0.5 mm and 0.02 s, respectively, in compliance with the test endpoints recommended by international professional societies. The source air kerma strength was verified experimentally within method uncertainties estimated as 1.44% (k = 1). The source positioning correction method employed did not introduce bias to experimental results of irradiations where source positioning was accurate. Development of a custom jig attachable to the chamber array for accurate and reproducible experimental set up would improve testing accuracy and obviate the need for source positioning correction in air kerma strength verification.Significance.Delivery of a single irradiation plan, optimized based on results of this work, to a 2D ion chamber array can be used for concurrent testing of source position, dwell time and air kerma strength, and the procedure can be expedited through automation. Chamber arrays merit further study in treatment planning QA and real time,in vivodose verification.
Collapse
Affiliation(s)
- Aristea Lekatou
- Medical Physics Laboratory, Department of Medicine, National and Kapodistrian University of Athens, 75 Mikras Asias, 11527 Athens, Greece
| | - Vasiliki Peppa
- Medical Physics Laboratory, Department of Medicine, National and Kapodistrian University of Athens, 75 Mikras Asias, 11527 Athens, Greece
| | - Pantelis Karaiskos
- Medical Physics Laboratory, Department of Medicine, National and Kapodistrian University of Athens, 75 Mikras Asias, 11527 Athens, Greece
| | - Evangelos Pantelis
- Medical Physics Laboratory, Department of Medicine, National and Kapodistrian University of Athens, 75 Mikras Asias, 11527 Athens, Greece
| | - Panagiotis Papagiannis
- Medical Physics Laboratory, Department of Medicine, National and Kapodistrian University of Athens, 75 Mikras Asias, 11527 Athens, Greece
| |
Collapse
|
6
|
Jørgensen EB, Buus S, Bentzen L, Hokland SB, Rylander S, Kertzscher G, Beddar S, Tanderup K, Johansen JG. 3D dose reconstruction based on in vivo dosimetry for determining the dosimetric impact of geometric variations in high-dose-rate prostate brachytherapy. Radiother Oncol 2022; 171:62-68. [PMID: 35033604 DOI: 10.1016/j.radonc.2022.01.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2021] [Revised: 12/30/2021] [Accepted: 01/04/2022] [Indexed: 11/16/2022]
Abstract
INTRODUCTION In vivo dosimetry (IVD) can be used for source tracking (ST), i.e., estimating source positions, during brachytherapy. The aim of this study was to exploit IVD-based ST to perform 3D dose reconstruction for high-dose-rate prostate brachytherapy and to evaluate the robustness of the treatments against observed geometric variations. MATERIALS AND METHODS Twenty-three fractions of high-dose-rate prostate brachytherapy were analysed. The treatment planning was based on MRI. Time-resolved IVD was performed using a fibre-coupled scintillator. ST was retrospectively performed using the IVD measurements. The ST identified 2D positional shifts of each treatment catheter and thereby inferred updated source positions. For each fraction, the dose was recalculated based on the source-tracked catheter positions and compared with the original plan dose using differences in dose volume histogram indices. RESULTS Of 352 treatment catheters, 344 had shifts of less than 5 mm. Shifts between 5 and 10 mm were observed for 3 catheters, and shifts greater than 10 mm for 2 catheters. The ST failed for 3 catheters. The maximum relative difference in clinical target volume (prostate + 3 mm isotropic margin) D90% was 5%. In one fraction, the bladder D2cm3 dose increased by 18% (1.4Gy) due to a single source position being inside the bladder rather than nearby as planned. The max increase in urethra dose was 1.5Gy (15%). CONCLUSION IVD-based 3D dose reconstruction for high-dose-rate prostate brachytherapy is feasible. The dosimetric impact of the observed catheter shifts was limited. Dose reconstruction can therefore aid in determining the dosimetric impact of geometric variations and errors in brachytherapy.
Collapse
Affiliation(s)
- Erik B Jørgensen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark.
| | - Simon Buus
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Lise Bentzen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Susanne Rylander
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Sam Beddar
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States
| | - Kari Tanderup
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Jacob G Johansen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| |
Collapse
|
7
|
Linares Rosales HM, Couture G, Archambault L, Beddar S, Després P, Beaulieu L. On the use of machine learning methods for mPSD calibration in HDR brachytherapy. Phys Med 2021; 91:73-79. [PMID: 34717139 DOI: 10.1016/j.ejmp.2021.10.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Revised: 07/15/2021] [Accepted: 10/02/2021] [Indexed: 11/17/2022] Open
Abstract
We sought to evaluate the feasibility of using machine learning (ML) algorithms for multipoint plastic scintillator detector (mPSD) calibration in high-dose-rate (HDR) brachytherapy. Dose measurements were conducted under HDR brachytherapy conditions. The dosimetry system consisted of an optimized 1-mm-core mPSD and a compact assembly of photomultiplier tubes coupled with dichroic mirrors and filters. An 192Ir source was remotely controlled and sent to various positions in a homemade PMMA holder, ensuring 0.1-mm positional accuracy. Dose measurements covering a range of 0.5 to 12 cm of source displacement were carried out according to TG-43 U1 recommendations. Individual scintillator doses were decoupled using a linear regression model, a random forest estimator, and artificial neural network algorithms. The dose predicted by the TG-43U1 formalism was used as the reference for system calibration and ML algorithm training. The performance of the different algorithms was evaluated using different sample sizes and distances to the source for the mPSD system calibration. We found that the calibration conditions influenced the accuracy in predicting the measured dose. The decoupling methods' deviations from the expected TG-43 U1 dose generally remained below 20%. However, the dose prediction with the three algorithms was accurate to within 7% relative to the dose predicted by the TG-43 U1 formalism when measurements were performed in the same range of distances used for calibration. In such cases, the predictions with random forest exhibited minimal deviations (<2%). However, the performance random forest was compromised when the predictions were done beyond the range of distances used for calibration. Because the linear regression algorithm can extrapolate the data, the dose prediction by the linear regression was less influenced by the calibration conditions than random forest. The linear regression algorithm's behavior along the distances to the source was smoother than those for the random forest and neural network algorithms, but the observed deviations were more significant than those for the neural network and random forest algorithms. The number of available measurements for training purposes influenced the random forest and neural network models the most. Their accuracy tended to converge toward deviation values close to 1% from a number of dwell positions greater than 100. In performing HDR brachytherapy dose measurements with an optimized mPSD system, ML algorithms are good alternatives for precise dose reporting and treatment assessment during this kind of cancer treatment.
Collapse
Affiliation(s)
- Haydee M Linares Rosales
- Département de physique, de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec, Canada; Département de radio-oncologie et Axe Oncologie du CRCHU de Québec, CHU de Québec - Université Laval, QC, Canada.
| | - Gabriel Couture
- Département de physique, de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec, Canada; Département de radio-oncologie et Axe Oncologie du CRCHU de Québec, CHU de Québec - Université Laval, QC, Canada
| | - Louis Archambault
- Département de physique, de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec, Canada; Département de radio-oncologie et Axe Oncologie du CRCHU de Québec, CHU de Québec - Université Laval, QC, Canada
| | - Sam Beddar
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, United States; The University of Texas MD Anderson UTHealth Graduate School of Biomedical Sciences, Houston, TX, United States
| | - Philippe Després
- Département de physique, de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec, Canada; Département de radio-oncologie et Axe Oncologie du CRCHU de Québec, CHU de Québec - Université Laval, QC, Canada
| | - Luc Beaulieu
- Département de physique, de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec, Canada; Département de radio-oncologie et Axe Oncologie du CRCHU de Québec, CHU de Québec - Université Laval, QC, Canada
| |
Collapse
|
8
|
Jørgensen EB, Johansen JG, Overgaard J, Piché-Meunier D, Tho D, Rosales HML, Tanderup K, Beaulieu L, Kertzscher G, Beddar S. A high-Z inorganic scintillator-based detector for time-resolved in vivo dosimetry during brachytherapy. Med Phys 2021; 48:7382-7398. [PMID: 34586641 DOI: 10.1002/mp.15257] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Revised: 08/08/2021] [Accepted: 09/09/2021] [Indexed: 11/11/2022] Open
Abstract
PURPOSE High-dose rate (HDR) and pulsed-dose rate (PDR) brachytherapy would benefit from an independent treatment verification system to monitor treatment delivery and to detect errors in real time. This paper characterizes and provides an uncertainty budget for a detector based on a fiber-coupled high-Z inorganic scintillator capable of performing time-resolved in vivo dosimetry during HDR and PDR brachytherapy. METHOD The detector was composed of a detector probe and an optical reader. The detector probe consisted of either a 0.5 × 0.4 × 0.4 mm3 (HDR) or a 1.0 × 0.4 × 0.4 mm3 (PDR) cuboid ZnSe:O crystal glued onto an optical-fiber cable. The outer diameter of the detector probes was 1 mm, and fit inside standard brachytherapy catheters. The signal from the detector probe was read out at 20 Hz by a photodiode and a data acquisition device inside the optical reader. In order to construct an uncertainty budget for the detector, six characteristics were determined: (1) temperature dependence of the detector probe, (2) energy dependence as a function of the probe-to-source position in 2D (determined with 2 mm resolution using a robotic arm), (3) the signal-to-noise ratio (SNR), (4) short-term stability over 8 h, and (5) long-term stability of three optical readers and four probes used for in vivo monitoring in HDR and PDR treatments over 21 months (196 treatments and 189 detector calibrations, and (6) dose-rate dependence. RESULTS The total uncertainty of the detector at a 20 mm probe-to-source distance was < 5.1% and < 5.8% for the HDR and PDR versions, respectively. Regarding the above characteristics, (1) the sensitivity of the detector decreased by an average of 1.4%/°C for detector probe temperatures varying from 22 to 37°C; (2) the energy dependence of the detector was nonlinear and depended on both probe-to-source distance and the angle between the probe and the brachytherapy source; (3) the median SNRs were 187 and 34 at a 20 mm probe-to-source distance for the HDR and PDR versions, respectively (corresponding median source activities of 4.8 and 0.56 Ci, respectively); (4) the detector response varied by 0.6% in 11 identical irradiations over 8 h; (5) the sensitivity of the four detector probes decreased systematically by 0-1.2%/100 Gy of dose delivered to the probes, and random fluctuations of 4.8% in the sensitivity were observed for the three probes used in PDR and 1.9% for the probe used in HDR; and (6) the detector response was linear with dose rate. CONCLUSION ZnSe:O detectors can be used effectively for in vivo dosimetry and with high accuracy for HDR and PDR brachytherapy applications.
Collapse
Affiliation(s)
- Erik B Jørgensen
- Health Graduate School, Aarhus University, Aarhus, Denmark.,Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Jacob G Johansen
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Dominique Piché-Meunier
- Département de physique-de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec City, Quebec, Canada.,Département de radio-oncologie et Axe Oncologie, CHU de Québec-Université Laval, Québec City, Quebec, Canada
| | - Daline Tho
- Département de physique-de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec City, Quebec, Canada.,Département de radio-oncologie et Axe Oncologie, CHU de Québec-Université Laval, Québec City, Quebec, Canada
| | - Haydee M L Rosales
- Département de physique-de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec City, Quebec, Canada.,Département de radio-oncologie et Axe Oncologie, CHU de Québec-Université Laval, Québec City, Quebec, Canada
| | - Kari Tanderup
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Luc Beaulieu
- Département de physique-de génie physique et d'optique et Centre de recherche sur le cancer, Université Laval, Québec City, Quebec, Canada.,Département de radio-oncologie et Axe Oncologie, CHU de Québec-Université Laval, Québec City, Quebec, Canada
| | - Gustavo Kertzscher
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark.,Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Sam Beddar
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| |
Collapse
|
9
|
Quantifying clinical severity of physics errors in high-dose rate prostate brachytherapy using simulations. Brachytherapy 2021; 20:1062-1069. [PMID: 34193362 DOI: 10.1016/j.brachy.2021.05.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 05/10/2021] [Accepted: 05/17/2021] [Indexed: 11/22/2022]
Abstract
PURPOSE To quantitatively evaluate through automated simulations the clinical significance of potential high-dose rate (HDR) prostate brachytherapy (HDRPB) physics errors selected from our internal failure-modes and effect analysis (FMEA). METHODS AND MATERIALS A list of failure modes was compiled and scored independently by 8 brachytherapy physicists on a one-to-ten scale for severity (S), occurrence (O), and detectability (D), with risk priority number (RPN) = SxOxD. Variability of RPNs across observers (standard deviation/average) was calculated. Six idealized HDRPB plans were generated, and error simulations were performed: single (N = 1722) and systematic (N = 126) catheter shifts (craniocaudal; -1cm:1 cm); single catheter digitization errors (tip and connector needle-tips displaced independently in random directions; 0.1 cm:0.5 cm; N = 44,318); and swaps (two catheters swapped during digitization or connection; N = 528). The deviations due to each error in prostate D90%, urethra D20%, and rectum D1cm3 were analyzed using two thresholds: 5-20% (possible clinical impact) and >20% (potentially reportable events). RESULTS Twenty-nine relevant failure modes were described. Overall, RPNs ranged from 6 to 108 (average ± 1 standard deviation, 46 ± 23), with responder variability ranging from 19% to 184% (average 75% ± 30%). Potentially reportable events were observed in the simulations for systematic shifts >0.4 cm for prostate and digitization errors >0.3 cm for the urethra and >0.4 cm for rectum. Possible clinical impact was observed for catheter swaps (all organs), systematic shifts >0.2 cm for prostate and >0.4 cm for rectum, and digitization errors >0.2 cm for prostate and >0.1 cm for urethra and rectum. CONCLUSIONS A high variability in RPN scores was observed. Systematic simulations can provide insight in the severity scoring of multiple failure modes, supplementing typical FMEA approaches.
Collapse
|
10
|
Fonseca GP, van Wagenberg T, Voncken R, Podesta M, van Beveren C, van Limbergen E, Lutgens L, Vanneste B, Berbee M, Reniers B, Verhaegen F. Brachytherapy treatment verification using gamma radiation from the internal treatment source combined with an imaging panel-a phantom study. Phys Med Biol 2021; 66. [PMID: 33831856 DOI: 10.1088/1361-6560/abf605] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2021] [Accepted: 04/08/2021] [Indexed: 12/15/2022]
Abstract
Brachytherapy has an excellent clinical outcome for different treatment sites. However,in vivotreatment verification is not performed in the majority of hospitals due to the lack of proper monitoring systems. This study investigates the use of an imaging panel (IP) and the photons emitted by a high dose rate (HDR)192Ir source to track source motion and obtain some information related to the patient anatomy. The feasibility of this approach was studied by monitoring the treatment delivery to a 3D printed phantom that mimicks a prostate patient. A 3D printed phantom was designed with a template for needle insertion, a cavity ('rectum') to insert an ultrasound probe, and lateral cavities used to place tissue-equivalent materials. CT images were acquired to create HDR192Ir treatment plans with a range of dwell times, interdwell distances and needle arrangements. Treatment delivery was verified with an IP placed at several positions around the phantom using radiopaque markers on the outer surface to register acquired IP images with the planning CT. All dwell positions were identified using acquisition times ≤0.11 s (frame rates ≥ 9 fps). Interdwell distances and dwell positions (in relation to the IP) were verified with accuracy better than 0.1 cm. Radiopaque markers were visible in the acquired images and could be used for registration with CT images. Uncertainties for image registration (IP and planning CT) between 0.1 and 0.4 cm. The IP is sensitive to tissue-mimicking insert composition and showed phantom boundaries that could be used to improve treatment verification. The IP provided sufficient time and spatial resolution for real-time source tracking and allows for the registration of the planning CT and IP images. The results obtained in this study indicate that several treatment errors could be detected including swapped catheters, incorrect dwell times and dwell positions.
Collapse
Affiliation(s)
- G P Fonseca
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - T van Wagenberg
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - R Voncken
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - M Podesta
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - C van Beveren
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - E van Limbergen
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - L Lutgens
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - B Vanneste
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - M Berbee
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| | - B Reniers
- Research group NuTeC, Centre for Environmental Sciences, Hasselt University, Diepenbeek, Belgium
| | - F Verhaegen
- Department of Radiation Oncology (MAASTRO), GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre+, Doctor Tanslaan 12, 6229 ET Maastricht, The Netherlands
| |
Collapse
|
11
|
Linares Rosales HM, Johansen JG, Kertzscher G, Tanderup K, Beaulieu L, Beddar S. 3D source tracking and error detection in HDR using two independent scintillator dosimetry systems. Med Phys 2021; 48:2095-2107. [PMID: 33222208 DOI: 10.1002/mp.14607] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2020] [Revised: 09/22/2020] [Accepted: 11/01/2020] [Indexed: 01/22/2023] Open
Abstract
PURPOSE The aim of this study is to perform three-dimensional (3D) source position reconstruction by combining in vivo dosimetry measurements from two independent detector systems. METHODS Time-resolved dosimetry was performed in a water phantom during HDR brachytherapy irradiation with 192 Ir source using two detector systems. The first was based on three plastic scintillator detectors and the second on a single inorganic crystal (CsI:Tl). Brachytherapy treatments were simulated in water under TG-43U1 conditions, including a HDR prostate plan. Treatment needles were placed in distances covering a range of source movement of 120 mm around the detectors. The distance from each dwell position to each scintillator was determined based on the measured dose rates. The three distances given by the mPSD were recalculated to a position along the catheter (z) and a distance radially away from the mPSD (xy) for each dwell position (a circumference around the mPSD). The source x, y, and z coordinates were derived from the intersection of the mPSD's circumference with the sphere around the ISD based on the distance to this detector. We evaluated the accuracy of the source position reconstruction as a function of the distance to the source, the most likely location for detector positioning within a prostate volume, as well as the capacity to detect positioning errors. RESULTS Approximately 4000 source dwell positions were tracked for eight different HDR plans. An intersection of the mPSD torus and the ISD sphere was observed in 77.2% of the dwell positions, assuming no uncertainty in the dose rate determined distance. This increased to 100% if 1σ search regions were added. However, only 73(96)% of the expected dwell positions were found within the intersection band for 1(2) σ uncertainties. The agreement between the source's reconstructed and expected positions was within 3 mm for a range of distances to the source up to 50 mm. The experiments on a HDR prostate plan, showed that by having at least one of the detectors located in the middle of the prostate volume, reduces the measurement deviations considerably compared to scenarios where the detectors were located outside of the prostate volume. The analysis showed a detection probability that, in most cases, is far from the random detection threshold. Errors of 1(2) mm can be detected in ranges of 5-25 (25-50) mm from the source, with a true detection probability rate higher than 80%, while the false probability rate is kept below 20%. CONCLUSIONS By combining two detector responses, we enabled the determination of the absolute source coordinates. The combination of the mPSD and the ISD in vivo dosimetry constitutes a promising alternative for real-time 3D source tracking in HDR brachytherapy.
Collapse
Affiliation(s)
| | - Jacob G Johansen
- Department of Oncology, Aarhus University Hospital, Aarhus C, Denmark
| | | | - Kari Tanderup
- Department of Oncology, Aarhus University Hospital, Aarhus C, Denmark
| | - Luc Beaulieu
- CHU de Quebec-Université Laval, Quebec, Canada.,Université Laval, Quebec, Canada
| | - Sam Beddar
- The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| |
Collapse
|
12
|
Jørgensen EB, Kertzscher G, Buus S, Bentzen L, Hokland SB, Rylander S, Tanderup K, Johansen JG. Accuracy of an in vivo dosimetry-based source tracking method for afterloading brachytherapy - A phantom study. Med Phys 2021; 48:2614-2623. [PMID: 33655555 DOI: 10.1002/mp.14812] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 01/07/2021] [Accepted: 02/04/2021] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To report on the accuracy of an in vivo dosimetry (IVD)-based source tracking (ST) method for high dose rate (HDR) prostate brachytherapy (BT). METHODS The ST was performed on a needle-by-needle basis. A least square fit of the expected to the measured dose rate was performed using the active dwell positions in the given needle. Two fitting parameters were used to determine the position of each needle relative to the IVD detector: radial (away or toward the detector) and longitudinal (along the axis of the treatment needle). The accuracy of the ST was assessed in a phantom where the geometries of five HDR prostate BT treatments previously treated at our clinic were reproduced. For each of the five treatment geometries, one irradiation was performed with the detector placed in the middle of the implant. Furthermore, four additional irradiations were performed for one of the geometries where the detector was retracted caudally in four steps of 10-15 mm and up to 12 mm inferior of the most inferior active dwell position, which represented the prostate apex. The time resolved dose measurements were retrieved at a rate of 20 Hz using a detector based on an Al2 O3 :C radioluminescence crystal, which was placed inside a standard BT needle. Individual calibrations of the detector were performed prior to each of the nine irradiations. RESULTS Source tracking could be applied in all needles across all nine irradiations. For irradiations with the detector located in the middle region of the implant (a total of 89 needles), the mean ± standard deviation (SD, k = 1) accuracy was -0.01 mm ± 0.38 mm and 0.30 mm ± 0.38 mm in the radial and longitudinal directions, respectively. Caudal retraction of the detector did not lead to reduced accuracy as long as the detector was located superior to the most inferior active dwell positions in all needles. However, reduced accuracy was observed for detector positions inferior to the most inferior active dwell positions which corresponded to detector positions in and inferior to the prostate apex region. Detector positions in the prostate apex and 12 mm inferior to the prostate resulted in mean ± SD (k = 1) ST accuracy of 0.7 mm ± 1 mm and 2.8 mm ± 1.6 mm, respectively, in radial direction, and -1.7 mm ± 1 mm and -2.1 mm ± 1.1 mm, respectively, in longitudinal direction. The largest deviations for the configurations with those detector positions were 2.6 and 5.4 mm, respectively, in the radial direction and -3.5 and -3.8 mm, respectively, in the longitudinal direction. CONCLUSION This phantom study demonstrates that ST based on IVD during prostate BT is adequately accurate for clinical use. The ST yields submillimeter accuracy on needle positions as long as the IVD detector is positioned superior to at least one active dwell position in all needles. Locations of the detector inferior to the prostate apex result in decreased ST accuracy while detector locations in the apex region and above are advantageous for clinical applications.
Collapse
Affiliation(s)
- Erik B Jørgensen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | | | - Simon Buus
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Lise Bentzen
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | | | - Susanne Rylander
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Kari Tanderup
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Jacob G Johansen
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| |
Collapse
|
13
|
Debnath SBC, Ferre M, Tonneau D, Fauquet C, Tallet A, Goncalves A, Darreon J. High resolution small-scale inorganic scintillator detector: HDR brachytherapy application. Med Phys 2021; 48:1485-1496. [PMID: 33476399 DOI: 10.1002/mp.14727] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Revised: 12/21/2020] [Accepted: 01/07/2021] [Indexed: 12/15/2022] Open
Abstract
PURPOSE Brachytherapy (BT) deals with high gradient internal dose irradiation made up of a complex system where the source is placed nearby the tumor to destroy cancerous cells. A primary concern of clinical safety in BT is quality assurance to ensure the best matches between the delivered and prescribed doses targeting small volume tumors and sparing surrounding healthy tissues. Hence, the purpose of this study is to evaluate the performance of a point size inorganic scintillator detector (ISD) in terms of high dose rate brachytherapy (HDR-BT) treatment. METHODS A prototype of the dose verification system has been developed based on scintillating dosimetry to measure a high dose rate while using an 192 Ir BT source. The associated dose rate is measured in photons/s employing a highly sensitive photon counter (design data: 20 photons/s). Dose measurement was performed as a function of source-to-detector distance according to TG43U1 recommendations. Overall measurements were carried out inside water phantoms keeping the ISD along the BT needle; a minimum of 0.1 cm distance was maintained between each measurement point. The planned dwell times were measured accurately from the difference of two adjacent times of transit. The ISD system performances were also evaluated in terms of dose linearity, energy dependency, scintillation stability, signal-to-noise ratio (SNR), and signal-to-background ratio (SBR). Finally, a comparison was presented between the ISD measurements and results obtained from TG43 reference dataset. RESULTS The detection efficiency of the ISD was verified by measuring the planned dwell times at different dwell positions. Measurements demonstrated that the ISD has a perfectly linear behavior with dose rate (R2 = 1) and shows high SNR (>35) and SBR (>36) values even at the lowest dose rate investigated at around 10 cm from the source. Standard deviation (1σ) remains within 0.03% of signal magnitude, and less than 0.01% STEM signal was monitored at 0.1 cm source-to-detector distance. Stability of 0.54% is achieved, and afterglow stays less than 1% of the total signal in all the irradiations. Excellent symmetrical behavior of the dose rate regarding source position was observed at different radiation planes. Finally, a comparison with TG-43 reference dataset shows that corrected measurements agreed with simulation data within 1.2% and 1.3%, and valid for the source-to-detector distance greater than 0.25 cm. CONCLUSION The proposed ISD in this study anticipated that the system could be promoted to validate with further clinical investigations. It allows an appropriate dose verification with dwell time estimation during source tracking and suitable dose measurement with a high spatial resolution both nearby (high dose gradient) and far (low dose gradient) from the source position.
Collapse
Affiliation(s)
| | | | - Didier Tonneau
- Aix Marseille Université, CNRS, CINaM UMR 7325, Marseille, 13288, France
| | - Carole Fauquet
- Aix Marseille Université, CNRS, CINaM UMR 7325, Marseille, 13288, France
| | - Agnes Tallet
- Institut Paoli-Calmettes, Marseille, 13009, France
| | - Anthony Goncalves
- Institut Paoli-Calmettes, Marseille, 13009, France.,Aix Marseille Université, CNRS UMR 7258, INSERM UMR 1068, CRCM, Marseille, 13009, France
| | | |
Collapse
|
14
|
Muenkel J, Xu Z, Traughber BJ, Baig T, Xu K, Langmack C, Harris E, Podder TK. Feasibility of improving patient's safety with in vivo dose tracking in intracavitary and interstitial HDR brachytherapy. Brachytherapy 2020; 20:353-360. [PMID: 33187822 DOI: 10.1016/j.brachy.2020.09.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2020] [Revised: 08/29/2020] [Accepted: 09/14/2020] [Indexed: 11/19/2022]
Abstract
PURPOSE The in vivo dosimetric monitoring in HDR brachytherapy is important for improving patient safety. However, there are very limited options available for clinical application. In this study, we present a new in vivo dose measurement system with a plastic scintillating detector (PSD) for GYN HDR brachytherapy. METHODS An FDA approved PSD system, called OARtrac (AngioDynamics, Latham, NY), was used with various applicators for in vivo dose measurements for GYN patients. An institutional workflow was established for the clinical implementation of the dosimetric system. Action levels were proposed based on the measurement and system uncertainty for measurement deviations. From October 2018 to September 2019, a total of 75 measurements (48 fractions) were acquired from 14 patients who underwent HDR brachytherapy using either a multichannel cylinder, Venezia applicator, or Syed-Neblett template. The PSDs were placed in predetermined catheters/channels. A planning CT was acquired for treatment planning in Oncentra (Elekta, Version-4.5.2) TPS. The PSDs were contoured on the CT images, and the PSD D90% values were used as the expected doses for comparison with the measured doses. RESULTS The mean difference from patient measurements was -0.22% ± 5.98%, with 26% being the largest deviation from the expected value (Syed case). Large deviations were observed when detectors were placed in the area where dose rates were less than 1 cGy/s. CONCLUSIONS The establishment of clinical workflow for the in vivo dosimetry for both the intracavitary and interstitial GYN HDR brachytherapy will potentially improve the safety of the patient treatment.
Collapse
Affiliation(s)
- Jessica Muenkel
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH
| | - Zhengzheng Xu
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH.
| | - Bryan J Traughber
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH; School of Medicine, Case Western Reserve University, Cleveland, OH
| | - Tanvir Baig
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH
| | - Keying Xu
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH
| | - Christian Langmack
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH
| | - Eleanor Harris
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH; School of Medicine, Case Western Reserve University, Cleveland, OH
| | - Tarun K Podder
- Department of Radiation Oncology, University Hospitals Cleveland Medical Center, Cleveland, OH; School of Medicine, Case Western Reserve University, Cleveland, OH
| |
Collapse
|
15
|
Fonseca GP, Johansen JG, Smith RL, Beaulieu L, Beddar S, Kertzscher G, Verhaegen F, Tanderup K. In vivo dosimetry in brachytherapy: Requirements and future directions for research, development, and clinical practice. PHYSICS & IMAGING IN RADIATION ONCOLOGY 2020; 16:1-11. [PMID: 33458336 PMCID: PMC7807583 DOI: 10.1016/j.phro.2020.09.002] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/14/2020] [Revised: 08/24/2020] [Accepted: 09/17/2020] [Indexed: 12/19/2022]
Abstract
Brachytherapy can deliver high doses to the target while sparing healthy tissues due to its steep dose gradient leading to excellent clinical outcome. Treatment accuracy depends on several manual steps making brachytherapy susceptible to operational mistakes. Currently, treatment delivery verification is not routinely available and has led, in some cases, to systematic errors going unnoticed for years. The brachytherapy community promoted developments in in vivo dosimetry (IVD) through research groups and small companies. Although very few of the systems have been used clinically, it was demonstrated that the likelihood of detecting deviations from the treatment plan increases significantly with time-resolved methods. Time–resolved methods could interrupt a treatment avoiding gross errors which is not possible with time-integrated dosimetry. In addition, lower experimental uncertainties can be achieved by using source-tracking instead of direct dose measurements. However, the detector position in relation to the patient anatomy remains a main source of uncertainty. The next steps towards clinical implementation will require clinical trials and systematic reporting of errors and near-misses. It is of utmost importance for each IVD system that its sensitivity to different types of errors is well understood, so that end-users can select the most suitable method for their needs. This report aims to formulate requirements for the stakeholders (clinics, vendors, and researchers) to facilitate increased clinical use of IVD in brachytherapy. The report focuses on high dose-rate IVD in brachytherapy providing an overview and outlining the need for further development and research.
Collapse
Affiliation(s)
- Gabriel P Fonseca
- Department of Radiation Oncology (Maastro), GROW School for Oncology, Maastricht University Medical Centre+, Maastricht, Doctor Tanslaan 12, 6229 ET Maastricht, the Netherlands
| | - Jacob G Johansen
- Department of Oncology, Aarhus University Hospital, Palle Juul-Jensens Boulevard 99, DK-8200 Aarhus, Denmark
| | - Ryan L Smith
- Alfred Health Radiation Oncology, Alfred Health, 55 Commercial Rd, Melbourne, VIC 3004, Australia
| | - Luc Beaulieu
- Department of Physics, Engineering Physics & Optics and Cancer Research Center, Université Laval, Quebec City, QC, Canada.,Department of Radiation Oncology, Research Center of CHU de Québec, Université Laval, Quebec City, QC, Canada
| | - Sam Beddar
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Blvd., Unit 1420, Houston, TX 77030, United States
| | - Gustavo Kertzscher
- Department of Oncology, Aarhus University Hospital, Palle Juul-Jensens Boulevard 99, DK-8200 Aarhus, Denmark
| | - Frank Verhaegen
- Department of Radiation Oncology (Maastro), GROW School for Oncology, Maastricht University Medical Centre+, Maastricht, Doctor Tanslaan 12, 6229 ET Maastricht, the Netherlands
| | - Kari Tanderup
- Department of Oncology, Aarhus University Hospital, Palle Juul-Jensens Boulevard 99, DK-8200 Aarhus, Denmark
| |
Collapse
|
16
|
Abstract
Many brachytherapy (BT) errors could be detected with real-time in vivo dosimetry technology. Inorganic scintillation detectors (ISDs) have demonstrated promising capabilities for BT, because some ISD materials can generate scintillation signals large enough that (a) the background signal emitted in the fiber-optic cable (stem signal) is insignificant, and (b) small detector volumes can be used to avoid volume averaging effects in steep dose gradients near BT sources. We investigated the characteristics of five ISD materials to identify one that is appropriate for BT. ISDs consisting of a 0.26 to 1.0 mm3 volume of ruby (Al2O3:Cr), a mixture of Y2O3:Eu and YVO4:Eu, ZnSe:O, or CsI:Tl coupled to a fiber-optic cable were irradiated in a water-equivalent phantom using a high-dose-rate 192Ir BT source. Detectors based on plastic scintillators BCF-12 and BCF-60 (0.8 mm3 volume) were used as a reference. Measurements demonstrated that the ruby, Y2O3:Eu+YVO4:Eu, ZnSe:O, and CsI:Tl ISDs emitted scintillation signals that were up to 19, 19, 250, and 880 times greater, respectively, than that of the BCF-12 detector. While the total signals of the plastic scintillation detectors were dominated by the stem signal for source positions 0.5 cm from the fiber-optic cable and >3.5 cm from the scintillator volume, the stem signal for the ruby and Y2O3:Eu+YVO4:Eu ISDs were <1% of the total signal for source positions <3.4 and <4.4 cm from the scintillator, respectively, and <0.7% and <0.5% for the ZnSe:O and CsI:Tl ISDs, respectively, for positions ⩽8.0 cm. In contrast to the other ISDs, the Y2O3:Eu+YVO4:Eu ISD exhibited unstable scintillation and significant afterglow. All ISDs exhibited significant energy dependence, i.e. their dose response to distance-dependent 192Ir energy spectra differed significantly from the absorbed dose in water. Provided that energy dependence is accounted for, ZnSe:O ISDs are promising for use in error detection and patient safety monitoring during BT.
Collapse
Affiliation(s)
- Gustavo Kertzscher
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, United States of America
| | | |
Collapse
|
17
|
Tariq M, Gomez C, Riegel AC. Dosimetric impact of placement errors in optically stimulated luminescent in vivo dosimetry in radiotherapy. PHYSICS & IMAGING IN RADIATION ONCOLOGY 2019; 11:63-68. [PMID: 33458280 PMCID: PMC7807655 DOI: 10.1016/j.phro.2019.08.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Revised: 08/27/2019] [Accepted: 08/28/2019] [Indexed: 11/02/2022]
Abstract
Introduction Studies have suggested that optically stimulated luminescent dosimeters (OSLDs) can be used for in vivo dosimetry of intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT). Clinical uncertainties such as placement error have not been thoroughly investigated. The purpose of this work was to measure OSLD placement error in a clinical sample and analyze its dosimetric impact. Methods and materials The analysis consisted of three parts: first, quantification of placement error in a clinical sample of 128 patients yielding 293 cone-beam CT (CBCT) with visible OSLDs registered to the treatment plan; Second, correlation of placement error and clinical OSLD measurements; third, simulation of dosimeter placement in the treatment plan and correlation of recalculated dose with placement error. Results In the first analysis, average placement error was 9.7 ± 9.5 mm. In the second analysis, placement error and measured-to-planned dose agreement yielded no correlation (R2 = 0.02) for a subsample of 77 CBCTs of 55 head-and-neck patients. Average placement error was 7.0 ± 6.0 mm. Several factors, including image-guided shifts, introduced uncharacterized uncertainty to the measured-to-planned dose agreement. The third analysis isolated placement error from these other effects. Average dosimetric error was -2.4 ± 19.3%. Simulated dosimetric impact was weakly correlated with placement error (R2 = 0.39). Removing outliers reduced the average dosimetric error to -2.1 ± 10.9%, marginally improving the correlation (R2 = 0.44). Conclusion Placement error can substantially impact measured-to-planned dose agreement of OSLDs in high gradient regions, demonstrating the criticality of accurate dosimeter placement for IMRT and VMAT treatments.
Collapse
Affiliation(s)
- Mahin Tariq
- St. Peter's Health Partners, Albany, NY 12208, United States.,Hofstra University, Hempstead, NY 11549, United States.,Northwell Health, New Hyde Park, NY 11402, United States
| | - Cindy Gomez
- Hofstra University, Hempstead, NY 11549, United States
| | - Adam C Riegel
- Hofstra University, Hempstead, NY 11549, United States.,Northwell Health, New Hyde Park, NY 11402, United States.,Donald and Barbara Zucker School of Medicine at Hofstra/Northwell, Hempstead, NY 11549, United States
| |
Collapse
|
18
|
Johansen J, Kertzscher G, Jørgensen E, Rylander S, Bentzen L, Hokland S, Søndergaard C, With A, Buus S, Tanderup K. Dwell time verification in brachytherapy based on time resolved in vivo dosimetry. Phys Med 2019; 60:156-161. [DOI: 10.1016/j.ejmp.2019.03.031] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 02/22/2019] [Accepted: 03/29/2019] [Indexed: 10/27/2022] Open
|
19
|
Johansen JG, Rylander S, Buus S, Bentzen L, Hokland SB, Søndergaard CS, With AKM, Kertzscher G, Tanderup K. Time-resolved in vivo dosimetry for source tracking in brachytherapy. Brachytherapy 2018; 17:122-132. [DOI: 10.1016/j.brachy.2017.08.009] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2017] [Revised: 08/11/2017] [Accepted: 08/15/2017] [Indexed: 11/25/2022]
|
20
|
Tanderup K, Kirisits C, Damato AL. Treatment delivery verification in brachytherapy: Prospects of technology innovation. Brachytherapy 2018; 17:1-6. [PMID: 29406123 DOI: 10.1016/j.brachy.2017.12.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Kari Tanderup
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark.
| | - Christian Kirisits
- Department of Radiotherapy, Medical University of Vienna, Comprehensive Cancer Center, Vienna, Austria
| | - Antonio L Damato
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY
| |
Collapse
|
21
|
Smith RL, Hanlon M, Panettieri V, Millar JL, Matheson B, Haworth A, Franich RD. An integrated system for clinical treatment verification of HDR prostate brachytherapy combining source tracking with pretreatment imaging. Brachytherapy 2018; 17:111-121. [DOI: 10.1016/j.brachy.2017.08.004] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2017] [Revised: 07/29/2017] [Accepted: 08/09/2017] [Indexed: 10/18/2022]
|
22
|
Persson M, Nilsson J, Carlsson Tedgren Å. Experience of using MOSFET detectors for dose verification measurements in an end-to-end 192Ir brachytherapy quality assurance system. Brachytherapy 2017; 17:227-233. [PMID: 29110967 DOI: 10.1016/j.brachy.2017.09.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2016] [Revised: 09/22/2017] [Accepted: 09/22/2017] [Indexed: 11/26/2022]
Abstract
PURPOSE Establishment of an end-to-end system for the brachytherapy (BT) dosimetric chain could be valuable in clinical quality assurance. Here, the development of such a system using MOSFET (metal oxide semiconductor field effect transistor) detectors and experience gained during 2 years of use are reported with focus on the performance of the MOSFET detectors. METHODS AND MATERIALS A bolus phantom was constructed with two implants, mimicking prostate and head & neck treatments, using steel needles and plastic catheters to guide the 192Ir source and house the MOSFET detectors. The phantom was taken through the BT treatment chain from image acquisition to dose evaluation. During the 2-year evaluation-period, delivered doses were verified a total of 56 times using MOSFET detectors which had been calibrated in an external 60Co beam. An initial experimental investigation on beam quality differences between 192Ir and 60Co is reported. RESULTS The standard deviation in repeated MOSFET measurements was below 3% in the six measurement points with dose levels above 2 Gy. MOSFET measurements overestimated treatment planning system doses by 2-7%. Distance-dependent experimental beam quality correction factors derived in a phantom of similar size as that used for end-to-end tests applied on a time-resolved measurement improved the agreement. CONCLUSIONS MOSFET detectors provide values stable over time and function well for use as detectors for end-to-end quality assurance purposes in 192Ir BT. Beam quality correction factors should address not only distance from source but also phantom dimensions.
Collapse
Affiliation(s)
- Maria Persson
- Department of Medical Radiation Physics and Nuclear Medicine, Section for Radiotherapy Physics and Engineering, Karolinska University Hospital, Stockholm, Sweden.
| | - Josef Nilsson
- Department of Medical Radiation Physics and Nuclear Medicine, Section for Radiotherapy Physics and Engineering, Karolinska University Hospital, Stockholm, Sweden
| | - Åsa Carlsson Tedgren
- Department of Medical Radiation Physics and Nuclear Medicine, Section for Radiotherapy Physics and Engineering, Karolinska University Hospital, Stockholm, Sweden; Department of Medical and Health Sciences and Center for Medical Image Science and Visualization, Radiation Physics, Linköping University, Linköping, Sweden
| |
Collapse
|
23
|
Smith RL, Haworth A, Panettieri V, Millar JL, Franich RD. 3D catheter reconstruction in HDR prostate brachytherapy for pre-treatment verification using a flat panel detector. Phys Med 2017. [PMID: 28625472 DOI: 10.1016/j.ejmp.2017.06.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
PURPOSE High dose rate prostate brachytherapy is a widely-practiced treatment, delivering large conformal doses in relatively few treatment fractions. Inter- and intra-fraction catheter displacements have been reported. Unrecognized displacement can have a significant impact on dosimetry. Knowledge of the implant geometry at the time of treatment is important for ensuring safe and effective treatment. In this work we demonstrate a method to reconstruct the catheter positions pre-treatment, using a 'shift' imaging technique, and perform registration with the treatment plan for verification relative to the prostate. METHODS Two oblique 'shift' images were acquired of a phantom containing brachytherapy catheters, representing the patient immediately pre-treatment. Using a back projection approach, the catheter paths were reconstructed in 3D and registered with the planned catheter paths. The robustness of the reconstruction and registration process was investigated as a function of phantom rotation. Catheter displacement detection was performed and compared to known applied displacements. RESULTS Reconstruction of the implant geometry in 3D immediately prior to treatment was achieved. A mean reconstruction uncertainty of 0.8mm was determined for all catheters with a mean registration uncertainty of 0.5mm. A catheter displacement detection threshold of 2.2mm was demonstrated. Catheter displacements were all detected to within 0.5mm of the applied displacements. CONCLUSION This technique is robust and sensitive to assess catheter displacements throughout the implant volume. This approach provides a method to detect, in 3D, changes in catheter positions relative to the prostate. The method has sufficient sensitivity to enable clinically significant decisions immediately prior to treatment delivery.
Collapse
Affiliation(s)
- Ryan L Smith
- Alfred Health Radiation Oncology, The Alfred Hospital, Melbourne 3004, VIC, Australia; School of Science, RMIT University, Melbourne 3000, VIC, Australia.
| | - Annette Haworth
- School of Science, RMIT University, Melbourne 3000, VIC, Australia; School of Physics, University of Sydney, Camperdown 2006, NSW, Australia
| | - Vanessa Panettieri
- Alfred Health Radiation Oncology, The Alfred Hospital, Melbourne 3004, VIC, Australia
| | - Jeremy L Millar
- Alfred Health Radiation Oncology, The Alfred Hospital, Melbourne 3004, VIC, Australia; School of Science, RMIT University, Melbourne 3000, VIC, Australia; Central Clinical School, Monash University, Melbourne 3004, Victoria, Australia
| | - Rick D Franich
- School of Science, RMIT University, Melbourne 3000, VIC, Australia
| |
Collapse
|
24
|
Kertzscher G, Beddar S. Ruby-based inorganic scintillation detectors for 192Ir brachytherapy. Phys Med Biol 2016; 61:7744-7764. [PMID: 27740947 DOI: 10.1088/0031-9155/61/21/7744] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We tested the potential of ruby inorganic scintillation detectors (ISDs) for use in brachytherapy and investigated various unwanted luminescence properties that may compromise their accuracy. The ISDs were composed of a ruby crystal coupled to a poly(methyl methacrylate) fiber-optic cable and a charge-coupled device camera. The ISD also included a long-pass filter that was sandwiched between the ruby crystal and the fiber-optic cable. The long-pass filter prevented the Cerenkov and fluorescence background light (stem signal) induced in the fiber-optic cable from striking the ruby crystal, which generates unwanted photoluminescence rather than the desired radioluminescence. The relative contributions of the radioluminescence signal and the stem signal were quantified by exposing the ruby detectors to a high-dose-rate brachytherapy source. The photoluminescence signal was quantified by irradiating the fiber-optic cable with the detector volume shielded. Other experiments addressed time-dependent luminescence properties and compared the ISDs to commonly used organic scintillator detectors (BCF-12, BCF-60). When the brachytherapy source dwelled 0.5 cm away from the fiber-optic cable, the unwanted photoluminescence was reduced from >5% to <1% of the total signal as long as the ISD incorporated the long-pass filter. The stem signal was suppressed with a band-pass filter and was <3% as long as the source distance from the scintillator was <7 cm. Some ruby crystals exhibited time-dependent luminescence properties that altered the ruby signal by >5% within 10 s from the onset of irradiation and after the source had retracted. The ruby-based ISDs generated signals of up to 20 times that of BCF-12-based detectors. The study presents solutions to unwanted luminescence properties of ruby-based ISDs for high-dose-rate brachytherapy. An optic filter should be sandwiched between the ruby crystal and the fiber-optic cable to suppress the photoluminescence. Furthermore, we recommend avoiding ruby crystals that exhibit significant time-dependent luminescence.
Collapse
Affiliation(s)
- Gustavo Kertzscher
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | | |
Collapse
|
25
|
|
26
|
Smith RL, Haworth A, Panettieri V, Millar JL, Franich RD. A method for verification of treatment delivery in HDR prostate brachytherapy using a flat panel detector for both imaging and source tracking. Med Phys 2016; 43:2435. [DOI: 10.1118/1.4946820] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
27
|
Kertzscher G, Rosenfeld A, Beddar S, Tanderup K, Cygler JE. In vivo dosimetry: trends and prospects for brachytherapy. Br J Radiol 2014; 87:20140206. [PMID: 25007037 DOI: 10.1259/bjr.20140206] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
The error types during brachytherapy (BT) treatments and their occurrence rates are not well known. The limited knowledge is partly attributed to the lack of independent verification systems of the treatment progression in the clinical workflow routine. Within the field of in vivo dosimetry (IVD), it is established that real-time IVD can provide efficient error detection and treatment verification. However, it is also recognized that widespread implementations are hampered by the lack of available high-accuracy IVD systems that are straightforward for the clinical staff to use. This article highlights the capabilities of the state-of-the-art IVD technology in the context of error detection and quality assurance (QA) and discusses related prospects of the latest developments within the field. The article emphasizes the main challenges responsible for the limited practice of IVD and provides descriptions on how they can be overcome. Finally, the article suggests a framework for collaborations between BT clinics that implemented IVD on a routine basis and postulates that such collaborations could improve BT QA measures and the knowledge about BT error types and their occurrence rates.
Collapse
Affiliation(s)
- G Kertzscher
- 1 Centre for Nuclear Technologies, Technical University of Denmark, Roskilde, Denmark
| | | | | | | | | |
Collapse
|