Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11

Artificial Neural Network as the FPGA Trigger in the Cyclone V based Front-End for a Detection of Neutrino-Origin Showers

Neutrinos play a fundamental role in the understanding of the origin of ultra-high-energy cosmic rays. They interact through charged and neutral currents in the atmosphere generating extensive air showers. However, their a very low rate of events potentially generated by neutrinos is a significant challenge for a detection technique and requires both sophisticated algorithms and high-resolution hardware. A trigger based on a artificial neural network was implemented into the Cyclone{sup R} V E FPGA 5CEFA9F31I7 - the heart of the prototype Front-End boards developed for tests of new algorithms in the Pierre Auger surface detectors. Showers for muon and tau neutrino initiating particles on various altitudes, angles and energies were simulated in CORSICA and Offline platforms giving pattern of ADC traces in Auger water Cherenkov detectors. The 3-layer 12-8-1 neural network was taught in MATLAB by simulated ADC traces according the Levenberg-Marquardt algorithm. Results show that a probability of a ADC traces generation is very low due to a small neutrino cross-section. Nevertheless, ADC traces, if occur, for 1-10 EeV showers are relatively short and can be analyzed by 16-point input algorithm. We optimized the coefficients from MATLAB to get a maximal range of potentially registered events and for fixed-point FPGA processing to minimize calculation errors. New sophisticated triggers implemented in Cyclone{sup R} V E FPGAs with large amount of DSP blocks, embedded memory running with 120 - 160 MHz sampling may support a discovery of neutrino events in the Pierre Auger Observatory. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Blockchain-Enabled Cyber-Secure Microgrid Control Using Consensus Algorithm

This paper proposes a distributed control method based on a consensus algorithm for distributed energy resources (DERs) using blockchain as a secure communication medium for cyber resilience. Each DER communicates with a local blockchain server that is maintained by the distributed ledger technology to allow for securely sharing local measurements among neighbouring assets to achieve the global control objectives, i.e., voltage and frequency regulation as well as accurate power sharing among the DERs, including collective grid-forming capability. To prove that the distributed control can retain system stability under a blockchain-induced variable communication delay, Lyapunov function-based stability analysis is carried out. This paper demonstrates the concept on an 11-bus test case developed in MATLAB Simulink, which has been modified from the IEEE 9-bus test case, to study dynamic operations of the five inverter-based DERs working with a blockchain-induced variable delay. The results validate the superior performance of the proposed control method, compared to heavily compromised operations of the other test cases without the security measure, affected by heavy communication delays and communication interruptions.

blockchains↗

Numerical Simulations of Cathodoluminescence Measurements in Thin-Film Solar Cells

Cathodoluminescence (CL) measurements can be applied to assess grain-boundary (GB) and grain-interior (GI) recombination in thin-film solar cell materials and made quantitative if we can develop CL models that account for material and measurement complexities. Recently, we developed a three-dimensional numerical CL model, based in MATLAB, that simulates the GI CL intensity as a function of four parameters: grain size, GI lifetime, and GB and surface recombination velocities. The model assumes that GB electrostatic potentials are screened by the high excesscarrier densities used in the CL measurement such that transport is governed by ambipolar diffusion. Here, we develop models to address directly GB potentials and their effects on these measurements. First, we transfer the MATLAB-based model to COMSOL software, and then introduce shallow donors to the GBs to produce potentials. We also develop a two-dimensional model in COMSOL to simulate CL GB contrast with GB potentials. Simulations indicate that GB potentials can increase or decrease CL intensities relative to the zero-potential case. However, the high electron-beam currents typically applied in CL measurements minimize the impact of GB potentials.

14 SOLAR ENERGY↗

SU-D-9A-03: STAMP: Simulator for Texture Analysis in MRI/PET

