A low-cost phantom design for evaluating spine SABR calculations in the presence of prosthetic vertebral stabilization.

Dose-perturbation characteristics are important to consider during the calculation of radiation therapy protocols for patients who are going to receive high doses that would reach the tolerance limits of the spinal cord [1]. Several studies have investigated dose perturbations introduced by metal implants in close proximity to spine SABR treatments [2-7]. However, there is a lack of work assessing this effect using the RayStation TPS [8]. We present an initial design for a low-cost phantom to evaluate spine stereotactic ablative radiotherapy (SABR) in the presence of prosthetic vertebral stabilization. The phantom is modular, allowing the prosthetic at the centre of the phantom to be removed by exchanging the central block. It also includes space to insert ion chamber and film. The agreement of the RayStation TPS (v8.0B) collapsed cone convolution (CCC) calculation and measurement was determined for phantom versions with and without prosthetic. There was little to no change in the agreement between the measured and calculated dose when introducing metallic hardware. This suggests that our Raystation-based SABR planning approach for patients with spinal hardware meets clinical expectations. Departments without access to anthropomorphic phantoms may find this design useful but should test their phantom design in typical clinical settings to ensure it is robust to real world situations.

Scatter imaging during lung stereotactic body radiation therapy characterized with phantom studies.

By collecting photons scattered out of the therapy beam, scatter imaging creates images of the treated volume. Two phantoms were used to assess the possible application of scatter imaging for markerless tracking of lung tumors during stereotactic body radiation therapy (SBRT) treatment. A scatter-imaging camera was assembled with a CsI flat-panel detector and a 5 mm diameter pinhole collimator. Scatter images were collected during the irradiation of phantoms with megavoltage photons. To assess scatter image quality, spherical phantom lung tumors of 2.1-2.8 cm diameters were placed inside a static, anthropomorphic phantom. To show the efficacy of the technique with a moving target (3 cm diameter), the position of a simulated tumor was tracked in scatter images during sinusoidal motion (15 mm amplitude, 0.25 Hz frequency) in a dynamic lung phantom in open-field, dynamic conformal arc (DCA), and volumetric modulated arc therapy (VMAT) deliveries. Anatomical features are identifiable on static phantom scatter images collected with 10 MU of delivered dose (2.1 cm diameter lung tumor contrast-to-noise ratio of 4.4). The contrast-to-noise ratio increases with tumor size and delivered dose. During dynamic motion, the position of the 3.0 cm diameter lung tumor was identified with a root-mean-square error of 0.8, 1.2, and 2.9 mm for open field (0.3 s frame integration), DCA (0.5 s), and VMAT (0.5 s), respectively. Based on phantom studies, scatter imaging is a potential technique for markerless lung tumor tracking during SBRT without additional imaging dose. Quality scatter images may be collected at low, clinically relevant doses (10 MU). Scatter images are capable of sub-millimeter tracking precision, but modulation decreases accuracy.

Compton scatter imaging: A promising modality for image guidance in lung stereotactic body radiation therapy.

PURPOSE: Lung stereotactic body radiation therapy (SBRT) requires delivering large radiation doses with millimeter accuracy, making image guidance essential. An approach to forming images of patient anatomy from Compton-scattered photons during lung SBRT is presented. METHODS: To investigate the potential of scatter imaging, a pinhole collimator and flat-panel detector are used for spatial localization and detection of photons scattered during external beam therapy using lung SBRT treatment conditions (6 MV FFF beam). MCNP Monte Carlo software is used to develop a model to simulate scatter images. This model is validated by comparing experimental and simulated phantom images. Patient scatter images are then simulated from 4DCT data. RESULTS: Experimental lung tumor phantom images have sufficient contrast-to-noise to visualize the tumor with as few as 10 MU (0.5 s temporal resolution). The relative signal intensity from objects of different composition as well as lung tumor contrast for simulated phantom images agree quantitatively with experimental images, thus validating the Monte Carlo model. Scatter images are shown to display high contrast between different materials (lung, water, bone). Simulated patient images show superior (~double) tumor contrast compared to MV transmission images. CONCLUSIONS: Compton scatter imaging is a promising modality for directly imaging patient anatomy during treatment without additional radiation, and it has the potential to complement existing technologies and aid tumor tracking and lung SBRT image guidance.

