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

SU-G-BRA-02: Development of a Learning Based Block Matching Algorithm for Ultrasound Tracking in Radiotherapy

Purpose: To develop an ultrasound learning-based tracking algorithm with the potential to provide real-time motion traces of anatomy-based fiducials that may aid in the effective delivery of external beam radiation. Methods: The algorithm was developed in Matlab R2015a and consists of two main stages: reference frame selection, and localized block matching. Immediately following frame acquisition, a normalized cross-correlation (NCC) similarity metric is used to determine a reference frame most similar to the current frame from a series of training set images that were acquired during a pretreatment scan. Segmented features in the reference frame provide the basis for the localized block matching to determine the feature locations in the current frame. The boundary points of the reference frame segmentation are used as the initial locations for the block matching and NCC is used to find the most similar block in the current frame. The best matched block locations in the current frame comprise the updated feature boundary. The algorithm was tested using five features from two sets of ultrasound patient data obtained from MICCAI 2014 CLUST. Due to the lack of a training set associated with the image sequences, the first 200 frames of the image sets were considered a valid training set for preliminary testing, and tracking was performed over the remaining frames. Results: Tracking of the five vessel features resulted in an average tracking error of 1.21 mm relative to predefined annotations. The average analysis rate was 15.7 FPS with analysis for one of the two patients reaching real-time speeds. Computations were performed on an i5-3230M at 2.60 GHz. Conclusion: Preliminary tests show tracking errors comparable with similar algorithms at close to real-time speeds. Extension of the work onto a GPU platform has the potential to achieve real-time performance, making tracking for therapy applications a feasible option. This work is partially funded by NIH grant R01CA190298.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-IeP1-08: MR Geometric Distortion Dependency On Imaging Sequence, Acquisition Orientation and Receiver Bandwidth of a Dedicated 1.5T MR-Simulator

Purpose: To investigate the 3D geometric distortion of four potential MR sequences for radiotheraptic applications, and its dependency on sequence-type, acquisition-orientation and receiver-bandwidth from a dedicated 1.5T 700mm-wide bore MR-simulator (Magnetom-Aera, Sienmens Healthcare, Erlangen, Germany), using a large customized geometric accuracy phantom. Methods: This work studied 3D gradient-echo (VIBE) and spin-echo (SPACE) sequences for anatomical imaging; a specific ultra-short-TE sequence (PETRA) potentially for bone imaging and MR-based dosimetry; and a motion-insensitive sequence (BLADE) for dynamic applications like 4D-MRI. Integrated geometric-correction was employed, three orthogonal acquisition-orientations and up to three receiver-bandwidths were used, yielding 27 acquisitions for testing (Table 1a).A customized geometric accuracy phantom (polyurethane, MR/CT invisible, W×L×H:55×55×32.5cm3) was constructed and filled with 3892 spherical markers (6mm diameter, MR/CT visible) arranged on a 25mm-interval 3D isotropic-grid (Fig.1). The marker positions in MR images were quantitatively calculated and compared against those in the CT-reference using customized MatLab scripts. Results: The average distortion within various diameter-of-spherical-volumes (DSVs) and the usable DSVs under various distortion limits were measured (Tables 1b-c). It was observed that distortions fluctuated when sequence-type, acquisition-orientation or receiver-bandwidth changed (e.g. within 300mm-DSV, the lowest/highest average distortions of VIBE were 0.40mm/0.59mm, a 47.5% difference). According to AAPM-TG66 (<1mm distortion, left-most column of Table 1c), PETRA (Largest-DSV:253.9mm) has the potential on brain treatment, while BLADE (Largest-DSV:207.2mm) may need improvement for thoracic/abdominal applications. The results of VIBE (Largest-DSVs:294.3mm, the best among tested acquisitions) and SPACE (Largest-DSVs:267.7mm) suggests their potentials on head and neck applications. These Largest-DSVs were attained on different acquisition-orientations and receiver-bandwidths. Conclusion: Geometric distortion was shown to be dependent on sequence-type, acquisition-orientation and receiver-bandwidth. In the experiment, no configuration in any one of these factors could consistently reduce distortion while the others were varying. The distortion analysis result is a valuable guideline for sequence selection and optimization for MR-aided radiotherapy applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-JeP2-13: Spatial Accuracy Evaluation for Real-Time MR Guided Radiation Therapy Using a Novel Large-Field MRI Distortion Phantom

