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

Potential benefits of dosimetric VMAT tracking verified with 3D film measurements

Purpose: To evaluate three different plan adaptation strategies using 3D film-stack dose measurements of both focal boost and hypofractionated prostate VMAT treatments. The adaptation strategies (a couch shift, geometric tracking, and dosimetric tracking) were applied for three realistic intrafraction prostate motions. Methods: A focal boost (35 × 2.2 and 35 × 2.7 Gy) and a hypofractionated (5 × 7.25 Gy) prostate VMAT plan were created for a heterogeneous phantom that allows for internal prostate motion. For these plans geometric tracking and dosimetric tracking were evaluated by ionization chamber (IC) point dose measurements (zero-D) and measurements using a stack of EBT3 films (3D). The geometric tracking applied translations, rotations, and scaling of the MLC aperture in response to realistic prostate motions. The dosimetric tracking additionally corrected the monitor units to resolve variations due to difference in depth, tissue heterogeneity, and MLC-aperture. The tracking was based on the positions of four fiducial points only. The film measurements were compared to the gold standard (i.e., IC measurements) and the planned dose distribution. Additionally, the 3D measurements were converted to dose volume histograms, tumor control probability, and normal tissue complication probability parameters (DVH/TCP/NTCP) as a direct estimate of clinical relevance of the proposed tracking. Results: Compared to the planned dose distribution, measurements without prostate motion and tracking showed already a reduced homogeneity of the dose distribution. Adding prostate motion further blurs the DVHs for all treatment approaches. The clinical practice (no tracking) delivered the dose distribution inside the PTV but off target (CTV), resulting in boost dose errors up to 10%. The geometric and dosimetric tracking corrected the dose distribution’s position. Moreover, the dosimetric tracking could achieve the planned boost DVH, but not the DVH of the more homogeneously irradiated prostate. A drawback of both the geometric and dosimetric tracking was a reduced MLC blocking caused by the rotational component of the MLC aperture corrections. Because of the used CTV to PTV margins and the high doses in the considered fractionation schemes, the TCP differed less than 0.02 from the planned value for all targets and all correction methods. The rectal NTCP constraints, however, could not be realized using any of these methods. Conclusions: The geometric and dosimetric tracking use only a limited input, but they deposit the dose distribution with higher geometric accuracy than the clinical practice. The latter case has boost dose errors up to 10%. The increased accuracy has a modest impact [Δ(NT)CP < 0.02] because of the applied margins and the high dose levels used. To allow further margin reduction tracking methods are vital. The proposed methodology could further be improved by implementing a rotational correction using collimator rotations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Technical Note: Statistical dependences between channels in radiochromic film readings. Implications in multichannel dosimetry

Purpose: This note studies the statistical relationships between color channels in radiochromic film readings with flatbed scanners. The same relationships are studied for noise. Finally, their implications for multichannel film dosimetry are discussed. Methods: Radiochromic films exposed to wedged fields of 6 MV energy were read in a flatbed scanner. The joint histograms of pairs of color channels were used to obtain the joint and conditional probability density functions between channels. Then, the conditional expectations and variances of one channel given another channel were obtained. Noise was extracted from film readings by means of a multiresolution analysis. Two different dose ranges were analyzed, the first one ranging from 112 to 473 cGy and the second one from 52 to 1290 cGy. Results: For the smallest dose range, the conditional expectations of one channel given another channel can be approximated by linear functions, while the conditional variances are fairly constant. The slopes of the linear relationships between channels can be used to simplify the expression that estimates the dose by means of the multichannel method. The slopes of the linear relationships between each channel and the red one can also be interpreted as weights in the final contribution to dose estimation. However, for the largest dose range, the conditional expectations of one channel given another channel are no longer linear functions. Finally, noises in different channels were found to correlate weakly. Conclusions: Signals present in different channels of radiochromic film readings show a strong statistical dependence. By contrast, noise correlates weakly between channels. For the smallest dose range analyzed, the linear behavior between the conditional expectation of one channel given another channel can be used to simplify calculations in multichannel film dosimetry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Nanoscale measurements of proton tracks using fluorescent nuclear track detectors

