Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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 73 records · Page 4

Cloud detection using satellite measurements of infrared and visible radiances for ISCCP

This paper, the first of three, describes the cloud detection part of the International Satellite Cloud Climatology Project (ISCCP) analysis. Key features of the cloud detection alogrithm are (1) use of space and time radiance variation tests over several different space and time domains to account for the global variety of cloudy and clear characteristics, (2) estimation of clear radiance values for every time and place, and, (3) use of radiance thresholds that vary with the type of surface and climate regime. Design of the detection algorithm was supported by global, multiyear surveys of the statistical behavior of satellite-measured infrared and visible radiances to determine those characteristics that differentiate cloudy and clear scenes and how these characteristics vary among climate regimes. A summary of these statistical results is presented to illustrate how the cloud detection method works in a variety of circumstances. The sensitivity of the results to changing test parameter values is determined to provide a first estimate of the uncertainty of ISCCP cloud amounts. These test results (which exclude polar regions) suggest detection uncertainties of about 10% with possible negative biases of 5% (especially at night).

Rossow, Illiam B.↗

Target Detection via Cognitive Radars Using Change-Point Detection, Learning, and Adaptation

Many radar detection algorithms that assume a stationary environment (clutter) have been proposed and analyzed over the years. However, in practice, changes in the nonstationary environment can perturb the parameters of the clutter distribution, or even alter the clutter distribution family, which can greatly deteriorate the target detection capability. To avoid such potential performance degradation, cognitive radar systems are envisioned which are required to rapidly realize the nonstationarity, accurately learn the new characteristics of the environments, and adaptively update the detector. In this paper, aiming to develop a fully cognitive radar for target detection in nonstationary environments, we propose a unifying framework that integrates (i) change-point detection of clutter distributions by using a data-driven cumulative sum (CUSUM) algorithm and its extended version, (ii) learning/identification of clutter distribution by applying sparse theory and kernel density estimation methods, and (iii) adaptive target detection by automatically modifying the likelihood-ratio test and corresponding detection threshold. Further, with extensive numerical examples, we demonstrate the achieved improvements in detection performance due to the proposed framework in comparison to a nonadaptive case, an adaptive matched filter (AMF) method, and the clairvoyant case. Herein, we also use Wilcoxon rank-sum tests to evaluate the statistical significance of the performance improvements

42 ENGINEERING↗

When Will the Antarctic Ozone Hole Recover?

The Antarctic ozone hole demonstrates large-scale, man-made affects on our atmosphere. Surface observations now show that human produced ozone depleting substances (ODSs) are declining. The ozone hole should soon start to diminish because of this decline. In this talk we will demonstrate an ozone hole parametric model. This model is based upon: 1) a new algorithm for estimating 61 and Br levels over Antarctica and 2) late-spring Antarctic stratospheric temperatures. This parametric model explains 95% of the ozone hole area's variance. We use future ODS levels to predict ozone hole recovery. Full recovery to 1980 levels will occur in approximately 2068. The ozone hole area will very slowly decline over the next 2 decades. Detection of a statistically significant decrease of area will not occur until approximately 2024. We further show that nominal Antarctic stratospheric greenhouse gas forced temperature change should have a small impact on the ozone hole.

Newman, Paul A.↗

When Will the Antarctic Ozone Hole Recover?

The Antarctic ozone hole demonstrates large-scale, man-made affects on our atmosphere. Surface observations now show that human produced ozone depleting substances (ODSs) are declining. The ozone hole should soon start to diminish because of this decline. Herein we demonstrate an ozone hole parametric model. This model is based upon: 1) a new algorithm for estimating C1 and Br levels over Antarctica and 2) late-spring Antarctic stratospheric temperatures. This parametric model explains 95% of the ozone hole area s variance. We use future ODS levels to predict ozone hole recovery. Full recovery to 1980 levels will occur in approximately 2068. The ozone hole area will very slowly decline over the next 2 decades. Detection of a statistically significant decrease of area will not occur until approximately 2024. We further show that nominal Antarctic stratospheric greenhouse gas forced temperature change should have a small impact on the ozone hole.

Newman, Paul A.↗

Audibility of Multiple, Low-Frequency Tonal Signals in Noise

