D Nazareth

Optimization of beam angles for intensity modulated radiation therapy treatment planning using genetic algorithm on a distributed computing platform.

Planning intensity modulated radiation therapy (IMRT) treatment involves selection of several angle parameters as well as specification of structures and constraints employed in the optimization process. Including these parameters in the combinatorial search space vastly increases the computational burden, and therefore the parameter selection is normally performed manually by a clinician, based on clinical experience. We have investigated the use of a genetic algorithm (GA) and distributed-computing platform to optimize the gantry angle parameters and provide insight into additional structures, which may be necessary, in the dose optimization process to produce optimal IMRT treatment plans. For an IMRT prostate patient, we produced the first generation of 40 samples, each of five gantry angles, by selecting from a uniform random distribution, subject to certain adjacency and opposition constraints. Dose optimization was performed by distributing the 40-plan workload over several machines running a commercial treatment planning system. A score was assigned to each resulting plan, based on how well it satisfied clinically-relevant constraints. The second generation of 40 samples was produced by combining the highest-scoring samples using techniques of crossover and mutation. The process was repeated until the sixth generation, and the results compared with a clinical (equally-spaced) gantry angle configuration. In the sixth generation, 34 of the 40 GA samples achieved better scores than the clinical plan, with the best plan showing an improvement of 84%. Moreover, the resulting configuration of beam angles tended to cluster toward the patients sides, indicating where the inclusion of additional structures in the dose optimization process may avoid dose hot spots. Additional parameter selection in IMRT leads to a large-scale computational problem. We have demonstrated that the GA combined with a distributed-computing platform can be applied to optimize gantry angle selection within a reasonable amount of time.

Solving Beam-Angle Selection and Dose Optimization Simultaneously via High-Throughput Computing

We provide a framework for integrating two stages of radiation treatment planning (RTP): beam-angle selection (BAS) and dose optimization (DO). The framework is applied to both classical three-dimensional conformal radiotherapy and advanced intensity-modulated radiation therapy. Automated BAS and improved dose distribution are achieved within the framework. A metaheuristic approach, nested partitions, is applied. Alternative BAS and DO algorithms or commercial RTP software and clinical experience can be embedded within the framework to provide new methods for warm starts and evaluations of the quality of beam-angle set samples. Computational efficiency is achieved by utilizing high-throughput computing via the Condor system. Computational results show that our framework has led to a significant improvement in terms of solution quality and delivery time compared with current clinical practice.

Effect of Ultrasound Probe on Dose Delivery During Real-time Ultrasound-Guided Tumor Tracking

Ultrasound is a noninvasive and less costly modality for real-time imaging of soft tissues. It has the capability of tracking soft tissue at levels of submillimeter precision even in the presence of radiation beams. The effect of a transducer on radiation dose is not fully known. The best imaging location for an ultrasound transducer happens to coincide with the path of an anterior-posterior beam in intensity modulated radiation therapy (IMRT). This study indicates a significant change in dose when this juxtaposition occurs. If the anterior-posterior beam is avoided in IMRT planning, however, the effect of the transducer on radiotherapy is found to be negligible

Vessel size measurements in angiograms: Manual measurements

Vessel size measurement is perhaps the most often performed quantitative analysis in diagnostic and interventional angiography. Although automated vessel sizing techniques are generally considered to have good accuracy and precision, we have observed that clinicians rarely use these techniques in standard clinical practice, choosing to indicate the edges of vessels and catheters to determine sizes and calibrate magnifications, i.e., manual measurements. Thus, we undertook an investigation of the accuracy and precision of vessel sizes calculated from manually indicated edges of vessels. Manual measurements were performed by three neuroradiologists and three physicists. Vessel sizes ranged from 0.1-3.0 mm in simulation studies and 0.3-6.4 mm in phantom studies. Simulation resolution functions had full-widths-at-half-maximum (FWHM) ranging from 0.0 to 0.5 mm. Phantom studies were performed with 4.5 in., 6 in., 9 in., and 12 in. image intensifier modes, magnification factor = 1, with and without zooming. The accuracy and reproducibility of the measurements ranged from 0.1 to 0.2 mm, depending on vessel size, resolution, and pixel size, and zoom. These results indicate that manual measurements may have accuracies comparable to automated techniques for vessels with sizes greater than 1 mm, but that automated techniques which take into account the resolution function should be used for vessels with sizes smaller than 1 mm.

Quantification of incidental mediastinal and hilar irradiation delivered during definitive stereotactic body radiation therapy for peripheral non-small cell lung cancer

To determine the amount of incidental radiation dose received by the mediastinal and hilar nodes for patients with non-small cell lung cancer (NSCLC) treated with stereotactic body radiation therapy (SBRT). Fifty consecutive patients with NSCLC, treated using an SBRT technique, were identified. Of these patients, 38 had a prescription dose of 60 Gy in 20-Gy fractions and were eligible for analysis. For each patient, ipsilateral upper (level 2) and lower (level 4) paratracheal, and hilar (level 10) nodal regions were contoured on the planning computed tomography (CT) images. Using the clinical treatment plan, dose and volume calculations were performed retrospectively for each nodal region. SBRT to upper lobe tumors resulted in an average total ipsilateral mean dose of between 5.2 and 7.8 Gy for the most proximal paratracheal nodal stations (2R and 4R for right upper lobe lesions, 2L and 4L for left upper lobe lesions). SBRT to lower lobe tumors resulted in an average total ipsilateral mean dose of between 15.6 and 21.5 Gy for the most proximal hilar nodal stations (10R for right lower lobe lesions, 10 l for left lower lobe lesions). Doses to more distal nodes were substantially lower than 5 Gy. The often substantial incidental irradiation, delivered during SBRT for peripheral NSCLC of the lower lobes to the most proximal hilar lymph nodes may be therapeutic for low-volume, subclinical nodal disease. Treatment of peripheral upper lobe lung tumors delivers less incidental irradiation to the paratracheal lymph nodes with lower likelihood of therapeutic benefit.