Purpose: The authors describe a method in which fluorescence nuclear track detectors (FNTDs), novel track detectors with nanoscale spatial resolution, are used to determine the linear energy transfer (LET) of individual proton tracks from proton therapy beams by allowing visualization and 3D reconstruction of such tracks. Methods: FNTDs were exposed to proton therapy beams with nominal energies ranging from 100 to 250 MeV. Proton track images were then recorded by confocal microscopy of the FNTDs. Proton tracks in the FNTD images were fit by using a Gaussian function to extract fluorescence amplitudes. Histograms of fluorescence amplitudes were then compared with LET spectra. Results: The authors successfully used FNTDs to register individual proton tracks from high-energy proton therapy beams, allowing reconstruction of 3D images of proton tracks along with delta rays. The track amplitudes from FNTDs could be used to parameterize LET spectra, allowing the LET of individual proton tracks from therapeutic proton beams to be determined. Conclusions: FNTDs can be used to directly visualize proton tracks and their delta rays at the nanoscale level. Because the track intensities in the FNTDs correlate with LET, they could be used further to measure LET of individual proton tracks. This method may be useful for measuring nanoscale radiation quantities and for measuring the LET of individual proton tracks in radiation biology experiments.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology

Purpose: Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE. Methods: Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of {sup 99m}Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a “liver’s segment,” or distributing multiple uptaking nodules (“nonuniform liver”), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and {sup 99m}Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with {sup 90}Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the “gold standard,” GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again. Results: Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%–14% for a liver segment, and by 7%–12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from −10%/40% for a liver segment, and −10%/20% for a nonuniform liver, to −13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, −38%/−13% in liver’s segment, −61%/−40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: −7%/13%). Conclusions: Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Four-dimensional dose reconstruction through in vivo phase matching of cine images of electronic portal imaging device

Purpose: A method is proposed to reconstruct a four-dimensional (4D) dose distribution using phase matching of measured cine images to precalculated images of electronic portal imaging device (EPID). Methods: (1) A phantom, designed to simulate a tumor in lung (a polystyrene block with a 3 cm diameter embedded in cork), was placed on a sinusoidally moving platform with an amplitude of 1 cm and a period of 4 s. Ten-phase 4D computed tomography (CT) images of the phantom were acquired. A planning target volume (PTV) was created by adding a margin of 1 cm around the internal target volume of the tumor. (2) Three beams were designed, which included a static beam, a theoretical dynamic beam, and a planning-optimized dynamic beam (PODB). While the theoretical beam was made by manually programming a simplistic sliding leaf motion, the planning-optimized beam was obtained from treatment planning. From the three beams, three-dimensional (3D) doses on the phantom were calculated; 4D dose was calculated by means of the ten phase images (integrated over phases afterward); serving as “reference” images, phase-specific EPID dose images under the lung phantom were also calculated for each of the ten phases. (3) Cine EPID images were acquired while the beams were irradiated to the moving phantom. (4) Each cine image was phase-matched to a phase-specific CT image at which common irradiation occurred by intercomparing the cine image with the reference images. (5) Each cine image was used to reconstruct dose in the phase-matched CT image, and the reconstructed doses were summed over all phases. (6) The summation was compared with forwardly calculated 4D and 3D dose distributions. Accounting for realistic situations, intratreatment breathing irregularity was simulated by assuming an amplitude of 0.5 cm for the phantom during a portion of breathing trace in which the phase matching could not be performed. Intertreatment breathing irregularity between the time of treatment and the time of planning CT was considered by utilizing the same reduced amplitude when the phantom was irradiated. To examine the phase matching in a humanoid environment, the matching was also performed in a digital phantom (4D XCAT phantom). Results: For the static, the theoretical, and the planning-optimized dynamic beams, the 4D reconstructed doses showed agreement with the forwardly calculated 4D doses within the gamma pass rates of 92.7%, 100%, and 98.1%, respectively, at the isocenter plane given by 3%/3 mm criteria. Excellent agreement in dose volume histogram of PTV and lung-PTV was also found between the two 4D doses, while substantial differences were found between the 3D and the 4D doses. The significant breathing irregularities modeled in this study were found not to be noticeably affecting the reconstructed dose. The phase matching was performed equally well in a digital phantom. Conclusions: The method of retrospective phase determination of a moving object under irradiation provided successful 4D dose reconstruction. This method will provide accurate quality assurance and facilitate adaptive therapy when distinguishable objects such as well-defined tumors, diaphragm, and organs with markers (pancreas and liver) are covered by treatment beam apertures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Generic and robust method for automatic segmentation of PET images using an active contour model