Purpose: Evaluate a large-field MRI phantom for assessment of geometric distortion in whole-body MRI for real-time MR guided radiation therapy. Methods: A prototype CIRS large-field MRI distortion phantom consisting of a PMMA cylinder (33 cm diameter, 30 cm length) containing a 3D-printed orthogonal grid (3 mm diameter rods, 20 mm apart), was filled with 6 mM NiCl{sub 2} and 30 mM NaCl solution. The phantom was scanned at 1.5T and 3.0T on a GE HDxt and Discovery MR750, respectively, and at 0.35T on a ViewRay system. Scans were obtained with and without 3D distortion correction to demonstrate the impact of such corrections. CT images were used as a reference standard for analysis of geometric distortion, as determined by a fully automated gradient-search method developed in Matlab. Results: 1,116 grid points distributed throughout a cylindrical volume 28 cm in diameter and 16 cm in length were identified and analyzed. With 3D distortion correction, average/maximum displacements for the 1.5, 3.0, and 0.35T systems were 0.84/2.91, 1.00/2.97, and 0.95/2.37 mm, respectively. The percentage of points with less than (1.0, 1.5, 2.0 mm) total displacement were (73%, 92%, 97%), (54%, 85%, 97%), and (55%, 90%, 99%), respectively. A reduced scan volume of 20 × 20 × 10 cm{sup 3} (representative of a head and neck scan volume) consisting of 420 points was also analyzed. In this volume, the percentage of points with less than (1.0, 1.5, 2.0 mm) total displacement were (90%, 99%, 100%), (63%, 95%, 100%), and (75%, 96%, 100%), respectively. Without 3D distortion correction, average/maximum displacements were 1.35/3.67, 1.67/4.46, and 1.51/3.89 mm, respectively. Conclusion: The prototype large-field MRI distortion phantom and developed software provide a thorough assessment of 3D spatial distortions in MRI. The distortions measured were acceptable for RT applications, both for the high field strengths and the system configuration developed by ViewRay.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-JeP3-09: Tumor Location Prediction Using Natural Respiratory Volume for Respiratory Gated Radiation Therapy (RGRT): System Verification Study

Purpose: Respiratory gated radiation therapy (RGRT) gives accurate results when a patient’s breathing is stable and regular. Thus, the patient should be fully aware during respiratory pattern training before undergoing the RGRT treatment. In order to bypass the process of respiratory pattern training, we propose a target location prediction system for RGRT that uses only natural respiratory volume, and confirm its application. Methods: In order to verify the proposed target location prediction system, an in-house phantom set was used. This set involves a chest phantom including target, external markers, and motion generator. Natural respiratory volume signals were generated using the random function in MATLAB code. In the chest phantom, the target takes a linear motion based on the respiratory signal. After a four-dimensional computed tomography (4DCT) scan of the in-house phantom, the motion trajectory was derived as a linear equation. The accuracy of the linear equation was compared with that of the motion algorithm used by the operating motion generator. In addition, we attempted target location prediction using random respiratory volume values. Results: The correspondence rate of the linear equation derived from the 4DCT images with the motion algorithm of the motion generator was 99.41%. In addition, the average error rate of target location prediction was 1.23% for 26 cases. Conclusion: We confirmed the applicability of our proposed target location prediction system for RGRT using natural respiratory volume. If additional clinical studies can be conducted, a more accurate prediction system can be realized without requiring respiratory pattern training.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-TeP4-02: A Method for Evaluating the Direct Impact of Failed IMRT QAs On Patient Dose

