Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATLAB”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-J-09: Radioactive Seed Localization for Breast Lumpectomy - Towards Optimization

Purpose: The aim of this work is to propose a method to optimize radioactive source localization (RSL) for non-palpable breast cancer surgery. RSL is commonly used as a guiding technique during surgery for excision of non-palpable tumors. A collimated hand-held detector is used to localize radioactive sources implanted in tumors. Incisions made by the surgeon are based on maximum observed detector counts, and tumors are subsequently resected based on an arbitrary estimate of the counts expected at the surgical margin boundary. This work focuses on building a framework to predict detector counts expected throughout the procedure to improve surgical margins. Methods: A gamma detection system called the Neoprobe GDS was used for this work. The probe consists of a cesium zinc telluride crystal and a collimator. For this work, an I-125 Best Medical model 2301 source was used. The source was placed in three different phantoms, a PMMA, a Breast (25%- glandular tissue/75%- adipose tissue) and a Breast (75-25) phantom with a backscatter thickness of 6 cm. Counts detected by the probe were recorded with varying amounts of phantom thicknesses placed on top of the source. A calibration curve was generated using MATLAB based on the counts recorded for the calibration dataset acquired with the PMMA phantom. Results: The observed detector counts data used as the validation set was accurately predicted to within ±3.2%, ±6.9%, ±8.4% for the PMMA, Breast (75-25), Breast (25–75) phantom respectively. The average difference between predicted and observed counts was −0.4%, 2.4%, 1.4% with a standard deviation of 1.2 %, 1.8%, 3.4% for the PMMA, Breast (75-25), Breast (25–75) phantom respectively. Conclusion: The results of this work provide a basis for characterization of a detector used for RSL. Counts were predicted to within ±9% for three different phantoms without the application of a density correction factor.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-99: Dose Accumulation and Evaluation in Lung SBRT Among All Phases of Respiration

Purpose: To calculate the total planning dose on lung tumors (GTV) by accumulating the dose received in all respiration phases. Methods: A patient 4D planning CT (phase-binned, from a Siemens Somatom CT) was used to locate the GTV of a lung tumor in all respiratory phases with Pinnacle (v9.10). GTV contours defined in all phases were projected to the reference phase, where the ITV was defined. Centroids were calculated for all the GTV projections. No deformation or rotation was taken into account. The only GTV contour as defined in the reference phase was voxelized to track each voxel individually. We accumulated the absorbed dose in different phases on each voxel. A 3DCRT and a VMAT plan were designed on the reference phase fulfilling the ITV dosimetric requirements, using the 10MV FFF photon model from an Elekta Versa linac. ITV-to-PTV margins were set to 5mm. In-house developed MATLAB code was used for tumor voxeling and dose accumulation, assuming that the dose distribution planned in the reference phase behaved as a “dose-cloud” during patient breathing. Results: We tested the method on a patient 4DCT set of images exhibiting limited tumor motion (<5mm). For the 3DCRT plan, D95 was calculated for the GTV with motion and for the ITV, showing an agreement of 0.04%. For the VMAT plan, we calculated the D95 for every phase as if the GTV in that phase had received the whole treatment. Differences in D95 for all phases are within 1%, and estimate the potential interplay effect during delivery. Conclusion: A method for dose accumulation and assessment was developed that can compare GTV motion with ITV dosage, and estimate the potential interplay effect for VMAT plans. Work in progress includes the incorporation of deformable image registration and 4D CBCT dose calculation for dose reconstruction and assessment during treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-102: Lower Esophagus Margin Implications Based On Rapid Computational Algorithm for SBRT

