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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-E-06: Teaching About the Gamma Camera and Ultrasound Imaging

Purpose: Instructional modules on applications of physics in medicine are being developed. The target audience consists of students who have had an introductory undergraduate physics course. This presentation will concentrate on an active learning approach to teach the principles of the gamma camera. There will also be a description of an apparatus to teach ultrasound imaging. Methods: Since a real gamma camera is not feasible in the undergraduate classroom, we have developed two types of optical apparatus that teach the main principles. To understand the collimator, LEDS mimic gamma emitters in the body, and the photons pass through an array of tubes. The distance, spacing, diameter, and length of the tubes can be varied to understand the effect upon the resolution of the image. To determine the positions of the gamma emitters, a second apparatus uses a movable green laser, fluorescent plastic in lieu of the scintillation crystal, acrylic rods that mimic the PMTs, and a photodetector to measure the intensity. The position of the laser is calculated with a centroid algorithm.To teach the principles of ultrasound imaging, we are using the sound head and pulser box of an educational product, variable gain amplifier, rotation table, digital oscilloscope, Matlab software, and phantoms. Results: Gamma camera curriculum materials have been implemented in the classroom at Loyola in 2014 and 2015. Written work shows good knowledge retention and a more complete understanding of the material. Preliminary ultrasound imaging materials were run in 2015. Conclusion: Active learning methods add another dimension to descriptions in textbooks and are effective in keeping the students engaged during class time. The teaching apparatus for the gamma camera and ultrasound imaging can be expanded to include more cases, and could potentially improve students’ understanding of artifacts and distortions in the images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-E-07: An Adaptable Approach for Education On Medical Physics at Undergraduate and Postgraduate Levels

Purpose: To contribute to the professional profile of future medical physicists, technologists and physicians, and implement an adaptable educational strategy at both undergraduate and postgraduate levels. Methods: The Medical Physics Block of Electives (MPBE) designed was adapted to the Program of B.S. in Physics. The conferences and practical activities were developed with participatory methods, with interdisciplinary collaboration from research institutions and hospitals engaged on projects of Research, Development and Innovation (RDI). The scientific education was implemented by means of critical analysis of scientific papers and seminars where students debated on solutions for real research problems faced by medical physicists. This approach included courses for graduates not associated to educational programs of Medical Physics (MP). Results: The implementation of the MPBE began in September 2014, with the electives of Radiation MP and Introduction to Nuclear Magnetic Resonance. The students of second year received an Introduction to MP. This initiative was validated by the departmental Methodological Workshop, which promoted the full implementation of the MPBE. Both postgraduated and undergraduate trainees participated in practices with our DICOM viewer system, a local prototype for photoplethysmography and a home-made interface for ROC analysis, built with MATLAB. All these tools were designed and constructed in previous RDI projects. The collaborative supervision of University’s researchers with clinical medical physicists will allow to overcome the limitations of residency in hospitals, to reduce the workload for clinical supervisors and develop appropriate educational activities. Conclusion: We demonstrated the feasibility of adaptable educational strategies, considering available resources. This provides an innovative way for prospective medical physicists, technologists and radiation oncologists. This strategy can be implemented in several regions without formal programs of MP, like most of developing countries. Starting with undergraduate students would allow to reach appropriate certification faster than most of traditional or alternative approaches for education on MP. The authors acknowledge Radiation Consulting Group, LLC, an Arizona Corporation which promotes the use of ionizing radiation in the healing arts, for the “Oscar Luis Caballero” travel grant. The authors thanks to professors Meisbel Daudinot, David Adame and Alexander Pascau for the practices through imagis, imageROC and ANGIODIN PD3000 respectively.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-P-26: Oncospace: A Shared Radiation Oncology Database System Designed for Personalized Medicine, Decision Support, and Research

