Real-time catheter tracking for high-dose-rate prostate brachytherapy using an electromagnetic 3D-guidance device: a preliminary performance study

US-guided EM tracked system for HDR brachytherapy: A first in-men randomized study for whole prostate treatment

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.

Image-guided Raman spectroscopy navigation system to improve transperineal prostate cancer detection. Part 1: Raman spectroscopy fiber-optics system and in situ tissue characterization

SIGNIFICANCE

The diagnosis of prostate cancer (PCa) and focal treatment by brachytherapy are limited by the lack of precise intraoperative information to target tumors during biopsy collection and radiation seed placement. Image-guidance techniques could improve the safety and diagnostic yield of biopsy collection as well as increase the efficacy of radiotherapy.

AIM

To estimate the accuracy of PCa detection using in situ Raman spectroscopy (RS) in a pilot in-human clinical study and assess biochemical differences between in vivo and ex vivo measurements.

APPROACH

A new miniature RS fiber-optics system equipped with an electromagnetic (EM) tracker was guided by trans-rectal ultrasound-guided imaging, fused with preoperative magnetic resonance imaging to acquire 49 spectra in situ (in vivo) from 18 PCa patients. In addition, 179 spectra were acquired ex vivo in fresh prostate samples from 14 patients who underwent radical prostatectomy. Two machine-learning models were trained to discriminate cancer from normal prostate tissue from both in situ and ex vivo datasets.

RESULTS

A support vector machine (SVM) model was trained on the in situ dataset and its performance was evaluated using leave-one-patient-out cross validation from 28 normal prostate measurements and 21 in-tumor measurements. The model performed at 86% sensitivity and 72% specificity. Similarly, an SVM model was trained with the ex vivo dataset from 152 normal prostate measurements and 27 tumor measurements showing reduced cancer detection performance mostly attributable to spatial registration inaccuracies between probe measurements and histology assessment. A qualitative comparison between in situ and ex vivo measurements demonstrated a one-to-one correspondence and similar ratios between the main Raman bands (e.g., amide I-II bands, phenylalanine).

CONCLUSIONS

PCa detection can be achieved using RS and machine learning models for image-guidance applications using in situ measurements during prostate biopsy procedures.

Performance of an enhanced afterloader with electromagnetic tracking capabilities for channel reconstruction and error detection

PURPOSE

To assess catheter reconstruction and error detection performance of an afterloader (Elekta Brachytherapy, Veenendaal, The Netherlands) equipped with electromagnetic (EM) tracking capabilities.

MATERIALS/METHODS

The Flexitron research unit used was equipped with a special check cable integrating an EM sensor (NDI Aurora V3) that enables tracking and reconstruction capability. The reconstructions of a 24-cm long catheter were performed using two methods: continuous fixed-speed check cable backward stepping (at 1, 2.5, 5, 10, 25 and 50 cm/s) and stepping through each dwell position every 1 mm. The ability of the system to differentiate between two closely located (parallel) catheters was investigated by connecting catheters to the afterloader and moving it from its axis with an increment of 1 mm. A robotic arm (Meca500, Mecademic, Montreal) with an accuracy of 0.01 mm was used to move the catheter between each reconstruction. Reconstructions were obtained with a locally weighted scatterplot smoothing algorithm. To quantify the reconstruction accuracy, distances between two catheters were computed along the reconstruction track with a 5 mm step. The reconstructions of curve catheter paths were assessed through parallel and perpendicular phantom configuration to the EM field generator. Indexer length and lateral errors were simulated and a ROC analysis was made.

RESULTS

Using a 50 cm/s check cable speed does not allow for accurate reconstructions. A slower check cable speed results in better reconstruction performance and smaller standard deviations. At 1 cm/s, a catheter can be shifted laterally down to 1 mm and all paths can be uniquely identified. The optimum operating distance from the field generator (50 to 300 mm) resulted in a lower absolute mean deviation from the expected value (0.2 ± 0.1 mm) versus being positioned on the edge of the electromagnetic sensitive detection volume (0.6 ±0.3 mm). The reconstructions of curved catheters with a check cable speed under 5 cm/s gave a 0.8 mm ±0.3 mm error, or better. All indexer and lateral shifts of 1 mm were detected with a check cable speed of 2.5 cm/s or lower.

CONCLUSIONS

The EM-equipped Flexitron afterloader is able to track and reconstruct catheters with high accuracy. A speed under 5 cm/s is recommended for straight and curved catheter reconstructions. It allows catheter identification down to 1 mm inter-catheter distance shift. The check cable can also be used to detect common shift errors.

Commissioning of an intra-operative US guided prostate HDR system integrating an EM tracking technology

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.

Laser-Generated Leaky Acoustic Wave Imaging for Interventional Guidewire Guidance

