Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Filter error”

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 1,225 records · Page 68

Model-based VQ for image data archival, retrieval and distribution

An ideal image compression technique for image data archival, retrieval and distribution would be one with the asymmetrical computational requirements of Vector Quantization (VQ), but without the complications arising from VQ codebooks. Codebook generation and maintenance are stumbling blocks which have limited the use of VQ as a practical image compression algorithm. Model-based VQ (MVQ), a variant of VQ described here, has the computational properties of VQ but does not require explicit codebooks. The codebooks are internally generated using mean removed error and Human Visual System (HVS) models. The error model assumed is the Laplacian distribution with mean, lambda-computed from a sample of the input image. A Laplacian distribution with mean, lambda, is generated with uniform random number generator. These random numbers are grouped into vectors. These vectors are further conditioned to make them perceptually meaningful by filtering the DCT coefficients from each vector. The DCT coefficients are filtered by multiplying by a weight matrix that is found to be optimal for human perception. The inverse DCT is performed to produce the conditioned vectors for the codebook. The only image dependent parameter used in the generation of codebook is the mean, lambda, that is included in the coded file to repeat the codebook generation process for decoding.

Manohar, Mareboyana↗

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗

A Spatially Coherent Global Soil Moisture Product with Improved Temporal Resolution

Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.

spatial↗

Estimating Bare Earth in Sparse Boreal Forests With WorldView Stereo Imagery

Circumboreal forests are currently experiencing rapid climate warming which is altering their structure, productivity, and status as a carbon sink. Very-high resolution (VHR, <2 m) stereo derived digital surface models are available to monitor these forests, but, a similar resolution digital terrain model (DTM) is required to extract information about tree height, which is often used to estimate carbon content. To our knowledge, no openly available VHR DTM currently exists. To address this need, we developed approaches to extract DTMs by filtering VHR stereo point clouds (PCs) in sparse canopies of Alaska. Our evaluation consisted of two stereo processing methods with three PC search radii, at six different tree canopy cover (TCC) intervals. We found VHR DTMs were robust for estimating bare ground at TCC intervals less than 40% with vertical errors <1.6 m using airborne small footprint light detection and ranging (LiDAR) as reference.

Christopher S R Neigh↗

Estimating Bare Earth in Sparse Boreal Forests With WorldView Stereo Imagery

Circumboreal forests are currently experiencing rapid climate warming which is altering their structure, productivity, and status as a carbon sink. Very-high resolution (VHR, <2 m) stereo derived digital surface models are available to monitor these forests, but, a similar resolution digital terrain model (DTM) is required to extract information about tree height, which is often used to estimate carbon content. To our knowledge, no openly available VHR DTM currently exists. To address this need, we developed approaches to extract DTMs by filtering VHR stereo point clouds (PCs) in sparse canopies of Alaska. Our evaluation consisted of two stereo processing methods with three PC search radii, at six different tree canopy cover (TCC) intervals. We found VHR DTMs were robust for estimating bare ground at TCC intervals less than 40% with vertical errors <1.6 m using airborne small footprint light detection and ranging (LiDAR) as reference.

Christopher S. R. Neigh↗

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗

Iterative Learning Control for Video-rate Atomic Force Microscopy

We present a control scheme for video-rate atomic force microscopy with rosette pattern. The controller structure involves a feedback internal-model-based controller and a feedforward iterative learning controller. The iterative learning controller is designed to improve tracking performance of the feedback-controlled scanner by rejecting the repetitive disturbances arising from the system nonlinearities. We investigate the performance of two inversion techniques for constructing the learning filter. We conduct tracking experiments using a two-degree-of-freedom microelectromechanical system (MEMS) nanopositioner at frame rates ranging from 5 to 20 frames per second. Furthermore, the results reveal that the algorithm converges rapidly and the iterative learning controller significantly reduces both the transient and steady-state tracking errors. We acquire and report a series of high-resolution time-lapsed video-rate AFM images with the rosette pattern.

42 ENGINEERING↗

Identification of Possible Stellar Companions via Speckle Interferometry in a Sample of Be Stars

We present speckle interferometric measurements of Be stars made with the 1.60 m Perkin-Elmer telescope at the Pico dos Dias Observatory of the Laboratório Nacional de Astrofísica in Brazil. Using three different filters (V, R, and I), we obtained astrometric measurements for 23 stars, during three observing runs allocated in 2017 March, 2018 April, and 2019 July. We identified pairs with angular separation ranging from 0.″10 to 6.″69, including six systems with separations smaller than 1.″0. We estimated a mean error in separation of 15 mas, 1.°2 in position angle, and 0.29 mag for the magnitude difference. We discovered three new possible double systems that are not identified in the literature: HD 174886, HD 159032, and HD 168135 with an average angular separation of 0.″20, 3.″23, and 5.″41, respectively. In addition, we identified four new pairs in multiple systems previously known in the Washington Double Star Catalog. However, for 16 stars of our sample, we did not find evidence of binarity, based on the diffraction-limited resolution, the field of view, and/or the sensitivity of our observations. From the cross-match with data from Gaia Data Release 2, we verified the possibility of some of these systems be gravitationally bound or not.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation in nonlinear systems with transport delay.

The problem of estimation of state in nonlinear dynamical systems containing time delays is studied. The plant is specified by a set of nonlinear differential-difference equations. Observations are a nonlinear function of current and/or delayed states. Both contain additive disturbances. The criterion used for the optimal estimates is the integral of the weighted squared error. Using the theory of the calculus of variations, equations are developed for the estimation. They are first expressed in the form of a split boundary value problem, which is then converted to an initial value problem for on-line estimation. The result yields a sequential estimation scheme in which filtered and smoothed estimates are computed in a sequential manner. The applicability of the procedure is demonstrated by a practical example.