Purpose: Although positron emission tomography (PET) images have shown potential to improve the accuracy of targeting in radiation therapy planning and assessment of response to treatment, the boundaries of tumors are not easily distinguishable from surrounding normal tissue owing to the low spatial resolution and inherent noisy characteristics of PET images. The objective of this study is to develop a generic and robust method for automatic delineation of tumor volumes using an active contour model and to evaluate its performance using phantom and clinical studies. Methods: MASAC, a method for automatic segmentation using an active contour model, incorporates the histogram fuzzy C-means clustering, and localized and textural information to constrain the active contour to detect boundaries in an accurate and robust manner. Moreover, the lattice Boltzmann method is used as an alternative approach for solving the level set equation to make it faster and suitable for parallel programming. Twenty simulated phantom studies and 16 clinical studies, including six cases of pharyngolaryngeal squamous cell carcinoma and ten cases of nonsmall cell lung cancer, were included to evaluate its performance. Besides, the proposed method was also compared with the contourlet-based active contour algorithm (CAC) and Schaefer’s thresholding method (ST). The relative volume error (RE), Dice similarity coefficient (DSC), and classification error (CE) metrics were used to analyze the results quantitatively. Results: For the simulated phantom studies (PSs), MASAC and CAC provide similar segmentations of the different lesions, while ST fails to achieve reliable results. For the clinical datasets (2 cases with connected high-uptake regions excluded) (CSs), CAC provides for the lowest mean RE (−8.38% ± 27.49%), while MASAC achieves the best mean DSC (0.71 ± 0.09) and mean CE (53.92% ± 12.65%), respectively. MASAC could reliably quantify different types of lesions assessed in this work with good accuracy, resulting in a mean RE of −13.35% ± 11.87% and −11.15% ± 23.66%, a mean DSC of 0.89 ± 0.05 and 0.71 ± 0.09, and a mean CE of 19.19% ± 7.89% and 53.92% ± 12.65%, for PSs and CSs, respectively. Conclusions: The authors’ results demonstrate that the developed novel PET segmentation algorithm is applicable to various types of lesions in the authors’ study and is capable of producing accurate and consistent target volume delineations, potentially resulting in reduced intraobserver and interobserver variabilities observed when using manual delineation and improved accuracy in treatment planning and outcome evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Backscatter dose effects for high atomic number materials being irradiated in the presence of a magnetic field: A Monte Carlo study for the MRI linac

Purpose: To quantify and explain the backscatter dose effects for clinically relevant high atomic number materials being irradiated in the presence of a 1.5 T transverse magnetic field. Methods: Interface effects were investigated using Monte Carlo simulation techniques. We used GPUMCD (v5.1) and GEANT4 (v10.1) for this purpose. GPUMCD is a commercial software written for the Elekta AB, MRI linac. Dose was scored using GPUMCD in cubic voxels of side 1 and 0.5 mm, in two different virtual phantoms of dimensions 20 × 20 × 20 cm and 5 × 5 × 13.3 cm, respectively. A photon beam was generated from a point 143.5 cm away from the isocenter with energy distribution sampled from a histogram representing the true Elekta, MRI linac photon spectrum. A slab of variable thickness and position containing either bone, aluminum, titanium, stainless steel, or one of the two different dental filling materials was inserted as an inhomogeneity in the 20 × 20 × 20 cm phantom. The 5 × 5 × 13.3 cm phantom was used as a clinical test case in order to explain the dose perturbation effects for a head and neck cancer patient. The back scatter dose factor (BSDF) was defined as the ratio of the doses at a given depth with and without the presence of the inhomogeneity. Backscattered electron fluence was calculated at the inhomogeneity interface using GEANT4. A 1.5 T magnetic field was applied perpendicular to the direction of the beam in both phantoms, identical to the geometry in the Elekta MRI linac. Results: With the application of a 1.5 T magnetic field, all the BSDF’s were reduced by 12%–47%, compared to the no magnetic field case. The corresponding backscattered electron fluence at the interface was also reduced by 45%–64%. The reduction in the BSDF at the interface, due to the application of the magnetic field, is manifested in a different manner for each material. In the case of bone, the dose drops at the interface contrary to the expected increase when no magnetic field is applied. In the case of aluminum, the dose at the interface is the same with and without the presence of the aluminum. For all of the other materials the dose increases at the interface. Conclusions: The reduction in dose at the interface, in the presence of the magnetic field, is directly related to the reduction in backscattered electron fluence. This reduction occurs due to two different reasons. First, the electron spectrum hitting the interface is changed when the magnetic field is turned on, which results in changes in the electron scattering probability. Second, some electrons that have curved trajectories due to the presence of the magnetic field are absorbed by the higher density side of the interface and no longer contribute to the backscattered electron fluence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Modification and validation of an analytical source model for external beam radiotherapy Monte Carlo dose calculations