Purpose: To develop a convenient simulation platform to facilitate PET/MR image analysis with the prospect of gaining a better understanding of the influence of acquisition parameters on PET/MRI textural features. The simulation platform is demonstrated by showing textural variations of a representative case study using different image acquisition parameters. Methods: The simulation platform is composed of MRI simulators JEMRIS and SIMRI to achieve simulations of customized MRI sequences on sample tumor models. The PET simulator GATE is used to get 2D and 3D Monte Carlo acquisitions of voxelized PET sources using a phantom geometry and a customized scanner architecture. The platform incorporates a series of graphical user interfaces written in Matlab. Two GUIs are used to facilitate communication with the simulation executables installed on a computer cluster. A third GUI is used to collect and display the clinical and simulated images, as well as fused PET/MRI images, and perform computation of textural features.To illustrate the capabilities of this platform, one FDG-PET and T1-weighted (T1w) digitized tumor models were generated from clinical images of a soft-tissue sarcoma patient. Numerically simulated MR images were produced using 3 different echo times (TE) and 5 different repetition times (TR). PET 2D images were simulated using an OSEM algorithm with 1 to 32 iterations and a post-reconstruction Gaussian filter of 0, 2, 4 or 6 mm width. Results: STAMP was successfully used to produce numerically simulated FDG-PET and MRI images, and to calculate their corresponding textures. Three typical textures (GLCM-Contrast, GLSZM-ZSV and NGTDM-Coarseness) were found to vary by a range of 45% on average compared to reference scanning conditions in the case of FDG-PET, and by a range of 40% in the case of T1w MRI. Conclusion: We have successfully developed a Matlab-based simulation platform to facilitate PET/MRI texture image analysis for outcome prediction.

60 APPLIED LIFE SCIENCES↗

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↗

Sci—Thur AM: YIS - 11: Estimation of Bladder-Wall Cumulative Dose in Multi-Fraction Image-Based Gynaecological Brachytherapy Using Deformable Point Set Registration

Dosimetric parameters based on dose-volume histograms (DVH) of contoured structures are routinely used to evaluate dose delivered to target structures and organs at risk. However, the DVH provides no information on the spatial distribution of the dose in situations of repeated fractions with changes in organ shape or size. The aim of this research was to develop methods to more accurately determine geometrically localized, cumulative dose to the bladder wall in intracavitary brachytherapy for cervical cancer. The CT scans and treatment plans of 20 cervical cancer patients were used. Each patient was treated with five high-dose-rate (HDR) brachytherapy fractions of 600cGy prescribed dose. The bladder inner and outer surfaces were delineated using MIM Maestro software (MIM Software Inc.) and were imported into MATLAB (MathWorks) as 3-dimensional point clouds constituting the “bladder wall”. A point-set registration toolbox for MATLAB, Coherent Point Drift (CPD), was used to non-rigidly transform the bladder-wall points from four of the fractions to the coordinate system of the remaining (reference) fraction, which was chosen to be the emptiest bladder for each patient. The doses were accumulated on the reference fraction and new cumulative dosimetric parameters were calculated. The LENT-SOMA toxicity scores of these patients were studied against the cumulative dose parameters. Based on this study, there was no significant correlation between the toxicity scores and the determined cumulative dose parameters.

60 APPLIED LIFE SCIENCES↗

IBEX: An open infrastructure software platform to facilitate collaborative work in radiomics

Purpose: Radiomics, which is the high-throughput extraction and analysis of quantitative image features, has been shown to have considerable potential to quantify the tumor phenotype. However, at present, a lack of software infrastructure has impeded the development of radiomics and its applications. Therefore, the authors developed the imaging biomarker explorer (IBEX), an open infrastructure software platform that flexibly supports common radiomics workflow tasks such as multimodality image data import and review, development of feature extraction algorithms, model validation, and consistent data sharing among multiple institutions. Methods: The IBEX software package was developed using the MATLAB and C/C++ programming languages. The software architecture deploys the modern model-view-controller, unit testing, and function handle programming concepts to isolate each quantitative imaging analysis task, to validate if their relevant data and algorithms are fit for use, and to plug in new modules. On one hand, IBEX is self-contained and ready to use: it has implemented common data importers, common image filters, and common feature extraction algorithms. On the other hand, IBEX provides an integrated development environment on top of MATLAB and C/C++, so users are not limited to its built-in functions. In the IBEX developer studio, users can plug in, debug, and test new algorithms, extending IBEX’s functionality. IBEX also supports quality assurance for data and feature algorithms: image data, regions of interest, and feature algorithm-related data can be reviewed, validated, and/or modified. More importantly, two key elements in collaborative workflows, the consistency of data sharing and the reproducibility of calculation result, are embedded in the IBEX workflow: image data, feature algorithms, and model validation including newly developed ones from different users can be easily and consistently shared so that results can be more easily reproduced between institutions. Results: Researchers with a variety of technical skill levels, including radiation oncologists, physicists, and computer scientists, have found the IBEX software to be intuitive, powerful, and easy to use. IBEX can be run at any computer with the windows operating system and 1GB RAM. The authors fully validated the implementation of all importers, preprocessing algorithms, and feature extraction algorithms. Windows version 1.0 beta of stand-alone IBEX and IBEX’s source code can be downloaded. Conclusions: The authors successfully implemented IBEX, an open infrastructure software platform that streamlines common radiomics workflow tasks. Its transparency, flexibility, and portability can greatly accelerate the pace of radiomics research and pave the way toward successful clinical translation.