Ultrasound (US) is widely used to visualize both tissue and the positions of surgical instruments in real time during surgery. Previously we proposed a new method to exploit US imaging and laser-generated leaky acoustic waves (LAWs) for needle visualization. Although successful, that method only detects the position of a needle tip, with the location of the entire needle deduced from knowing that the needle is straight. The purpose of the current study was to develop a beamforming-based method for the direct visualization of objects. The approach can be applied to objects with arbitrary shapes, such as the guidewires that are commonly used in interventional guidance. With this method, illumination by a short laser pulse generates photoacoustic waves at the top of the guidewire that propagate down its metal surface. These waves then leak into the surrounding tissue, which can be detected by a US array transducer. The time of flight consists of two parts: 1) the propagation time of the guided waves on the guidewire and 2) the propagation time of the US that leaks into the tissue. In principle, an image of the guidewire can be formed based on array beamforming by taking the propagation time on the metal into consideration. Furthermore, we introduced directional filtering and a matched filter to compress the dispersion signal associated with long propagation times. The results showed that guidewires could be detected at depths of at least 70 mm. The maximum detectable angle was 56.3°. LAW imaging with a 1268-mm-long guidewire was also demonstrated. The proposed method has considerable potential in new clinical applications.

A review of brachytherapy physical phantoms developed over the last 20 years: clinical purpose and future requirements

Within the brachytherapy community, many phantoms are constructed in-house, and less commercial development is observed as compared to the field of external beam. Computational or virtual phantom design has seen considerable growth; however, physical phantoms are beneficial for brachytherapy, in which quality is dependent on physical processes, such as accuracy of source placement. Focusing on the design of physical phantoms, this review paper presents a summary of brachytherapy specific phantoms in published journal articles over the last twenty years (January 1, 2000 – December 31, 2019). The papers were analyzed and tabulated by their primary clinical purpose, which was deduced from their associated publications.

A substantial body of work has been published on phantom designs from the brachytherapy community, but a standardized method of reporting technical aspects of the phantoms is lacking. In-house phantom development demonstrates an increasing interest in magnetic resonance (MR) tissue mimicking materials, which is not yet reflected in commercial phantoms available for brachytherapy. The evaluation of phantom design provides insight into the way, in which brachytherapy practice has changed over time, and demonstrates the customised and broad nature of treatments offered.

Positional and angular tracking of HDR 192 Ir source for brachytherapy quality assurance using radiochromic film dosimetry

PURPOSE

To quantify and verify the dosimetric impact of high-dose rate (HDR) source positional uncertainty in brachytherapy, and to introduce a model for three-dimensional (3D) position tracking of the HDR source based on a two-dimensional (2D) measurement. This model has been utilized for the development of a comprehensive source quality assurance (QA) method using radiochromic film (RCF) dosimetry including assessment of different digitization uncertainties.

METHODS

An algorithm was developed and verified to generate 2D dose maps of the mHDR-V2 ¹⁹² Ir source (Elekta, Veenendaal, Netherlands) based on the AAPM TG-43 formalism. The limits of the dosimetric error associated with source (0.9 mm diameter) positional uncertainty were evaluated and experimentally verified with EBT3 film measurements for 6F (2.0 mm diameter) and 4F (1.3 mm diameter) size catheters at the surface (4F, 6F) and 10 mm further (4F only). To quantify this uncertainty, a source tracking model was developed to incorporate the unique geometric features of all isodose lines (IDLs) within any given 2D dose map away from the source. The tracking model normalized the dose map to its maximum, then quantified the IDLs using blob analysis based on features such as area, perimeter, weighted centroid, elliptic orientation, and circularity. The Pearson correlation coefficients (PCCs) between these features and source coordinates (x, y, z, θ_y , θ_z ) were calculated. To experimentally verify the accuracy of the tracking model, EBT3 film pieces were positioned within a Solid Water® (SW) phantom above and below the source and they were exposed simultaneously.

RESULTS

The maximum measured dosimetric variations on the 6F and 4F catheter surfaces were 39.8% and 36.1%, respectively. At 10 mm further, the variation reduced to 2.6% for the 4F catheter which is in agreement with the calculations. The source center (x, y) was strongly correlated with the low IDL-weighted centroid (PCC = 0.99), while the distance to source (z) was correlated with the IDL areas (PCC = 0.96) and perimeters (PCC = 0.99). The source orientation θ_y was correlated with the difference between high and low IDL-weighted centroids (PCC = 0.98), while θ_z was correlated with the elliptic orientation of the 60-90% IDLs (PCC = 0.97) for a maximum distance of z = 5 mm. Beyond 5 mm, IDL circularity was significant, therefore limiting the determination of θ_z (PCC ≤ 0.48). The measured positional errors from the film sets above and below the source indicated a source position at the bottom of the catheter (-0.24 ± 0.07 mm).