Purpose: A dose calculation tool, which combines the accuracy of the dose planning method (DPM) Monte Carlo code and the versatility of a practical analytical multisource model, which was previously reported has been improved and validated for the Varian 6 and 10 MV linear accelerators (linacs). The calculation tool can be used to calculate doses in advanced clinical application studies. One shortcoming of current clinical trials that report dose from patient plans is the lack of a standardized dose calculation methodology. Because commercial treatment planning systems (TPSs) have their own dose calculation algorithms and the clinical trial participant who uses these systems is responsible for commissioning the beam model, variation exists in the reported calculated dose distributions. Today’s modern linac is manufactured to tight specifications so that variability within a linac model is quite low. The expectation is that a single dose calculation tool for a specific linac model can be used to accurately recalculate dose from patient plans that have been submitted to the clinical trial community from any institution. The calculation tool would provide for a more meaningful outcome analysis. Methods: The analytical source model was described by a primary point source, a secondary extra-focal source, and a contaminant electron source. Off-axis energy softening and fluence effects were also included. The additions of hyperbolic functions have been incorporated into the model to correct for the changes in output and in electron contamination with field size. A multileaf collimator (MLC) model is included to facilitate phantom and patient dose calculations. An offset to the MLC leaf positions was used to correct for the rudimentary assumed primary point source. Results: Dose calculations of the depth dose and profiles for field sizes 4 × 4 to 40 × 40 cm agree with measurement within 2% of the maximum dose or 2 mm distance to agreement (DTA) for 95% of the data points tested. The model was capable of predicting the depth of the maximum dose within 1 mm. Anthropomorphic phantom benchmark testing of modulated and patterned MLCs treatment plans showed agreement to measurement within 3% in target regions using thermoluminescent dosimeters (TLD). Using radiochromic film normalized to TLD, a gamma criteria of 3% of maximum dose and 2 mm DTA was applied with a pass rate of least 85% in the high dose, high gradient, and low dose regions. Finally, recalculations of patient plans using DPM showed good agreement relative to a commercial TPS when comparing dose volume histograms and 2D dose distributions. Conclusions: A unique analytical source model coupled to the dose planning method Monte Carlo dose calculation code has been modified and validated using basic beam data and anthropomorphic phantom measurement. While this tool can be applied in general use for a particular linac model, specifically it was developed to provide a singular methodology to independently assess treatment plan dose distributions from those clinical institutions participating in National Cancer Institute trials.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-C-202-05: Pilot Study of Online Treatment Evaluation and Adaptive Re-Planning for Laryngeal SBRT

Purpose: We have instigated a phase I trial of 5-fraction stereotactic body radiotherapy (SBRT) for advanced-stage laryngeal cancer. We conducted this pilot dosimetric study to confirm the potential utility of online adaptive re-planning to preserve treatment quality. Methods: Ten cases of larynx cancer were evaluated. Baseline and daily SBRT treatment plans were generated per trial protocol. Daily volumetric images were acquired prior to every fraction of treatment. Reference simulation CT images were deformably registered to daily volumetric images using Eclipse. Planning contours were then deformably propagated to daily images. Reference SBRT plans were directly copied to calculate delivered dose distributions on deformed reference CT images. In-house software platform has been developed to calculate cumulative dose over a course of treatment in four steps: 1) deforming delivered dose grid to reference CT images using deformation information exported from Eclipse; 2) generating tetrahedrons using deformed dose grid as vertices; 3) resampling dose to a high resolution within every tetrahedron; 4) calculating dose-volume histograms. Our inhouse software was benchmarked with a commercial software, Mirada. Results: In all ten cases including 49 fractions of treatments, delivered daily doses were completely evaluated and treatment could be re-planned within 10 minutes. Prescription dose coverage of PTV was less than intended in 53% of fractions of treatment (mean: 94%, range: 84%–98%) while minimum coverage of CTV and GTV was 94% and 97%, respectively. Maximum bystander point dose limits to arytenoids, parotids, and spinal cord remained respected in all cases, although variances in carotid artery doses were observed in a minority of cases. Conclusion: Although GTV and CTV coverage is preserved by in-room 3D image guidance of larynx SBRT, PTV coverage can vary significantly from intended plans. Online adaptive treatment evaluation and re-planning is potentially necessary and our procedure is clinically applicable to fully preserve treatment quality. This project is supported by CPRIT Individual Investigator Research Award RP150386.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-D-207B-07: Development of a CT-Radiomics Based Early Response Prediction Model During Delivery of Chemoradiation Therapy for Pancreatic Cancer