Purpose: We developed a method to calculate patient doses corresponding to IMRT QA measurements in order to determine and assess the actual dose delivered for plans with failed (or borderline) IMRT QA. This work demonstrates the feasibility of automatically computing delivered patient dose from portal dosimetry measurements in the Varian TPS system, which would provide a valuable and clinically viable IMRT QA tool for physicists and physicians. Methods: IMRT QA fluences were measured using portal dosimetry, processed using in-house matlab software, and imported back into Eclipse to calculate dose on the planning CT. To validate the proposed workflow, the Eclipse calculated portal dose for a 5-field sliding window prostate boost plan was processed as described above. The resulting dose was compared to the planned dose and found to be within 0.5 Gy. Two IMRT QA results for the prostate boost plan (one that failed and one that passed) were processed and the resulting patient doses were evaluated. Results: The max dose difference between IMRT QA #1 and the original planned and approved dose is 4.5 Gy, while the difference between the planned and IMRT QA #2 dose is 4.0 Gy. The inferior portion of the PTV is slightly underdosed in both plans, and the superior portion is slightly overdosed. The patient dose resulting from IMRT QA #1 and #2 differs by only 0.5 Gy. With this new information, it may be argued that the evaluated plan alteration to obtain passing gamma analysis produced clinically irrelevant differences. Conclusion: Evaluation of the delivered QA dose on the planning CT provides valuable information about the clinical relevance of failed or borderline IMRT QAs. This particular workflow demonstrates the feasibility of pushing the measured IMRT QA portal dosimetry results directly back onto the patient planning CT within the Varian system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-TeP4-07: Automatic EPID-Based 2D Measurement of MLC Leaf Offset as a Quality Control Tool

Purpose: The MLC dosimetric leaf gap (DLG) and transmission are measured parameters which impact the dosimetric accuracy of IMRT and VMAT plans. This investigation aims to develop an efficient and accurate routine constancy check of the physical DLG in two dimensions. Methods: The manufacturer’s recommended DLG measurement method was modified by using 5 fields instead of 11 and by utilizing the Electronic Portal Imaging Device (EPID). Validations were accomplished using an ion chamber (IC) in solid water and a 2D IC array. EPID data was collected for 6 months on multiple TrueBeam linacs using both Millennium and HD MLCs at 5 different clinics in an international consortium. Matlab code was written to automatically analyze the images and calculate the 2D results. Sensitivity was investigated by introducing deliberate leaf position errors. MLC calibration and initialization history was recorded to allow quantification of their impact. Results were analyzed using statistical process control (SPC). Results: The EPID method took approximately 5 minutes. Due to detector response, the EPID measured DLG and transmission differed from the IC values but were reproducible and consistent with changes measured using the ICs. For the Millennium MLC, the EPID measured DLG and transmission were both consistently lower than IC results. The EPID method was implemented as leaf offset and transmission constancy tests (LOC and TC). Based on 6 months of measurements, the initial leaf-specific action thresholds for changes from baseline were set to 0.1 mm. Upper and lower control limits for variation were developed for each machine. Conclusion: Leaf offset and transmission constancy tests were implemented on Varian HD and Millennium MLCs using an EPID and found to be efficient and accurate. The test is effective for monitoring MLC performance using dynamic delivery and performing process control on the DLG in 2D, thus enhancing dosimetric accuracy. This work was supported by a grant from Varian Medical Systems.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-G-TeP4-13: Interfraction Treatment Monitoring Using Integrated Invivo EPID Images

