Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “histogram”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An Entropy-Based Test and Development Framework for Uncertainty Modeling in Level-Set Visualizations

We present a simple comparative framework for testing and developing uncertainty modeling in uncertain marching cubes implementations. The selection of a model to represent the probability distribution of uncertain values directly influences the memory use, run time, and accuracy of an uncertainty visualization algorithm. We use an entropy calculation directly on ensemble data to establish an expected result and then compare the entropy from various probability models, including uniform, Gaussian, histogram, and quantile models. Our results verify that models matching the distribution of the ensemble indeed match the entropy. We further show that fewer bins in nonparametric histogram models are more effective whereas large numbers of bins in quantile models approach data accuracy.

Sisneros, Robert↗
DOE OSTI2012

Radiobiology for eye plaque brachytherapy and evaluation of implant duration and radionuclide choice using an objective function

Purpose: Clinical optimization of Collaborative Ocular Melanoma Study (COMS) eye plaque brachytherapy is currently limited to tumor coverage, consensus prescription dosage, and dose calculations to ocular structures. The biologically effective dose (BED) of temporary brachytherapy treatments is a function of both chosen radionuclide R and implant duration T. This study endeavored to evaluate BED delivered to the tumor volume and surrounding ocular structures as a function of plaque position P, prescription dose, R, and T. Methods: Plaque-heterogeneity-corrected dose distributions were generated with MCNP5 for the range of currently available COMS plaques loaded with sources using three available low-energy radionuclides. These physical dose distributions were imported into the PINNACLE{sup 3} treatment planning system using the TG-43 hybrid technique and used to generate dose volume histograms for a T = 7 day implant within a reference eye geometry including the ciliary body, cornea, eyelid, foveola, lacrimal gland, lens, optic disc, optic nerve, retina, and tumor at eight standard treatment positions. The equation of Dale and Jones was employed to create biologically effective dose volume histograms (BEDVHs), allowing for BED volumetric analysis of all ROIs. Isobiologically effective prescription doses were calculated for T = 5 days down to 0.01 days, with BEDVHs subsequently generated for all ROIs using correspondingly reduced prescription doses. Objective functions were created to evaluate the BEDVHs as a function of R and T. These objective functions are mathematically accessible and sufficiently general to be applied to temporary or permanent brachytherapy implants for a variety of disease sites. Results: Reducing T from 7 to 0.01 days for a 10 mm plaque produced an average BED benefit of 26%, 20%, and 17% for {sup 103}Pd, {sup 125}I, and {sup 131}Cs, respectively, for all P; 16 and 22 mm plaque results were more position-dependent. {sup 103}Pd produced a 16%-35% BED benefit over {sup 125}I, whereas {sup 131}Cs produced a 3%-7% BED detriment, independent of P, T, and plaque size. Additionally, corresponding organ at risk physical doses were lowest using {sup 103}Pd in all circumstances. Conclusions: The results suggest that shorter implant durations may correlate with more favorable outcomes compared to 7 day implants when treating small or medium intraocular lesions. The data also indicate that implant duration may be safely reduced if the prescription physical dose is likewise diminished and that {sup 103}Pd offers a substantial radiobiological benefit over {sup 125}I and {sup 131}Cs irrespective of plaque position, implant duration, and tumor size.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Effectiveness of robust optimization in intensity-modulated proton therapy planning for head and neck cancers

Purpose: Intensity-modulated proton therapy (IMPT) is highly sensitive to uncertainties in beam range and patient setup. Conventionally, these uncertainties are dealt using geometrically expanded planning target volume (PTV). In this paper, the authors evaluated a robust optimization method that deals with the uncertainties directly during the spot weight optimization to ensure clinical target volume (CTV) coverage without using PTV. The authors compared the two methods for a population of head and neck (H and N) cancer patients. Methods: Two sets of IMPT plans were generated for 14 H and N cases, one being PTV-based conventionally optimized and the other CTV-based robustly optimized. For the PTV-based conventionally optimized plans, the uncertainties are accounted for by expanding CTV to PTV via margins and delivering the prescribed dose to PTV. For the CTV-based robustly optimized plans, spot weight optimization was guided to reduce the discrepancy in doses under extreme setup and range uncertainties directly, while delivering the prescribed dose to CTV rather than PTV. For each of these plans, the authors calculated dose distributions under various uncertainty settings. The root-mean-square dose (RMSD) for each voxel was computed and the area under the RMSD-volume histogram curves (AUC) was used to relatively compare plan robustness. Data derived from the dose volume histogram in the worst-case and nominal doses were used to evaluate the plan optimality. Then the plan evaluation metrics were averaged over the 14 cases and were compared with two-sided paired t tests. Results: CTV-based robust optimization led to more robust (i.e., smaller AUCs) plans for both targets and organs. Under the worst-case scenario and the nominal scenario, CTV-based robustly optimized plans showed better target coverage (i.e., greater D{sub 95%}), improved dose homogeneity (i.e., smaller D{sub 5%}- D{sub 95%}), and lower or equivalent dose to organs at risk. Conclusions: CTV-based robust optimization provided significantly more robust dose distributions to targets and organs than PTV-based conventional optimization in H and N using IMPT. Eliminating the use of PTV and planning directly based on CTV provided better or equivalent normal tissue sparing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Dose-shaping using targeted sparse optimization