Purpose: Advancement in Radiation Oncology (RO) practice develops through evidence based medicine and clinical trial. Knowledge usable for treatment planning, decision support and research is contained in our clinical data, stored in an Oncospace database. This data store and the tools for populating and analyzing it are compatible with standard RO practice and are shared with collaborating institutions. The question is - what protocol for system development and data sharing within an Oncospace Consortium? We focus our example on the technology and data meaning necessary to share across the Consortium. Methods: Oncospace consists of a database schema, planning and outcome data import and web based analysis tools.1) Database: The Consortium implements a federated data store; each member collects and maintains its own data within an Oncospace schema. For privacy, PHI is contained within a single table, accessible to the database owner.2) Import: Spatial dose data from treatment plans (Pinnacle or DICOM) is imported via Oncolink. Treatment outcomes are imported from an OIS (MOSAIQ).3) Analysis: JHU has built a number of webpages to answer analysis questions. Oncospace data can also be analyzed via MATLAB or SAS queries.These materials are available to Consortium members, who contribute enhancements and improvements. Results: 1) The Oncospace Consortium now consists of RO centers at JHU, UVA, UW and the University of Toronto. These members have successfully installed and populated Oncospace databases with over 1000 patients collectively.2) Members contributing code and getting updates via SVN repository. Errors are reported and tracked via Redmine. Teleconferences include strategizing design and code reviews.3) Successfully remotely queried federated databases to combine multiple institutions’ DVH data for dose-toxicity analysis (see below – data combined from JHU and UW Oncospace). Conclusion: RO data sharing can and has been effected according to the Oncospace Consortium model: http://oncospace.radonc.jhmi.edu/ . John Wong - SRA from Elekta; Todd McNutt - SRA from Elekta; Michael Bowers - funded by Elekta.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-P-36: Evaluation of MLC Positioning Errors in Dynamic IMRT Treatments by Analyzing Dynalog Files

Purpose: To assess the accuracy of MLC positioning in Varian linear accelerator, in dynamic IMRT technique, from the analysis of dynalog files generated by the MLC controller. Methods: In Clinac accelerators (pre-TrueBeam technology), control system has an approximately 50ms delay (one control cycle time). Then, the system compares the measured position to the planned position corresponding to the next control cycle. As it has been confirmed by Varian technical support, this effect causes that measured positions appear in dynalogs one cycle out of phase with respect to the planned positions. Around 9000 dynalogs have been analyzed, coming from the three linear accelerators of our center (one Trilogy and two Clinac 21EX) equipped with a Millennium 120 MLC. In order to compare our results to recent publications, leaf positioning errors (RMS and 95th percentile) are calculated with and without delay effect. Dynalogs have been analyzed using a in-house Matlab software. Results: The RMS errors were 0.341, 0.339 and 0.348mm for each Linac; being the average error 0.343 mm. The 95th percentiles of the error were 0.617, 0.607 and 0.625; with an average of 0.617mm. A recent multi-institution study carried out by Kerns et al. found a mean leaf RMS error of 0.32mm and a 95th percentile error value of 0.64mm.Without delay effect, mean leaf RMS errors obtained were 0.040, 0.042 and 0.038mm for each treatment machine; being the average 0.040mm. The 95th percentile error values obtained were 0.057, 0.058 and 0.054 mm, with an average of 0.056mm. Conclusion: Results obtained for the mean leaf RMS error and the mean 95th percentile were consistent with the multi-institution study. Calculated error statistics with delay effect are significantly larger due to the speed proportional and systematic leaf offset. Consequently it is proposed to correct this effect in dynalogs analysis to determine the MLC performance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-P-59: A Graphical Interface for XCAT Phantom Configuration, Generation and Processing