Purpose: Pilot study developing a CT-texture based model for early assessment of treatment response during the delivery of chemoradiation therapy (CRT) for pancreatic cancer. Methods: Daily CT data acquired for 24 pancreatic head cancer patients using CT-on-rails, during the routine CT-guided CRT delivery with a radiation dose of 50.4 Gy in 28 fractions, were analyzed. The pancreas head was contoured on each daily CT. Texture analysis was performed within the pancreas head contour using a research tool (IBEX). Over 1300 texture metrics including: grey level co-occurrence, run-length, histogram, neighborhood intensity difference, and geometrical shape features were calculated for each daily CT. Metric-trend information was established by finding the best fit of either a linear, quadratic, or exponential function for each metric value verses accumulated dose. Thus all the daily CT texture information was consolidated into a best-fit trend type for a given patient and texture metric. Linear correlation was performed between the patient histological response vector (good, medium, poor) and all combinations of 23 patient subgroups (statistical jackknife) determining which metrics were most correlated to response and repeatedly reliable across most patients. Control correlations against CT scanner, reconstruction kernel, and gated/nongated CT images were also calculated. Euclidean distance measure was used to group/sort patient vectors based on the data of these trend-response metrics. Results: We found four specific trend-metrics (Gray Level Coocurence Matrix311-1InverseDiffMomentNorm, Gray Level Coocurence Matrix311-1InverseDiffNorm, Gray Level Coocurence Matrix311-1 Homogeneity2, and Intensity Direct Local StdMean) that were highly correlated with patient response and repeatedly reliable. Our four trend-metric model successfully ordered our pilot response dataset (p=0.00070). We found no significant correlation to our control parameters: gating (p=0.7717), scanner (p=0.9741), and kernel (p=0.8586). Conclusion: We have successfully created a CT-texture based early treatment response prediction model using the CTs acquired during the delivery of chemoradiation therapy for pancreatic cancer. Future testing is required to validate the model with more patient data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-P-35: A Multi-Institutional Plan Quality Checking Tool Built On Oncospace: A Shared Radiation Oncology Database System

Purpose: Late toxicity from radiation to critical structures limits the possible dose in Radiation Therapy. Perfectly conformal treatment of a target is not realizable, so the clinician must accept a certain level of collateral radiation to nearby OARs. But how much? General guidelines exist for healthy tissue sparing which guide RT treatment planning, but are these guidelines good enough to create the optimal plan given the individualized patient anatomy? We propose a means to evaluate the planned dose level to an OAR using a multi-institutional data-store of previously treated patients, so a clinician might reconsider planning objectives. Methods: The tool is built on Oncospace, a federated data-store system, which consists of planning data import, web based analysis tools, and a database containing:1) DVHs: dose by percent volume delivered to each ROI for each patient previously treated and included in the database.2) Overlap Volume Histograms (OVHs): Anatomical measure defined as the percent volume of an ROI within a given distance to target structures.Clinicians know what OARs are important to spare. For any ROI, Oncospace knows for which patients’ anatomy that ROI was harder to plan in the past (the OVH is less). The planned dose should be close to the least dose of previous patients. The tool displays the dose those OARs were subjected to, and the clinician can make a determination about the planning objectives used.Multiple institutions contribute to the Oncospace Consortium, and their DVH and OVH data are combined and color coded in the output. Results: The Oncospace website provides a plan quality display tool which identifies harder to treat patients, and graphically displays the dose delivered to them for comparison with the proposed plan. Conclusion: The Oncospace Consortium manages a data-store of previously treated patients which can be used for quality checking new plans. Grant funding by Elekta.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-P-54: Guidelines to Check Image Registration QA of a Clinical Deformation Registration Software: A Single Institution Preliminary Study