Purpose: Dose volume histograms (DVHs) are common tools in radiation therapy treatment planning to characterize plan quality. As statistical metrics, DVHs provide a compact summary of the underlying plan at the cost of losing spatial information: the same or similar dose-volume histograms can arise from substantially different spatial dose maps. This is exactly the reason why physicians and physicists scrutinize dose maps even after they satisfy all DVH endpoints numerically. However, up to this point, little has been done to control spatial phenomena, such as the spatial distribution of hot spots, which has significant clinical implications. To this end, the authors propose a novel objective function that enables a more direct tradeoff between target coverage, organ-sparing, and planning target volume (PTV) homogeneity, and presents our findings from four prostate cases, a pancreas case, and a head-and-neck case to illustrate the advantages and general applicability of our method.Methods: In designing the energy minimization objective (E{sub tot}{sup sparse}), the authors utilized the following robust cost functions: (1) an asymmetric linear well function to allow differential penalties for underdose, relaxation of prescription dose, and overdose in the PTV; (2) a two-piece linear function to heavily penalize high dose and mildly penalize low and intermediate dose in organs-at risk (OARs); and (3) a total variation energy, i.e., the L{sub 1} norm applied to the first-order approximation of the dose gradient in the PTV. By minimizing a weighted sum of these robust costs, general conformity to dose prescription and dose-gradient prescription is achieved while encouraging prescription violations to follow a Laplace distribution. In contrast, conventional quadratic objectives are associated with a Gaussian distribution of violations, which is less forgiving to large violations of prescription than the Laplace distribution. As a result, the proposed objective E{sub tot}{sup sparse} improves tradeoff between planning goals by 'sacrificing' voxels that have already been violated to improve PTV coverage, PTV homogeneity, and/or OAR-sparing. In doing so, overall plan quality is increased since these large violations only arise if a net reduction in E{sub tot}{sup sparse} occurs as a result. For example, large violations to dose prescription in the PTV in E{sub tot}{sup sparse}-optimized plans will naturally localize to voxels in and around PTV-OAR overlaps where OAR-sparing may be increased without compromising target coverage. The authors compared the results of our method and the corresponding clinical plans using analyses of DVH plots, dose maps, and two quantitative metrics that quantify PTV homogeneity and overdose. These metrics do not penalize underdose since E{sub tot}{sup sparse}-optimized plans were planned such that their target coverage was similar or better than that of the clinical plans. Finally, plan deliverability was assessed with the 2D modulation index.Results: The proposed method was implemented using IBM's CPLEX optimization package (ILOG CPLEX, Sunnyvale, CA) and required 1-4 min to solve with a 12-core Intel i7 processor. In the testing procedure, the authors optimized for several points on the Pareto surface of four 7-field 6MV prostate cases that were optimized for different levels of PTV homogeneity and OAR-sparing. The generated results were compared against each other and the clinical plan by analyzing their DVH plots and dose maps. After developing intuition by planning the four prostate cases, which had relatively few tradeoffs, the authors applied our method to a 7-field 6 MV pancreas case and a 9-field 6MV head-and-neck case to test the potential impact of our method on more challenging cases. The authors found that our formulation: (1) provided excellent flexibility for balancing OAR-sparing with PTV homogeneity; and (2) permitted the dose planner more control over the evolution of the PTV's spatial dose distribution than conventional objective functions. In particular, E{sub tot}{sup sparse}-optimized plans for the pancreas case and head-and-neck case exhibited substantially improved sparing of the spinal cord and parotid glands, respectively, while maintaining or improving sparing for other OARs and markedly improving PTV homogeneity. Plan deliverability for E{sub tot}{sup sparse}-optimized plans was shown to be better than their associated clinical plans, according to the two-dimensional modulation index.Conclusions: These results suggest that our formulation may be used to improve dose-shaping and OAR-sparing for complicated disease sites, such as the pancreas or head and neck. Furthermore, our objective function and constraints are linear and constitute a linear program, which converges to the global minimum quickly, and can be easily implemented in treatment planning software. Thus, the authors expect fast translation of our method to the clinic where it may have a positive impact on plan quality for challenging disease sites.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Novel multimodality segmentation using level sets and Jensen-Rényi divergence

