Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “histogram”

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

Factorization-based texture segmentation

This study introduces a factorization-based approach that efficiently segments textured images. We use local spectral histograms as features, and construct an M × N feature matrix using M-dimensional feature vectors in an N-pixel image. Based on the observation that each feature can be approximated by a linear combination of several representative features, we factor the feature matrix into two matrices-one consisting of the representative features and the other containing the weights of representative features at each pixel used for linear combination. The factorization method is based on singular value decomposition and nonnegative matrix factorization. The method uses local spectral histograms to discriminate region appearances in a computationally efficient way and at the same time accurately localizes region boundaries. Finally, the experiments conducted on public segmentation data sets show the promise of this simple yet powerful approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Simulations of cloud‐radiation interaction using large‐scale forcing derived from the CINDY/DYNAMO northern sounding array

Abstract The recently completed CINDY/DYNAMO field campaign observed two Madden‐Julian oscillation (MJO) events in the equatorial Indian Ocean from October to December 2011. Prior work has indicated that the moist static energy anomalies in these events grew and were sustained to a significant extent by radiative feedbacks. We present here a study of radiative fluxes and clouds in a set of cloud‐resolving simulations of these MJO events. The simulations are driven by the large‐scale forcing data set derived from the DYNAMO northern sounding array observations, and carried out in a doubly periodic domain using the Weather Research and Forecasting (WRF) model. Simulated cloud properties and radiative fluxes are compared to those derived from the S‐PolKa radar and satellite observations. To accommodate the uncertainty in simulated cloud microphysics, a number of single‐moment (1M) and double‐moment (2M) microphysical schemes in the WRF model are tested. The 1M schemes tend to underestimate radiative flux anomalies in the active phases of the MJO events, while the 2M schemes perform better, but can overestimate radiative flux anomalies. All the tested microphysics schemes exhibit biases in the shapes of the histograms of radiative fluxes and radar reflectivity. Histograms of radiative fluxes and brightness temperature indicate that radiative biases are not evenly distributed; the most significant bias occurs in rainy areas with OLR less than 150 W/m 2 in the 2M schemes. Analysis of simulated radar reflectivities indicates that this radiative flux uncertainty is closely related to the simulated stratiform cloud coverage. Single‐moment schemes underestimate stratiform cloudiness by a factor of 2, whereas 2M schemes simulate much more stratiform cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Using the apparent diffusion coefficient to identifying MGMT promoter methylation status early in glioblastoma: importance of analytical method

Accurate knowledge of O{sup 6}-methylguanine methyltransferase (MGMT) gene promoter subtype in patients with glioblastoma (GBM) is important for treatment. However, this test is not always available. Pre-operative diffusion MRI (dMRI) can be used to probe tumour biology using the apparent diffusion coefficient (ADC); however, its ability to act as a surrogate to predict MGMT status has shown mixed results. We investigated whether this was due to variations in the method used to analyse ADC. We undertook a retrospective study of 32 patients with GBM who had MGMT status measured. Matching pre-operative MRI data were used to calculate the ADC within contrast enhancing regions of tumour. The relationship between ADC and MGMT was examined using two published ADC methods. A strong trend between a measure of ‘minimum ADC’ and methylation status was seen. An elevated minimum ADC was more likely in the methylated compared to the unmethylated MGMT group (U = 56, P = 0.0561). In contrast, utilising a two-mixture model histogram approach, a significant reduction in mean measure of the ‘low ADC’ component within the histogram was associated with an MGMT promoter methylation subtype (P < 0.0246). This study shows that within the same patient cohort, the method selected to analyse ADC measures has a significant bearing on the use of that metric as a surrogate marker of MGMT status. Thus for dMRI data to be clinically useful, consistent methods of data analysis need to be established prior to establishing any relationship with genetic or epigenetic profiling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Juxtatumoral perinephric fat analysis in clear cell renal cell carcinoma

