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 145 records · Page 8

Time filtering of event based neutron scattering data: A pathway to study the dynamic structural responses of materials

Time-resolved diffraction has become a vital tool for probing dynamic responses to an applied stimulus. Such experiments traditionally use hardware solutions to histogram measured data into their respective bin. We will show that a major advantage of event-based data acquisition, which time-stamps measured diffraction data with 100 ns accuracy, is much preferred over hardware histogramming of the data by enabling postprocessing for advanced custom binning using a software solution. This approach is made even more powerful by coupling measured diffraction data with metadata about the applied stimuli and material response. In this work, we present a time-filter approach that leverages the power of event-based diffraction collection to reduce stroboscopic data measured over many hours into equally weighted segments that represent subsets of the response to a single cycle of the applied stimulus. We demonstrate this approach by observing ferroelectric/ferroelastic domain wall motion during electric field cycling of BaTiO 3 . Furthermore, the developed approach can readily be expanded to investigate other dynamic phenomena using complex sample environments.

47 OTHER INSTRUMENTATION↗

Fluctuation statistics in the scrape-off layer of Alcator C-Mod

We study long time series of the ion saturation current and floating potential, sampled by Langmuir probes dwelled in the outboard mid-plane scrape-off layer and embedded in the lower divertor baffle of Alcator C-Mod. A series of ohmically heated L-mode plasma discharges is investigated with line-averaged plasma density ranging from ${{\bar{n}}_{\text{e}}}/{{n}_{\text{G}}}=0.15$ to 0.42, where $n_ G$ is the Greenwald density. All ion saturation current time series that are sampled in the far scrape-off layer are characterized by large-amplitude burst events. Coefficients of skewness and excess kurtosis of the time series obey a quadratic relationship and their histograms coincide partially upon proper normalization. Histograms of the ion saturation current time series are found to agree well with a prediction of a stochastic model for the particle density fluctuations in scrape-off layer plasmas. The distribution of the waiting times between successive large-amplitude burst events and of the burst amplitudes are approximately described by exponential distributions. The average waiting time and burst amplitude are found to vary weakly with the line-averaged plasma density. Conditional averaging reveals that the radial blob velocity, estimated from floating potential measurements, increases with the normalized burst amplitude in the outboard mid-plane scrape-off layer. For low density discharges, the conditionally averaged waveform of the floating potential associated with large amplitude bursts at the divertor probes has a dipolar shape. In detached divertor conditions the average waveform is random, indicating electrical disconnection of blobs from the sheaths at the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An implementation of neural simulation-based inference for parameter estimation in ATLAS

Neural simulation-based inference (NSBI) is a powerful class of machine-learning-based methods for statistical inference that naturally handles high-dimensional parameter estimation without the need to bin data into low-dimensional summary histograms. Such methods are promising for a range of measurements, including at the Large Hadron Collider, where no single observable may be optimal to scan over the entire theoretical phase space under consideration, or where binning data into histograms could result in a loss of sensitivity. This work develops a NSBI framework for statistical inference, using neural networks to estimate probability density ratios, which enables the application to a full-scale analysis. It incorporates a large number of systematic uncertainties, quantifies the uncertainty due to the finite number of events in training samples, develops a method to construct confidence intervals, and demonstrates a series of intermediate diagnostic checks that can be performed to validate the robustness of the method. As an example, the power and feasibility of the method are assessed on simulated data for a simplified version of an off-shell Higgs boson couplings measurement in the four-lepton final states. This approach represents an extension to the standard statistical methodology used by the experiments at the Large Hadron Collider, and can benefit many physics analyses.

frequentist statistics↗

Discriminative versus generative approaches to simulation-based inference