Purpose: Positron emission tomography (PET) is playing an increasing role in radiotherapy treatment planning. However, despite progress, robust algorithms for PET and multimodal image segmentation are still lacking, especially if the algorithm were extended to image-guided and adaptive radiotherapy (IGART). This work presents a novel multimodality segmentation algorithm using the Jensen-Rényi divergence (JRD) to evolve the geometric level set contour. The algorithm offers improved noise tolerance which is particularly applicable to segmentation of regions found in PET and cone-beam computed tomography. Methods: A steepest gradient ascent optimization method is used in conjunction with the JRD and a level set active contour to iteratively evolve a contour to partition an image based on statistical divergence of the intensity histograms. The algorithm is evaluated using PET scans of pharyngolaryngeal squamous cell carcinoma with the corresponding histological reference. The multimodality extension of the algorithm is evaluated using 22 PET/CT scans of patients with lung carcinoma and a physical phantom scanned under varying image quality conditions. Results: The average concordance index (CI) of the JRD segmentation of the PET images was 0.56 with an average classification error of 65%. The segmentation of the lung carcinoma images had a maximum diameter relative error of 63%, 19.5%, and 14.8% when using CT, PET, and combined PET/CT images, respectively. The estimated maximal diameters of the gross tumor volume (GTV) showed a high correlation with the macroscopically determined maximal diameters, with aR{sup 2} value of 0.85 and 0.88 using the PET and PET/CT images, respectively. Results from the physical phantom show that the JRD is more robust to image noise compared to mutual information and region growing. Conclusions: The JRD has shown improved noise tolerance compared to mutual information for the purpose of PET image segmentation. Presented is a flexible framework for multimodal image segmentation that can incorporate a large number of inputs efficiently for IGART.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Multicriteria optimization of the spatial dose distribution

Purpose: Treatment planning for radiation therapy involves trade-offs with respect to different clinical goals. Typically, the dose distribution is evaluated based on few statistics and dose–volume histograms. Particularly for stereotactic treatments, the spatial dose distribution represents further criteria, e.g., when considering the gradient between subregions of volumes of interest. The authors have studied how to consider the spatial dose distribution using a multicriteria optimization approach.Methods: The authors have extended a stepwise multicriteria optimization approach to include criteria with respect to the local dose distribution. Based on a three-dimensional visualization of the dose the authors use a software tool allowing interaction with the dose distribution to map objectives with respect to its shape to a constrained optimization problem. Similarly, conflicting criteria are highlighted and the planner decides if and where to relax the shape of the dose distribution.Results: To demonstrate the potential of spatial multicriteria optimization, the tool was applied to a prostate and meningioma case. For the prostate case, local sparing of the rectal wall and shaping of a boost volume are achieved through local relaxations and while maintaining the remaining dose distribution. For the meningioma, target coverage is improved by compromising low dose conformality toward noncritical structures. A comparison of dose–volume histograms illustrates the importance of spatial information for achieving the trade-offs.Conclusions: The results show that it is possible to consider the location of conflicting criteria during treatment planning. Particularly, it is possible to conserve already achieved goals with respect to the dose distribution, to visualize potential trade-offs, and to relax constraints locally. Hence, the proposed approach facilitates a systematic exploration of the optimal shape of the dose distribution.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Improvement of dose distribution in breast radiotherapy using a reversible transverse magnetic field Linac-MR unit

Purpose: To investigate the improvement in dose distribution in tangential breast radiotherapy using a reversible transverse magnetic field that maintains the same direction of Lorentz force between two fields. The investigation has a potential application in future Linac-MR units. Methods: Computed tomography images of four patients and magnetic fields of 0.25–1.5 Tesla (T) were used for Monte Carlo simulation. Two patients had intact breast while the other two had mastectomy. Simulations of planning and chest wall irradiation were similar to the actual clinical process. The direction of superior-inferior magnetic field for the medial treatment beam was reversed for the lateral beam. Results: For the ipsilateral lung and heart mean doses were reduced by a mean (range) of 45.8% (27.6%–58.6%) and 26.0% (20.2%–38.9%), respectively, depending on various treatment plan setups. The mean V{sub 20} for ipsilateral lung was reduced by 55.0% (43.6%–77.3%). In addition acceptable results were shown after simulation of 0.25 T magnetic field demonstrated in dose-volume reductions of the heart, ipsilateral lung, and noninvolved skin. Conclusions: Applying a reversible magnetic field during breast radiotherapy, not only reduces the dose to the lung and heart but also produces a sharp drop dose volume histogram for planning target volume, because of bending of the path of secondary charged particles toward the chest wall by the Lorentz force. The simulations have shown that use of the magnetic field at 1.5 T is not feasible for clinical applications due to the increase of ipsilateral chest wall skin dose in comparison to the conventional planning while 0.25 T is suitable for all patients due to dose reduction to the chest wall skin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Metallic artifact mitigation and organ-constrained tissue assignment for Monte Carlo calculations of permanent implant lung brachytherapy