Purpose: To investigate inter-fraction differences of dose delivery by analyzing portal images acquired during treatment and implement an automated system to generate a report for each fraction. Large differences in images between fractions can alert the physicist of possible machine performance issues or patient set-up errors. Methods: A Varian Novalis Tx equipped with a HD120 MLC and aS1000 electronic portal imaging device (EPID) was used in our study. EPID images are acquired in continuous acquisition mode for 32 volumetric arc therapy (VMAT) patients. The images are summed to create an image for each arc and a single image for each fraction. The first fraction is designated as the reference unless a machine error prevented acquisition of all images. The images for each beam as well as the fraction image are compared using gamma analysis at 1%/1mm, 2%/2mm and 3%/3mm. A report is then generated using an in house MatLab program containing the comparison for the current fraction as well as a history of previous fractions. The reports are automatically sent via email to the physicist for review. Fractions in which the total number of images was not within 5% of the reference number of images were not included in the results. Results: 91 of the 182 fractions recorded an image count within 5% of the reference. Gamma averages over all fractions and patients were 96.2% ±0.8% at 3%/3mm, 92.9% ±1% at 2%/2mm and 80.6% ±1.8% at 1%/1mm. The SD between fractions for each patient ranged from .004% to 10.4%. Of the 91 fractions 3 flagged due to low gamma values. After further investigation no significant errors were found. Conclusion: This toolkit can be used for in-vivo monitoring of treatment plan delivery an alert the physics staff of any inter-fraction discrepancies that may require further investigation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-CAMPUS-TeP1-03: Pre-Treatment Surface Imaging Based Collision Detection

Purpose: Modern radiotherapy increasingly employs large immobilization devices, gantry attachments, and couch rotations for treatments. All of which raise the risk of collisions between the patient and the gantry / couch. Collision detection is often achieved by manually checking each couch position in the treatment room and sometimes results in extraneous imaging if collisions are detected after image based setup has begun. In the interest of improving efficiency and avoiding extra imaging, we explore the use of a surface imaging based collision detection model. Methods: Surfaces acquired from AlignRT (VisionRT, London, UK) were transferred in wavefront format to a custom Matlab (Mathworks, Natick, MA) software package (CCHECK). Computed tomography (CT) scans acquired at the same time were sent to CCHECK in DICOM format. In CCHECK, binary maps of the surfaces were created and overlaid on the CT images based on the fixed relationship of the AlignRT and CT coordinate systems. Isocenters were added through a graphical user interface (GUI). CCHECK then compares the inputted surfaces to a model of the linear accelerator (linac) to check for collisions at defined gantry and couch positions. Note, CCHECK may be used with or without a CT. Results: The nominal surface image field of view is 650 mm × 900 mm, with variance based on patient position and size. The accuracy of collision detections is primarily based on the linac model and the surface mapping process. The current linac model and mapping process yield detection accuracies on the order of 5 mm, assuming no change in patient posture between surface acquisition and treatment. Conclusions: CCHECK provides a non-ionizing method to check for collisions without the patient in the treatment room. Collision detection accuracy may be improved with more robust linac modeling. Additional gantry attachments (e.g. conical collimators) can be easily added to the model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-CAMPUS-TeP2-01: A Graph Form ADMM Algorithm for Constrained Quadratic Radiation Treatment Planning

Purpose: In radiation therapy optimization the constraints can be either hard constraints which must be satisfied or soft constraints which are included but do not need to be satisfied exactly. Currently the voxel dose constraints are viewed as soft constraints and included as a part of the objective function and approximated as an unconstrained problem. However in some treatment planning cases the constraints should be specified as hard constraints and solved by constrained optimization. The goal of this work is to present a computation efficiency graph form alternating direction method of multipliers (ADMM) algorithm for constrained quadratic treatment planning optimization and compare it with several commonly used algorithms/toolbox. Method: ADMM can be viewed as an attempt to blend the benefits of dual decomposition and augmented Lagrangian methods for constrained optimization. Various proximal operators were first constructed as applicable to quadratic IMRT constrained optimization and the problem was formulated in a graph form of ADMM. A pre-iteration operation for the projection of a point to a graph was also proposed to further accelerate the computation. Result: The graph form ADMM algorithm was tested by the Common Optimization for Radiation Therapy (CORT) dataset including TG119, prostate, liver, and head & neck cases. Both unconstrained and constrained optimization problems were formulated for comparison purposes. All optimizations were solved by LBFGS, IPOPT, Matlab built-in toolbox, CVX (implementing SeDuMi) and Mosek solvers. For unconstrained optimization, it was found that LBFGS performs the best, and it was 3–5 times faster than graph form ADMM. However, for constrained optimization, graph form ADMM was 8 – 100 times faster than the other solvers. Conclusion: A graph form ADMM can be applied to constrained quadratic IMRT optimization. It is more computationally efficient than several other commercial and noncommercial optimizers and it also used significantly less computer memory.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-AB-BRC-12: Optimized Parallel MonteCarlo Dose Calculations for Secondary MU Checks