Most of the fundamental, emergent, and phenomenological parameters of particle and nuclear physics are determined through parametric template fits. Simulations are used to populate histograms which are then matched to data. This approach is inherently lossy, since histograms are binned and low-dimensional. Deep learning has enabled unbinned and high-dimensional parameter estimation through neural likelihood(-ratio) estimation. We compare two approaches for neural simulation-based inference (NSBI): one based on discriminative learning (classification) and one based on generative modeling. These two approaches are directly evaluated on the same datasets, with a similar level of hyperparameter optimization in both cases. In addition to a Gaussian dataset, we study NSBI using a Higgs boson dataset from the FAIR Universe Challenge. We find that both the direct likelihood and likelihood ratio estimation are able to effectively extract parameters with reasonable uncertainties. For the numerical examples and within the set of hyperparameters studied, we found that the likelihood ratio method is more accurate and/or precise. Both methods have a significant spread from the network training and would require ensembling or other mitigation strategies in practice.

high energy physics↗

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗

Geometrical-Based Generative Adversarial Network to Enhance Digital Rock Image Quality

X-ray microcomputed tomography (micro-CT) is a common tool for the study of porous media structures and properties. High-quality micro-CT data are required to accurately capture pore structures. Acquiring high-quality micro-CT data, however, is not always possible, owing to application limitations and experimental constraints. Therefore, we propose a geometrical-based generative adversarial network (GAN) to rapidly restore noisy micro-CT images to their clean counterparts. The training data and related ground-truth (GT) data are scanned for 7 min and 9.5 h, respectively. To evaluate the performance of the geometrical-based GAN, a 6003 voxel image that has never been used for training is reconstructed and compared with the corresponding GT image. Histogram matching and linear normalization are implemented to adjust the histogram of the reconstructed image to that of the GT image. A watershed-based segmentation method is then applied to delineate pore and solid phases. Lastly, we measure the Minkowski functionals and petrophysical properties, including absolute permeability, pore size distribution, drainage capillary pressure-saturation curve, and imbibition relative permeability, to estimate the physical accuracy of the denoised image. The results show that the proposed geometrical-based GAN can accurately restore noisy micro-CT data. By reducing the scanning time from 9.5 h to 7 min, the expenditure of collecting micro-CT can be decreased significantly. This is particularly important for applications where time-lapse images of a dynamic process are required, high-throughput imaging is necessary for real-time data analysis, or where the quantification of large sample volumes is required.

58 GEOSCIENCES↗

MILK : a Python scripting interface to MAUD for automation of Rietveld analysis

Modern diffraction experiments ( e.g. in situ parametric studies) present scientists with many diffraction patterns to analyze. Interactive analyses via graphical user interfaces tend to slow down obtaining quantitative results such as lattice parameters and phase fractions. Furthermore, Rietveld refinement strategies ( i.e. the parameter turn-on-off sequences) tend to be instrument specific or even specific to a given dataset, such that selection of strategies can become a bottleneck for efficient data analysis. Managing multi-histogram datasets such as from multi-bank neutron diffractometers or caked 2D synchrotron data presents additional challenges due to the large number of histogram-specific parameters. To overcome these challenges in the Rietveld software Material Analysis Using Diffraction ( MAUD ), the MAUD Interface Language Kit ( MILK ) is developed along with an updated text batch interface for MAUD . The open-source software MILK is computer-platform independent and is packaged as a Python library that interfaces with MAUD . Using MILK , model selection ( e.g. various texture or peak-broadening models), Rietveld parameter manipulation and distributed parallel batch computing can be performed through a high-level Python interface. A high-level interface enables analysis workflows to be easily programmed, shared and applied to large datasets, and external tools to be integrated with MAUD . Through modification to the MAUD batch interface, plot and data exports have been improved. The resulting hierarchical folders from Rietveld refinements with MILK are compatible with Cinema: Debye–Scherrer , a tool for visualizing and inspecting the results of multi-parameter analyses of large quantities of diffraction data. In this manuscript, the combined Python scripting and visualization capability of MILK is demonstrated with a quantitative texture and phase analysis of data collected at the HIPPO neutron diffractometer.

97 MATHEMATICS AND COMPUTING↗