Purpose: To investigate methods of generating accurate patient-specific computational phantoms for the Monte Carlo calculation of lung brachytherapy patient dose distributions. Methods: Four metallic artifact mitigation methods are applied to six lung brachytherapy patient computed tomography (CT) images: simple threshold replacement (STR) identifies high CT values in the vicinity of the seeds and replaces them with estimated true values; fan beam virtual sinogram replaces artifact-affected values in a virtual sinogram and performs a filtered back-projection to generate a corrected image; 3D median filter replaces voxel values that differ from the median value in a region of interest surrounding the voxel and then applies a second filter to reduce noise; and a combination of fan beam virtual sinogram and STR. Computational phantoms are generated from artifact-corrected and uncorrected images using several tissue assignment schemes: both lung-contour constrained and unconstrained global schemes are considered. Voxel mass densities are assigned based on voxel CT number or using the nominal tissue mass densities. Dose distributions are calculated using the EGSnrc user-code BrachyDose for{sup 125}I, {sup 103}Pd, and {sup 131}Cs seeds and are compared directly as well as through dose volume histograms and dose metrics for target volumes surrounding surgical sutures. Results: Metallic artifact mitigation techniques vary in ability to reduce artifacts while preserving tissue detail. Notably, images corrected with the fan beam virtual sinogram have reduced artifacts but residual artifacts near sources remain requiring additional use of STR; the 3D median filter removes artifacts but simultaneously removes detail in lung and bone. Doses vary considerably between computational phantoms with the largest differences arising from artifact-affected voxels assigned to bone in the vicinity of the seeds. Consequently, when metallic artifact reduction and constrained tissue assignment within lung contours are employed in generated phantoms, this erroneous assignment is reduced, generally resulting in higher doses. Lung-constrained tissue assignment also results in increased doses in regions of interest due to a reduction in the erroneous assignment of adipose to voxels within lung contours. Differences in dose metrics calculated for different computational phantoms are sensitive to radionuclide photon spectra with the largest differences for{sup 103}Pd seeds and smallest but still considerable differences for {sup 131}Cs seeds. Conclusions: Despite producing differences in CT images, dose metrics calculated using the STR, fan beam + STR, and 3D median filter techniques produce similar dose metrics. Results suggest that the accuracy of dose distributions for permanent implant lung brachytherapy is improved by applying lung-constrained tissue assignment schemes to metallic artifact corrected images.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Model-based dose calculations for COMS eye plaque brachytherapy using an anatomically realistic eye phantom

Purpose : To investigate the effects of the composition and geometry of ocular media and tissues surrounding the eye on dose distributions for COMS eye plaque brachytherapy with{sup 125}I, {sup 103}Pd, or {sup 131}Cs seeds, and to investigate doses to ocular structures. Methods : An anatomically and compositionally realistic voxelized eye model with a medial tumor is developed based on a literature review. Mass energy absorption and attenuation coefficients for ocular media are calculated. Radiation transport and dose deposition are simulated using the EGSnrc Monte Carlo user-code BrachyDose for a fully loaded COMS eye plaque within a water phantom and our full eye model for the three radionuclides. A TG-43 simulation with the same seed configuration in a water phantom neglecting the plaque and interseed effects is also performed. The impact on dose distributions of varying tumor position, as well as tumor and surrounding tissue media is investigated. Each simulation and radionuclide is compared using isodose contours, dose volume histograms for the lens and tumor, maximum, minimum, and average doses to structures of interest, and doses to voxels of interest within the eye. Results : Mass energy absorption and attenuation coefficients of the ocular media differ from those of water by as much as 12% within the 20–30 keV photon energy range. For all radionuclides studied, average doses to the tumor and lens regions in the full eye model differ from those for the plaque in water by 8%–10% and 13%–14%, respectively; the average doses to the tumor and lens regions differ between the full eye model and the TG-43 simulation by 2%–17% and 29%–34%, respectively. Replacing the surrounding tissues in the eye model with water increases the maximum and average doses to the lens by 2% and 3%, respectively. Substituting the tumor medium in the eye model for water, soft tissue, or an alternate melanoma composition affects tumor dose compared to the default eye model simulation by up to 16%. In the full eye model simulations, the average dose to the lens is larger by 7%–9% than the dose to the center of the lens, and the maximum dose to the optic nerve is 17%–22% higher than the dose to the optic disk for all radionuclides. In general, when normalized to the same prescription dose at the tumor apex, doses delivered to all structures of interest in the full eye model are lowest for{sup 103}Pd and highest for {sup 131}Cs, except for the tumor where the average dose is highest for {sup 103}Pd and lowest for {sup 131}Cs. Conclusions : The eye is not radiologically water-equivalent, as doses from simulations of the plaque in the full eye model differ considerably from doses for the plaque in a water phantom and from simulated TG-43 calculated doses. This demonstrates the importance of model-based dose calculations for eye plaque brachytherapy, for which accurate elemental compositions of ocular media are necessary.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

