Search NASA⌕ Search

Engineering topics

Gu, X

Publications and source records attributed to Gu, X.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-J-12: An Image-Guided Soft Robotic Patient Positioning System for Maskless Head-And-Neck Cancer Radiotherapy: A Proof-Of-Concept Study

Purpose: We propose a surface-image-guided soft robotic patient positioning system for maskless head-and-neck radiotherapy. The ultimate goal of this project is to utilize a soft robot to realize non-rigid patient positioning and real-time motion compensation. In this proof-of-concept study, we design a position-based visual servoing control system for an air-bladder-based soft robot and investigate its performance in controlling the flexion/extension cranial motion on a mannequin head phantom. Methods: The current system consists of Microsoft Kinect depth camera, an inflatable air bladder (IAB), pressured air source, pneumatic valve actuators, custom-built current regulators, and a National Instruments myRIO microcontroller. The performance of the designed system was evaluated on a mannequin head, with a ball joint fixed below its neck to simulate torso-induced head motion along flexion/extension direction. The IAB is placed beneath the mannequin head. The Kinect camera captures images of the mannequin head, extracts the face, and measures the position of the head relative to the camera. This distance is sent to the myRIO, which runs control algorithms and sends actuation commands to the valves, inflating and deflating the IAB to induce head motion. Results: For a step input, i.e. regulation of the head to a constant displacement, the maximum error was a 6% overshoot, which the system then reduces to 0% steady-state error. In this initial investigation, the settling time to reach the regulated position was approximately 8 seconds, with 2 seconds of delay between the command start of motion due to capacitance of the pneumatics, for a total of 10 seconds to regulate the error. Conclusion: The surface image-guided soft robotic patient positioning system can achieve accurate mannequin head flexion/extension motion. Given this promising initial Result, the extension of the current one-dimensional soft robot control to multiple IABs for non-rigid positioning control will be pursued.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-06: 4D Particle Swarm Optimization to Enable Lung SBRT in Patients with Central And/or Large Tumors