Bimodal Visualization of Industrial X-Ray and Neutron Computed Tomography Data

Advanced manufacturing creates increasingly complex objects with material compositions that are often difficult to characterize by a single modality. Our collaborating domain scientists are going beyond traditional methods by employing both X-ray and neutron computed tomography to obtain complementary representations expected to better resolve material boundaries. However, the use of two modalities creates its own challenges for visualization, requiring either complex adjustments of bimodal transfer functions or the need for multiple views. Together with experts in nondestructive evaluation, we designed a novel interactive bimodal visualization approach to create a combined view of the co-registered X-ray and neutron acquisitions of industrial objects. Using an automatic topological segmentation of the bivariate histogram of X-ray and neutron values as a starting point, the system provides a simple yet effective interface to easily create, explore, and adjust a bimodal visualization. Here, we propose a widget with simple brushing interactions that enables the user to quickly correct the segmented histogram results. Our semiautomated system enables domain experts to intuitively explore large bimodal datasets without the need for either advanced segmentation algorithms or knowledge of visualization techniques. We demonstrate our approach using synthetic examples, industrial phantom objects created to stress bimodal scanning techniques, and real-world objects, and we discuss expert feedback.

image segmentation↗

An Entropy-Based Test and Development Framework for Uncertainty Modeling in Level-Set Visualizations

We present a simple comparative framework for testing and developing uncertainty modeling in uncertain marching cubes implementations. The selection of a model to represent the probability distribution of uncertain values directly influences the memory use, run time, and accuracy of an uncertainty visualization algorithm. We use an entropy calculation directly on ensemble data to establish an expected result and then compare the entropy from various probability models, including uniform, Gaussian, histogram, and quantile models. Our results verify that models matching the distribution of the ensemble indeed match the entropy. We further show that fewer bins in nonparametric histogram models are more effective whereas large numbers of bins in quantile models approach data accuracy.

Sisneros, Robert↗

Dynamic PET simulator via tomographic emission projection for kinetic modeling and parametric image studies

Purpose: To develop and evaluate a fast and simple tool called dPETSTEP (Dynamic PET Simulator of Tracers via Emission Projection), for dynamic PET simulations as an alternative to Monte Carlo (MC), useful for educational purposes and evaluation of the effects of the clinical environment, postprocessing choices, etc., on dynamic and parametric images. Methods: The tool was developed in MATLAB using both new and previously reported modules of PETSTEP (PET Simulator of Tracers via Emission Projection). Time activity curves are generated for each voxel of the input parametric image, whereby effects of imaging system blurring, counting noise, scatters, randoms, and attenuation are simulated for each frame. Each frame is then reconstructed into images according to the user specified method, settings, and corrections. Reconstructed images were compared to MC data, and simple Gaussian noised time activity curves (GAUSS). Results: dPETSTEP was 8000 times faster than MC. Dynamic images from dPETSTEP had a root mean square error that was within 4% on average of that of MC images, whereas the GAUSS images were within 11%. The average bias in dPETSTEP and MC images was the same, while GAUSS differed by 3% points. Noise profiles in dPETSTEP images conformed well to MC images, confirmed visually by scatter plot histograms, and statistically by tumor region of interest histogram comparisons that showed no significant differences (p < 0.01). Compared to GAUSS, dPETSTEP images and noise properties agreed better with MC. Conclusions: The authors have developed a fast and easy one-stop solution for simulations of dynamic PET and parametric images, and demonstrated that it generates both images and subsequent parametric images with very similar noise properties to those of MC images, in a fraction of the time. They believe dPETSTEP to be very useful for generating fast, simple, and realistic results, however since it uses simple scatter and random models it may not be suitable for studies investigating these phenomena. dPETSTEP can be downloaded free of cost from https://github.com/CRossSchmidtlein/dPETSTEP.

60 APPLIED LIFE SCIENCES↗

SU-C-BRA-06: Automatic Brain Tumor Segmentation for Stereotactic Radiosurgery Applications