Purpose: To quantify inter-fraction esophagus-variation. Methods: Computed tomography and daily on-treatment 0.3-T MRI data sets for 7 patients were analyzed using a novel Matlab-based (Mathworks, Natick, MA) rapid computational method. Rigid registration was performed from the cricoid to the gastro-esophageal junction. CT and MR-based contours were compared at slice intervals of 3mm. Variation was quantified by “expansion,” defined as additional length in any radial direction from CT contour to MR contour. Expansion computations were performed with 360° of freedom in each axial slice. We partitioned expansions into left anterior, right anterior, right posterior, and left posterior quadrants (LA, RA, RP, and LP, respectively). Sample means were compared by analysis of variance (ANOVA) and Fisher’s Protected Least Significant Difference test. Results: Fifteen fractions and 1121 axial slices from 7 patients undergoing SBRT for primary lung cancer (3) and metastatic lung disease (4) were analyzed, generating 41,970 measurements. Mean LA, RA, RP, and LP expansions were 4.30±0.05 mm, 3.71±0.05mm, 3.17±0.07, and 3.98±0.06mm, respectively. 50.13% of all axial slices showed variation > 5 mm in one or more directions. Variation was greatest in lower esophagus with mean LA, RA, RP, and LP expansion (5.98±0.09 mm, 4.59±0.09 mm, 4.04±0.16 mm, and 5.41±0.16 mm, respectively). The difference was significant compared to mid and upper esophagus (p<.0001). The 95th percentiles of expansion for LA, RA, RP, LP were 13.36 mm, 9.97 mm, 11.29 mm, and 12.19 mm, respectively. Conclusion: Analysis of on-treatment MR imaging of the lower esophagus during thoracic SBRT suggests margin expansions of 13.36 mm LA, 9.97 mm RA, 11.29 mm RP, 12.19 mm LP would account for 95% of measurements. Our novel algorithm for rapid assessment of margin expansion for critical structures with 360° of freedom in each axial slice enables continuously adaptive patient-specific margins which may reduce overlap with high-dose regions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-174: A Series of Computational Human Phantoms in DICOM-RT Format for Normal Tissue Dose Reconstruction in Epidemiological Studies

Purpose: Epidemiological studies of second cancer risk in radiotherapy patients often require individualized dose estimates of normal tissues. Prior to 3D conformal radiation therapy planning, patient anatomy information was mostly limited to 2D radiological images or not even available. Generic patient CT images are often used in commercial radiotherapy treatment planning system (TPS) to reconstruct normal tissue doses. The objective of the current work was to develop a series of reference size computational human phantoms in DICOM-RT format for direct use in dose reconstruction in TPS. Methods: Contours of 93 organs and tissues were extracted from a series of pediatric and adult hybrid computational human phantoms (newborn, 1-, 5-, 10-, 15-year-old, and adult males and females) using Rhinoceros software. A MATLAB script was created to convert the contours into the DICOM-RT structure format. The simulated CT images with the resolution of 1×1×3 mm3 were also generated from the binary phantom format and coupled with the DICOM-structure files. Accurate volumes of the organs were drawn in the format using precise delineation of the contours in converted format. Due to complex geometry of organs, higher resolution (1×1×1 mm3) was found to be more efficient in the conversion of newborn and 1-year-old phantoms. Results: Contour sets were efficiently converted into DICOM-RT structures in relatively short time (about 30 minutes for each phantom). A good agreement was observed in the volumes between the original phantoms and the converted contours for large organs (NRMSD<1.0%) and small organs (NRMSD<7.7%). Conclusion: A comprehensive series of computational human phantoms in DICOM-RT format was created to support epidemiological studies of second cancer risks in radiotherapy patients. We confirmed the DICOM-RT phantoms were successfully imported into the TPS programs of major vendors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-180: A Reference Data Set for Testing Two Dimension Registration Algorithms