Purpose: The purpose of the study was to evaluate the feasibility of using contrast-enhanced computed tomography (CECT)-based texture analysis (CTTA) metrics to differentiate between juxtatumoral perinephric fat (JPF) surrounding low-grade (ISUP 1–2) versus high-grade (ISUP 3–4) clear cell renal cell carcinoma (ccRCC). Methods: In this IRB-approved study, we retrospectively queried the surgical database between June 2009 and April 2016 and identified 83 patients with pathologically confirmed ccRCC (low grade: n = 54, mean age = 61.5 years, 18F/35M; high grade n = 30, mean age = 61.7 years, 8F/22M) who also had pre-operative multiphase CT acquisitions. CT images were transferred to a 3D workstation, and nephrographic phase JPF regions were manually segmented. Using an in-house developed Matlab program, a CTTA panel comprising of texture metrics extracted using six different methods, histogram, 2D- and 3D-Gray-level co-occurrence matrix (GLCM) and Gray-level difference matrix (GLDM), and 2D-Fast Fourier Transform (FFT) analyses, was applied to the segmented images to assess JPF textural heterogeneity in low- versus high-grade ccRCC. Univariate analysis and receiver-operator characteristics (ROC) analysis were used to assess interclass differences in texture metrics and their prediction accuracy, respectively. Results: All methods except GLCM consistently revealed increased heterogeneity in the JPF surrounding high- versus low-grade ccRCC. FFT showed increased complexity index (p < 0.01). Histogram analysis showed increased kurtosis and positive skewness in (p < 0.03), and GLDM analysis showed decreased measure of correlation coefficient (MCC) (p < 0.04). Several of the GLCM metrics showed statistically significant (p < 0.04) textural differences between the two groups, but with no consistent trend. ROC analysis showed that MCC in GLCM analysis had an area under the curve of 0.75. Conclusions: Our study suggests that CTTA of ccRCC shows statistically significant textural differences in JPF surrounding high- versus low-grade ccRCC.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Measurement uncertainty of rossi-alpha neutron experiments

Rossi-alpha neutron experiments are used to estimate the prompt neutron decay constant of a fissile assembly, a quantity of widespread interest in applications including in nuclear nonproliferation and criticality safety. This work develops a mathematical model to efficiently estimate measurement uncertainty of Rossi-alpha neutron experiments inferred from a two-exponential fit model with histogram binning. The derived uncertainty estimates were validated using repeated Rossi-alpha measurements of a subcritical, 4.5-kg sphere of weapons-grade, alpha-phase plutonium with nickel, copper, tungsten, and polyethylene reflectors. The estimates of uncertainty for the histogram data produced by the model were conservative and agree with the reference uncertainties within noise. The estimates of the prompt neutron decay constant uncertainty agreed with the reference uncertainties within one standard deviation. In conclusion, the proposed model will reduce total measurement times, ultimately reducing operational and procedural costs in application.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Probability Density Function Control for Stochastic Nonlinear Systems using Monte Carlo Simulation

This paper presents a fast implementation framework of output probability density function (PDF) control for a class of stochastic nonlinear systems which are subjected to non-Gaussian noises. The statistical properties of the system outputs can be adjusted by shaping the dynamic output probability density function to track the reference stochastic distribution. However, the dynamic probability density function evolution is very difficult to obtain analytically even if the system model and the stochastic distributions of the noises are known. Motivated by Monte Carlo simulation, the dynamic probability density function can be estimated by sampling data which forms the contribution of this paper. In particular, the sampling points are generated following the stochastic distribution of the noise for each instant. These points go through the system would generate the histogram for system outputs, then the dynamic model can be established based on the dynamic histogram which reflects the randomness and the nonlinear dynamics of the investigated system. Based on the established model, the output probability density function tracking can be achieved and the simulation results and discussions show the effectiveness and benefits of the presented framework.

Zhang, Qichun↗
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

Decreased Risk of Radiation Pneumonitis With Incidental Concurrent Use of Angiotensin-Converting Enzyme Inhibitors and Thoracic Radiation Therapy