Purpose: Stereotactic radiosurgery (SRS), which delivers a potent dose of highly conformal radiation to the target in a single fraction, requires accurate tumor delineation for treatment planning. We present an automatic segmentation strategy, that synergizes intensity histogram thresholding, super-voxel clustering, and level-set based contour evolving methods to efficiently and accurately delineate SRS brain tumors on contrast-enhance T1-weighted (T1c) Magnetic Resonance Images (MRI). Methods: The developed auto-segmentation strategy consists of three major steps. Firstly, tumor sites are localized through 2D slice intensity histogram scanning. Then, super voxels are obtained through clustering the corresponding voxels in 3D with reference to the similarity metrics composited from spatial distance and intensity difference. The combination of the above two could generate the initial contour surface. Finally, a localized region active contour model is utilized to evolve the surface to achieve the accurate delineation of the tumors. The developed method was evaluated on numerical phantom data, synthetic BRATS (Multimodal Brain Tumor Image Segmentation challenge) data, and clinical patients’ data. The auto-segmentation results were quantitatively evaluated by comparing to ground truths with both volume and surface similarity metrics. Results: DICE coefficient (DC) was performed as a quantitative metric to evaluate the auto-segmentation in the numerical phantom with 8 tumors. DCs are 0.999±0.001 without noise, 0.969±0.065 with Rician noise and 0.976±0.038 with Gaussian noise. DC, NMI (Normalized Mutual Information), SSIM (Structural Similarity) and Hausdorff distance (HD) were calculated as the metrics for the BRATS and patients’ data. Assessment of BRATS data across 25 tumor segmentation yield DC 0.886±0.078, NMI 0.817±0.108, SSIM 0.997±0.002, and HD 6.483±4.079mm. Evaluation on 8 patients with total 14 tumor sites yield DC 0.872±0.070, NMI 0.824±0.078, SSIM 0.999±0.001, and HD 5.926±6.141mm. Conclusion: The developed automatic segmentation strategy, which yields accurate brain tumor delineation in evaluation cases, is promising for its application in SRS treatment planning.

60 APPLIED LIFE SCIENCES↗

SU-D-BRC-05: Effects of Motion and Variable RBE in a Lung Patient Treated with Passively Scattered Protons

Purpose: To evaluate effects of motion and variable relative biological effectiveness (RBE) in a lung cancer patient treated with passively scattered proton therapy using dose volume histograms associated with patient dose computed using three different methods. Methods: A proton treatment plan of a lung cancer patient optimized using clinical treatment planning system (TPS) was used to construct a detailed Monte Carlo (MC) model of the beam delivery system and the patient specific aperture and compensator. A phase space file containing all particles transported through the beam line was collected at the distal surface of the range compensator and subsequently transported through two different patient models. The first model was based on the average CT used by the TPS and the second model included all 10 phases of the corresponding 4DCT. The physical dose and proton linear energy transfer (LET) were computed in each voxel of two models and used to compute constant and variable RBE MC dose on average CT and 4D CT. The MC computed doses were compared to the TPS dose using dose volume histograms for relevant structures. Results: The results show significant differences in doses to the target and critical structures suggesting the need for more accurate proton dose computation methods. In particular, the 4D dose shows reduced coverage of the target and higher dose to the spinal cord, while variable RBE dose shows higher lung dose. Conclusion: The methodology developed in this pilot study is currently used for the analysis of a cohort of ∼90 lung patients from a clinical trial comparing proton and photon therapy for lung cancer. The results from this study will help us in determining the clinical significance of more accurate dose computation models in proton therapy.

60 APPLIED LIFE SCIENCES↗

SU-D-207B-04: Morphological Features of MRI as a Correlate of Capsular Contracture in Breast Cancer Patients with Implant-Based Reconstructions

