Search NASA⌕ Search

SEARCH · Search NASA

Results for “LINEAR FILTER”

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 19 records

The Use of Photon Beams of a Flattening Filter-free Linear Accelerator for Hypofractionated Volumetric Modulated Arc Therapy in Localized Prostate Cancer

Purpose: To evaluate the potential usage of flattening filter-free (FFF) photon beams in the treatment of prostate cancer. Methods and Materials: Volumetric-modulated arc therapy (VMAT) treatment planning was performed for 7 patients using TrueBeam{sup Registered-Sign} linear accelerator and photon beams with (X6, X10) and without (X6FFF, X10FFF) flattening filter. Prescribed dose was 19 Multiplication-Sign 3 Gy = 57 Gy. One or two 360 Degree-Sign arcs with dose rate of 600 MU/min for flattened beams, and 1,200 MU/min for FFF beams were used. Results: No difference was detected between the four beams in PTV coverage, conformity, and homogeneity. Mean body dose and body volume receiving 50% of the prescribed dose decreased with increasing mean energy (r{sup 2} = 0.8275, p < 0.01). X6FFF delivered 3.6% more dose compared with the X6 (p < 0.01). X10FFF delivered 3.0% (p < 0.01), and the X10 5.8% (p < 0.01) less mean body dose compared with X6. There was a significant increase in the mean dose to the rectum for the X10 compared with X6 (2.6%, p < 0.01). Mean dose to the bladder increased by 1.3% for X6FFF and decreased by 2.3% for X10FFF. Using a single arc and FFF, treatment time was reduced by 35% (2 SD = 10%). Conclusion: FFF beams resulted in dose distributions similar to flattened beams. X10FFF beam provided the best solution, sparing rectum and bladder and minimizing whole-body dose. FFF beams lead to a time efficient treatment delivery, particularly when combined with hypofractionated VMAT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Adaptive linear predictor FIR filter based on the Cyclone V FPGA with HPS to reduce narrow band RFI in AERA radio detection of cosmic rays

We present the new approach to a filtering of radio frequency interferences (RFI) in the Auger Engineering Radio Array (AERA) which study the electromagnetic part of the Extensive Air Showers. The radio stations can observe radio signals caused by coherent emissions due to geomagnetic radiation and charge excess processes. AERA observes frequency band from 30 to 80 MHz. This range is highly contaminated by human-made RFI. In order to improve the signal to noise ratio RFI filters are used in AERA to suppress this contamination. The first kind of filter used by AERA was the Median one, based on the Fast Fourier Transform (FFT) technique. The second one, which is currently in use, is the infinite impulse response (IIR) notch filter. The proposed new filter is a finite impulse response (FIR) filter based on a linear prediction (LP). A periodic contamination hidden in a registered signal (digitized in the ADC) can be extracted and next subtracted to make signal cleaner. The FIR filter requires a calculation of n=32, 64 or even 128 coefficients (dependent on a required speed or accuracy) by solving of n linear equations with coefficients built from the covariance Toeplitz matrix. This matrix can be solved by the Levinson recursion, which is much faster than the Gauss procedure. The filter has been already tested in the real AERA radio stations on Argentinean pampas with a very successful results. The linear equations were solved either in the virtual soft-core NIOSR processor (implemented in the FPGA chip as a net of logic elements) or in the external Voipac PXA270M ARM processor. The NIOS processor is relatively slow (50 MHz internal clock), calculations performed in an external processor consume a significant amount of time for data exchange between the FPGA and the processor. Test showed a very good efficiency of the RFI suppression for stationary (long-term) contaminations. However, we observed a short-time contaminations, which could not be suppressed either by the IIR-notch filter or by the FIR filter based on the linear predictions. For the LP FIR filter the refreshment time of the filter coefficients was to long and filter did not keep up with the changes of a contamination structure, mainly due to a long calculation time in a slow processors. We propose to use the Cyclone V SE chip with embedded micro-controller operating with 925 MHz internal clock to significantly reduce a refreshment time of the FIR coefficients. The lab results are promising. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

