Search NASA⌕ Search

Engineering topics

Zhou, F

Publications and source records attributed to Zhou, F.

Fast neural network training on a cluster of GPUs for action recognition with high accuracy

In this work, we propose algorithms and techniques to accelerate training of deep neural networks for action recognition on a cluster of GPUs. The convergence analysis of our algorithm shows it is possible to reduce communication cost and at the same time minimize the number of iterations needed for convergence. We customize the Adam optimizer for our distributed algorithm to improve efficiency. In addition, we employ transfer-learning to further reduce training time while improving validation accuracy. For the UCF101 and HMDB51 datasets, the validation accuracies achieved are 93.1% and 67.9% respectively. With an additional end-to-end trained temporal stream, the validation accuracies achieved for UCF101 and HMDB51 are 93.47% and 81.24% respectively. As far as we know, these are the highest accuracies achieved with the two-stream approach using ResNet that does not involve computationally expensive 3D convolutions or pretraining on much larger datasets.

97 MATHEMATICS AND COMPUTING↗

High Thermoelectric Power Factor in Intermetallic Co Si Arising from Energy Filtering of Electrons by Phonon Scattering

Intermetallic compounds are usually poor thermoelectrics due to the high electronic densities of states at the Fermi level and concomitantly low Seebeck coefficients. However, intermetallic B20-type CoSi has been experimentally shown to attain remarkably large negative values of the Seebeck coefficient. We provide a theoretical explanation for this surprising fact using state-of-the-art first-principles calculations with explicit treatment of electron-phonon interactions. We find that the unique band structure of CoSi, which has both massless and heavy fermion bands near the Fermi level, facilitates strong scattering of the low-energy electrons by phonons into the heavy band, resulting in effective energy filtering and high thermal voltage. Our study predicts that a very large thermoelectric power factor of 80 μ W cm -1 K -2 or higher is experimentally accessible in the 300–600 K range and highlights a general principle for identifying intermetallic compounds with large thermopower.

36 MATERIALS SCIENCE↗

SU-F-J-98: Improvement and Evaluation Of Deformation Image Registration On Parotid Glands During Radiation Therapy for Nasopharyngeal Cancer

Purpose: To quantitatively evaluate the strategic innovation and accuracy variation of deformation registration algorithm for parotid glands using the similarity Dice coefficient during the course of radiation therapy (RT) for nasopharyngeal cancer (NPC). Methods: Daily MVCT data for 10 patients with pathologically proven nasopharyngeal cancers were analyzed. The data were acquired using tomotherapy (TomoTherapy, Accuray) at the PLA General Hospital. The prescription dose to the primary target was 70Gy in 33 fractions. Two kinds of contours for parotid glands on daily MVCTs were obtained by populating these contours from planning CTs to the daily CTs via rigid-body registration with or without the rotation shifts using the in-house tools and the Adaptive plan software (Adaptive Plan, TomoTherapy), and were edited manually if necessary. The diffeomorphic Demons algorithm developed in the in-house tool was used to propagate the parotid structures from the daily CTs to planning CTs. The differences of the mapped parotid contours in two methods were evaluated using Dice similarity index (DSI). Two-tailed t-test analysis was carried out to compare the DSI changes during the course of RT. Results: For 10 patient plans, the accuracy of deformation image registration (DIR) with the rotation shift was obviously better than those without the rotation shift. The Dice scores of the ipsi- and contra-lateral parotids for with and without the rotation shifts were found to be correlated with each other [0.904±0.031 vs 0.919±0.030 (p<0.001); 0.900±0.031 vs 0.910±0.032 (p<0.001)]. The Dice scores for the parotids have shown the reduction with the changes of parotid volumes during RT. The DSI values between the first and last fraction were 0.932±0.020 vs 0.899±0.030 in 10 patient plans. Conclusion: DIR was successfully improved using the strategic innovation for ART. And the decrease of DIR accuracy has also been found during the delivery of fractionated radiotherapy. This work was supported in part by the grant from Chinese Natural Science Foundation (Grant No. 11105225).

60 APPLIED LIFE SCIENCES↗

SU-F-T-302: Implementation of the RadCalc Image Analysis Tool for IMRT QA