CONCLUSIONS

Isodose line features of a 2D dose map away from the HDR source can reveal its spatial coordinates. RCF was shown to be a suitable dosimeter for source tracking and dosimetry. This technique offers a novel source QA method and has the potential to be used for QA of commercial and customized applicators.

Electromagnetic (EM) catheter path tracking in ultrasound-guided brachytherapy of the breast

PURPOSE

To evaluate a novel navigation system for breast brachytherapy, based on ultrasound (US)-guided catheter needle implantations followed by electromagnetic (EM) tracking of catheter paths.

METHODS

Breast phantoms were produced, containing US-visible tumors. Ultrasound was used to localize the tumor pose and volume within the phantom, followed by planning an optimal catheter pattern through the tumor using navigation software. An electromagnetic (EM)-tracked catheter needle was used to insert the catheters in the desired pattern. The inserted catheters were visualized on a post-implant CT, serving as ground truth. Electromagnetic (EM) tracking and reconstruction of the inserted catheter paths were performed by pulling a flexible EM guidewire through each catheter, performed in two clinical brachytherapy suites. The accuracy of EM catheter tracking was evaluated by calculating the Hausdorff distance between the EM-tracked and CT-based catheter paths. The accuracy and clinical feasibility of EM catheter tracking were also evaluated in three breast cancer patients, performed in a separate experiment room.

RESULTS

A total of 71 catheter needles were implanted into 12 phantoms using US guidance and EM navigation, in an average ± SD time of 8.1 ± 2.9 min. The accuracy of EM catheter tracking was dependent on the brachytherapy suite: 2.0 ± 1.2 mm in suite 1 and 0.6 ± 0.2 mm in suite 2. EM catheter tracking was successfully performed in three breast brachytherapy patients. Catheter tracking typically took less than 5 min and had an average accuracy of 1.7 ± 0.3 mm.

CONCLUSION

Our preliminary results show a potential role for US guidance and EM needle navigation for implantation of catheters for breast brachytherapy. EM catheter tracking can accurately assess the implant geometry in breast brachytherapy patients. This methodology has the potential to evaluate catheter positions directly after the implantation and during the several fractions of the treatment.

Feasibility study of navigated endoscopy for the placement of high dose rate brachytherapy applicators in the esophagus and lung

PURPOSE

To evaluate the electromagnetic (EM) tracking of endoscopes and applicators as a method of positioning a high dose rate (HDR) luminal applicator.

METHOD

An anatomical phantom consisting of a rigid trachea and flexible esophagus was used to compare applicator placement measurements using EM tracking vs the traditional method using two-dimensional (2D) fluoroscopy and surface skin markers. The phantom included a tumor in the esophagus and several pairs of optically visible points inside the lumen that were used to simulate proximal and distal ends of tumors of varying lengths. The esophagus tumor and lung points were visible on a computed tomography (CT) image of the phantom, which was used as ground truth for the measurements. The EM tracking system was registered to the CT image using fiducial markers. A flexible endoscope was tracked using the EM system and the locations of the proximal and distal ends of the tumor identified and this position recorded. An EM-tracked applicator was then inserted and positioned relative to the tumor markings. The applicator path was mapped using the EM tracking. The gross tumor length (GTL) and the distance between the first dwell position and distal edge of tumor (offset) were measured using the EM tracking and 2D fluoroscopy methods and compared to the same measurements on the CT image.

RESULTS

The errors in GTL using EM tracking were on average -0.5 ± 1.7 mm and 0.7 ± 3.6 mm for esophagus and lung measurements, similar to errors measured using the 2D fluoroscopy method of -0.9 ± 1.2 mm and 3.4 ± 4.4 mm. Offset measurements were slightly larger while using EM tracking relative to the fluoroscopy method but these were not statistically significant.

CONCLUSIONS

Electromagnetic tracking for placement of lumen applicators is feasible and accurate. Tracking of the endoscope that is used to identify the proximal and distal ends of the tumor and of the applicator during insertion generates accurate three-dimensional measurements of the applicator path, GTL and offset. Guiding the placement of intraluminal applicators using EM navigation is potentially attractive for cases with complex insertions, such as those with nonlinear paths or multiple applicator insertions.

Dynamic Modulated Brachytherapy (DMBT) Balloon Applicator for Accelerated Partial Breast Irradiation

PURPOSE

To propose a novel high-dose-rate brachytherapy applicator for balloon-based dynamic modulated brachytherapy (DMBT) for accelerated partial breast irradiation (APBI) and to demonstrate its dosimetric advantage compared to the widely used Contura applicator.

METHODS AND MATERIALS