Purpose: To design a comprehensive open-source, publicly available, graphical user interface (GUI) to facilitate the configuration, generation, processing and use of the 4D Extended Cardiac-Torso (XCAT) phantom. Methods: The XCAT phantom includes over 9000 anatomical objects as well as respiratory, cardiac and tumor motion. It is widely used for research studies in medical imaging and radiotherapy. The phantom generation process involves the configuration of a text script to parameterize the geometry, motion, and composition of the whole body and objects within it, and to generate simulated PET or CT images. To avoid the need for manual editing or script writing, our MATLAB-based GUI uses slider controls, drop-down lists, buttons and graphical text input to parameterize and process the phantom. Results: Our GUI can be used to: a) generate parameter files; b) generate the voxelized phantom; c) combine the phantom with a lesion; d) display the phantom; e) produce average and maximum intensity images from the phantom output files; f) incorporate irregular patient breathing patterns; and f) generate DICOM files containing phantom images. The GUI provides local help information using tool-tip strings on the currently selected phantom, minimizing the need for external documentation. The DICOM generation feature is intended to simplify the process of importing the phantom images into radiotherapy treatment planning systems or other clinical software. Conclusion: The GUI simplifies and automates the use of the XCAT phantom for imaging-based research projects in medical imaging or radiotherapy. This has the potential to accelerate research conducted with the XCAT phantom, or to ease the learning curve for new users. This tool does not include the XCAT phantom software itself. We would like to acknowledge funding from MRA, Varian Medical Systems Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-I-45: Feasibility for Using Iodine Quantification to Assist Diagnosis in Dual Energy Contrast-Enhanced Digital Mammography

Purpose: The objective of this study is to develop quantitative calibration between image quality indexes and iodine concentration with dual-energy (DE) contrast-enhanced digital mammography (CEDM) techniques and further serve as the assistance for diagnosis. Methods: Custom-made acrylic phantom with dimensions of 24×30 cm{sup 2} simulated breast thickness from 2 to 6 cm was used in the calibration. The phantom contained matrix of four times four holes of 3 mm deep with a diameter of 15 mm for filling contrast agent with area density ranged from 0.1 to 10 mg/cm{sup 2}. All the image acquisitions were performed on a full-field digital mammography system (Senographe Essential, GE) with dual energy acquisitions. Mean pixel value (MPV), and contrast-to-noise ratio (CNR) were used for evaluating the relationship between image quality indexes and iodine concentration. Iodine map and CNR map could further be constructed with these calibration curves applied pixel by pixel utilized MATLAB software. Minimum iodine concentration could also be calculated with the visibility threshold of CNR=5 according the Rose model. Results: When evaluating the DE subtraction images, MPV increased linearly as the iodine concentration increased with all the phantom thickness surveyed (R{sup 2} between 0.989 and 0.992). Lesions with increased iodine uptake could thus be enhanced in the color-encoded iodine maps, and the mean iodine concentration could be obtained through the ROI measurements. As for investigating CNR performance, linear relationships were also shown between the iodine concentration and CNR (R{sup 2} between 0.983 and 0.990). Minimum iodine area density of 1.45, 1.73, 1.80, 1.73 and 1.72 mg/cm{sup 2} for phantom thickness of 2, 3, 4, 5, 6 cm were calculated based on Rose’s visualization criteria. Conclusion: Quantitative calibration between image quality indexes and iodine concentrations may further serving as the assistance for analyzing contrast enhancement for patient participating the dual energy CEDM procedures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-14: A Comparison of a 2.5MV Imaging Beam to KV and 6MV Imaging Beams