Purpose: To assess RadCalc as a means for analyzing fluence maps of IMRT QA patient plans. Eclipse generated fluence maps were compared to RadCalc fluence maps generated from both Eclipse parameters and from dynalog trajectory files. Methods: Six IMRT plans consisting of fifty individual fields were compared both field-by-field and as composite plans. These plans were exported from Eclipse to RadCalc. Each plan was then delivered in QA mode on a Varian Clinac iX, while recording and saving the dynalog trajectory files. These files were then imported into RadCalc, providing three sources of generating fluence maps: Eclipse, RadCalc, and dynalog files. Using the image analysis tool in RadCalc, a gamma analysis was then performed as a point-by-point comparison for each IMRT field. The Eclipse fluence map was used as the baseline for each comparison. Results: Based on the manufacturer’s recommendations, all fields were normalized to the maximum pixel value within each Eclipse field. All data was analyzed using a gamma index of 3mm/3% with passing criteria of 90%. For both dynalog analysis and RadCalc fluence analysis, 48/50 created fluence maps passed at greater than 90% when compared with the Eclipse baseline. In analyzing the composite of each of the six patients’ plans, all plans passed over 90%. Conclusion: Using the dynalog trajectory files in combination with RadCalc Version 6.3 image analysis is a promising metric for verifying IMRT QA passing rates. Notably, if a field failed, it did so for both dynalog and RadCalc compared to Eclipse. This suggests RadCalc can accurately simulate fluence maps, with similar results to dynalog comparisons. Limitations within the RadCalc software were discovered in the analytical process such as pixel resolution and the inability to set minimum threshold values. Comparisons could be extended to include dose map and distance to agreement analysis by expanding software capabilities.

60 APPLIED LIFE SCIENCES↗

SU-F-T-347: An Absolute Dose-Volume Constraint Based Deterministic Optimization Framework for Multi-Co60 Source Focused Radiotherapy

Purpose: Treatment plan optimization in multi-Co60 source focused radiotherapy with multiple isocenters is challenging, because dose distribution is normalized to maximum dose during optimization and evaluation. The objective functions are traditionally defined based on relative dosimetric distribution. This study presents an alternative absolute dose-volume constraint (ADC) based deterministic optimization framework (ADC-DOF). Methods: The initial isocenters are placed on the eroded target surface. Collimator size is chosen based on the area of 2D contour on corresponding axial slice. The isocenter spacing is determined by adjacent collimator sizes. The weights are optimized by minimizing the deviation from ADCs using the steepest descent technique. An iterative procedure is developed to reduce the number of isocenters, where the isocenter with lowest weight is removed without affecting plan quality. The ADC-DOF is compared with the genetic algorithm (GA) using the same arbitrary shaped target (254cc), with a 15mm margin ring structure representing normal tissues. Results: For ADC-DOF, the ADCs imposed on target and ring are (D100>10Gy, D50,10, 0<12Gy, 15Gy and 20Gy) and (D40<10Gy). The resulting D100, 50, 10, 0 and D40 are (9.9Gy, 12.0Gy, 14.1Gy and 16.2Gy) and (10.2Gy). The objectives of GA are to maximize 50% isodose target coverage (TC) while minimize the dose delivered to the ring structure, which results in 97% TC and 47.2% average dose in ring structure. For ADC-DOF (GA) techniques, 20 out of 38 (10 out of 12) initial isocenters are used in the final plan, and the computation time is 8.7s (412.2s) on an i5 computer. Conclusion: We have developed a new optimization technique using ADC and deterministic optimization. Compared with GA, ADC-DOF uses more isocenters but is faster and more robust, and achieves a better conformity. For future work, we will focus on developing a more effective mechanism for initial isocenter determination.

60 APPLIED LIFE SCIENCES↗

SU-F-T-620: Development of a Convolution/Superposition Dose Engine for CyberKnife System

Purpose: Current CyberKnife treatment planning system (TPS) provided two dose calculation algorithms: Ray-tracing and Monte Carlo. Ray-tracing algorithm is fast, but less accurate, and also can’t handle irregular fields since a multi-leaf collimator system was recently introduced to CyberKnife M6 system. Monte Carlo method has well-known accuracy, but the current version still takes a long time to finish dose calculations. The purpose of this paper is to develop a GPU-based fast C/S dose engine for CyberKnife system to achieve both accuracy and efficiency. Methods: The TERMA distribution from a poly-energetic source was calculated based on beam’s eye view coordinate system, which is GPU friendly and has linear complexity. The dose distribution was then computed by inversely collecting the energy depositions from all TERMA points along 192 collapsed-cone directions. EGSnrc user code was used to pre-calculate energy deposition kernels (EDKs) for a series of mono-energy photons The energy spectrum was reconstructed based on measured tissue maximum ratio (TMR) curve, the TERMA averaged cumulative kernels was then calculated. Beam hardening parameters and intensity profiles were optimized based on measurement data from CyberKnife system. Results: The difference between measured and calculated TMR are less than 1% for all collimators except in the build-up regions. The calculated profiles also showed good agreements with the measured doses within 1% except in the penumbra regions. The developed C/S dose engine was also used to evaluate four clinical CyberKnife treatment plans, the results showed a better dose calculation accuracy than Ray-tracing algorithm compared with Monte Carlo method for heterogeneous cases. For the dose calculation time, it takes about several seconds for one beam depends on collimator size and dose calculation grids. Conclusion: A GPU-based C/S dose engine has been developed for CyberKnife system, which was proven to be efficient and accurate for clinical purpose, and can be easily implemented in TPS.

60 APPLIED LIFE SCIENCES↗

SU-G-TeP1-08: LINAC Head Geometry Modeling for Cyber Knife System