60 APPLIED LIFE SCIENCES↗

Novel spirometry based on optical surface imaging

Purpose: To evaluate the feasibility of using optical surface imaging (OSI) to measure the dynamic tidal volume (TV) of the human torso during free breathing. Methods: We performed experiments to measure volume or volume change in geometric and deformable phantoms as well as human subjects using OSI. To assess the accuracy of OSI in volume determination, we performed experiments using five geometric phantoms and two deformable body phantoms and compared the values with those derived from geometric calculations and computed tomography (CT) measurements, respectively. To apply this technique to human subjects, an institutional review board protocol was established and three healthy volunteers were studied. In the human experiment, a high-speed image capture mode of OSI was applied to acquire torso images at 4–5 frames per second, which was synchronized with conventional spirometric measurements at 5 Hz. An in-house MATLAB program was developed to interactively define the volume of interest (VOI), separate the thorax and abdomen, and automatically calculate the thoracic and abdominal volumes within the VOIs. The torso volume change (TV C = ΔV{sub torso} = ΔV{sub thorax} + ΔV{sub abdomen}) was automatically calculated using full-exhalation phase as the reference. The volumetric breathing pattern (BP{sub v} = ΔV{sub thorax}/ΔV{sub torso}) quantifying thoracic and abdominal volume variations was also calculated. Under quiet breathing, TVC should equal the tidal volume measured concurrently by a spirometer with a conversion factor (1.08) accounting for internal and external differences of temperature and moisture. Another MATLAB program was implemented to control the conventional spirometer that was used as the standard. Results: The volumes measured from the OSI imaging of geometric phantoms agreed with the calculated volumes with a discrepancy of 0.0% ± 1.6% (range −1.9% to 2.5%). In measurements from the deformable torso/thorax phantoms, the volume differences measured using OSI imaging and CT imaging were 1.2% ± 2.1% (range −0.5% to 3.6%), with a linear regression fitting (slope = 1.02 and R{sup 2} = 0.999). In volunteers, the relative error in OSI tidal volume measurement was −2.2% ± 4.9% (range −9.2% to 4.8%) and a correlation of r = 0.98 was found with spirometric measurement. The breathing pattern values of the three volunteers were substantially different from each other (BP{sub v} = 0.15, 0.45, and 0.32). Conclusions: This study demonstrates the feasibility of using OSI to measure breathing tidal volumes and breathing patterns with adequate accuracy. This is the first time that dynamic breathing tidal volume as well as breathing patterns is measured using optical surface imaging. The OSI-observed movement of the entire torso could serve as a new respiratory surrogate in the treatment room during radiation therapy.

60 APPLIED LIFE SCIENCES↗

SU-D-207-07: Implementation of Full/half Bowtie Filter Model in a Commercial Treatment Planning System for Kilovoltage X-Ray Imaging Dose Estimation