Purpose: To create and characterize a reference data set for testing image registration algorithms that transform portal image (PI) to digitally reconstructed radiograph (DRR). Methods: Anterior-posterior (AP) and Lateral (LAT) projection and DRR image pairs from nine cases representing four different anatomical sites (head and neck, thoracic, abdominal, and pelvis) were selected for this study. Five experts will perform manual registration by placing landmarks points (LMPs) on the DRR and finding their corresponding points on the PI using computer assisted manual point selection tool (CAMPST), a custom-made MATLAB software tool developed in house. The landmark selection process will be repeated on both the PI and the DRR in order to characterize inter- and -intra observer variations associated with the point selection process. Inter and an intra observer variation in LMPs was done using Bland-Altman (B&A) analysis and one-way analysis of variance. We set our limit such that the absolute value of the mean difference between the readings should not exceed 3mm. Later on in this project we will test different two dimension (2D) image registration algorithms and quantify the uncertainty associated with their registration. Results: Using one-way analysis of variance (ANOVA) there was no variations within the readers. When Bland-Altman analysis was used the variation within the readers was acceptable. The variation was higher in the PI compared to the DRR.ConclusionThe variation seen for the PI is because although the PI has a much better spatial resolution the poor resolution on the DRR makes it difficult to locate the actual corresponding anatomical feature on the PI. We hope this becomes more evident when all the readers complete the point selection. The reason for quantifying inter- and -intra observer variation tells us to what degree of accuracy a manual registration can be done. Research supported by William Beaumont Hospital Research Start Up Fund.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-209: Quantification of Image-Guidance Benefit in Image-Guided Radiotherapy of Cancers

Purpose: Image-guidance has been widely used in radiation oncology for accurate radiotherapy. The goal of this study is to quantify the benefit of image-guidance in image-guided radiotherapy (IGRT) of cancers. Methods: In this study, a new index termed image-guidance benefit (IGB), was proposed to quantify the benefit of image-guidance in cancer radiotherapy. It is calculated as a ratio of the square sum of dose differences between planning dose matrix and actual delivery dose matrix post image-guidance to the square sum of dose summation between the two matrixes, summing over all dose scoring voxels. Ranging from 0 to 1, larger IGB values indicate larger benefit out of image-guidance. With IRB approval, the DICOM RT files and 3D couch shifts applied during IGRT of 2219 patients were collected, based on which patient-specific IGB values were calculated with an in-house MATLAB code. Results: In this study, the mean IGB value was found to be 0.0398 (0.000583–0.999) with a positive correlation between IGB value and 3D couch shift vector at 0.0435 per cm (P<0.0001). With 2 mm shift as a threshold above which an image-guidance is deemed clinically necessary, the corresponding mean IGB value was 0.00457, much less than 0.0398 (P<0.001). However, the IGB values of 56 cases based on couch shifts were less than those based on 2 mm shift. Conclusion: The IGB values were patient-specific and site-dependent. Using 2 mm shifts as criterion for applying image-guidance, the applied image-guidance procedures were found clinically necessary and highly beneficial in 97.5% of cancer patients. However, image-guidance procedures were found over-used in about 2.5% of cancer patients in our current practices of IGRT, which should be avoided as they added no benefit in improving delivery accuracy while increasing cancer risk.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-20: Novel Catheter Lumen Recognition Algorithm for Rapid Digitization

Purpose: Manual catheter recognition remains a time-consuming aspect of high-dose-rate brachytherapy (HDR) treatment planning. In this work, a novel catheter lumen recognition algorithm was created for accurate and rapid digitization. Methods: MatLab v8.5 was used to create the catheter recognition algorithm. Initially, the algorithm searches the patient CT dataset using an intensity based k-means filter designed to locate catheters. Once the catheters have been located, seed points are manually selected to initialize digitization of each catheter. From each seed point, the algorithm searches locally in order to automatically digitize the remaining catheter. This digitization is accomplished by finding pixels with similar image curvature and divergence parameters compared to the seed pixel. Newly digitized pixels are treated as new seed positions, and hessian image analysis is used to direct the algorithm toward neighboring catheter pixels, and to make the algorithm insensitive to adjacent catheters that are unresolvable on CT, air pockets, and high Z artifacts. The algorithm was tested using 11 HDR treatment plans, including the Syed template, tandem and ovoid applicator, and multi-catheter lung brachytherapy. Digitization error was calculated by comparing manually determined catheter positions to those determined by the algorithm. Results: he digitization error was 0.23 mm ± 0.14 mm axially and 0.62 mm ± 0.13 mm longitudinally at the tip. The time of digitization, following initial seed placement was less than 1 second per catheter. The maximum total time required to digitize all tested applicators was 4 minutes (Syed template with 15 needles). Conclusion: This algorithm successfully digitizes HDR catheters for a variety of applicators with or without CT markers. The minimal axial error demonstrates the accuracy of the algorithm, and its insensitivity to image artifacts and challenging catheter positioning. Future work to automatically place initial seed positions would improve the algorithm speed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-46: The Effect of Inter-Seed Attenuation and Tissue Composition in Prostate 125I Brachytherapy Dose Calculations