Purpose: The objective of this study is to verify and analyze the accuracy of a clinical deformable image registration (DIR) software. Methods: To test clinical DIR software qualitatively and quantitatively, we focused on lung radiotherapy and analyzed a single (Lung) patient CT scan. Artificial anatomical changes were applied to account for daily variations during the course of treatment including the planning target volume (PTV) and organs at risk (OAR). The primary CT (pCT) and the structure set (pST) was deformed with commercial tool (ImSimQA-Oncology Systems Limited) and after artificial deformation (dCT and dST) sent to another commercial tool (VelocityAI-Varian Medical Systems). In Velocity, the deformed CT and structures (dCT and dST) were inversely deformed back to original primary CT (dbpCT and dbpST). We compared the dbpST and pST structure sets using similarity metrics. Furthermore, a binary deformation field vector (BDF) was created and sent to ImSimQA software for comparison with known “ground truth” deformation vector fields (DVF). Results: An image similarity comparison was made by using “ground truth” DVF and “deformed output” BDF with an output of normalized “cross correlation (CC)” and “mutual information (MI)” in ImSimQA software. Results for the lung case were MI=0.66 and CC=0.99. The artificial structure deformation in both pST and dbpST was analyzed using DICE coefficient, mean distance to conformity (MDC) and deformation field error volume histogram (DFEVH) by comparing them before and after inverse deformation. We have noticed inadequate structure match for CTV, ITV and PTV due to close proximity of heart and overall affected by lung expansion. Conclusion: We have seen similarity between pCT and dbpCT but not so well between pST and dbpST, because of inadequate structure deformation in clinical DIR system. This system based quality assurance test will prepare us for adopting the guidelines of upcoming AAPM task group 132 protocol.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-18: Updates to the Computational Environment for Radiological Research for Image Analysis

Purpose: To present new tools in CERR for Texture Analysis and Visualization. Method: (1) Quantitative Image Analysis: We added the ability to compute Haralick texture features based on local neighbourhood. The Texture features depend on many parameters used in their derivation. For example: (a) directionality, (b) quantization of image, (c) patch-size for the neighborhood, (d) handling of the edge voxels within the region of interest, (e) Averaging co-occurance matrix vs texture features for different directions etc. A graphical user interface was built to set these parameters and then visualize their impact on the resulting texture maps. The entire functionality was written in Matlab. Array indexing was used to speed up the texture calculation. The computation speed is very competitive with the ITK library. Moreover, our implementation works with multiple CPUs and the computation time can be further reduced by using multiple processor threads. In order to reduce the Haralick texture maps into scalar features, we propose the use of Texture Volume Histograms. This lets users make use of the entire distribution of texture values within the region of interest rather than using just the mean and the standard deviations. (2) Qualitative/Visualization tools: The derived texture maps are stored as a new scan (derived) within CERR’s planC data structure. A display that compares various scans was built to show the raw image and the derived texture maps side-by-side. These images are positionally linked and can be navigated together. CERR’s graphics handling was updated and sped-up to be compatible with the newer Matlab versions. As a result, the users can use (a) different window levels and colormaps for different viewports, (b) click-and-drag or use mouse scroll-wheel to navigate slices. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics and Outcomes modeling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-24: Identifying Prognostic Imaging Biomarkers in Early Stage Lung Cancer Using Radiomics

Purpose: Patients diagnosed with early stage lung cancer have favorable outcomes when treated with surgery or stereotactic radiotherapy. However, a significant proportion (∼20%) of patients will develop metastatic disease and eventually die of the disease. The purpose of this work is to identify quantitative imaging biomarkers from CT for predicting overall survival in early stage lung cancer. Methods: In this institutional review board-approved HIPPA-compliant retrospective study, we retrospectively analyzed the diagnostic CT scans of 110 patients with early stage lung cancer. Data from 70 patients were used for training/discovery purposes, while those of remaining 40 patients were used for independent validation. We extracted 191 radiomic features, including statistical, histogram, morphological, and texture features. Cox proportional hazard regression model, coupled with the least absolute shrinkage and selection operator (LASSO), was used to predict overall survival based on the radiomic features. Results: The optimal prognostic model included three image features from the Law’s feature and wavelet texture. In the discovery cohort, this model achieved a concordance index or CI=0.67, and it separated the low-risk from high-risk groups in predicting overall survival (hazard ratio=2.72, log-rank p=0.007). In the independent validation cohort, this radiomic signature achieved a CI=0.62, and significantly stratified the low-risk and high-risk groups in terms of overall survival (hazard ratio=2.20, log-rank p=0.042). Conclusion: We identified CT imaging characteristics associated with overall survival in early stage lung cancer. If prospectively validated, this could potentially help identify high-risk patients who might benefit from adjuvant systemic therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-35: Repeatability of Texture Features in T1- and T2-Weighted MR Images