Purpose: Secondary MU checks are an important tool used during a physics review of a treatment plan. Commercial software packages offer varying degrees of theoretical dose calculation accuracy, depending on the modality involved. Dose calculations of VMAT plans are especially prone to error due to the large approximations involved. Monte Carlo (MC) methods are not commonly used due to their long run times. We investigated two methods to increase the computational efficiency of MC dose simulations with the BEAMnrc code. Distributed computing resources, along with optimized code compilation, will allow for accurate and efficient VMAT dose calculations. Methods: The BEAMnrc package was installed on a high performance computing cluster accessible to our clinic. MATLAB and PYTHON scripts were developed to convert a clinical VMAT DICOM plan into BEAMnrc input files. The BEAMnrc installation was optimized by running the VMAT simulations through profiling tools which indicated the behavior of the constituent routines in the code, e.g. the bremsstrahlung splitting routine, and the specified random number generator. This information aided in determining the most efficient compiling parallel configuration for the specific CPU’s available on our cluster, resulting in the fastest VMAT simulation times. Our method was evaluated with calculations involving 10{sup 8} – 10{sup 9} particle histories which are sufficient to verify patient dose using VMAT. Results: Parallelization allowed the calculation of patient dose on the order of 10 – 15 hours with 100 parallel jobs. Due to the compiler optimization process, further speed increases of 23% were achieved when compared with the open-source compiler BEAMnrc packages. Conclusion: Analysis of the BEAMnrc code allowed us to optimize the compiler configuration for VMAT dose calculations. In future work, the optimized MC code, in conjunction with the parallel processing capabilities of BEAMnrc, will be applied to provide accurate and efficient secondary MU checks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-D-209-03: Alignment of the Patient Graphic Model Using Fluoroscopic Images for Skin Dose Mapping

Purpose: The Dose Tracking System (DTS) was developed to provide realtime feedback of skin dose and dose rate during interventional fluoroscopic procedures. A color map on a 3D graphic of the patient represents the cumulative dose distribution on the skin. Automated image correlation algorithms are described which use the fluoroscopic procedure images to align and scale the patient graphic for more accurate dose mapping. Methods: Currently, the DTS employs manual patient graphic selection and alignment. To improve the accuracy of dose mapping and automate the software, various methods are explored to extract information about the beam location and patient morphology from the procedure images. To match patient anatomy with a reference projection image, preprocessing is first used, including edge enhancement, edge detection, and contour detection. Template matching algorithms from OpenCV are then employed to find the location of the beam. Once a match is found, the reference graphic is scaled and rotated to fit the patient, using image registration correlation functions in Matlab. The algorithm runs correlation functions for all points and maps all correlation confidences to a surface map. The highest point of correlation is used for alignment and scaling. The transformation data is saved for later model scaling. Results: Anatomic recognition is used to find matching features between model and image and image registration correlation provides for alignment and scaling at any rotation angle with less than onesecond runtime, and at noise levels in excess of 150% of those found in normal procedures. Conclusion: The algorithm provides the necessary scaling and alignment tools to improve the accuracy of dose distribution mapping on the patient graphic with the DTS. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-D-209-05: Automatic Calculation of Organ and Effective Dose for CBCT and Interventional Fluoroscopic Procedures