Purpose: To compare image quality metrics and dose of TrueBeam V2.0’s 2.5MV imaging beam and kV and 6MV images. Methods: To evaluate the MV image quality, the Standard Imaging QC-3 and Varian Las Vegas (LV) phantoms were imaged using the ‘quality’ and ‘low dose’ modes and then processed using RIT113 V6.3. The LEEDS phantom was used to evaluate the kV image quality. The signal to noise ratio (SNR) was also evaluated in patient images using Matlab. In addition, dose per image was evaluated at a depth of 5cm using solid water for a 28.6 cm × 28.6 cm field size, which is representative of the largest jaw settings at an SID of 150cm. Results: The 2.5MV images had lower dose than the 6 MV images and a contrast to noise ratio (CNR) about 1.4 times higher, when evaluated using the QC-3. When energy was held constant but dose varied, the different modes, ‘low dose’ and ‘quality’, showed less than an 8% difference in CNR. The ‘quality’ modes demonstrated better spatial resolution than the ‘low dose’; however, even with the ‘low dose’ all line pairs were distinct except for the 0.75lp/mm on the 2.5MV. The LV phantom was used to measure low contrast detectability and showed similar results to the QC-3. Several patient images all confirmed that SNR were highest in kV images followed by 2.5MV and then 6MV. Qualitatively, for anatomical areas with large variability in thickness, like lateral head and necks, 2.5MV images show more anatomy, such as shoulder position, than kV images. Conclusions: The kV images clearly provide the best image metrics per unit dose. The 2.5MV beam showed excellent contrast at a lower dose than 6MV and may be superior to kV for difficult to image areas that include large changes in anatomical thickness. P Balter: Varian, Sun Nuclear, Philips, CPRIT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-61: Monitoring Tumor Motion in Real-Time with EPID Imaging During Cervical Cancer Treatment

Purpose: Traditional external beam radiotherapy for cervical cancer requires setup by external skin marks. In order to improve treatment accuracy and reduce planning margin for more conformal therapy, it is essential to monitor tumor positions interfractionally and intrafractionally. We demonstrate feasibility of monitoring cervical tumor motion online using EPID imaging from Beam’s Eye View. Methods: Prior to treatment, 1∼2 cylindrical radio opaque markers were implanted into inferior aspect of cervix tumor. During external beam treatments on a Varian 2100C by 4-field 3D plans, treatment beam images were acquired continuously by an EPID. A Matlab program was developed to locate internal markers on MV images. Based on 2D marker positions obtained from different treatment fields, their 3D positions were estimated for every treatment fraction. Results: There were 398 images acquired during different treatment fractions of three cervical cancer patients. Markers were successfully located on every frame of image at an analysis speed of about 1 second per frame. Intrafraction motions were evaluated by comparing marker positions relative to the position on the first frame of image. The maximum intrafraction motion of the markers was 1.6 mm. Interfraction motions were evaluated by comparing 3D marker positions at different treatment fractions. The maximum interfraction motion was up to 10 mm. Careful comparison found that this is due to patient positioning since the bony structures shifted with the markers. Conclusion: This method provides a cost-free and simple solution for online tumor tracking for cervical cancer treatment since it is feasible to acquire and export EPID images with fast analysis in real time. This method does not need any extra equipment or deliver extra dose to patients. The online tumor motion information will be very useful to reduce planning margins and improve treatment accuracy, which is particularly important for SBRT treatment with long delivery time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-93: Development of Pre-Contoured Human Model Library in DICOM-RT Format for the Epidemiological Study of the Radiotherapy Patients

Purpose: Prior to 3D conformal radiation therapy planning, patient anatomy information was mostly limited to 2D beams-eye-view from the conventional simulator. To analyze the outcomes of such treatments for radiation late effects, 3D computational human models are often used in commercial treatment planning systems (TPSs). However, several underlying difficulties such as time-consuming manual delineation procedures of a large number of structures in the model have always limited its applications. Primary objective of this work was to develop a human model library for the epidemiological study by converting 3D-surface model organs to DICOM-RT format (DICOM-RT structure) using an in-house built software. We converted the ICRP reference human models to DICOM-RT models, which can be readily adopted for various dose calculations. Methods: MATLAB based code were utilized to convert the contour drawings extracted in text-format from the 3D graphic-tool, Rhinoceros into DICOM-RT structure format for 50 different organs of each model using a 16GB dual-core processor. The conversion periods were measured for each DICOM-RT models, and the reconstructed structure volumes were validated against the original 3D-surface models in the TPS. Ten reference hybrid whole-body models (8-pediatric and 2-adults) were automatically processed to create DICOM-RT computational human model library. Results: Mean contour conversion period was found to be 580 (N=2) and 394.5 (N=8) seconds for 50 organs in the adult and pediatric models respectively. A good agreement for large organs (NRMSD <1.0%) and small organs (NRMSD <7.7%) was also observed between the original volumes and corresponding DICOM-RT structure volumes of the organs. Conclusion: The ICRP reference human models were converted into DICOM-RT format to support the epidemiological study using a large cohort of conventional radiotherapy patients. Due to its DICOM-compatibility, the library may be implemented to many other different applications. We also expect to develop the library by including additional models in future.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-111: Finite Element-Based Deformable Image Registration of Pleural Cavity for Photodynamic Therapy