The first clinical implementation of electromagnetic transponder-guided MLC tracking

Purpose: We report on the clinical process, quality assurance, and geometric and dosimetric results of the first clinical implementation of electromagnetic transponder-guided MLC tracking which occurred on 28 November 2013 at the Northern Sydney Cancer Centre. Methods: An electromagnetic transponder-based positioning system (Calypso) was modified to send the target position output to in-house-developed MLC tracking code, which adjusts the leaf positions to optimally align the treatment beam with the real-time target position. Clinical process and quality assurance procedures were developed and performed. The first clinical implementation of electromagnetic transponder-guided MLC tracking was for a prostate cancer patient being treated with dual-arc VMAT (RapidArc). For the first fraction of the first patient treatment of electromagnetic transponder-guided MLC tracking we recorded the in-room time and transponder positions, and performed dose reconstruction to estimate the delivered dose and also the dose received had MLC tracking not been used. Results: The total in-room time was 21 min with 2 min of beam delivery. No additional time was needed for MLC tracking and there were no beam holds. The average prostate position from the initial setup was 1.2 mm, mostly an anterior shift. Dose reconstruction analysis of the delivered dose with MLC tracking showed similar isodose and target dose volume histograms to the planned treatment and a 4.6% increase in the fractional rectal V{sub 60}. Dose reconstruction without motion compensation showed a 30% increase in the fractional rectal V{sub 60} from that planned, even for the small motion. Conclusions: The real-time beam-target correction method, electromagnetic transponder-guided MLC tracking, has been translated to the clinic. This achievement represents a milestone in improving geometric and dosimetric accuracy, and by inference treatment outcomes, in cancer radiotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Impact of dose size in single fraction spatially fractionated (grid) radiotherapy for melanoma

Purpose: To evaluate the impact of dose size in single fraction, spatially fractionated (grid) radiotherapy for selectively killing infiltrated melanoma cancer cells of different tumor sizes, using different radiobiological models. Methods: A Monte Carlo technique was employed to calculate the 3D dose distribution of a commercially available megavoltage grid collimator in a 6 MV beam. The linear-quadratic (LQ) and modified linear quadratic (MLQ) models were used separately to evaluate the therapeutic outcome of a series of single fraction regimens that employed grid therapy to treat both acute and late responding melanomas of varying sizes. The dose prescription point was at the center of the tumor volume. Dose sizes ranging from 1 to 30 Gy at 100% dose line were modeled. Tumors were either touching the skin surface or having their centers at a depth of 3 cm. The equivalent uniform dose (EUD) to the melanoma cells and the therapeutic ratio (TR) were defined by comparing grid therapy with the traditional open debulking field. The clinical outcomes from recent reports were used to verify the authors’ model. Results: Dose profiles at different depths and 3D dose distributions in a series of 3D melanomas treated with grid therapy were obtained. The EUDs and TRs for all sizes of 3D tumors involved at different doses were derived through the LQ and MLQ models, and a practical equation was derived. The EUD was only one fifth of the prescribed dose. The TR was dependent on the prescribed dose and on the LQ parameters of both the interspersed cancer and normal tissue cells. The results from the LQ model were consistent with those of the MLQ model. At 20 Gy, the EUD and TR by the LQ model were 2.8% higher and 1% lower than by the MLQ, while at 10 Gy, the EUD and TR as defined by the LQ model were only 1.4% higher and 0.8% lower, respectively. The dose volume histograms of grid therapy for a 10 cm tumor showed different dosimetric characteristics from those of conventional radiotherapy. A significant portion of the tumor volume received a very large dose in grid therapy, which ensures significant tumor cell killing in these regions. Conversely, some areas received a relatively small dose, thereby sparing interspersed normal cells and increasing radiation tolerance. The radiobiology modeling results indicated that grid therapy could be useful for treating acutely responding melanomas infiltrating radiosensitive normal tissues. The theoretical model predictions were supported by the clinical outcomes. Conclusions: Grid therapy functions by selectively killing infiltrating tumor cells and concomitantly sparing interspersed normal cells. The TR depends on the radiosensitivity of the cell population, dose, tumor size, and location. Because the volumes of very high dose regions are small, the LQ model can be used safely to predict the clinical outcomes of grid therapy. When treating melanomas with a dose of 15 Gy or higher, single fraction grid therapy is clearly advantageous for sparing interspersed normal cells. The existence of a threshold fraction dose, which was found in the authors’ theoretical simulations, was confirmed by clinical observations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Dose impact in radiographic lung injury following lung SBRT: Statistical analysis and geometric interpretation