Purpose: Patients presenting with large and/or centrally-located lung tumors are currently considered ineligible for highly potent regimens such as SBRT due to concerns of toxicity to normal tissues and organs-at-risk (OARs). We present a particle swarm optimization (PSO)-based 4D planning technique, designed for MLC tracking delivery, that exploits the temporal dimension as an additional degree of freedom to significantly improve OAR-sparing and reduce toxicity to levels clinically considered as acceptable for SBRT administration. Methods: Two early-stage SBRT-ineligible NSCLC patients were considered, presenting with tumors of maximum dimensions of 7.4cm (central-right lobe; 1.5cm motion) and 11.9cm (upper-right lobe; 1cm motion). In each case, the target and normal structures were manually contoured on each of the ten 4DCT phases. Corresponding ten initial 3D-conformal plans (Pt#1: 7-beams; Pt#2: 9-beams) were generated using the Eclipse planning system. Using 4D-PSO, fluence weights were optimized over all beams and all phases (70 and 90 apertures for Pt1&2, respectively). Doses to normal tissues and OARs were compared with clinicallyestablished lung SBRT guidelines based on RTOG-0236. Results: The PSO-based 4D SBRT plan yielded tumor coverage and dose—sparing for parallel and serial OARs within the SBRT guidelines for both patients. The dose-sparing compared to the clinically-delivered conventionallyfractionated plan for Patient 1 (Patient 2) was: heart Dmean = 11% (33%); lung V20 = 16% (21%); lung Dmean = 7% (20%); spinal cord Dmax = 5% (16%); spinal cord Dmean = 7% (33%); esophagus Dmax = 0% (18%). Conclusion: Truly 4D planning can significantly reduce dose to normal tissues and OARs. Such sparing opens up the possibility of using highly potent and effective regimens such as lung SBRT for patients who were conventionally considered SBRT non-eligible. Given the large, non-convex solution space, PSO represents an attractive, parallelizable tool to successfully achieve a globally optimal solution for this problem. This work was supported through funding from the National Institutes of Health and Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-175: Clinical Evaluations of Monte Carlo-Based Inverse Treatment Plan Optimization for Intensity Modulated Radiotherapy

Purpose: Pencil-beam or superposition-convolution type dose calculation algorithms are routinely used in inverse plan optimization for intensity modulated radiation therapy (IMRT). However, due to their limited accuracy in some challenging cases, e.g. lung, the resulting dose may lose its optimality after being recomputed using an accurate algorithm, e.g. Monte Carlo (MC). It is the objective of this study to evaluate the feasibility and advantages of a new method to include MC in the treatment planning process. Methods: We developed a scheme to iteratively perform MC-based beamlet dose calculations and plan optimization. In the MC stage, a GPU-based dose engine was used and the particle number sampled from a beamlet was proportional to its optimized fluence from the previous step. We tested this scheme in four lung cancer IMRT cases. For each case, the original plan dose, plan dose re-computed by MC, and dose optimized by our scheme were obtained. Clinically relevant dosimetric quantities in these three plans were compared. Results: Although the original plan achieved a satisfactory PDV dose coverage, after re-computing doses using MC method, it was found that the PTV D95% were reduced by 4.60%–6.67%. After re-optimizing these cases with our scheme, the PTV coverage was improved to the same level as in the original plan, while the critical OAR coverages were maintained to clinically acceptable levels. Regarding the computation time, it took on average 144 sec per case using only one GPU card, including both MC-based beamlet dose calculation and treatment plan optimization. Conclusion: The achieved dosimetric gains and high computational efficiency indicate the feasibility and advantages of the proposed MC-based IMRT optimization method. Comprehensive validations in more patient cases are in progress.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-184: Clinical VMAT QA Practice Using LINAC Delivery Log Files

Purpose: To evaluate the accuracy of volumetric modulated arc therapy (VMAT) treatment delivery dose clouds by comparing linac log data to doses measured using an ionization chamber and film. Methods: A commercial IMRT quality assurance (QA) process utilizing a DICOM-RT framework was tested for clinical practice using 30 prostate and 30 head and neck VMAT plans. Delivered 3D VMAT dose distributions were independently checked using a PinPoint ionization chamber and radiographic film in a solid water phantom. DICOM RT coordinates were used to extract the corresponding point and planar doses from 3D log file dose distributions. Point doses were evaluated by computing the percent error between log file and chamber measured values. A planar dose evaluation was performed for each plan using a 2D gamma analysis with 3% global dose difference and 3 mm isodose point distance criteria. The same analysis was performed to compare treatment planning system (TPS) doses to measured values to establish a baseline assessment of agreement. Results: The mean percent error between log file and ionization chamber dose was 1.0%±2.1% for prostate VMAT plans and −0.2%±1.4% for head and neck plans. The corresponding TPS calculated and measured ionization chamber values agree within 1.7%±1.6%. The average 2D gamma passing rates for the log file comparison to film are 98.8%±1.0% and 96.2%±4.2% for the prostate and head and neck plans, respectively. The corresponding passing rates for the TPS comparison to film are 99.4%±0.5% and 93.9%±5.1%. Overall, the point dose and film data indicate that log file determined doses are in excellent agreement with measured values. Conclusion: Clinical VMAT QA practice using LINAC treatment log files is a fast and reliable method for patient-specific plan evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-217: Comprehensive Dosimetric Evaluation On 3D-CRT, IMRT and Non-Coplanar Arc Treatment for Prone Accelerated Partial Breast Irradiation (APBI)

Purpose: Accelerated partial breast irradiation (APBI) is an effective treatment for early stage breast-cancer. Irradiation in a prone position can mitigate breast motion and spare heart and lung. In this study, a comprehensive study is performed to evaluate various treatment techniques for prone APBI treatment including: 3D-CRT, IMRT, co-planar and non-coplanar partial arcs treatment. Methods: In this treatment planning study, a left breast patient treated in prone position in our clinic was imported into Varian Eclipse TPS. Six beams tangential to chest wall were used in both 3D-CRT and IMRT plans. These six beams were coplanar in a transactional plane achieved by both gantry and couch rotation. A 60-beam IMRT plan was also created to explore the maximum benefit of co-planar IMRT. Within deliverable couch rotation range (±30°), partial arc treatment plans with one and up to ten couch positions were generated for comparison. For each plan, 30Gy in 6 fractions was prescribed to 95% PTV volume. Critical dosimetric parameters, such as conformity index, mean, maximum, and volume dose of organ at risk, are evaluated. Results: The conformity indexes (CI) are 3.53, 3.17, 2.21 and 1.08 respectively to 3D-CRT, 6-beam IMRT, 60-beam IMRT, and two-partial-arcs coplanar plans. However, arc plans increase heart dose. CI for non-coplanar arc plans decreases from 1.19 to 1.10 when increases couch positions. Maximum dose in ipsilateral lung (1.98 to 1.13 Gy), and heart (0.62 to 0.43 Gy) are steadily decreased with the increased number of non-coplanar arcs. Conclusions: The dosimetric evaluation results show that partial arc plans have improved CIs compared to conventional 3D-CRT and IMRT plans. Increasing number of partial arcs decreases lung and heart dose. The dosimetric benefit obtained from non-coplanar arcs should be considered with treatment delivery time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-253: Development of a GDPM Monte Carlo Based Quality Assurance Tool for Cyberknife

Purpose: To develop a 3D quality assurance (QA) tool for Cyberknife plans using a GPU-based Monte Carlo dose calculation package, gDPM. The developed QA tool will be used as a secondary check for Cyberknife treatment planning software (TPS) reported dose distributions. Methods: Phase space files of the 5, 7.5, 10, 15, and 60 mm iris cones were downloaded from IAEA phase-space database. From the 60mm iris cone phase space file, we were able to generate raw phase space files for all other cones using a particle rejection method. The particles in the raw phase space files were binned with respect to the radial position and energy for commissioning. During the commissioning process, gDPM calculated percent depth dose and off-center ratios which were compared to those measured in the Blue Phantom. The off-axis data was measured with an SAD setup at depths of 15, 50, 100, 200, and 300 mm. At the commissioning stage, an optimization problem was solved to adjust the binned particles weights to minimize the difference between the calculated and measured ones. Commissioning validations will be performed by measuring point dose data of patient specific plans delivered to water phantoms. With commissioned phase space files, 3D patient-specific dose distributions will be calculated and compared against TPS reported dose. Results: Beyond the initial buildup region, the root mean square difference between the calculated and measured percent depth dose differences was less than 0.5%. The full-width half-max data from the off-axis ratio calculations was found to be within 0.1 mm of the measured data. Conclusion: Our current work showed excellent agreement between the gDPM calculated dose and the measured data for all cone sizes and types for the Cyberknife system. Implementation of a fast and accurate QA tool for patient specific plans will be feasible with this tool.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-334: Dosimetric Impacts of Organ at Risk (OAR) Contour Errors in Prostate Volumetric Modulated Arc Therapy (VMAT) Planning: A Sensitivity Study

Purpose: To investigate the dosimetric impacts of OAR contour errors in prostate VMAT planning with a sensitivity analysis method. Methods: Ten randomly selected prostate VMAT patients were used to simulate OAR contour error in rectum and bladder. For each OAR 12 even spaced contour control points were selected in every three slices. For each simulation, only one of the control points was simulated to move inwards or outwards by up to 6 voxels and coordinates of adjacent voxels to that moving control point were moved accordingly with a 3D-spline smooth function. An in-house software was used to predict OARs dosimetric endpoints based on geometric relationship between PTV and OAR. In this study, the V75, V70, V65, V60 of rectum and V80, V75, V70, V65 of bladder with and without perturbed contours were calculated and compared. Results: The percentage of OAR dose volume difference around the reference OAR contours were plotted as iso-error lines overlaid on CT images. The significant difference was shown in OAR contours adjacent to the PTV where high dose gradient exists. When one of the 12 points adjacent to the PTV moved up to 6 voxels, maximum errors of 1.99%, 2.46%, 2.89%, and 3.19% were found in V75, V70, V65, and V60 of the rectum, and 1.59%, 1.67%, 1.87%, and 1.90% were found in V80, V75, V70, and V65 of the bladder, respectively. Conclusion: We can quantify the dosimetric impact of OAR contouring error by evaluating the percentage of OAR dose volume changes in prostate VMAT planning. For the adaptive radiation therapy, the iso-error lines can provide the planner a guideline that which portion of contour needed to be checked carefully because dosimetric sensitivity of contour errors are different for different part of contours.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-500: Initial Implementation of GPU-Based Particle Swarm Optimization for 4D IMRT Planning in Lung SBRT

Purpose 4D-IMRT planning, combined with dynamic MLC tracking delivery, utilizes the temporal dimension as an additional degree of freedom to achieve improved OAR-sparing. The computational complexity for such optimization increases exponentially with increase in dimensionality. In order to accomplish this task in a clinically-feasible time frame, we present an initial implementation of GPU-based 4D-IMRT planning based on particle swarm optimization (PSO). Methods The target and normal structures were manually contoured on ten phases of a 4DCT scan of a NSCLC patient with a 54cm3 right-lower-lobe tumor (1.5cm motion). Corresponding ten 3D-IMRT plans were created in the Eclipse treatment planning system (Ver-13.6). A vendor-provided scripting interface was used to export 3D-dose matrices corresponding to each control point (10 phases × 9 beams × 166 control points = 14,940), which served as input to PSO. The optimization task was to iteratively adjust the weights of each control point and scale the corresponding dose matrices. In order to handle the large amount of data in GPU memory, dose matrices were sparsified and placed in contiguous memory blocks with the 14,940 weight-variables. PSO was implemented on CPU (dual-Xeon, 3.1GHz) and GPU (dual-K20 Tesla, 2496 cores, 3.52Tflops, each) platforms. NiftyReg, an open-source deformable image registration package, was used to calculate the summed dose. Results The 4D-PSO plan yielded PTV coverage comparable to the clinical ITV-based plan and significantly higher OAR-sparing, as follows: lung Dmean=33%; lung V20=27%; spinal cord Dmax=26%; esophagus Dmax=42%; heart Dmax=0%; heart Dmean=47%. The GPU-PSO processing time for 14940 variables and 7 PSO-particles was 41% that of CPU-PSO (199 vs. 488 minutes). Conclusion Truly 4D-IMRT planning can yield significant OAR dose-sparing while preserving PTV coverage. The corresponding optimization problem is large-scale, non-convex and computationally rigorous. Our initial results indicate that GPU-based PSO with further software optimization can make such planning clinically feasible. This work was supported through funding from the National Institutes of Health and Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-812: Volumetric Modulated Arc Therapy-Total Body Irradiation (VMAT-TBI) V.s. Conventional Extended SSD-TBI (cTBI): A Dosimetric Comparisom

Purpose: To perform a dosimetric evaluation on a new developed volumetric modulated arc therapy based total body irradiation (VMAT-TBI). Methods: Three patients were CT scanned with an indexed rotatable body frame to get whole body CT images. Concatenated CT images were imported in Pinnacle treatment planning system and whole body and lung were contoured as PTV and organ at risk, respectively. Treatment plans were generated by matching multiple isocenter volumetric modulated arc (VMAT) fields of the upper body and multiple isocenter parallel-opposed fields of the lower body. For each plan, 1200 cGy in 8 fractions was prescribed to the whole body volume and the lung dose was constrained to a mean dose of 750 cGy. Such a two-level dose plan was achieved by inverse planning of the torso VMAT fields. For comparison, conventional standing TBI (cTBI) plans were generated on the same whole body CT images at an extended SSD (550cm).The shape of compensators and lung blocks are simulated using body segments and lung contours Compensation was calculated based on the patient CT images, in mimic of the standing TBI treatment. The whole body dose distribution of cTBI plans were calculated with a home-developed GPU Monte Carlo dose engine. Calculated cTBI dose distribution was prescribed to the mid-body point at umbilical level. Results: The VMAT-TBI treatment plans of three patients’ plans achieved 80.2%±5.0% coverage of the total body volume within ±10% of the prescription dose, while cTBI treatment plans achieved 72.2%±4.0% coverage of the total body volume. The averaged mean lung dose of all three patients is lower for VMAT-TBI (7.48 cGy) than for cTBI (8.96 cGy). Conclusion: The proposed patient comfort-oriented VMAT-TBI technique provides for a uniform dose distribution within the total body while reducing the dose to the lungs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MO-DE-210-06: Development of a Supercompounded 3D Volumetric Ultrasound Image Guidance System for Prone Accelerated Partial Breast Irradiation (APBI)

Purpose: For early-stage breast cancer, accelerated partial breast irradiation (APBI) is a cost-effective breast-conserving treatment. Irradiation in a prone position can mitigate respiratory induced breast movement and achieve maximal sparing of heart and lung tissues. However, accurate dose delivery is challenging due to breast deformation and lumpectomy cavity shrinkage. We propose a 3D volumetric ultrasound (US) image guidance system for accurate prone APBI Methods: The designed system, set beneath the prone breast board, consists of a water container, an US scanner, and a two-layer breast immobilization cup. The outer layer of the breast cup forms the inner wall of water container while the inner layer is attached to patient breast directly to immobilization. The US transducer scans is attached to the outer-layer of breast cup at the dent of water container. Rotational US scans in a transverse plane are achieved by simultaneously rotating water container and transducer, and multiple transverse scanning forms a 3D scan. A supercompounding-technique-based volumetric US reconstruction algorithm is developed for 3D image reconstruction. The performance of the designed system is evaluated with two custom-made gelatin phantoms containing several cylindrical inserts filled in with water (11% reflection coefficient between materials). One phantom is designed for positioning evaluation while the other is for scaling assessment. Results: In the positioning evaluation phantom, the central distances between the inserts are 15, 20, 30 and 40 mm. The distances on reconstructed images differ by −0.19, −0.65, −0.11 and −1.67 mm, respectively. In the scaling evaluation phantom, inserts are 12.7, 19.05, 25.40 and 31.75 mm in diameter. Measured inserts’ sizes on images differed by 0.23, 0.19, −0.1 and 0.22 mm, respectively. Conclusion: The phantom evaluation results show that the developed 3D volumetric US system can accurately localize target position and determine target volume, and is a promising image-guidance tool for prone APBI.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-F-BRF-08: Conformal Mapping-Based 3D Surface Matching and Registration

Purpose: Recently, non-rigid 3D surface matching and registration has been used extensively in engineering and medicine. However, matching 3D surfaces undergoing non-rigid deformation accurately is still a challenging mathematical problem. In this study, we present a novel algorithm to address this issue by introducing intrinsic symmetry to the registration Methods: Our computational algorithm for symmetric conformal mapping is divided into three major steps: 1) Finding the symmetric plane; 2) Finding feature points; and 3) Performing cross registration. The key strategy is to preserve the symmetry during the conformal mapping, such that the image on the parameter domain is symmetric and the area distortion factor on the parameter image is also symmetric. Several novel algorithms were developed using different conformal geometric tools. One was based on solving Riemann-Cauchy equation and the other one employed curvature flow Results: Our algorithm was implemented using generic C++ on Windows XP and used conjugate gradient search optimization for acceleration. The human face 3D surface images were acquired using a high speed 3D scanner based on the phase-shifting method. The scanning speed was 30 frames/sec. The image resolution for each frame was 640 × 480. For 3D human face surfaces with different expressions, postures, and boundaries, our algorithms were able to produce consistent result on the texture pattern on the overlapping region Conclusion: We proposed a novel algorithm to improve the robustness of conformal geometric methods by incorporating the symmetric information into the mapping process. To objectively evaluate its performance, we compared it with most existing techniques. Experimental results indicated that our method outperformed all the others in terms of robustness. The technique has a great potential in real-time patient monitoring and tracking in image-guided radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-F-BRF-09: A Non-Rigid Point Matching Method for Accurate Bladder Dose Summation in Cervical Cancer HDR Brachytherapy