Purpose: The pleural volumes will deform during surgery portion of the pleural photodynamic therapy (PDT) of lung cancer when the pleural cavity is opened. This impact the delivered dose when using highly conformal treatment techniques. In this study, a finite element-based (FEM) deformable image registration is used to quantify the anatomical variation between the contours for the pleural cavities obtained in the operating room and those determined from pre-surgery computed tomography (CT) scans. Methods: An infrared camera-based navigation system (NDI) is used during PDT to track the anatomical changes and contour the lung and chest cavity. A series of CTs of the lungs, in the same patient, are also acquired before the surgery. The structure contour of lung and the CTs are processed and contoured in Matlab and MeshLab. Then, the contours are imported into COMSOL Multiphysics 5.0, where the FEM-based deformable image registration is obtained using the deformed mesh - moving mesh (ALE) model. The NDI acquired lung contour is considered as the reference contour, and the CT contour is used as the target one, which will be deformed. Results: The reconstructed three-dimensional contours from both NDI and CT can be converted to COMSOL so that a three-dimensional ALE model can be developed. The contours can be registered using COMSOL ALE moving mesh model, which takes into account the deformation along x, y and z-axes. The deformed contour has good matches to the reference contour after the dynamic matching process. The resulting 3D deformation map can be used to obtain the locations of other critical anatomic structures, e.g., heart, during surgery. Conclusion: Deformable image registration can fuse images acquired by different modalities. It provides insights into the development of phenomenon and variation in normal anatomical structures over time. The initial assessments of three-dimensional registration show good agreement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-163: A Biomechanical Lung Model for Respiratory Motion Study

Purpose: This work presents a biomechanical model to investigate the complex respiratory motion for the lung tumor tracking in radiosurgery by computer simulation. Methods: The models include networked massspring-dampers to describe the tumor motion, different types of surrogate signals, and the force generated by the diaphragm. Each mass-springdamper has the same mechanical structure and each model can have different numbers of mass-spring-dampers. Both linear and nonlinear stiffness parameters were considered, and the damping ratio was tuned in a range so that the tumor motion was over-damped (no natural tumor oscillation occurs without force from the diaphragm). The simulation was run by using ODE45 (ordinary differential equations by Runge-Kutta method) in MATLAB, and all time courses of motions and inputs (force) were generated and compared. Results: The curvature of the motion time courses around their peaks was sensitive to the damping ratio. Therefore, the damping ratio can be determined based on the clinical data of a high sampling rate. The peak values of different signals and the time the peaks occurred were compared, and it was found that the diaphragm force had a time lead over the tumor motion, and the lead time (0.1–0.4 seconds) depended on the distance between the tumor and the diaphragm. Conclusion: We reported a model based analysis approach for the spatial and temporal relation between the motion of the lung tumor and the surrogate signals. Due to the phase lead of the diaphragm in comparing with the lung tumor motion, the measurement of diaphragm motion (or its electromyography signal) can be used as a beam gating signal in radiosurgery, and it can also be an additional surrogate signal for better tumor motion tracking. The research is funded by the American Cancer Society (ACS) grant. The grant name is: Frameless SRS Based on Robotic Head Motion Cancellation. The grant number is: RSG-13-313-01-CCE.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-164: Estimation of DVH Variation for PTV Due to Interfraction Organ Motion in Prostate VMAT Using Gaussian Error Function