Purpose: Capsular contracture (CC) is a serious complication in patients receiving implant-based reconstruction for breast cancer. The goal of this study was to identify image-based correlates of CC using MRI imaging in breast cancer patients who received both MRI and clinical evaluation following reconstructive surgery. Methods: We analyzed a retrospective dataset of 50 patients who had both a diagnostic MR and a plastic surgeon’s evaluations of CC score (Baker’s score) within a six month period following mastectomy and reconstructive surgery. T2w sagittal MRIs (TR/TE = 3500/102 ms, slice thickness = 4 mm) were used for morphological shape features (roundness, eccentricity, solidity, extent and ratio-length) and histogram features (median, skewness and kurtosis) of the implant and the pectoralis muscle overlying the implant. Implant and pectoralis muscles were segmented in 3D using Computation Environment for Radiological Research (CERR) and shape and histogram features were calculated as a function of Baker’s score. Results: Shape features such as roundness and eccentricity were statistically significant in differentiating grade 1 and grade 2 (p = 0.009; p = 0.06) as well as grade 1 and grade 3 CC (p = 0.001; p = 0.006). Solidity and extent were statistically significant in differentiating grade 1 and grade 3 CC (p = 0.04; p = 0.04). Ratio-length was statistically significant in differentiating all grades of CC except grade 2 and grade 3 that showed borderline significance (p = 0.06). The muscle thickness, median intensity and kurtosis were significant in differentiating between grade 1 and grade 3 (p = 0.02), grade 1 and grade 2 (p = 0.03) and grade 1 and grade 3 (p = 0.01) respectively. Conclusion: Morphological shape features described on MR images were associated with the severity of CC. MRI may be important in objectively evaluating outcomes in breast cancer patients who undergo implant reconstruction.

60 APPLIED LIFE SCIENCES↗

SU-F-R-31: Identification of Robust Normal Lung CT Texture Features for the Prediction of Radiation-Induced Lung Disease

Purpose: Normal lung CT texture features have been used for the prediction of radiation-induced lung disease (radiation pneumonitis and radiation fibrosis). For these features to be clinically useful, they need to be relatively invariant (robust) to tumor size and not correlated with normal lung volume. Methods: The free-breathing CTs of 14 lung SBRT patients were studied. Different sizes of GTVs were simulated with spheres placed at the upper lobe and lower lobe respectively in the normal lung (contralateral to tumor). 27 texture features (9 from intensity histogram, 8 from grey-level co-occurrence matrix [GLCM] and 10 from grey-level run-length matrix [GLRM]) were extracted from [normal lung-GTV]. To measure the variability of a feature F, the relative difference D=|Fref -Fsim|/Fref*100% was calculated, where Fref was for the entire normal lung and Fsim was for [normal lung-GTV]. A feature was considered as robust if the largest non-outlier (Q3+1.5*IQR) D was less than 5%, and considered as not correlated with normal lung volume when their Pearson correlation was lower than 0.50. Results: Only 11 features were robust. All first-order intensity-histogram features (mean, max, etc.) were robust, while most higher-order features (skewness, kurtosis, etc.) were unrobust. Only two of the GLCM and four of the GLRM features were robust. Larger GTV resulted greater feature variation, this was particularly true for unrobust features. All robust features were not correlated with normal lung volume while three unrobust features showed high correlation. Excessive variations were observed in two low grey-level run features and were later identified to be from one patient with local lung diseases (atelectasis) in the normal lung. There was no dependence on GTV location. Conclusion: We identified 11 robust normal lung CT texture features that can be further examined for the prediction of radiation-induced lung disease. Interestingly, low grey-level run features identified normal lung diseases. This work was supported in part by the National Cancer Institute Grants R01CA172638.

60 APPLIED LIFE SCIENCES↗

SU-F-R-56: Early Assessment of Treatment Response During Radiation Therapy Delivery for Esophageal Cancer Using Quantitative CT