Purpose: To evaluate repeatability of lung tumor texture features from inspiration/expiration MR image pairs for potential use in patient specific care models and applications. Repeatability is a desirable and necessary characteristic of features included in such models. Methods: T1-weighted Volumetric Interpolation Breath-Hold Examination (VIBE) and/or T2-weighted MRI scans were acquired for 15 patients with non-small cell lung cancer before and during radiotherapy for a total of 32 and 34 same session inspiration-expiration breath-hold image pairs respectively. Bias correction was applied to the VIBE (VIBE-BC) and T2-weighted (T2-BC) images. Fifty-nine texture features at five wavelet decomposition ratios were extracted from the delineated primary tumor including: histogram(HIST), gray level co-occurrence matrix(GLCM), gray level run length matrix(GLRLM), gray level size zone matrix(GLSZM), and neighborhood gray tone different matrix (NGTDM) based features. Repeatability of the texture features for VIBE, VIBE-BC, T2-weighted, and T2-BC image pairs was evaluated by the concordance correlation coefficient (CCC) between corresponding image pairs, with a value greater than 0.90 indicating repeatability. Results: For the VIBE image pairs, the percentage of repeatable texture features by wavelet ratio was between 20% and 24% of the 59 extracted features; the T2-weighted image pairs exhibited repeatability in the range of 44–49%. The percentage dropped to 10–20% for the VIBE-BC images, and 12–14% for the T2-BC images. In addition, five texture features were found to be repeatable in all four image sets including two GLRLM, two GLZSM, and one NGTDN features. No single texture feature category was repeatable among all three image types; however, certain categories performed more consistently on a per image type basis. Conclusion: We identified repeatable texture features on T1- and T2-weighted MRI scans. These texture features should be further investigated for use in specific applications such as tissue classification and changes during radiation therapy utilizing a standard imaging protocol. Authors have the following disclosures: a research agreement with Philips Medical systems (Hugo, Weiss), a license agreement with Varian Medical Systems (Hugo, Weiss), research grants from the National Institute of Health (Hugo, Weiss), UpToDate royalties (Weiss), and none(Mahon, Ford, Karki). Authors have no potential conflicts of interest to disclose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-36: Validating Quantitative Radiomic Texture Features for Oncologic PET: A Digital Phantom Study

Purpose: Radiomic texture features derived from the oncologic PET have recently been brought under intense investigation within the context of patient stratification and treatment outcome prediction in a variety of cancer types; however, their validity has not yet been examined. This work is aimed to validate radiomic PET texture metrics through the use of realistic simulations in the ground truth setting. Methods: Simulation of FDG-PET was conducted by applying the Zubal phantom as an attenuation map to the SimSET software package that employs Monte Carlo techniques to model the physical process of emission imaging. A total of 15 irregularly-shaped lesions featuring heterogeneous activity distribution were simulated. For each simulated lesion, 28 texture features in relation to the intensity histograms (GLIH), grey-level co-occurrence matrices (GLCOM), neighborhood difference matrices (GLNDM), and zone size matrices (GLZSM) were evaluated and compared with their respective values extracted from the ground truth activity map. Results: In reference to the values from the ground truth images, texture parameters appearing on the simulated data varied with a range of 0.73–3026.2% for GLIH-based, 0.02–100.1% for GLCOM-based, 1.11–173.8% for GLNDM-based, and 0.35–66.3% for GLZSM-based. For majority of the examined texture metrics (16/28), their values on the simulated data differed significantly from those from the ground truth images (P-value ranges from <0.0001 to 0.04). Features not exhibiting significant difference comprised of GLIH-based standard deviation, GLCO-based energy and entropy, GLND-based coarseness and contrast, and GLZS-based low gray-level zone emphasis, high gray-level zone emphasis, short zone low gray-level emphasis, long zone low gray-level emphasis, long zone high gray-level emphasis, and zone size nonuniformity. Conclusion: The extent to which PET imaging disturbs texture appearance is feature-dependent and could be substantial. It is thus advised that use of PET texture parameters for predictive and prognostic measurements in oncologic setting awaits further systematic and critical evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-47: Quantitative Shape Relationship Analysis of PTV Modification for Critical Anatomy Sparing and Its Impact On Pathologic Response for Neoadjuvant Stereotactic Radiotherapy for Pancreatic Cancer

