EMVision: An electromagnetic tracking guidance system for accurate placement of interstitial brachytherapy applicators.

PURPOSE: To design, construct, and evaluate a system for image-guided placement of brachytherapy applicators using electromagnetic tracking (EMT) technology for improved procedure quality and efficiency. The system, named EMVision, provides a three-dimensional display of planned needle sites, visibility of the anatomy and needle position during placement, and reference tracking to account for generator or target anatomy shifts. METHODS: The EMVision EMT guidance system registers CT and EMT reference frames using an automated point-based algorithm according to tandem/ovoid locations. Subsequent needle placement and digitization is guided with an intraluminal EMT sensor and user interface that displays the needle's location in axial, coronal, sagittal, and 3D-volumetric CT views. The interface can overlay contours, preplanned needles, and provide 6DOF compensation for tandem movement during needle placement/manipulation. EMVision accuracy was evaluated in phantom and human cadaver by comparing EMVision's DICOM needle positions with a ground-truth, postimplant CT. RESULTS: Proof of concept was demonstrated for EMT-assisted placement of brachytherapy needles. EMVision accuracy in phantom (mean ± standard deviation) on a brachysuite CT table was 0.76 ± 0.13 mm for needle tips placed up to 75 mm from the tandem/ovoids and 0.52 ± 0.27 mm for needle shafts at distances up to 100 mm from the tandem/ovoids. Performance in human cadaver was similar, with tip and shaft accuracies of 0.77 ± 0.14 mm and 0.40 ± 0.21 mm, respectively. CONCLUSION: EMVision provides sub-millimeter accuracy for the placement of brachytherapy needles without repeated or continuous imaging. The technology can be used to reduce brachytherapy procedure times, improve the correspondence between intended and actual needle positions, or decrease the trainee learning curve.

Color VISION for improved ultrasound visualization of brachytherapy needles.