Purpose: To investigate the feasibility of assessing treatment response using CTs during delivery of radiation therapy (RT) for esophageal cancer. Methods: Daily CTs acquired using a CT-on-Rails during the routine CT-guided RT for 20 patients with stage II to IV esophageal cancers were analyzed. All patients were treated with combined chemotherapy and IMRT of 45–50 Gy in 25 fractions, and were followed up for two years. Contours of GTV, spinal cord, and non-specified tissue (NST) irradiated with low dose were generated on each daily CT. A series of CT-texture metrics including Hounsfield Unit (HU) histogram, mean HU, standard derivation (STD), entropy, and energy were obtained in these contours on each daily CT. The changes of these metrics and GTV volume during RT delivery were calculated and correlated with treatment outcome. Results: Changes in CT texture (e.g., HU histogram) in GTV and spinal cord (but not in NST) were observed during RT delivery and were consistently increased with radiation dose. For the 20 cases studied, the mean HU in GTV was reduced on average by 4.0HU from the first to the last fractions, while 8 patients (responders) had larger reductions in GTV mean HU (average 7.8 HU) with an average GTV reduction of 51% and had increased consistently in GTV STD and entropy with radiation dose. The rest of 12 patients (non-responders) had lower reductions in GTV mean HU (average 1.5HU) and almost no change in STD and entropy. For the 8 responders, 2 experienced complete response, 7 (88%) survived and 1 died. In contrast, for the 12 non-responders, 4 (33%) survived and 8 died. Conclusion: Radiation can induce changes in CT texture in tumor (e.g., mean HU) during the delivery of RT for esophageal cancer. If validated with more data, such changes may be used for early prediction of RT response for esophageal cancer.

60 APPLIED LIFE SCIENCES↗

SU-F-J-207: Non-Small Cell Lung Cancer Patient Survival Prediction with Quantitative Tumor Textures Analysis in Baseline CT

Purpose: Chemo-radiation therapy (CRT) is widely used in treating patients with locally advanced non-small cell lung cancer (NSCLC). Determination of the likelihood of patient response to treatment and optimization of treatment regime is of clinical significance. Up to date, no imaging biomarker has reliably correlated to NSCLC patient survival rate. This pilot study is to extract CT texture information from tumor regions for patient survival prediction. Methods: Thirteen patients with stage II-III NSCLC were treated using CRT with a median dose of 6210 cGy. Non-contrast-enhanced CT images were acquired for treatment planning and retrospectively collected for this study. Texture analysis was applied in segmented tumor regions using the Local Binary Pattern method (LBP). By comparing its HU with neighboring voxels, the LBPs of a voxel were measured in multiple scales with different group radiuses and numbers of neighbors. The LBP histograms formed a multi-dimensional texture vector for each patient, which was then used to establish and test a Support Vector Machine (SVM) model to predict patients’ one year survival. The leave-one-out cross validation strategy was used recursively to enlarge the training set and derive a reliable predictor. The predictions were compared with the true clinical outcomes. Results: A 10-dimensional LBP histogram was extracted from 3D segmented tumor region for each of the 13 patients. Using the SVM model with the leave-one-out strategy, only 1 out of 13 patients was misclassified. The experiments showed an accuracy of 93%, sensitivity of 100%, and specificity of 86%. Conclusion: Within the framework of a Support Vector Machine based model, the Local Binary Pattern method is able to extract a quantitative imaging biomarker in the prediction of NSCLC patient survival. More patients are to be included in the study.

60 APPLIED LIFE SCIENCES↗

SU-F-T-662: Feasibility Study of Fe3O4/TaOx Nano Particles as a Radiosensitizer for Radiation Therapy