Purpose: To demonstrate a new method of evaluating dose response of treatment-induced lung radiographic injury post-SBRT (stereotactic body radiotherapy) treatment and the discovery of bimodal dose behavior within clinically identified injury volumes. Methods: Follow-up CT scans at 3, 6, and 12 months were acquired from 24 patients treated with SBRT for stage-1 primary lung cancers or oligometastic lesions. Injury regions in these scans were propagated to the planning CT coordinates by performing deformable registration of the follow-ups to the planning CTs. A bimodal behavior was repeatedly observed from the probability distribution for dose values within the deformed injury regions. Based on a mixture-Gaussian assumption, an Expectation-Maximization (EM) algorithm was used to obtain characteristic parameters for such distribution. Geometric analysis was performed to interpret such parameters and infer the critical dose level that is potentially inductive of post-SBRT lung injury. Results: The Gaussian mixture obtained from the EM algorithm closely approximates the empirical dose histogram within the injury volume with good consistency. The average Kullback-Leibler divergence values between the empirical differential dose volume histogram and the EM-obtained Gaussian mixture distribution were calculated to be 0.069, 0.063, and 0.092 for the 3, 6, and 12 month follow-up groups, respectively. The lower Gaussian component was located at approximately 70% prescription dose (35 Gy) for all three follow-up time points. The higher Gaussian component, contributed by the dose received by planning target volume, was located at around 107% of the prescription dose. Geometrical analysis suggests the mean of the lower Gaussian component, located at 35 Gy, as a possible indicator for a critical dose that induces lung injury after SBRT. Conclusions: An innovative and improved method for analyzing the correspondence between lung radiographic injury and SBRT treatment dose has been demonstrated. Bimodal behavior was observed in the dose distribution of lung injury after SBRT. Novel statistical and geometrical analysis has shown that the systematically quantified low-dose peak at approximately 35 Gy, or 70% prescription dose, is a good indication of a critical dose for injury. The determined critical dose of 35 Gy resembles the critical dose volume limit of 30 Gy for ipsilateral bronchus in RTOG 0618 and results from previous studies. The authors seek to further extend this improved analysis method to a larger cohort to better understand the interpatient variation in radiographic lung injury dose response post-SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-J-93: Parametrisation of Dose to the Mucosa of the Anterior Rectal Wall in Transrectal Ultrasound Guided High-Dose-Rate Brachytherapy of the Prostate

Purpose: In high-dose-rate (HDR) brachytherapy of the prostate, radiation is delivered from a number of radioactive sources which are inserted via catheter into the target volume. The rectal mucosa also receives dose during the treatment, which may lead to late toxicity effects. To allow possible links between rectal dose and toxicity to be investigated, suitable methods of parametrising the rectal dose are needed. Methods: During treatment of a series of 95 patients, anatomy and catheter locations were monitored by transrectal ultrasound, and target volume positions were contoured on the ultrasound scan by the therapist. The anterior rectal mucosal wall was identified by contouring the transrectal ultrasound balloon within the ultrasound scan. Source positions and dwell times, along with the dose delivered to the patient were computed using the Oncentra Prostate treatment planning system (TPS). Data for the series of patients were exported from the TPS in Dicom format, and a series of parametrisation methods were developed in a Matlab environment to assess the rectal dose. Results: Contours of the anterior rectal mucosa were voxelised within Matlab to allow the dose to the rectal mucosa to be analysed directly from the 3D dose grid. Dose parametrisations based on dose-surface (DSH) and dose-line (DLH) histograms were obtained. Both lateral and longitudinal extents of the mucosal dose were parametrised using dose-line histograms in the relevant directions. Conclusion: We have developed a series of dose parametrisations for quantifying the dose to the rectal mucosa during HDR prostate brachytherapy which are suitable for future studies investigating potential associations between mucosal dose and late toxicity effects. The geometry of the transrectal probe standardises the rectal anatomy, making this treatment technique particularly suited to studies of this nature.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-275: Radiobiological Evaluation of Intensity Modulated Radiotherapy Treatment for Locally Advanced Head and Neck Squamous Cell Carcinomas