Purpose: To implement full/half bowtie filter models in a commercial treatment planning system (TPS) to calculate kilovoltage (kV) x-ray imaging dose of Varian On-Board Imager (OBI) cone beam CT (CBCT) system. Methods: Full/half bowtie filters of Varian OBI were created as compensator models in Pinnacle TPS (version 9.6) using Matlab software (version 2011a). The profiles of both bowtie filters were acquired from the manufacturer, imported into the Matlab system and hard coded in binary file format. A Pinnacle script was written to import each bowtie filter data into a Pinnacle treatment plan as a compensator. A kV x-ray beam model without including the compensator model was commissioned per each bowtie filter setting based on percent depth dose and lateral profile data acquired from Monte Carlo simulations. To validate the bowtie filter models, a rectangular water phantom was generated in the planning system and an anterior/posterior beam with each bowtie filter was created. Using the Pinnacle script, each bowtie filter compensator was added to the treatment plan. Lateral profile at the depth of 3cm and percent depth dose were measured using an ion chamber and compared with the data extracted from the treatment plans. Results: The kV x-ray beams for both full and half bowtie filter have been modeled in a commercial TPS. The difference of lateral and depth dose profiles between dose calculations and ion chamber measurements were within 6%. Conclusion: Both full/half bowtie filter models provide reasonable results in kV x-ray dose calculations in the water phantom. This study demonstrates the possibility of using a model-based treatment planning system to calculate the kV imaging dose for both full and half bowtie filter modes. Further study is to be performed to evaluate the models in clinical situations.

60 APPLIED LIFE SCIENCES↗

SU-E-T-651: Quantification of Dosimetric Accuracy of Respiratory Gated Stereotactic Body Radiation Therapy

Purpose: To quantify the dosimetric accuracy of respiratory gated stereotactic body radiation therapy delivery using dynamic thorax phantom. Methods: Three patients with mobile target (2 lung, 1liver) were chosen. Retrospective 4DCT image sets were acquired for using Varian RPM system. An in-house MATLAB program was designed for MIP, MinIP and AvgIP generation. ITV was contoured on MIP image set for lung patients and on MinIP for liver patient. Dynamic IMRT plans were generated on selected phase bin image set in Eclipse (v10.0) planning system. CIRS dynamic thorax phantom was used to perform the dosimetric quality assurance. Patient breathing pattern file from RPM system was converted to phantom compatible file by an in-house MATLAB program. This respiratory pattern fed to the CIRS dynamic thorax phantom. 4DCT image set was acquired for this phantom using patient breathing pattern. Verification plans were generated using patient gating window and delivered on the phantom. Measurements were carried out using with ion chamber and EBT2 film. Exposed films were analyzed and evaluated in FilmQA software. Results: The stability of gated output in comparison with un-gated output was within 0.5%. The Ion chamber measured and TPS calculated dose compared for all the patients. The difference observed was 0.45%, −0.52% and −0.54 for Patient 1, Patient2 and Patient 3 respectively.Gamma value evaluated from EBT film shows pass rates from 92.41% to 99.93% for 3% dose difference and 3mm distance to agreement criteria. Conclusion: Dosimetric accuracy of respiratory gated SBRT delivery for lung and liver was dosimetrically acceptable. The Ion chamber measured dose was within 0.203±0.5659% of the expected dose. Gamma pass rates were within 96.63±3.84% of the expected dose.

61 RADIATION PROTECTION AND DOSIMETRY↗

TU-AB-303-11: Predict Parotids Deformation Applying SIS Epidemiological Model in H&N Adaptive RT

Purpose: The aim is to investigate the use of epidemiological models to predict morphological variations in patients undergoing radiation therapy (RT). The susceptible-infected-susceptible (SIS) deterministic model was applied to simulate warping within a focused region of interest (ROI). Hypothesis is to consider each voxel like a single subject of the whole sample and to treat displacement vector fields like an infection. Methods: Using Raystation hybrid deformation algorithms and automatic re-contouring based on mesh grid, we post-processed 360 MVCT images of 12 H&N patients treated with Tomotherapy. Study focused on parotid glands, identified by literature and previous analysis, as ROI more susceptible to warping in H&N region. Susceptible (S) and infectious (I) cases were identified in voxels with inter-fraction movement respectively under and over a set threshold. IronPython scripting allowed to export positions and displacement data of surface voxels for every fraction. A MATLAB homemade toolbox was developed to model the SIS. Results: SIS model was validated simulating organ motion on QUASAR phantom. Applying model in patients, within a [0–1cm] range, a single voxel movement of 0.4cm was selected as displacement threshold. SIS indexes were evaluated by MATLAB simulations. Dynamic time warping algorithm was used to assess matching between model and parotids behavior days of treatments. The best fit of the model was obtained with contact rate of 7.89±0.94 and recovery rate of 2.36±0.21. Conclusion: SIS model can follow daily structures evolutions, making possible to compare warping conditions and highlighting challenges due to abnormal variation and set-up errors. By epidemiology approach, organ motion could be assessed and predicted not in terms of average of the whole ROI, but in a voxel-by-voxel deterministic trend. Identifying anatomical region subjected to variations, would be possible to focus clinic controls within a cohort of pre-selected patients eligible for adaptive RT. The research is partially co-funded by the Italian Research Grant: Dose warping methods for IGRT and Adaptive RT: dose accumulation based on organ motion and anatomical variations of the patients during radiation therapy treatments,MoH (GR-2010-2318757) and Tecnologie Avanzate S.r.l.(Italy)