Purpose: To investigate the difference of dose distributions with/without the effect of inter-seed attenuation and tissue compositions in prostate {sup 125}I brachytherapy dose calculations, using Monte Carlo simulations of Particle and Heavy Ion Transport code System (PHITS). Methods: The dose distributions in {sup 125}I prostate brachytherapy were calculated using PHITS for non-simultaneous and simultaneous alignments of STM1251 sources in water or prostate phantom for six patients. The PHITS input file was created from DICOM-RT file which includes source coordinates and structures for clinical target volume (CTV) and organs at risk (OARs) of urethra and rectum, using in-house Matlab software. Photon and electron cutoff energies were set to 1 keV and 100 MeV, respectively. The dose distributions were calculated with the kerma approximation and the voxel size of 1 × 1 × 1 mm{sup 3}. The number of incident photon was set to be the statistical uncertainty (1σ) of less than 1%. The effect of inter-seed attenuation and prostate tissue compositions was evaluated from dose volume histograms (DVHs) for each structure, by comparing to results of the AAPM TG-43 dose calculation (without the effect of inter-seed attenuation and prostate tissue compositions). Results: The dose reduction due to the inter-seed attenuation by source capsules was approximately 2% for CTV and OARs compared to those of TG-43. In additions, by considering prostate tissue composition, the D{sub 90} and V{sub 100} of CTV reduced by 6% and 1%, respectively. Conclusion: It needs to consider the dose reduction due to the inter-seed attenuation and tissue composition in prostate {sup 125}I brachytherapy dose calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-73: Experimental Determination of the Effective Point of Measurement in Electron Beams Using a Commercial Scintillation Detector

Purpose: The aim of this work was to determine experimentally the effective point of measurement (EPOM) in clinical electron beams for three cylindrical ionization chambers using a commercial scintillation detector as a reference detector. Methods: Percent depth dose (PDD) curves were measured using an Exradin W1 scintillation detector and were used as a representative PDD to water. Depth dose curves were measured with the Exradin A18, A1SL, and A28 ionization chambers. The raw ionization chamber curve data were corrected by the chamber fluence perturbation correction factor and restricted mass collisional stopping power ratio at each depth to obtain a percent depth dose curve to the gas volume (PDDGV) of the detector. Ratios of the W1 PDD to the ion chamber PDDGV were calculated for each measurement depth. The W1 PDD curve was shifted by small depth increments, Δz, until the ratio of the W1 PDD to the ion chamber PDDGV was depth-independent (optimal Δz). A MATLAB routine was developed to determine the optimal Δz value. Results: The optimal Δz shift was used as an estimate of the EPOM for each chamber. The average calculated EPOM shifts (expressed as a fraction of the chamber cavity radius) for the A18, A1SL, and A28 ionization chambers were 0.21 ± 0.04, 0.10 ± 0.05, and 0.22 ± 0.03, respectively. Conclusion: The experimentally determined EPOM values for the A18 and A1SL in this work agreed with the simulated values of Muir and Rogers (MedPhys 2014). The results also indicate that the Exradin W1 scintillator is water equivalent for electron energies of 6 MeV, 9 MeV, 12 MeV, and 16 MeV. In addition, we confirmed that the AAPM TG51 recommended EPOM shift of 0.5 times the cavity radius is not accurate for the A18 and A1SL chambers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-F-T-121: Abdominal Compression Effectively Reduces the Interplay Effect and Enables Pencil Beam Scanning Proton Therapy of Liver Tumors