Dosimetric effects of saline- versus water-filled balloon applicators for IORT using the model S700 electronic brachytherapy source.

PURPOSE: The Xoft Axxent Electronic Brachytherapy System (Xoft, Inc., San Jose, CA) is a viable option for intraoperative radiation therapy (IORT) treatment of early-stage breast cancer. The low-energy (50-kVp) X-ray source simplifies shielding and increases relative biological effectiveness but increases dose distribution sensitivity to medium composition. Treatment planning systems typically assume homogenous water for brachytherapy dose calculations, including precalculated atlas plans for Xoft IORT. However, Xoft recommends saline for balloon applicator filling. This study investigates dosimetric differences due to increased effective atomic number (Zeff) for saline (Zeff = 7.56) versus water (Zeff = 7.42). METHODS: Balloon applicator diameters range from 3 to 6 cm. Monte Carlo N-Particle software is used to calculate dose at the surface (Ds) of and 1 cm away (D1cm) from the water-/saline-filled balloon applicator using a single dwell at the applicator center as a simple estimation of the dosimetry and multiple dwells simulating the clinical dose distributions for the atlas plans. RESULTS: Single-dwell plans show a 4.4-6.1% decrease in Ds for the 3- to 6-cm diameter applicators due to the saline. Multidwell plans show similar results: 4.9% and 6.4% Ds decrease, for 4-cm and 6-cm diameter applicators, respectively. For the single-dwell plans, D1cm decreases 3.6-5.2% for the 3- to 6-cm diameter applicators. For the multidwell plans, D1cm decreases 3.3% and 5.3% for the 4-cm and 6-cm applicators, respectively. CONCLUSIONS: The dosimetric effect introduced by saline versus water filling for Xoft balloon applicator-based IORT treatments is ∼5%. Users should be aware of this in the context of both treatment planning and patient outcome studies.

Characterization of Compton-scatter imaging with an analytical simulation method.

By collimating the photons scattered when a megavoltage therapy beam interacts with the patient, a Compton-scatter image may be formed without the delivery of an extra dose. To characterize and assess the potential of the technique, an analytical model for simulating scatter images was developed and validated against Monte Carlo (MC). For three phantoms, the scatter images collected during irradiation with a 6 MV flattening-filter-free therapy beam were simulated. Images, profiles, and spectra were compared for different phantoms and different irradiation angles. The proposed analytical method simulates accurate scatter images up to 1000 times faster than MC. Minor differences between MC and analytical simulated images are attributed to limitations in the isotropic superposition/convolution algorithm used to analytically model multiple-order scattering. For a detector placed at 90° relative to the treatment beam, the simulated scattered photon energy spectrum peaks at 140-220 keV, and 40-50% of the photons are the result of multiple scattering. The high energy photons originate at the beam entrance. Increasing the angle between source and detector increases the average energy of the collected photons and decreases the relative contribution of multiple scattered photons. Multiple scattered photons cause blurring in the image. For an ideal 5 mm diameter pinhole collimator placed 18.5 cm from the isocenter, 10 cGy of deposited dose (2 Hz imaging rate for 1200 MU min-1 treatment delivery) is expected to generate an average 1000 photons per mm2 at the detector. For the considered lung tumor CT phantom, the contrast is high enough to clearly identify the lung tumor in the scatter image. Increasing the treatment beam size perpendicular to the detector plane decreases the contrast, although the scatter subject contrast is expected to be greater than the megavoltage transmission image contrast. With the analytical method, real-time tumor tracking may be possible through comparison of simulated and acquired patient images.

The sensitivity of ArcCHECK-based gamma analysis to manufactured errors in helical tomotherapy radiation delivery.

Three-dimensional measurement arrays are an efficient means of acquiring a distribution of data for patient plan delivery QA. However, the tie between plan integrity and traditional gamma-based analysis of these data are not clear. This study explores the sensitivity of such analysis by creating errors in Helical Tomotherapy delivery and measuring the passing rates with an ArcCHECK cylindrical diode array. Errors were introduced in each of the couch speed, leaf open time, and gantry starting position in increasing magnitude while the resulting gamma passing rates were tabulated. The error size required to degrade the gamma passing rate to 90% or below was on average a 3% change in couch speed, 5° in gantry synchronization, or a 5 ms in leaf closing speed for a 3%/3 mm Van Dyk gamma analysis. This varied with plan type, with prostate plans exhibiting less sensitivity than head and neck plans and with gamma analysis criteria, but in all cases the error magnitudes were large compared to actual machine tolerances. These findings suggest that the sensitivity of ArcCHECK-based gamma analysis to single-mode errors in tomotherapy plans is dependent upon plan and analysis type and at traditional passing thresholds unable to detect small defects in the plan.

