Engineering Papers⌕ Search

Engineering topics

Liu, C

Publications and source records attributed to Liu, C.

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

SU-F-T-566: Absolute Film Dosimetry for Stereotactic Radiosurgery and Stereotactic Body Radiotherapy Quality Assurance Using Gafchromic EBT3 Films

Purpose: To evaluate the dosimetric uncertainty associated with Gafchromic (EBT3) films and establish an absolute dosimetry protocol for Stereotactic Radiosurgery (SRS) and Stereotactic Body Radiotherapy (SBRT). Methods: EBT3 films were irradiated at each of seven different dose levels between 1 and 15 Gy with open fields, and standard deviations of dose maps were calculated at each color channel for evaluation. A scanner non-uniform response correction map was built by registering and comparing film doses to the reference diode array-based dose map delivered with the same doses. To determine the temporal dependence of EBT3 films, the average correction factors of different dose levels as a function of time were evaluated up to four days after irradiation. An integrated film dosimetry protocol was developed for dose calibration, calibration curve fitting, dose mapping, and profile/gamma analysis. Patient specific quality assurance (PSQA) was performed for 93 SRS/SBRT treatment plans. Results: The scanner response varied within 1% for the field sizes less than 5 × 5 cm{sup 2}, and up to 5% for the field sizes of 10 × 10 cm{sup 2}. The scanner correction method was able to remove visually evident, irregular detector responses found for larger field sizes. The dose response of the film changed rapidly (∼10%) in the first two hours and plateaued afterwards, ∼3% change between 2 and 24 hours. The mean uncertainties (mean of the standard deviations) were <0.5% over the dose range 1∼15Gy for all color channels for the OD response curves. The percentage of points passing the 3%/1mm gamma criteria based on absolute dose analysis, averaged over all tests, was 95.0 ± 4.2. Conclusion: We have developed an absolute film dose dosimetry protocol using EBT3 films. The overall uncertainty has been established to be approximately 1% for SRS and SBRT PSQA. The work was supported by a Research Scholar Grant, RSG-15-137-01-CCE from the American Cancer Society.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRB-05: Automation of the Photon Dosimetric Quality Assurance Program of a Linear Accelerator

Purpose: To develop an automated method to calculate a linear accelerator (LINAC) photon radiation field size, flatness, symmetry, output and beam quality in a single delivery for flattened (FF) and flattening-filter-free (FFF) beams using an ionization chamber array. Methods: The proposed method consists of three control points that deliver 30×30, 10×10 and 5×5cm{sup 2} fields (FF or FFF) in a step-and-shoot sequence where the number of monitor units is weighted for each field size. The IC Profiler (Sun Nuclear Inc.) with 5mm detector spacing was used for this study. The corrected counts (CCs) were calculated and the locations of the maxima and minima values of the first-order gradient determined data of each sub field. Then, all CCs for each field size are summed in order to obtain the final profiles. For each profile, the radiation field size, symmetry, flatness, output factor and beam quality were calculated. For field size calculation, a parameterized gradient method was used. For method validation, profiles were collected in the detector array both, individually and as part of the step-and-shoot plan, with 9.9cm buildup for FF and FFF beams at 90cm source-to-surface distance. The same data were collected with the device (plus buildup) placed on a movable platform to achieve a 1mm resolution. Results: The differences between the dosimetric quantities calculated from both deliveries, individually and step-and-shoot, were within 0.31±0.20% and 0.04±0.02mm. The differences between the calculated field sizes with 5mm and 1mm resolution were ±0.1mm. Conclusion: The proposed single delivery method proved to be simple and efficient in automating the photon dosimetric monthly and annual quality assurance.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-G-TeP1-10: Predicting Photon Percentage Depth Dose From Profile Due to Changes in Beam Quality