60 APPLIED LIFE SCIENCES↗

Analytic reconstruction algorithms for triple-source CT with horizontal data truncation

Purpose: This paper explores a triple-source imaging method with horizontal data truncation to enlarge the field of view (FOV) for big objects. Methods: The study is conducted by using theoretical analysis, mathematical deduction, and numerical simulations. The proposed algorithms are implemented in c + + and MATLAB. While the basic platform is constructed in MATLAB, the computationally intensive segments are coded in c + +, which are linked via a MEX interface. Results: A triple-source circular scanning configuration with horizontal data truncation is developed, where three pairs of x-ray sources and detectors are unevenly distributed on the same circle to cover the whole imaging object. For this triple-source configuration, a fan-beam filtered backprojection-type algorithm is derived for truncated full-scan projections without data rebinning. The algorithm is also extended for horizontally truncated half-scan projections and cone-beam projections in a Feldkamp-type framework. Using their method, the FOV is enlarged twofold to threefold to scan bigger objects with high speed and quality. The numerical simulation results confirm the correctness and effectiveness of the developed algorithms. Conclusions: The triple-source scanning configuration with horizontal data truncation cannot only keep most of the advantages of a traditional multisource system but also cover a larger FOV for big imaging objects. In addition, because the filtering is shift-invariant, the proposed algorithms are very fast and easily parallelized on graphic processing units.

60 APPLIED LIFE SCIENCES↗

SU-F-R-18: Updates to the Computational Environment for Radiological Research for Image Analysis

Purpose: To present new tools in CERR for Texture Analysis and Visualization. Method: (1) Quantitative Image Analysis: We added the ability to compute Haralick texture features based on local neighbourhood. The Texture features depend on many parameters used in their derivation. For example: (a) directionality, (b) quantization of image, (c) patch-size for the neighborhood, (d) handling of the edge voxels within the region of interest, (e) Averaging co-occurance matrix vs texture features for different directions etc. A graphical user interface was built to set these parameters and then visualize their impact on the resulting texture maps. The entire functionality was written in Matlab. Array indexing was used to speed up the texture calculation. The computation speed is very competitive with the ITK library. Moreover, our implementation works with multiple CPUs and the computation time can be further reduced by using multiple processor threads. In order to reduce the Haralick texture maps into scalar features, we propose the use of Texture Volume Histograms. This lets users make use of the entire distribution of texture values within the region of interest rather than using just the mean and the standard deviations. (2) Qualitative/Visualization tools: The derived texture maps are stored as a new scan (derived) within CERR’s planC data structure. A display that compares various scans was built to show the raw image and the derived texture maps side-by-side. These images are positionally linked and can be navigated together. CERR’s graphics handling was updated and sped-up to be compatible with the newer Matlab versions. As a result, the users can use (a) different window levels and colormaps for different viewports, (b) click-and-drag or use mouse scroll-wheel to navigate slices. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics and Outcomes modeling.

60 APPLIED LIFE SCIENCES↗

SU-F-T-167: Measurement of Proton Beam Spot Size Using a 2D Ion Chamber Array for Daily QA