Purpose: Angiotensin-converting enzyme (ACE) inhibitors have been shown to mitigate radiation-induced lung injury in preclinical models. The aim of this study was to evaluate whether ACE inhibitors decrease the risk of radiation pneumonitis in lung cancer patients receiving thoracic irradiation. Methods and Materials: Patients with Stage I through III small-cell and non-small-cell lung cancer treated definitively with radiation from 2004-2009 at the Clement J. Zablocki Veterans Affairs Medical Center were retrospectively reviewed. Acute pulmonary toxicity was quantified within 6 months of completion of treatment according to the Common Terminology Criteria for Adverse Events version 4. The use of ACE inhibitors, nonsteroidal anti-inflammatory drugs, inhaled glucocorticosteroids, statins, and angiotensin receptor blockers; dose-volume histogram parameters; and patient factors were assessed for association with Grade 2 or higher pneumonitis. Results: A total of 162 patients met the criteria for inclusion. The majority of patients had Stage III disease (64%) and received concurrent chemotherapy (61%). Sixty-two patients were identified as ACE inhibitor users (38%). All patients had acceptable radiation plans based on dose-volume histogram constraints (V20 [volume of lung receiving at least 20 Gy] {<=}37% and mean lung dose {<=}20 Gy) with the exception of 2 patients who did not meet both criteria. Grade 2 or higher pulmonary toxicity occurred in 12 patients (7.4%). The rate of Grade 2 or higher pneumonitis was lower in ACE inhibitor users vs. nonusers (2% vs. 11%, p = 0.032). Rates of Grade 2 or higher pneumonitis were significantly increased in patients aged greater than 70 years (16% vs. 2%, p = 0.005) or in whom V5 (volume of lung receiving at least 5 Gy) was 50% or greater (13% vs. 4%, p = 0.04). V10 (volume of lung receiving at least 10 Gy), V20, V30 (volume of lung receiving at least 30 Gy), and mean lung dose were not independently associated with Grade 2 or higher pneumonitis. Conclusion: ACE inhibitors may decrease the incidence of radiation pneumonitis in patients receiving thoracic radiation for lung cancer. These findings are consistent with preclinical evidence and should be prospectively evaluated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Parotid Gland Sparing With Helical Tomotherapy in Head-and-Neck Cancer

Purpose: This study evaluated the ability of helical tomotherapy to spare the function of the parotid glands in patients with head-and-neck cancer by analyzing dose-volume histograms, salivary gland scintigraphy, and quality of life assessment. Methods and Materials: Data from 76 consecutive patients treated with helical tomotherapy (Hi-Art Tomotherapy) at University Hospital Brussel were analyzed. During planning, priority was given to planning target volume (PTV) coverage: {>=}95% of the dose must be delivered to {>=}95% of the PTV. Elective nodal regions received 54 Gy (1.8 Gy/fraction). A dose of 70.5 Gy (2.35 Gy/fraction) was prescribed to the primary tumor and pathologic lymph nodes (simultaneous integrated boost scheme). Objective scoring of salivary excretion was performed by salivary gland scintigraphy. Subjective scoring of salivary gland function was evaluated by the European Organization for Research and Treatment of Cancer quality of life questionnaires Quality of Life Questionnaire-C30 (QLQ-C30) and Quality of Life Questionnaire-Head and Neck 35 (H and N35). Results: Analysis of dose-volume histograms (DVHs) showed excellent coverage of the PTV. The volume of PTV receiving 95% of the prescribed dose (V95%) was 99.4 (range, 96.3-99.9). DVH analysis of parotid gland showed a median value of the mean parotid dose of 32.1 Gy (range, 17.5-70.3 Gy). The median parotid volume receiving a dose <26 Gy was 51.2%. Quality of life evaluation demonstrated an initial deterioration of almost all scales and items in QLQ-C30 and QLQ-H and N35. Most items improved in time, and some reached baseline values 18 months after treatment. Conclusion: DVH analysis, scintigraphic evaluation of parotid function, and quality of life assessment of our patient group showed that helical tomotherapy makes it possible to preserve parotid gland function without compromising disease control. We recommend mean parotid doses of <34 Gy and doses <26 Gy to a maximum 47% of the parotid volume as planning goals. Intensity-modulated radiotherapy should be considered as standard treatment in patients with head-and-neck cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Dose to the Developing Dentition During Therapeutic Irradiation: Organ at Risk Determination and Clinical Implications