Purpose: To predict photon percentage depth dose (PDD) from profile due to a change in flattened (FF) and flattening-filter-free (FFF) beam quality. Methods: 6MV photon beam PDDs and profiles in a 3D water tank (3DW) and profiles in an ionization chamber array (ICP) were collected for different field sizes and depths with FF and FFF beams in a Versa HD (Elekta Ltd.). The energy was adjusted by changing the bending magnet current (BMC) ±15% from the clinical beam (6MV) in 5% increments. For baseline establishment, PDDs(depth≥3cm) were parameterized with bi-exponential functions and the PDD 20 to 10cm ratios (PDD{sub 20,10}) were calculated. Then, the FF profile at 10cm from the central axis (Pr{sub 10}) and the slope of the FFF central linear region (SFFF) were calculated. Calibration curves were established: (1) change in Pr{sub 10} and SFFF as functions of the change in PDD{sub 20,10} and (2) change in PDD(depth=3, 15 and 30cm) as function of the change in PDD{sub 20,10}. The differences between Pr{sub 10} and SFFF from baseline were calculated and, from calibration curves, changes in PDD{sub 20,10} and PDD(depth=3, 15 and 30cm) were obtained. Then, absolute PDD(depth=3, 15 and 30cm) values were input into a least-square-optimization algorithm to calculate the bi-exponential function’s optimal coefficients and generate the PDD(depths≥3cm). Results: The change in PDD{sub 20,10} relative to baseline increased (<±4%) with BMC. Pr{sub 10} increased (±6%) and SFFF decreased (±11%) with BMC. Relative differences between measured and calculated (i.e. PDD calculation from Pr{sub 10} and SFFF) PDDs were less than 1%. Results apply to FF and FFF beams measured in 3DW and ICP. Conclusion: Pr{sub 10} and SFFF are more sensitive than PDD to changes in beam energy and PDD information can be accurately generated from them. With known 3DW and ICP profile relationship, ICP can be used to obtain PDD for current photon beam.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-G-TeP4-03: A Multileaf Collimator Calibration and Quality Assurance Technique Using An Electronic Portal Imaging Device

Purpose: To develop an accurate and quick multileaf collimator (MLC) calibration and quality assurance technique using an electronic portal imaging device (EPID) Methods: The MLC models used include the MLCi and Agility (Elekta Ltd). This technique consists of two 22(L)x10(W) cm{sup 2} fields with 0{sup 0} and 180{sup 0} collimator angles centered to an offset EPID. The MLC opening is estimated by calculating the profile at the image’s center in the image’s horizontal direction. Scans in the image’s vertical direction were calculated every 20 pixels in the inner 70% of estimated MLC opening. The profiles’ edges were fitted with linear equations to determine the image’s rotation angle. Then, crossline profiles were scanned at the center of each leaf taking into account the leaf’s width at isocenter and the rotation angle. The profiles’ edges determine the location of the leaves’ edges and these were subtracted from the reference leaf’s position in order to determine the relative leaf offsets. The edge location of all profiles was determined by using the parameterized gradient of the penumbra region. The technique was tested against an established diode array-based method, and for different MLC systems, patterns, gantry angles, days, energies, beam modalities and MLC openings. Results: The differences between the proposed and established methods were 0.26±0.19mm. The leaf offsets’ deviation was <0.3mm (5 months period). For pattern fields, the differences between predetermined and calculated offsets were 0.18±0.18mm. The leaf offset deviation of measurements with different energies and MLC openings were <0.1mm and <0.3mm, respectively. The differences between offsets of FF and FFF beams were 0.01±0.02mm (<0.07mm). The differences between the offsets at different gantry angles were 0.08±0.15mm. Conclusion: The proposed method proved to be accurate and efficient in calculating the relative leaf offsets. Parameterized field edge is essential to obtain accurate result by eliminating the noise from EPID.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-BRC-01: A New Fabrication Method for Secondary Skin Collimation Using 3D Scanner and 3D Printer