F3D Image Processing and Analysis for Many - and Multi-core Platforms

F3D is written in OpenCL, so it achieve[sic] platform-portable parallelism on modern mutli-core CPUs and many-core GPUs. The interface and mechanims to access F3D core are written in Java as a plugin for Fiji/ImageJ to deliver several key image-processing algorithms necessary to remove artifacts from micro-tomography data. The algorithms consist of data parallel aware filters that can efficiently utilizes[sic] resources and can work on out of core datasets and scale efficiently across multiple accelerators. Optimizing for data parallel filters, streaming out of core datasets, and efficient resource and memory and data managements over complex execution sequence of filters greatly expedites any scientific workflow with image processing requirements. F3D performs several different types of 3D image processing operations, such as non-linear filtering using bilateral filtering and/or median filtering and/or morphological operators (MM). F3D gray-level MM operators are one-pass constant time methods that can perform morphological transformations with a line-structuring element oriented in discrete directions. Additionally, MM operators can be applied to gray-scale images, and consist of two parts: (a) a reference shape or structuring element, which is translated over the image, and (b) a mechanism, or operation, that defines the comparisons to be performed between the image and the structuring element. This tool provides a critical component within many complex pipelines such as those for performing automated segmentation of image stacks. F3D is also called a "descendent" of Quant-CT, another software we developed in the past. These two modules are to be integrated in a next version. Further details were reported in: D.M. Ushizima, T. Perciano, H. Krishnan, B. Loring, H. Bale, D. Parkinson, and J. Sethian. Structure recognition from high-resolution images of ceramic composites. IEEE International Conference on Big Data, October 2014.

Source record↗

F3D Image Processing and Analysis for Many - and Multi-core Platforms

F3D is written in OpenCL, so it achieve[sic] platform-portable parallelism on modern mutli-core CPUs and many-core GPUs. The interface and mechanims to access F3D core are written in Java as a plugin for Fiji/ImageJ to deliver several key image-processing algorithms necessary to remove artifacts from micro-tomography data. The algorithms consist of data parallel aware filters that can efficiently utilizes[sic] resources and can work on out of core datasets and scale efficiently across multiple accelerators. Optimizing for data parallel filters, streaming out of core datasets, and efficient resource and memory and data managements over complex execution sequence of filters greatly expedites any scientific workflow with image processing requirements. F3D performs several different types of 3D image processing operations, such as non-linear filtering using bilateral filtering and/or median filtering and/or morphological operators (MM). F3D gray-level MM operators are one-pass constant time methods that can perform morphological transformations with a line-structuring element oriented in discrete directions. Additionally, MM operators can be applied to gray-scale images, and consist of two parts: (a) a reference shape or structuring element, which is translated over the image, and (b) a mechanism, or operation, that defines the comparisons to be performed between the image and the structuring element. This tool provides a critical component within many complex pipelines such as those for performing automated segmentation of image stacks. F3D is also called a "descendent" of Quant-CT, another software we developed in the past. These two modules are to be integrated in a next version. Further details were reported in: D.M. Ushizima, T. Perciano, H. Krishnan, B. Loring, H. Bale, D. Parkinson, and J. Sethian. Structure recognition from high-resolution images of ceramic composites. IEEE International Conference on Big Data, October 2014.

Source record↗

Metamaterials-based sensor to detect and locate nonlinear elastic sources

In recent years, acoustic metamaterials have attracted increasing scientific interest for very diverse technological applications ranging from sound abatement to ultrasonic imaging, mainly due to their ability to act as band-stop filters. At the same time, the concept of chaotic cavities has been recently proposed as an efficient tool to enhance the quality of nonlinear signal analysis, particularly in the ultrasonic/acoustic case. The goal of the present paper is to merge the two concepts in order to propose a metamaterial-based device that can be used as a natural and selective linear filter for the detection of signals resulting from the propagation of elastic waves in nonlinear materials, e.g., in the presence of damage, and as a detector for the damage itself in time reversal experiments. Numerical simulations demonstrate the feasibility of the approach and the potential of the device in providing improved signal-to-noise ratios and enhanced focusing on the defect locations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prediction of Power Measurements Using Adaptive Filters