Purpose: Irradiation of pediatric facial structures can cause severe impairment of permanent teeth later in life. We therefore focused on primary and permanent teeth as organs at risk, investigating the ability to identify individual teeth in children and infants and to correlate dose distributions with subsequent dental toxicity. Methods and Materials: We retrospectively reviewed 14 pediatric patients who received a maximum dose >20 Gy(relative biological effectiveness, RBE) to 1 or more primary or permanent teeth between 2003 and 2009. The patients (aged 1-16 years) received spot-scanning proton therapy with 46 to 66 Gy(RBE) in 23 to 33 daily fractions for a variety of tumors, including rhabdomyosarcoma (n=10), sarcoma (n=2), teratoma (n=1), and carcinoma (n=1). Individual teeth were contoured on axial slices from planning computed tomography (CT) scans. Dose-volume histogram data were retrospectively obtained from total calculated delivered treatments. Dental follow-up information was obtained from external care providers. Results: All primary teeth and permanent incisors, canines, premolars, and first and second molars were identifiable on CT scans in all patients as early as 1 year of age. Dose-volume histogram analysis showed wide dose variability, with a median 37 Gy(RBE) per tooth dose range across all individuals, and a median 50 Gy(RBE) intraindividual dose range across all teeth. Dental follow-up revealed absence of significant toxicity in 7 of 10 patients but severe localized toxicity in teeth receiving >20 Gy(RBE) among 3 patients who were all treated at <4 years of age. Conclusions: CT-based assessment of dose distribution to individual teeth is feasible, although delayed calcification may complicate tooth identification in the youngest patients. Patterns of dental dose exposure vary markedly within and among patients, corresponding to rapid dose falloff with protons. Severe localized dental toxicity was observed in a few patients receiving the largest doses of radiation at the youngest ages; however, multiple factors including concurrent chemotherapy confounded the dose-effect relationship. Further studies with larger cohorts and appropriate controls will be required.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Long-term Results of Carbon Ion Radiation Therapy for Locally Advanced or Unfavorably Located Choroidal Melanoma: Usefulness of CT-based 2-Port Orthogonal Therapy for Reducing the Incidence of Neovascular Glaucoma

Purpose: To determine the long-term results of carbon ion radiation therapy (C-ion RT) in patients with choroidal melanoma, and to assess the usefulness of CT-based 2-port irradiation in reducing the risk of neovascular glaucoma (NVG). Methods and Materials: Between January 2001 and February 2012, a total of 116 patients with locally advanced or unfavorably located choroidal melanoma received CT-based C-ion RT. Of these patients, 114 were followed up for more than 6 months and their data analyzed. The numbers of T3 and T2 patients (International Union Against Cancer [UICC], 5th edition) were 106 and 8, respectively. The total dose of C-ion RT varied from 60 to 85 GyE, with each dose given in 5 fractions. Since October 2005, 2-port therapy (51 patients) has been used in an attempt to reduce the risk of NVG. A dose-volume histogram analysis was also performed in 106 patients. Results: The median follow-up was 4.6 years (range, 0.5-10.6 years). The 5-year overall survival, cause-specific survival, local control, distant metastasis-free survival, and eye retention rates were 80.4% (95% confidence interval 89.0%-71.8%), 82.2% (90.6%-73.8%), 92.8% (98.5%-87.1%), 72.1% (81.9%-62.3%), and 92.8% (98.1%-87.5%), respectively. The overall 5-year NVG incidence rate was 35.9% (25.9%-45.9%) and that of 1-port group and 2-port group were 41.6% (29.3%-54.0%) and 13.9% (3.2%-24.6%) with statistically significant difference (P<.001). The dose-volume histogram analysis showed that the average irradiated volume of the iris-ciliary body was significantly lower in the non-NVG group than in the NVG group at all dose levels, and significantly lower in the 2-port group than in the 1-port group at high dose levels. Conclusions: The long-term results of C-ion RT for choroidal melanoma are satisfactory. CT-based 2-port C-ion RT can be used to reduce the high-dose irradiated volume of the iris-ciliary body and the resulting risk of NVG.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A Novel Method for Predicting Late Genitourinary Toxicity After Prostate Radiation Therapy and the Need for Age-Based Risk-Adapted Dose Constraints