Purpose: Using secondary skin collimation (SSC) to protect the critical organ surrounding the tumors is always desirable for electron and/or ortho-voltage treatments. However, the time-consuming and labor-intensive fabrication processes of current method (manually fabricate the lead mask) restrict the general use of SSC. The aim of this study is to develop a new SSC fabrication method using 3D scanning and 3D printing technologies in order to largely decrease the human labor involvement and fabrication time, also improve the fabrication accuracy. Methods: First, the patient surface was scanned with a 3D scanner (Structure Sensor, Occipital, CO). The raw scan data was subsequently transferred to a 3D modeling software (Rhinoceros, Ver. 5.0, McNeel North America, Seattle, WA). The tumor contour was then digitized and shielding region was determined by clinicians in the same software. The corresponding SCC conformed to the skin surface was then automatically generated by the software with the proper shielding thickness. The shell of the SCC (with hollow inside) was consequently printed by a 3D printer (Lulzbot TAZ, Aleph Objects, CO) using plastic material. Finally, the hollow mold of SCC was filled up with a melted cerrobend alloy. Once the cerrobend alloy cooled down, the fabrication of SCC was accomplished. Results: The results indicated the proposed method can achieve a much shorter time on making a SCC compared with tradition fabrication method. The processes of making a skin contour model for patients have been eliminated with the new method. SCC created by the new method possessed better accuracy and better conformality to patient’s skin contours. Conclusion: In this study, we have demonstrated a new method for the SCC fabrication. It is anticipated that our method can be an alternative way to replace conventional manual-based methods for electron and/or ortho-voltage SCC fabrication. This research was supported by the Global Ph.D. Fellowship Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (Grant No. 2015H1A2A1034071), Republic of Korea.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207B-05: Correlation of Normal Lung Density Changes with Dose After Stereotactic Body Radiotherapy (SBRT) for Early Stage Lung Cancer

Purpose: To investigate correlation of normal lung CT density changes with dose accuracy and outcome after SBRT for patients with early stage lung cancer. Methods: Dose distributions for patients originally planned and treated using a 1-D pencil beam-based (PB-1D) dose algorithm were retrospectively recomputed using algorithms: 3-D pencil beam (PB-3D), and model-based Methods: AAA, Acuros XB (AXB), and Monte Carlo (MC). Prescription dose was 12 Gy × 4 fractions. Planning CT images were rigidly registered to the followup CT datasets at 6–9 months after treatment. Corresponding dose distributions were mapped from the planning to followup CT images. Following the method of Palma et al .(1–2), Hounsfield Unit (HU) changes in lung density in individual, 5 Gy, dose bins from 5–45 Gy were assessed in the peri-tumor region, defined as a uniform, 3 cm expansion around the ITV(1). Results: There is a 10–15% displacement of the high dose region (40–45 Gy) with the model-based algorithms, relative to the PB method, due to the electron scattering of dose away from the tumor into normal lung tissue (Fig.1). Consequently, the high-dose lung region falls within the 40–45 Gy dose range, causing an increase in HU change in this region, as predicted by model-based algorithms (Fig.2). The patient with the highest HU change (∼110) had mild radiation pneumonitis, and the patient with HU change of ∼80–90 had shortness of breath. No evidence of pneumonitis was observed for the 3 patients with smaller CT density changes (<50 HU). Changes in CT densities, and dose-response correlation, as computed with model-based algorithms, are in excellent agreement with the findings of Palma et al. (1–2). Conclusion: Dose computed with PB (1D or 3D) algorithms was poorly correlated with clinically relevant CT density changes, as opposed to model-based algorithms. A larger cohort of patients is needed to confirm these results. This work was supported in part by a grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207B-09: Margin Reduction for Planning Target Volume (PTV) in Patients with Localized Prostate Cancer: Impact On Delivered Dose and Quality of Life