The main purpose of this study is to examine the audibility of multiple, low-frequency tones that are placed in distinct auditory channels. Three experiments are described, the goals of which are to determine if the presence of sound in multiple channels results in enhanced audibility and to assess the applicability of the Statistical Summation Model (SSM) to this frequency range. This model predicts that for the case of multiple signals that are in separate auditory channels, implying statistical independence, each with sensitivity value d prime of i, the resulting total sensitivity is given by the square root of the sum of the squares of the individual d prime of i values. In common with previous studies conducted at higher frequencies, the signals are pure tones and the maskers are broadband noise. The requirement that low frequency tones be placed in separate auditory filters limited the number of tones to a maximum of three. The first of the three experiments measured the change in masked thresholds for two- and three-tone signals relative to the level of the equally-detectable single tones. The multiple tone signals were composed of combinations of 55, 120 and 200 Hz tones. The measured changes in thresholds exceeded those predicted by the SSM, although they did not differ statistically from the model predictions. The second experiment employed the same overall approach but acquired more data and concentrated on the three-tone signal. Once again, the measured changes in masked threshold exceeded the model predictions, this time to a statistically-significant degree. Two issues were postulated with the potential to yield inflated changes in masked threshold: interaction between tones resulting in perceptible intermodulation/difference tones, and the assumption that the tones were in distinct auditory filters and statistically independent of one another. The third experiment used two sets of three-tone signals to address these latter concerns. The first set of three tones was composed of harmonically related tone frequencies of 55, 110 and 165 Hz, which was an attempt to reduce effects of intermodulation difference tones. The second set of three tones was chosen to be 110, 220 and 330 Hz, again reducing effects of difference tones, but also providing greater separation between tones. Results for the first set of three tones compared to those of the earlier experiments indicated that intermodulation was not an important effect. The second set of three tones (110, 220, 330 Hz) yielded changes in masked thresholds that, on average, were in good agreement with the SSM, although intersubject variability was large and prohibited a definitive conclusion regarding the concern that tone spacing was inadequate. The results of the three experiments showed that the masked threshold of sounds with multiple (two or three) equally-detectable low frequency tones was lower than those of the single tones. In other words, it is clear that audibility is enhanced by the presence of signals in multiple auditory filters. This finding is consistent with most previous research conducted at higher frequencies. In contrast with previous research, test subjects were, on average, able to detect multitone sounds at lower levels than those predicted using the SSM. Analyses that included Monte Carlo simulations showed that normally distributed errors in the single tone thresholds result in biased estimates of the thresholds of multitone sounds. This phenomenon is likely responsible for at least a substantial fraction of the unexpected deviation of measurements from SSM predictions.

Rafaelof, Menachem↗

In Situ Investigation of Upper Airway Occlusion in Sleep Disordered Breathing Using Ultrasonic Transducer Arrays

This work presents a novel application of ultrasound for the real-time, non-invasive investigation of occlusion of the upper airway during events of obstructive sleep apnea/hypopnea syndrome. It is hypothesized that ultrasonic pulses applied to the neck during apneic events produce spectral and temporal features that can detect apnea occurrence. Theoretical models of ultrasound propagation and an in vitro test were conducted to test this hypothesis in both transmission and reflection modes. Complete specifications and technical details of the system design and fabrication, which is mounted on each subject’s neck, are presented, including the methodology. Nine patients (seven male and two female, mean age of 42 years, with a range of 25 to 56 years, and body mass index 37.6 ± 6.6 kg/m2) were recruited for a full night study, which included simultaneous nocturnal polysomnography for the validation of the results. Nine temporal features and four spectral features were extracted from the envelope of the received pulse waveform. These were used to compute 26 metrics to quantify the changes in the ultrasonic waveforms between normal breathing and apneic events. The statistical analysis of the collected ultrasonic data showed that at least two or more of the proposed features could detect apneic events in all subjects. The findings establish the feasibility of the proposed method as a cost-effective and non-invasive OSAHS screening tool.

60 APPLIED LIFE SCIENCES↗

Compression of digital images over local area networks. Appendix 1: Item 3