Purpose: We developed a method to predict the change of DVH for PTV due to interfraction organ motion in prostate VMAT without repeating the CT scan and treatment planning. The method is based on a pre-calculated patient database with DVH curves of PTV modelled by the Gaussian error function (GEF). Methods: For a group of 30 patients with different prostate sizes, their VMAT plans were recalculated by shifting their PTVs 1 cm with 10 increments in the anterior-posterior, left-right and superior-inferior directions. The DVH curve of PTV in each replan was then fitted by the GEF to determine parameters describing the shape of curve. Information of parameters, varying with the DVH change due to prostate motion for different prostate sizes, was analyzed and stored in a database of a program written by MATLAB. Results: To predict a new DVH for PTV due to prostate interfraction motion, prostate size and shift distance with direction were input to the program. Parameters modelling the DVH for PTV were determined based on the pre-calculated patient dataset. From the new parameters, DVH curves of PTVs with and without considering the prostate motion were plotted for comparison. The program was verified with different prostate cases involving interfraction prostate shifts and replans. Conclusion: Variation of DVH for PTV in prostate VMAT can be predicted using a pre-calculated patient database with DVH curve fitting. The computing time is fast because CT rescan and replan are not required. This quick DVH estimation can help radiation staff to determine if the changed PTV coverage due to prostate shift is tolerable in the treatment. However, it should be noted that the program can only consider prostate interfraction motions along three axes, and is restricted to prostate VMAT plan using the same plan script in the treatment planning system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-169: The Dosimetric and Temporal Effects of Respiratory-Gated Radiation Therapy in Lung Cancer Patients

Purpose: To evaluate the dosimetric and temporal effects of high dose rate treatment mode for respiratory-gated radiation therapy in lung cancer patients. Methods: Treatment plans from five lung cancer patients (3 nongated (Group 1), 2 gated at 80EX-80IN (Group 2)) were retrospectively evaluated. The maximum tumor motions range from 6–12 mm. Using the same planning criteria, four new treatment plans, corresponding to four gating windows (20EX–20IN, 40EX–40IN, 60EX–60IN, and 80EX–80IN), were generated for each patient. Mean tumor dose (MTD), mean lung dose (MLD), and lung V20 were used to assess the dosimetric effects. A MATLAB algorithm was developed to compute treatment time by considering gantry rotation time, time to position collimator leaves, dose delivery time (scaled relative to the gating window), and communication overhead. Treatment delivery time for each plan was estimated using a 500 MU/min dose rate for the original plans and a 1500 MU/min dose rate for the gated plans. Results: Differences in MTD were less than 1Gy across plans for all five patients. MLD and lung V20 were on average reduced between −16.1% to −6.0% and −20.0% to −7.2%, respectively for non-gated plans when compared with the corresponding gated plans, and between − 5.8% to −4.2% and −7.0% to −5.4%, respectively for plans originally gated at 80EX–80IN when compared with the corresponding 20EX-20IN to 60EX– 60IN gated plans. Treatment delivery times of gated plans using high dose rate were reduced on average between −19.7% (−1.9min) to −27.2% (−2.7min) for originally non-gated plans and −15.6% (−0.9min) to −20.3% (−1.2min) for originally 80EX-80IN gated plans. Conclusion: Respiratory-gated radiation therapy in lung cancer patients can reduce lung toxicity, while maintaining tumor dose. Using a gated high-dose-rate treatment, delivery time comparable to non-gated normal-dose-rate treatment can be achieved. This research is supported by Siemens Medical Solutions USA, Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-188: Theoretical Estimation of Margin Necessary for Markerless Motion Tracking