Purpose: To estimate the delivered (cumulative) dose to targets and organs at risk for localized prostate cancer patients treated with reduced PTV margins and to evaluate preliminary patient reported quality-of-life (QOL). Methods: Under an IRB-approved protocol, 20 prostate cancer patients (including 11 control patients) were treated with reduced planning margins (5 mm uniform with 4 mm at prostate/rectum interface). Control patients had standard margin (10/6 mm)-based treatments. A parameter-optimized Elastix algorithm along with energy-mass mapping was used to deform and resample dose of the day onto the planning CT for each fraction to estimate the delivered dose over all fractions. QOL data were collected via Expanded Prostate cancer Index Composite (EPIC-26) questionnaires at time points pre-treatment, post-treatment, and at 2, 6, 12, 18 month follow-ups. Standardized QOL scores [range: 0–100] were determined and baseline-corrected by subtracting pre-treatment QOL data. Mean QOL differences between the margin reduced group and control group (QOLmr-QOLcontrol) were calculated for first 18 months. Results: The difference between the cumulative mean dose (Dmean) and the planned mean dose (±SD) for PTV, prostate, bladder, and rectum were −2.2±1.0, 0.3±0.5, −0.7±2.6, and −2.1±1.3 Gy respectively for the margin-reduced group, and −0.8±2.0, 0.9±1.4, - 0.7±3.1 and −1.0±2.4 Gy for the control group. Difference between the two groups was statistically insignificant (p=0.1). Standardized and baseline corrected QOLmr-QOLcontrol for EPIC domains categorized as “Urinary Incontinence”, “Urinary Irritative/Obstructive”, “Bowel”, “Sexual”, and “Hormonal” were 0.6, 12.1, 9.1, 13.3, and −0.9 for the 18 months following radiation therapy (higher values better). Delivered dose to rectum showed a weak correlation to “Bowel” domain (Pearson’s coefficient −0.24, p<0.001), while bladder dose did not correlate to Urinary Incontinence/Irritative/Obstructive QOL domains. Conclusion: The margin-reduced group exhibited clinically meaningful improvement of QOL without compromising the PTV dose. A larger number of patients and greater follow-up is needed to draw unequivocal conclusions. This work was supported in part by a research grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-C-304-01: Investigation of Various Detector Response Functions and Their Geometry Dependence in a Novel Method to Address Ion Chamber Volume Averaging Effect

Purpose: A novel convolution-based approach has been proposed to address ion chamber (IC) volume averaging effect (VAE) for the commissioning of commercial treatment planning systems (TPS). We investigate the use of various convolution kernels and its impact on the accuracy of beam models. Methods: Our approach simulates the VAE by iteratively convolving the calculated beam profiles with a detector response function (DRF) while optimizing the beam model. At convergence, the convolved profiles match the measured profiles, indicating the calculated profiles match the “true” beam profiles. To validate the approach, beam profiles of an Elekta LINAC were repeatedly collected with ICs of various volumes (CC04, CC13 and SNC 125) to obtain clinically acceptable beam models. The TPS-calculated profiles were convolved externally with the DRF of respective IC. The beam model parameters were reoptimized using Nelder-Mead method by forcing the convolved profiles to match the measured profiles. We evaluated three types of DRFs (Gaussian, Lorentzian, and parabolic) and the impact of kernel dependence on field geometry (depth and field size). The profiles calculated with beam models were compared with SNC EDGE diode-measured profiles. Results: The method was successfully implemented with Pinnacle Scripting and Matlab. The reoptimization converged in ∼10 minutes. For all tested ICs and DRFs, penumbra widths of the TPS-calculated profiles and diode-measured profiles were within 1.0 mm. Gaussian function had the best performance with mean penumbra width difference within 0.5 mm. The use of geometry dependent DRFs showed marginal improvement, reducing the penumbra width differences to less than 0.3 mm. Significant increase in IMRT QA passing rates was achieved with the optimized beam model. Conclusion: The proposed approach significantly improved the accuracy of the TPS beam model. Gaussian functions as the convolution kernel performed consistently better than Lorentzian and parabolic functions. Further improvement of using geometry-dependent kernels was marginal.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-66: Significant Anatomical and Dosimetric Changes Observed with the Pharyngeal Constrictor During Head and Neck Radiotherapy Elicited From Daily Deformable Image Registration and Dose Accumulation