Purpose: To study if abdominal compression can reduce breathing motion and mitigate interplay effect in pencil beam scanning proton therapy (PBSPT) treatment of liver tumors in order to better spare healthy liver volumes compared with photon therapy. Methods: Ten patients, six having large tumors initially treated with IMRT and four having small tumors treated with SBRT, were replanned for PBSPT. ITV and beam-specific PTVs based on 4D-CT were used to ensure target coverage in PBSPT. The use of an abdominal compression belt and volumetric repainting was investigated to mitigate the interplay effect between breathing motion and PBSPT dynamic delivery. An in-house Matlab script has been developed to simulate this interplay effect. The dose is computed on each phase individually by sorting all spots according to their simulated delivery timing. The final dose distribution is then obtained by accumulating all dose maps to a reference phase. Results: For equivalent target coverage PBSPT reduced average healthy liver dose by 9.5% of the prescription dose compared with IMRT/SBRT. Abdominal compression of 113.2±42.2 mmHg was effective for all 10 patients and reduced average motion by 2.25 mm. As a result, the average ITV volume decreased from 128.2% to 123.1% of CTV volume. Similarly, the average beam-specific PTV volume decreased from 193.2% to 183.3%. For 8 of the 10 patients, the average motion was reduced below 5 mm, and up to 3 repainting were sufficient to mitigate interplay. For the other two patients with larger residual motion, 4–5 repainting were needed. Conclusion: We recommend evaluation of the 4DCT motion histogram following simulation and the interplay effect following treatment planning in order to personalize the use of compression and volumetric repainting for each patient. Abdominal compression enables safe and more effective PBS treatment of liver tumors by reduction of motion and interplay effect. Kevin Souris is supported by IBA and Televie Grant from F.R.S.-FNRS. Liyong Lin is partially supported by Varian.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-138: Commissioning and Evaluating Dose Computation Models for a Dedicated Proton Line Scanning Beam Nozzle in Eclipse Treatment Planning System

Purpose: In this study, we present an effective method to derive low dose envelope of the proton in-air spot fluence at beam positions other than the isocenter to reduce amount of measurements required for planning commission. Also, we demonstrate commissioning and validation results of this method to the Eclipse treatment planning system (version 13.0.29) for a Sumitomo dedicated proton line scanning beam nozzle. Methods: The in-air spot profiles at five beam-axis positions (±200, ±100 and 0 mm) were obtained in trigger mode using a MP3 Water tank (PTW-Freiburg) and a pinpoint ionization chamber (model 31014, PTW-Freiburg). Low dose envelope (below 1% of the center dose) of the spot profile at isocenter was obtained by repeated point measurements to minimize dosimetry uncertainty. The double Gaussian (DG) model was used to fit and obtain optimal σ1, σ2 and their corresponding weightings through our in-house MATLAB (Mathworks) program. σ1, σ2 were assumed to expand linearly along the beam axis from a virtual source position calculated by back projecting fitted sigmas from the single Gaussian (SG) model. Absolute doses in water were validated using an Advanced Markus chamber at the depth of 2cm with Pristine Peak (BP) R90d ranging from 5–32 cm for 10×10 cm2 scanned fields. The field size factors were verified with square fields from 2 to 20 cm at 2cm and before BP depth. Results: The absolute dose outputs were found to be within ±3%. For field size factor, the agreement between calculated and measurement were within ±2% at 2cm and ±3% before BP, except for the field size below 2×2 cm2. Conclusion: The double Gaussian model was found to be sufficient for characterizing the Sumitomo dedicated proton line scanning nozzle. With our effective double Gaussian fitting method, we are able to save significant proton beam time with acceptable output accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-145: Modeling of a Proton Pencil Beam Partially Blocked by a Collimator