Stoller, R. L.↗

Attitude determination techniques for Nimbus 6 and Nimbus G

Techniques for measuring attitude deviations in the Nimbus 6 and Nimbus G 3-axis stabilized, zero-momentum, earth observation satellites are presented. The on-orbit telemetered data include scanner roll and pitch error signals, an illuminated DSAS which can read the azimuth and elevation of the sunline when the spacecraft is not in the umbra, the signal from the single gyro in the RMP, and the momenta of the reaction wheels. Mathematical methods for interpreting the measured data are described with attention to filter equations, especially the Kalman filter because of its ability to eliminate high frequency noise.

Das, A.↗

Combined optimization of image-gathering and image-processing systems for scene feature detection

The relationship between the image gathering and image processing systems for minimum mean squared error estimation of scene characteristics is investigated. A stochastic optimization problem is formulated where the objective is to determine a spatial characteristic of the scene rather than a feature of the already blurred, sampled and noisy image data. An analytical solution for the optimal characteristic image processor is developed. The Wiener filter for the sampled image case is obtained as a special case, where the desired characteristic is scene restoration. Optimal edge detection is investigated using the Laplacian operator x G as the desired characteristic, where G is a two dimensional Gaussian distribution function. It is shown that the optimal edge detector compensates for the blurring introduced by the image gathering optics, and notably, that it is not circularly symmetric. The lack of circular symmetry is largely due to the geometric effects of the sampling lattice used in image acquisition. The optimal image gathering optical transfer function is also investigated and the results of a sensitivity analysis are shown.

Halyo, Nesim↗

An application of fuzzy sets in real time filtering problems

The human decision-making task is modeled with fuzzy sets and the Kalman filter updates to the state vector are weighted using fuzzy functions. Results of the study show that the use of fuzzy set models gives results comparable to those requiring human assistance, e.g. lock on to false targets, and better results when the problem is caused by noise and/or bias, or unexpected errors in the state vector at initial acquisition.

Lea, Robert N.↗

A compressed marine data set for geomagnetic field modeling

Some 13 million scalar magnetic field data points that have been collected from the world's ocean areas reside in the collection of the National Geophysical Data Center. In order to derive a suitable data set for modeling the geomagnetic field of the earth, each ship track is divided into 220 km segments. The distribution of the reduced data in position, time and local time is discussed. The along-track filtering process described has proved to be an effective method of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for field modeling. This process also provides the benefits of smoothing short-wavelength crystal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be utilized in the least squares modeling.

Langel, R. A.↗

Open-Loop Acquisition Of Frequency In BPSK

Proposed open-loop analog/digital signal-processing system would be Costas-loop error detector functioning in closed-loop manner overall. Detector estimates difference between frequency of input signal and internal reference oscillator. Estimate used to close frequency-control loop. Precise symbol timing not necessary. Performance better than systems that effect open-loop acquisition using integrators instead of low-pass filters in arms of Costas loops and in which performance varies with symbol timing.

Shah, Biren N.↗

Observations of the Neuton Star RXJ0002.9 + 6246

The High Resolution Imager (HRI) observation of the CTB1 region requested in the original proposal was performed and data received. The primary target, the candidate neutron star, was detected at the predicted intensity. As proposed, both timing and spatial extent analysis has been performed on the extracted photons. The preliminary timing analysis showed no evidence for the faint pulsations seen in the PSPC data. As the pulse fraction was rather small, this was not an unexpected result. We have, however, recently developed more sophisticated timing analysis capabilities and are in the process of reanalyzing the data to more accurately constrain the pulse fraction, if indeed any evidence for pulsation is found. We have supplemented the HRI observations with optical observations. Filter images, in BVR & I, have been obtained as well as spectra for candidate counterparts. As of this writing, the original analysis holds, i.e. the x-ray flux is not consistent with the L(sub x)/L(sub opt) ratios expected from the candidate stars within the 3 sigma positional error circles. HRI confirmation of the flux as well as the refinement of the position have allowed us to limit the candidate list and refine the original L(sub x)/L(sub opt) analysis.

Hailey, Charles J.↗

Optimized FPGA Implementation of Multi-Rate FIR Filters Through Thread Decomposition

Multirate (decimation/interpolation) filters are among the essential signal processing components in spaceborne instruments where Finite Impulse Response (FIR) filters are often used to minimize nonlinear group delay and finite-precision effects. Cascaded (multi-stage) designs of Multi-Rate FIR (MRFIR) filters are further used for large rate change ratio, in order to lower the required throughput while simultaneously achieving comparable or better performance than single-stage designs. Traditional representation and implementation of MRFIR employ polyphase decomposition of the original filter structure, whose main purpose is to compute only the needed output at the lowest possible sampling rate. In this paper, an alternative representation and implementation technique, called TD-MRFIR (Thread Decomposition MRFIR), is presented. The basic idea is to decompose MRFIR into output computational threads, in contrast to a structural decomposition of the original filter as done in the polyphase decomposition. Each thread represents an instance of the finite convolution required to produce a single output of the MRFIR. The filter is thus viewed as a finite collection of concurrent threads. The technical details of TD-MRFIR will be explained, first showing its applicability to the implementation of downsampling, upsampling, and resampling FIR filters, and then describing a general strategy to optimally allocate the number of filter taps. A particular FPGA design of multi-stage TD-MRFIR for the L-band radar of NASA's SMAP (Soil Moisture Active Passive) instrument is demonstrated; and its implementation results in several targeted FPGA devices are summarized in terms of the functional (bit width, fixed-point error) and performance (time closure, resource usage, and power estimation) parameters.

DSP↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA PACE Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