Purpose: To investigate the feasibility of using multifunctional Fe{sub 3}O{sub 4}/TaOx(core / shell) nano particles developed for CT and MRI contrast agent as dose enhancing radiosensitizers. Methods: Firstly, to verify the imaging detectability of Fe{sub 3}O{sub 4}/TaOx nano particles, in-vivo tests were conducted. Approximately 600 mg/kg of 19 nm diameter Fe{sub 3}O{sub 4}/TaOx nano particles dispersed in phosphate buffered saline(PBS) were injected to ten nude Balb/c mice through the tail vein. Difference between pre- and post-injection images was analyzed by computing the pixel histogram and correlation coefficient factor using MATLAB in the user defined ROI. Secondly, to quantify the potential therapeutic enhancement with nano materials, DER (Dose Enhancement Ratio) and number of SER (Secondary Electron Ratio) were computed using TOPAS(ver.2.0 P-03) MC simulation. Results: In CT, MRI imaging, the aorta, the blood vessel, and the liver were clearly visualized after intravenous injection of Fe{sub 3}O{sub 4}/TaOx nano particles. There was large different between pre and post-injection images of Histogram data and Coefficients of correlation factor in CT and MR are 0.006, 0.060, respectively. When 70 MeV protons were irradiated for a Gold, Tantalum, TaOx, Fe{sub 3}O{sub 4}/TaOx, Fe{sub 3}O{sub 4} nano particle, DER was 2.59, 2.41, 1.68, 1.54 and 1.36 respectively. Similarly, SER increment was 2.31, 2.15, 1.56, 1.46, and 1.27 for Gold, Tantalum, TaOx, Fe{sub 3}O{sub 4}/TaOx, Fe{sub 3}O{sub 4} nano particle, respectively. Conclusion: Fe{sub 3}O{sub 4}/TaOx nano particles have potential as a multifunctional agent which enhances the accuracy in cancer detection through visualization of developed tumor lesion and increases the therapeutic effect in proton therapy. The dose enhancement with Fe{sub 3}O{sub 4}/TaOx was estimated as half of the Gold. However, tumor targeting such as combined with magnetic field may overcome the low DER. This research was supported by the NRF funded by the Ministry of Science, ICT & Future Planning (2012M3A9B6055201 and 2012R1A1A2042414), Samsung Medical Center grant[GFO1130081].

60 APPLIED LIFE SCIENCES↗

SU-G-TeP1-09: Modality-Specific Dose Gradient Modeling for Prostate IMRT Using Spherical Distance Maps of PTV and Isodose Contours

Purpose: Overlapping volume histogram (OVH) and distance-to-target histogram (DTH) calculations rely on the assumption that dose gradients are symmetric with respect to primary target volume (PTV) expansion and minimum distance to PTV surface, respectively. It is desirable to lift this assumption and instead account for achievable modality-specific dose gradients (MSDG) for a given PTV shape. Methods: From a library of 96 prostate 7-beam IMRT plans, we computed spherical distance maps (SDMs) for PTVs and 3 iso-dose contours. Each SDM contains the minimum distances between plan isocenter and the object’s surface for a fixed set of azimuthal and polar angles. We performed principal component analysis (PCA)-based missing data recovery with PTV SDM as input and a single iso-dose contour SDM as output. Repeating this process for the set of iso-dose contours sparsely reconstructed the MSDG for a given PTV. DVH points were computed from the MSDG for bladder and rectum (OARs) as a natural way of casting patient-specific geometric information into modality-specific dose space (vs. OVH and DTH where data is purely geometric). For comparison, we implemented a cumulative (c)DTH-based prediction algorithm, and produced DVH for both OARs separately. We then computed linear regressions between each method and the DVH of the original patient plan for three different OAR DVH points. Results: R-squared for MSD-Gcomputed DVH vs. original plan DVH at V90%, V80% and V60% of max dose were 0.86, 0.83, and 0.78, for bladder and 0.36, 0.52, and 0.55 for rectum, respectively; R-squared for cDTH-based predicted DVH vs. original plan DVH were 0.67, 0.81, and 0.83, for bladder and 0.44, 0.50, and 0.51, for rectum, respectively. Conclusion: By simply modeling the MSDG about a given PTV, we are able to reproduce DVHs consistent with a validated cDTH-based DVH prediction. Regression analysis suggests MSDG could further improve predictions.

60 APPLIED LIFE SCIENCES↗