Purpose: To develop a mathematical model of the fluence distribution of a proton pencil beam that is partially blocked by a collimator. Methods: In one dimension, an exact analytical solution of the proton distribution behind a collimator edge was derived, taking the beam’s phase space, the edge position, and the scattering downstream of the collimator as input. Extension to two dimensions was achieved using a Clarkson-like Method: determine the radial distances from each calculation point to the collimator outline, calculate the beam’s phase space in the direction of each radial line, apply the one-dimensional model along each radial line, and sum the contributions. A Matlab implementation of the full model produced high-resolution distributions with relatively short calculation times (of about 2 ms per calculation point). The model was validated by comparison to Monte Carlo simulations for select pencil beam phase space distributions and collimator shapes. Results: Agreement between model and simulation confirm the exact mathematical solution of the one-dimensional model. In two dimensions the agreement is reasonable, but not exact, due to the approximations introduced in the Clarkson-like method. Conclusion: A model has been developed and validated that gives the lateral fluence distribution of a pencil beam beyond a collimator. Notably, this model does not include free parameters and the final fluence is completely determined by the phase space parameters of the beam, already existing in the treatment planning system, as well as both axial and lateral position of the aperture edge. As such, it can be incorporated in a full pencil-beam dose calculation algorithm allowing for the modeling of field-specific collimators in proton pencil beam scanning and, potentially, the optimization of collimator shape and/or pencil phase space in addition to the customary spot weight optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SU-F-T-154: An Evaluation and Quantification of Secondary Neutron Radiation Dose Due to Double Scatter and Pencil Beam Scanning Proton Therapy

Proton radiation therapy can deliver high radiation doses to tumors while sparing normal tissue. However, protons yield secondary neutron and gamma radiation that is difficult to detect, small in comparison to the prescribed dose, and not accounted for in most treatment planning systems. The risk for secondary malignancies after proton therapy may be dependent on the quality of this dose. Consequently, there is interest in characterizing the secondary radiation. Previously, we used the dual ionization chamber method to measure the separate absorbed dose from gamma-rays and neutrons secondary to the proton beam1, relying on characterization of ionization chamber response in the unknown neutron spectrum from Monte Carlo simulation. We developed a procedure to use Shieldwerx activation foils, with neutron activation energies ranging from 0.025 eV to 13.5 MeV, to measure the neutron energy spectrum from double scattering (DS) and pencil beam scanning (PBS) protons outside of the treatment volume in a water tank. The activated foils are transferred to a NaI well chamber for gamma-ray spectroscopy and activity measurement. Since PBS treats in layers, the switching time between layers is used to correct for the decay of the activated foils and the relative dose per layer is assumed to be proportional to the neutron fluence per layer. MATLAB code was developed to incorporate the layer delivery and switching time into a calculation of foil activity, which is then used to determine the neutron energy fluence from tabulated foil activation energy thresholds.1. Diffenderfer et. al., Med. Phys., 38(11) 2011.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-156: Monte Carlo Simulation Using TOPAS for Synchrotron Based Proton Discrete Spot Scanning System