Differential Pulse Code Modulation (DPCM) has been used with speech for many years. It has not been as successful for images because of poor edge performance. The only corruption in DPC is quantizer error but this corruption becomes quite large in the region of an edge because of the abrupt changes in the statistics of the signal. We introduce two improved DPCM schemes; Edge correcting DPCM and Edge Preservation Differential Coding. These two coding schemes will detect the edges and take action to correct them. In an Edge Correcting scheme, the quantizer error for an edge is encoded using a recursive quantizer with entropy coding and sent to the receiver as side information. In an Edge Preserving scheme, when the quantizer input falls in the overload region, the quantizer error is encoded and sent to the receiver repeatedly until the quantizer input falls in the inner levels. Therefore these coding schemes increase the bit rate in the region of an edge and require variable rate channels. We implement these two variable rate coding schemes on a token wing network. Timed token protocol supports two classes of messages; asynchronous and synchronous. The synchronous class provides a pre-allocated bandwidth and guaranteed response time. The remaining bandwidth is dynamically allocated to the asynchronous class. The Edge Correcting DPCM is simulated by considering the edge information under the asynchronous class. For the simulation of the Edge Preserving scheme, the amount of information sent each time is fixed, but the length of the packet or the bit rate for that packet is chosen depending on the availability capacity. The performance of the network, and the performance of the image coding algorithms, is studied.

Gorjala, Bhargavi↗

Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology

Changing trends in ecosystem productivity can be quantified using satellite observations of Normalized Difference Vegetation Index (NDVI). However, the estimation of trends from NDVI time series differs substantially depending on analyzed satellite dataset, the corresponding spatiotemporal resolution, and the applied statistical method. Here we compare the performance of a wide range of trend estimation methods and demonstrate that performance decreases with increasing inter-annual variability in the NDVI time series. Trend slope estimates based on annual aggregated time series or based on a seasonal-trend model show better performances than methods that remove the seasonal cycle of the time series. A breakpoint detection analysis reveals that an overestimation of breakpoints in NDVI trends can result in wrong or even opposite trend estimates. Based on our results, we give practical recommendations for the application of trend methods on long-term NDVI time series. Particularly, we apply and compare different methods on NDVI time series in Alaska, where both greening and browning trends have been previously observed. Here, the multi-method uncertainty of NDVI trends is quantified through the application of the different trend estimation methods. Our results indicate that greening NDVI trends in Alaska are more spatially and temporally prevalent than browning trends. We also show that detected breakpoints in NDVI trends tend to coincide with large fires. Overall, our analyses demonstrate that seasonal trend methods need to be improved against inter-annual variability to quantify changing trends in ecosystem productivity with higher accuracy.

Inter-Annual↗

Detecting hidden transient events in noisy nonlinear time-series

The information impulse function (IIF), running Variance, and local Hölder Exponent are three conceptually different time-series evaluation techniques. These techniques examine time-series for local changes in information content, statistical variation, and point-wise smoothness, respectively. Using simulated data emulating a randomly excited nonlinear dynamical system, this study interrogates the utility of each method to correctly differentiate a transient event from the background while simultaneously locating it in time. Computational experiments are designed and conducted to evaluate the efficacy of each technique by varying pulse size, time location, and noise level in time-series. Our findings reveal that, in most cases, the first instance of a transient event is more easily observed with the information-based approach of IIF than with the Variance and local Hölder Exponent methods. While our study highlights the unique strengths of each technique, the results suggest that very robust and reliable event detection for nonlinear systems producing noisy time-series data can be obtained by incorporating the IIF into the analysis.

42 ENGINEERING↗

Identification of Human-Induced Changes in Atmospheric Moisture Content

Data from the satellite-based Special Sensor Microwave Imager (SSM/I) show that the total atmospheric moisture content over oceans has increased by 0.41 kg/sq m per decade since 1988. Results from current climate models indicate that water vapor increases of this magnitude cannot be explained by climate noise alone. In a formal detection and attribution analysis using the pooled results from 22 different climate models, the simulated "fingerprint" pattern of anthropogenically caused changes in water vapor is identifiable with high statistical confidence in the SSM/I data. Experiments in which forcing factors are varied individually suggest that this fingerprint "match" is primarily due to human-caused increases in greenhouse gases and not to solar forcing or recovery from the eruption of Mount Pinatubo. Our findings provide preliminary evidence of an emerging anthropogenic signal in the moisture content of earth's atmosphere.

Santer, B.D.↗

Passive/active Microwave Soil Moisture Change Disaggregation Using Smapvex12 Data