The DMBT balloon device consists of a fixed central channel enabling real-time, in vivo dosimetry and an outer motion-dynamic, adjustable-radius channel capable of moving to any angular position within the balloon. This design allows placement of dwell positions anywhere within the balloon volume, guaranteeing optimal placement and generation of the applicator and treatment plan, respectively. Thirteen clinical treatment plans for patients with early-stage breast cancer receiving APBI after lumpectomy using Contura were retrospectively obtained under institutional review board approval. New treatment plans were created by replacing the Contura with the DMBT device. DMBT plans were limited to 4 angular positions and an outer channel radius of 1.5 cm. The new plans were optimized to limit dose to ribs and skin while maintaining target coverage similar to that of the clinical plan.

RESULTS

Similar target coverage was obtained for the DMBT plans compared with clinical Contura plans. Across all patients the mean (standard deviation) reductions in D0.1 cc to the ribs and skin were 6.70% (6.28%) and 5.13% (6.54%), respectively. A threshold separation distance between the balloon surface and the organ at risk (OAR), below which dosimetric changes of greater than 5% were obtained, was observed to be 12 mm for ribs and skin. When both OARs were far from the balloon, DMBT plans were of similar quality to Contura plans, as expected.

CONCLUSIONS

This study demonstrates the superior ability of the APBI DMBT applicator to spare OARs while achieving target coverage comparable to current treatment plans, especially when in close proximity. The DMBT balloon may enable new modes of dynamic high-dose-rate treatment delivery and allow for ultrahypofractionated dose regimens to be safely used.

Technical Note: A comparison of point set registration methods for electromagnetic tracking

PURPOSE

High dose rate brachytherapy applies intense and destructive radiation. A treatment plan defines radiation source dwell positions to avoid irradiating healthy tissue. The study discusses methods to quantify any positional changes of source locations along the various treatment sessions.

METHODS

Electromagnetic tracking (EMT) localizes the radiation source during the treatment sessions. But in each session the relative position of the patient relative to the filed generator is changed. Hence, the measured dwell point sets need to be registered onto each other to render them comparable. Two point set registration techniques are compared: a probabilistic method called coherent point drift (CPD) and a multidimensional scaling (MDS) technique.

RESULTS

Both enable using EMT without external registration and achieve very similar results with respect to dwell position determination of the radiation source. Still MDS achieves smaller grand average deviations (CPD-rPSR: MD = 2.55 mm, MDS-PSR: MD = 2.15 mm) between subsequent dwell position determinations, which also show less variance (CPD-rPSR: IQR = 4 mm, MDS-PSR: IQR = 3 mm). Furthermore, MDS is not based on approximations and does not need an iterative procedure to track sensor positions inside the implanted catheters.

CONCLUSION

Although both methods achieve similar results, MDS is to be preferred over rigid CPD while nonrigid CPD is unsuitable as it does not preserve topology.

Review of strategies for MRI based reconstruction of endocavitary and interstitial applicators in brachytherapy of cervical cancer

Brachytherapy plays an essential role in the curative intent management of locally advanced cervical cancer. The introduction of the magnetic resonance (MR) as a preferred image modality and the development of new type of applicators with interstitial components have further improved its benefits. The aim of this work is to review the current status of one important aspect in the cervix cancer brachytherapy procedure, namely catheter reconstruction. MR compatible intracavitary and interstitial applicators are described. Considerations about the use of MR imaging (MRI) regarding appropriate strategies for applicator reconstruction, technical requirements, MR sequences, patient preparation and applicator commissioning are included. It is recommendable to perform the reconstruction process in the same image study employed by the physician for contouring, that is, T2 weighted (T2W) sequences. Nevertheless, a clear identification of the source path inside the catheters and the applicators is a challenge when using exclusively T2W sequences. For the intracavitary component of the implant, sometimes the catheters may be filled with some substance that produces a high intensity signal on MRI. However, this strategy is not feasible for plastic tubes or titanium needles, which, moreover, induce magnetic susceptibility artifacts. In these situations, the use of applicator libraries available in the treatment planning system (TPS) is useful, since they not only include accurate geometrical models of the intracavitary applicators, but also recent developments have made possible the implementation of the interstitial component. Another strategy to improve the reconstruction process is based on the incorporation of MR markers, such as small pellets, to be used as anchor points. Many institutions employ computed tomography (CT) as a supporting image modality. The registration of CT and MR image sets should be carefully performed, and its uncertainty previously assessed. Besides, an important research work is being carried out regarding the use of ultrasound and electromagnetic tracking technologies.

MR-based source localization for MR-guided HDR brachytherapy