Purpose: To propose and validate a deformable point matching scheme for surface deformation to facilitate accurate bladder dose summation for fractionated HDR cervical cancer treatment. Method: A deformable point matching scheme based on the thin plate spline robust point matching (TPSRPM) algorithm is proposed for bladder surface registration. The surface of bladders segmented from fractional CT images is extracted and discretized with triangular surface mesh. Deformation between the two bladder surfaces are obtained by matching the two meshes' vertices via the TPS-RPM algorithm, and the deformation vector fields (DVFs) characteristic of this deformation is estimated by B-spline approximation. Numerically, the algorithm is quantitatively compared with the Demons algorithm using five clinical cervical cancer cases by several metrics: vertex-to-vertex distance (VVD), Hausdorff distance (HD), percent error (PE), and conformity index (CI). Experimentally, the algorithm is validated on a balloon phantom with 12 surface fiducial markers. The balloon is inflated with different amount of water, and the displacement of fiducial markers is benchmarked as ground truth to study TPS-RPM calculated DVFs' accuracy. Results: In numerical evaluation, the mean VVD is 3.7(±2.0) mm after Demons, and 1.3(±0.9) mm after TPS-RPM. The mean HD is 14.4 mm after Demons, and 5.3mm after TPS-RPM. The mean PE is 101.7% after Demons and decreases to 18.7% after TPS-RPM. The mean CI is 0.63 after Demons, and increases to 0.90 after TPS-RPM. In the phantom study, the mean Euclidean distance of the fiducials is 7.4±3.0mm and 4.2±1.8mm after Demons and TPS-RPM, respectively. Conclusions: The bladder wall deformation is more accurate using the feature-based TPS-RPM algorithm than the intensity-based Demons algorithm, indicating that TPS-RPM has the potential for accurate bladder dose deformation and dose summation for multi-fractional cervical HDR brachytherapy. This work is supported in part by the National Natural ScienceFoundation of China (no 30970866 and no 81301940)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-F-BRF-12: Investigating Dosimetric Effects of Inter-Fraction Deformation in Lung Cancer Stereotactic Body Radiotherapy (SBRT)