Purpose: To compare PCXMC and EGSnrc calculated organ and effective radiation doses from cone-beam computed tomography (CBCT) and interventional fluoroscopically-guided procedures using automatic exposure-event grouping. Methods: For CBCT, we used PCXMC20Rotation.exe to automatically calculate the doses and compared the results to those calculated using EGSnrc with the Zubal patient phantom. For interventional procedures, we use the dose tracking system (DTS) which we previously developed to produce a log file of all geometry and exposure parameters for every x-ray pulse during a procedure, and the data in the log file is input into PCXMC and EGSnrc for dose calculation. A MATLAB program reads data from the log files and groups similar exposures to reduce calculation time. The definition files are then automatically generated in the format used by PCXMC and EGSnrc. Processing is done at the end of the procedure after all exposures are completed. Results: For the Toshiba Infinix CBCT LCI-Middle-Abdominal protocol, most organ doses calculated with PCXMC20Rotation closely matched those calculated with EGSnrc. The effective doses were 33.77 mSv with PCXMC20Rotation and 32.46 mSv with EGSnrc. For a simulated interventional cardiac procedure, similar close agreement in organ dose was obtained between the two codes; the effective doses were 12.02 mSv with PCXMC and 11.35 mSv with EGSnrc. The calculations can be completed on a PC without manual intervention in less than 15 minutes with PCXMC and in about 10 hours with EGSnrc, depending on the level of data grouping and accuracy desired. Conclusion: Effective dose and most organ doses in CBCT and interventional radiology calculated by PCXMC closely match those calculated by EGSnrc. Data grouping, which can be done automatically, makes the calculation time with PCXMC on a standard PC acceptable. This capability expands the dose information that can be provided by the DTS. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-201-05: Varian MPC as a Statistical Process Control Tool

Purpose: Quality assurance in radiotherapy requires the measurement of various machine parameters to ensure they remain within permitted values over time. In Truebeam release 2.0 the Machine Performance Check (MPC) was released allowing beam output and machine axis movements to be assessed in a single test. We aim to evaluate the Varian Machine Performance Check (MPC) as a tool for Statistical Process Control (SPC). Methods: Varian’s MPC tool was used on three Truebeam and one EDGE linac for a period of approximately one year. MPC was commissioned against independent systems. After this period the data were reviewed to determine whether or not the MPC was useful as a process control tool. Analyses on individual tests were analysed using Shewhart control plots, using Matlab for analysis. Principal component analysis was used to determine if a multivariate model was of any benefit in analysing the data. Results: Control charts were found to be useful to detect beam output changes, worn T-nuts and jaw calibration issues. Upper and lower control limits were defined at the 95% level. Multivariate SPC was performed using Principal Component Analysis. We found little evidence of clustering beyond that which might be naively expected such as beam uniformity and beam output. Whilst this makes multivariate analysis of little use it suggests that each test is giving independent information. Conclusion: The variety of independent parameters tested in MPC makes it a sensitive tool for routine machine QA. We have determined that using control charts in our QA programme would rapidly detect changes in machine performance. The use of control charts allows large quantities of tests to be performed on all linacs without visual inspection of all results. The use of control limits alerts users when data are inconsistent with previous measurements before they become out of specification. A. Carver has received a speaker’s honorarium from Varian.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

TU-FG-201-07: Development of SRS Conical Collimator Collision Prediction Software for Radiation Treatment Safety

Purpose: Stereotactic radiosurgery (SRS) treatments using conical collimators can potentially result in gantry collision with treatment table due to limited collision-clear spaces. An in-house software was developed to help the SRS treatment planner mitigate potential SRS conical collimator (Varian Medical System, Palo Alto, CA) collisions with the treatment table. This software was designed to remove treatment re-planning secondary to unexpected collisions. Methods: A BrainLAB SRS ICT Frameless Extension used for SRS treatments in our clinic was mathematically modelled using surface points registered to the 3D co-ordinate space of the couch extension. The surface points are transformed based on the treatment isocenter point and potential collisions are determined in 3D space for couch and gantry angle combinations. The distance between the SRS conical collimators and LINAC isocenter is known. The collision detection model was programmed in MATLAB (Mathwork, Natick, MA) to display graphical plots of the calculations, and the plotted data is used to avoid the gantry and couch angle combinations that would likely result in a collision. We have utilized the cone collision tool for 23 SRS cone treatment plans (8 retrospective and 15 prospective for 10 patients). Results: Twenty one plans strongly agreed with the software tool prediction for collision. However, in two plans, a collision was observed with a 0.5 cm margin when the software predicted no collision. Therefore, additional margins were added to the clearance criteria in the program to achieve a lower risk of actual collisions. Conclusion: Our in-house developed collision check software successfully avoided SRS cone re-planning by 91.3% due to a reduction in cone collisions with the treatment table. Future developments to our software will include a CT image data set based collision prediction model as well as a beam angle optimization tool to avoid normal critical tissues as well as previously treated lesions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-209-04: Testing of Digital Image Receptors Using AAPM TG-150’s Draft Recommendations - Investigating the Impact of Different Processing Parameters