With the advent of smart grid concept, Internet of Things (IoT) and the deployment of smart meters, the cyberattack threats on power networks have increased due to the use of communication systems that can be accessed by adversaries. Attackers will have the ability to manipulate the outcomes of smart meters which in turn influence the core application of Energy Management System 9EMS): State Estimation (SE). Bad data analytic tools may fail to detect some attacks into measurements. Meanwhile, Machine Learning (ML) solutions have been proposed for detecting False Data Injection (FDI) attacks. However, there is a lack of ML time-series solutions presented in the state-of-the-art that is yet to be complex. In signal processing, time-series solutions do not only consider the signal, but also the statistics of the signal over time. Therefore, in this paper, a machine learning for time-series solutions is presented as an application to model the measurements of the power grid that are used in SE. The presented model takes into account adaptive linear and non-linear filters: Finite Impulse Response (FIR), and Infinite Impulse Response (IIR). The presented models are implemented and performed on the IEEE-118 bus system. The results indicate the advantage of applying those filters over the state-of-the-art machine learning solutions.

Hamad, Khaled↗

Minimum feature size control in level set topology optimization via density fields

A level set topology optimization approach that uses an auxiliary density field to nucleate holes during the optimization process and achieves minimum feature size control in optimized designs is explored. The level set field determines the solid-void interface and the density field describes the distribution of a fictitious porous material using the solid isotropic material with penalization. These fields are governed by two sets of independent optimization variables which are initially coupled using a penalty for hole nucleation. The strength of the density field penalization and projection is gradually increased during the optimization process to promote a 0-1 density distribution. In addition, a second penalty regulates the evolution of the density field in the void phase. The treatment of the density field combined with the second penalty mitigate the appearance of small design features. The minimum feature size of optimized designs is controlled by the radius of the linear filter applied to the density optimization variables. The structural response is predicted by the extended finite element method, the sensitivities by the adjoint method, and the optimization variables are updated by a gradient-based optimization algorithm. Numerical examples investigate the robustness of this approach with respect to algorithmic parameters and mesh refinement. The results show the applicability of the combined density level set topology optimization approach for both optimal hole nucleation and for minimum feature size control in 2D and 3D. This comes, however, at the cost of a more complex problem formulation and additional computational cost due to an increased number of optimization variables.

42 ENGINEERING↗

Methodology and performance comparison of statistical learning pulse shape classifiers as demonstrated with organic liquid scintillator [Performance comparison of statistical learning pulse shape classifiers]

In this study, we present novel methods for automated pulse shape discrimination. The classifiers are trained using simple radionuclide sources and do not require ground truth labeling. We test their performance using labels derived from time of flight experiments and present the results in terms of energy-dependent Receiver Operating Characteristic (ROC) curves. In addition, we also train and test standard pulse shape discrimination methods on the same data for comparison. We find multiple methods that can yield similar false neutron and true neutron rates at 24 keVee as tail-to-total or Gatti’s optimal linear filter yield at 54 keVee.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Convolutional Dictionary Learning: A Comparative Review and New Algorithms

Convolutional sparse representations are a form of sparse representation with a dictionary that has a structure that is equivalent to convolution with a set of linear filters. Additionally, while effective algorithms have recently been developed for the convolutional sparse coding problem, the corresponding dictionary learning problem is substantially more challenging. Furthermore, although a number of different approaches have been proposed, the absence of thorough comparisons between them makes it difficult to determine which of them represents the current state of the art. The present work both addresses this deficiency and proposes some new approaches that outperform existing ones in certain contexts. A thorough set of performance comparisons indicates a very wide range of performance differences among the existing and proposed methods, and clearly identifies those that are the most effective.

