1
|
Lavallee MC, Cantin A, Aubin S, Lefebvre M, Marrier AS, Bouchard I, Fiset C, Villeneuve-Gauthier A, Foster W, Martin AG, Carignan D, Beaulieu L, Vigneault E. US-guided EM tracked system for HDR brachytherapy: A first in-men randomized study for whole prostate treatment. Brachytherapy 2024; 23:64-72. [PMID: 37806788 DOI: 10.1016/j.brachy.2023.08.003] [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/03/2023] [Revised: 07/10/2023] [Accepted: 08/09/2023] [Indexed: 10/10/2023]
Abstract
PURPOSE An electromagnetic tracking device (EMT) has been integrated in an HDR 3D ultrasound guidance system for prostate HDR. The aim of this study was to compare the efficiency of HDR workflows with and without EM tracking. METHODS AND MATERIALS A total of 58 patients with a 15 Gy HDR prostate boost were randomized in two arms and two operation room (OR) procedures using: (1) the EMT investigational device, and (2) the Oncentra prostate system (OCP). OR times were compared for both techniques. RESULTS The overall procedure median time was about 20% shorter for EMT (63 min) compared to OCP (79 min). The US acquisition and contouring was longer for OCP compared to EMT (23 min vs. 16 min). The catheter reconstruction's median times were 23 min and 13 min for OCP and EMT respectively. For the automatic reconstruction with EMT, 62% of cases required no or few manual corrections. Using the EM technology in an OR environment was challenging. In some cases, interferences or the stiffness of the stylet introduced errors in the reconstruction of catheters. The last step was the dosimetry with median times of 11 min (OCP) and 15.5 min (EMT). Finally, it was observed that there was no learning curve associated with the introduction of this new technology. CONCLUSIONS The EMT device offers an efficient solution for automatic catheter reconstruction for HDR prostate while reducing the possibility of mis-reconstructed catheters caused by issues of visualization in the US images. Because of that, the overall OR times was shorter when using the EMT system.
Collapse
Affiliation(s)
- Marie-Claude Lavallee
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada; Département de physique, de génie physique et d'optique, Université Laval, Québec, Canada
| | - Audrey Cantin
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada
| | - Sylviane Aubin
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada
| | - Martine Lefebvre
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada
| | | | | | - Cedric Fiset
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada
| | | | - William Foster
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada
| | - Andre-Guy Martin
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada; CRCHU de Québec and Centre de recherche sur le cancer de l'Université Laval, Québec, Canada
| | - Damien Carignan
- CRCHU de Québec and Centre de recherche sur le cancer de l'Université Laval, Québec, Canada
| | - Luc Beaulieu
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada; CRCHU de Québec and Centre de recherche sur le cancer de l'Université Laval, Québec, Canada; Département de physique, de génie physique et d'optique, Université Laval, Québec, Canada
| | - Eric Vigneault
- CHU de Québe - Université Laval, Radiation Oncology, Québec, Canada; CRCHU de Québec and Centre de recherche sur le cancer de l'Université Laval, Québec, Canada.
| |
Collapse
|
2
|
Decabooter E, Hilgers GC, De Rouck J, Salvo K, Van Wingerden J, Bosmans H, van der Heyden B, Qamhiyeh S, Papalazarou C, Kaatee R, Pittomvils G, Bogaert E. Survey on fan-beam computed tomography for radiotherapy: Imaging for dose calculation and delineation. Phys Imaging Radiat Oncol 2024; 29:100522. [PMID: 38152701 PMCID: PMC10750173 DOI: 10.1016/j.phro.2023.100522] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Revised: 11/26/2023] [Accepted: 11/28/2023] [Indexed: 12/29/2023] Open
Abstract
Background and purpose To obtain an understanding of current practice, professional needs and future directions in the field of fan-beam CT in RT, a survey was conducted. This work presents the collected information regarding the use of CT imaging for dose calculation and structure delineation. Materials and methods An online institutional survey was distributed to medical physics experts employed at Belgian and Dutch radiotherapy institutions to assess the status, challenges, and future directions of QA practices for fan-beam CT. A maximum of 143 questions covered topics such as CT scanner availability, CT scanner specifications, QA protocols, treatment simulation workflow, and radiotherapy dose calculation. Answer forms were collected between 1-Sep-2022 and 10-Oct-2022. Results A 66 % response rate was achieved, yielding data on a total of 58 CT scanners. For MV photon therapy, all single-energy CT scans are reconstructed in Hounsfield Units for delineation or dose calculation, and a direct- or stoichiometric method was used to convert CT numbers for dose calculation. Limited use of dual-energy CT is reported for photon (N = 3) and proton dose calculations (N = 1). For brachytherapy, most institutions adopt water-based dose calculation, while approximately 26 % of the institutions take tissue heterogeneity into account. Commissioning and regular QA include eleven tasks, which are performed by two or more professions (29/31) with varying frequencies. Conclusions Dual usage of a planning CT limits protocol optimization for both tissue characterization and delineation. DECT has been implemented only gradually. A variation of QA testing frequencies and tests are reported.
Collapse
Affiliation(s)
- Esther Decabooter
- Department of Radiation Oncology (Maastro Clinic), GROW School for Oncology, Maastricht University Medical Centre+, Maastricht, The Netherlands
| | | | - Joke De Rouck
- Department of Radiotherapy, AZ Sint Lucas, Ghent, Belgium
| | - Koen Salvo
- Department of Radiotherapy, AZ Sint-Maarten, Mechelen, Belgium
| | - Jacobus Van Wingerden
- Department of Medical Physics, Haaglanden Medical Centre, Leidschendam, The Netherlands
| | - Hilde Bosmans
- Department of Medical Radiation Physics, University Hospital Leuven, Belgium
| | - Brent van der Heyden
- IBiTech-MEDISIP, Department of Electronics and Information Systems, Ghent University, Ghent, Belgium
- Department of Oncology, Laboratory of Experimental Radiotherapy, KU Leuven, Leuven, Belgium
| | - Sima Qamhiyeh
- University Hospitals Leuven, Department of Radiation Oncology, Leuven, Belgium
| | - Chrysi Papalazarou
- Department of Radiotherapy, Leiden University Medical Center, Leiden, The Netherlands
| | - Robert Kaatee
- Radiotherapy Institute Friesland, Leeuwarden, The Netherlands
| | - Geert Pittomvils
- Department of Radiation-Oncology, Ghent University Hospital, Ghent, Belgium
| | - Evelien Bogaert
- Department of Radiation-Oncology, Ghent University Hospital, Ghent, Belgium
| |
Collapse
|
3
|
Automating implant reconstruction in interstitial brachytherapy of the breast: A hybrid approach combining electromagnetic tracking and image segmentation. Radiother Oncol 2022; 176:172-178. [PMID: 36181920 DOI: 10.1016/j.radonc.2022.09.019] [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: 06/26/2022] [Revised: 08/24/2022] [Accepted: 09/23/2022] [Indexed: 12/14/2022]
Abstract
BACKGROUND AND PURPOSE To develop a method for automatic reconstruction of catheter implants in interstitial brachytherapy (iBT) of the breast by means of electromagnetic tracking (EMT) with the goal of making treatment planning as time-effective and accurate as possible. MATERIALS AND METHODS The implant geometry of 64 patients was recorded using an afterloader prototype with EMT functionality immediately after the planning CT. EMT data were transferred to the CT image space by rigidly registering the catheter fixation buttons as landmarks. To further improve reconstruction accuracy, the EMT reconstruction points were used as starting points to define small regions of interest (ROI) in the CT image. Within these ROIs, the catheter track was segmented in the CT using image processing operations such as thresholding and blob detection, thus refining the reconstruction. The perpendicular distance between the refined EMT implant reconstruction points and the manually reconstructed catheters by an experienced treatment planner was calculated as a measure of their geometric agreement. RESULTS Geometrically, the refined EMT based implant reconstruction shows excellent agreement with the manual reconstruction. The median distance across all patients is 0.25 mm and the 95th percentile is 1 mm. Refinement takes approximately 0.5 s per reconstruction point and typically does not exceed 3 min per implant at no user interaction. CONCLUSION The refined EMT based implant reconstruction proved to be extremely accurate and fast compared to manual reconstruction. The presented procedure can in principle be easily transferred to clinical routine and therefore has enormous potential to provide significant time savings in iBT treatment planning whilst improving reconstruction accuracy.
Collapse
|
4
|
Lavallée MC, Cantin A, Monéger F, Lefebvre M, Foster W, Vigneault E, Beaulieu L. Commissioning of an intra-operative US guided prostate HDR system integrating an EM tracking technology. Brachytherapy 2021; 20:1296-1304. [PMID: 34272132 DOI: 10.1016/j.brachy.2021.05.163] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Revised: 05/05/2021] [Accepted: 05/31/2021] [Indexed: 11/17/2022]
Abstract
PURPOSE Ultrasound-based planning for high-dose-rate prostate brachytherapy is commonly used in the clinic, mainly because it offers fast real-time image-guided capability at a relatively low cost. The main difficulty with US planning is the catheter reconstruction due to artefacts (from multiple catheters) and echogenicity. Electromagnetic tracking (EMT) system offers a fast and accurate solution for automatic reconstruction of catheters using the EMT technology. In this study, the commissioning and performance evaluation of the new real-time prostate high-dose-rate brachytherapy investigational system from Philips Disease Management Solutions integrating EMT was performed before its clinical integration. METHOD AND MATERIALS The Philips' clinical investigational system includes a treatment planning software (TPS) that was commissioned based on AAPM TG53 and TG56 recommendations for the use of TPS in brachytherapy. First, the CIRS - model 045A - QA phantom was used to evaluate the ultrasound (US) image quality and 3D image handling. Distances, volumes, and dimensions of the structures inside the phantom were measured and compared to the actual values. The calibration reproducibility and accuracy of the electromagnetic (EM) sensor used to track the US probe (rotation and translation) were performed using a specifically designed QA tool mounted on the probe and immersed in a salted water tank. This was performed for 3 different B&K 8848 US probes to evaluate the sensitivity of EM calibration to the probe geometric properties (manufacturing process). The new TPS performance was compared to that in OncentraBrachy (OcB) V4.5.5 (Elekta) using 30 clinical cases as part of a retrospective study. Following the system commissioning, clinical workflows were explored; tests were performed with the brachytherapy team on phantoms and finally implemented in the clinic. RESULTS US image quality evaluation showed a mean difference with actual dimensions (lengths, widths and distances) of 0.4 mm (±0.3 mm) and mean difference in volume sizes of 0.2 cc (±0.2 cc). Then, the calibration of the US-to-EM coordinate system was performed for 3 different probes. For each probe, 3 measurements were acquired for every position of the calibration tool and measurements were repeated 3 times for a total of 27 measurements per probe per plane. The error was slightly higher in transverse mode compared to sagittal mode with mean values of 0.6 ± 0.2 mm and 0.3 ± 0.1 mm respectively. 30 clinical cases were used to compare the new TPS performance to OcB (IPSA). Optimized plans obtained with both systems were all clinically acceptable, but the plans from the Philips system have slightly higher V150% values, V200% values and dose to organs at risk. In the case of organs at risk, plans could have been manually modified to reduce the dose. Philips' system had a larger number of active dwell positions and longer treatment times. CONCLUSIONS The first clinical version of Philips' system was proven to be stable, accurate and precise. The fully integrated EM tracking technology opens the way for automated catheter reconstruction and on-the-fly dynamical replanning.
Collapse
Affiliation(s)
| | | | - Fanny Monéger
- CRCHU de Quebec and Centre de recherche en Cancérologie Université Laval, Québec, Canada; Polytech Clermont-Ferrand, Département de génie physique, Aubière, France
| | | | | | - Eric Vigneault
- CHU de Québec, Radiation Oncology, Québec, Canada; CRCHU de Quebec and Centre de recherche en Cancérologie Université Laval, Québec, Canada
| | - Luc Beaulieu
- CHU de Québec, Radiation Oncology, Québec, Canada; CRCHU de Quebec and Centre de recherche en Cancérologie Université Laval, Québec, Canada; Département de physique, de génie physique et d'optique, Université Laval, Québec, Canada
| |
Collapse
|
5
|
Deufel CL, Tian S, Yan BB, Vaishnav BD, Haddock MG, Petersen IA. Automated applicator digitization for high-dose-rate cervix brachytherapy using image thresholding and density-based clustering. Brachytherapy 2019; 19:111-118. [PMID: 31594729 DOI: 10.1016/j.brachy.2019.09.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Revised: 08/13/2019] [Accepted: 09/09/2019] [Indexed: 10/25/2022]
Abstract
PURPOSE The purpose of the study was to develop and evaluate an automated digitization algorithm for high-dose-rate cervix brachytherapy, with the goal of reducing the duration of treatment planning, staff resources, variability, and potential for human error. METHODS An automated digitization algorithm was developed and retrospectively evaluated using treatment planning data from 10 patients with cervix cancer who were treated with a titanium tandem and ovoids applicator set. Applicators were segmented, without human interaction, by thresholding CT images to isolate high-density voxels and assigning the voxels to applicator and nonapplicator structures using HDBSCAN, a density-based linkage clustering algorithm. The applicator contours were determined from the centroid of the clustered voxels on each image slice and written to a treatment plan file. Automated contours were evaluated against manual digitization using distance and dosimetric metrics. RESULTS A close agreement between automatic and manual digitization was observed. The mean magnitude of contour disagreement for 10 patients equaled 0.3 mm. Hausdorff distances were ≤1.0 mm. The applicator tip coordinates had submillimeter agreement. The median and mean dose volume histogram parameter differences were less than or equal to 1% for high-risk clinical target volume D90, high-risk clinical target volume D95, bladder D2cc, rectum D2cc, large bowel D2cc, and small bowel D2cc. The average execution time for the automated algorithm was less than 30 s. CONCLUSION The digitization of titanium tandem and ovoids applicators for high-dose-rate brachytherapy treatment planning can be automated using straightforward thresholding and clustering algorithms. The adoption of automated digitization is expected to improve the consistency of treatment plans and reduce the duration of treatment planning.
Collapse
Affiliation(s)
| | - Shulan Tian
- Division of Biomedical Statistics and Informatics, Department of Health Sciences Research, Mayo Clinic, Rochester, MN
| | - Benjamin B Yan
- Department of Radiation Oncology, Mayo Clinic, Rochester, MN
| | | | | | - Ivy A Petersen
- Department of Radiation Oncology, Mayo Clinic, Rochester, MN
| |
Collapse
|
6
|
Kallis K, Kaltsas T, Kreppner S, Lotter M, Strnad V, Fietkau R, Bert C. Impact of inter- and intra-observer variabilities of catheter reconstruction on multi-catheter interstitial brachytherapy of breast cancer patients. Radiother Oncol 2019; 135:25-32. [PMID: 31015167 DOI: 10.1016/j.radonc.2019.02.017] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 02/20/2019] [Accepted: 02/21/2019] [Indexed: 10/27/2022]
Abstract
PURPOSE The aim of this study was to evaluate inter- and intra-observer variabilities of catheter reconstruction and its dosimetric impact for multi-catheter interstitial breast cancer patients. METHODS AND MATERIALS In order to evaluate inter-observer variabilities (IOV) three medical physicists reconstructed the catheter traces of 13 patients. These manual reconstructions were further compared to the automatic reconstruction algorithm integrated into the planning system and one on purpose imprecise manual reconstruction. For intra-observer variabilities (IAV) repeated reconstructions of two physicists were compared for 13 patients. In total 426 catheters were considered. Keeping dwell times, dwell positions, the optimization and the normalization relative points constant the geometrical deviations between the corresponding dwell positions of the reference data set and the investigated reconstructions were evaluated. Also, the effect on the quality indices, such as coverage index (CI), dose non-uniformity ratio (DNR) or conformal index (COIN), and the exposure of the organs at risk were analyzed. RESULTS Over all patients and all different catheter reconstructions considered for IOV a mean deviation between the corresponding dwell positions of 0.60 ± 0.35 mm was detected. The first observer had a mean deviation of 0.54 ± 0.32 mm, whereas the second observer yielded a mean deviation of 0.58 ± 0.37 mm. The length of the catheter traces varied in the mean by 0.51 ± 0.45 mm. The mean relative deviation of the CI, DNR, COIN, mean heart dose and mean lung dose varied by 0.27 ± 0.31%, 0.0027 ± 0.0025, 0.0036 ± 0.0033, 0.024 ± 0.019%, 0.05 ± 0.11%, respectively. The skin dose (D0.2ccm) changed in the maximum 8.52%. On average IAV reached a deviation between the corresponding dwell positions of 0.49 ± 0.30 mm. IOVs and IAVs proved to be significantly different (Wilcoxon's test p < 0.01). CONCLUSIONS The study proved that a repeated reconstruction of the catheter traces does not lead to large geometrical deviations or to a significant change in the dose exposure. But the lack of ground truth makes the estimation of the quality of the reconstruction challenging. A precise reconstruction mapping the reality is a necessity for the planned dose delivery. With all considered reconstruction techniques reliable quality indices for the target and the organs at risk could be obtained.
Collapse
Affiliation(s)
- Karoline Kallis
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Theodoros Kaltsas
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Stephan Kreppner
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Michael Lotter
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Vratislav Strnad
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Rainer Fietkau
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Christoph Bert
- Department of Radiation Oncology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.
| |
Collapse
|
7
|
Wang W, Viswanathan AN, Damato AL, Chen Y, Tse Z, Pan L, Tokuda J, Seethamraju RT, Dumoulin CL, Schmidt EJ, Cormack RA. Evaluation of an active magnetic resonance tracking system for interstitial brachytherapy. Med Phys 2016; 42:7114-21. [PMID: 26632065 DOI: 10.1118/1.4935535] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE In gynecologic cancers, magnetic resonance (MR) imaging is the modality of choice for visualizing tumors and their surroundings because of superior soft-tissue contrast. Real-time MR guidance of catheter placement in interstitial brachytherapy facilitates target coverage, and would be further improved by providing intraprocedural estimates of dosimetric coverage. A major obstacle to intraprocedural dosimetry is the time needed for catheter trajectory reconstruction. Herein the authors evaluate an active MR tracking (MRTR) system which provides rapid catheter tip localization and trajectory reconstruction. The authors assess the reliability and spatial accuracy of the MRTR system in comparison to standard catheter digitization using magnetic resonance imaging (MRI) and CT. METHODS The MRTR system includes a stylet with microcoils mounted on its shaft, which can be inserted into brachytherapy catheters and tracked by a dedicated MRTR sequence. Catheter tip localization errors of the MRTR system and their dependence on catheter locations and orientation inside the MR scanner were quantified with a water phantom. The distances between the tracked tip positions of the MRTR stylet and the predefined ground-truth tip positions were calculated for measurements performed at seven locations and with nine orientations. To evaluate catheter trajectory reconstruction, fifteen brachytherapy catheters were placed into a gel phantom with an embedded catheter fixation framework, with parallel or crossed paths. The MRTR stylet was then inserted sequentially into each catheter. During the removal of the MRTR stylet from within each catheter, a MRTR measurement was performed at 40 Hz to acquire the instantaneous stylet tip position, resulting in a series of three-dimensional (3D) positions along the catheter's trajectory. A 3D polynomial curve was fit to the tracked positions for each catheter, and equally spaced dwell points were then generated along the curve. High-resolution 3D MRI of the phantom was performed followed by catheter digitization based on the catheter's imaging artifacts. The catheter trajectory error was characterized in terms of the mean distance between corresponding dwell points in MRTR-generated catheter trajectory and MRI-based catheter digitization. The MRTR-based catheter trajectory reconstruction process was also performed on three gynecologic cancer patients, and then compared with catheter digitization based on MRI and CT. RESULTS The catheter tip localization error increased as the MRTR stylet moved further off-center and as the stylet's orientation deviated from the main magnetic field direction. Fifteen catheters' trajectories were reconstructed by MRTR. Compared with MRI-based digitization, the mean 3D error of MRTR-generated trajectories was 1.5 ± 0.5 mm with an in-plane error of 0.7 ± 0.2 mm and a tip error of 1.7 ± 0.5 mm. MRTR resolved ambiguity in catheter assignment due to crossed catheter paths, which is a common problem in image-based catheter digitization. In the patient studies, the MRTR-generated catheter trajectory was consistent with digitization based on both MRI and CT. CONCLUSIONS The MRTR system provides accurate catheter tip localization and trajectory reconstruction in the MR environment. Relative to the image-based methods, it improves the speed, safety, and reliability of the catheter trajectory reconstruction in interstitial brachytherapy. MRTR may enable in-procedural dosimetric evaluation of implant target coverage.
Collapse
Affiliation(s)
- Wei Wang
- Department of Radiology, Brigham and Women's Hospital, Boston, Massachusetts 02115 and Department of Radiation Oncology, Dana-Farber Cancer Institute/Brigham and Women's Hospital, Boston, Massachusetts 02115
| | - Akila N Viswanathan
- Department of Radiation Oncology, Dana-Farber Cancer Institute/Brigham and Women's Hospital, Boston, Massachusetts 02115
| | - Antonio L Damato
- Department of Radiation Oncology, Dana-Farber Cancer Institute/Brigham and Women's Hospital, Boston, Massachusetts 02115
| | - Yue Chen
- Department of Engineering, The University of Georgia, Athens, Georgia 30602
| | - Zion Tse
- Department of Engineering, The University of Georgia, Athens, Georgia 30602
| | - Li Pan
- Siemens Healthcare USA, Baltimore, Maryland 21287
| | - Junichi Tokuda
- Department of Radiology, Brigham and Women's Hospital, Boston, Massachusetts 02115
| | | | - Charles L Dumoulin
- Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio 45229
| | - Ehud J Schmidt
- Department of Radiology, Brigham and Women's Hospital, Boston, Massachusetts 02115
| | - Robert A Cormack
- Department of Radiation Oncology, Dana-Farber Cancer Institute/Brigham and Women's Hospital, Boston, Massachusetts 02115
| |
Collapse
|
8
|
Poulin E, Racine E, Binnekamp D, Beaulieu L. Fast, automatic, and accurate catheter reconstruction in HDR brachytherapy using an electromagnetic 3D tracking system. Med Phys 2015; 42:1227-32. [DOI: 10.1118/1.4908011] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
9
|
Pokhrel D, Murphy MJ, Todor DA, Weiss E, Williamson JF. Localizing intracavitary brachytherapy applicators from cone-beam CT x-ray projections via a novel iterative forward projection matching algorithm. Med Phys 2011; 38:1070-80. [PMID: 21452744 DOI: 10.1118/1.3544661] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To present a novel method for reconstructing the 3D pose (position and orientation) of radio-opaque applicators of known but arbitrary shape from a small set of 2D x-ray projections in support of intraoperative brachytherapy planning. METHODS The generalized iterative forward projection matching (gIFPM) algorithm finds the six degree-of-freedom pose of an arbitrary rigid object by minimizing the sum-of-squared-intensity differences (SSQD) between the computed and experimentally acquired autosegmented projection of the objects. Starting with an initial estimate of the object's pose, gIFPM iteratively refines the pose parameters (3D position and three Euler angles) until the SSQD converges. The object, here specialized to a Fletcher-Weeks intracavitary brachytherapy (ICB) applicator, is represented by a fine mesh of discrete points derived from complex combinatorial geometric models of the actual applicators. Three pairs of computed and measured projection images with known imaging geometry are used. Projection images of an intrauterine tandem and colpostats were acquired from an ACUITY cone-beam CT digital simulator. An image postprocessing step was performed to create blurred binary applicators only images. To quantify gIFPM accuracy, the reconstructed 3D pose of the applicator model was forward projected and overlaid with the measured images and empirically calculated the nearest-neighbor applicator positional difference for each image pair. RESULTS In the numerical simulations, the tandem and colpostats positions (x,y,z) and orientations (alpha, beta, gamma) were estimated with accuracies of 0.6 mm and 2 degrees, respectively. For experimentally acquired images of actual applicators, the residual 2D registration error was less than 1.8 mm for each image pair, corresponding to about 1 mm positioning accuracy at isocenter, with a total computation time of less than 1.5 min on a 1 GHz processor. CONCLUSIONS This work describes a novel, accurate, fast, and completely automatic method to localize radio-opaque applicators of arbitrary shape from measured 2D x-ray projections. The results demonstrate approximately 1 mm accuracy while compared against the measured applicator projections. No lateral film is needed. By localizing the applicator internal structure as well as radioactive sources, the effect of intra-applicator and interapplicator attenuation can be included in the resultant dose calculations. Further validation tests using clinically acquired tandem and colpostats images will be performed for the accurate and robust applicator/sources localization in ICB patients.
Collapse
Affiliation(s)
- Damodar Pokhrel
- Department of Radiation Oncology, School of Medicine, Virginia Commonwealth University, Richmond, Virginia 23298, USA
| | | | | | | | | |
Collapse
|
10
|
Kolkman-Deurloo IKK, Deleye XGJ, Jansen PP, Koper PCM. Anatomy based inverse planning in HDR prostate brachytherapy. Radiother Oncol 2005; 73:73-7. [PMID: 15465149 DOI: 10.1016/j.radonc.2004.08.014] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2004] [Revised: 08/02/2004] [Accepted: 08/16/2004] [Indexed: 10/26/2022]
Abstract
The purpose of this study is to evaluate anatomy based inverse planning as implemented in PLATO BPS 14.2 for planning of HDR prostate implants. Six patients were analysed. The dose distributions were optimized using geometric optimization followed by graphical optimization (GO), anatomy based inverse planning or standard inverse optimization (SIO), tuned inverse optimization (TIO) and tuned inverse optimization followed by graphical optimization (GOTIO). The mean target coverage was 93+/-4%, 53+/-11%, 74+/-8%, 90+/-3%, respectively, for GO, SIO, TIO and GOTIO. The conformal index COIN was 0.74+/-0.02, 0.43+/-0.15 and 0.77+/-0.07, respectively, for GO, SIO and GOTIO. Improved dose homogeneity was found when comparing GOTIO with GO.
Collapse
Affiliation(s)
- Inger-Karine K Kolkman-Deurloo
- Department of Radiation Oncology, Erasmus MC-Daniel Den Hoed Cancer Center, Groene Hilledijk 301, 3075 EA Rotterdam, The Netherlands
| | | | | | | |
Collapse
|
11
|
Chang JN, Suh TS, Park SY, Cho KH, Kim S. A dual-energy technique for enhanced localization accuracy in intracavitary brachytherapy. Med Phys 2005; 32:376-9. [PMID: 15789582 DOI: 10.1118/1.1843352] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
The orthogonal imaging method is commonly used for source localization in brachytherapy. In some cases, however, difficulty is encountered in determining the dummy sources because of the presence of either contrast materials or bony structures. We here offer a novel method for source localization utilizing a dual-energy, radiographic technique. In this approach, two sets of orthogonal radiographic images (anterior-posterior and lateral views) are obtained using two different x-ray energies. Image processing (i.e., subtraction between two image sets) is carried out to enhance the source image. In a study performed using a laboratory developed pelvic phantom, it was demonstrated that the dual-energy method could significantly enhance the image quality of the dummy sources, and improve the achievable precision and relative accuracy in localization of source positions. When directly combined with digital imaging modalities, the dual-energy method can be a useful technique to improve the accuracy in brachytherapy source localization from planar radiographs.
Collapse
Affiliation(s)
- Ji-Na Chang
- Department of Biomedical Engineering, College of Medicine, Catholic University, Seoul, Korea
| | | | | | | | | |
Collapse
|
12
|
Kim Y, Hsu ICJ, Lessard E, Pouliot J, Vujic J. Dose uncertainty due to computed tomography (CT) slice thickness in CT-based high dose rate brachytherapy of the prostate cancer. Med Phys 2004; 31:2543-8. [PMID: 15487736 DOI: 10.1118/1.1785454] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
In computed tomography (CT)-based high dose rate (HDR) brachytherapy, the uncertainty in the localization of the longitudinal catheter-tip positions due to the discrete CT slice thickness, results in a delivered dose uncertainty. Catheter coordinates were extracted from five patients treated for prostate cancer, and three simulation scenarios were followed to mimic the longitudinal imprecision of the catheter tips, hence the dwell positions. All catheters were displaced (1) forward, (2) backward, or (3) randomly distributed within the space defined by one CT slice thickness, for thicknesses ranging from 2 to 5 mm. Average and standard deviation values of the relative dose variations are reported for the various catheter displacement scenarios. Also, the dose points were grouped according to their relative position in the prostate, inner, peripheral and outer area of prostate and base, median and apex zones, in order to estimate the spatial sensitivity of the dose errors. For scenarios (1) and (2), the dose uncertainties due to the finite slice thickness increase linearly with the slice spacing, from 3% to 8% for the slice thickness values ranging from 2 to 5 mm, respectively. The more realistic scenario (3) yields average errors ranging from 0.7% to 1.7%. The apex and the base show larger dose errors and variability of dose errors than the median of the prostate. No statistical difference was observed among different transversal sections of the prostate. A CT slice thickness of 3 mm appears to be a good compromise showing an acceptable average dose uncertainty of 1%, without unduly increasing the number of slices.
Collapse
Affiliation(s)
- Yongbok Kim
- Department of Radiation Oncology, University of California, San Francisco, Mount Zion Comprehensive Cancer Center, 1600 Divisadero Street, San Francisco, California 94143-1708, USA.
| | | | | | | | | |
Collapse
|
13
|
Pantelis E, Papagiannis P, Anagnostopoulos G, Baltas D, Karaiskos P, Sandilos P, Sakelliou L. Evaluation of a TG-43 compliant analytical dosimetry model in clinical192Ir HDR brachytherapy treatment planning and assessment of the significance of source position and catheter reconstruction uncertainties. Phys Med Biol 2003; 49:55-67. [PMID: 14971772 DOI: 10.1088/0031-9155/49/1/004] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A simple, time efficient, analytical model incorporating heterogeneities and body dimensions around a point 192Ir source is generalized for accurate dosimetry around commercially available 192Ir brachytherapy sources. The generalized model was verified in dosimetry of a clinical 192Ir high dose rate prostate monotherapy application, involving 16 catheters and 83 source dwell positions, through comparison with corresponding treatment planning system data. The computational time efficiency and accuracy of the proposed model allowed the assessment of the impact that uncertainties in source dwell positions and catheter reconstruction may have on dose distributions, and how these could potentially affect the clinical outcome. Results revealed that a 0.1 cm catheter reconstruction uncertainty and a 0.15 cm source position uncertainty along the catheter lead to a dose uncertainty of less than 2% for doses lower than 200% of the prescribed dose, reaching up to 5% for points lying in close proximity to the catheters. These uncertainties were found to have no impact (less than 1%) on dose volume histogram results of both the planning target volume and the urethra. A catheter reconstruction uncertainty as high as 0.2 cm results in a dose uncertainty greater than 2%, reaching up to 9%, only for points inside the 150% contour. However, even in this case, the impact on dose volume histogram calculations is less than 3%.
Collapse
Affiliation(s)
- E Pantelis
- Nuclear and Particle Physics Section, Physics Department, University of Athens, Panepistimioupolis, Ilisia, 157 71 Athens, Greece
| | | | | | | | | | | | | |
Collapse
|
14
|
Anagnostopoulos G, Baltas D, Geretschlaeger A, Martin T, Papagiannis P, Tselis N, Zamboglou N. In vivo thermoluminescence dosimetry dose verification of transperineal 192Ir high-dose-rate brachytherapy using CT-based planning for the treatment of prostate cancer. Int J Radiat Oncol Biol Phys 2003; 57:1183-91. [PMID: 14575851 DOI: 10.1016/s0360-3016(03)00762-4] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
PURPOSE To evaluate the potential of in vivo thermoluminescence dosimetry to estimate the accuracy of dose delivery in conformal high-dose-rate brachytherapy of prostate cancer. METHODS AND MATERIALS A total of 50 LiF, TLD-100 cylindrical rods were calibrated in the dose range of interest and used as a batch for all fractions. Fourteen dosimeters for every treatment fraction were loaded in a plastic 4F catheter that was fixed in either one of the 6F needles implanted for treatment purposes or in an extra needle implanted after consulting with the patient. The 6F needles were placed either close to the urethra or in the vicinity of the median posterior wall of the prostate. Initial results are presented for 18 treatment fractions in 5 patients and compared to corresponding data calculated using the commercial treatment planning system used for the planning of the treatments based on CT images acquired postimplantation. RESULTS The maximum observed mean difference between planned and delivered dose within a single treatment fraction was 8.57% +/- 2.61% (root mean square [RMS] errors from 4.03% to 9.73%). Corresponding values obtained after averaging results over all fractions of a patient were 6.88% +/- 4.93% (RMS errors from 4.82% to 7.32%). Experimental results of each fraction corresponding to the same patient point were found to agree within experimental uncertainties. CONCLUSIONS Experimental results indicate that the proposed method is feasible for dose verification purposes and suggest that dose delivery in transperineal high-dose-rate brachytherapy after CT-based planning can be of acceptable accuracy.
Collapse
Affiliation(s)
- G Anagnostopoulos
- Department of Medical Physics and Engineering, Strahlenklinik, Klinikum Offenbach, Offenbach, Germany
| | | | | | | | | | | | | |
Collapse
|
15
|
Kipouros P, Papagiannis P, Sakelliou L, Karaiskos P, Sandilos P, Baras P, Seimenis I, Kozicki M, Anagnostopoulos G, Baltas D. 3D dose verification in 192Ir HDR prostate monotherapy using polymer gels and MRI. Med Phys 2003; 30:2031-9. [PMID: 12945969 DOI: 10.1118/1.1590437] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
VIPAR polymer gels and 3D MRI techniques were evaluated for their ability to provide experimental verification of 3D dose distributions in a simulation of a 192Ir prostate monotherapy clinical application. A real clinical treatment plan was utilized, generated by post-irradiation, CT based calculations derived from Plato BPS and Swift treatment planning systems. The simulated treatment plan involved the use of 10 catheters and 39 source positions within a glass vessel of appropriate dimensions, homogeneously filled with the VIPAR gel. 3D high resolution MR scanning of the gel produced T2 relaxation time maps, from which 3D dose distributions were derived via an appropriate calibration procedure. Results were compared to corresponding dose distributions obtained from the Plato and Swift treatment planning systems. Quantitative comparison, on a point by point basis, was based on user adopted acceptance criteria of 5% dose-difference and 3 mm distance-to-agreement. Significant deviations between experimental and calculated dose distributions were found for doses lower than 50% due to the reduced dose resolution of the method in the low dose, low dose gradient region. Measurement errors were observed at 1.0-1.5 mm around each catheter due to MR imaging susceptibility artifacts. For most remaining points the acceptance criteria were fulfilled. Systematic offsets of the order of 1-2 mm, observed between measured and corresponding calculated isocontours at specific segments, are attributed to the 1 mm uncertainty in catheter reconstruction and 1 mm uncertainty in the alignment of the MR and CT imaging planes.
Collapse
Affiliation(s)
- P Kipouros
- Nuclear and Particle Physics Section, Physics Department, University of Athens, Panepistimioupolis, Ilisia 157 71, Athens, Greece.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
16
|
Milickovic N, Lahanas M, Papagiannopoulo M, Zamboglou N, Baltas D. Multiobjective anatomy-based dose optimization for HDR-brachytherapy with constraint free deterministic algorithms. Phys Med Biol 2002; 47:2263-80. [PMID: 12164586 DOI: 10.1088/0031-9155/47/13/306] [Citation(s) in RCA: 41] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In high dose rate (HDR) brachytherapy, conventional dose optimization algorithms consider multiple objectives in the form of an aggregate function that transforms the multiobjective problem into a single-objective problem. As a result, there is a loss of information on the available alternative possible solutions. This method assumes that the treatment planner exactly understands the correlation between competing objectives and knows the physical constraints. This knowledge is provided by the Pareto trade-off set obtained by single-objective optimization algorithms with a repeated optimization with different importance vectors. A mapping technique avoids non-feasible solutions with negative dwell weights and allows the use of constraint free gradient-based deterministic algorithms. We compare various such algorithms and methods which could improve their performance. This finally allows us to generate a large number of solutions in a few minutes. We use objectives expressed in terms of dose variances obtained from a few hundred sampling points in the planning target volume (PTV) and in organs at risk (OAR). We compare two- to four-dimensional Pareto fronts obtained with the deterministic algorithms and with a fast-simulated annealing algorithm. For PTV-based objectives, due to the convex objective functions, the obtained solutions are global optimal. If OARs are included, then the solutions found are also global optimal, although local minima may be present as suggested.
Collapse
Affiliation(s)
- N Milickovic
- Department of Medical Physics and Engineering, Strahlenklinik, Klinikum Offenbach, Offenbach am Main, Germany.
| | | | | | | | | |
Collapse
|
17
|
Giannouli S, Baltas D, Milickovic N, Lahanas M, Kolotas C, Zamboglou N, Uzunoglu N. Autoactivation of source dwell positions for HDR brachytherapy treatment planning. Med Phys 2000; 27:2517-20. [PMID: 11128303 DOI: 10.1118/1.1315315] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
The most accurate classical dose optimization algorithms in HDR brachytherapy strongly depend on an appropriate selection of source dwell positions which fulfill user-defined geometrical boundary conditions which are relative to patient anatomy. Most anatomical situations, such as for prostate and head and neck tumors, are complex and can require geometries with 5-15 catheters with 48 possible dwell positions per catheter depending on the tumor volume. The manual selection of dwell positions using visual checks by trial and error is very time consuming. This can only be improved by the use of a technique which automatically recognizes and selects the optimum dwell positions for each catheter. We have developed an algorithm, termed an autoactivation algorithm, which improves implant planning by providing a facility for the necessary automatic recognition of HDR source dwell positions.
Collapse
Affiliation(s)
- S Giannouli
- Department of Medical Physics and Engineering, Strahlenklinik, Klinikum Offenbach, Germany.
| | | | | | | | | | | | | |
Collapse
|
18
|
Milickovic N, Baltast D, Giannouli S, Lahanas M, Zamboglou N. CT imaging based digitally reconstructed radiographs and their application in brachytherapy. Phys Med Biol 2000; 45:2787-800. [PMID: 11049172 DOI: 10.1088/0031-9155/45/10/305] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
The aim of our study was to develop an algorithm to simulate the digitally reconstructed radiograph (DRR) calculation process for different beam qualities (photon energies) in the range 50 keV to 12 MeV. This was achieved using volumetric anatomical data for the patient obtained from three-dimensional diagnostic CT images. These DRR images can be used in three-dimensional treatment planning for external beam radiotherapy as well as for brachytherapy in the same way as conventional radiographic films. The advantages of using such DRRs in modern 3D brachytherapy treatment planning are shown. A number of tools are described, illustrating that the application of DRRs in brachytherapy is very useful.
Collapse
Affiliation(s)
- N Milickovic
- Department of Medical Physics & Engineering, Strahlenklinik, Klinikum Offenbach, Germany.
| | | | | | | | | |
Collapse
|