The SMAPVEX12 (Soil Moisture Active Passive (SMAP) Validation Experiment 2012) experiment was conducted during June-July 2012 in Manitoba, Canada with the goal of collecting remote sensing data and ground measurements for the development and testing of soil moisture retrieval algorithms under varying vegetation and soil conditions for the SMAP satellite. The aircraft based soil moisture data provided by the passive/active microwave sensor PALS (Passive and Active L-band System) has a nominal spatial resolution of 1600 m. However, this resolution is not compatible with agricultural, meteorological and hydrological studies that require high spatial resolutions and this issue can be solved by soil moisture disaggregation. The soil moisture disaggregation algorithm integrates radiometer soil moisture retrievals and high-resolution radar observations and it can provide soil moisture estimates at a finer scale than the radiometer data alone. In this study, a change detection algorithm was used for disaggregation of coarse resolution passive microwave soil moisture retrievals with radar backscatter coefficients obtained from the higher spatial resolution UAVSAR (Unmanned Air Vehicle Synthetic Aperture Radar) at crop field scale. The accuracy of the disaggregated change in soil moisture was evaluated using ground based soil moisture measurements collected during SMAPVEX12 campaign. The results showed that soil moisture spatial variabilities were better characterized by the disaggregated change in soil moisture estimates at 5 m / 800 m resolution as well as good agreement with in situ measurements. It also showed that VWC (Vegetation Water Content) did not have a big impact on disaggregation algorithm performance, with R2 of the disaggregated results ranging 0.628-0.794. The 5 m and 800m resolution disaggregated soil moisture did no show significant difference in statistical performance variables.

Radar Backscatter↗

Robust Measurement of Stellar Streams around the Milky Way: Correcting Spatially Variable Observational Selection Effects in Optical Imaging Surveys

Observations of density variations in stellar streams are a promising probe of low-mass dark matter substructure in the Milky Way. However, survey systematics such as variations in seeing and sky brightness can also induce artificial fluctuations in the observed densities of known stellar streams. These variations arise because survey conditions affect both object detection and star–galaxy misclassification rates. To mitigate these effects, we use Balrog synthetic source injections in the Dark Energy Survey (DES) Y3 data to calculate detection rate variations and classification rates as functions of survey properties. We show that these rates are nearly separable with respect to survey properties and can be estimated with sufficient statistics from the synthetic catalogs. Applying these corrections reduces the standard deviation of relative detection rates across the DES footprint by a factor of 5, and our corrections significantly change the inferred linear density of the Phoenix stream when including faint objects. Additionally, for artificial streams with DES-like survey properties we are able to recover density power spectra with reduced bias. We also find that uncorrected power-spectrum results for Legacy Survey of Space and Time (LSST)-like data can be around 5 times more biased, highlighting the need for such corrections in future ground-based surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Ensembles of radial basis function networks for spectroscopic detection of cervical precancer

The mortality related to cervical cancer can be substantially reduced through early detection and treatment. However, current detection techniques, such as Pap smear and colposcopy, fail to achieve a concurrently high sensitivity and specificity. In vivo fluorescence spectroscopy is a technique which quickly, noninvasively and quantitatively probes the biochemical and morphological changes that occur in precancerous tissue. A multivariate statistical algorithm was used to extract clinically useful information from tissue spectra acquired from 361 cervical sites from 95 patients at 337-, 380-, and 460-nm excitation wavelengths. The multivariate statistical analysis was also employed to reduce the number of fluorescence excitation-emission wavelength pairs required to discriminate healthy tissue samples from precancerous tissue samples. The use of connectionist methods such as multilayered perceptrons, radial basis function (RBF) networks, and ensembles of such networks was investigated. RBF ensemble algorithms based on fluorescence spectra potentially provide automated and near real-time implementation of precancer detection in the hands of nonexperts. The results are more reliable, direct, and accurate than those achieved by either human experts or multivariate statistical algorithms.

Cervix Neoplasms/pathology↗