Purpose: We studied dosimetric effects of inter-fraction deformation in lung stereotactic body radiotherapy (SBRT), in order to investigate the necessity of adaptive re-planning for lung SBRT treatments. Methods: Six lung cancer patients with different treatment fractions were retrospectively investigated. All the patients were immobilized and localized with a stereotactic body frame and were treated under cone-beam CT (CBCT) image guidance at each fraction. We calculated the actual delivered dose of the treatment plan using the up-to-date patient geometry of each fraction, and compared the dose with the intended plan dose to investigate the dosimetric effects of the inter-fraction deformation. Deformable registration was carried out between the treatment planning CT and the CBCT of each fraction to obtain deformed planning CT for more accurate dose calculations of the delivered dose. The extent of the inter-fraction deformation was also evaluated by calculating the dice similarity coefficient between the delineated structures on the planning CT and those on the deformed planning CT. Results: The average dice coefficients for PTV, spinal cord, esophagus were 0.87, 0.83 and 0.69, respectively. The volume of PTV covered by prescription dose was decreased by 23.78% on average for all fractions and all patients. For spinal cord and esophagus, the volumes covered by the constraint dose were increased by 4.57% and 3.83%. The maximum dose was also increased by 4.11% for spinal cord and 4.29% for esophagus. Conclusion: Due to inter-fraction deformation, large deterioration was found in both PTV coverage and OAR sparing, which demonstrated the needs for adaptive re-planning of lung SBRT cases to improve target coverage while reducing radiation dose to nearby normal tissues.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