Purpose: Knowledge of the LINAC head information is critical for model based dose calculation algorithms. However, the geometries are difficult to measure precisely. The purpose of this study is to develop linac head models for Cyber Knife system (CKS). Methods: For CKS, the commissioning data were measured in water at 800mm SAD. The measured full width at half maximum (FWHM) for each cone was found greater than the nominal value, this was further confirmed by additional film measurement in air. Diameter correction, cone shift and source shift models (DCM, CSM and SSM) are proposed to account for the differences. In DCM, a cone-specific correction is applied. For CSM and SSM, a single shift is applied to the cone or source physical position. All three models were validated with an in-house developed pencil beam dose calculation algorithm, and further evaluated by the collimator scatter factor (Sc) correction. Results: The mean square error (MSE) between nominal diameter and the FWHM derived from commissioning data and in-air measurement are 0.54mm and 0.44mm, with the discrepancy increasing with cone size. Optimal shift for CSM and SSM is found to be 9mm upward and 18mm downward, respectively. The MSE in FWHM is reduced to 0.04mm and 0.14mm for DCM and CSM (SSM). Both DCM and CSM result in the same set of Sc values. Combining all cones at SAD 600–1000mm, the average deviation from 1 in Sc of DCM (CSM) and SSM is 2.6% and 2.2%, and reduced to 0.9% and 0.7% for the cones with diameter greater than 15mm. Conclusion: We developed three geometrical models for CKS. All models can handle the discrepancy between vendor specifications and commissioning data. And SSM has the best performance for Sc correction. The study also validated that a point source can be used in CKS dose calculation algorithms.

60 APPLIED LIFE SCIENCES↗

SU-E-T-181: Clinical Implementation of Task Group 176

Purpose: The implementation of TG-176 with regards to immobilization devices and couch tops as they effect dosimetric treatment planning. Methods: The external devices used clinically were scanned to measure their HU values. Plans were created in the Eclipse treatment planning system (TPS) using these devices, one that accounted for the correct HU value of the each device and another that did not include the device as a structure. A dose subtraction was performed between the two plans to evaluate the dosimetric differences. The metrics used for evaluation included attenuation and surface dose. Plan parameters were varied to evaluate the impact of the devices in different clinical scenarios. Results: While the exact HU values of our results are clinic-dependent, the protocol being implemented is widely applicable. We recommend a four step process for implementing this task group. First, physics should scan each treatment device to determine accurate HU values. Second, CT therapists should include in the setup note which table top was used during patient CT simulation and are asked to make immobilization devices as uniform in thickness as possible. Therapists should also index the devices whenever possible so beam will traverse the same area of the device. Third, the dosimetrist should manually correct the HU value for any external device, including the couch. For H&N cases, the rails must be removed from the couch structure. When rails are used during treatments, it is important to make note of their exact position in the setup notes. Finally, physicians should be made aware that there could be changes in surface doses depending on whether or not immobilization devices or couch tops are in the beam path. Conclusion: The protocol outlined above was implemented to reduce the errors that arise from ignoring effects of external devices, thus ensuring safer, more accurate patient treatments.

60 APPLIED LIFE SCIENCES↗

SU-E-T-252: Developing a Pencil Beam Dose Calculation Algorithm for CyberKnife System

Purpose: Currently there are two dose calculation algorithms available in the Cyberknife planning system: ray-tracing and Monte Carlo, which is either not accurate or time-consuming for irregular field shaped by the MLC that was recently introduced. The purpose of this study is to develop a fast and accurate pencil beam dose calculation algorithm which can handle irregular field. Methods: A pencil beam dose calculation algorithm widely used in Linac system is modified. The algorithm models both primary (short range) and scatter (long range) components with a single input parameter: TPR{sub 20}/{sub 10}. The TPR{sub 20}/{sub 20}/{sub 10} value was first estimated to derive an initial set of pencil beam model parameters (PBMP). The agreement between predicted and measured TPRs for all cones were evaluated using the root mean square of the difference (RMSTPR), which was then minimized by adjusting PBMPs. PBMPs are further tuned to minimize OCR RMS (RMSocr) by focusing at the outfield region. Finally, an arbitrary intensity profile is optimized by minimizing RMSocr difference at infield region. To test model validity, the PBMPs were obtained by fitting to only a subset of cones (4) and applied to all cones (12) for evaluation. Results: With RMS values normalized to the dmax and all cones combined, the average RMSTPR at build-up and descending region is 2.3% and 0.4%, respectively. The RMSocr at infield, penumbra and outfield region is 1.5%, 7.8% and 0.6%, respectively. Average DTA in penumbra region is 0.5mm. There is no trend found in TPR or OCR agreement among cones or depths. Conclusion: We have developed a pencil beam algorithm for Cyberknife system. The prediction agrees well with commissioning data. Only a subset of measurements is needed to derive the model. Further improvements are needed for TPR buildup region and OCR penumbra. Experimental validations on MLC shaped irregular field needs to be performed. This work was partially supported by the National Natural Science Foundation of China (61171005) and the China Scholarship Council (CSC)

60 APPLIED LIFE SCIENCES↗