Purpose: To evaluate the radiobiological outcome of Intensity Modulated Radiotherapy Treatment (IMRT) for locally advanced head and neck squamous cell carcinomas using HART (Histogram Analysis in Radiation Therapy; J Appl Clin Med Phys 11(1): 137–157, 2010) program and compare with the clinical outcomes. Methods: We have treated 20 patients of stage III and IV HNSCC Oropharynx and hypopharynx with accelerated IMRT technique and concurrent chemotherapy. Delineation of tumor and normal tissues were done using Danish Head and Neck Cancer Group (DAHANCA) contouring guidelines and radiotherapy was delivered to a dose of 70Gy in 35 fractions to the primary and involved lymph nodes, 63Gy to intermediate risk areas and 56 Gy to lower risk areas, Monday to Saturday, 6 Days/week using 6 MV Photons with an expected overall treatment time of 6 weeks. The TCP and NTCP's were calculated from the dose-volume histogram (DVH) statistics using the Poisson Statistics (PS) and JT Lyman models respectively and the Resultwas correlated with clinical outcomes of the patients with mean follow up of 24 months. Results: Using HART program, the TCP (0.89± 0.01) of primary tumor and the NTCP for parotids (0.20±0.12), spinal cord (0.05±0.01), esophagus (0.30±0.2), mandible (0.35±0.21), Oral cavity (0.37±0.18), Larynx (0.30±0.15) were estimated and correlated with clinical outcome of the patients. Conclusion: Accelerated IMRT with Chemotherapy is a clinical feasible option in the treatment of locally advanced HNSCC with encouraging initial tumour response and acceptable acute toxicities. The correlation between the clinical outcomes and radiobiological model estimated parameters using HART programs are found to be satisfactory.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-375: Passive Scattering to Pencil-Beam-Scanning Comparison for Medulloblastoma Proton Therapy: LET Distributions and Radiobiological Implications

Purpose: To compare the linear energy transfer (LET) distributions between passive scattering and pencil beam scanning proton radiation therapy techniques for medulloblastoma patients and study the potential radiobiological implications. Methods: A group of medulloblastoma patients, previously treated with passive scattering (PS) proton craniospinal irradiation followed by prosterior fossa or involved field boost, were selected from the patient database of our institution. Using the beam geometry and planning computed tomography (CT) image sets of the original treatment plans, pencil beam scanning (PBS) treatment plans were generated for the cranial treatment for each patient, with average beam spot size of 8mm (sigma in air at isocenter). 3-dimensional dose and LET distributions were calculated by Monte Carlo methods (TOPAS) both for the original passive scattering and new pencil beam scanning treatment plans. LET volume histograms were calculated for the target and OARs and compared for the two delivery methods. Variable RBE weighted dose distributions and volume histograms were also calculated using a variable dose and LET-based model. Results: Better dose conformity was achieved with PBS planning compared to PS, leading to increased dose coverage for the boost target area and decreased average dose to the structures adjacent to it and critical structures outside the whole brain treatment field. LET values for the target were lower for PBS plans. Elevated LET values for OARs close to the boosted target areas were noticed, due to end of range of proton beams falling inside these structures, resulting in higher RBE weighted dose for these structures compared to the clinical RBE value of 1.1. Conclusion: Transitioning from passive scattering to pencil beam scanning proton radiation treatment can be dosimetrically beneficial for medulloblastoma patients. LET–guided treatment planning could contribute to better decision making for these cases, especially for critical structures at close proximity to the boosted target area.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-604: Dosimetric Dependence On the Collimator Angle in Prostate Volumetric Modulated Arc Therapy