MO-C-17A-11: A Segmentation and Point Matching Enhanced Deformable Image Registration Method for Dose Accumulation Between HDR CT Images

Purpose: To propose and validate a novel and accurate deformable image registration (DIR) scheme to facilitate dose accumulation among treatment fractions of high-dose-rate (HDR) gynecological brachytherapy. Method: We have developed a method to adapt DIR algorithms to gynecologic anatomies with HDR applicators by incorporating a segmentation step and a point-matching step into an existing DIR framework. In the segmentation step, random walks algorithm is used to accurately segment and remove the applicator region (AR) in the HDR CT image. A semi-automatic seed point generation approach is developed to obtain the incremented foreground and background point sets to feed the random walks algorithm. In the subsequent point-matching step, a feature-based thin-plate spline-robust point matching (TPS-RPM) algorithm is employed for AR surface point matching. With the resulting mapping, a DVF characteristic of the deformation between the two AR surfaces is generated by B-spline approximation, which serves as the initial DVF for the following Demons DIR between the two AR-free HDR CT images. Finally, the calculated DVF via Demons combined with the initial one serve as the final DVF to map doses between HDR fractions. Results: The segmentation and registration accuracy are quantitatively assessed by nine clinical HDR cases from three gynecological cancer patients. The quantitative results as well as the visual inspection of the DIR indicate that our proposed method can suppress the interference of the applicator with the DIR algorithm, and accurately register HDR CT images as well as deform and add interfractional HDR doses. Conclusions: We have developed a novel and robust DIR scheme that can perform registration between HDR gynecological CT images and yield accurate registration results. This new DIR scheme has potential for accurate interfractional HDR dose accumulation. This work is supported in part by the National Natural ScienceFoundation of China (no 30970866 and no 81301940)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