Purpose: To develop an algorithm which can calculate the Full Width Half Maximum (FWHM) of a Proton Pencil Beam from a 2D dimensional ion chamber array (IBA Matrixx) with limited spatial resolution ( 7.6 mm inter chamber distance). The algorithm would allow beam FWHM measurements to be taken during daily QA without an appreciable time increase. Methods: Combinations of 147 MeV single spot beams were delivered onto an IBA Matrixx and concurrently on EBT3 films for a standard. Data were collected around the Bragg Peak region and evaluated by a custom MATLAB script based on our algorithm using a least squared analysis. A set of artificial data, modified with random noise, was also processed to test for robustness. Results: The Matlab script processed Matixx data shows acceptable agreement (within 5%) with film measurements with no single measurement differing by more than 1.8 mm. In cases where the spots show some degree of asymmetry, the algorithm is able to resolve the differences. The algorithm was able to process artificial data with noise up to 15% of the maximum value. Time assays of each measurement took less than 3 minutes to perform, indicating that such measurements may be efficiently added to daily QA treatment. Conclusion: The developed algorithm can be implemented in daily QA program for Proton Pencil Beam scanning beams (PBS) with Matrixx to extract spot size and position information. The developed algorithm may be extended to small field sizes in photon clinic.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-F-T-458: Tracking Trends of TG-142 Parameters Via Analysis of Data Recorded by 2D Chamber Array

Purpose: With increasing QA demands of medical physicists in clinical radiation oncology, the need for an effective method of tracking clinical data has become paramount. A tool was produced which scans through data automatically recorded by a 2D chamber array and extracts relevant information recommended by TG-142. Using this extracted information a timely and comprehensive analysis of QA parameters can be easily performed enabling efficient monthly checks on multiple linear accelerators simultaneously. Methods: A PTW STARCHECK chamber array was used to record several months of beam outputs from two Varian 2100 series linear accelerators and a Varian NovalisTx−. In conjunction with the chamber array, a beam quality phantom was used to simultaneously to determine beam quality. A minimalist GUI was created in MatLab that allows a user to set the file path of the data for each modality to be analyzed. These file paths are recorded to a MatLab structure and then subsequently accessed by a script written in Python (version 3.5.1) which then extracts values required to perform monthly checks as outlined by recommendations from TG-142. The script incorporates calculations to determine if the values recorded by the chamber array fall within an acceptable threshold. Results: Values obtained by the script are written to a spreadsheet where results can be easily viewed and annotated with a “pass” or “fail” and saved for further analysis. In addition to creating a new scheme for reviewing monthly checks, this application allows for able to succinctly store data for follow up analysis. Conclusion: By utilizing this tool, parameters recommended by TG-142 for multiple linear accelerators can be rapidly obtained and analyzed which can be used for evaluation of monthly checks.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-T-531: Determination of Site-Specific Dynamic-Jaw Versus Static-Jaw RapidArc Delivery

Purpose: Dynamic-jaw tracking maximizes the area blocked by both jaw and MLC in RapidArc. We developed a method to quantify jaw tracking. Methods: An Eclipse Scripting API (ESAPI) was used to export beam parameters for each arc’s control points. The specific beam parameters extracted were: gantry angle, control point number, meterset, x-jaw positions, y-jaw positions, MLC bank-number, MLC leaf-number, and MLC leaf-position. Each arc contained 178 control points with 120 MLC positions. MATLAB routines were written to process these parameters in order to calculate both the beam aperture (unblocked) size for each control point. An average aperture size was weighted by meterset. Jaw factor was defined as the ratio between dynamic-jaw to static-jaw aperture size. Jaw factor was determined for forty retrospectively replanned patients treated with static-jaw delivery sites including lung, brain, prostate, H&N, rectum, and bladder. Results: Most patients had multiple arcs and reduced-field boosts, resulting in 151 fields. Of these, the lowest (0.4722) and highest (0.9622) jaw factor was observed in prostate and rectal cases, respectively. The median jaw factor was 0.7917 meaning there is the potential unincreased blocking by 20%. Clinically, the dynamic-jaw tracking represents an area surrounding the target which would receive MLC-only leakage transmission of 1.68% versus 0.1% with jaws. Jaw-tracking was more pronounced at areas farther from the target. In prostate patients, the rectum and bladder had 5.5% and 6.3% lower mean dose, respectively; the structures closer to the prostate such as the rectum and bladder both had 1.4% lower mean dose. Conclusion: A custom ESAPI script was coupled with a MATLAB routine in order to extract beam parameters from static-jaw plans and their replanned dynamic-jaw deliveries. The effects were quantified using jaw factor which is the ratio between the meterset weighted aperture size for dynamic-jaw fields versus static-jaw fields.