Purpose: To evaluate implementation of AAPM TG-150’s draft recommendations via a parameter study for testing the performance of digital image receptors. Methods: Flat field images were acquired from 9 calibrated digital image receptors associated with 9 new portable digital radiography systems (Carestream Health, Inc.) based on the draft recommendations and manufacturer-specified calibration conditions (set of 4 images at input detector air kerma ranging from 1 to 25 µGy). Effects of exposure response function (linearized and logarithmic), ‘Presentation Intent Type’ (‘For Processing’ and ‘For Presentation’), detector orientation with respect to the anode-cathode axis (4 orientations; 900 rotations per iteration), different ROI sizes (5×5–40×40 mm{sup 2}) and elimination of varying dimensions of image border (0 mm i.e., without boundary elimination to 150 mm) on signal, noise, signal-to-noise ratio (SNR) and the associated nonuniformities were evaluated. Images were analyzed in Matlab and quantities were compared using ANOVA. Results: Signal, noise and SNR values averaged over 9 systems with default parameter values in draft recommendations were 4837.2±139.4, 19.7±0.9 and 246.4±10.1 (mean ± standard deviation), respectively (at input detector air kerma: 12.5 µGy). Signal, noise and SNR showed characteristic dependency on exposure response function and on ‘Presentation Intent Type’. These values were not affected by ROI size and detector orientation, but analysis showed that eliminating the edge pixels along the boundary was required for the noise parameter (coefficient of variation range for noise: 72%–106% and 3%–4% without and with boundary elimination; respectively). Local and global nonuniformities showed a similar dependence on the need for boundary elimination. Interestingly, computed non-uniformities showed agreement with manufacturer-reported values except for noise non-uniformities in two units; artifacts were seen in images from these two units highlighting the importance of independent evaluations. Conclusion: The effect of different parameters on performance characterization of digital image receptors was evaluated based on TG-150’s draft recommendations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-BRA-04: A Novel Superconducting Magnet Design for Optimized Patient Access and Minimal SSD for Use in a Linac-MR Hybrid

Purpose: A prototype rotating hybrid MR imaging system and linac has been developed to allow for simultaneous imaging and radiation delivery parallel to B{sub 0}. However, the design of a compact magnet capable of rotation in a small vault with sufficient patient access and a typical clinical source-to-surface distance (SSD) is challenging. This work presents a novel superconducting magnet design that allows for a reduced SSD and ample patient access by moving the superconducting coils to the side of the yoke. The yoke and pole-plate structures are shaped to direct the magnetic flux appropriately. Methods: The surface of the pole plate for the magnet assembly is optimized. The magnetic field calculations required in this work are performed with the 3D finite element method software package Opera-3D. Each tentative design strategy is virtually modeled in this software package and externally controlled by MATLAB, with its key geometries defined as variables. The particle swarm optimization algorithm is used to optimize the variables subject to the minimization of a cost function. At each iteration, Opera-3D will solve the magnetic field solution over a field-of-view suitable for MR imaging and the degree of field uniformity will be assessed to calculate the value of the cost function associated with that iteration. Results: An optimized magnet assembly that generates a homogenous 0.2T magnetic field over an ellipsoid with large axis of 30 cm and small axes of 20 cm is obtained. Conclusion: The distinct features of this model are the minimal distance between the yoke’s top and the isocentre and the improved patient access. On the other hand, having homogeneity over an ellipsoid give us a larger field-of-view, essential for geometric accuracy of the MRI system. The increase of B{sub 0} from 0.2T in the present model to 0.5T is the subject of future work. Funding Sources: Alberta Innovates - Health Solutions (AIHS)| Disclosure and Conflict of Interest: B. Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta biplanar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-206-04: An Effective Homomorphic Unsharp Mask Filtering Method to Correct Intensity Inhomogeneity in Daily Treatment MR Images