For the purpose of MR-guided high-dose-rate (HDR) brachytherapy, a method for real-time localization of an HDR brachytherapy source was developed, which requires high spatial and temporal resolutions. MR-based localization of an HDR source serves two main aims. First, it enables real-time treatment verification by determination of the HDR source positions during treatment. Second, when using a dummy source, MR-based source localization provides an automatic detection of the source dwell positions after catheter insertion, allowing elimination of the catheter reconstruction procedure. Localization of the HDR source was conducted by simulation of the MR artifacts, followed by a phase correlation localization algorithm applied to the MR images and the simulated images, to determine the position of the HDR source in the MR images. To increase the temporal resolution of the MR acquisition, the spatial resolution was decreased, and a subpixel localization operation was introduced. Furthermore, parallel imaging (sensitivity encoding) was applied to further decrease the MR scan time. The localization method was validated by a comparison with CT, and the accuracy and precision were investigated. The results demonstrated that the described method could be used to determine the HDR source position with a high accuracy (0.4-0.6 mm) and a high precision (⩽0.1 mm), at high temporal resolutions (0.15-1.2 s per slice). This would enable real-time treatment verification as well as an automatic detection of the source dwell positions.

Assessment of the implant geometry in fractionated interstitial HDR breast brachytherapy using an electromagnetic tracking system

PURPOSE

During the partial-breast treatment course by interstitial brachytherapy, electromagnetic tracking (EMT) was applied to measure the implant geometry. Implant-geometry variation, choice of reference data, and three registration methods were assessed.

METHODS AND MATERIALS

The implant geometry was measured in 28 patients after catheter implantation (EMT_bed), during CT imaging (EMT_CT), and in each of up to n = 9 treatment fractions (EMT_{F(k), k = 1, 2,… n}). EMT_F(k) were registered to the planned implant reconstruction (CT_plan) by using all dwell positions (DPs), the button centers, or three fiducial sensors on the patient's skin. Variation in implant geometry obtained from EMT_F(k) was assessed for EMT_bed, EMT_CT, and CT_plan.

RESULTS

EMT was used to measure 3932 catheters. A duration of 6.5 ± 1.7 min was needed for each implant measurement (mean, 17 catheters) plus setup of the EMT system. Data registration based on the DP deviated significantly lower than registration on button centers or fiducial sensors. Within a registration group, there was a <0.5-mm difference in the choice of reference data. Using CT_plan as reference for registration, the mean residual distance of DPs on EMT-derived DPs was found at 2.1 ± 1.6 mm (EMT_bed), 1.3 ± 0.9 mm (EMT_CT), and 2.5 ± 1.5 mm (EMT_F(k)).

CONCLUSIONS

EMT can assess the implant geometry in high-dose-rate interstitial brachytherapy breast treatments. EMT_bed, EMT_CT, and CT_plan data can serve as reference for assessment of implant changes.

On the use of multi-dimensional scaling and electromagnetic tracking in high dose rate brachytherapy

High dose rate brachytherapy affords a frequent reassurance of the precise dwell positions of the radiation source. The current investigation proposes a multi-dimensional scaling transformation of both data sets to estimate dwell positions without any external reference. Furthermore, the related distributions of dwell positions are characterized by uni-or bi-modal heavy-tailed distributions. The latter are well represented by α-stable distributions. The newly proposed data analysis provides dwell position deviations with high accuracy, and, furthermore, offers a convenient visualization of the actual shapes of the catheters which guide the radiation source during the treatment.

Electromagnetic tracking (EMT) technology for improved treatment quality assurance in interstitial brachytherapy

Electromagnetic Tracking (EMT) is a novel technique for error detection and quality assurance (QA) in interstitial high dose rate brachytherapy (HDR-iBT). The purpose of this study is to provide a concept for data acquisition developed as part of a clinical evaluation study on the use of EMT during interstitial treatment of breast cancer patients. The stability, accuracy, and precision of EMT-determined dwell positions were quantified. Dwell position reconstruction based on EMT was investigated on CT table, HDR table and PDR bed to examine the influence on precision and accuracy in a typical clinical workflow. All investigations were performed using a precise PMMA phantom. The track of catheters inserted in that phantom was measured by manually inserting a 5 degree of freedom (DoF) sensor while recording the position of three 6DoF fiducial sensors on the phantom surface to correct motion influences. From the corrected data, dwell positions were reconstructed along the catheter's track. The accuracy of the EMT-determined dwell positions was quantified by the residual distances to reference dwell positions after using a rigid registration. Precision and accuracy were investigated for different phantom-table and sensor-field generator (FG) distances. The measured precision of the EMT-determined dwell positions was ≤ 0.28 mm (95th percentile). Stability tests showed a drift of 0.03 mm in the first 20 min of use. Sudden shaking of the FG or (large) metallic objects close to the FG degrade the precision. The accuracy with respect to the reference dwell positions was on all clinical tables < 1 mm at 200 mm FG distance and 120 mm phantom-table distance. Phantom measurements showed that EMT-determined localization of dwell positions in HDR-iBT is stable, precise, and sufficiently accurate for clinical assessment. The presented method may be viable for clinical applications in HDR-iBT, like implant definition, error detection or quantification of uncertainties. Further clinical investigations are needed.