Purpose: To evaluate the anatomical changes and associated dosimetric consequences to the pharyngeal constrictor (PC) that occurs during head and neck radiotherapy (H&N RT). Methods: A cohort of 13 oro-pharyngeal cancer patients, who had daily CBCT’s for localization, was retrospectively studied. On every 5th CBCT, PC was manually delineated by a radiation oncologist. The anterior-posterior PC thickness was measured at the C3 level. Delivered dose to PC was estimated by calculating daily doses on CBCT’s, and accumulating to corresponding planning CT images. For accumulation, a parameter-optimized B- spline-based deformable image registration algorithm (Elastix) was used, in conjunction with an energy-mass mapping dose transfer algorithm. Mean and maximum dose (Dmean, Dmax) to PC was determined and compared with corresponding planned quantities. Results: The mean (±standard deviation) volume increase (ΔV) and thickness increase (Δt) over the course of 35 total fractions were 54±33% (11.9±7.6 cc), and 63±39% (2.9±1.9 mm), respectively. The resultant cumulative mean dose increase from planned dose to PC (ΔDmean) was 1.4±1.3% (0.9±0.8 Gy), while the maximum dose increase (ΔDmax) was 0.0±1.6% (0.0±1.1 Gy). Patients with adaptive replanning (n=6) showed a smaller mean dose increase than those without (n=7); 0.5±0.2% (0.3±0.1 Gy) vs. 2.2±1.4% (1.4±0.9 Gy). There was a statistically significant (p<0.0001) strong correlation between ΔDmean and Δt (Pearson coefficient r=0.78), and a moderate-to-strong correlation (r=0.52) between ΔDmean and ΔV. Correlation between ΔDmean and weight loss ΔW (r=0.1), as well as ΔV and ΔW (r=0.2) were negligible. Conclusion: Patients were found to undergo considerable anatomical changes to pharyngeal constrictor during H&N RT, resulting in non-negligible dose deviations from intended dose. Results are indicative that pharyngeal constrictor thickness, measured at C3 level, is a good predictor for the dose change to the organ. Daily deformable registration and dose accumulation provide a reliable means to assess important anatomical and dosimetric changes to pharyngeal constrictor occurring during treatment. This work was supported in part by a research grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-114: A Practical Hybrid Method for Improving the Quality of CT-CBCT Deformable Image Registration for Head and Neck Radiotherapy

Purpose: Accurate deformable image registration (DIR) between CT and CBCT in H&N is challenging. In this study, we propose a practical hybrid method that uses not only the pixel intensities but also organ physical properties, structure volume of interest (VOI), and interactive local registrations. Methods: Five oropharyngeal cancer patients were selected retrospectively. For each patient, the planning CT was registered to the last fraction CBCT, where the anatomy difference was largest. A three step registration strategy was tested; Step1) DIR using pixel intensity only, Step2) DIR with additional use of structure VOI and rigidity penalty, and Step3) interactive local correction. For Step1, a public-domain open-source DIR algorithm was used (cubic B-spline, mutual information, steepest gradient optimization, and 4-level multi-resolution). For Step2, rigidity penalty was applied on bony anatomies and brain, and a structure VOI was used to handle the body truncation such as the shoulder cut-off on CBCT. Finally, in Step3, the registrations were reviewed on our in-house developed software and the erroneous areas were corrected via a local registration using level-set motion algorithm. Results: After Step1, there were considerable amount of registration errors in soft tissues and unrealistic stretching in the posterior to the neck and near the shoulder due to body truncation. The brain was also found deformed to a measurable extent near the superior border of CBCT. Such errors could be effectively removed by using a structure VOI and rigidity penalty. The rest of the local soft tissue error could be corrected using the interactive software tool. The estimated interactive correction time was approximately 5 minutes. Conclusion: The DIR using only the image pixel intensity was vulnerable to noise and body truncation. A corrective action was inevitable to achieve good quality of registrations. We found the proposed three-step hybrid method efficient and practical for CT/CBCT registrations in H&N. My department receives grant support from Industrial partners: (a) Varian Medical Systems, Palo Alto, CA, and (b) Philips HealthCare, Best, Netherlands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-257: A PCA Model to Predict Adaptive Changes for Head&neck Patients Based On Extraction of Geometric Features From Daily CBCT Datasets