Level-set segmentation of pulmonary nodules in megavolt electronic portal images using a CT prior

PURPOSE Pulmonary nodules present unique problems during radiation treatment due to nodule position uncertainty that is caused by respiration. The radiation field has to be enlarged to account for nodule motion during treatment. The purpose of this work is to provide a method of locating a pulmonary nodule in a megavolt portal image that can be used to reduce the internal target volume (ITV) during radiation therapy. A reduction in the ITV would result in a decrease in radiation toxicity to healthy tissue. METHODS Eight patients with nonsmall cell lung cancer were used in this study. CT scans that include the pulmonary nodule were captured with a GE Healthcare LightSpeed RT 16 scanner. Megavolt portal images were acquired with a Varian Trilogy unit equipped with an AS1000 electronic portal imaging device. The nodule localization method uses grayscale morphological filtering and level-set segmentation with a prior. The treatment-time portion of the algorithm is implemented on a graphical processing unit. RESULTS The method was retrospectively tested on eight cases that include a total of 151 megavolt portal image frames. The method reduced the nodule position uncertainty by an average of 40% for seven out of the eight cases. The treatment phase portion of the method has a subsecond execution time that makes it suitable for near-real-time nodule localization. CONCLUSIONS A method was developed to localize a pulmonary nodule in a megavolt portal image. The method uses the characteristics of the nodule in a prior CT scan to enhance the nodule in the portal image and to identify the nodule region by level-set segmentation. In a retrospective study, the method reduced the nodule position uncertainty by an average of 40% for seven out of the eight cases studied.

Applying graphics processor units to Monte Carlo dose calculation in radiation therapy

We investigate the potential in using of using a graphics processor unit (GPU) for Monte-Carlo (MC)-based radiation dose calculations. The percent depth dose (PDD) of photons in a medium with known absorption and scattering coefficients is computed using a MC simulation running on both a standard CPU and a GPU. We demonstrate that the GPUs capability for massive parallel processing provides a significant acceleration in the MC calculation, and offers a significant advantage for distributed stochastic simulations on a single computer. Harnessing this potential of GPUs will help in the early adoption of MC for routine planning in a clinical environment.

Region of interest (ROI) microangiography: imager development

A new high spatial resolution micro-angiographic camera will enable routine viewing within a region of interest of detailed vascular structure unable to be seen with current full field of view (FOV) angiographic detectors. Such details include perforator vessels, vessel contractility or compliance, and condition and location of 50 micron or smaller stent wires. Although the basic CsI(Tl) phosphor-optical taper-CCD design of the new ROI micro-angiographic camera is essentially the same as that of the pre-clinical prototype, many of the physical parameters are much improved. The FOV is 5 cm X 5 cm vs. the previous 1 cm X 1 cm; the phosphor thickness is 350 - 400 micron vs. the previous 100 micron; the taper ratio is now 1.8 rather than 3.0 (2.8X improvement in light collection). The pixel size is either 25 or 50 micron. Additionally, detector noise may now be carefully considered in the camera design as may mechanical supporting mechanisms, methods to synchronize image acquisition with exposure and the effects of other physical factors such as exposure parameters, tube loading, focal spot size and geometric unsharpness. It is expected that this new capability should allow improved treatments and further development of smaller interventional devices and catheter delivery systems.

Comparison of hybrid volumetric modulated arc therapy (VMAT) technique and double arc VMAT technique in the treatment of prostate cancer

Abstract Background. Volumetric modulated arc therapy (VMAT) has quickly become accepted as standard of care for the treatment of prostate cancer based on studies showing it is able to provide faster delivery with adequate target coverage and reduced monitor units while maintaining organ at risk (OAR) sparing. This study aims to demonstrate the potential to increase dose conformality with increased planner control and OAR sparing using a hybrid treatment technique compared to VMAT. Methods. Eleven patients having been previously treated for prostate cancer with VMAT techniques were replanned with a hybrid technique on Varian Treatment Planning System. Multiple static IMRT fields (2 to 3) were planned initially based on critical OAR to reduce dose but provide some planning treatment volume (PTV) coverage. This was used as a base dose plan to provide 30-35% coverage for a single arc VMAT plan. Results. The clinical VMAT plan was used as a control for the purposes of comparison. Average of all OAR sparing between the hybrid technique and VMAT showed the hybrid plan delivering less dose in almost all cases except for V80 of the bladder and maximum dose to right femoral head. PTV coverage was superior with the VMAT technique. Monitor unit differences varied, with the hybrid plan able to deliver fewer units 37% of the time, similar results 18% of the time, and higher units 45% of the time. On average, the hybrid plan delivered 10% more monitor units. Conclusions. The hybrid plan can be delivered in a single gantry rotation combining aspects of VMAT with regions of dynamic intensity modulated radiation therapy (IMRT) within the treatment arc.

Image-Guided Treatment Planning and Dosimetry for Interstitial Photodynamic Therapy

The treatment scenario for photodynamic therapy (PDT) can vary between patients. We present a method that combines clinical images, radiation therapy planning contours, and Monte Carlo simulation in order to design individualized PDT treatment plans.