60 APPLIED LIFE SCIENCES↗

SU-G-BRB-04: Automated Output Factor Measurements Using Continuous Data Logging for Linac Commissioning

Purpose: Linac commissioning is a time consuming and labor intensive process, the streamline of which is highly desirable. In particular, manual measurement of output factors for a variety of field sizes and energy greatly hinders the commissioning efficiency. In this study, automated measurement of output factors was demonstrated as ‘one-click’ using data logging of an electrometer. Methods: Beams to be measured were created in the recording and verifying (R&V) system and configured for continuous delivery. An electrometer with an automatic data logging feature enabled continuous data collection for all fields without human intervention. The electrometer saved data into a spreadsheet every 0.5 seconds. A Matlab program was developed to analyze the excel data to monitor and check the data quality. Results: For each photon energy, output factors were measured for five configurations, including open field and four wedges. Each configuration includes 72 fields sizes, ranging from 4×4 to 20×30 cm{sup 2}. Using automation, it took 50 minutes to complete the measurement of 72 field sizes, in contrast to 80 minutes when using the manual approach. The automation avoided the necessity of redundant Linac status checks between fields as in the manual approach. In fact, the only limiting factor in such automation is Linac overheating. The data collection beams in the R&V system are reusable, and the simplified process is less error-prone. In addition, our Matlab program extracted the output factors faithfully from data logging, and the discrepancy between the automatic and manual measurement is within ±0.3%. For two separate automated measurements 30 days apart, consistency check shows a discrepancy within ±1% for 6MV photon with a 60 degree wedge. Conclusion: Automated output factor measurements can save time by 40% when compared with conventional manual approach. This work laid ground for further improvement for the automation of Linac commissioning.

43 PARTICLE ACCELERATORS↗

SU-G-TeP1-01: A Simulation Study to Investigate Maximum Allowable Deformations of Implant Geometry Before Plan Objectives Are Violated in Prostate HDR Brachytherapy

Purpose: In prostate HDR brachytherapy dose distributions are highly sensitive to changes in prostate volume and catheter displacements. We investigate the maximum deformations in implant geometry before planning objectives are violated. Methods: A typical prostate Ir-192 HDR brachytherapy reference plan was calculated on the Oncentra planning system, which used CT images from a tissue equivalent prostate phantom (CIRS Model 053S) embedded inside a pelvis wax phantom. The prostate was deformed and catheters were displaced in simulations using a code written in MATLAB. For each deformation dose distributions were calculated, based on TG43 methods, using the MATLAB code. The calculations were validated through comparison with Oncentra calculations for the reference plan, and agreed within 0.12%SD and 0.3%SD for dose and volume, respectively. Isotropic prostate volume deformations of up to +34% to −27% relative to its original volume, and longitudinal catheter displacements of 7.5 mm in superior and inferior directions were simulated. Planning objectives were based on American Brachytherapy Society guidelines for prostate and urethra volumes. A plan violated the planning objectives when less than 90% of the prostate volume received the prescribed dose or higher (V{sub 100}), or the urethral volume receiving 125% of prescribed dose or higher was more than 1 cc (U{sub 125}). Lastly, the dose homogeneity index (DHI=1-V{sub 150}/V{sub 100}) was evaluated; a plan was considered sub-optimal when the DHI fell below 0.62. Results and Conclusion: Planning objectives were violated when the prostate expanded by 10.7±0.5% or contracted by 11.0±0.2%; objectives were also violated when catheters were displaced by 4.15±0.15 mm and 3.70±0.15 mm in the superior and inferior directions, respectively. The DHI changes did not affect the plan optimality, except in the case of prostate compression. In general, catheter displacements have a significantly larger impact on plan optimality than prostate volume changes.

60 APPLIED LIFE SCIENCES↗