BACKGROUND: High dose rate brachytherapy is commonly used in the treatment of prostate cancer. Treatment planning is often performed under transrectal ultrasound (US) guidance, but brachytherapy needles can be challenging to digitize due to the presence of poor US conspicuity and imaging artifacts. The plan accuracy and quality, however, are dependent on the proper visualization of the needles with millimeter accuracy. PURPOSE: This work describes a technique for generating a color overlay of needle locations atop the grayscale US image. Prototype devices were developed to produce vibrations in the brachytherapy needles that generate a high contrast color Doppler (CD) signal that highlights the needle locations with superior contrast and reduced artifacts. Denoted by the acronym color VISION (Vibrationally Induced Shimmering for Identifying an Object's Nature), the technology has the potential to improve applicator conspicuity and facilitate automated applicator digitization. METHODS: Three prototype vibrational devices with frequencies between 200-450 Hz were designed in-house and evaluated with needle implants in a phantom and cadaveric male pelvis using: (1) an actuator attached to the front of a prostate needle template; (2) an actuator attached to the top of the needle template; and (3) a hand-held actuator with a stylet, inserted directly into a needle's inner lumen. Acquired images were postprocessed in MATLAB to evaluate the potential for automated digitization. RESULTS: All prototype devices produced localized shimmering in implanted brachytherapy needles in both the axial and sagittal planes. The template mounted actuators provided better vibrational coupling and ease of operation than the stylet prototype. The Michelson contrast, or visibility, of the shimmering CD signal was 100% compared with ≤40% for B-mode imaging of a single needle. Proof-of-principle for automated applicator digitization using only the CD signal was demonstrated. CONCLUSIONS: The color VISION prototype devices successfully coupled mechanical vibrations into brachytherapy needles to generate US CD shimmering and accurately highlight brachytherapy needle locations. The high contrast and natively registered signal are promising for future work to automate the needle digitization and provide a real-time visual overlay of the applicator on the B-mode US image.

A script-enabled interactive checklist document for efficient management of electronic devices in a busy multimodality radiotherapy clinic.

PURPOSE: Develop an efficient, interactive, and instructive checklist document for the management of implanted electronic medical devices in a multimodality radiotherapy clinic. METHODS: The built-in scripting and interactivity of a popular commercial word processor was used to develop an interactive document that changes the information presented to the user based on drop-down selections. The interactivity and scripting were compatible with the radiation oncology information system (ROIS) which allows the document to be accessible by all team members and serve as a permanent record in a patient's electronic chart. RESULTS: The final interactive document, which was clinically deployed after beta testing with a group consisting of nurses and medical physicists, presents information and action plans to the user based on multiple departmental medical device decision trees that are specific to the combination of device, treatment modality, rhythm-pacing dependence for cardiac devices, and distance from the device to the treatment volume. CONCLUSION: A script-enabled interactive document was developed for a busy multimodality clinic, condensing multiple comprehensive departmental guidelines spanning multiple device types and treatment modalities into a single interactive checklist accessible within the ROIS. Given the wide accessibility of the commercial word processor, this approach could be adopted by other clinics to streamline their own respective workflows.

A new way to visualize prostate brachytherapy needles using ultrasound color Doppler and needle surface modifications.

PURPOSE: Suboptimal ultrasound conspicuity of the brachytherapy applicator can lead to inaccurate image reconstructions of the applicator resulting in decreased tumor control or increased normal tissue dose. This feasibility study aims to improve ultrasound conspicuity of high-dose rate (HDR) brachytherapy needles by modifying the surface of the needles to produce a color Doppler twinkling signature. MATERIALS AND METHODS: Surface modifications of standard 17-gauge titanium HDR brachytherapy needles included laser-scribing, application of polymethyl methacrylate (PMMA), and coating with a commercially available echogenic coating. Laser-scribing was performed with variable widths (0.1-1 mm) and depths (10-100 µm). The echogenic coating was applied with 3 different thicknesses (27, 40, and 64 µm). Unmodified and modified needles were imaged under B-mode and color Doppler ultrasound in phantom and cadaver, and the signal strength was recorded. RESULTS: Laser-scribed, PMMA-coated, and echogenic-coated brachytherapy needles produced a twinkling signature along the needle shaft on color Doppler ultrasound. Twinkling was observed with laser-scribe depths >20 µm and widths >0.1 mm and from echogenic coatings 40 µm and 64 µm thick. Twinkling was not observed with unmodified needles. The twinkling signature had a spectral composition with a uniform magnitude between the velocities of 2 to 16 cm/s. CONCLUSIONS: Color Doppler ultrasound of surface-modified brachytherapy applicators may improve applicator conspicuity aiding applicator placement and digitization. HDR brachytherapy needles may be modified to produce the twinkling signature via laser-scribing, PMMA rings, or applying an echogenic coating.

Improving ultrasound-based brachytherapy needle conspicuity by applying an echogenic coating.

BACKGROUND: Applicator conspicuity in ultrasound-guided brachytherapy procedures is commonly impaired by imaging artifacts or non-ideal imaging geometry, which can slow down applicator position digitization and increase the geometric uncertainty of the delivered dose distribution. PURPOSE: The purpose of this study was to improve the conspicuity of high-dose rate (HDR) brachytherapy needles under B-mode ultrasound imaging by applying an echogenic surface coating. Our hypothesis was that an echogenic coating would reduce artifacts and improve needle visualization within regions of signal degradation. METHODS: In this study, 17-gauge, 25-cm long titanium HDR brachytherapy needles were coated with acoustically reflective microspheres over a 2.5 cm region starting from the needle tip. Three coating thicknesses (27 µm, 40 µm, 64 µm) were compared against an uncoated control needle. The coated and uncoated needles were imaged using B-mode ultrasound in a tissue-equivalent prostate phantom and in a cadaverous male pelvis using a transrectal probe. Needle conspicuity was assessed under multiple conditions: a single needle implant, an implant with multiple needles between the probe and the needle of interest, and an angled needle implant. All images were assessed qualitatively for needle conspicuity and the presence of artifacts and quantitatively using grey-scale image intensity values. RESULTS: The 64 µm echogenic coating reduced the magnitude of reverberation artifacts by 31 ± 14% and comet tail artifacts by 40%-70%. The echogenic coating also improved needle contrast, measured by the relative differences in signal intensity compared with the adjacent environment, when needles were angled up to 30° with respect to the transducer probe in the cadaver. The improvements in conspicuity and artifact reduction increased with increasing coating thickness. The performance of the needles coated with the 64 µm thickness was qualitatively superior and yielded high-contrast, well-circumscribed signals in the cadaverous male pelvis, even under situations where a needle was acoustically shadowed by multiple other needles. CONCLUSIONS: An echogenic surface coating reduced imaging artifacts and improved needle conspicuity under realistic clinical conditions for ultrasound-based prostate or gynecological brachytherapy. The improved conspicuity has the potential to improve the efficiency of needle placement and the accuracy of needle position digitization during brachytherapy procedures.

A fast GPU-accelerated Monte Carlo engine for calculation of MLC-collimated electron fields.

BACKGROUND: Although intensity-modulated radiation therapy and volumetric arc therapy have revolutionized photon external beam therapies, the technological advances associated with electron beam therapy have fallen behind. Modern linear accelerators contain technologies that would allow for more advanced forms of electron treatments, such as beam collimation, using the conventional photon multi-leaf collimator (MLC); however, no commercial solutions exist that calculate dose from such beam delivery modes. Additionally, for clinical adoption to occur, dose calculation times would need to be on par with that of modern dose calculation algorithms. PURPOSE: This work developed a graphics processing unit (GPU)-accelerated Monte Carlo (MC) engine incorporating the Varian TrueBeam linac head geometry for a rapid calculation of electron beams collimated using the conventional photon MLC. METHODS: A compute unified device architecture framework was created for the following: (1) transport of electrons and photons through the linac head geometry, considering multiple scattering, Bremsstrahlung, Møller, Compton, and pair production interactions; (2) electron and photon propagation through the CT geometry, considering all interactions plus the photoelectric effect; and (3) secondary particle cascades through the linac head and within the CT geometry. The linac head collimating geometry was modeled according to the specifications provided by the vendor, who also provided phase-space files. The MC was benchmarked against EGSnrc/DOSXYZnrc/GEANT by simulating individual interactions with simple geometries, pencil, and square beam dose calculations in various phantoms. MC-calculated dose distributions for MLC and jaw-collimated electron fields were compared to measurements in a water phantom and with radiochromic film. RESULTS: Pencil and square beam dose distributions are in good agreement with DOSXYZnrc. Angular and spatial distributions for multiple scattering and secondary particle production in thin slab geometries are in good agreement with EGSnrc and GEANT. Dose profiles for MLC and jaw-collimated 6-20-MeV electron beams showed an average absolute difference of 1.1 and 1.9 mm for the FWHM and 80%-20% penumbra from measured profiles. Percent depth doses showed differences of <5% for as compared to measurement. The computation time on an NVIDIA Tesla V100 card was 2.5 min to achieve a dose uncertainty of <1%, which is ∼300 times faster than published results in a similar geometry using a single-CPU core. CONCLUSIONS: The GPU-based MC can quickly calculate dose for electron fields collimated using the conventional photon MLC. The fast calculation times will allow for a rapid calculation of electron fields for mixed photon and electron particle therapy.

Evaluation of SrAl2 O4 :Eu, Dy phosphor for potential applications in thermoluminescent dosimetry.

PURPOSE: To evaluate the use of commercial-grade strontium aluminate phosphorescent powder as a thermoluminescent (TL) dosimeter for clinical radiotherapy beams. MATERIALS AND METHOD: Commercially available Eu2+ , Dy3+ co-doped strontium aluminate powder (SrAl2 O4 :Eu, Dy) was annealed and then irradiated using 20 × 20 cm2 field size, with 6-MV (PDD10  = 70.7) and 18-MV (PDD10  = 79.4) photon beams and and 9-MeV (R50  = 3.6), 15 MeV (R50  = 5.9) and 18-MeV (R50  = 7.2) electron beams. To calibrate the relationship between the TL readings and the irradiated doses, TL glow curves were acquired for doses up to 600 cGy at all beam energies. For the percentage depth dose (PDD) measurement, the SrAl2 O4 :Eu, Dy powder was sandwiched by solid water phantoms, with varying thickness of solid water placed above to determine the depth. PDDs were measured at four representative depths and compared against the commissioning depth dose data for each beam energy. RESULTS: Linear dose response models of doses up to 200 cGy were created for all beam energies. Superlinearity was observed with doses greater than 200 cGy. The PDD measurement acquired experimentally agrees well with the commissioning data of the medical linear accelerator. Trapping parameters such as order of kinetics, activation energy and frequency factor have been obtained via TL glow curve analysis. CONCLUSION: The linear dose response demonstrates that SrAl2 O4 :Eu, Dy is a potential TLD dosimeter for both electron beams and photon beams at different beam energies. The PDD measurements further support its potential use in quality assurance and radiation dosimetry.

Use of Reduced Activity Seeds in Breast Radioactive Seed Localization.

Radioactive seed localization procedures, using I seeds of typical activity 3.7 MBq and higher, are performed to localize nonpalpable lesions in the breast for surgical excision and pathology analysis. This study evaluated the use and dosimetry of I seeds of activity <3.7 MBq in radioactive seed localization procedures through retrospective health record review, Monte Carlo simulation, and experimental detection. An average seed strength at the time of specimen excision of 2.48 ± 0.629 MBq was used in 295 radioactive seed localization procedures at Gundersen Health System in La Crosse, Wisconsin, US. The average explanted seed activity served as a basis for Monte Carlo simulation of an I IsoAid Advantage seed embedded in soft tissue, which scored the dose deposited to soft tissue. Tabulated values of the dose to postsurgical residual tissue as a function of explanted tumor radius were shown and compared with previously published results. Use of seeds of activity from 1.44 to 3.7 MBq at the time of excision did not adversely affect seed detection and excision. The absorbed dose to residual tissue calculated using Monte Carlo was an average of 1.4 times larger than previously published results when scaled to identical seed strengths. This study demonstrates that seeds of activity <3.7 MBq can be used for radioactive seed localization procedures with no loss in efficacy and a benefit of reduced radiation dose to patients. This is important because the estimated radiation dose to residual tissue is approximately 1.4 times higher than previously indicated.

Space-variant deconvolution of Cerenkov light images acquired from a curved surface.

PURPOSE: Cerenkov photons are generated by high-energy radiation used in external beam radiation therapy (EBRT). This study expands upon the Cerenkov light dosimetry formula previously developed to relate an image of Cerenkov photons to the primary beam fluence. Extension of this formulation allows for deconvolution to be performed on images acquired from curved geometries. METHODS: The integral equation, which represented the formation of Cerenkov photon image from an incident high-energy photon beam, was expanded to allow for space-variance of the convolution kernel called as the Cerenkov scatter function (CSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo (MC) particle simulation software was used to obtain the CSF for different incident beam angles. The image of a curved surface was first projected to a flat plane by using a perspective correction method. Then, the planar image was partitioned into small segments (or blocks), where a CSF corresponding to a specific beam incident angle was applied for deconvolution. The block size and the margin around the block were optimized by studying the effects of those parameters on the deconvolution accuracy for a test image. We evaluated three deconvolution techniques: Richardson-Lucy, Blind, and Total Variation minimization (TV/L2) algorithms, to select the most accurate method for the current applications. RESULTS: Analysis of deconvolution algorithms showed that the TV/L2 method provided the most accurate solution to the deconvolution problem for Cerenkov imaging. Optimization of space-variant deconvolution parameters showed that including a margin that is at least 42.9% of the image width provided the most accurate product image. There was no optimal size for the deconvolution area and should be chosen based on the presence of unique CSF kernels within an image. Space-variant deconvolution improved measured field size in Cerenkov photon images by 7.37%, as compared with 1.74% by space-invariant deconvolution. Space-variant deconvolution improved measured penumbra by 99.3%, as compared with 76.7% by space-invariant deconvolution. Space-variant deconvolution introduced artifacts in flat regions of the beam. Artifacts were avoided through selective space-variant deconvolution in only the penumbra region. CONCLUSIONS: Primary photon fluence distributions of a curved surface can be obtained by using space-variant deconvolution methods in Cerenkov light dosimetry. The TV/L2 algorithm is the best method for deconvolution of Cerenkov photon images from an open-field beam derived from either a flat or curved surface. The partition size chosen for space-variant deconvolution should be at least six times the full width at half maximum (FWHM) of the corresponding scatter kernel used in deconvolution. Space-variant deconvolution is necessary if the incident beam angle difference is larger than 6 ∘ between regions of an image.

Characterization of the Cerenkov scatter function: a convolution kernel for Cerenkov light dosimetry.

Cerenkov light is created in clinical applications involving high-energy radiation such as in radiation therapy. There is considerable interest in using Cerenkov light as a means to perform in vivo dosimetry during radiation therapy; however, a better understanding of the light-to-dose relationship is needed. One such method to solve this relationship is that of a deconvolution formulation, which relies on the Cerenkov scatter function (CSF). The CSF describes the creation of Cerenkov photons by a pencil beam of high-energy radiation, and the subsequent scattering that occurs before emission from the irradiated medium surface. This study investigated the dependence of the CSF on common radiation beam parameters (beam energy and incident angle) and the type of irradiated medium. An analytical equation with fitting coefficients of the CSF was obtained for common beam energies in a stratified skin model and optical phantom. Perturbation analysis was performed to investigate the dependence of the deconvolved Cerenkov images on the full-width at half-maximum and amplitude of the CSF. The irradiated material and beam angle had a large impact on the deconvolution process, whereas the beam energy had little effect.

A mathematical deconvolution formulation for superficial dose distribution measurement by Cerenkov light dosimetry.

PURPOSE: Cerenkov photons are created by high-energy radiation beams used for radiation therapy. In this study, we developed a Cerenkov light dosimetry technique to obtain a two-dimensional dose distribution in a superficial region of medium from the images of Cerenkov photons by using a deconvolution method. METHODS: An integral equation was derived to represent the Cerenkov photon image acquired by a camera for a given incident high-energy photon beam by using convolution kernels. Subsequently, an equation relating the planar dose at a depth to a Cerenkov photon image using the well-known relationship between the incident beam fluence and the dose distribution in a medium was obtained. The final equation contained a convolution kernel called the Cerenkov dose scatter function (CDSF). The CDSF function was obtained by deconvolving the Cerenkov scatter function (CSF) with the dose scatter function (DSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo particle simulation software was used to obtain the CSF and DSF. The dose distribution was calculated from the Cerenkov photon intensity data using an iterative deconvolution method with the CDSF. The theoretical formulation was experimentally evaluated by using an optical phantom irradiated by high-energy photon beams. RESULTS: The intensity of the deconvolved Cerenkov photon image showed linear dependence on the dose rate and the photon beam energy. The relative intensity showed a field size dependence similar to the beam output factor. Deconvolved Cerenkov images showed improvement in dose profiles compared with the raw image data. In particular, the deconvolution significantly improved the agreement in the high dose gradient region, such as in the penumbra. Deconvolution with a single iteration was found to provide the most accurate solution of the dose. Two-dimensional dose distributions of the deconvolved Cerenkov images agreed well with the reference distributions for both square fields and a multileaf collimator (MLC) defined, irregularly shaped field. CONCLUSIONS: The proposed technique improved the accuracy of the Cerenkov photon dosimetry in the penumbra region. The results of this study showed initial validation of the deconvolution method for beam profile measurements in a homogeneous media. The new formulation accounted for the physical processes of Cerenkov photon transport in the medium more accurately than previously published methods.