Purpose: Stereotactic body radiation therapy (SBRT) may be used to increase surgery candidacy in borderline resectable (BRPC) and locally advanced (LAPC) pancreatic cancer. However, the planning target volume (PTV) may need to be limited to avoid toxicity when the gross tumor volume (GTV) is anatomically involved with surrounding critical structures. Our study aims to characterize the coverage of GTV and investigate the association between modified PTV and pathologic (pCR) or near pathologic (npCR) complete response rates determined from the surgical specimen. Methods: Patients treated with neoadjuvant pancreas SBRT followed by surgery from 2010–2015 were selected from Oncospace. Overlap volume histogram (OVH) analysis was performed to determine the extent of compromise of the PTV from both the GTV and a standard target (GTV+3mm). Subsequently, normalized overlap volume (%) was calculated for: (1) GTV-PTV, and (2) GTV+3mm expansion-PTV. A logistic regression model was used to identify the association between the overlap ratios and ≥ npCR(pCR/npCR) stratified by active breathing control (ABC) versus free-breathing status. Results: Eighty-one (BRPC: n=42, LAPC: n=39) patients were available for analysis. Nearly 40% (31/81) had ≥npCR and 75% (61/81) were able to complete ABC. Mean coverage of the GTV-PTV was 92.6% (range, 59.9%–100%, SD = 8.68) and coverage of the GTV+3mm expansion-PTV was 85. 2% (range, 59.9% −100.0%, SD= 8.67). Among the patients with ABC, every 10% increase in GTV coverage doubled the odds to have ≥npCR (OR = 1.82, p=0.06). Coverage of GTV+3mm expansion was not associated with ≥npCR regardless of ABC status. Conclusion: Preferential sparing of critical anatomy over GTV-PTV coverage with ABC management suggests worse ≥npCR rates for neoadjuvant SBRT in BRPC and LAPC. Limiting the GTV and GTV+3mm expansion in free-breathing patients was not associated with pathologic response perhaps due to larger GTV definitions as a result of motion artifacts in free-breathing CT scans. Collaboration with Toshiba, Elekta, and Phillips.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-R-54: CT-Texture Based Early Tumor Treatment Response Assessment During Radiation Therapy Delivery: Small Cell Versus Non-Small Cell Lung Cancers

Purpose: Tumor treatment response may potentially be assessed during radiation therapy (RT) by analyzing changes in CT-textures. We investigated the different early RT-responses between small cell (SCLC) and non-small cell lung cancer (NSCLC) as assessed by CT-texture. Methods: Daily diagnostic-quality CT acquired during routine CT-guided RT using a CT-on-Rails for 13-NSCLC and 5-SCLC patients were analyzed. These patient had ages ranging from 45–78 and 38–63 years, respectively, for NSCLC and SCLC groups, and tumor-stages ranging from T2-T4, and were treated with either RT or chemotherapy and RT with 45–66Gy/ 20–34 fractions. Gross-tumor volume (GTV) contour was generated on each daily CT by populating GTV contour from simulation to daily CTs with manual editing if necessary. CT-texture parameters, such as Hounsfield Unit (HU) histogram, mean HU, skewness, kurtosis, entropy, and short-run high-gray level emphasis (SRHGLE), were calculated in GTV from each daily CT-set using an in house software tool. Difference in changes of these texture parameters during RT between NSCLC and SCLC was analyzed and compared with GTV volume changes. Results: Radiation-induced changes in CT-texture were different between SCLC and NSCLC. Average changes from first to the last fractions for NSCLC and SCLC in GTV were 28±10(12–44) and 30±15(11–47) HU (mean HU reduction), 12.7% and 18.3% (entropy), 50% and 55% (SRHGLE), 19% and 22% (kurtosis), and 5.2% and 22% (skewness), respectively. Good correlation in kurtosis changes and GTV was seen (R{sup 2}=0.8923) for SCLC, but not for NSCLC (R{sup 2}=0.4748). SCLC had better correlations between GTV volume reduction and entropy (SCLC R{sup 2}=0.847; NSCLC R{sup 2}=0.6485), skewness (SCLC R{sup 2}=0.935; NSCLC R{sup 2}=0.7666), or SRHGLE (SCLC R{sup 2}=0.9619; NSCLC R{sup 2}=0.787). Conclusion: NSCLC and SCLC exhibited different early RT-responses as assessed by CT-texture changes during RT-delivery. The observed larger changes in various CT-texture parameters for SCLC indicate that SCLC may respond to RT more rapid than NSCLC.

60 APPLIED LIFE SCIENCES↗