Purpose: The purpose of this study is to investigate the dose-volume variations of planning target volume (PTV) and organs-at-risk (OARs) in prostate volumetric modulated arc therapy (VMAT) when using different collimator angles. It is because collimator angle awareness is essential for planner to produce an optimal prostate VMAT plan in a rational time. Methods: Single-arc VMAT plans at different collimator angles (0o, 15o, 30o, 45o, 60o, 75o and 90o) were created systematically using a Harold heterogeneous pelvis phantom. For each change of collimator angle, a new plan was re-optimized for that angle. The prescription dose was 78 Gy per 39 fractions. Conformity index (CI), homogeneity index (HI), gradient index, machine monitor unit, dose-volume histogram, the mean and maximum doses of the PTV were calculated and analyzed. On the other hand, dose-volume histogram, the mean and maximum doses of the OARs such as bladder, rectum and femoral heads for different collimator angles were determined from the plans. Results: There was no significance difference, based on the plan dose-volume evaluation criteria, found in the VMAT optimizations for all studied collimator angles. Higher CI and lower HI were found for the 45o collimator angle. In addition, the 15o collimator angle provided lower HI similar to the 45o collimator angle. The 75o and 90o collimator angle were found good for the rectum sparing, and the 75o and 30o collimator angle were found good for the right and left femur sparing, respectively. The PTV dose coverage for each plan was comparatively independent of the collimator angle. Conclusion: The dosimetric results in this study are useful to the planner to select different collimator angles to improve the PTV coverage and OAR sparing in prostate VMAT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-632: A Dosimetric Comparison of the 3D-CRT Planning of Chest Wall in Post-Mastectomy Breast Cancer Patients, with and Without Breast Board Setup

Purpose: Breast boards are used in breast radiation which increases normal lung and heart doses, when supraclavicular field is included. Therefore, in this study through dose volume histogram (DVHs), lung and heart doses comparison was done between two different setups i.e. with and without breast board, for the treatment of left chest wall and supraclavicular fossa in postmastectomy left breast cancer. Methods: In this study, CT-Simulation scans of ten breast cancer patients were done with and without breast board, at Shifa International Hospitals Islamabad, to investigate the differences between the two different setups of the irradiation of left chest wall in terms of lung and heart doses. For immobilization, support under the neck, shoulders and arms was used. Precise PLAN 2.15 treatment planning system (TPS) was used for 3D-CRT planning. The total prescribed dose for both the plans was 5000 cGy/25 fractions. The chest wall was treated with a pair of tangential photon fields and the upper supraclavicular nodal regions were treated with an anterior photon field. A mono-isocentric technique was used to match the tangential fields with the anterior field at the isocentre. The dose volume histogram was used to compare the doses of heart and ipsilateral lung. Results: Both the plans of each patient were generated and compared. DVH results showed that for the same PTV dose coverage, plans without breast board resulted in a reduction of lung and heart doses compared with the plans with breast board. There was significant reductions in V20, V<25 and mean doses for lung and V<9 and mean doses for heart. Conclusion: In comparison of both the plans, setup without breast board significantly reduced the dose-volume of the ipsilateral lung and heart in left chest wall patients. Waived registration request has been submitted.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-634: Analysis of Volume Based GYN HDR Brachytherapy Plans for Dose Calculation to Organs At Risk(OAR)

Purpose: We have analyzed the dose volume histogram of 140 CT based HDR brachytherapy plans and evaluated the dose received to OAR ; rectum, bladder and sigmoid colon based on recommendations from ICRU and Image guided brachytherapy working group for cervical cancer . Methods: Our treatment protocol consist of XRT to whole pelvis with 45 Gy at 1.8Gy/fraction followed by 30 Gy at 6 Gy per fraction by HDR brachytherapy in 2 weeks . The CT compatible tandem and ovoid applicators were used and stabilized with radio opaque packing material. The patient was stabilized using special re-locatable implant table and stirrups for reproducibility of the geometry during treatment. The CT scan images were taken at 3mm slice thickness and exported to the treatment planning computer. The OAR structures, bladder, rectum and sigmoid colon were outlined on the images along with the applicators. The prescription dose was targeted to A left and A right as defined in Manchester system and optimized on geometry . The dosimetry was compared on all plans using the parameter Ci.sec.cGy-1 . Using the Dose Volume Histogram (DVH) obtained from the plans the doses to rectum, sigmoid colon and bladder for ICRU defined points and 2cc volume were analyzed and reported. The following criteria were used for limiting the tolerance dose by volume (D2cc) were calculated. The rectum and sigmoid colon doses were limited to <75Gy. The bladder dose was limited to < 90Gy from both XRT and HDR brachytherapy. Results: The average total (XRT+HDRBT) BED values to prescription volume was 120 Gy. Dose 2cc to rectum was 70Gy +/− 17Gy, dose to 2cc bladder was 82+/−32 Gy. The average Ci.sec.cGy-1 calculated for the HDR plans was 6.99 +/− 0.5 Conclusion: The image based treatment planning enabled to evaluati volume based dose to critical structures for clinical interpretation.

60 APPLIED LIFE SCIENCES↗