97 MATHEMATICS AND COMPUTING↗

Stochastic analysis of surface metrology

The design and evaluation of the expected performance of optical systems require sophisticated and reliable information about the surface topography for planned optical elements before they are fabricated. Modern x-ray source facilities are reliant upon the availability of optics with unprecedented quality (surface slope accuracy <0.1 μrad). The problem is especially complex in the case of x-ray optics, particularly for the X-ray Surveyor under development and other missions. The high angular resolution and throughput of future x-ray space observatories requires hundreds of square meters of high-quality optics. The uniqueness of the optics and limited number of proficient vendors makes the fabrication extremely time consuming and expensive, mostly due to the limitations in accuracy and measurement rate of metrology used in fabrication. Here, we discuss improvements in metrology efficacy via comprehensive statistical analysis of a compact volume of metrology data. The data are considered stochastic, and a statistical model called invertible time-invariant linear filter (InTILF) is developed now for two-dimensional (2-D) surface profiles to provide compact description of the 2-D data in addition to one-dimensional data treated so far. The InTILF model captures stochastic patterns in the data and can be used as a quality metric and feedback to polishing processes, avoiding high-resolution metrology measurements over the entire optical surface. The modeling, implemented in our BeatMark™ software, allows simulating metrology data for optics made by the same vendor and technology. The data are vital for reliable specification for optical fabrication, to be exactly adequate for the required system performance.

42 ENGINEERING↗

SU-C-207B-05: Tissue Segmentation of Computed Tomography Images Using a Random Forest Algorithm: A Feasibility Study

Purpose: Develop an automated Random Forest algorithm for tissue segmentation of CT examinations. Methods: Seven materials were classified for segmentation: background, lung/internal gas, fat, muscle, solid organ parenchyma, blood/contrast, and bone using Matlab and the Trainable Weka Segmentation (TWS) plugin of FIJI. The following classifier feature filters of TWS were investigated: minimum, maximum, mean, and variance each evaluated over a pixel radius of 2n, (n = 0–4). Also noise reduction and edge preserving filters, Gaussian, bilateral, Kuwahara, and anisotropic diffusion, were evaluated. The algorithm used 200 trees with 2 features per node. A training data set was established using an anonymized patient’s (male, 20 yr, 72 kg) chest-abdomen-pelvis CT examination. To establish segmentation ground truth, the training data were manually segmented using Eclipse planning software, and an intra-observer reproducibility test was conducted. Six additional patient data sets were segmented based on classifier data generated from the training data. Accuracy of segmentation was determined by calculating the Dice similarity coefficient (DSC) between manual and auto segmented images. Results: The optimized autosegmentation algorithm resulted in 16 features calculated using maximum, mean, variance, and Gaussian blur filters with kernel radii of 1, 2, and 4 pixels, in addition to the original CT number, and Kuwahara filter (linear kernel of 19 pixels). Ground truth had a DSC of 0.94 (range: 0.90–0.99) for adult and 0.92 (range: 0.85–0.99) for pediatric data sets across all seven segmentation classes. The automated algorithm produced segmentation with an average DSC of 0.85 ± 0.04 (range: 0.81–1.00) for the adult patients, and 0.86 ± 0.03 (range: 0.80–0.99) for the pediatric patients. Conclusion: The TWS Random Forest auto-segmentation algorithm was optimized for CT environment, and able to segment seven material classes over a range of body habitus and CT protocol parameters with an average DSC of 0.86 ± 0.04 (range: 0.80–0.99).

60 APPLIED LIFE SCIENCES↗

NA-22 Particle Classifier Quarterly Report FY18 Q1 [Automated Particle Classifier]