Purpose: Using daily cone beam CTs (CBCTs) to develop principal component analysis (PCA) models of anatomical changes in head and neck (H&N) patients and to assess the possibility of using these prospectively in adaptive radiation therapy (ART). Methods: Planning CT (pCT) images of 4 H&N patients were deformed to model several different systematic changes in patient anatomy during the course of the radiation therapy (RT). A Pinnacle plugin was used to linearly interpolate the systematic change in patient for the 35 fraction RT course and to generate a set of 35 synthetic CBCTs. Each synthetic CBCT represents the systematic change in patient anatomy for each fraction. Deformation vector fields (DVFs) were acquired between the pCT and synthetic CBCTs with random fraction-to-fraction changes were superimposed on the DVFs. A patient-specific PCA model was built using these DVFs containing systematic plus random changes. It was hypothesized that resulting eigenDVFs (EDVFs) with largest eigenvalues represent the major anatomical deformations during the course of treatment. Results: For all 4 patients, the PCA model provided different results depending on the type and size of systematic change in patient’s body. PCA was more successful in capturing the systematic changes early in the treatment course when these were of a larger scale with respect to the random fraction-to-fraction changes in patient’s anatomy. For smaller scale systematic changes, random changes in patient could completely “hide” the systematic change. Conclusion: The leading EDVF from the patientspecific PCA models could tentatively be identified as a major systematic change during treatment if the systematic change is large enough with respect to random fraction-to-fraction changes. Otherwise, leading EDVF could not represent systematic changes reliably. This work is expected to facilitate development of population-based PCA models that can be used to prospectively identify significant anatomical changes early in treatment. This work is supported in part by a grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-584: Optical Tracking Guided Patient-Specific VMAT QA with ArcCHECK

Purpose: To investigate the novel use of an in-house optical tracking system (OTS) to improve the efficacy of VMAT QA with a cylindrical dosimeter (ArcCHECK™). Methods: The translational and rotational setup errors of ArcCHECK are convoluted which makes it challenging to position the device efficiently and accurately. We first aligned the ArcCHECK to the machine cross-hair at three cardinal angles (0°, 90°, and 270°) to establish a reference position. Four infrared reflective markers were attached to the back of the device. An OTS with 0.2mm/0.2° accuracy was used to control its setup uncertainty. Translational uncertainties of 1 mm and 2 mm in three directions (in, right, and up) were applied on the device. Four open beams of various field sizes and six clinical VMAT arcs were delivered and measured for all simulated setup errors. The measurements were compared with Pinnacle™ calculations using Gamma analysis to evaluate the impact of setup uncertainty. This study also evaluated the improvement in setup reproducibility and efficiency with the aid of the OTS. Results: For open beams, with 3%/3mm, the mean passing rates dropped by less than 5% for all shifts; with 2%/2mm, two significant drops(>5%) were observed: 15.38±6.75% for 2 mm lateral shift and 9.35±4.88% for 2 mm longitudinal shift. For VMAT arcs, the mean passing rates using 2%/2mm dropped by 10.47±7.46% and 22.02±11.39% for 1 and 2 mm shift, respectively. With 3%/3mm, significant drop only occurred with 2 mm longitudinal shift (13.73±8.30%). Setup time could be reduced by >15 min with the aid of the OTS. Conclusion: OTS is an effective tool for separating translational and rotational uncertainties. The current VMAT QA solution was not strongly sensitive to translation errors of 2mm with widely accepted criterion (3%/3mm). This finding raises concerns regarding the efficacy of such QA system in detecting errors in the dynamic VMAT delivery.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

SU-E-T-597: Parameterization of the Photon Beam Dosimetry for a Commercial Linear Accelerator

Purpose: In radiation therapy, accurate data acquisition of photon beam dosimetric quantities is important for (1) beam modeling data input into a treatment planning system (TPS), (2) comparing measured and TPS modelled data, (3) a linear accelerator’s (linac) beam characteristics quality assurance process, and (4) establishing a standard data set for data comparison, etcetera. Parameterization of the photon beam dosimetry creates a portable data set that is easy to implement for different applications such as those previously mentioned. The aim of this study is to develop methods to parameterize photon percentage depth doses(PDD), profiles, and total scatter output factors(Scp). Methods: Scp, PDDs and profiles for different field sizes (from 2×2 to 40×40cm{sup 2}), depths and energies were measured in a linac using a three-dimensional water tank. All data were smoothed and profile data were also centered, symmetrized and geometrically scaled. The Scp and PDD data were analyzed using exponential functions. For modelling of open and wedge field profiles, each side was divided into three regions described by exponential, sigmoid and Gaussian equations. The model’s equations were chosen based on the physical principles described by these dosimetric quantities. The equations’ parameters were determined using a least square optimization method with the minimal amount of measured data necessary. The model’s accuracy was then evaluated via the calculation of absolute differences and distance–to–agreement analysis in low gradient and high gradient regions, respectively. Results: All differences in the PDDs’ buildup and the profiles’ penumbra regions were less than 2 mm and 0.5 mm, respectively. Differences in the low gradient regions were 0.20 ± 0.20% and 0.50 ± 0.35% for PDDs and profiles, respectively. For Scp data, all differences were less than 0.5%. Conclusion: This novel analytical model with minimum measurement requirements proved to accurately calculate PDDs, profiles, and Scp for different field sizes, depths and energies.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