Background: There are no well-established normal tissue sparing dose–volume histogram (DVH) criteria that limit the risk of urinary toxicity from prostate radiation therapy (RT). The aim of this study was to determine which criteria predict late toxicity among various DVH parameters when contouring the entire solid bladder and its contents versus the bladder wall. The area under the histogram curve (AUHC) was also analyzed. Methods and Materials: From 1993 to 2000, 503 men with prostate cancer received 3-dimensional conformal RT (median follow-up time, 71 months). The whole bladder and the bladder wall were contoured in all patients. The primary endpoint was grade ≥2 genitourinary (GU) toxicity occurring ≥3 months after completion of RT. Cox regressions of time to grade ≥2 toxicity were estimated separately for the entire bladder and bladder wall. Concordance probability estimates (CPE) assessed model discriminative ability. Before training the models, an external random test group of 100 men was set aside for testing. Separate analyses were performed based on the mean age (≤ 68 vs >68 years). Results: Age, pretreatment urinary symptoms, mean dose (entire bladder and bladder wall), and AUHC (entire bladder and bladder wall) were significant (P<.05) in multivariable analysis. Overall, bladder wall CPE values were higher than solid bladder values. The AUHC for bladder wall provided the greatest discrimination for late bladder toxicity when compared with alternative DVH points, with CPE values of 0.68 for age ≤68 years and 0.81 for age >68 years. Conclusion: The AUHC method based on bladder wall volumes was superior for predicting late GU toxicity. Age >68 years was associated with late grade ≥2 GU toxicity, which suggests that risk-adapted dose constraints based on age should be explored.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Statistical Assessment of Proton Treatment Plans Under Setup and Range Uncertainties

Purpose: To evaluate a method for quantifying the effect of setup errors and range uncertainties on dose distribution and dose–volume histogram using statistical parameters; and to assess existing planning practice in selected treatment sites under setup and range uncertainties. Methods and Materials: Twenty passively scattered proton lung cancer plans, 10 prostate, and 1 brain cancer scanning-beam proton plan(s) were analyzed. To account for the dose under uncertainties, we performed a comprehensive simulation in which the dose was recalculated 600 times per given plan under the influence of random and systematic setup errors and proton range errors. On the basis of simulation results, we determined the probability of dose variations and calculated the expected values and standard deviations of dose–volume histograms. The uncertainties in dose were spatially visualized on the planning CT as a probability map of failure to target coverage or overdose of critical structures. Results: The expected value of target coverage under the uncertainties was consistently lower than that of the nominal value determined from the clinical target volume coverage without setup error or range uncertainty, with a mean difference of −1.1% (−0.9% for breath-hold), −0.3%, and −2.2% for lung, prostate, and a brain cases, respectively. The organs with most sensitive dose under uncertainties were esophagus and spinal cord for lung, rectum for prostate, and brain stem for brain cancer. Conclusions: A clinically feasible robustness plan analysis tool based on direct dose calculation and statistical simulation has been developed. Both the expectation value and standard deviation are useful to evaluate the impact of uncertainties. The existing proton beam planning method used in this institution seems to be adequate in terms of target coverage. However, structures that are small in volume or located near the target area showed greater sensitivity to uncertainties.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Risk Factors for Neovascular Glaucoma After Proton Beam Therapy of Uveal Melanoma: A Detailed Analysis of Tumor and Dose–Volume Parameters

Purpose: To determine neovascular glaucoma (NVG) incidence and identify contributing tumor and dosing factors in uveal melanoma patients treated with proton beam radiation therapy (PBRT). Methods and Materials: A total of 704 PBRT patients treated by a single surgeon (DHC) for uveal melanoma (1996-2010) were reviewed for NVG in our prospectively maintained database. All patients received 56 GyE in 4 fractions. Median follow-up was 58.3 months. Analyses included the Kaplan-Meier method to estimate NVG distributions, univariate log–rank tests, and Cox's proportional hazards multivariate analysis using likelihood ratio tests to identify independent risk factors of NVG among patient, tumor, and dose–volume histogram parameters. Results: The 5-year PBRT NVG rate was 12.7% (95% confidence interval [CI] 10.2%-15.9%). The 5-year rate of enucleation due to NVG was 4.9% (95% CI 3.4%-7.2%). Univariately, the NVG rate increased significantly with larger tumor diameter (P<.0001), greater height (P<.0001), higher T stage (P<.0001), and closer proximity to the disc (P=.002). Dose–volume histogram analysis revealed that if >30% of the lens or ciliary body received ≥50% dose (≥28 GyE), there was a higher probability of NVG (P<.0001 for both). Furthermore, if 100% of the disc or macula received ≥28 GyE, the NVG rate was higher (P<.0001 and P=.03, respectively). If both anterior and posterior doses were above specified cut points, NVG risk was highest (P<.0001). Multivariate analysis confirmed significant independent risk factors to include tumor height (P<.0001), age (P<.0001), %disc treated to ≥50% Dose (<100% vs 100%) (P=.0007), larger tumor diameter (P=.01), %lens treated to ≥90% Dose (0 vs >0%-30% vs >30%) (P=.01), and optic nerve length treated to ≥90% Dose (≤1 mm vs >1 mm) (P=.02). Conclusions: Our current PBRT patients experience a low rate of NVG and resultant enucleation compared with historical data. The present analysis shows that tumor height, diameter, and anterior as well as posterior critical structure dose–volume parameters may be used to predict NVG risk.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Effects of Respiratory Motion on Passively Scattered Proton Therapy Versus Intensity Modulated Photon Therapy for Stage III Lung Cancer: Are Proton Plans More Sensitive to Breathing Motion?