Seismic Excitation of the Polar Motion, 1977-1993

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed by CHAO and GROSS (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0-1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards approximately 140deg E, away from the actual observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field by CHAO and GROSS (1987), manifests some geodynamic behavior yet to be explored.

Chao, Benjamin Fong↗

Impact of Glass Irradiation on Laser-Induced Breakdown Spectroscopy Data Analysis

Increased absorption of optical materials arising from exposure to ionizing radiation must be accounted for to accurately analyze laser-induced breakdown spectroscopy (LIBS) data retrieved from high-radiation environments. We evaluate this effect on two examples that mimic the diagnostics placed within novel nuclear reactor designs. The analysis is performed on LIBS data measured with 1% Xe gas in an ambient He environment and 1% Eu in a molten LiCl-KCl matrix, along with the measured optical absorption from the gamma- and neutron-irradiated low-OH fused silica and sapphire glasses. Significant changes in the number of laser shots required to reach a 3σ detection level are observed for the Eu data, increasing by two orders of magnitude after exposure to a 1.7 × 1017 n/cm2 neutron fluence. For all cases examined, the spectral dependence of absorption results in the introduction of systematic errors. Moreover, if lines from different spectral regions are used to create Boltzmann plots, this attenuation leads to statistically significant changes in the temperatures calculated from the Xe II lines and Eu II lines, lowering them from 8000 ± 610 K to 6900 ± 810 K and from 15,800 ± 400 K to 7200 ± 800 K, respectively, for exposure to the 1.7 × 1017 n/cm2 fluence. The temperature range required for a 95% confidence interval for the calculated temperature is also broadened. In the case of measuring the Xe spectrum, these effects may be mitigated using only the longer-wavelength spectral region, where radiation attenuation is relatively small, or through analysis using the iterative Saha–Boltzmann method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Seismic Excitation of the Polar Motion

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed by Chao and Gross (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0-1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards approx. 140 deg E, away from the actually observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field by Chao and Gross (1987), manifests some geodynamic behavior yet to be explored.

Chao, Benjamin Fong↗

Variability in the power spectrum of solar five-minute oscillations

It is noted that two-dimensional power spectra of solar five-minute oscillations exhibit prominent ridge structures in (k, omega) space, k being the horizontal wavenumber and omega the temporal frequency. The positions of these ridges in k and omega can be used in probing temperature and velocity structures in the subphotosphere. The progress made in a program to establish whether power spectra taken on separate days show shifts in ridge locations is surveyed, noting that these may arise from different velocity and temperature patterns having been brought into the sampling region by solar rotation. Power spectra are obtained for six days of observations of Doppler velocities using the Mg I 5173 A and Fe I 5434 A spectral lines. Each data set encompasses 8 to 11 hr in time and samples a region 256 x 1024 arcsec in spatial extent, with a spatial resolution of 2 arcsec and temporal sampling of 65 s. Statistically significant shifts are detected in ridge locations between certain data sets. When analyzed in terms of eastward and westward propagating waves, the character of these displacements is seen as implying that changes have occurred in both temperature and horizontal velocity fields underlying the observing window. The velocity changes are estimated to be of the order of 100 m/s.

Hill, F.↗

Downsampling Photodetector Array with Windowing

In a photon counting detector array, each pixel in the array produces an electrical pulse when an incident photon on that pixel is detected. Detection and demodulation of an optical communication signal that modulated the intensity of the optical signal requires counting the number of photon arrivals over a given interval. As the size of photon counting photodetector arrays increases, parallel processing of all the pixels exceeds the resources available in current application-specific integrated circuit (ASIC) and gate array (GA) technology; the desire for a high fill factor in avalanche photodiode (APD) detector arrays also precludes this. Through the use of downsampling and windowing portions of the detector array, the processing is distributed between the ASIC and GA. This allows demodulation of the optical communication signal incident on a large photon counting detector array, as well as providing architecture amenable to algorithmic changes. The detector array readout ASIC functions as a parallel-to-serial converter, serializing the photodetector array output for subsequent processing. Additional downsampling functionality for each pixel is added to this ASIC. Due to the large number of pixels in the array, the readout time of the entire photodetector is greater than the time between photon arrivals; therefore, a downsampling pre-processing step is done in order to increase the time allowed for the readout to occur. Each pixel drives a small counter that is incremented at every detected photon arrival or, equivalently, the charge in a storage capacitor is incremented. At the end of a user-configurable counting period (calculated independently from the ASIC), the counters are sampled and cleared. This downsampled photon count information is then sent one counter word at a time to the GA. For a large array, processing even the downsampled pixel counts exceeds the capabilities of the GA. Windowing of the array, whereby several subsets of pixels are designated for processing, is used to further reduce the computational requirements. The grouping of the designated pixel frame as the photon count information is sent one word at a time to the GA, the aggregation of the pixels in a window can be achieved by selecting only the designated pixel counts from the serial stream of photon counts, thereby obviating the need to store the entire frame of pixel count in the gate array. The pixel count se quence from each window can then be processed, forming lower-rate pixel statistics for each window. By having this processing occur in the GA rather than in the ASIC, future changes to the processing algorithm can be readily implemented. The high-bandwidth requirements of a photon counting array combined with the properties of the optical modulation being detected by the array present a unique problem that has not been addressed by current CCD or CMOS sensor array solutions.

Patawaran, Ferze D.↗