Optimization for high-dose-rate brachytherapy of cervical cancer with adaptive simulated annealing and gradient descent.

PURPOSE: To validate an in-house optimization program that uses adaptive simulated annealing (ASA) and gradient descent (GD) algorithms and investigate features of physical dose and generalized equivalent uniform dose (gEUD)-based objective functions in high-dose-rate (HDR) brachytherapy for cervical cancer. METHODS: Eight Syed/Neblett template-based cervical cancer HDR interstitial brachytherapy cases were used for this study. Brachytherapy treatment plans were first generated using inverse planning simulated annealing (IPSA). Using the same dwell positions designated in IPSA, plans were then optimized with both physical dose and gEUD-based objective functions, using both ASA and GD algorithms. Comparisons were made between plans both qualitatively and based on dose-volume parameters, evaluating each optimization method and objective function. A hybrid objective function was also designed and implemented in the in-house program. RESULTS: The ASA plans are higher on bladder V75% and D2cc (p=0.034) and lower on rectum V75% and D2cc (p=0.034) than the IPSA plans. The ASA and GD plans are not significantly different. The gEUD-based plans have higher homogeneity index (p=0.034), lower overdose index (p=0.005), and lower rectum gEUD and normal tissue complication probability (p=0.005) than the physical dose-based plans. The hybrid function can produce a plan with dosimetric parameters between the physical dose-based and gEUD-based plans. The optimized plans with the same objective value and dose-volume histogram could have different dose distributions. CONCLUSIONS: Our optimization program based on ASA and GD algorithms is flexible on objective functions, optimization parameters, and can generate optimized plans comparable with IPSA.

Thermal effusivity: a promising imaging biomarker to predict radiation-induced skin injuries.

An effective screening technology is needed to triage individuals at the time of radiation incidents involving a large population. Three-dimensional thermal tomography is a relatively new development in active thermal imaging technology that produces cross-sectional images based on the subject's ability to transfer heat-thermal effusivity-at the voxel level. This noninvasive imaging modality has been used successfully in nondestructive examination of complex materials; also it has been shown to predict the severity of radiation-induced skin injuries several days before the manifestation of severe moist desquamations or blister formation symptoms in mice at 40 Gy. If these results are confirmed at lower dose levels in human subjects, a thermal tomography imaging device may be an ideal screening tool in radiation emergencies. This imaging method is non-invasive, relatively simple, easily adaptable for field use, and when properly deployed, it will enhance public emergency preparedness for incidents involving unexpected radiation exposure.

Thermal effusivity changes as a precursor to moist desquamation.

Skin toxicity is a ubiquitous side effect in radiotherapy and can be difficult to predict. Moist desquamation in cancer patients can decrease quality of life and occasionally demand unplanned treatment breaks thus worsening outcome. In breast cancer patients, moist desquamation occurs approximately one-third of the time, and while avenues such as intensity-modulated radiation therapy exist to decrease skin side effects, they may be prohibitively expensive to distribute widely. To selectively target patients who are at risk for high skin toxicity, toxicity prediction beyond heuristics is required. This study presents 3D thermal tomography, a translation technology that employs active thermal imaging to map the thermal effusivity of skin. Irradiated mice were imaged throughout reaction development to establish a correlation between effusivity changes and eventual toxicity severity. Female hairless mice (n = 11) were anesthetized and irradiated to 40 Gy in one fraction using a 1 cm Leipzig brachytherapy applicator with an Ir-192 source. After irradiation, thermal imaging was conducted daily with a flash lamp and infrared camera. Effusivity was calculated using custom software and tracked within irradiated and contralateral control regions. Mice were retrospectively grouped into high-grade (moist desquamation present, n = 6) and low-grade (n = 5). All mice showed an increase in the relative average effusivity difference among the treated and control regions between irradiation and peak reaction between 12 and 15 days after irradiation. The high-grade group showed an earlier increase in relative average effusivity difference (mean 1.7 days after irradiation versus 4.4 days after irradiation) than the low-grade group, and had a significantly greater relative average effusivity difference between 2-5 days after irradiation. We concluded that 3D thermal tomography is quick, non-invasive, non-ionizing and exhibited a correlative difference between mice that eventually developed moist desquamation and those that only presented dry desquamation. With further development, it may prove to be a useful tool in the clinic for differentiating patients who require preventative measures to reduce skin toxicity.