MO-G-BRE-01: A Real-Time Virtual Delivery System for Photon Radiotherapy Delivery Monitoring

Purpose: Treatment delivery monitoring is important for radiotherapy, which enables catching dosimetric error at the earliest possible opportunity. This project develops a virtual delivery system to monitor the dose delivery process of photon radiotherapy in real-time using GPU-based Monte Carlo (MC) method. Methods: The simulation process consists of 3 parallel CPU threads. A thread T1 is responsible for communication with a linac, which acquires a set of linac status parameters, e.g. gantry angles, MLC configurations, and beam MUs every 20 ms. Since linac vendors currently do not offer interface to acquire data in real time, we mimic this process by fetching information from a linac dynalog file at the set frequency. Instantaneous beam fluence map (FM) is calculated. A FM buffer is also created in T1 and the instantaneous FM is accumulated to it. This process continues, until a ready signal is received from thread T2 on which an inhouse developed MC dose engine executes on GPU. At that moment, the accumulated FM is transferred to T2 for dose calculations, and the FM buffer in T1 is cleared. Once the calculation finishes, the resulting 3D dose distribution is directed to thread T3, which displays it in three orthogonal planes overlaid on the CT image for treatment monitoring. This process continues to monitor the 3D dose distribution in real-time. Results: An IMRT and a VMAT cases used in our patient-specific QA are studied. Maximum dose differences between our system and treatment planning system are 0.98% and 1.58% for the two cases, respectively. The average time per MC calculation is 0.1sec with <2% relative uncertainty. The update frequency of ∼10Hz is considered as real time. Conclusion: By embedding a GPU-based MC code in a novel data/work flow, it is possible to achieve real-time MC dose calculations to monitor delivery process.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-D-9A-02: Automated Landmark-Guided CT to Cone-Beam CT Deformable Image Registration