Purpose: To quantify and compare the effects of respiratory motion on paired passively scattered proton therapy (PSPT) and intensity modulated photon therapy (IMRT) plans; and to establish the relationship between the magnitude of tumor motion and the respiratory-induced dose difference for both modalities. Methods and Materials: In a randomized clinical trial comparing PSPT and IMRT, radiation therapy plans have been designed according to common planning protocols. Four-dimensional (4D) dose was computed for PSPT and IMRT plans for a patient cohort with respiratory motion ranging from 3 to 17 mm. Image registration and dose accumulation were performed using grayscale-based deformable image registration algorithms. The dose–volume histogram (DVH) differences (4D-3D [3D = 3-dimensional]) were compared for PSPT and IMRT. Changes in 4D-3D dose were correlated to the magnitude of tumor respiratory motion. Results: The average 4D-3D dose to 95% of the internal target volume was close to zero, with 19 of 20 patients within 1% of prescribed dose for both modalities. The mean 4D-3D between the 2 modalities was not statistically significant (P<.05) for all dose–volume histogram indices (mean ± SD) except the lung V5 (PSPT: +1.1% ± 0.9%; IMRT: +0.4% ± 1.2%) and maximum cord dose (PSPT: +1.5 ± 2.9 Gy; IMRT: 0.0 ± 0.2 Gy). Changes in 4D-3D dose were correlated to tumor motion for only 2 indices: dose to 95% planning target volume, and heterogeneity index. Conclusions: With our current margin formalisms, target coverage was maintained in the presence of respiratory motion up to 17 mm for both PSPT and IMRT. Only 2 of 11 4D-3D indices (lung V5 and spinal cord maximum) were statistically distinguishable between PSPT and IMRT, contrary to the notion that proton therapy will be more susceptible to respiratory motion. Because of the lack of strong correlations with 4D-3D dose differences in PSPT and IMRT, the extent of tumor motion was not an adequate predictor of potential dosimetric error caused by breathing motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Functional Data Analysis in NTCP Modeling: A New Method to Explore the Radiation Dose-Volume Effects

Purpose/Objective(s): To describe a novel method to explore radiation dose-volume effects. Functional data analysis is used to investigate the information contained in differential dose-volume histograms. The method is applied to the normal tissue complication probability modeling of rectal bleeding (RB) for patients irradiated in the prostatic bed by 3-dimensional conformal radiation therapy. Methods and Materials: Kernel density estimation was used to estimate the individual probability density functions from each of the 141 rectum differential dose-volume histograms. Functional principal component analysis was performed on the estimated probability density functions to explore the variation modes in the dose distribution. The functional principal components were then tested for association with RB using logistic regression adapted to functional covariates (FLR). For comparison, 3 other normal tissue complication probability models were considered: the Lyman-Kutcher-Burman model, logistic model based on standard dosimetric parameters (LM), and logistic model based on multivariate principal component analysis (PCA). Results: The incidence rate of grade ≥2 RB was 14%. V{sub 65Gy} was the most predictive factor for the LM (P=.058). The best fit for the Lyman-Kutcher-Burman model was obtained with n=0.12, m = 0.17, and TD50 = 72.6 Gy. In PCA and FLR, the components that describe the interdependence between the relative volumes exposed at intermediate and high doses were the most correlated to the complication. The FLR parameter function leads to a better understanding of the volume effect by including the treatment specificity in the delivered mechanistic information. For RB grade ≥2, patients with advanced age are significantly at risk (odds ratio, 1.123; 95% confidence interval, 1.03-1.22), and the fits of the LM, PCA, and functional principal component analysis models are significantly improved by including this clinical factor. Conclusion: Functional data analysis provides an attractive method for flexibly estimating the dose-volume effect for normal tissues in external radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Impact of Fractionation and Dose in a Multivariate Model for Radiation-Induced Chest Wall Pain