Purpose: The daily treatment MRIs acquired on MR-IGRT systems, like diagnostic MRIs, suffer from intensity inhomogeneity issue, associated with B1 and B0 inhomogeneities. An improved homomorphic unsharp mask (HUM) filtering method, automatic and robust body segmentation, and imaging field-of-view (FOV) detection methods were developed to compute the multiplicative slow-varying correction field and correct the intensity inhomogeneity. The goal is to improve and normalize the voxel intensity so that the images could be processed more accurately by quantitative methods (e.g., segmentation and registration) that require consistent image voxel intensity values. Methods: HUM methods have been widely used for years. A body mask is required, otherwise the body surface in the corrected image would be incorrectly bright due to the sudden intensity transition at the body surface. In this study, we developed an improved HUM-based correction method that includes three main components: 1) Robust body segmentation on the normalized image gradient map, 2) Robust FOV detection (needed for body segmentation) using region growing and morphologic filters, and 3) An effective implementation of HUM using repeated Gaussian convolution. Results: The proposed method was successfully tested on patient images of common anatomical sites (H/N, lung, abdomen and pelvis). Initial qualitative comparisons showed that this improved HUM method outperformed three recently published algorithms (FCM, LEMS, MICO) in both computation speed (by 50+ times) and robustness (in intermediate to severe inhomogeneity situations). Currently implemented in MATLAB, it takes 20 to 25 seconds to process a 3D MRI volume. Conclusion: Compared to more sophisticated MRI inhomogeneity correction algorithms, the improved HUM method is simple and effective. The inhomogeneity correction, body mask, and FOV detection methods developed in this study would be useful as preprocessing tools for many MRI-related research and clinical applications in radiotherapy. Authors have received research grants from ViewRay and Varian.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-FG-207B-06: Plaque Composition Measurement with Dual Energy Computed Tomography

Purpose: To investigate the feasibility of characterizing arterial plaque composition in terms of water, lipid and protein or calcium using dual energy computed tomography. Characterization of plaque composition can potentially help distinguish vulnerable from stable plaques. Methods: Simulations studies were performed by the CT simulator based on ASTRA tomography toolbox. The beam energy for dual energy images was selected to be 80 kVp and 135 kVp. The radiation dose and energy spectrum for the CT simulator were carefully calibrated with respect to a 320-slice CT scanner. A digital chest phantom was constructed using Matlab for calibration and plaque measurement. Pure water, lipid, protein or calcium was used for calibration and a mixture of different volume percentages of these materials were used for validation purposes. Non-calcified plaque was simulated using water, lipid and protein with volumetric percentage range of 35%∼65%, 5%∼60% and 5%∼40%, respectively. Calcified plaque was simulated using water, lipid and calcium with volumetric percentage range of 50%∼80%, 8%∼45% and 3%∼13%, respectively. We employed iterative sinogram processing (ISP) to reduce the beam hardening effect in the simulation to improve the decomposition results. Results: The simulated known composition and dual energy decomposition results were in good agreement. Water, lipid and protein (calcium) mixtures were decomposed into water, lipid and protein (calcium) contents. The RMS errors of volumetric percentage for the water, lipid and protein (non-calcified plaque) decomposition, as compared to known values, were estimated to be approximately 5.74%, 2.54%, and 0.95% respectively. The RMS errors of volumetric percentage for the water, lipid and Calcium (calcified plaque) decomposition, as compared to known values, were estimated to be approximately 7.4%, 8.64%, and 0.08% respectively. Conclusion: The results of this study suggest that the dual energy decomposition can potentially be used to quantify the water, lipid, and protein or calcium composition of a plaque with relatively good accuracy. Grant funding from Toshiba Medical Systems and Philips Medical Systems.

60 APPLIED LIFE SCIENCES↗