Purpose: The anatomical changes that occur between the simulation CT and daily cone-beam CT (CBCT) are investigated using an automated landmark-guided deformable image registration (LDIR) algorithm with simultaneous intensity correction. LDIR was designed to be accurate in the presence of tissue intensity mismatch and heavy noise contamination. Method: An auto-landmark generation algorithm was used in conjunction with a local small volume (LSV) gradient matching search engine to map corresponding landmarks between the CBCT and planning CT. The LSVs offsets were used to perform an initial deformation, generate landmarks, and correct local intensity mismatch. The landmarks act as stabilizing controlpoints in the Demons objective function. The accuracy of the LDIR algorithm was evaluated on one synthetic case with ground truth and data of ten head and neck cancer patients. The deformation vector field (DVF) accuracy was accessed using a synthetic case. The Root mean square error of the 3D canny edge (RMSECE), mutual information (MI), and feature similarity index metric (FSIM) were used to access the accuracy of LDIR on the patient data. The quality of the corresponding deformed contours was verified by an attending physician. Results: The resulting 90 percentile DVF error for the synthetic case was within 5.63mm for the original demons algorithm, 2.84mm for intensity correction alone, 2.45mm using controlpoints without intensity correction, and 1.48 mm for the LDIR algorithm. For the five patients the mean RMSECE of the original CT, Demons deformed CT, intensity corrected Demons CT, control-point stabilized deformed CT, and LDIR CT was 0.24, 0.26, 0.20, 0.20, and 0.16 respectively. Conclusion: LDIR is accurate in the presence of multimodal intensity mismatch and CBCT noise contamination. Since LDIR is GPU based it can be implemented with minimal additional strain on clinical resources. This project has been supported by a CPRIT individual investigator award RP11032.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-A-9A-02: BEST IN PHYSICS (THERAPY) - 4D IMRT Planning Using Highly- Parallelizable Particle Swarm Optimization