Comparison of computed tomography scout based reference point localization to conventional film and axial computed tomography.

Identification of source positions after implantation is an important step in brachytherapy planning. Reconstruction is traditionally performed from films taken by conventional simulators, but these are gradually being replaced in the clinic by computed tomography (CT) simulators. The present study explored the use of a scout image-based reconstruction algorithm that replaces the use of traditional film, while exhibiting low sensitivity to metal-induced artifacts that can appear in 3D CT methods. In addition, the accuracy of an in-house graphical software implementation of scout-based reconstruction was compared with seed location reconstructions for 2 phantoms by conventional simulator and CT measurements. One phantom was constructed using a planar fixed grid of 1.5-mm diameter ball bearings (BBs) with 40-mm spacing. The second was a Fletcher-Suit applicator embedded in Styrofoam (Dow Chemical Co., Midland, MI) with one 3.2-mm-diameter BB inserted into each of 6 surrounding holes. Conventional simulator, kilovoltage CT (kVCT), megavoltage CT, and scout-based methods were evaluated by their ability to calculate the distance between seeds (40 mm for the fixed grid, 30-120 mm in Fletcher-Suit). All methods were able to reconstruct the fixed grid distances with an average deviation of <1%. The worst single deviations (approximately 6%) were exhibited in the 2 volumetric CT methods. In the Fletcher-Suit phantom, the intermodality agreement was within approximately 3%, with the conventional sim measuring marginally larger distances, with kVCT the smallest. All of the established reconstruction methods exhibited similar abilities to detect the distances between BBs. The 3D CT-based methods, with lower axial resolution, showed more variation, particularly with the smaller BBs. With a software implementation, scout-based reconstruction is an appealing approach because it simplifies data acquisition over film-based reconstruction without requiring any specialized equipment and does not carry risk of misreads caused by artifacts.

Updating a 3-D vertebral body finite element model using 2-D images.

In osteoporotic patients, vertebral strength is often evaluated in the clinical setting using bone densitometry methods, such as dual energy X-ray absorptiometry. Finite element models based on computed tomography (CT) have been shown to outperform such methods in predicting fracture strength, but repetitive use of CT scans may be impractical due to cost, availability, and radiation exposure. We propose a method of generating a vertebral model at an initial time point which can subsequently be updated using several digital radiographs by using an algebraic reconstruction technique (ART) to find the three-dimensional density distribution. The effectiveness of the algorithm was measured by comparison of the error of the reconstructed model to the error incurred by using the outdated model. Application of the ART was able to reduce density errors from 30% to under 7% and to reduce errors in calculated vertebral stiffness from 125% to under 10%. This preliminary study demonstrates that the method is valid and could potentially improve fracture risk diagnostics drastically.

A hierarchical approach to finite element modeling of the human spine.

Finite element modeling has become an important part of biomechanics research in the last 30 years. In particular, virtual experimentation on the spine bypasses several limitations on mechanical testing: any load and constraint may be applied to a finite element model, and experiments may be repeated with slight alterations on the same specimen. Because computer power is still limited, each virtual model is designed with a particular scale in mind. To analyze spinal motion, a motion segment consisting of two or more vertebrae and the connective tissue is modeled. For vertebral strength, a single vertebra, with or without posterior elements, is analyzed. To understand the biomechanics of the cancellous bone and to examine bone remodeling, the individual trabeculae comprising the cancellous bone are examined. These three areas represent the bulk of current research. Scientists also explore whole-spine mechanics and ultrastructural dynamics, but these suffer greatly from a lack of physiological data on the observed phenomena. To a large extent, the spinal physiological loading scenario and the load distribution at the different structural levels are unknown. As in vivo measurement techniques and computational power grow, so will the usefulness of finite element modeling.