MO-F-CAMPUS-T-01: Radiosurgery of Multiple Brain Metastases with Single-Isocenter VMAT: Optimizing Treatment Geometry to Reduce Normal Brain Dose

Purpose: To develop an optimization algorithm to reduce normal brain dose by optimizing couch and collimator angles for single isocenter multiple targets treatment of stereotactic radiosurgery. Methods: Three metastatic brain lesions were retrospectively planned using single-isocenter volumetric modulated arc therapy (VMAT). Three matrices were developed to calculate the projection of each lesion on Beam’s Eye View (BEV) by the rotating couch, collimator and gantry respectively. The island blocking problem was addressed by computing the total area of open space between any two lesions with shared MLC leaf pairs. The couch and collimator angles resulting in the smallest open areas were the optimized angles for each treatment arc. Two treatment plans with and without couch and collimator angle optimization were developed using the same objective functions and to achieve 99% of each target volume receiving full prescription dose of 18Gy. Plan quality was evaluated by calculating each target’s Conformity Index (CI), Gradient Index (GI), and Homogeneity index (HI), and absolute volume of normal brain V8Gy, V10Gy, V12Gy, and V14Gy. Results: Using the new couch/collimator optimization strategy, dose to normal brain tissue was reduced substantially. V8, V10, V12, and V14 decreased by 2.3%, 3.6%, 3.5%, and 6%, respectively. There were no significant differences in the conformity index, gradient index, and homogeneity index between two treatment plans with and without the new optimization algorithm. Conclusion: We have developed a solution to the island blocking problem in delivering radiation to multiple brain metastases with shared isocenter. Significant reduction in dose to normal brain was achieved by using optimal couch and collimator angles that minimize total area of open space between any of the two lesions with shared MLC leaf pairs. This technique has been integrated into Eclipse treatment system using scripting API.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-G-207-03: Mask Guided Image Reconstruction (MGIR): A Novel Method for Ultra-Low-Dose 3D and Enhanced 4D Cone-Beam Computer-Tomography

Purpose: Recently, compressed sensing (CS) based iterative reconstruction (IR) method is receiving attentions to reconstruct high quality cone beam computed tomography (CBCT) images using sparsely sampled or noisy projections. The aim of this study is to develop a novel baseline algorithm called Mask Guided Image Reconstruction (MGIR), which can provide superior image quality for both low-dose 3DCBCT and 4DCBCT under single mathematical framework. Methods: In MGIR, the unknown CBCT volume was mathematically modeled as a combination of two regions where anatomical structures are 1) within the priori-defined mask and 2) outside the mask. Then we update each part of images alternatively thorough solving minimization problems based on CS type IR. For low-dose 3DCBCT, the former region is defined as the anatomically complex region where it is focused to preserve edge information while latter region is defined as contrast uniform, and hence aggressively updated to remove noise/artifact. In 4DCBCT, the regions are separated as the common static part and moving part. Then, static volume and moving volumes were updated with global and phase sorted projection respectively, to optimize the image quality of both moving and static part simultaneously. Results: Examination of MGIR algorithm showed that high quality of both low-dose 3DCBCT and 4DCBCT images can be reconstructed without compromising the image resolution and imaging dose or scanning time respectively. For low-dose 3DCBCT, a clinical viable and high resolution head-and-neck image can be obtained while cutting the dose by 83%. In 4DCBCT, excellent quality 4DCBCT images could be reconstructed while requiring no more projection data and imaging dose than a typical clinical 3DCBCT scan. Conclusion: The results shown that the image quality of MGIR was superior compared to other published CS based IR algorithms for both 4DCBCT and low-dose 3DCBCT. This makes our MGIR algorithm potentially useful in various on-line clinical applications. Provisional Patent: UF#15476; WGS Ref. No. U1198.70067US00.

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↗