Real-time electromagnetic tracking-based treatment platform for high-dose-rate prostate brachytherapy: Clinical workflows and end-to-end validation

PURPOSE

New technologies were integrated into a novel treatment platform combining electromagnetically (EM) tracked catheters, a 3D ultrasound (3DUS) imaging device, and a new treatment planning system to provide a real-time prostate high-dose-rate (HDR) brachytherapy treatment system. This work defines workflows for offline CT and online 3DUS planning scenarios and preclinical end-to-end validation of the platform.

METHODS AND MATERIALS

The platform is composed of an EM-tracked stylet, a EM-tracked 3DUS probe, and an EM-tracked template guide, all used with the NDI Aurora field generator (NDI, Ontario, Canada). The treatment planning system performs continuous position and angular readings from all three EM sensors into a streamlined environment that allows for (1) contouring; (2) planning; (3) catheter insertion guidance and reconstruction; (4) QA of catheter path and tip position; and (5) exporting to an afterloader. Data were gathered on the times required for the various key steps of the 3DUS-based workflow.

RESULTS

The complete 3DUS-based workflow on 16-catheter implant phantoms took approximately 15 min. This time is expected to increase for actual patients. Plan generation is fast (7.6 ± 2.5s) and the initial catheter reconstruction with updated dose distribution is obtained at no (time) cost as part of the insertion process. Subsequent catheter reconstruction takes on average 10.5 ± 3.1s per catheter, representing less than 3 min for a 16-catheter implant.

CONCLUSIONS

This preclinical study suggests that EM technology could help to significantly streamline real-time US-based high-dose-rate prostate brachytherapy.

Improved electromagnetic tracking for catheter path reconstruction with application in high-dose-rate brachytherapy

PURPOSE

Electromagnetic (EM) catheter tracking has recently been introduced in order to enable prompt and uncomplicated reconstruction of catheter paths in various clinical interventions. However, EM tracking is prone to measurement errors which can compromise the outcome of the procedure. Minimizing catheter tracking errors is therefore paramount to improve the path reconstruction accuracy.

METHODS

An extended Kalman filter (EKF) was employed to combine the nonlinear kinematic model of an EM sensor inside the catheter, with both its position and orientation measurements. The formulation of the kinematic model was based on the nonholonomic motion constraints of the EM sensor inside the catheter. Experimental verification was carried out in a clinical HDR suite. Ten catheters were inserted with mean curvatures varying from 0 to [Formula: see text] in a phantom. A miniaturized Ascension (Burlington, Vermont, USA) trakSTAR EM sensor (model 55) was threaded within each catheter at various speeds ranging from 7.4 to [Formula: see text]. The nonholonomic EKF was applied on the tracking data in order to statistically improve the EM tracking accuracy. A sample reconstruction error was defined at each point as the Euclidean distance between the estimated EM measurement and its corresponding ground truth. A path reconstruction accuracy was defined as the root mean square of the sample reconstruction errors, while the path reconstruction precision was defined as the standard deviation of these sample reconstruction errors. The impacts of sensor velocity and path curvature on the nonholonomic EKF method were determined. Finally, the nonholonomic EKF catheter path reconstructions were compared with the reconstructions provided by the manufacturer's filters under default settings, namely the AC wide notch and the DC adaptive filter.

RESULTS

With a path reconstruction accuracy of 1.9 mm, the nonholonomic EKF surpassed the performance of the manufacturer's filters (2.4 mm) by 21% and the raw EM measurements (3.5 mm) by 46%. Similarly, with a path reconstruction precision of 0.8 mm, the nonholonomic EKF surpassed the performance of the manufacturer's filters (1.0 mm) by 20% and the raw EM measurements (1.7 mm) by 53%. Path reconstruction accuracies did not follow an apparent trend when varying the path curvature and sensor velocity; instead, reconstruction accuracies were predominantly impacted by the position of the EM field transmitter ([Formula: see text]).

CONCLUSION

The advanced nonholonomic EKF is effective in reducing EM measurement errors when reconstructing catheter paths, is robust to path curvature and sensor speed, and runs in real time. Our approach is promising for a plurality of clinical procedures requiring catheter reconstructions, such as cardiovascular interventions, pulmonary applications (Bender et al. in medical image computing and computer-assisted intervention-MICCAI 99. Springer, Berlin, pp 981-989, 1999), and brachytherapy.

Electromagnetic tracking for treatment verification in interstitial brachytherapy

Electromagnetic tracking (EMT) is used in several medical fields to determine the position and orientation of dedicated sensors, e.g., attached to surgical tools. Recently, EMT has been introduced to brachytherapy for implant reconstruction and error detection. The manuscript briefly summarizes the main issues of EMT and error detection in brachytherapy. The potential and complementarity of EMT as treatment verification technology will be discussed in relation to in vivo dosimetry and imaging.