Purpose: To determine the role of patient/tumor characteristics, radiation dose, and fractionation using the linear-quadratic (LQ) model to predict stereotactic body radiation therapy–induced grade ≥2 chest wall pain (CWP2) in a larger series and develop clinically useful constraints for patients treated with different fraction numbers. Methods and Materials: A total of 316 lung tumors in 295 patients were treated with stereotactic body radiation therapy in 3 to 5 fractions to 39 to 60 Gy. Absolute dose–absolute volume chest wall (CW) histograms were acquired. The raw dose-volume histograms (α/β = ∞ Gy) were converted via the LQ model to equivalent doses in 2-Gy fractions (normalized total dose, NTD) with α/β from 0 to 25 Gy in 0.1-Gy steps. The Cox proportional hazards (CPH) model was used in univariate and multivariate models to identify and assess CWP2 exposed to a given physical and NTD. Results: The median follow-up was 15.4 months, and the median time to development of CWP2 was 7.4 months. On a univariate CPH model, prescription dose, prescription dose per fraction, number of fractions, D83cc, distance of tumor to CW, and body mass index were all statistically significant for the development of CWP2. Linear-quadratic correction improved the CPH model significance over the physical dose. The best-fit α/β was 2.1 Gy, and the physical dose (α/β = ∞ Gy) was outside the upper 95% confidence limit. With α/β = 2.1 Gy, V{sub NTD99Gy} was most significant, with median V{sub NTD99Gy} = 31.5 cm{sup 3} (hazard ratio 3.87, P<.001). Conclusion: There were several predictive factors for the development of CWP2. The LQ-adjusted doses using the best-fit α/β = 2.1 Gy is a better predictor of CWP2 than the physical dose. To aid dosimetrists, we have calculated the physical dose equivalent corresponding to V{sub NTD99Gy} = 31.5 cm{sup 3} for the 3- to 5-fraction groups.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Exploratory Study of 4D versus 3D Robust Optimization in Intensity Modulated Proton Therapy for Lung Cancer

Purpose: The purpose of this study was to compare the impact of uncertainties and interplay on 3-dimensional (3D) and 4D robustly optimized intensity modulated proton therapy (IMPT) plans for lung cancer in an exploratory methodology study. Methods and Materials: IMPT plans were created for 11 nonrandomly selected non-small cell lung cancer (NSCLC) cases: 3D robustly optimized plans on average CTs with internal gross tumor volume density overridden to irradiate internal target volume, and 4D robustly optimized plans on 4D computed tomography (CT) to irradiate clinical target volume (CTV). Regular fractionation (66 Gy [relative biological effectiveness; RBE] in 33 fractions) was considered. In 4D optimization, the CTV of individual phases received nonuniform doses to achieve a uniform cumulative dose. The root-mean-square dose-volume histograms (RVH) measured the sensitivity of the dose to uncertainties, and the areas under the RVH curve (AUCs) were used to evaluate plan robustness. Dose evaluation software modeled time-dependent spot delivery to incorporate interplay effect with randomized starting phases of each field per fraction. Dose-volume histogram (DVH) indices comparing CTV coverage, homogeneity, and normal tissue sparing were evaluated using Wilcoxon signed rank test. Results: 4D robust optimization plans led to smaller AUC for CTV (14.26 vs 18.61, respectively; P=.001), better CTV coverage (Gy [RBE]) (D{sub 95%} CTV: 60.6 vs 55.2, respectively; P=.001), and better CTV homogeneity (D{sub 5%}-D{sub 95%} CTV: 10.3 vs 17.7, resspectively; P=.002) in the face of uncertainties. With interplay effect considered, 4D robust optimization produced plans with better target coverage (D{sub 95%} CTV: 64.5 vs 63.8, respectively; P=.0068), comparable target homogeneity, and comparable normal tissue protection. The benefits from 4D robust optimization were most obvious for the 2 typical stage III lung cancer patients. Conclusions: Our exploratory methodology study showed that, compared to 3D robust optimization, 4D robust optimization produced significantly more robust and interplay-effect-resistant plans for targets with comparable dose distributions for normal tissues. A further study with a larger and more realistic patient population is warranted to generalize the conclusions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