Purpose: This study provides an overview of the design and commissioning of the Monte Carlo (MC) model of the spot-scanning proton therapy nozzle and its implementation for the patient plan simulation. Methods: The Hitachi PROBEAT V scanning nozzle was simulated based on vendor specifications using the TOPAS extension of Geant4 code. FLUKA MC simulation was also utilized to provide supporting data for the main simulation. Validation of the MC model was performed using vendor provided data and measurements collected during acceptance/commissioning of the proton therapy machine. Actual patient plans using CT based treatment geometry were simulated and compared to the dose distributions produced by the treatment planning system (Varian Eclipse 13.6), and patient quality assurance measurements. In-house MATLAB scripts are used for converting DICOM data into TOPAS input files. Results: Comparison analysis of integrated depth doses (IDDs), therapeutic ranges (R90), and spot shape/sizes at different distances from the isocenter, indicate good agreement between MC and measurements. R90 agreement is within 0.15 mm across all energy tunes. IDDs and spot shapes/sizes differences are within statistical error of simulation (less than 1.5%). The MC simulated data, validated with physical measurements, were used for the commissioning of the treatment planning system. Patient geometry simulations were conducted based on the Eclipse produced DICOM plans. Conclusion: The treatment nozzle and standard option beam model were implemented in the TOPAS framework to simulate a highly conformal discrete spot-scanning proton beam system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-187: Quantifying Normal Tissue Sparing with 4D Robust Optimization of Intensity Modulated Proton Therapy

Purpose: To report an approach to quantify the normal tissue sparing for 4D robustly-optimized versus PTV-optimized IMPT plans. Methods: We generated two sets of 90 DVHs from a patient’s 10-phase 4D CT set; one by conventional PTV-based optimization done in the Eclipse treatment planning system, and the other by an in-house robust optimization algorithm. The 90 DVHs were created for the following scenarios in each of the ten phases of the 4DCT: ± 5mm shift along x, y, z; ± 3.5% range uncertainty and a nominal scenario. A Matlab function written by Gay and Niemierko was modified to calculate EUD for each DVH for the following structures: esophagus, heart, ipsilateral lung and spinal cord. An F-test determined whether or not the variances of each structure’s DVHs were statistically different. Then a t-test determined if the average EUDs for each optimization algorithm were statistically significantly different. Results: T-test results showed each structure had a statistically significant difference in average EUD when comparing robust optimization versus PTV-based optimization. Under robust optimization all structures except the spinal cord received lower EUDs than PTV-based optimization. Using robust optimization the average EUDs decreased 1.45% for the esophagus, 1.54% for the heart and 5.45% for the ipsilateral lung. The average EUD to the spinal cord increased 24.86% but was still well below tolerance. Conclusion: This work has helped quantify a qualitative relationship noted earlier in our work: that robust optimization leads to plans with greater normal tissue sparing compared to PTV-based optimization. Except in the case of the spinal cord all structures received a lower EUD under robust optimization and these results are statistically significant. While the average EUD to the spinal cord increased to 25.06 Gy under robust optimization it is still well under the TD50 value of 66.5 Gy from Emami et al. Supported in part by the NCI U19 CA021239.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-210: The Variable Virtual Source-To-Axis Distance Effect On A Compact Proton Pencil Beam Scanning System

Purpose: We investigate the spot characteristic and dose profiles properties from a compact gantry proton therapy system. This compact design features a dedicated pencil beam scanning nozzle with the scanning magnet located upstream of the final 60 degree bending magnet. Due to the unique beam line design, uncertainty has been raised in the virtual source-to-axis distance (SAD). We investigate its potential clinical impact through measurements and simulation. Methods: A scintillator camera based detector was used to measure spot characteristics and position accuracy. An ion chamber array device was used to measure planar dose profile. Dose profile in-air simulation was performed using in-house built MATLAB program based on additional spot parameters directly from measurements. Spot characteristics such as position and in-air sigma values were used to general simulated 2D elliptical Gaussian spots. The virtual SAD distance changes in the longitudinal direction were also simulated. Planar dose profiles were generated by summation of simulated spots at the isocenter, 15 cm above the isocenter, and 15 cm below the isocenter for evaluation of potential clinical dosimetric impact. Results: We found that the virtual SAD varies depending on the spot location on the longitudinal axis. Measurements have shown that the variable SAD changes from 7 to 12 meters from one end to the other end of the treatment field in the longitudinal direction. The simulation shows that the planer dose profiles differences between the fixed SAD and variable SAD are within 3% from the isocenter profile and the lateral penumbras are within 1 mm difference. Conclusion: Our measurements and simulations show that there are minimum effects on the spot characteristics and dose profiles for this up-stream scanning compact system proton system. Further treatment planning study is needed with the variable virtual SAD accounted for in the planning system to show minimum dosimetric impact.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-229: A Novel Method for EPID-Based In-Vivo Exit Dose Verification for Intensity Modulated Radiotherapy