This is the first quarter of the third and final year of this project. Much of this quarter was taken up with preparing the comprehensive paper that details the collection of appropriate data, the pre-processing of data, the forms and development of the algorithms, and the scoring of performance. These four processes were still a bit unpolished among the different algorithms, so by synergy the paper raised the quality of performance of each algorithm. For completeness, the group implemented two more kinds of algorithms for performance comparison: Gatti-DeMartini’s optimal linear filter, and a more straightforward (and more effective) implementation of principal component analysis (PCA) than we had previously built.

47 OTHER INSTRUMENTATION↗

Autonomous Correction of Sensor Data Applied to Building Technologies Using Filtering Methods

Sensor data validity is extremely important in a number of applications, particularly building technologies where collected data are used to determine performance. An example of this is Oak Ridge National Laboratory s ZEBRAlliance research project, which consists of four single-family homes located in Oak Ridge, TN. The homes are outfitted with a total of 1,218 sensors to determine the performance of a variety of different technologies integrated within each home. Issues arise with such a large amount of sensors, such as missing or corrupt data. This paper aims to eliminate these problems using: (1) Kalman filtering and (2) linear prediction filtering techniques. Five types of data are the focus of this paper: (1) temperature; (2) humidity; (3) energy consumption; (4) pressure; and (5) airflow. Simulations show the Kalman filtering method performed best in predicting temperature, humidity, pressure, and airflow data, while the linear prediction filtering method performed best with energy consumption data.

Castello, Charles C↗

SU-G-IeP2-08: Investigation On Signal Detectability in Volumetric Cone Beam CT Images with Anatomical Background

Purpose: To investigate the slice direction dependent detectability in cone beam CT images with anatomical background. Methods: We generated 3D anatomical background images using breast anatomy model. To generate 3D breast anatomy, we filtered 3D Gaussian noise with a square root of 1/f{sup 3}, and then assigned the attenuation coefficient of glandular (0.8cm{sup −1}) and adipose (0.46 cm{sup −1}) tissues based on voxel values. Projections were acquired by forward projection, and quantum noise was added to the projection data. The projection data were reconstructed by FDK algorithm. We compared the detectability of a 3 mm spherical signal in the image reconstructed from four different backprojection Methods: Hanning weighted ramp filter with linear interpolation (RECON1), Hanning weighted ramp filter with Fourier interpolation (RECON2), ramp filter with linear interpolation (RECON3), and ramp filter with Fourier interpolation (RECON4), respectively. We computed task SNR of the spherical signal in transverse and longitudinal planes using channelized Hotelling observer with Laguerre-Gauss channels. Results: Transverse plane has similar task SNR values for different backprojection methods, while longitudinal plane has a maximum task SNR value in RECON1. For all backprojection methods, longitudinal plane has higher task SNR than transverse plane. Conclusion: In this work, we investigated detectability for different slice direction in cone beam CT images with anatomical background. Longitudinal plane has a higher task SNR than transverse plane, and backprojection with hanning weighted ramp filter with linear interpolation method (i.e., RECON1) produced the highest task SNR among four different backprojection methods. This research was supported by the MSIP (Ministry of Science, ICT and Future Planning), Korea, under the IT Consilience Creative Programs(IITP-2015-R0346-15-1008) supervised by the IITP (Institute for Information & Communications Technology Promotion), Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the MSIP (2015R1C1A1A01052268) and framework of international cooperation program managed by NRF (NRF-2015K2A1A2067635).

60 APPLIED LIFE SCIENCES↗

The feasibility of a piecewise-linear dynamic bowtie filter