Performance and suitability assessment of a real-time 3D electromagnetic needle tracking system for interstitial brachytherapy

PURPOSE

Accurate insertion and overall needle positioning are key requirements for effective brachytherapy treatments. This work aims at demonstrating the accuracy performance and the suitability of the Aurora(®) V1 Planar Field Generator (PFG) electromagnetic tracking system (EMTS) for real-time treatment assistance in interstitial brachytherapy procedures.

MATERIAL AND METHODS

The system's performance was characterized in two distinct studies. First, in an environment free of EM disturbance, the boundaries of the detection volume of the EMTS were characterized and a tracking error analysis was performed. Secondly, a distortion analysis was conducted as a means of assessing the tracking accuracy performance of the system in the presence of potential EM disturbance generated by the proximity of standard brachytherapy components.

RESULTS

The tracking accuracy experiments showed that positional errors were typically 2 ± 1 mm in a zone restricted to the first 30 cm of the detection volume. However, at the edges of the detection volume, sensor position errors of up to 16 mm were recorded. On the other hand, orientation errors remained low at ± 2° for most of the measurements. The EM distortion analysis showed that the presence of typical brachytherapy components in vicinity of the EMTS had little influence on tracking accuracy. Position errors of less than 1 mm were recorded with all components except with a metallic arm support, which induced a mean absolute error of approximately 1.4 mm when located 10 cm away from the needle sensor.

CONCLUSIONS

The Aurora(®) V1 PFG EMTS possesses a great potential for real-time treatment assistance in general interstitial brachytherapy. In view of our experimental results, we however recommend that the needle axis remains as parallel as possible to the generator surface during treatment and that the tracking zone be restricted to the first 30 cm from the generator surface.

Review of advanced catheter technologies in radiation oncology brachytherapy procedures

The development of new catheter and applicator technologies in recent years has significantly improved treatment accuracy, efficiency, and outcomes in brachytherapy. In this paper, we review these advances, focusing on the performance of catheter imaging and reconstruction techniques in brachytherapy procedures using magnetic resonance images and electromagnetic tracking. The accuracy of catheter reconstruction, imaging artifacts, and other notable properties of plastic and titanium applicators in gynecologic treatments are reviewed. The accuracy, noise performance, and limitations of electromagnetic tracking for catheter reconstruction are discussed. Several newly developed applicators for accelerated partial breast irradiation and gynecologic treatments are also reviewed. New hypofractionated high dose rate treatment schemes in prostate cancer and accelerated partial breast irradiation are presented.

A system to use electromagnetic tracking for the quality assurance of brachytherapy catheter digitization

PURPOSE

To investigate the use of a system using electromagnetic tracking (EMT), post-processing and an error-detection algorithm for detecting errors and resolving uncertainties in high-dose-rate brachytherapy catheter digitization for treatment planning.

METHODS

EMT was used to localize 15 catheters inserted into a phantom using a stepwise acquisition technique. Five distinct acquisition experiments were performed. Noise associated with the acquisition was calculated. The dwell location configuration was extracted from the EMT data. A CT scan of the phantom was performed, and five distinct catheter digitization sessions were performed. No a priori registration of the CT scan coordinate system with the EMT coordinate system was performed. CT-based digitization was automatically extracted from the brachytherapy plan DICOM files (CT), and rigid registration was performed between EMT and CT dwell positions. EMT registration error was characterized in terms of the mean and maximum distance between corresponding EMT and CT dwell positions per catheter. An algorithm for error detection and identification was presented. Three types of errors were systematically simulated: swap of two catheter numbers, partial swap of catheter number identification for parts of the catheters (mix), and catheter-tip shift. Error-detection sensitivity (number of simulated scenarios correctly identified as containing an error/number of simulated scenarios containing an error) and specificity (number of scenarios correctly identified as not containing errors/number of correct scenarios) were calculated. Catheter identification sensitivity (number of catheters correctly identified as erroneous across all scenarios/number of erroneous catheters across all scenarios) and specificity (number of catheters correctly identified as correct across all scenarios/number of correct catheters across all scenarios) were calculated. The mean detected and identified shift was calculated.

RESULTS

The maximum noise ±1 standard deviation associated with the EMT acquisitions was 1.0 ± 0.1 mm, and the mean noise was 0.6 ± 0.1 mm. Registration of all the EMT and CT dwell positions was associated with a mean catheter error of 0.6 ± 0.2 mm, a maximum catheter error of 0.9 ± 0.4 mm, a mean dwell error of 1.0 ± 0.3 mm, and a maximum dwell error of 1.3 ± 0.7 mm. Error detection and catheter identification sensitivity and specificity of 100% were observed for swap, mix and shift (≥2.6 mm for error detection; ≥2.7 mm for catheter identification) errors. A mean detected shift of 1.8 ± 0.4 mm and a mean identified shift of 1.9 ± 0.4 mm were observed.