Purpose: Electronic portal imaging device (EPID) can be used to acquire a two-dimensional exit dose distribution during treatment delivery, thus allowing the in-vivo verification of the dose delivery through a comparison of measured portal images to predicted portal dose images (PDI). The aim of this study was to present a novel method to easily and accurately predict PDI, and to establish an EPID-based in-vivo dose verification method during IMRT treatments. Methods: We developed a model to determine the predicted portal dose at the same plane of the EPID detector location. The Varian EPID (aS1000) positions at 150cm source-to-detector-distance (SDD), and can be used to acquire in-vivo exit dose using Portal Dosimetry (PD) function. Our model was generated to make an equivalent water thickness represent the buildup plate of EPID. The exit dose at extend SDD plane with patient CT data in the beam can be calculated as the predicted PDI in the treatment planning system (TPS). After that, the PDI was converted to the fluence at SDD of 150cm using the inverse square law coded in MATLAB. Five head-and-neck and prostate IMRT patient plans contain 32 fields were investigated to evaluate the feasibility of this new method. The measured EPID image was compared with PDI using the gamma analysis. Results: The average results for cumulative dose comparison were 81.9% and 91.6% for 3%, 3mm and 4%, 4mm gamma criteria, respectively. Results indicate that the patient transit dosimetry predicted algorithm compares well with EPID measured PD doses for test situations. Conclusion: Our new method can be used as an easy and feasible tool for online EPID-based in-vivo dose delivery verification for IMRT treatments. It can be implemented for fast detecting those obvious treatment delivery errors for individual field and patient quality assurance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-259: GPR Tables for the Estimation of Mid-Plane Dose Using EPID

Purpose: To develop a simple method for estimating the mid-plane dose (MPD) of a patient using Electronic Portal imaging Device (EPID). Methods: A Varian TrueBeam with aSi100 EPID was used in this study. The EPID images were acquired for a 30 cm × 30 cm homogeneous slab phantom and a 30 cm diameter 20 cm thick cylindrical phantom in the continuous dosimetry mode. The acquired EPID images in XIM format were imported into in-house MATLAB program for the data analysis. First, the dosimetric characteristics of EPID were studied for dose-response linearity, dose-rate dependence, and field size dependence. Next, the average pixels values of the EPID images were correlated with the MPD measured by an ionisation chamber for various thicknesses of the slab phantom (8 cm – 30 cm) and for various square field sizes (3×3 cm{sup 2} – 25×25 cm{sup 2} at the isocenter). Look-up tables called as GPR tables were then generated for both SSD and SAD setup by taking the ratio of MPD measured by the ionisation chamber and the corresponding EPID pixel values. The accuracy of the GPR tables was evaluated by varying the field size, phantom thickness, and wedge angles with the slab and cylindrical phantoms. Results: The dose response of EPID was linear from 20 MU to 300 MU. The EPID response for different dose rates from 40 MU/min to 600 MU/min was within ±1%. The difference in the doses from the GPR tables and the doses measured by the ionization chambers were within 2% for slab phantoms, and 3% for the cylindrical phantom for various field sizes, phantom thickness, and wedge angles. Conclusion: GPR tables are a ready reckoner for in-vivo dosimetry and it can be used to quickly estimate the MPD value from the EPID images with an accuracy of ±3% for common clinical treatment. project work funded by Union for International cancer control(UICC) under ICRETT fellowship.

60 APPLIED LIFE SCIENCES↗