Purpose: The prepatient attenuator (or 'bowtie filter') in CT is used to modulate the flux as a function of fan angle of the x-ray beam incident on the patient. Traditional, static bowtie filters are tailored only for very generic scans and for the average patient. The authors propose a design for a dynamic bowtie that can produce a time-dependent piecewise-linear attenuation profile. This dynamic bowtie may reduce dynamic range, dose or scatter, but in this work they focus on its ability to reduce dynamic range, which may be particularly important for systems employing photon-counting detectors. Methods: The dynamic bowtie is composed of a set of triangular wedges. Each wedge is independently moved in order to produce a time-dependent piecewise-linear attenuation profile. Simulations of the bowtie are conducted to estimate the dynamic range reduction in six clinical datasets. The control of the dynamic bowtie is determined by solving a convex optimization problem, and the dose is estimated using Monte Carlo techniques. Beam hardening artifacts are also simulated. Results: The dynamic range is reduced by factors ranging from 2.4 to 27 depending on the part of the body studied. With a dynamic range minimization objective, the dose to the patient can be reduced from 6% to 33% while maintaining peak image noise. Further reduction in dose may be possible with a specific dose reduction objective. Beam hardening artifacts are suppressed with a two-pass algorithm. Conclusions: A dynamic bowtie producing a time-dependent, piecewise-linear attenuation profile is possible and can be used to modulate the flux of the scanner to the imaging task. Initial simulations show a large reduction in dynamic range. Several other applications are possible.

60 APPLIED LIFE SCIENCES↗

SU-E-I-10: Investigation On Detectability of a Small Target for Different Slice Direction of a Volumetric Cone Beam CT Image

Purpose: To investigate the detectability of a small target for different slice direction of a volumetric cone beam CT image and its impact on dose reduction. Methods: Analytic projection data of a sphere object (1 mm diameter, 0.2/cm attenuation coefficient) were generated and reconstructed by FDK algorithm. In this work, we compared the detectability of the small target from four different backprojection Methods: hanning weighted ramp filter with linear interpolation (RECON 1), hanning weighted ramp filter with Fourier interpolation (RECON2), ramp filter with linear interpolation (RECON 3), and ramp filter with Fourier interpolation (RECON4), respectively. For noise simulation, 200 photons per measurement were used, and the noise only data were reconstructed using FDK algorithm. For each reconstructed volume, axial and coronal slice were extracted and detection-SNR was calculated using channelized Hotelling observer (CHO) with dense difference-of-Gaussian (D-DOG) channels. Results: Detection-SNR of coronal images varies for different backprojection methods, while axial images have a similar detection-SNR. Detection-SNR{sup 2} ratios of coronal and axial images in RECON1 and RECON2 are 1.33 and 1.15, implying that the coronal image has a better detectability than axial image. In other words, using coronal slices for the small target detection can reduce the patient dose about 33% and 15% compared to using axial slices in RECON 1 and RECON 2. Conclusion: In this work, we investigated slice direction dependent detectability of a volumetric cone beam CT image. RECON 1 and RECON 2 produced the highest detection-SNR, with better detectability in coronal slices. These results indicate that it is more beneficial to use coronal slice to improve detectability of a small target in a volumetric cone beam CT image. This research was supported by the MSIP (Ministry of Science, ICT and Future Planning), Korea, under the IT Consilience Creative Program (NIPA-2014-H0201-14-1002) supervised by the NIPA (National IT Industry Promotion Agency). Authors declares that s/he has no conflict of Interest in relation to the work in this abstract.

60 APPLIED LIFE SCIENCES↗

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SIMULATIONS OF EARLY BARYONIC STRUCTURE FORMATION WITH STREAM VELOCITY. II. THE GAS FRACTION

Understanding the gas content of high-redshift halos is crucial for studying the formation of the first generation of galaxies and reionization. Recently, Tseliakhovich and Hirata showed that the relative 'stream' velocity between the dark matter and baryons at the time of recombination-formally a second-order effect, but an unusually large one-can influence the later structure formation history of the universe. We quantify the effect of the stream velocity on the so-called characteristic mass-the minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation epoch-using three different high-resolution sets of cosmological simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between the dark matter and baryons. In order to understand this effect theoretically, we generalize the linear theory filtering mass to properly account for the difference between the dark matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical ansatzes for the characteristic mass are substantially less precise.

79 ASTRONOMY AND ASTROPHYSICS↗