Purpose: To estimate the margin necessary to adequately cover the target using markerless motion tracking (MMT) of lung lesions given the uncertainty in tracking and the size of the target. Methods: Simulations were developed in Matlab to determine the effect of tumor size and tracking uncertainty on the margin necessary to achieve adequate coverage of the target. For simplicity, the lung tumor was approximated by a circle on a 2D radiograph. The tumor was varied in size from a diameter of 0.1 − 30 mm in increments of 0.1 mm. From our previous studies using dual energy markerless motion tracking, we estimated tracking uncertainties in x and y to have a standard deviation of 2 mm. A Gaussian was used to simulate the deviation between the tracked location and true target location. For each size tumor, 100,000 deviations were randomly generated, the margin necessary to achieve at least 95% coverage 95% of the time was recorded. Additional simulations were run for varying uncertainties to demonstrate the effect of the tracking accuracy on the margin size. Results: The simulations showed an inverse relationship between tumor size and margin necessary to achieve 95% coverage 95% of the time using the MMT technique. The margin decreased exponentially with target size. An increase in tracking accuracy expectedly showed a decrease in margin size as well. Conclusion: In our clinic a 5 mm expansion of the internal target volume (ITV) is used to define the planning target volume (PTV). These simulations show that for tracking accuracies in x and y better than 2 mm, the margin required is less than 5 mm. This simple simulation can provide physicians with a guideline estimation for the margin necessary for use of MMT clinically based on the accuracy of their tracking and the size of the tumor.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-253: The Radiomics Toolbox in the Computational Environment for Radiological Research (CERR)

Purpose: To present an open source and free platform to facilitate radiomics research — The “Radiomics toolbox” in CERR. Method: There is scarcity of open source tools that support end-to-end modeling of image features to predict patient outcomes. The “Radiomics toolbox” strives to fill the need for such a software platform. The platform supports (1) import of various kinds of image modalities like CT, PET, MR, SPECT, US. (2) Contouring tools to delineate structures of interest. (3) Extraction and storage of image based features like 1st order statistics, gray-scale co-occurrence and zonesize matrix based texture features and shape features and (4) Statistical Analysis. Statistical analysis of the extracted features is supported with basic functionality that includes univariate correlations, Kaplan-Meir curves and advanced functionality that includes feature reduction and multivariate modeling. The graphical user interface and the data management are performed with Matlab for the ease of development and readability of code and features for wide audience. Open-source software developed with other programming languages is integrated to enhance various components of this toolbox. For example: Java-based DCM4CHE for import of DICOM, R for statistical analysis. Results: The Radiomics toolbox will be distributed as an open source, GNU copyrighted software. The toolbox was prototyped for modeling Oropharyngeal PET dataset at MSKCC. The analysis will be presented in a separate paper. Conclusion: The Radiomics Toolbox provides an extensible platform for extracting and modeling image features. To emphasize new uses of CERR for radiomics and image-based research, we have changed the name from the “Computational Environment for Radiotherapy Research” to the “Computational Environment for Radiological Research”.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-262: Variability in Texture Analysis of Gynecological Tumors in the Context of An 18F-FDG PET Adaptive Protocol

Purpose: This study examines the effect on texture analysis due to variable reconstruction of PET images in the context of an adaptive FDG PET protocol for node positive gynecologic cancer patients. By measuring variability in texture features from baseline and intra-treatment PET-CT, we can isolate unreliable texture features due to large variation. Methods: A subset of seven patients with node positive gynecological cancers visible on PET was selected for this study. Prescribed dose varied between 45–50.4Gy, with a 55–70Gy boost to the PET positive nodes. A baseline and intratreatment (between 30–36Gy) PET-CT were obtained on a Siemens Biograph mCT. Each clinical PET image set was reconstructed 6 times using a TrueX+TOF algorithm with varying iterations and Gaussian filter. Baseline and intra-treatment primary GTVs were segmented using PET Edge (MIM Software Inc., Cleveland, OH), a semi-automatic gradient-based algorithm, on the clinical PET and transferred to the other reconstructed sets. Using an in-house MATLAB program, four 3D texture matrices describing relationships between voxel intensities in the GTV were generated: co-occurrence, run length, size zone, and neighborhood difference. From these, 39 textural features characterizing texture were calculated in addition to SUV histogram features. The percent variability among parameters was first calculated. Each reconstructed texture feature from baseline and intra-treatment per patient was normalized to the clinical baseline scan and compared using the Wilcoxon signed-rank test in order to isolate variations due to reconstruction parameters. Results: For the baseline scans, 13 texture features showed a mean range greater than 10%. For the intra scans, 28 texture features showed a mean range greater than 10%. Comparing baseline to intra scans, 25 texture features showed p <0.05. Conclusion: Variability due to different reconstruction parameters increased with treatment, however, the majority of texture features showed significant changes during treatment independent of reconstruction effects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-268: Change of CT Number During the Course of Chemoradiation Therapy for Pancreatic Cancer