Purpose: We present a particle swarm optimization (PSO)-based 4D IMRT planning technique designed for dynamic MLC tracking delivery to lung tumors. The key idea is to utilize the temporal dimension as an additional degree of freedom rather than a constraint in order to achieve improved sparing of organs at risk (OARs). Methods: The target and normal structures were manually contoured on each of the ten phases of a 4DCT scan acquired from a lung SBRT patient who exhibited 1.5cm tumor motion despite the use of abdominal compression. Corresponding ten IMRT plans were generated using the Eclipse treatment planning system. These plans served as initial guess solutions for the PSO algorithm. Fluence weights were optimized over the entire solution space i.e., 10 phases × 12 beams × 166 control points. The size of the solution space motivated our choice of PSO, which is a highly parallelizable stochastic global optimization technique that is well-suited for such large problems. A summed fluence map was created using an in-house B-spline deformable image registration. Each plan was compared with a corresponding, internal target volume (ITV)-based IMRT plan. Results: The PSO 4D IMRT plan yielded comparable PTV coverage and significantly higher dose—sparing for parallel and serial OARs compared to the ITV-based plan. The dose-sparing achieved via PSO-4DIMRT was: lung Dmean = 28%; lung V20 = 90%; spinal cord Dmax = 23%; esophagus Dmax = 31%; heart Dmax = 51%; heart Dmean = 64%. Conclusion: Truly 4D IMRT that uses the temporal dimension as an additional degree of freedom can achieve significant dose sparing of serial and parallel OARs. Given the large solution space, PSO represents an attractive, parallelizable tool to achieve globally optimal solutions for such problems. This work was supported through funding from the National Institutes of Health and Varian Medical Systems. Amit Sawant has research funding from Varian Medical Systems, VisionRT Ltd. and Elekta.

60 APPLIED LIFE SCIENCES↗