CONCLUSIONS

Registration of the EMT dwell positions to the CT dwell positions was possible with a residual mean error per catheter of 0.6 ± 0.2 mm and a maximum error for any dwell of 1.3 ± 0.7 mm. These low residual registration errors show that quality assurance of the general characteristics of the catheters and of possible errors affecting one specific dwell position is possible. The sensitivity and specificity of the catheter digitization verification algorithm was 100% for swap and mix errors and for shifts ≥2.6 mm. On average, shifts ≥1.8 mm were detected, and shifts ≥1.9 mm were detected and identified.

Augmented reality assisted laparoscopic partial nephrectomy

Computer assisted navigation is a widely adopted technique in neurosurgery and orthopedics. However, it is rarely used for surgeries on abdominal organs. In this paper, we propose a novel, noninvasive method based on electromagnetic tracking to determine the pose of the kidney. As a clinical use case, we show a complete surgical navigation system for augmented reality assisted laparoscopic partial nephrectomy. Experiments were performed ex vivo on pig kidneys and the evaluation showed an excellent augmented reality alignment error of 2.1 mm ± 1.2 mm.

Electromagnetic tracking for catheter reconstruction in ultrasound-guided high-dose-rate brachytherapy of the prostate

PURPOSE

The accurate delivery of high-dose-rate brachytherapy is dependent on the correct identification of the position and shape of the treatment catheters. In many brachytherapy clinics, transrectal ultrasound (TRUS) imaging is used to identify the catheters. However, manual catheter identification on TRUS images can be time consuming, subjective, and operator dependent because of calcifications and distal shadowing artifacts. We report the use of electromagnetic (EM) tracking technology to map the position and shape of catheters inserted in a tissue-mimicking phantom.

METHODS AND MATERIALS

The accuracy of the EM system was comprehensively quantified using a three-axis robotic system. In addition, EM tracks acquired from catheters in a phantom were compared with catheter positions determined from TRUS and CT images to compare EM system performance to standard clinical imaging modalities. The tracking experiments were performed in a controlled laboratory environment and also in a typical brachytherapy operating room to test for potential EM distortions.

RESULTS

The robotic validation of the EM system yielded a mean accuracy of <0.5 mm for a clinically acceptable field of view in a nondistorting environment. The EM-tracked catheter representations were found to have an accuracy of <1 mm when compared with TRUS- and CT-identified positions, both in the laboratory environment and in the brachytherapy operating room. The achievable accuracy depends to a large extent on the calibration of the TRUS probe, geometry of the tracked devices relative to the EM field generator, and locations of surrounding clinical equipment. To address the issue of variable accuracy, a robust calibration algorithm has been developed and integrated into the workflow. The proposed mapping technique was also found to improve the workflow efficiency of catheter identification.

CONCLUSIONS

The high baseline accuracy of the EM system, the consistent agreement between EM-tracked, TRUS- and CT-identified catheters, and the improved workflow efficiency illustrate the potential value of using EM tracking for catheter mapping in high-dose-rate brachytherapy.

Electromagnetic tracking system with reduced distortion using quadratic excitation

PURPOSE

Electromagnetic tracking systems, frequently used in minimally invasive surgery, are affected by conductive distorters. The influence of conductive distorters on electromagnetic tracking system accuracy can be reduced through magnetic field modifications. This approach was developed and tested.

METHODS

The voltage induced directly by the emitting coil in the sensing coil without additional influence by the conductive distorter depends on the first derivative of the voltage on the emitting coil. The voltage which is induced indirectly by the emitting coil across the conductive distorter in the sensing coil, however, depends on the second derivative of the voltage on the emitting coil. The electromagnetic tracking system takes advantage of this difference by supplying the emitting coil with a quadratic excitation voltage. The method is adaptive relative to the amount of distortion cause by the conductive distorters. This approach is evaluated with an experimental setup of the electromagnetic tracking system.

RESULTS

In vitro testing showed that the maximal error decreased from 10.9 to 3.8 mm when the quadratic voltage was used to excite the emitting coil instead of the sinusoidal voltage. Furthermore, the root mean square error in the proximity of the aluminum disk used as a conductive distorter was reduced from 3.5 to 1.6 mm when the electromagnetic tracking system used the quadratic instead of sinusoidal excitation.

CONCLUSIONS

Electromagnetic tracking with quadratic excitation is immune to the effects of a conductive distorter, especially compared with sinusoidal excitation of the emitting coil. Quadratic excitation of electromagnetic tracking for computer-assisted surgery is promising for clinical applications.