Purpose: It has been observed radiation can induce changes in CT number (CTN) inside tumor during the course of radiation therapy (RT) for several tumor sites including lung and head and neck, suggesting that the CTN change may be potentially used to assess RT response. In this study, we investigate the CTN changes inside tumor during the course of chemoradiation therapy (CRT) for pancreatic cancer. Methods: Daily diagnostic-quality CT data acquired during IGRT for 17 pancreatic head cancer patients using an in-room CT (CTVision, Siemens) were analyzed. All patients were treated with a radiation dose of 50.4 in 1.8 Gy per fraction. On each daily CT set, The contour of the pancreatic head, included in the treatment target, was generated by populating the pancreatic head contour from the planning CT or MRI using an auto-segmentation tool based on deformable registration (ABAS, Elekta) with manual editing if necessary. The CTN at each voxel in the pancreatic head contour was extracted and the 3D distribution of the CTNs was processed using MATLAB. The mean value of CTN distribution was used to quantify the daily CTN change in the pancreatic head. Results: Reduction of CTN in pancreatic head was observed during the CRT delivery in 14 out the 17 (82%) patients studied. Although the average reduction is only 3.5 Houncefield Unit (HU), this change is significant (p<0.01). Among them, there are 7 patients who had a CTN drop larger than 5 HU, ranging from 6.0 to 11.8 HU. In contrast to this trend, CTN was increased in 3 patients. Conclusion: Measurable changes in the CT number in tumor target were observed during the course of chemoradiation therapy for the pancreas cancer patients, indicating this radiation-induced CTN change may be used to assess treatment response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-20: A Correlation Study of 2D and 3D Gamma Passing Rates for Prostate IMRT Plans

Purpose: To investigate the correlation between the two-dimensional gamma passing rate (2D %GP) and three-dimensional gamma passing rate (3D %GP) in prostate IMRT quality assurance. Methods: Eleven prostate IMRT plans were randomly selected from the clinical database and were used to obtain dose distributions in the phantom and patient. Three types of delivery errors (MLC bank sag errors, central MLC errors and monitor unit errors) were intentionally introduced to modify the clinical plans through an in-house Matlab program. This resulted in 187 modified plans. The 2D %GP and 3D %GP were analyzed using different dose-difference and distance-toagreement (1%-1mm, 2%-2mm and 3%-3mm) and 20% dose threshold. The 2D %GP and 3D %GP were then compared not only for the whole region, but also for the PTVs and critical structures using the statistical Pearson’s correlation coefficient (γ). Results: For different delivery errors, the average comparison of 2D %GP and 3D %GP showed different conclusions. The statistical correlation coefficients between 2D %GP and 3D %GP for the whole dose distribution showed that except for 3%/3mm criterion, 2D %GP and 3D %GP of 1%/1mm criterion and 2%/2mm criterion had strong correlations (Pearson’s γ value >0.8). Compared with the whole region, the correlations of 2D %GP and 3D %GP for PTV were better (the γ value for 1%/1mm, 2%/2mm and 3%/3mm criterion was 0.959, 0.931 and 0.855, respectively). However for the rectum, there was no correlation between 2D %GP and 3D %GP. Conclusion: For prostate IMRT, the correlation between 2D %GP and 3D %GP for the PTV is better than that for normal structures. The lower dose-difference and DTA criterion shows less difference between 2D %GP and 3D %GP. Other factors such as the dosimeter characteristics and TPS algorithm bias may also influence the correlation between 2D %GP and 3D %GP.

60 APPLIED LIFE SCIENCES↗