Taylor, Brian A; Hwang, Ken-Pin; Hazle, John D; Stafford, R Jason
2009-03-01
The authors investigated the performance of the iterative Steiglitz-McBride (SM) algorithm on an autoregressive moving average (ARMA) model of signals from a fast, sparsely sampled, multiecho, chemical shift imaging (CSI) acquisition using simulation, phantom, ex vivo, and in vivo experiments with a focus on its potential usage in magnetic resonance (MR)-guided interventions. The ARMA signal model facilitated a rapid calculation of the chemical shift, apparent spin-spin relaxation time (T2*), and complex amplitudes of a multipeak system from a limited number of echoes (< or equal 16). Numerical simulations of one- and two-peak systems were used to assess the accuracy and uncertainty in the calculated spectral parameters as a function of acquisition and tissue parameters. The measured uncertainties from simulation were compared to the theoretical Cramer-Rao lower bound (CRLB) for the acquisition. Measurements made in phantoms were used to validate the T2* estimates and to validate uncertainty estimates made from the CRLB. We demonstrated application to real-time MR-guided interventions ex vivo by using the technique to monitor a percutaneous ethanol injection into a bovine liver and in vivo to monitor a laser-induced thermal therapy treatment in a canine brain. Simulation results showed that the chemical shift and amplitude uncertainties reached their respective CRLB at a signal-to-noise ratio (SNR) > or =5 for echo train lengths (ETLs) > or =4 using a fixed echo spacing of 3.3 ms. T2* estimates from the signal model possessed higher uncertainties but reached the CRLB at larger SNRs and/or ETLs. Highly accurate estimates for the chemical shift (<0.01 ppm) and amplitude (<1.0%) were obtained with > or =4 echoes and for T2*(<1.0%) with > or =7 echoes. We conclude that, over a reasonable range of SNR, the SM algorithm is a robust estimator of spectral parameters from fast CSI acquisitions that acquire < or =16 echoes for one- and two-peak systems. Preliminary ex vivo
Taylor, Brian A.; Hwang, Ken-Pin; Hazle, John D.; Stafford, R. Jason
2009-01-01
The authors investigated the performance of the iterative Steiglitz–McBride (SM) algorithm on an autoregressive moving average (ARMA) model of signals from a fast, sparsely sampled, multiecho, chemical shift imaging (CSI) acquisition using simulation, phantom, ex vivo, and in vivo experiments with a focus on its potential usage in magnetic resonance (MR)-guided interventions. The ARMA signal model facilitated a rapid calculation of the chemical shift, apparent spin-spin relaxation time (T2*), and complex amplitudes of a multipeak system from a limited number of echoes (≤16). Numerical simulations of one- and two-peak systems were used to assess the accuracy and uncertainty in the calculated spectral parameters as a function of acquisition and tissue parameters. The measured uncertainties from simulation were compared to the theoretical Cramer–Rao lower bound (CRLB) for the acquisition. Measurements made in phantoms were used to validate the T2* estimates and to validate uncertainty estimates made from the CRLB. We demonstrated application to real-time MR-guided interventions ex vivo by using the technique to monitor a percutaneous ethanol injection into a bovine liver and in vivo to monitor a laser-induced thermal therapy treatment in a canine brain. Simulation results showed that the chemical shift and amplitude uncertainties reached their respective CRLB at a signal-to-noise ratio (SNR)≥5 for echo train lengths (ETLs)≥4 using a fixed echo spacing of 3.3 ms. T2* estimates from the signal model possessed higher uncertainties but reached the CRLB at larger SNRs and∕or ETLs. Highly accurate estimates for the chemical shift (<0.01 ppm) and amplitude (<1.0%) were obtained with ≥4 echoes and for T2* (<1.0%) with ≥7 echoes. We conclude that, over a reasonable range of SNR, the SM algorithm is a robust estimator of spectral parameters from fast CSI acquisitions that acquire ≤16 echoes for one- and two-peak systems. Preliminary ex vivo and in vivo
NASA Astrophysics Data System (ADS)
Bousi, Evgenia; Pitris, Costas
2015-03-01
Fourier Domain (FD) Optical Coherence Tomography (OCT) interferograms require a Fourier transformation in order to be converted to A-Scans representing the backscattering intensity from the different depths of the tissue microstructure. Most often, this transformation is performed using a discrete Fourier transform, i.e. the well-known Fast Fourier Transform (FFT). However, there are many alternatives for performing the necessary spectral conversion. Autoregressive (AR) spectral estimation techniques are one such example. The parametric nature of AR techniques offers several advantages, compared to the commonly-used FFT, including better convergence and less susceptibility to noise. They can also be adjusted to represent more or less of the signal detail depending on the order of the autoregression. These features make them uniquely suited for processing the FD OCT data. The advantages of the proposed methodology are illustrated on in vivo skin imaging data and the resolution is verified on single back-reflections from a glass surface. AR spectral estimation can be used to convert the interferograms to A-Scans while at the same time reducing the artifacts caused by high intensity back-reflections (by -20 dB) and diminishing the speckle (by -12 dB) all without the degradation in the resolution associated with other techniques.
Application of multivariate autoregressive spectrum estimation to ULF waves
NASA Technical Reports Server (NTRS)
Ioannidis, G. A.
1975-01-01
The estimation of the power spectrum of a time series by fitting a finite autoregressive model to the data has recently found widespread application in the physical sciences. The extension of this method to the analysis of vector time series is presented here through its application to ULF waves observed in the magnetosphere by the ATS 6 synchronous satellite. Autoregressive spectral estimates of the power and cross-power spectra of these waves are computed with computer programs developed by the author and are compared with the corresponding Blackman-Tukey spectral estimates. The resulting spectral density matrices are then analyzed to determine the direction of propagation and polarization of the observed waves.
Chen, L. Leon; Madhavan, Radhika; Rapoport, Benjamin I.; Anderson, William S.
2012-01-01
Neural oscillations are important features in a working central nervous system, facilitating efficient communication across large networks of neurons. They are implicated in a diverse range of processes such as synchronization and synaptic plasticity, and can be seen in a variety of cognitive processes. For example, hippocampal theta oscillations are thought to be a crucial component of memory encoding and retrieval. To better study the role of these oscillations in various cognitive processes, and to be able to build clinical applications around them, accurate and precise estimations of the instantaneous frequency and phase are required. Here, we present methodology based on autoregressive modeling to accomplish this in real time. This allows the targeting of stimulation to a specific phase of a detected oscillation. We first assess performance of the algorithm on two signals where the exact phase and frequency are known. Then, using intracranial EEG recorded from two patients performing a Sternberg memory task, we characterize our algorithm’s phase-locking performance on physiologic theta oscillations: optimizing algorithm parameters on the first patient using a genetic algorithm, we carried out cross-validation procedures on subsequent trials and electrodes within the same patient, as well as on data recorded from the second patient. PMID:21292589
MAXIMUM LIKELIHOOD ESTIMATION FOR PERIODIC AUTOREGRESSIVE MOVING AVERAGE MODELS.
Vecchia, A.V.
1985-01-01
A useful class of models for seasonal time series that cannot be filtered or standardized to achieve second-order stationarity is that of periodic autoregressive moving average (PARMA) models, which are extensions of ARMA models that allow periodic (seasonal) parameters. An approximation to the exact likelihood for Gaussian PARMA processes is developed, and a straightforward algorithm for its maximization is presented. The algorithm is tested on several periodic ARMA(1, 1) models through simulation studies and is compared to moment estimation via the seasonal Yule-Walker equations. Applicability of the technique is demonstrated through an analysis of a seasonal stream-flow series from the Rio Caroni River in Venezuela.
NASA Astrophysics Data System (ADS)
Tary, J. B.; Herrera, R. H.; van der Baan, M.
2014-01-01
Recent studies show that the frequency content of continuous passive recordings contains useful information for the study of hydraulic fracturing experiments as well as longstanding applications in volcano and global seismology. The short-time Fourier transform (STFT) is usually used to obtain the time-frequency representation of a seismic trace. Yet, this transform has two main disadvantages, namely its fixed time-frequency resolution and spectral leakage. Here, we describe two methods based on autoregressive (AR) models: the short-time autoregressive method (ST-AR) and the Kalman smoother (KS). These two methods allow for the AR coefficients to vary over time in order to follow time-varying frequency contents. The outcome of AR methods depends mainly on the number of AR coefficients. We use a robust approach to estimate the optimum order of the AR methods that best matches the spectral comparison between Fourier and AR spectra. Comparing the outcomes of the three methods on a synthetic signal, a long-period volcanic event, and microseismic data, we show that the STFT and both AR methods are able to track fast changes in frequency content. The STFT provides reasonable results even for noisy data using a simple and effective algorithm. The coefficients of the AR filter are defined at all time in the case of the KS. However, its better time resolution is slightly offset by a lower frequency resolution and its computational complexity. The ST-AR has a high spectral resolution and the lowest sensitivity to background noises, facilitating the identification of the various frequency components. The estimated AR coefficients can also be used to extract parts of the signal. The study of long-term phenomena, such as resonance frequencies, or transient events, such as long-period events, could help to gain further insight on reservoir deformation during hydraulic fracturing experiments as well as global or volcano seismological signals.
Least squares estimation of Generalized Space Time AutoRegressive (GSTAR) model and its properties
NASA Astrophysics Data System (ADS)
Ruchjana, Budi Nurani; Borovkova, Svetlana A.; Lopuhaa, H. P.
2012-05-01
In this paper we studied a least squares estimation parameters of the Generalized Space Time AutoRegressive (GSTAR) model and its properties, especially in consistency and asymptotic normality. We use R software to estimate the GSTAR parameter and apply the model toward real phenomena of data, such as an oil production data at volcanic layer.
Piazza, C; Cantiani, C; Tacchino, G; Molteni, M; Reni, G; Bianchi, A M
2014-01-01
The ability to process rapidly-occurring auditory stimuli plays an important role in the mechanisms of language acquisition. For this reason, the research community has begun to investigate infant auditory processing, particularly using the Event Related Potentials (ERP) technique. In this paper we approach this issue by means of time domain and time-frequency domain analysis. For the latter, we propose the use of Adaptive Autoregressive (AAR) identification with spectral power decomposition. Results show EEG delta-theta oscillation enhancement related to the processing of acoustic frequency and duration changes, suggesting that, as expected, power modulation encodes rapid auditory processing (RAP) in infants and that the time-frequency analysis method proposed is able to identify this modulation. PMID:25571014
Spectral estimators in elastography.
Konofagou, E E; Varghese, T; Ophir, J
2000-03-01
Like velocity, strain induces a time delay and a time scaling to the received signal. Elastography typically uses time delay techniques to indirectly (i.e. via the displacement estimate) measure tissue strain induced by an applied compression, and considers time scaling as a source of distortion. More recently, we have shown that the time scaling factor can also be spectrally estimated and used as a direct measure of strain. Strain causes a Doppler-like frequency shift and a change in bandwidth of the bandpass power spectrum of the echo signal. Two frequency shift strain estimators are described that have been proven to be more robust but less precise when compared to time delay estimators, both in simulations and experiments. The increased robustness is due to the insensitivity of the spectral techniques to phase decorrelation noise. In this paper we discuss and compare the theoretical and experimental findings obtained with traditional time delay estimators and with the newly proposed spectral methods. PMID:10829698
Power spectral estimation algorithms
NASA Technical Reports Server (NTRS)
Bhatia, Manjit S.
1989-01-01
Algorithms to estimate the power spectrum using Maximum Entropy Methods were developed. These algorithms were coded in FORTRAN 77 and were implemented on the VAX 780. The important considerations in this analysis are: (1) resolution, i.e., how close in frequency two spectral components can be spaced and still be identified; (2) dynamic range, i.e., how small a spectral peak can be, relative to the largest, and still be observed in the spectra; and (3) variance, i.e., how accurate the estimate of the spectra is to the actual spectra. The application of the algorithms based on Maximum Entropy Methods to a variety of data shows that these criteria are met quite well. Additional work in this direction would help confirm the findings. All of the software developed was turned over to the technical monitor. A copy of a typical program is included. Some of the actual data and graphs used on this data are also included.
NASA Astrophysics Data System (ADS)
Yi, Guo-Sheng; Wang, Jiang; Deng, Bin; Wei, Xi-Le; Han, Chun-Xiao
2013-02-01
To investigate whether and how manual acupuncture (MA) modulates brain activities, we design an experiment where acupuncture at acupoint ST36 of the right leg is used to obtain electroencephalograph (EEG) signals in healthy subjects. We adopt the autoregressive (AR) Burg method to estimate the power spectrum of EEG signals and analyze the relative powers in delta (0 Hz-4 Hz), theta (4 Hz-8 Hz), alpha (8 Hz-13 Hz), and beta (13 Hz-30 Hz) bands. Our results show that MA at ST36 can significantly increase the EEG slow wave relative power (delta band) and reduce the fast wave relative powers (alpha and beta bands), while there are no statistical differences in theta band relative power between different acupuncture states. In order to quantify the ratio of slow to fast wave EEG activity, we compute the power ratio index. It is found that the MA can significantly increase the power ratio index, especially in frontal and central lobes. All the results highlight the modulation of brain activities with MA and may provide potential help for the clinical use of acupuncture. The proposed quantitative method of acupuncture signals may be further used to make MA more standardized.
Ueda, Masashi; Tomobe, Katsuma; Setoguchi, Keiichi; Endou, Akira
2002-02-15
The response of a sensor depends on its operating conditions, and thus it is desirable to develop an in-service method for response time estimation. The applicability of the autoregressive (AR) model for this purpose was examined in the case of the fuel subassembly outlet coolant thermocouples and the primary circuit electromagnetic flowmeter (EMF) of Monju, the prototype fast breeder reactor in Japan.The use of an AR model with exogenous input (ARX model) is possible when the physical variable to be sensed can be observed by an alternative means with a faster response time than that of the sensor in question. In the case of the subassembly outlet thermocouple, the temperature output from an eddy-current sensor, during pseudorandom reactor power variation, served as the exogenous input.In respect to the thermocouple response, AR and ARX modeling were shown to be applicable, and the transient responses thus derived agreed well with each other and with the results measured by means of a step change in sodium temperature. However, the primary circuit EMF response time, estimated using the AR model, decreased with increasing flow rate even when approaching the rated flow, demonstrating that the method was not completely applicable. Nevertheless, it can be concluded that the response is faster than that estimated in the rated condition.
Chon, K H; Cohen, R J; Holstein-Rathlou, N H
1997-01-01
A linear and nonlinear autoregressive moving average (ARMA) identification algorithm is developed for modeling time series data. The algorithm uses Laguerre expansion of kernals (LEK) to estimate Volterra-Wiener kernals. However, instead of estimating linear and nonlinear system dynamics via moving average models, as is the case for the Volterra-Wiener analysis, we propose an ARMA model-based approach. The proposed algorithm is essentially the same as LEK, but this algorithm is extended to include past values of the output as well. Thus, all of the advantages associated with using the Laguerre function remain with our algorithm; but, by extending the algorithm to the linear and nonlinear ARMA model, a significant reduction in the number of Laguerre functions can be made, compared with the Volterra-Wiener approach. This translates into a more compact system representation and makes the physiological interpretation of higher order kernels easier. Furthermore, simulation results show better performance of the proposed approach in estimating the system dynamics than LEK in certain cases, and it remains effective in the presence of significant additive measurement noise. PMID:9236985
Correlation analysis for long time series by robustly estimated autoregressive stochastic processes
NASA Astrophysics Data System (ADS)
Schuh, Wolf-Dieter; Brockmann, Jan-Martin; Kargoll, Boris
2015-04-01
Modern sensors and satellite missions deliver huge data sets and long time series of observations. These data sets have to be handled with care because of changing correlations, conspicuous data and possible outliers. Tailored concepts for data selection and robust techniques to estimate the correlation characteristics allow for a better/optimal exploitation of the information of these measurements. In this presentation we give an overview of standard techniques for estimating correlations occurring in long time series in the time domain as well as in the frequency domain. We discuss the pros and cons especially with the focus on the intensified occurrence of conspicuous data and outliers. We present a concept to classify the measurements and isolate conspicuous data. We propose to describe the varying correlation behavior of the measurement series by an autoregressive stochastic process and give some hints how to construct adaptive filters to decorrelate the measurement series and to handle the huge covariance matrices. As study object we use time series from gravity gradient data collected during the GOCE low orbit operation campaign (LOOC). Due to the low orbit these data from 13-Jun-2014 to 21-Oct-2014 have more or less the same potential to recover the Earth gravity field with the same accuracy than all the data from the rest of the entire mission. Therefore these data are extraordinarily valuable but hard to handle, because of conspicuous data due to maneuvers during the orbit lowering phases, overall increase in drag, saturation of ion thrusters and other (currently) unexplained effects.
Estimation of source function and medium response function by autoregressive method
NASA Astrophysics Data System (ADS)
Nair, G. J.
1983-04-01
By modelling seismograms as "low" and "high" order autoregressive (AR) processes, the source function and the medium response function are separated from a single channel seismogram. Akaike's final prediction error is used as a statistic to select the appropriate "low" and "high" AR order of the process. Case studies of synthetic data show that the recovered source and reflectivity functions compare very well with the input functions. Using this method, arrivals of the surface reflected P phases of five explosions from the Soviet region and of two earthquakes from Kamchatka, recorded at Gauribidanur Seismic Array (GBA), India, are identified. Certain features of the source and source region of these events are also inferred.
Spectral procedures for estimating crop biomass
Wanjura, D.F.; Hatfield, J.L.
1985-05-01
Spectral reflectance was measured semi-weekly and used to estimate leaf area and plant dry weight accumulation in cotton, soybeans, and sunflower. Integration of spectral crop growth cycle curves explained up to 95 and 91%, respectively, of the variation in cotton lint yield and dry weight. A theoretical relationship for dry weight accumulation, in which only intercepted radiation or intercepted radiation and solar energy to biomass conversion efficiency were spectrally estimated, explained 99 and 96%, respectively, of the observed plant dry weight variation of the three crops. These results demonstrate the feasibility of predicting crop biomass from spectral measurements collected frequently during the growing season. 15 references.
Sava, H; Durand, L G; Cloutier, G
1999-05-01
To achieve an accurate estimation of the instantaneous turbulent velocity fluctuations downstream of prosthetic heart valves in vivo, the variability of the spectral method used to measure the mean frequency shift of the Doppler signal (i.e. the Doppler velocity) should be minimised. This paper investigates the performance of various short-time spectral parametric methods such as the short-time Fourier transform, autoregressive modelling based on two different approaches, autoregressive moving average modelling based on the Steiglitz-McBride method, and Prony's spectral method. A simulated Doppler signal was used to evaluate the performance of the above mentioned spectral methods and Gaussian noise was added to obtain a set of signals with various signal-to-noise ratios. Two different parameters were used to evaluate the performance of each method in terms of variability and accurate matching of the theoretical Doppler mean instantaneous frequency variation within the cardiac cycle. Results show that autoregressive modelling outperforms the other investigated spectral techniques for window lengths varying between 1 and 10 ms. Among the autoregressive algorithms implemented, it is shown that the maximum entropy method based on a block data processing technique gives the best results for a signal-to-noise ratio of 20 dB. However, at 10 and 0 dB, the Levinson-Durbin algorithm surpasses the performance of the maximum entropy method. It is expected that the intrinsic variance of the spectral methods can be an important source of error for the estimation of the turbulence intensity. The range of this error varies from 0.38% to 24% depending on the parameters of the spectral method and the signal-to-noise ratio. PMID:10505377
Spectral moment estimation in MST radars
NASA Technical Reports Server (NTRS)
Woodman, R. F.
1983-01-01
Signal processing techniques used in Mesosphere-Stratosphere-Troposphere (MST) radars are reviewed. Techniques which produce good estimates of the total power, frequency shift, and spectral width of the radar power spectra are considered. Non-linear curve fitting, autocovariance, autocorrelation, covariance, and maximum likelihood estimators are discussed.
NASA Astrophysics Data System (ADS)
Foffani, Guglielmo; Bianchi, Anna M.; Priori, Alberto; Baselli, Giuseppe
2004-09-01
We propose a method that combines adaptive autoregressive (AAR) identification and spectral power decomposition for the study of movement-related spectral changes in scalp EEG signals and basal ganglia local field potentials (LFPs). This approach introduces the concept of movement-related poles, allowing one to study not only the classical event-related desynchronizations (ERD) and synchronizations (ERS), which correspond to modulations of power, but also event-related modulations of frequency. We applied the method to analyze movement-related EEG signals and LFPs contemporarily recorded from the sensorimotor cortex, the globus pallidus internus (GPi) and the subthalamic nucleus (STN) in a patient with Parkinson's disease who underwent stereotactic neurosurgery for the implant of deep brain stimulation (DBS) electrodes. In the AAR identification we compared the whale and the exponential forgetting factors, showing that the whale forgetting provides a better disturbance rejection and it is therefore more suitable to investigate movement-related brain activity. Movement-related power modulations were consistent with previous studies. In addition, movement-related frequency modulations were observed from both scalp EEG signals and basal ganglia LFPs. The method therefore represents an effective approach to the study of movement-related brain activity.
A two dimensional power spectral estimate for some nonstationary processes. M.S. Thesis
NASA Technical Reports Server (NTRS)
Smith, Gregory L.
1989-01-01
A two dimensional estimate for the power spectral density of a nonstationary process is being developed. The estimate will be applied to helicopter noise data which is clearly nonstationary. The acoustic pressure from the isolated main rotor and isolated tail rotor is known to be periodically correlated (PC) and the combined noise from the main and tail rotors is assumed to be correlation autoregressive (CAR). The results of this nonstationary analysis will be compared with the current method of assuming that the data is stationary and analyzing it as such. Another method of analysis is to introduce a random phase shift into the data as shown by Papoulis to produce a time history which can then be accurately modeled as stationary. This method will also be investigated for the helicopter data. A method used to determine the period of a PC process when the period is not know is discussed. The period of a PC process must be known in order to produce an accurate spectral representation for the process. The spectral estimate is developed. The bias and variability of the estimate are also discussed. Finally, the current method for analyzing nonstationary data is compared to that of using a two dimensional spectral representation. In addition, the method of phase shifting the data is examined.
Small Sample Properties of Bayesian Multivariate Autoregressive Time Series Models
ERIC Educational Resources Information Center
Price, Larry R.
2012-01-01
The aim of this study was to compare the small sample (N = 1, 3, 5, 10, 15) performance of a Bayesian multivariate vector autoregressive (BVAR-SEM) time series model relative to frequentist power and parameter estimation bias. A multivariate autoregressive model was developed based on correlated autoregressive time series vectors of varying…
Statistics of the Spectral Kurtosis Estimator
NASA Astrophysics Data System (ADS)
Nita, Gelu M.; Gary, Dale E.
2010-05-01
Spectral kurtosis (SK) is a statistical approach for detecting and removing radio-frequency interference (RFI) in radio astronomy data. In this article, the statistical properties of the SK estimator are investigated and all moments of its probability density function are analytically determined. These moments provide a means to determine the tail probabilities of the estimator that are essential to defining the thresholds for RFI discrimination. It is shown that, for a number of accumulated spectra M>=24, the first SK standard moments satisfy the conditions required by a Pearson type IV probability density function (pdf), which is shown to accurately reproduce the observed distributions. The cumulative function (CF) of the Pearson type IV is then found, in both analytical and numerical forms, suitable for accurate estimation of the tail probabilities of the SK estimator. This same framework is also shown to be applicable to the related time-domain kurtosis (TDK) estimator, whose pdf corresponds to Pearson type IV when the number of time-domain samples is M>=46. The pdf and CF also are determined for this case.
Alegana, Victor A.; Atkinson, Peter M.; Wright, Jim A.; Kamwi, Richard; Uusiku, Petrina; Katokele, Stark; Snow, Robert W.; Noor, Abdisalan M.
2013-01-01
As malaria transmission declines, it becomes increasingly important to monitor changes in malaria incidence rather than prevalence. Here, a spatio-temporal model was used to identify constituencies with high malaria incidence to guide malaria control. Malaria cases were assembled across all age groups along with several environmental covariates. A Bayesian conditional-autoregressive model was used to model the spatial and temporal variation of incidence after adjusting for test positivity rates and health facility utilisation. Of the 144,744 malaria cases recorded in Namibia in 2009, 134,851 were suspected and 9893 were parasitologically confirmed. The mean annual incidence based on the Bayesian model predictions was 13 cases per 1000 population with the highest incidence predicted for constituencies bordering Angola and Zambia. The smoothed maps of incidence highlight trends in disease incidence. For Namibia, the 2009 maps provide a baseline for monitoring the targets of pre-elimination. PMID:24238079
Alegana, Victor A; Atkinson, Peter M; Wright, Jim A; Kamwi, Richard; Uusiku, Petrina; Katokele, Stark; Snow, Robert W; Noor, Abdisalan M
2013-12-01
As malaria transmission declines, it becomes increasingly important to monitor changes in malaria incidence rather than prevalence. Here, a spatio-temporal model was used to identify constituencies with high malaria incidence to guide malaria control. Malaria cases were assembled across all age groups along with several environmental covariates. A Bayesian conditional-autoregressive model was used to model the spatial and temporal variation of incidence after adjusting for test positivity rates and health facility utilisation. Of the 144,744 malaria cases recorded in Namibia in 2009, 134,851 were suspected and 9893 were parasitologically confirmed. The mean annual incidence based on the Bayesian model predictions was 13 cases per 1000 population with the highest incidence predicted for constituencies bordering Angola and Zambia. The smoothed maps of incidence highlight trends in disease incidence. For Namibia, the 2009 maps provide a baseline for monitoring the targets of pre-elimination. PMID:24238079
NASA Astrophysics Data System (ADS)
Liu, Xinyu; Chen, Si; Luo, Yuemei; Bo, En; Wang, Nanshuo; Yu, Xiaojun; Liu, Linbo
2016-02-01
The evaluation of the endothelium coverage on the vessel wall is most wanted by cardiologists. Arterial endothelial cells play a crucial role in keeping low-density lipoprotein and leukocytes from entering into the intima. The damage of endothelial cells is considered as the first step of atherosclerosis development and the presence of endothelial cells is an indicator of arterial healing after stent implantation. Intravascular OCT (IVOCT) is the highest-resolution coronary imaging modality, but it is still limited by an axial resolution of 10-15 µm. This limitation in axial resolution hinders our ability to visualize cellular level details associated with coronary atherosclerosis. Spectral estimation optical coherence tomography (SE-OCT) uses modern spectral estimation techniques and may help reveal the microstructures underlying the resolution limit. In this presentation, we conduct an ex vivo study using SE-OCT to image the endothelium cells on the fresh swine aorta. We find that in OCT images with an axial resolution of 10 µm, we may gain the visibility of individual endothelium cells by applying the autoregressive spectral estimation techniques to enhance the axial resolution. We believe the SE-OCT can provide a potential to evaluate the coverage of endothelium cells using current IVOCT with a 10-µm axial resolution.
Digital spectral estimation and modeling of Space Shuttle flight data
NASA Technical Reports Server (NTRS)
Spanos, P. D.; Mushung, L. J.; Nelson, D. A. R., Jr.; Hamilton, D. A.
1988-01-01
Application of the digital signal processing technique of autoregressive-moving-average (ARMA) modeling to the estimation of power spectra and shock spectra from Space Shuttle lift-off flight accelerograms is described in this paper. The background for application to ARMA of lift-off accelerograms which are non-stationary in nature is exemplified through a step-by-step discussion of actual numerical results. Included is a discussion of pertinent mathematical background for the ARMA approximations. Potential areas for application of ARMA modeling in payload integration activities are suggested.
Jacob, Benjamin G; Griffith, Daniel A; Muturi, Ephantus J; Caamano, Erick X; Githure, John I; Novak, Robert J
2009-01-01
Background Autoregressive regression coefficients for Anopheles arabiensis aquatic habitat models are usually assessed using global error techniques and are reported as error covariance matrices. A global statistic, however, will summarize error estimates from multiple habitat locations. This makes it difficult to identify where there are clusters of An. arabiensis aquatic habitats of acceptable prediction. It is therefore useful to conduct some form of spatial error analysis to detect clusters of An. arabiensis aquatic habitats based on uncertainty residuals from individual sampled habitats. In this research, a method of error estimation for spatial simulation models was demonstrated using autocorrelation indices and eigenfunction spatial filters to distinguish among the effects of parameter uncertainty on a stochastic simulation of ecological sampled Anopheles aquatic habitat covariates. A test for diagnostic checking error residuals in an An. arabiensis aquatic habitat model may enable intervention efforts targeting productive habitats clusters, based on larval/pupal productivity, by using the asymptotic distribution of parameter estimates from a residual autocovariance matrix. The models considered in this research extends a normal regression analysis previously considered in the literature. Methods Field and remote-sampled data were collected during July 2006 to December 2007 in Karima rice-village complex in Mwea, Kenya. SAS 9.1.4® was used to explore univariate statistics, correlations, distributions, and to generate global autocorrelation statistics from the ecological sampled datasets. A local autocorrelation index was also generated using spatial covariance parameters (i.e., Moran's Indices) in a SAS/GIS® database. The Moran's statistic was decomposed into orthogonal and uncorrelated synthetic map pattern components using a Poisson model with a gamma-distributed mean (i.e. negative binomial regression). The eigenfunction values from the spatial
Spectral abundance fraction estimation of materials using Kalman filters
NASA Astrophysics Data System (ADS)
Wang, Su; Chang, Chein; Jensen, Janet L.; Jensen, James O.
2004-12-01
Kalman filter has been widely used in statistical signal processing for parameter estimation. Although a Kalman filter approach has been recently developed for spectral unmixing, referred to as Kalman filter-based linear unmixing (KFLU), its applicability to spectral characterization within a single pixel vector has not been explored. This paper presents a new application of Kalman filtering in spectral estimation and quantification. It develops a Kalman filter-based spectral signature esimator (KFSSE) which is different from the KFLU in the sense that the former performs a Kalman filter wavelength by wavelength across a spectral signature as opposed to the latter which implements a Kalman filter pixel vector by pixel vector in an image cube. The idea of the KFSSE is to implement the state equation to characterize the true spectral signature, while the measurement equation is being used to describe the spectral signature to be processed. Additionally, since a Kalman filter can accurately estimate spectral abundance fraction of a signature, our proposed KFSSE can further used for spectral quantification for subpixel targets and mixed pixel vectors, called Kalman filter-based spectral quantifier (KFSQ). Such spectral quantification is particularly important for chemical/biological defense which requires quantification of detected agents for damage control assessment. Several different types of hyperspectral data are used for experiments to demonstrate the ability of the KFSSE in estimation of spectral signature and the utility of the KFSQ in spectral quantification.
NASA Astrophysics Data System (ADS)
Liu, Xinyu; Yu, Xiaojun; Wang, Nanshuo; Bo, En; Luo, Yuemei; Chen, Si; Cui, Dongyao; Liu, Linbo
2016-03-01
The sample depth reflectivity profile of Fourier domain optical coherence tomography (FD-OCT) is estimated from the inverse Fourier transform of the spectral interference signals (interferograms). As a result, the axial resolution is fundamentally limited by the coherence length of the light source. We demonstrate an axial resolution improvement method by using the autoregressive spectral estimation technique to instead of the inverse Fourier transform to analyze the spectral interferograms, which is named as spectral estimation OCT (SE-OCT). SE-OCT improves the axial resolution by a factor of up to 4.7 compared with the corresponding FD-OCT. Furthermore, SE-OCT provides a complete sidelobe suppression in the point-spread function. Using phantoms such as an air wedge and micro particles, we prove the ability of resolution improvement. To test SE-OCT for real biological tissue, we image the rat cornea and demonstrate that SE-OCT enables clear identification of corneal endothelium anatomical details ex vivo. We also find that the performance of SE-OCT is depended on SNR of the feature object. To evaluate the potential usage and define the application scope of SE-OCT, we further investigate the property of SNR dependence and the artifacts that may be caused. We find SE-OCT may be uniquely suited for viewing high SNR layer structures, such as the epithelium and endothelium in cornea, retina and aorta. Given that SE-OCT can be implemented in the FD-OCT devices easily, the new capabilities provided by SE-OCT are likely to offer immediate improvements to the diagnosis and management of diseases based on OCT imaging.
Reference data set for camera spectral sensitivity estimation.
Darrodi, Maryam Mohammadzadeh; Finlayson, Graham; Goodman, Teresa; Mackiewicz, Michal
2015-03-01
In this article, we describe a spectral sensitivity measurement procedure at the National Physical Laboratory, London, with the aim of obtaining ground truth spectral sensitivity functions for Nikon D5100 and Sigma SD1 Merill cameras. The novelty of our data is that the potential measurement errors are estimated at each wavelength. We determine how well the measured spectral sensitivity functions represent the actual camera sensitivity functions (as a function of wavelength). The second contribution of this paper is to test the performance of various leading sensor estimation techniques implemented from the literature using measured and synthetic data and also evaluate them based on ground truth data for the two cameras. We conclude that the estimation techniques tested are not sufficiently accurate when compared with our measured ground truth data and that there remains significant scope to improve estimation algorithms for spectral estimation. To help in this endeavor, we will make all our data available online for the community. PMID:26366649
Yield estimation from hyperspectral imagery using Spectral Angle Mapper (SAM)
Technology Transfer Automated Retrieval System (TEKTRAN)
Vegetation indices (VIs) derived from remotely sensed imagery are commonly used to estimate crop yields. Spectral angle mapper (SAM) provides an alternative approach to quantifying the spectral differences among all pixels in imagery and therefore has the potential for mapping yield variability. The...
Estimation of agronomic variables using spectral signatures
NASA Technical Reports Server (NTRS)
Goel, N. S.; Thompson, R. L.
1984-01-01
Techniques for the determination of leaf area index or leaf angle distribution from remote-sensing canopy-reflectance (CR) measurements are developed on the basis of empirical models relating CR to parameters such as soil and vegetation spectral properties, solar flux, and viewing angle. A general procedure for inverting CR models is presented and applied to the models of Suits (1972), Verhoef and Bunnik (1981), and Norman (1979) in the IR range. Numerical results for a soybean canopy are compared in a table, and the error sensitivity of the inverted models is shown to be relatively high, requiring the use of ancillary data such as soil reflectance, leaf reflectance, and leaf transmittance.
Schaffer, Thorsten; Hensel, Bernhard; Weigand, Christian; Schüttler, Jürgen; Jeleazcov, Christian
2014-10-01
Heart rate variability (HRV) analysis is increasingly used in anaesthesia and intensive care monitoring of spontaneous breathing and mechanical ventilated patients. In the frequency domain, different estimation methods of the power spectral density (PSD) of RR-intervals lead to different results. Therefore, we investigated the PSD estimates of fast Fourier transform (FFT), autoregressive modeling (AR) and Lomb-Scargle periodogram (LSP) for 25 young healthy subjects subjected to metronomic breathing. The optimum method for determination of HRV spectral parameters under paced respiration was identified by evaluating the relative error (RE) and the root mean square relative error (RMSRE) for each breathing frequency (BF) and spectral estimation method. Additionally, the sympathovagal balance was investigated by calculating the low frequency/high frequency (LF/HF) ratio. Above 7 breaths per minute, all methods showed a significant increase in LF/HF ratio with increasing BF. On average, the RMSRE of FFT was lower than for LSP and AR. Therefore, under paced respiration conditions, estimating RR-interval PSD using FFT is recommend. PMID:23508826
Spectral reflectance estimation using a six-color scanner
NASA Astrophysics Data System (ADS)
Tominaga, Shoji; Kohno, Satoshi; Kakinuma, Hirokazu; Nohara, Fuminori; Horiuchi, Takahiko
2009-01-01
A method is proposed for estimating the spectral reflectance function of an object surface by using a six-color scanner. The scanner is regarded as a six-band spectral imaging system, since it captures six color channels in total from two separate scans using two difference lamps. First, we describe the basic characteristics of the imaging systems for a HP color scanner and a multiband camera used for comparison. Second, we describe a computational method for recovering surface-spectral reflectances from the noisy sensor outputs. A LMMSE estimator is presented as an optimal estimator. We discuss the reflectance estimation for non-flat surfaces with shading effect. A solution method is presented for the reliable reflectance estimation. Finally, the performance of the proposed method is examined in detail on experiments using the Macbeth Color Checker and non-flat objects.
Spectral Irradiance Estimation in the Near Infrared
NASA Astrophysics Data System (ADS)
Woolley, R. C.
2003-12-01
Currently, photometric data in the near IR is publicly available in only a few select wavelength bands from sources such as IRAS, 2Mass, and the MSX point source catalog. Efforts have been made to construct composite spectra using existing photometry by Martin Cohen (1995) and Michael P. Egan (1996). Astrometric and photometric databases in various wavelength bands are needed as a tool for improving pointing accuracy and calibration of space-based infrared observations. Methods for extrapolating and interpolating existing photometry in order to produce these databases are not as straightforward as previously expected. Spectral irradiance in the IR region deviates greatly from the standard Planck function and is dominated by absorption and emission in the stellar atmosphere. Since most stars that are bright in the IR are late K and M class stars (often variable), this problem is compounded. In this poster a rough method is put forth to extrapolate irradiances from existing 2Mass magnitudes and B-V temperatures. A more rigorous analysis as well as observational data are needed to verify the results.
Comparison of spectral estimators for characterizing fractionated atrial electrograms
2013-01-01
Background Complex fractionated atrial electrograms (CFAE) acquired during atrial fibrillation (AF) are commonly assessed using the discrete Fourier transform (DFT), but this can lead to inaccuracy. In this study, spectral estimators derived by averaging the autocorrelation function at lags were compared to the DFT. Method Bipolar CFAE of at least 16 s duration were obtained from pulmonary vein ostia and left atrial free wall sites (9 paroxysmal and 10 persistent AF patients). Power spectra were computed using the DFT and three other methods: 1. a novel spectral estimator based on signal averaging (NSE), 2. the NSE with harmonic removal (NSH), and 3. the autocorrelation function average at lags (AFA). Three spectral parameters were calculated: 1. the largest fundamental spectral peak, known as the dominant frequency (DF), 2. the DF amplitude (DA), and 3. the mean spectral profile (MP), which quantifies noise floor level. For each spectral estimator and parameter, the significance of the difference between paroxysmal and persistent AF was determined. Results For all estimators, mean DA and mean DF values were higher in persistent AF, while the mean MP value was higher in paroxysmal AF. The differences in means between paroxysmals and persistents were highly significant for 3/3 NSE and NSH measurements and for 2/3 DFT and AFA measurements (p<0.001). For all estimators, the standard deviation in DA and MP values were higher in persistent AF, while the standard deviation in DF value was higher in paroxysmal AF. Differences in standard deviations between paroxysmals and persistents were highly significant in 2/3 NSE and NSH measurements, in 1/3 AFA measurements, and in 0/3 DFT measurements. Conclusions Measurements made from all four spectral estimators were in agreement as to whether the means and standard deviations in three spectral parameters were greater in CFAEs acquired from paroxysmal or in persistent AF patients. Since the measurements were consistent, use of
Spectral estimates of solar radiation intercepted by corn canopies
NASA Technical Reports Server (NTRS)
Bauer, M. E. (Principal Investigator); Daughtry, C. S. T.; Gallo, K. P.
1982-01-01
Reflectance factor data were acquired with a Landsat band radiometer throughout two growing seasons for corn (Zea mays L.) canopies differing in planting dates, populations, and soil types. Agronomic data collected included leaf area index (LAI), biomass, development stage, and final grain yields. The spectral variable, greenness, was associated with 78 percent of the variation in LAI over all treatments. Single observations of LAI or greenness have limited value in predicting corn yields. The proportions of solar radiation intercepted (SRI) by these canopies were estimated using either measured LAI or greenness. Both SRI estimates, when accumulated over the growing season, accounted for approximately 65 percent of the variation in yields. Models which simulated the daily effects of weather and intercepted solar radiation on growth had the highest correlations to grain yields. This concept of estimating intercepted solar radiation using spectral data represents a viable approach for merging spectral and meteorological data for crop yield models.
Estimation of Wheat Agronomic Parameters using New Spectral Indices
Jin, Xiu-liang; Diao, Wan-ying; Xiao, Chun-hua; Wang, Fang-yong; Chen, Bing; Wang, Ke-ru; Li, Shao-kun
2013-01-01
Crop agronomic parameters (leaf area index (LAI), nitrogen (N) uptake, total chlorophyll (Chl) content ) are very important for the prediction of crop growth. The objective of this experiment was to investigate whether the wheat LAI, N uptake, and total Chl content could be accurately predicted using spectral indices collected at different stages of wheat growth. Firstly, the product of the optimized soil-adjusted vegetation index and wheat biomass dry weight (OSAVI×BDW) were used to estimate LAI, N uptake, and total Chl content; secondly, BDW was replaced by spectral indices to establish new spectral indices (OSAVI×OSAVI, OSAVI×SIPI, OSAVI×CIred edge, OSAVI×CIgreen mode and OSAVI×EVI2); finally, we used the new spectral indices for estimating LAI, N uptake, and total Chl content. The results showed that the new spectral indices could be used to accurately estimate LAI, N uptake, and total Chl content. The highest R2 and the lowest RMSEs were 0.711 and 0.78 (OSAVI×EVI2), 0.785 and 3.98 g/m2 (OSAVI×CIred edge) and 0.846 and 0.65 g/m2 (OSAVI×CIred edge) for LAI, nitrogen uptake and total Chl content, respectively. The new spectral indices performed better than the OSAVI alone, and the problems of a lack of sensitivity at earlier growth stages and saturation at later growth stages, which are typically associated with the OSAVI, were improved. The overall results indicated that this new spectral indices provided the best approximation for the estimation of agronomic indices for all growth stages of wheat. PMID:24023639
Optimal estimation of spectral reflectance based on metamerism
NASA Astrophysics Data System (ADS)
Chou, Tzren-Ru; Lin, Wei-Ju
2012-01-01
In this paper, we proposed an accurate estimation method for spectral reflectance of objects captured in an image. The spectral reflectance is simply modeled by a linear combination of three basic spectrums of R, G, and B colors respectively, named as spectral reflective bases of objects, which are acquired by solving a linear system based on the principle of color metamerism. Some experiments were performed to evaluate the accuracy of the estimated spectral reflectance of objects. The average mean square error of 24 colors in Macbeth checker between we simulated and the measured is 0.0866, and the maximum is 0.310. In addition, the average color difference of the 24 colors is less than 1.5 under the D65 illuminant. There are 13 colors having their color difference values less than 1, and other 8 colors having the values during the range of 1 and 2. Only three colors are relatively larger, with the differences of 2.558, 4.130 and 2.569, from the colors of No. 2, No. 13, and No. 18 in Macbeth checker respectively. Furthermore, the computational cost of this spectral estimation is very low and suitable for many practical applications in real time.
[Optimized Spectral Indices Based Estimation of Forage Grass Biomass].
An, Hai-bo; Li, Fei; Zhao, Meng-li; Liu, Ya-jun
2015-11-01
As an important indicator of forage production, aboveground biomass will directly illustrate the growth of forage grass. Therefore, Real-time monitoring biomass of forage grass play a crucial role in performing suitable grazing and management in artificial and natural grassland. However, traditional sampling and measuring are time-consuming and labor-intensive. Recently, development of hyperspectral remote sensing provides the feasibility in timely and nondestructive deriving biomass of forage grass. In the present study, the main objectives were to explore the robustness of published and optimized spectral indices in estimating biomass of forage grass in natural and artificial pasture. The natural pasture with four grazing density (control, light grazing, moderate grazing and high grazing) was designed in desert steppe, and different forage cultivars with different N rate were conducted in artificial forage fields in Inner Mongolia. The canopy reflectance and biomass in each plot were measured during critical stages. The result showed that, due to the influence in canopy structure and biomass, the canopy reflectance have a great difference in different type of forage grass. The best performing spectral index varied in different species of forage grass with different treatments (R² = 0.00-0.69). The predictive ability of spectral indices decreased under low biomass of desert steppe, while red band based spectral indices lost sensitivity under moderate-high biomass of forage maize. When band combinations of simple ratio and normalized difference spectral indices were optimized in combined datasets of natural and artificial grassland, optimized spectral indices significant increased predictive ability and the model between biomass and optimized spectral indices had the highest R² (R² = 0.72) compared to published spectral indices. Sensitive analysis further confirmed that the optimized index had the lowest noise equivalent and were the best performing index in
Spectral estimation on a sphere in geophysics and cosmology
NASA Astrophysics Data System (ADS)
Dahlen, F. A.; Simons, Frederik J.
2008-09-01
We address the problem of estimating the spherical-harmonic power spectrum of a statistically isotropic scalar signal from noise-contaminated data on a region of the unit sphere. Three different methods of spectral estimation are considered: (i) the spherical analogue of the one-dimensional (1-D) periodogram, (ii) the maximum-likelihood method and (iii) a spherical analogue of the 1-D multitaper method. The periodogram exhibits strong spectral leakage, especially for small regions of area A << 4π, and is generally unsuitable for spherical spectral analysis applications, just as it is in 1-D. The maximum-likelihood method is particularly useful in the case of nearly-whole-sphere coverage, A ~ 4π, and has been widely used in cosmology to estimate the spectrum of the cosmic microwave background radiation from spacecraft observations. The spherical multitaper method affords easy control over the fundamental trade-off between spectral resolution and variance, and is easily implemented regardless of the region size, requiring neither non-linear iteration nor large-scale matrix inversion. As a result, the method is ideally suited for most applications in geophysics, geodesy or planetary science, where the objective is to obtain a spatially localized estimate of the spectrum of a signal from noisy data within a pre-selected and typically small region.
Estimating neugebauer primaries for multi-channel spectral printing modeling
NASA Astrophysics Data System (ADS)
Slavuj, Radovan; Coppel, Ludovic G.; Olen, Melissa; Hardeberg, Jon Yngve
2014-02-01
Multichannel printer modeling has been an active area of research in the field of spectral printing. The most commonly used models for characterization of such systems are the spectral Neugebauer (SN) and its extensions. This work addresses issues that can arise during calibration and testing of the SN model when modelling a 7-colorant printer. Since most substrates are limited in their capacity to take in large amount of ink, it is not always possible to print all colorant combinations necessary to determine the Neugebauer primaries (NP). A common solution is to estimate the nonprintable Neugebauer primaries from the single colorant primaries using the Kubelka-Munk (KM) optical model. In this work we test whether a better estimate can be obtained using general radiative transfer theory, which better represents the angular variation of the reflectance from highly absorbing media, and takes surface scattering into account. For this purpose we use the DORT2002 model. We conclude DORT2002 does not offer significant improvements over KM in the estimation of the NPs, but a significant improvement is obtained when using a simple surface scattering model. When the estimated primaries are used as inputs to the SN model instead of measured ones, it is found the SN model performs the same or better in terms of color difference and spectral error. If the mixed measured and estimated primaries are used as inputs to the SN model, it performs better than using either measured or estimated.
Induction machine faults detection using stator current parametric spectral estimation
NASA Astrophysics Data System (ADS)
El Bouchikhi, El Houssin; Choqueuse, Vincent; Benbouzid, Mohamed
2015-02-01
Current spectrum analysis is a proven technique for fault diagnosis in electrical machines. Current spectral estimation is usually performed using classical techniques such as periodogram (FFT) or its extensions. However, these techniques have several drawbacks since their frequency resolution is limited and additional post-processing algorithms are required to extract a relevant fault detection criterion. Therefore, this paper proposes a new parametric spectral estimator that fully exploits the faults sensitive frequencies. The proposed technique is based on the maximum likelihood estimator (MLE) and offers high-resolution capabilities. Based on this approach, a fault criterion is derived for detecting several fault types. The proposed technique is assessed using simulation signals, issued from a coupled electromagnetic circuits approach-based simulation tool for mechanical unbalance and electrical asymmetry faults detection. It is afterward validated using experiments on a 0.75-kW induction machine test bed for the particular case of bearing faults.
NASA Astrophysics Data System (ADS)
Baghi, Quentin; Métris, Gilles; Bergé, Joël; Christophe, Bruno; Touboul, Pierre; Rodrigues, Manuel
2016-06-01
We present a Gaussian regression method for time series with missing data and stationary residuals of unknown power spectral density (PSD). The missing data are efficiently estimated by their conditional expectation as in universal Kriging based on the circulant approximation of the complete data covariance. After initialization with an autoregressive fit of the noise, a few iterations of estimation/reconstruction steps are performed until convergence of the regression and PSD estimates, in a way similar to the expectation-conditional-maximization algorithm. The estimation can be performed for an arbitrary PSD provided that it is sufficiently smooth. The algorithm is developed in the framework of the MICROSCOPE space mission whose goal is to test the weak equivalence principle (WEP) with a precision of 10-15. We show by numerical simulations that the developed method allows us to meet three major requirements: to maintain the targeted precision of the WEP test in spite of the loss of data, to calculate a reliable estimate of this precision and of the noise level, and finally to provide consistent and faithful reconstructed data to the scientific community.
Yield estimation of sugarcane based on agrometeorological-spectral models
NASA Technical Reports Server (NTRS)
Rudorff, Bernardo Friedrich Theodor; Batista, Getulio Teixeira
1990-01-01
This work has the objective to assess the performance of a yield estimation model for sugarcane (Succharum officinarum). The model uses orbital gathered spectral data along with yield estimated from an agrometeorological model. The test site includes the sugarcane plantations of the Barra Grande Plant located in Lencois Paulista municipality in Sao Paulo State. Production data of four crop years were analyzed. Yield data observed in the first crop year (1983/84) were regressed against spectral and agrometeorological data of that same year. This provided the model to predict the yield for the following crop year i.e., 1984/85. The model to predict the yield of subsequent years (up to 1987/88) were developed similarly, incorporating all previous years data. The yield estimations obtained from these models explained 69, 54, and 50 percent of the yield variation in the 1984/85, 1985/86, and 1986/87 crop years, respectively. The accuracy of yield estimations based on spectral data only (vegetation index model) and on agrometeorological data only (agrometeorological model) were also investigated.
Respiratory impedance spectral estimation for digitally created random noise.
Davis, K A; Lutchen, K R
1991-01-01
Measurement of respiratory input mechanical impedance (Zrs) is noninvasive, requires minimal subject cooperation, and contains information related to mechanical lung function. A common approach to measure Zrs is to apply random noise pressure signals at the airway opening, measure the resulting flow variations, and then estimate Zrs using Fast-Fourier Transform (FFT) techniques. The goal of this study was to quantify how several signal processing issues affect the quality of a Zrs spectral estimate when the input pressure sequence is created digitally. Random noise driven pressure and flow time domain data were simulated for three models, which permitted predictions of Zrs characteristics previously reported from 0-4, 4-32, and 4-200 Hz. Then, the quality of the Zrs estimate was evaluated as a function of the number of runs ensemble averaged, windowing, flow signal-to-noise ratio (SNR), and pressure spectral magnitude shape magnitude of P(j omega). For a magnitude of P(j omega) with uniform power distribution and a SNR less than 100, the 0-4 Hz and 4-200 Hz Zrs estimates for 10 runs were poor (minimum coherence gamma 2 less than 0.75) particularly where Zrs is high. When the SNR greater than 200 and 10 runs were averaged, the minimum gamma 2 greater than 0.95. However, when magnitude of P(j omega) was matched to magnitude of Zrs, gamma 2 greater than 0.91 even for 5 runs and a SNR of 20. For data created digitally with equally spaced spectral content, the rectangular window was superior to the Hanning. Finally, coherence alone may not be a reliable measure of Zrs quality because coherence is only an estimate itself. We conclude that an accurate estimate of Zrs is best obtained by matching magnitude of P(j omega) to magnitude of Zin (subject and speaker) and using rectangular windowing. PMID:2048776
Estimation of spectral emissivity in the thermal infrared
NASA Technical Reports Server (NTRS)
Kryskowski, David; Maxwell, J. R.
1993-01-01
A number of algorithms are available in the literature that attempt to remove most of the effects of temperature from thermal multispectral data where the final goal is to extract emissivity differences. Early approaches include adjacent spectral band ratioing, broad band radiance normalization and the use of one band where emissivities are generally high (e.g., 11 to 12 micrometers) to determine the temperature. More recent work has produced two techniques that use data averaging to extract temperature to leave a quantity related to emissivity changes. These two techniques have been investigated and compared and appear to provide reasonable results. The analysis presented in this paper develops a thermal IR multispectral temperature/emissivity estimation procedure based on formal estimation theory, Gaussian statistics, and a stochastic radiance signal model including the effects of both temperature and emissivity. The importance of this work is that this is an optimal estimation procedure which will provide minimum variance estimates of temperature and emissivity changes directly. Section 2 discusses optimal linear spectral emissivity estimation and Section 3 is a summary.
Spectral estimators of absorbed photosynthetically active radiation in corn canopies
NASA Technical Reports Server (NTRS)
Gallo, K. P.; Daughtry, C. S. T.; Bauer, M. E.
1985-01-01
Most models of crop growth and yield require an estimate of canopy leaf area index (LAI) or absorption of radiation. Relationships between photosynthetically active radiation (PAR) absorbed by corn canopies and the spectral reflectance of the canopies were investigated. Reflectance factor data were acquired with a Landsat MSS band radiometer. From planting to silking, the three spectrally predicted vegetation indices examined were associated with more than 95 percent of the variability in absorbed PAR. The relationships developed between absorbed PAR and the three indices were evaluated with reflectance factor data acquired from corn canopies planted in 1979 through 1982. Seasonal cumulations of measured LAI and each of the three indices were associated with greater than 50 percent of the variation in final grain yields from the test years. Seasonal cumulations of daily absorbed PAR were associated with up to 73 percent of the variation in final grain yields. Absorbed PAR, cumulated through the growing season, is a better indicator of yield than cumulated leaf area index. Absorbed PAR may be estimated reliably from spectral reflectance data of crop canopies.
Spectral estimators of absorbed photosynthetically active radiation in corn canopies
NASA Technical Reports Server (NTRS)
Gallo, K. P.; Daughtry, C. S. T.; Bauer, M. E.
1984-01-01
Most models of crop growth and yield require an estimate of canopy leaf area index (LAI) or absorption of radiation. Relationships between photosynthetically active radiation (PAR) absorbed by corn canopies and the spectral reflectance of the canopies were investigated. Reflectance factor data were acquired with a LANDSAT MSS band radiometer. From planting to silking, the three spectrally predicted vegetation indices examined were associated with more than 95% of the variability in absorbed PAR. The relationships developed between absorbed PAR and the three indices were evaluated with reflectance factor data acquired from corn canopies planted in 1979 through 1982. Seasonal cumulations of measured LAI and each of the three indices were associated with greater than 50% of the variation in final grain yields from the test years. Seasonal cumulations of daily absorbed PAR were associated with up to 73% of the variation in final grain yields. Absorbed PAR, cumulated through the growing season, is a better indicator of yield than cumulated leaf area index. Absorbed PAR may be estimated reliably from spectral reflectance data of crop canopies.
Seeing red in cyclic stratigraphy: Spectral noise estimation for astrochronology
NASA Astrophysics Data System (ADS)
Meyers, Stephen R.
2012-09-01
Fundamental to the development of astronomical time scales is the recognition of oscillatory variability within stratigraphic data and its evaluation relative to a null "noise" hypothesis. In this study, Monte Carlo simulations are used to investigate the suitability of two commonly used noise hypotheses (the "conventional" and "robust" AR1 approaches), and the results highlight important limitations in both for cyclostratigraphic application. Perhaps most problematic, the robust AR1 method can result in inflated confidence level estimates and excessive clumping of false positives within the low frequency portion of the spectrum, especially when the underlying noise process has a high lag-1 autocorrelation. Given typical cyclostratigraphic records, this technique will often impose "significant" eccentricity band variability, even in the case of pure AR1 noise. An alternative spectral noise estimation method is proposed to overcome these problems, which simultaneously allows for departures from the AR1 assumption, and obtains high statistical power—that is, the ability to accurately identify astronomical signals when they are present in the data. We apply the method to un-tunedδ18O data from Miocene sediments of the Ceara Rise, indicating statistically significant spectral power at frequencies that are consistent with the published orbital interpretation of Weedon et al. (1997). Furthermore, evaluation of the frequency arrangement of the significant spatial bedding periods, using the average spectral misfit method for astrochronologic testing, reveals that the null hypothesis of no orbital influence can be rejected with a high degree of confidence (the 99.8% confidence level).
Convergence of oscillator spectral estimators for counted-frequency measurements.
NASA Technical Reports Server (NTRS)
Tausworthe, R. C.
1972-01-01
A common intermediary connecting frequency-noise calibration or testing of an oscillator to useful applications is the spectral density of the frequency-deviating process. In attempting to turn test data into predicts of performance characteristics, one is naturally led to estimation of statistical values by sample-mean and sample-variance techniques. However, sample means and sample variances themselves are statistical quantities that do not necessarily converge (in the mean-square sense) to actual ensemble-average means and variances, except perhaps for excessively large sample sizes. This is especially true for the flicker noise component of oscillators. This article shows, for the various types of noises found in oscillators, how sample averages converge (or do not converge) to their statistical counterparts. The convergence rate is shown to be the same for all oscillators of a given spectral type.
Efficient spectral estimation for time series with intermittent gaps
NASA Astrophysics Data System (ADS)
Smith, L. T.; Constable, C.
2009-12-01
Data from magnetic satellites like CHAMP, Ørsted, and Swarm can be used to study electromagnetic induction in Earth’s mantle. Time series of internal and external spherical harmonic coefficients (usually those associated with the predominantly dipolar structure of ring current variations) are used to determine Earth’s electromagnetic response as a function of frequency of the external variations. Inversion of this response can yield information about electrical conductivity variations in Earth’s mantle. The inductive response depends on frequency through skin depth, so it is desirable to work with the longest time series possible. Intermittent gaps in available data complicate attempts to estimate the power or cross spectra and thus the electromagnetic response for satellite records. Complete data series are most effectively analyzed using direct multi-taper spectral estimation, either with prolate multitapers that efficiently minimize broadband bias, or with a set designed to minimize local bias. The latter group have frequently been approximated by sine tapers. Intermittent gaps in data may be patched over using custom designed interpolation. We focus on a different approach, using sets of multitapers explicitly designed to accommodate gaps in the data. The optimization problems for the prolate and minimum bias tapers are altered to allow a specific arrangement of data samples, producing a modified eigenvalue-eigenfunction problem. We have shown that the prolate tapers with gaps and the minimum bias tapers with gaps provide higher resolution spectral estimates with less leakage than spectral averaging of data sections bounded by gaps. Current work is focused on producing efficient algorithms for spectral estimation of data series with gaps. A major limitation is the time and memory needed for the solution of large eigenvalue problems used to calculate the tapers for long time series. Fortunately only a limited set of the largest eigenvalues are needed, and
Alternative techniques for high-resolution spectral estimation of spectrally encoded endoscopy
NASA Astrophysics Data System (ADS)
Mousavi, Mahta; Duan, Lian; Javidi, Tara; Ellerbee, Audrey K.
2015-09-01
Spectrally encoded endoscopy (SEE) is a minimally invasive optical imaging modality capable of fast confocal imaging of internal tissue structures. Modern SEE systems use coherent sources to image deep within the tissue and data are processed similar to optical coherence tomography (OCT); however, standard processing of SEE data via the Fast Fourier Transform (FFT) leads to degradation of the axial resolution as the bandwidth of the source shrinks, resulting in a well-known trade-off between speed and axial resolution. Recognizing the limitation of FFT as a general spectral estimation algorithm to only take into account samples collected by the detector, in this work we investigate alternative high-resolution spectral estimation algorithms that exploit information such as sparsity and the general region position of the bulk sample to improve the axial resolution of processed SEE data. We validate the performance of these algorithms using bothMATLAB simulations and analysis of experimental results generated from a home-built OCT system to simulate an SEE system with variable scan rates. Our results open a new door towards using non-FFT algorithms to generate higher quality (i.e., higher resolution) SEE images at correspondingly fast scan rates, resulting in systems that are more accurate and more comfortable for patients due to the reduced image time.
Spectral estimation of temporal series at unequal intervals.
Abraira, V; Ibarz, J M
1986-06-01
Many biological variables present rhythmic oscillations at different frequencies. Most common techniques, which statistically characterize temporal series and permit the study of these rhythms, require equidistant sampling. However, it is not always possible to register at regular intervals many of the variables under study, either because of the nature of the phenomenon which generates them or because of the difficulty in obtaining the samples. This paper proposes a method for spectral estimation by means of fitting to the cosine functions of sampled variables using a nonuniform point process. PMID:3754801
How Accurately Do Spectral Methods Estimate Effective Elastic Thickness?
NASA Astrophysics Data System (ADS)
Perez-Gussinye, M.; Lowry, A. R.; Watts, A. B.; Velicogna, I.
2002-12-01
The effective elastic thickness, Te, is an important parameter that has the potential to provide information on the long-term thermal and mechanical properties of the the lithosphere. Previous studies have estimated Te using both forward and inverse (spectral) methods. While there is generally good agreement between the results obtained using these methods, spectral methods are limited because they depend on the spectral estimator and the window size chosen for analysis. In order to address this problem, we have used a multitaper technique which yields optimal estimates of the bias and variance of the Bouguer coherence function relating topography and gravity anomaly data. The technique has been tested using realistic synthetic topography and gravity. Synthetic data were generated assuming surface and sub-surface (buried) loading of an elastic plate with fractal statistics consistent with real data sets. The cases of uniform and spatially varying Te are examined. The topography and gravity anomaly data consist of 2000x2000 km grids sampled at 8 km interval. The bias in the Te estimate is assessed from the difference between the true Te value and the mean from analyzing 100 overlapping windows within the 2000x2000 km data grids. For the case in which Te is uniform, the bias and variance decrease with window size and increase with increasing true Te value. In the case of a spatially varying Te, however, there is a trade-off between spatial resolution and variance. With increasing window size the variance of the Te estimate decreases, but the spatial changes in Te are smeared out. We find that for a Te distribution consisting of a strong central circular region of Te=50 km (radius 600 km) and progressively smaller Te towards its edges, the 800x800 and 1000x1000 km window gave the best compromise between spatial resolution and variance. Our studies demonstrate that assumed stationarity of the relationship between gravity and topography data yields good results even in
NASA Astrophysics Data System (ADS)
Lopes, Sílvia R. C.; Prass, Taiane S.
2014-05-01
Here we present a theoretical study on the main properties of Fractionally Integrated Exponential Generalized Autoregressive Conditional Heteroskedastic (FIEGARCH) processes. We analyze the conditions for the existence, the invertibility, the stationarity and the ergodicity of these processes. We prove that, if { is a FIEGARCH(p,d,q) process then, under mild conditions, { is an ARFIMA(q,d,0) with correlated innovations, that is, an autoregressive fractionally integrated moving average process. The convergence order for the polynomial coefficients that describes the volatility is presented and results related to the spectral representation and to the covariance structure of both processes { and { are discussed. Expressions for the kurtosis and the asymmetry measures for any stationary FIEGARCH(p,d,q) process are also derived. The h-step ahead forecast for the processes {, { and { are given with their respective mean square error of forecast. The work also presents a Monte Carlo simulation study showing how to generate, estimate and forecast based on six different FIEGARCH models. The forecasting performance of six models belonging to the class of autoregressive conditional heteroskedastic models (namely, ARCH-type models) and radial basis models is compared through an empirical application to Brazilian stock market exchange index.
Spectral estimates of net radiation and soil heat flux
Daughtry, C.S.T.; Kustas, W.P.; Moran, M.S.; Pinter, P. J., Jr.; Jackson, R. D.; Brown, P.W.; Nichols, W.D.; Gay, L.W.
1990-01-01
Conventional methods of measuring surface energy balance are point measurements and represent only a small area. Remote sensing offers a potential means of measuring outgoing fluxes over large areas at the spatial resolution of the sensor. The objective of this study was to estimate net radiation (Rn) and soil heat flux (G) using remotely sensed multispectral data acquired from an aircraft over large agricultural fields. Ground-based instruments measured Rn and G at nine locations along the flight lines. Incoming fluxes were also measured by ground-based instruments. Outgoing fluxes were estimated using remotely sensed data. Remote Rn, estimated as the algebraic sum of incoming and outgoing fluxes, slightly underestimated Rn measured by the ground-based net radiometers. The mean absolute errors for remote Rn minus measured Rn were less than 7%. Remote G, estimated as a function of a spectral vegetation index and remote Rn, slightly overestimated measured G; however, the mean absolute error for remote G was 13%. Some of the differences between measured and remote values of Rn and G are associated with differences in instrument designs and measurement techniques. The root mean square error for available energy (Rn - G) was 12%. Thus, methods using both ground-based and remotely sensed data can provide reliable estimates of the available energy which can be partitioned into sensible and latent heat under nonadvective conditions. ?? 1990.
Optimizing spectral wave estimates with adjoint-based sensitivity maps
NASA Astrophysics Data System (ADS)
Orzech, Mark; Veeramony, Jay; Flampouris, Stylianos
2014-04-01
A discrete numerical adjoint has recently been developed for the stochastic wave model SWAN. In the present study, this adjoint code is used to construct spectral sensitivity maps for two nearshore domains. The maps display the correlations of spectral energy levels throughout the domain with the observed energy levels at a selected location or region of interest (LOI/ROI), providing a full spectrum of values at all locations in the domain. We investigate the effectiveness of sensitivity maps based on significant wave height ( H s ) in determining alternate offshore instrument deployment sites when a chosen nearshore location or region is inaccessible. Wave and bathymetry datasets are employed from one shallower, small-scale domain (Duck, NC) and one deeper, larger-scale domain (San Diego, CA). The effects of seasonal changes in wave climate, errors in bathymetry, and multiple assimilation points on sensitivity map shapes and model performance are investigated. Model accuracy is evaluated by comparing spectral statistics as well as with an RMS skill score, which estimates a mean model-data error across all spectral bins. Results indicate that data assimilation from identified high-sensitivity alternate locations consistently improves model performance at nearshore LOIs, while assimilation from low-sensitivity locations results in lesser or no improvement. Use of sub-sampled or alongshore-averaged bathymetry has a domain-specific effect on model performance when assimilating from a high-sensitivity alternate location. When multiple alternate assimilation locations are used from areas of lower sensitivity, model performance may be worse than with a single, high-sensitivity assimilation point.
A comparison of spectral estimation methods for the analysis of sibilant fricatives
Reidy, Patrick F.
2015-01-01
It has been argued that, to ensure accurate spectral feature estimates for sibilants, the spectral estimation method should include a low-variance spectral estimator; however, no empirical evaluation of estimation methods in terms of feature estimates has been given. The spectra of /s/ and /ʃ/ were estimated with different methods that varied the pre-emphasis filter and estimator. These methods were evaluated in terms of effects on two features (centroid and degree of sibilance) and on the detection of four linguistic contrasts within these features. Estimation method affected the spectral features but none of the tested linguistic contrasts. PMID:25920873
Smallwood, D. O.
1996-01-01
It is shown that the usual method for estimating the coherence functions (ordinary, partial, and multiple) for a general multiple-input! multiple-output problem can be expressed as a modified form of Cholesky decomposition of the cross-spectral density matrix of the input and output records. The results can be equivalently obtained using singular value decomposition (SVD) of the cross-spectral density matrix. Using SVD suggests a new form of fractional coherence. The formulation as a SVD problem also suggests a way to order the inputs when a natural physical order of the inputs is absent.
Informed spectral analysis: audio signal parameter estimation using side information
NASA Astrophysics Data System (ADS)
Fourer, Dominique; Marchand, Sylvain
2013-12-01
Parametric models are of great interest for representing and manipulating sounds. However, the quality of the resulting signals depends on the precision of the parameters. When the signals are available, these parameters can be estimated, but the presence of noise decreases the resulting precision of the estimation. Furthermore, the Cramér-Rao bound shows the minimal error reachable with the best estimator, which can be insufficient for demanding applications. These limitations can be overcome by using the coding approach which consists in directly transmitting the parameters with the best precision using the minimal bitrate. However, this approach does not take advantage of the information provided by the estimation from the signal and may require a larger bitrate and a loss of compatibility with existing file formats. The purpose of this article is to propose a compromised approach, called the 'informed approach,' which combines analysis with (coded) side information in order to increase the precision of parameter estimation using a lower bitrate than pure coding approaches, the audio signal being known. Thus, the analysis problem is presented in a coder/decoder configuration where the side information is computed and inaudibly embedded into the mixture signal at the coder. At the decoder, the extra information is extracted and is used to assist the analysis process. This study proposes applying this approach to audio spectral analysis using sinusoidal modeling which is a well-known model with practical applications and where theoretical bounds have been calculated. This work aims at uncovering new approaches for audio quality-based applications. It provides a solution for challenging problems like active listening of music, source separation, and realistic sound transformations.
Bayesian parameter estimation in spectral quantitative photoacoustic tomography
NASA Astrophysics Data System (ADS)
Pulkkinen, Aki; Cox, Ben T.; Arridge, Simon R.; Kaipio, Jari P.; Tarvainen, Tanja
2016-03-01
Photoacoustic tomography (PAT) is an imaging technique combining strong contrast of optical imaging to high spatial resolution of ultrasound imaging. These strengths are achieved via photoacoustic effect, where a spatial absorption of light pulse is converted into a measurable propagating ultrasound wave. The method is seen as a potential tool for small animal imaging, pre-clinical investigations, study of blood vessels and vasculature, as well as for cancer imaging. The goal in PAT is to form an image of the absorbed optical energy density field via acoustic inverse problem approaches from the measured ultrasound data. Quantitative PAT (QPAT) proceeds from these images and forms quantitative estimates of the optical properties of the target. This optical inverse problem of QPAT is illposed. To alleviate the issue, spectral QPAT (SQPAT) utilizes PAT data formed at multiple optical wavelengths simultaneously with optical parameter models of tissue to form quantitative estimates of the parameters of interest. In this work, the inverse problem of SQPAT is investigated. Light propagation is modelled using the diffusion equation. Optical absorption is described with chromophore concentration weighted sum of known chromophore absorption spectra. Scattering is described by Mie scattering theory with an exponential power law. In the inverse problem, the spatially varying unknown parameters of interest are the chromophore concentrations, the Mie scattering parameters (power law factor and the exponent), and Gruneisen parameter. The inverse problem is approached with a Bayesian method. It is numerically demonstrated, that estimation of all parameters of interest is possible with the approach.
Tracking closely spaced multiple sources via spectral-estimation techniques
NASA Astrophysics Data System (ADS)
Gabriel, W. F.
1982-06-01
Modern spectral-estimation techniques have achieved a level of performance that attracts interest in applications area such as the tracking of multiple spatial sources. In addition to the original "superresolution' capability, these techniques offer an apparent 'absence of sidelobes' characteristic and some reasonable solutions to the difficult radar coherent-source problem that involves a phase-dependent SNR (signal-to-noise ratio) penalty. This report reviews the situation briefly, and it discusses a few of the techniques that have been found useful, including natural or synthetic doppler shifts, non-Toeplitz forward-backward subaperture-shift processing, and recent eigenvalue/eigenvector analysis algorithms. The techniques are applied to multiple-source situations that include mixtures of coherent and noncoherent sources of unequal strengths, with either an 8-or a 12-element linear-array sampling aperture. The first test case involves the estimation of six sources, two of which are 95% correlated. The second test case involves a tracking-simulation display example of four moving sources: three are -10dB coherent sources 95% correlated, and the other is a strong 20-dB noncoherent source. These test cases demonstrate the remarkable improvements obtained with the recent estimation techniques, and they point to the possibilities for real-world applications.
NASA Astrophysics Data System (ADS)
Lana, X.; Burgueño, A.; Serra, C.; Martínez, M. D.
2015-09-01
Dry spell lengths, DSL, defined as the number of consecutive days with daily rain amounts below a given threshold, may provide relevant information about drought regimes. Taking advantage of a daily pluviometric database covering a great extension of Europe, a detailed analysis of the multifractality of the dry spell regimes is achieved. At the same time, an autoregressive process is applied with the aim of predicting DSL. A set of parameters, namely Hurst exponent, H, estimated from multifractal spectrum, f(α), critical Hölder exponent, α 0, for which f(α) reaches its maximum value, spectral width, W, and spectral asymmetry, B, permits a first clustering of European rain gauges in terms of the complexity of their DSL series. This set of parameters also allows distinguishing between time series describing fine- or smooth-structure of the DSL regime by using the complexity index, CI. Results of previous monofractal analyses also permits establishing comparisons between smooth-structures, relatively low correlation dimensions, notable predictive instability and anti-persistence of DSL for European areas, sometimes submitted to long droughts. Relationships are also found between the CI and the mean absolute deviation, MAD, and the optimum autoregressive order, OAO, of an ARIMA(p,d,0) autoregressive process applied to the DSL series. The detailed analysis of the discrepancies between empiric and predicted DSL underlines the uncertainty over predictability of long DSL, particularly for the Mediterranean region.
Estimation of canopy water content with MODIS spectral index
NASA Astrophysics Data System (ADS)
Jiang, Zuoning; Li, Lin; Ustin, Susan L.
2009-08-01
Canopy water content is an important variable for forestry and agriculture management. This study was aimed at building calibration models to estimate vegetation canopy (VC) equivalent water thickness (EWT) from high temporal resolution and large areal coverage MODIS images. The models were developed for a semi-arid area in Arizona (SMEX04) and the best one was applied to MODIS images covering a forest area in Southern Indiana. EWT derived from hyperspectral data in the process of atmospheric correction was used for calibrating MODIS spectral indices. Tested in this study were four vegetation indices: Normalized Difference Water Index (NDWI), Shortwave Infrared Water Stress Index (SIWSI), Normalized Difference Vegetation Index (NDVI), and Enhanced Vegetation Index (EVI), which were designed based on either water (NDWI and SIWSI) or chlorophyll absorptions (NDVI and EVI). Validating these indices on field measured EWT for the SMEX04 site resulted in R2 correlations of 0.7547, 0.7509, 0.7299 and 0.7547, respectively. According to regression equations, however, EWT estimated using NDWI and SIWSI shows a slope more close to 1 than those using NDVI and EVI when validated with ground measured EWT, thus showing a better prediction ability than the two chlorophyll indices. The SIWSI-EWT model was chosen to apply to a time series of MODIS images covering the Southern Indiana areas and the relationship of EWT derived from these images to precipitation was examined.
NASA Technical Reports Server (NTRS)
Scaife, Bradley James
1999-01-01
In any satellite communication, the Doppler shift associated with the satellite's position and velocity must be calculated in order to determine the carrier frequency. If the satellite state vector is unknown then some estimate must be formed of the Doppler-shifted carrier frequency. One elementary technique is to examine the signal spectrum and base the estimate on the dominant spectral component. If, however, the carrier is spread (as in most satellite communications) this technique may fail unless the chip rate-to-data rate ratio (processing gain) associated with the carrier is small. In this case, there may be enough spectral energy to allow peak detection against a noise background. In this thesis, we present a method to estimate the frequency (without knowledge of the Doppler shift) of a spread-spectrum carrier assuming a small processing gain and binary-phase shift keying (BPSK) modulation. Our method relies on an averaged discrete Fourier transform along with peak detection on spectral match filtered data. We provide theory and simulation results indicating the accuracy of this method. In addition, we will describe an all-digital hardware design based around a Motorola DSP56303 and high-speed A/D which implements this technique in real-time. The hardware design is to be used in NMSU's implementation of NASA's demand assignment, multiple access (DAMA) service.
MUSIC for Multidimensional Spectral Estimation: Stability and Super-Resolution
NASA Astrophysics Data System (ADS)
Liao, Wenjing
2015-12-01
This paper presents a performance analysis of the MUltiple SIgnal Classification (MUSIC) algorithm applied on $D$ dimensional single-snapshot spectral estimation while $s$ true frequencies are located on the continuum of a bounded domain. Inspired by the matrix pencil form, we construct a D-fold Hankel matrix from the measurements and exploit its Vandermonde decomposition in the noiseless case. MUSIC amounts to identifying a noise subspace, evaluating a noise-space correlation function, and localizing frequencies by searching the $s$ smallest local minima of the noise-space correlation function. In the noiseless case, $(2s)^D$ measurements guarantee an exact reconstruction by MUSIC as the noise-space correlation function vanishes exactly at true frequencies. When noise exists, we provide an explicit estimate on the perturbation of the noise-space correlation function in terms of noise level, dimension $D$, the minimum separation among frequencies, the maximum and minimum amplitudes while frequencies are separated by two Rayleigh Length (RL) at each direction. As a by-product the maximum and minimum non-zero singular values of the multidimensional Vandermonde matrix whose nodes are on the unit sphere are estimated under a gap condition of the nodes. Under the 2-RL separation condition, if noise is i.i.d. gaussian, we show that perturbation of the noise-space correlation function decays like $\\sqrt{\\log(\\#(\\mathbf{N}))/\\#(\\mathbf{N})}$ as the sample size $\\#(\\mathbf{N})$ increases. When the separation among frequencies drops below 2 RL, our numerical experiments show that the noise tolerance of MUSIC obeys a power law with the minimum separation of frequencies.
Chen, Szi-Wen; Chao, Shih-Chieh
2014-01-01
In this paper, a reweighted ℓ1-minimization based Compressed Sensing (CS) algorithm incorporating the Integral Pulse Frequency Modulation (IPFM) model for spectral estimation of HRV is introduced. Knowing as a novel sensing/sampling paradigm, the theory of CS asserts certain signals that are considered sparse or compressible can be possibly reconstructed from substantially fewer measurements than those required by traditional methods. Our study aims to employ a novel reweighted ℓ1-minimization CS method for deriving the spectrum of the modulating signal of IPFM model from incomplete RR measurements for HRV assessments. To evaluate the performance of HRV spectral estimation, a quantitative measure, referred to as the Percent Error Power (PEP) that measures the percentage of difference between the true spectrum and the spectrum derived from the incomplete RR dataset, was used. We studied the performance of spectral reconstruction from incomplete simulated and real HRV signals by experimentally truncating a number of RR data accordingly in the top portion, in the bottom portion, and in a random order from the original RR column vector. As a result, for up to 20% data truncation/loss the proposed reweighted ℓ1-minimization CS method produced, on average, 2.34%, 2.27%, and 4.55% PEP in the top, bottom, and random data-truncation cases, respectively, on Autoregressive (AR) model derived simulated HRV signals. Similarly, for up to 20% data loss the proposed method produced 5.15%, 4.33%, and 0.39% PEP in the top, bottom, and random data-truncation cases, respectively, on a real HRV database drawn from PhysioNet. Moreover, results generated by a number of intensive numerical experiments all indicated that the reweighted ℓ1-minimization CS method always achieved the most accurate and high-fidelity HRV spectral estimates in every aspect, compared with the ℓ1-minimization based method and Lomb's method used for estimating the spectrum of HRV from unevenly sampled RR
Stationary and integrated autoregressive neural network processes.
Trapletti, A; Leisch, F; Hornik, K
2000-10-01
We consider autoregressive neural network (AR-NN) processes driven by additive noise and demonstrate that the characteristic roots of the shortcuts-the standard conditions from linear time-series analysis-determine the stochastic behavior of the overall AR-NN process. If all the characteristic roots are outside the unit circle, then the process is ergodic and stationary. If at least one characteristic root lies inside the unit circle, then the process is transient. AR-NN processes with characteristic roots lying on the unit circle exhibit either ergodic, random walk, or transient behavior. We also analyze the class of integrated AR-NN (ARI-NN) processes and show that a standardized ARI-NN process "converges" to a Wiener process. Finally, least-squares estimation (training) of the stationary models and testing for nonstationarity is discussed. The estimators are shown to be consistent, and expressions on the limiting distributions are given. PMID:11032041
Biomass estimator for NIR image with a few additional spectral band images taken from light UAS
NASA Astrophysics Data System (ADS)
Pölönen, Ilkka; Salo, Heikki; Saari, Heikki; Kaivosoja, Jere; Pesonen, Liisa; Honkavaara, Eija
2012-05-01
A novel way to produce biomass estimation will offer possibilities for precision farming. Fertilizer prediction maps can be made based on accurate biomass estimation generated by a novel biomass estimator. By using this knowledge, a variable rate amount of fertilizers can be applied during the growing season. The innovation consists of light UAS, a high spatial resolution camera, and VTT's novel spectral camera. A few properly selected spectral wavelengths with NIR images and point clouds extracted by automatic image matching have been used in the estimation. The spectral wavelengths were chosen from green, red, and NIR channels.
Estimation of spectral distribution of sky radiance using a commercial digital camera.
Saito, Masanori; Iwabuchi, Hironobu; Murata, Isao
2016-01-10
Methods for estimating spectral distribution of sky radiance from images captured by a digital camera and for accurately estimating spectral responses of the camera are proposed. Spectral distribution of sky radiance is represented as a polynomial of the wavelength, with coefficients obtained from digital RGB counts by linear transformation. The spectral distribution of radiance as measured is consistent with that obtained by spectrometer and radiative transfer simulation for wavelengths of 430-680 nm, with standard deviation below 1%. Preliminary applications suggest this method is useful for detecting clouds and studying the relation between irradiance at the ground and cloud distribution. PMID:26835780
Estimation of Canopy Foliar Biomass with Spectral Reflectance Measurements
Technology Transfer Automated Retrieval System (TEKTRAN)
Canopy foliar biomass, defined as the product of leaf dry matter content and leaf area index, is an important measurement for global biogeochemical cycles. This study explores the potential for retrieving foliar biomass in green canopies using a spectral index, the Normalized Dry Matter Index (NDMI)...
Assessment of spectral indicies for crop residue cover estimation
Technology Transfer Automated Retrieval System (TEKTRAN)
The quantification of surficial crop residue cover is important for assessing agricultural tillage practices, rangeland health, and brush fire hazards. The Cellulose Absorption Index (CAI) and the Shortwave Infrared Normalized Difference Residue Index (SINDRI) are two spectral indices that have show...
Correlation autoregressive processes with application to helicopter noise
NASA Technical Reports Server (NTRS)
Hardin, J. C.; Miamee, A. G.
1990-01-01
This paper introduces a new class of random processes X(t), the autocorrelations R sub x (t1, t2) of which satisfy a linear relation for all t1 and t2 in some interval of the time axis. Such random processes are denoted as 'correlation-autoregressive'. This class is shown to include the familiar stationary and periodically correlated processes as well as many other, both harmonizable and nonharmonizable, nonstationary processes. When a process is correlation-autoregressive for all times and harmonizable, its two-dimensional power spectral density is shown to take a particularly simple form. The relationship of such processes to the class of stationary processes is examined. In addition, the application of such processes in the analysis of typical helicopter noise signals is described.
USE OF SPECTRAL ANGLE MAPPER (SAM) AND HYPERSPECTRAL IMAGERY FOR YIELD ESTIMATION
Technology Transfer Automated Retrieval System (TEKTRAN)
Vegetation indices (VIs) derived from remotely sensed imagery are commonly used to estimate crop yields. Spectral angle mapper (SAM) provides an alternative approach to quantifying the spectral differences among all pixels in imagery and therefore has the potential for mapping yield variability. The...
Xia, Peng; Shimozato, Yuki; Ito, Yasunori; Tahara, Tatsuki; Kakue, Takashi; Awatsuji, Yasuhiro; Nishio, Kenzo; Ura, Shogo; Kubota, Toshihiro; Matoba, Osamu
2011-12-01
We propose a color digital holography by using spectral estimation technique to improve the color reproduction of objects. In conventional color digital holography, there is insufficient spectral information in holograms, and the color of the reconstructed images depend on only reflectances at three discrete wavelengths used in the recording of holograms. Therefore the color-composite image of the three reconstructed images is not accurate in color reproduction. However, in our proposed method, the spectral estimation technique was applied, which has been reported in multispectral imaging. According to the spectral estimation technique, the continuous spectrum of object can be estimated and the color reproduction is improved. The effectiveness of the proposed method was confirmed by a numerical simulation and an experiment, and, in the results, the average color differences are decreased from 35.81 to 7.88 and from 43.60 to 25.28, respectively. PMID:22193005
ESTIMATION OF RESPONSE-SPECTRAL VALUES AS FUNCTIONS OF MAGNITUDE, DISTANCE, AND SITE CONDITIONS.
Joyner, W.B.; Boore, D.M.
1983-01-01
Horizontal pseudo-velocity response was analyzed for twelve shallow earthquakes in western North America. Estimation of response-spectral values was related to magnitude, distance and site conditions. Errors in the methods are analyzed.
Assessment of Spectral Indices for Crop Residue Cover Estimation
Technology Transfer Automated Retrieval System (TEKTRAN)
Agricultural soils are an important terrestrial carbon (C) stock, accounting for a significant portion of global C estimates. Soil tillage method is important in agricultural C sequestration models. Traditional intensive tillage systems greatly disturb the soil and have been shown to deplete the so...
Spectral Reflectance Estimates of Surface Soil Physical and Chemical Properties
Technology Transfer Automated Retrieval System (TEKTRAN)
Optical diffuse reflectance sensing in visible and near-infrared wavelength ranges is one approach to rapidly quantify soil properties for site-specific management. The objectives of this study were (1) to determine the accuracy of the reflectance approach for estimating physical and chemical proper...
Spectral responsivity estimation and noise effect analysis for digital imaging systems
NASA Astrophysics Data System (ADS)
Chang, Gao-Wei; Kuo, Hung-Zen; Tu, Chung-Fan
2004-02-01
The determination of spectral responsivities plays a significant role in analyzing and predicting the performance of digital imaging systems for remote sensing. For example, given the spectral response functions, we can readily obtain the colorimetric data from a camera corresponding to the remote illuminated objects. In this paper, we develop a filter-based optical system to estimate these functions. The design objective of this system is to effectively select a limited amount of spectral (or broadband) filters to characterize the spectral features of color imaging processes, which are contaminated with noise, so that the spectral response functions can be estimated with satisfactory accuracy. In our approach, a theoretical study is first presented to pave the way for this work, and then we propose a filter selection method based on the technique of orthogonal-triangular (QR) decomposition with column pivoting, called QRCP method. This method involves QR computations and a column permutation process, which determines a permutation matrix to conduct the subset (or filter) selection. Experimental results reveal that the proposed technique is truly consistent with the theoretical study on filter selections. As expected, the optical system with the filters selected from the QRCP method is much less sensitive to noise than those with other spectral filters from different selections. It turns out that our approach is an effective way to implement the optical system for estimating spectral responsivities of digital imaging systems.
NASA Technical Reports Server (NTRS)
Freedman, Ellis; Ryan, Robert; Pagnutti, Mary; Holekamp, Kara; Gasser, Gerald; Carver, David; Greer, Randy
2007-01-01
Spectral Dark Subtraction (SDS) provides good ground reflectance estimates across a variety of atmospheric conditions with no knowledge of those conditions. The algorithm may be sensitive to errors from stray light, calibration, and excessive haze/water vapor. SDS seems to provide better estimates than traditional algorithms using on-site atmospheric measurements much of the time.
On Spectral Classification and Astrophysical Parameter Estimation for Galactic Surveys
NASA Astrophysics Data System (ADS)
Re Fiorentin, Paola; Bailer-Jones, Coryn A. L.; Beers, Timothy C.; Zwitter, Tomaž
2008-12-01
We present several strategies that are being developed in order to classify and parameterize individual stars observed by Galactic surveys, and illustrate some results obtained from spectra obtained by the RAdial Velocity Experiment (RAVE) and the Sloan Digital Sky Survey (SDSS/SEGUE). We demonstrate the efficiency of our models for discrete source classification and stellar atmospheric parameter estimation (effective temperature, surface gravity, and metallicity), which use supervised machine learning algorithms along with a principal component analysis front-end compression phase that also enables knowledge discovery.
NASA Astrophysics Data System (ADS)
Funamizu, Hideki; Tokuno, Yuta; Aizu, Yoshihisa
2016-06-01
We investigate the estimation of spectral transmittance curves in color digital holographic microscopy using speckle illuminations. In color digital holography, it has the disadvantage in that the color-composite image gives poor color information due to the use of lasers with the two or three wavelengths. To overcome this disadvantage, the Wiener estimation method and an averaging process using multiple holograms are applied to color digital holographic microscopy. Estimated spectral transmittance and color-composite images are shown to indicate the usefulness of the proposed method.
[Estimation of Hunan forest carbon density based on spectral mixture analysis of MODIS data].
Yan, En-ping; Lin, Hui; Wang, Guang-xing; Chen, Zhen-xiong
2015-11-01
With the fast development of remote sensing technology, combining forest inventory sample plot data and remotely sensed images has become a widely used method to map forest carbon density. However, the existence of mixed pixels often impedes the improvement of forest carbon density mapping, especially when low spatial resolution images such as MODIS are used. In this study, MODIS images and national forest inventory sample plot data were used to conduct the study of estimation for forest carbon density. Linear spectral mixture analysis with and without constraint, and nonlinear spectral mixture analysis were compared to derive the fractions of different land use and land cover (LULC) types. Then sequential Gaussian co-simulation algorithm with and without the fraction images from spectral mixture analyses were employed to estimate forest carbon density of Hunan Province. Results showed that 1) Linear spectral mixture analysis with constraint, leading to a mean RMSE of 0.002, more accurately estimated the fractions of LULC types than linear spectral and nonlinear spectral mixture analyses; 2) Integrating spectral mixture analysis model and sequential Gaussian co-simulation algorithm increased the estimation accuracy of forest carbon density to 81.5% from 74.1%, and decreased the RMSE to 5.18 from 7.26; and 3) The mean value of forest carbon density for the province was 30.06 t · hm(-2), ranging from 0.00 to 67.35 t · hm(-2). This implied that the spectral mixture analysis provided a great potential to increase the estimation accuracy of forest carbon density on regional and global level. PMID:26915200
Effect of spectral shape on acoustic fatigue life estimates
NASA Technical Reports Server (NTRS)
Miles, R. N.
1992-01-01
Methods for estimating fatigue life due to random loading are briefly reviewed. These methods include a probabilistic approach in which the expected value of the rate of damage accumulation is computed by integrating over the probability density of damaging events and a method which consists of analyzing the response time history to count damaging events. It is noted that it is necessary to employ a time domain approach to perform Rainflow counting, while simple peak counting may be accomplished using the probabilistic method. Data obtained indicate that Rainflow counting produces significantly different fatigue life predictions than other methods that are commonly used in acoustic fatigue predictions. When low-frequency oscillations are present in a signal along with high-frequency components, peak counting will produce substantially shorter fatigue lives than Rainflow counting. It is concluded that Rainflow counting is capable of providing reliable fatigue life predictions for acoustic fatigue studies.
Intercepted photosynthetically active radiation estimated by spectral reflectance
NASA Technical Reports Server (NTRS)
Hatfield, J. L.; Asrar, G.; Kanemasu, E. T.
1984-01-01
Interception of photosynthetically active radiation (PAR) was evaluated relative to greenness and normalized difference (MSS (7-5)/(7+5) for five planting dates of wheat for 1978-79 and 1979-80 at Phoenix, Arizona. Intercepted PAR was calculated from leaf area index and stage of growth. Linear relatinships were found with greeness and normalized difference with separate relatinships describing growth and senescence of the crop. Normalized difference was significantly better than greenness for all planting dates. For the leaf area growth portion of the season the relation between PAR interception and normalized difference was the same over years and planting dates. For the leaf senescence phase the relationships showed more variability due to the lack of data on light interception in sparse and senescing canopies. Normalized difference could be used to estimate PAR interception throughout a growing season.
Spectral estimates of intercepted solar radiation by corn and soybean canopies
NASA Technical Reports Server (NTRS)
Gallo, K. P.; Brooks, C. C.; Daughtry, C. S. T.; Bauer, M. E.; Vanderbilt, V. C.
1982-01-01
Attention is given to the development of methods for combining spectral and meteorological data in crop yield models which are capable of providing accurate estimates of crop condition and yields throughout the growing season. The present investigation is concerned with initial tests of these concepts using spectral and agronomic data acquired in controlled experiments. The data were acquired at the Purdue University Agronomy Farm, 10 km northwest of West Lafayette, Indiana. Data were obtained throughout several growing seasons for corn and soybeans. Five methods or models for predicting yields were examined. On the basis of the obtained results, it is concluded that estimating intercepted solar radiation using spectral data is a viable approach for merging spectral and meteorological data in crop yield models.
Spectral estimation of made-up skin color under various conditions
NASA Astrophysics Data System (ADS)
Doi, Motonori; Ohtsuki, Rie; Tominaga, Shoji
2006-01-01
A method is proposed for estimating the spectral reflectance of made-up skin color under various conditions including the undesirable colored skin. The color of dark spot is caused by increasing the component of melanin. The reddish skin is caused by the increase of hemoglobin. Our method uses the Kubelka-Munk theory to calculate the surface spectral reflectance human skin. This theory calculates the reflectance and transmittance of the light passing through a turbid medium from the absorption and scattering of the medium. The spectral reflectance of made-up skin is estimated by adjusting parameters of the thickness of the makeup layer. The proposed estimation method is evaluated on an experiment in detail. First, we measure the spectral reflectance of facial skin under the three conditions of normal skin, undesirable skin, and made-up skin. The undesirable skin includes stain, suntan or ruddy skin. The made-up skin means the skin with foundation on the normal skin, the stain, the suntan and the ruddy skin. Second, we estimate the spectral reflectance of made-up skins from the reflectance of bare skins and optical characteristics of foundations. Good coincidence between the estimated reflectance and the direct measurement shows the feasibility of the proposed method.
NASA Astrophysics Data System (ADS)
Kira, Oz; Linker, Raphael; Gitelson, Anatoly
2015-06-01
Leaf pigment content provides valuable insight into the productivity, physiological and phenological status of vegetation. Measurement of spectral reflectance offers a fast, nondestructive method for pigment estimation. A number of methods were used previously for estimation of leaf pigment content, however, spectral bands employed varied widely among the models and data used. Our objective was to find informative spectral bands in three types of models, vegetation indices (VI), neural network (NN) and partial least squares (PLS) regression, for estimating leaf chlorophyll (Chl) and carotenoids (Car) contents of three unrelated tree species and to assess the accuracy of the models using a minimal number of bands. The bands selected by PLS, NN and VIs were in close agreement and did not depend on the data used. The results of the uninformative variable elimination PLS approach, where the reliability parameter was used as an indicator of the information contained in the spectral bands, confirmed the bands selected by the VIs, NN, and PLS models. All three types of models were able to accurately estimate Chl content with coefficient of variation below 12% for all three species with VI showing the best performance. NN and PLS using reflectance in four spectral bands were able to estimate accurately Car content with coefficient of variation below 14%. The quantitative framework presented here offers a new way of estimating foliar pigment content not requiring model re-parameterization for different species. The approach was tested using the spectral bands of the future Sentinel-2 satellite and the results of these simulations showed that accurate pigment estimation from satellite would be possible.
Spectral Approach to Optimal Estimation of the Global Average Temperature.
NASA Astrophysics Data System (ADS)
Shen, Samuel S. P.; North, Gerald R.; Kim, Kwang-Y.
1994-12-01
Making use of EOF analysis and statistical optimal averaging techniques, the problem of random sampling error in estimating the global average temperature by a network of surface stations has been investigated. The EOF representation makes it unnecessary to use simplified empirical models of the correlation structure of temperature anomalies. If an adjustable weight is assigned to each station according to the criterion of minimum mean-square error, a formula for this error can be derived that consists of a sum of contributions from successive EOF modes. The EOFs were calculated from both observed data and a noise-forced EBM for the problem of one-year and five-year averages. The mean square statistical sampling error depends on the spatial distribution of the stations, length of the averaging interval, and the choice of the weight for each station data stream. Examples used here include four symmetric configurations of 4 × 4, 6 × 4, 9 × 7, and 20 × 10 stations and the Angell-Korshover configuration. Comparisons with the 100-yr U.K. dataset show that correlations for the time series of the global temperature anomaly average between the full dataset and this study's sparse configurations are rather high. For example, the 63-station Angell-Korshover network with uniform weighting explains 92.7% of the total variance, whereas the same network with optimal weighting can lead to 97.8% explained total variance of the U.K. dataset.
Spectral approach to optimal estimation of the global average temperature
Shen, S.S.P.; North, G.R.; Kim, K.Y.
1994-12-01
Making use of EOF analysis and statistical optimal averaging techniques, the problem of random sampling error in estimating the global average temperature by a network of surface stations has been investigated. The EOF representation makes it unnecessary to use simplified empirical models of the correlation structure of temperature anomalies. If an adjustable weight is assigned to each station according to the criterion of minimum mean-square error, a formula for this error can be derived that consists of a sum of contributions from successive EOF modes. The EOFs were calculated from both observed data a noise-forced EBM for the problem of one-year and five-year averages. The mean square statistical sampling error depends on the spatial distribution of the stations, length of the averaging interval, and the choice of the weight for each station data stream. Examples used here include four symmetric configurations of 4 X 4, 5 X 4, 9 X 7, and 20 X 10 stations and the Angell-Korshover configuration. Comparisons with the 100-yr U.K. dataset show that correlations for the time series of the global temperature anomaly average between the full dataset and this study`s sparse configurations are rather high. For example, the 63-station Angell-Korshover network with uniform weighting explains 92.7% of the total variance, whereas the same network with optimal weighting can lead to 97.8% explained total variance of the U.K. dataset. 27 refs., 5 figs., 4 tabs.
Preliminary evaluation of spectral, normal and meteorological crop stage estimation approaches
NASA Technical Reports Server (NTRS)
Cate, R. B.; Artley, J. A.; Doraiswamy, P. C.; Hodges, T.; Kinsler, M. C.; Phinney, D. E.; Sestak, M. L. (Principal Investigator)
1980-01-01
Several of the projects in the AgRISTARS program require crop phenology information, including classification, acreage and yield estimation, and detection of episodal events. This study evaluates several crop calendar estimation techniques for their potential use in the program. The techniques, although generic in approach, were developed and tested on spring wheat data collected in 1978. There are three basic approaches to crop stage estimation: historical averages for an area (normal crop calendars), agrometeorological modeling of known crop-weather relationships agrometeorological (agromet) crop calendars, and interpretation of spectral signatures (spectral crop calendars). In all, 10 combinations of planting and biostage estimation models were evaluated. Dates of stage occurrence are estimated with biases between -4 and +4 days while root mean square errors range from 10 to 15 days. Results are inconclusive as to the superiority of any of the models and further evaluation of the models with the 1979 data set is recommended.
Cao, Dingcai; Barrionuevo, Pablo A
2015-03-01
The intrinsic circadian clock requires photoentrainment to synchronize the 24-hour solar day. Therefore, light stimulation is an important component of chronobiological research. Currently, the chronobiological research field overwhelmingly uses photopic illuminance that is based on the luminous efficiency function, V(λ), to quantify light levels. However, recent discovery of intrinsically photosensitive retinal ganglion cells (ipRGCs), which are activated by self-contained melanopsin photopigment and also by inputs from rods and cones, makes light specification using a one-dimensional unit inadequate. Since the current understanding of how different photoreceptor inputs contribute to the circadian system through ipRGCs is limited, it is recommended to specify light in terms of the excitations of five photoreceptors (S-, M-, L-cones, rods and ipRGCs; Lucas et al., 2014). In the current study, we assessed whether the spectral outputs from a commercially available spectral watch (i.e. Actiwatch Spectrum) could be used to estimate photoreceptor excitations. Based on the color sensor spectral sensitivity functions from a previously published work, as well as from our measurements, we computed spectral outputs in the long-wavelength range (R), middle-wavelength range (G), short-wavelength range (B) and broadband range (W) under 52 CIE illuminants (25 daylight illuminants, 27 fluorescent lights). We also computed the photoreceptor excitations for each illuminant using human photoreceptor spectral sensitivity functions. Linear regression analyses indicated that the Actiwatch spectral outputs could predict photoreceptor excitations reliably, under the assumption of linear responses of the Actiwatch color sensors. In addition, R, G, B outputs could classify illuminant types (fluorescent versus daylight illuminants) satisfactorily. However, the assessment of actual Actiwatch recording under several testing light sources showed that the spectral outputs were subject to
Power spectral density estimation by spline smoothing in the frequency domain
NASA Technical Reports Server (NTRS)
Defigueiredo, R. J. P.; Thompson, J. R.
1972-01-01
An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying FFT techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.
Power spectral density estimation by spline smoothing in the frequency domain.
NASA Technical Reports Server (NTRS)
De Figueiredo, R. J. P.; Thompson, J. R.
1972-01-01
An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying Fast Fourier Transform techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.-
Spectral estimation of gapped data and SAR imaging with angular diversity
NASA Astrophysics Data System (ADS)
Larsson, Erik G.; Li, Jian; Stoica, Peter; Liu, Guoqing; Williams, Robert L.
2001-08-01
The Amplitude and Phase EStimation (APES) approach to amplitude spectrum estimation has been receiving considerably attention recently. We develop an extension of APES for the spectral estimation of gapped (incomplete) data and apply it to synthetic aperture radar (SAR) imaging with angular diversity. It has recently been shown that APES minimizes a certain least-squares criterion with respect to the estimate of the spectrum. Our new algorithm is called gapped-data APES and is based on minimizing this criterion with respect to the missing data as well. Numerical results are presented to demonstrate the effectiveness of the proposed algorithm and its applicability to SAR imaging with angular diversity.
Daniell method for power spectral density estimation in atomic force microscopy.
Labuda, Aleksander
2016-03-01
An alternative method for power spectral density (PSD) estimation--the Daniell method--is revisited and compared to the most prevalent method used in the field of atomic force microscopy for quantifying cantilever thermal motion--the Bartlett method. Both methods are shown to underestimate the Q factor of a simple harmonic oscillator (SHO) by a predictable, and therefore correctable, amount in the absence of spurious deterministic noise sources. However, the Bartlett method is much more prone to spectral leakage which can obscure the thermal spectrum in the presence of deterministic noise. By the significant reduction in spectral leakage, the Daniell method leads to a more accurate representation of the true PSD and enables clear identification and rejection of deterministic noise peaks. This benefit is especially valuable for the development of automated PSD fitting algorithms for robust and accurate estimation of SHO parameters from a thermal spectrum. PMID:27036781
Time series autoregressive integrated moving average modeling of test-day milk yields of dairy ewes.
Macciotta, N P; Cappio-Borlino, A; Pulina, G
2000-05-01
Monthly test-day milk yields of 1200 dairy Sarda ewes were analyzed by time-series methods. Autocorrelation functions were calculated for lactations within parity classes and altitude of location of flocks. Spectral analysis of the successions of data was developed by Fourier transformation, and different Box-Jenkins autoregressive integrated moving average models were fitted. The separation of deterministic and stochastic components highlighted the autoregressive feature of milk production pattern. The forecasting power of autoregressive integrated moving average models was tested by predicting total milk production for a standardized lactation length of 225 d from only a few test-day records. Results indicated a greater forecasting capacity in comparison with standard methods and suggested further development of time-series analysis for studying lactation curves with more sophisticated methods, such as wavelet decomposition and neural network models. PMID:10821585
Spectral Estimation Model Construction of Heavy Metals in Mining Reclamation Areas
Dong, Jihong; Dai, Wenting; Xu, Jiren; Li, Songnian
2016-01-01
The study reported here examined, as the research subject, surface soils in the Liuxin mining area of Xuzhou, and explored the heavy metal content and spectral data by establishing quantitative models with Multivariable Linear Regression (MLR), Generalized Regression Neural Network (GRNN) and Sequential Minimal Optimization for Support Vector Machine (SMO-SVM) methods. The study results are as follows: (1) the estimations of the spectral inversion models established based on MLR, GRNN and SMO-SVM are satisfactory, and the MLR model provides the worst estimation, with R2 of more than 0.46. This result suggests that the stress sensitive bands of heavy metal pollution contain enough effective spectral information; (2) the GRNN model can simulate the data from small samples more effectively than the MLR model, and the R2 between the contents of the five heavy metals estimated by the GRNN model and the measured values are approximately 0.7; (3) the stability and accuracy of the spectral estimation using the SMO-SVM model are obviously better than that of the GRNN and MLR models. Among all five types of heavy metals, the estimation for cadmium (Cd) is the best when using the SMO-SVM model, and its R2 value reaches 0.8628; (4) using the optimal model to invert the Cd content in wheat that are planted on mine reclamation soil, the R2 and RMSE between the measured and the estimated values are 0.6683 and 0.0489, respectively. This result suggests that the method using the SMO-SVM model to estimate the contents of heavy metals in wheat samples is feasible. PMID:27367708
Spectral Estimation Model Construction of Heavy Metals in Mining Reclamation Areas.
Dong, Jihong; Dai, Wenting; Xu, Jiren; Li, Songnian
2016-01-01
The study reported here examined, as the research subject, surface soils in the Liuxin mining area of Xuzhou, and explored the heavy metal content and spectral data by establishing quantitative models with Multivariable Linear Regression (MLR), Generalized Regression Neural Network (GRNN) and Sequential Minimal Optimization for Support Vector Machine (SMO-SVM) methods. The study results are as follows: (1) the estimations of the spectral inversion models established based on MLR, GRNN and SMO-SVM are satisfactory, and the MLR model provides the worst estimation, with R² of more than 0.46. This result suggests that the stress sensitive bands of heavy metal pollution contain enough effective spectral information; (2) the GRNN model can simulate the data from small samples more effectively than the MLR model, and the R² between the contents of the five heavy metals estimated by the GRNN model and the measured values are approximately 0.7; (3) the stability and accuracy of the spectral estimation using the SMO-SVM model are obviously better than that of the GRNN and MLR models. Among all five types of heavy metals, the estimation for cadmium (Cd) is the best when using the SMO-SVM model, and its R² value reaches 0.8628; (4) using the optimal model to invert the Cd content in wheat that are planted on mine reclamation soil, the R² and RMSE between the measured and the estimated values are 0.6683 and 0.0489, respectively. This result suggests that the method using the SMO-SVM model to estimate the contents of heavy metals in wheat samples is feasible. PMID:27367708
An Overdetermined System for Improved Autocorrelation Based Spectral Moment Estimator Performance
NASA Technical Reports Server (NTRS)
Keel, Byron M.
1996-01-01
Autocorrelation based spectral moment estimators are typically derived using the Fourier transform relationship between the power spectrum and the autocorrelation function along with using either an assumed form of the autocorrelation function, e.g., Gaussian, or a generic complex form and applying properties of the characteristic function. Passarelli has used a series expansion of the general complex autocorrelation function and has expressed the coefficients in terms of central moments of the power spectrum. A truncation of this series will produce a closed system of equations which can be solved for the central moments of interest. The autocorrelation function at various lags is estimated from samples of the random process under observation. These estimates themselves are random variables and exhibit a bias and variance that is a function of the number of samples used in the estimates and the operational signal-to-noise ratio. This contributes to a degradation in performance of the moment estimators. This dissertation investigates the use autocorrelation function estimates at higher order lags to reduce the bias and standard deviation in spectral moment estimates. In particular, Passarelli's series expansion is cast in terms of an overdetermined system to form a framework under which the application of additional autocorrelation function estimates at higher order lags can be defined and assessed. The solution of the overdetermined system is the least squares solution. Furthermore, an overdetermined system can be solved for any moment or moments of interest and is not tied to a particular form of the power spectrum or corresponding autocorrelation function. As an application of this approach, autocorrelation based variance estimators are defined by a truncation of Passarelli's series expansion and applied to simulated Doppler weather radar returns which are characterized by a Gaussian shaped power spectrum. The performance of the variance estimators determined
NASA Astrophysics Data System (ADS)
Li, Zuchuan; Li, Lin; Song, Kaishan; Cassar, Nicolas
2013-03-01
Through its influence on the structure of pelagic ecosystems, phytoplankton size distribution (pico-, nano-, and micro-plankton) is believed to play a key role in "the biological pump." In this paper, an algorithm is proposed to estimate phytoplankton size fractions (PSF) for micro-, nano-, and pico-plankton (fm, fn, and fp, respectively) from the spectral features of remote-sensing data. From remote-sensing reflectance spectrum (Rrs(λ)), the algorithm constructs four types of spectral features: a normalized Rrs(λ), band ratios, continuum-removed spectra, and spectral curvatures. Using support vector machine recursive feature elimination, the algorithm ranks the constructed spectral features and Rrs(λ) according to their sensitivities to PSF which is then regressed against the sensitive spectral features through support vector regression. The algorithm is validated with (1) simulated Rrs(λ) and PSF, and (2) Rrs(λ) obtained by Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and PSF determined from High-Performance Liquid Chromatography (HPLC) pigments. The validation results show the overall effectiveness of the algorithm in estimating PSF, with R2 of (1) 0.938 (fm) for the simulated SeaWiFS data set; and (2) 0.617 (fm), 0.475 (fn), and 0.587 (fp) for the SeaWiFS satellite data set. The validation results also indicate that continuum-removed spectra and spectral curvatures are the dominant spectral features sensitive to PSF with their wavelengths mainly centered on the pigment-absorption domain. Global spatial distributions of fm, fn, and fp were mapped with monthly SeaWiFS images. Overall, their biogeographical distributions are consistent with our current understanding that pico-plankton account for a large proportion of total phytoplankton biomass in oligotrophic regions, nano-plankton in transitional areas, and micro-plankton in high-productivity regions.
NASA Technical Reports Server (NTRS)
Howell, Leonard W.; Whitaker, Ann F. (Technical Monitor)
2001-01-01
The maximum likelihood procedure is developed for estimating the three spectral parameters of an assumed broken power law energy spectrum from simulated detector responses and their statistical properties investigated. The estimation procedure is then generalized for application to real cosmic-ray data. To illustrate the procedure and its utility, analytical methods were developed in conjunction with a Monte Carlo simulation to explore the combination of the expected cosmic-ray environment with a generic space-based detector and its planned life cycle, allowing us to explore various detector features and their subsequent influence on estimating the spectral parameters. This study permits instrument developers to make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope.
On-line estimation of myoelectric signal spectral parameters and nonstationarities detection.
D'Alessio, T; Knaflitz, M; Balestra, G; Paggi, S
1993-09-01
In this communication, we present a method for detecting nonstationarities of random time series with an approximately Gaussian distribution of amplitudes. This method is suitable for real time implementation. Here we report some results obtained by applying them to a time series of spectral parameters of surface myoelectric signals, collected during voluntary isometric contractions of human muscles. Moreover, we describe the computerized system that we used to implement our detector of nonstationarity. This system is based on the TMS 320C25 DSP chip and realizes on-line estimation and display of spectral parameters, as well as detection of their nonstationarities, featuring a sampling frequency up to 20 k samples/s. A friendly user interface, fully menu driven, allows the user to select different options during the system's operation by means of hot keys. The accuracy of the system was tested by comparing its estimates with those of an off-line system, previously characterized, which we took as a reference. The estimates of spectral parameters obtained by means of the two systems were always consistent. The on-line stationarity detector was able to recognize rates of variation of the spectral parameters as small as 1%/s during contractions lasting 10-15 s. This sensitivity makes it suitable for clinical application. The set of results herein presented has been selected to highlight the main characteristics of the system. PMID:8288289
NASA Astrophysics Data System (ADS)
Michaelsen, Kelly E.; Krishnaswamy, Venkataramanan; Pogue, Brian W.; Poplack, Steven P.; Paulsen, Keith D.
2013-03-01
X-ray image pixel intensity and optical scattering are compared for 11 normal subjects to assess the feasibility of using X-ray imaging as a surrogate for optical scattering in NIR spectral tomography. Digital breast tomosynthesis exams, as well as twenty single point reflectance measurements of optical breast scattering are compared for a wide variety of breast sizes and densities to determine if scattering can be accurately predicted based on x-ray attenuation. If implemented, x-ray based scattering estimation will decrease exam time and cost as well as simplify the design of a newly developed integrated near infrared spectral tomography and digital breast tomosynthesis imaging system.
NASA Astrophysics Data System (ADS)
Shafian, S.; Maas, S. J.; Rajan, N.
2014-12-01
Water resources and agricultural applications require knowledge of crop water use (CWU) over a range of spatial and temporal scales. Due to the spatial density of meteorological stations, the resolution of CWU estimates based on these data is fairly coarse and not particularly suitable or reliable for water resources planning, irrigation scheduling and decision making. Various methods have been developed for quantifying CWU of agricultural crops. In this study, an improved version of the spectral crop coefficient which includes the effects of stomatal closure is applied. Raw digital count (DC) data in the red, near-infrared, and thermal infrared (TIR) spectral bands of Landsat-7 and Landsat-8 imaging sensors are used to construct the TIR-ground cover (GC) pixel data distribution and estimate the effects of stomatal closure. CWU is then estimated by combining results of the spectral crop coefficient approach and the stomatal closer effect. To test this approach, evapotranspiration was measured in 5 agricultural fields in the semi-arid Texas High Plains during the 2013 and 2014 growing seasons and compared to corresponding estimated values of CWU determined using this approach. The results showed that the estimated CWU from this approach was strongly correlated (R2 = 0.79) with observed evapotranspiration. In addition, the results showed that considering the stomatal closer effect in the proposed approach can improve the accuracy of the spectral crop coefficient method. These results suggest that the proposed approach is suitable for operational estimation of evapotranspiration and irrigation scheduling where irrigation is used to replace the daily CWU of a crop.
To center or not to center? Investigating inertia with a multilevel autoregressive model
Hamaker, Ellen L.; Grasman, Raoul P. P. P.
2015-01-01
Whether level 1 predictors should be centered per cluster has received considerable attention in the multilevel literature. While most agree that there is no one preferred approach, it has also been argued that cluster mean centering is desirable when the within-cluster slope and the between-cluster slope are expected to deviate, and the main interest is in the within-cluster slope. However, we show in a series of simulations that if one has a multilevel autoregressive model in which the level 1 predictor is the lagged outcome variable (i.e., the outcome variable at the previous occasion), cluster mean centering will in general lead to a downward bias in the parameter estimate of the within-cluster slope (i.e., the autoregressive relationship). This is particularly relevant if the main question is whether there is on average an autoregressive effect. Nonetheless, we show that if the main interest is in estimating the effect of a level 2 predictor on the autoregressive parameter (i.e., a cross-level interaction), cluster mean centering should be preferred over other forms of centering. Hence, researchers should be clear on what is considered the main goal of their study, and base their choice of centering method on this when using a multilevel autoregressive model. PMID:25688215
Daniell method for power spectral density estimation in atomic force microscopy
NASA Astrophysics Data System (ADS)
Labuda, Aleksander
2016-03-01
An alternative method for power spectral density (PSD) estimation—the Daniell method—is revisited and compared to the most prevalent method used in the field of atomic force microscopy for quantifying cantilever thermal motion—the Bartlett method. Both methods are shown to underestimate the Q factor of a simple harmonic oscillator (SHO) by a predictable, and therefore correctable, amount in the absence of spurious deterministic noise sources. However, the Bartlett method is much more prone to spectral leakage which can obscure the thermal spectrum in the presence of deterministic noise. By the significant reduction in spectral leakage, the Daniell method leads to a more accurate representation of the true PSD and enables clear identification and rejection of deterministic noise peaks. This benefit is especially valuable for the development of automated PSD fitting algorithms for robust and accurate estimation of SHO parameters from a thermal spectrum.
NASA Astrophysics Data System (ADS)
Warren, Russell E.; Vanderbeek, Richard G.
2004-08-01
Detection and estimation of materials in the atmosphere by lidar has heretofore required that the spectral dependence of the relevant cross section coefficients -- backscatter in the case of aerosols and absorptivity for vapors -- be known in advance. While this typically is a reasonable assumption in the case of vapor, the aerosol backscatter coefficients are complicated functions of particle size, shape, and refractive index, and are therefore usually not well characterized a priori. Using incorrect parameters will give biased concentration estimates and impair discrimination ability. This paper describes an approach for estimating both the spectral dependence of the aerosol backscatter and relative concentration range-dependence of a set of materials using multi-wavelength lidar. The approach is based on state-space filtering that applies a Kalman filter in range for concentration, and updates the backscatter spectral estimates through a sequential least-squares algorithm at each time step. The method is illustrated on aerosol-release data of the bio-simulant ovalbumin collected by ECBC during field tests in 2002, as well as synthetic data sets.
NASA Technical Reports Server (NTRS)
Garber, Donald P.
1993-01-01
A probability density function for the variability of ensemble averaged spectral estimates from helicopter acoustic signals in Gaussian background noise was evaluated. Numerical methods for calculating the density function and for determining confidence limits were explored. Density functions were predicted for both synthesized and experimental data and compared with observed spectral estimate variability.
[The Study of the Spectral Model for Estimating Pigment Contents of Tobacco Leaves in Field].
Ren, Xiao; Lao, Cai-lian; Xu, Zhao-li; Jin, Yan; Guo, Yan; Li, Jun-hui; Yang, Yu-hong
2015-06-01
Fast and non-destructive measurements of tobacco leaf pigment contents by spectroscopy in situ in the field has great significance in production guidance for nutrient diagnosis and growth monitoring of tobacco in vegetative growth stage, and it is also very important for the quality evaluation of tobacco leaves in mature stage. The purpose of this study is to estimate the chlorophyll and carotenoid contents of tobacco leaves using tobacco leaf spectrum collected in the field. Reflectance spectrum of tobacco leaves in vegetative growth stage and mature stage were collected in situ in the field and the pigment contents of tobacco leaf samples were measured in this study, taking the tobacco leaf samples collected in each and both stages as modeling sets respectively, and using the methods of support vector machine (SVM) and spectral indice to establish the pigment content estimation models, and then compare the prediction performance of the models built by different methods. The study results indicated that the difference of estimation performance by each stage or mixed stages is not significant. For chlorophyll content, SVM and spectral indice modeling methods can both have a well estimation performance, while for carotenoid content, SVM modeling method has a better estimation performance than spectral indice. The coefficient of determination and the root mean square error of SVM model for estimating tobacco leaf chlorophyll content by each stage were 0.867 6 and 0.014 7, while the coefficient of determination and the root mean square error of SVM model for estimating tobacco leaf chlorophyll content by mixed stages were 0.898 6 and 0.012 3; The coefficient of determination and the root mean square error for estimating tobacco leaf carotenoid content by each stage were 0.861 4 and 0.002 5, while the coefficient of determination and the root mean square error of SVM model for estimating tobacco leaf carotenoid content by mixed stages were 0.839 9 and 0.002 5. The
Performance Evaluation of the Spectral Centroid Downshift Method for Attenuation Estimation
Samimi, Kayvan; Varghese, Tomy
2015-01-01
Estimation of frequency-dependent ultrasonic attenuation is an important aspect of tissue characterization. Along with other acoustic parameters studied in quantitative ultrasound, the attenuation coefficient can be used to differentiate normal and pathological tissue. The spectral centroid downshift (CDS) method is one the most common frequency-domain approaches applied to this problem. In this study, a statistical analysis of this method’s performance was carried out based on a parametric model of the signal power spectrum in the presence of electronic noise. The parametric model used for the power spectrum of received RF data assumes a Gaussian spectral profile for the transmit pulse, and incorporates effects of attenuation, windowing, and electronic noise. Spectral moments were calculated and used to estimate second-order centroid statistics. A theoretical expression for the variance of a maximum likelihood estimator of attenuation coefficient was derived in terms of the centroid statistics and other model parameters, such as transmit pulse center frequency and bandwidth, RF data window length, SNR, and number of regression points. Theoretically predicted estimation variances were compared with experimentally estimated variances on RF data sets from both computer-simulated and physical tissue-mimicking phantoms. Scan parameter ranges for this study were electronic SNR from 10 to 70 dB, transmit pulse standard deviation from 0.5 to 4.1 MHz, transmit pulse center frequency from 2 to 8 MHz, and data window length from 3 to 17 mm. Acceptable agreement was observed between theoretical predictions and experimentally estimated values with differences smaller than 0.05 dB/cm/MHz across the parameter ranges investigated. This model helps predict the best attenuation estimation variance achievable with the CDS method, in terms of said scan parameters. PMID:25965681
[A spectral unmixing method of estimating main minerals abundance of lunar soils].
Yan, Bo-Kun; Li, Jian-Zhong; Gan, Fu-Ping; Yang, Su-Ming; Wang, Run-Sheng
2012-12-01
Estimating minerals abundance from reflectance spectra is one of the fundamental goals of remote sensing lunar exploration, and the main difficulties are the complicated mixing law of minerals spectrum and spectral features being sensitive to several kinds of factors such as topography, particle size and roughness etc. A method based on spectral unmixing was put forward and tested in the present paper. Before spectra are unmixed the spectral continuum is removed for clarifying and strengthening spectral features. The absorption features and reflectance features (the upward curving parts of spectra between absorption features) are integrated for unmixing to improve the unmixing performance. The Hapke model was used to correct unmixing error due to nonlinear mixing of minerals spectra. Forty three mixed spectra of olivine, clinopyroxene, hypersthene and plagioclase were used to validate the above method. The four minerals abundance was estimated under the conditions of being unaware of endmember spectra used to mix, granularity and chemical composition of minerals. Residual error, abundance error and correlation coefficient between retrieved and true abundance were 5.0 Vol%, 14.4 Vol% and 0.92 respectively. The method and result of this paper could be referred in the lunar minerals mapping of imaging spectrometer data such as M3. PMID:23427563
Software algorithm and hardware design for real-time implementation of new spectral estimator
2014-01-01
Background Real-time spectral analyzers can be difficult to implement for PC computer-based systems because of the potential for high computational cost, and algorithm complexity. In this work a new spectral estimator (NSE) is developed for real-time analysis, and compared with the discrete Fourier transform (DFT). Method Clinical data in the form of 216 fractionated atrial electrogram sequences were used as inputs. The sample rate for acquisition was 977 Hz, or approximately 1 millisecond between digital samples. Real-time NSE power spectra were generated for 16,384 consecutive data points. The same data sequences were used for spectral calculation using a radix-2 implementation of the DFT. The NSE algorithm was also developed for implementation as a real-time spectral analyzer electronic circuit board. Results The average interval for a single real-time spectral calculation in software was 3.29 μs for NSE versus 504.5 μs for DFT. Thus for real-time spectral analysis, the NSE algorithm is approximately 150× faster than the DFT. Over a 1 millisecond sampling period, the NSE algorithm had the capability to spectrally analyze a maximum of 303 data channels, while the DFT algorithm could only analyze a single channel. Moreover, for the 8 second sequences, the NSE spectral resolution in the 3-12 Hz range was 0.037 Hz while the DFT spectral resolution was only 0.122 Hz. The NSE was also found to be implementable as a standalone spectral analyzer board using approximately 26 integrated circuits at a cost of approximately $500. The software files used for analysis are included as a supplement, please see the Additional files 1 and 2. Conclusions The NSE real-time algorithm has low computational cost and complexity, and is implementable in both software and hardware for 1 millisecond updates of multichannel spectra. The algorithm may be helpful to guide radiofrequency catheter ablation in real time. PMID:24886214
NASA Astrophysics Data System (ADS)
Osman, Abdalla; Nourledin, Aboelamgd; El-Sheimy, Naser; Theriault, Jim; Campbell, Scott
2009-06-01
The problem of target detection and tracking in the ocean environment has attracted considerable attention due to its importance in military and civilian applications. Sonobuoys are one of the capable passive sonar systems used in underwater target detection. Target detection and bearing estimation are mainly obtained through spectral analysis of received signals. The frequency resolution introduced by current techniques is limited which affects the accuracy of target detection and bearing estimation at a relatively low signal-to-noise ratio (SNR). This research investigates the development of a bearing estimation method using fast orthogonal search (FOS) for enhanced spectral estimation. FOS is employed in this research in order to improve both target detection and bearing estimation in the case of low SNR inputs. The proposed methods were tested using simulated data developed for two different scenarios under different underwater environmental conditions. The results show that the proposed method is capable of enhancing the accuracy for target detection as well as bearing estimation especially in cases of a very low SNR.
A spectral-spatial-dynamic hierarchical Bayesian (SSD-HB) model for estimating soybean yield
NASA Astrophysics Data System (ADS)
Kazama, Yoriko; Kujirai, Toshihiro
2014-10-01
A method called a "spectral-spatial-dynamic hierarchical-Bayesian (SSD-HB) model," which can deal with many parameters (such as spectral and weather information all together) by reducing the occurrence of multicollinearity, is proposed. Experiments conducted on soybean yields in Brazil fields with a RapidEye satellite image indicate that the proposed SSD-HB model can predict soybean yield with a higher degree of accuracy than other estimation methods commonly used in remote-sensing applications. In the case of the SSD-HB model, the mean absolute error between estimated yield of the target area and actual yield is 0.28 t/ha, compared to 0.34 t/ha when conventional PLS regression was applied, showing the potential effectiveness of the proposed model.
NASA Astrophysics Data System (ADS)
Rao, Roshan
2016-04-01
Aerosol radiative forcing estimates with high certainty are required in climate change studies. The approach in estimating the aerosol radiative forcing by using the chemical composition of aerosols is not effective as the chemical composition data with radiative properties are not widely available. We look into the approach where ground based spectral radiation flux measurement is made and along with an Radtiative transfer (RT) model, radiative forcing is estimated. Measurements of spectral flux were made using an ASD spectroradiometer with 350 - 1050 nm wavelength range and a 3nm resolution during around 54 clear-sky days during which AOD range was around 0.01 to 0.7. Simultaneous measurements of black carbon were also made using Aethalometer (Magee Scientific) which ranged from around 1.5 ug/m3 to 8 ug/m3. The primary study involved in understanding the sensitivity of spectral flux due to change in individual aerosol species (Optical properties of Aerosols and Clouds (OPAC) classified aerosol species) using the SBDART RT model. This made us clearly distinguish the influence of different aerosol species on the spectral flux. Following this, a new technique has been introduced to estimate an optically equivalent mixture of aerosol species for the given location. The new method involves matching different combinations of aerosol species in OPAC model and RT model as long as the combination which gives the minimum root mean squared deviation from measured spectral flux is obtained. Using the optically equivalent aerosol mixture and RT model, aerosol radiative forcing is estimated. Also an alternate method to estimate the spectral SSA is discussed. Here, the RT model, the observed spectral flux and spectral AOD is used. Spectral AOD is input to RT model and SSA is varied till the minimum root mean squared difference between observed and simulated spectral flux from RT model is obtained. The methods discussed are limited to clear sky scenes and its accuracy to derive
Spectral estimation of received phase in the presence of amplitude scintillation
NASA Technical Reports Server (NTRS)
Vilnrotter, V. A.; Brown, D. H.; Hurd, W. J.
1988-01-01
A technique is demonstrated for obtaining the spectral parameters of the received carrier phase in the presence of carrier amplitude scintillation, by means of a digital phased locked loop. Since the random amplitude fluctuations generate time-varying loop characteristics, straightforward processing of the phase detector output does not provide accurate results. The method developed here performs a time-varying inverse filtering operation on the corrupted observables, thus recovering the original phase process and enabling accurate estimation of its underlying parameters.
A Real Valued Neural Network Based Autoregressive Energy Detector for Cognitive Radio Application.
Onumanyi, A J; Onwuka, E N; Aibinu, A M; Ugweje, O C; Salami, M J E
2014-01-01
A real valued neural network (RVNN) based energy detector (ED) is proposed and analyzed for cognitive radio (CR) application. This was developed using a known two-layered RVNN model to estimate the model coefficients of an autoregressive (AR) system. By using appropriate modules and a well-designed detector, the power spectral density (PSD) of the AR system transfer function was estimated and subsequent receiver operating characteristic (ROC) curves of the detector generated and analyzed. A high detection performance with low false alarm rate was observed for varying signal to noise ratio (SNR), sample number, and model order conditions. The proposed RVNN based ED was then compared to the simple periodogram (SP), Welch periodogram (WP), multitaper (MT), Yule-Walker (YW), Burg (BG), and covariance (CV) based ED techniques. The proposed detector showed better performance than the SP, WP, and MT while providing better false alarm performance than the YW, BG, and CV. Data provided here support the effectiveness of the proposed RVNN based ED for CR application. PMID:27379318
A Real Valued Neural Network Based Autoregressive Energy Detector for Cognitive Radio Application
Onumanyi, A. J.; Onwuka, E. N.; Aibinu, A. M.; Ugweje, O. C.; Salami, M. J. E.
2014-01-01
A real valued neural network (RVNN) based energy detector (ED) is proposed and analyzed for cognitive radio (CR) application. This was developed using a known two-layered RVNN model to estimate the model coefficients of an autoregressive (AR) system. By using appropriate modules and a well-designed detector, the power spectral density (PSD) of the AR system transfer function was estimated and subsequent receiver operating characteristic (ROC) curves of the detector generated and analyzed. A high detection performance with low false alarm rate was observed for varying signal to noise ratio (SNR), sample number, and model order conditions. The proposed RVNN based ED was then compared to the simple periodogram (SP), Welch periodogram (WP), multitaper (MT), Yule-Walker (YW), Burg (BG), and covariance (CV) based ED techniques. The proposed detector showed better performance than the SP, WP, and MT while providing better false alarm performance than the YW, BG, and CV. Data provided here support the effectiveness of the proposed RVNN based ED for CR application. PMID:27379318
Li, Ying; Wang, Hong; Li, Xiao Bing
2015-01-01
Vegetation is an important part of ecosystem and estimation of fractional vegetation cover is of significant meaning to monitoring of vegetation growth in a certain region. With Landsat TM images and HJ-1B images as data source, an improved selective endmember linear spectral mixture model (SELSMM) was put forward in this research to estimate the fractional vegetation cover in Huangfuchuan watershed in China. We compared the result with the vegetation coverage estimated with linear spectral mixture model (LSMM) and conducted accuracy test on the two results with field survey data to study the effectiveness of different models in estimation of vegetation coverage. Results indicated that: (1) the RMSE of the estimation result of SELSMM based on TM images is the lowest, which is 0.044. The RMSEs of the estimation results of LSMM based on TM images, SELSMM based on HJ-1B images and LSMM based on HJ-1B images are respectively 0.052, 0.077 and 0.082, which are all higher than that of SELSMM based on TM images; (2) the R2 of SELSMM based on TM images, LSMM based on TM images, SELSMM based on HJ-1B images and LSMM based on HJ-1B images are respectively 0.668, 0.531, 0.342 and 0.336. Among these models, SELSMM based on TM images has the highest estimation accuracy and also the highest correlation with measured vegetation coverage. Of the two methods tested, SELSMM is superior to LSMM in estimation of vegetation coverage and it is also better at unmixing mixed pixels of TM images than pixels of HJ-1B images. So, the SELSMM based on TM images is comparatively accurate and reliable in the research of regional fractional vegetation cover estimation. PMID:25905772
Proper orthogonal decomposition-based spectral higher-order stochastic estimation
NASA Astrophysics Data System (ADS)
Baars, Woutijn J.; Tinney, Charles E.
2014-05-01
A unique routine, capable of identifying both linear and higher-order coherence in multiple-input/output systems, is presented. The technique combines two well-established methods: Proper Orthogonal Decomposition (POD) and Higher-Order Spectra Analysis. The latter of these is based on known methods for characterizing nonlinear systems by way of Volterra series. In that, both linear and higher-order kernels are formed to quantify the spectral (nonlinear) transfer of energy between the system's input and output. This reduces essentially to spectral Linear Stochastic Estimation when only first-order terms are considered, and is therefore presented in the context of stochastic estimation as spectral Higher-Order Stochastic Estimation (HOSE). The trade-off to seeking higher-order transfer kernels is that the increased complexity restricts the analysis to single-input/output systems. Low-dimensional (POD-based) analysis techniques are inserted to alleviate this void as POD coefficients represent the dynamics of the spatial structures (modes) of a multi-degree-of-freedom system. The mathematical framework behind this POD-based HOSE method is first described. The method is then tested in the context of jet aeroacoustics by modeling acoustically efficient large-scale instabilities as combinations of wave packets. The growth, saturation, and decay of these spatially convecting wave packets are shown to couple both linearly and nonlinearly in the near-field to produce waveforms that propagate acoustically to the far-field for different frequency combinations.
Proper orthogonal decomposition-based spectral higher-order stochastic estimation
Baars, Woutijn J.; Tinney, Charles E.
2014-05-15
A unique routine, capable of identifying both linear and higher-order coherence in multiple-input/output systems, is presented. The technique combines two well-established methods: Proper Orthogonal Decomposition (POD) and Higher-Order Spectra Analysis. The latter of these is based on known methods for characterizing nonlinear systems by way of Volterra series. In that, both linear and higher-order kernels are formed to quantify the spectral (nonlinear) transfer of energy between the system's input and output. This reduces essentially to spectral Linear Stochastic Estimation when only first-order terms are considered, and is therefore presented in the context of stochastic estimation as spectral Higher-Order Stochastic Estimation (HOSE). The trade-off to seeking higher-order transfer kernels is that the increased complexity restricts the analysis to single-input/output systems. Low-dimensional (POD-based) analysis techniques are inserted to alleviate this void as POD coefficients represent the dynamics of the spatial structures (modes) of a multi-degree-of-freedom system. The mathematical framework behind this POD-based HOSE method is first described. The method is then tested in the context of jet aeroacoustics by modeling acoustically efficient large-scale instabilities as combinations of wave packets. The growth, saturation, and decay of these spatially convecting wave packets are shown to couple both linearly and nonlinearly in the near-field to produce waveforms that propagate acoustically to the far-field for different frequency combinations.
Li, Zenghui; Xu, Bin; Yang, Jian; Song, Jianshe
2015-01-01
This paper focuses on suppressing spectral overlap for sub-band spectral estimation, with which we can greatly decrease the computational complexity of existing spectral estimation algorithms, such as nonlinear least squares spectral analysis and non-quadratic regularized sparse representation. Firstly, our study shows that the nominal ability of the high-order analysis filter to suppress spectral overlap is greatly weakened when filtering a finite-length sequence, because many meaningless zeros are used as samples in convolution operations. Next, an extrapolation-based filtering strategy is proposed to produce a series of estimates as the substitutions of the zeros and to recover the suppression ability. Meanwhile, a steady-state Kalman predictor is applied to perform a linearly-optimal extrapolation. Finally, several typical methods for spectral analysis are applied to demonstrate the effectiveness of the proposed strategy. PMID:25609038
Li, Zenghui; Xu, Bin; Yang, Jian; Song, Jianshe
2015-01-01
This paper focuses on suppressing spectral overlap for sub-band spectral estimation, with which we can greatly decrease the computational complexity of existing spectral estimation algorithms, such as nonlinear least squares spectral analysis and non-quadratic regularized sparse representation. Firstly, our study shows that the nominal ability of the high-order analysis filter to suppress spectral overlap is greatly weakened when filtering a finite-length sequence, because many meaningless zeros are used as samples in convolution operations. Next, an extrapolation-based filtering strategy is proposed to produce a series of estimates as the substitutions of the zeros and to recover the suppression ability. Meanwhile, a steady-state Kalman predictor is applied to perform a linearly-optimal extrapolation. Finally, several typical methods for spectral analysis are applied to demonstrate the effectiveness of the proposed strategy. PMID:25609038
Spectral estimation from laser scanner data for accurate color rendering of objects
NASA Astrophysics Data System (ADS)
Baribeau, Rejean
2002-06-01
Estimation methods are studied for the recovery of the spectral reflectance across the visible range from the sensing at just three discrete laser wavelengths. Methods based on principal component analysis and on spline interpolation are judged based on the CIE94 color differences for some reference data sets. These include the Macbeth color checker, the OSA-UCS color charts, some artist pigments, and a collection of miscellaneous surface colors. The optimal three sampling wavelengths are also investigated. It is found that color can be estimated with average accuracy ΔE94 = 2.3 when optimal wavelengths 455 nm, 540 n, and 610 nm are used.
The use of spectral data in wheat yield estimation - An assessment of techniques explored in LACIE
NASA Technical Reports Server (NTRS)
Stuff, R. G.; Barnett, T. L.
1979-01-01
The object of the paper is to assess the results of the Large Area Crop Inventory Experiment (LACIE) and closely related research on yield estimation techniques based on remote sensing variables. The exploratory research conducted during LACIE substantiated the hypothesis of yield related information contained in Landsat multispectral scanner data and indicated some of its empirical characteristics. It is noted that leaf area and possibly other foliage features can be derived from spectral data for yield estimation through agrometeorological models and that multiple vegetative and grain related features may be discernable by Landsat derived wheat spectra at different points in the crop development.
NASA Technical Reports Server (NTRS)
Vukovich, Fred M.; Toll, David L.; Kennard, Ruth L.
1989-01-01
Surface biophysical estimates were derived from analysis of NOAA Advanced Very High Spectral Resolution (AVHRR) spectral data of the Senegalese area of west Africa. The parameters derived were of solar albedo, spectral visible and near-infrared band reflectance, spectral vegetative index, and ground temperature. Wet and dry linked AVHRR scenes from 1981 through 1985 in Senegal were analyzed for a semi-wet southerly site near Tambacounda and a predominantly dry northerly site near Podor. Related problems were studied to convert satellite derived radiance to biophysical estimates of the land surface. Problems studied were associated with sensor miscalibration, atmospheric and aerosol spatial variability, surface anisotropy of reflected radiation, narrow satellite band reflectance to broad solar band conversion, and ground emissivity correction. The middle-infrared reflectance was approximated with a visible AVHRR reflectance for improving solar albedo estimates. In addition, the spectral composition of solar irradiance (direct and diffuse radiation) between major spectral regions (i.e., ultraviolet, visible, near-infrared, and middle-infrared) was found to be insensitive to changes in the clear sky atmospheric optical depth in the narrow band to solar band conversion procedure. Solar albedo derived estimates for both sites were not found to change markedly with significant antecedent precipitation events or correspondingly from increases in green leaf vegetation density. The bright soil/substrate contributed to a high albedo for the dry related scenes, whereas the high internal leaf reflectance in green vegetation canopies in the near-infrared contributed to high solar albedo for the wet related scenes. The relationship between solar albedo and ground temperature was poor, indicating the solar albedo has little control of the ground temperature. The normalized difference vegetation index (NDVI) and the derived visible reflectance were more sensitive to antecedent
[Vegetation index estimation by chlorophyll content of grassland based on spectral analysis].
Xiao, Han; Chen, Xiu-Wan; Yang, Zhen-Yu; Li, Huai-Yu; Zhu, Han
2014-11-01
Comparing the methods of existing remote sensing research on the estimation of chlorophyll content, the present paper confirms that the vegetation index is one of the most practical and popular research methods. In recent years, the increasingly serious problem of grassland degradation. This paper, firstly, analyzes the measured reflectance spectral curve and its first derivative curve in the grasslands of Songpan, Sichuan and Gongger, Inner Mongolia, conducts correlation analysis between these two spectral curves and chlorophyll content, and finds out the regulation between REP (red edge position) and grassland chlorophyll content, that is, the higher the chlorophyll content is, the higher the REIP (red-edge inflection point) value would be. Then, this paper constructs GCI (grassland chlorophyll index) and selects the most suitable band for retrieval. Finally, this paper calculates the GCI by the use of satellite hyperspectral image, conducts the verification and accuracy analysis of the calculation results compared with chlorophyll content data collected from field of twice experiments. The result shows that for grassland chlorophyll content, GCI has stronger sensitivity than other indices of chlorophyll, and has higher estimation accuracy. GCI is the first proposed to estimate the grassland chlorophyll content, and has wide application potential for the remote sensing retrieval of grassland chlorophyll content. In addition, the grassland chlorophyll content estimation method based on remote sensing retrieval in this paper provides new research ideas for other vegetation biochemical parameters' estimation, vegetation growth status' evaluation and grassland ecological environment change's monitoring. PMID:25752061
White dwarf mass estimation with a new comprehensive X-ray spectral model of intermediate polars
NASA Astrophysics Data System (ADS)
Hayashi, Takayuki; Ishida, Manabu
A white dwarf (WD) mass is important astrophysical quantity because the WD explodes as a type Ia supernova when its mass reaches the Chandrasekhar mass limit of 1.4 solar mass. Many WD masses in intermediate polars (IPs) were measured with their X-ray spectra emitted from plasma flows channeled by strong magnetic fields of the WDs. For the WD mass estimation, multi-temperature X-ray spectral models have been used which made by summing up X-ray spectra emitted from the top to the bottom of the plasma flow. However, in previous studies, distributions of physical quantities such as temperature and density etc., which are base of the X-ray spectral model, were calculated with assumptions of accretion rate per unit area (call "specific accretion rate") a = 1 g cm(-2) s(-1) and cylindrical geometry for the plasma flows. In fact, a part of the WD masses estimated with the X-ray spectral model is not consistent with that dynamically measured. Therefore, we calculated the physical quantity distributions with the dipolar geometry and the wide range of the specific accretion rate a = 0.0001 - 100 g cm(-2) s(-1) . The calculations showed that the geometrical difference changes the physical quantity distributions and the lower specific accretion rate leads softer X-ray spectrum under a critical specific accretion rate. These results clearly indicate that the previous assumptions are not good approximation for low accretion IPs. We made a new spectral model of the plasma flow with our physical quantity distributions and applied that to Suzaku observations of high and low accretion rate IPs V1223 Sagittarii and EX Hydrae. As a results, our WD masses are almost consistent with the those dynamically measured. We will present the summary of our theoretical calculation and X-ray spectral model, and application to the {¥it Suzaku} observations.
NASA Astrophysics Data System (ADS)
Behrangi, Ali
In respond to the community demands, combining microwave (MW) and infrared (IR) estimates of precipitation has been an active area of research since past two decades. The anticipated launching of NASA's Global Precipitation Measurement (GPM) mission and the increasing number of spectral bands in recently launched geostationary platforms will provide greater opportunities for investigating new approaches to combine multi-source information towards improved global high resolution precipitation retrievals. After years of the communities' efforts the limitations of the existing techniques are: (1) Drawbacks of IR-only techniques to capture warm rainfall and screen out no-rain thin cirrus clouds; (2) Grid-box- only dependency of many algorithms with not much effort to capture the cloud textures whether in local or cloud patch scale; (3) Assumption of indirect relationship between rain rate and cloud-top temperature that force high intensity precipitation to any cold cloud; (4) Neglecting the dynamics and evolution of cloud in time; (5) Inconsistent combination of MW and IR-based precipitation estimations due to the combination strategies and as a result of above described shortcomings. This PhD dissertation attempts to improve the combination of data from Geostationary Earth Orbit (GEO) and Low-Earth Orbit (LEO) satellites in manners that will allow consistent high resolution integration of the more accurate precipitation estimates, directly observed through LEO's PMW sensors, into the short-term cloud evolution process, which can be inferred from GEO images. A set of novel approaches are introduced to cope with the listed limitations and is consist of the following four consecutive components: (1) starting with the GEO part and by using an artificial-neural network based method it is demonstrated that inclusion of multi-spectral data can ameliorate existing problems associated with IR-only precipitating retrievals; (2) through development of Precipitation Estimation
Seevers, P.M.; Sadowski, F.C.; Lauer, D.T.
1990-01-01
Retrospective satellite image data were evaluated for their ability to demonstrate the influence of center-pivot irrigation development in western Nebraska on spectral change and climate-related factors for the region. Periodic images of an albedo index and a normalized difference vegetation index (NDVI) were generated from calibrated Landsat multispectral scanner (MSS) data and used to monitor spectral changes associated with irrigation development from 1972 through 1986. The albedo index was not useful for monitoring irrigation development. For the NDVI, it was found that proportions of counties in irrigated agriculture, as discriminated by a threshold, were more highly correlated with reported ground estimates of irrigated agriculture than were county mean greenness values. A similar result was achieved when using coarse resolution Advanced Very High Resolution Radiometer (AVHRR) image data for estimating irrigated agriculture. The NDVI images were used to evaluate a procedure for making areal estimates of actual evapotranspiration (ET) volumes. Estimates of ET volumes for test counties, using reported ground acreages and corresponding standard crop coefficients, were correlated with the estimates of ET volume using crop coefficients scaled to NDVI values and pixel counts of crop areas. These county estimates were made under the assumption that soil water availability was unlimited. For nonirrigated vegetation, this may result in over-estimation of ET volumes. Ground information regarding crop types and acreages are required to derive the NDVI scaling factor. Potential ET, estimated with the Jensen-Haise model, is common to both methods. These results, achieved with both MSS and AVHRR data, show promise for providing climatologically important land surface information for regional and global climate models. ?? 1990 Kluwer Academic Publishers.
Moisture estimation in power transformer oil using acoustic signals and spectral kurtosis
NASA Astrophysics Data System (ADS)
Leite, Valéria C. M. N.; Veloso, Giscard F. C.; Borges da Silva, Luiz Eduardo; Lambert-Torres, Germano; Borges da Silva, Jonas G.; Onofre Pereira Pinto, João
2016-03-01
The aim of this paper is to present a new technique for estimating the contamination by moisture in power transformer insulating oil based on the spectral kurtosis analysis of the acoustic signals of partial discharges (PDs). Basically, in this approach, the spectral kurtosis of the PD acoustic signal is calculated and the correlation between its maximum value and the moisture percentage is explored to find a function that calculates the moisture percentage. The function can be easily implemented in DSP, FPGA, or any other type of embedded system for online moisture monitoring. To evaluate the proposed approach, an experiment is assembled with a piezoelectric sensor attached to a tank, which is filled with insulating oil samples contaminated by different levels of moisture. A device generating electrical discharges is submerged into the oil to simulate the occurrence of PDs. Detected acoustic signals are processed using fast kurtogram algorithm to extract spectral kurtosis values. The obtained data are used to find the fitting function that relates the water contamination to the maximum value of the spectral kurtosis. Experimental results show that the proposed method is suitable for online monitoring system of power transformers.
The Performance of Multilevel Growth Curve Models under an Autoregressive Moving Average Process
ERIC Educational Resources Information Center
Murphy, Daniel L.; Pituch, Keenan A.
2009-01-01
The authors examined the robustness of multilevel linear growth curve modeling to misspecification of an autoregressive moving average process. As previous research has shown (J. Ferron, R. Dailey, & Q. Yi, 2002; O. Kwok, S. G. West, & S. B. Green, 2007; S. Sivo, X. Fan, & L. Witta, 2005), estimates of the fixed effects were unbiased, and Type I…
NASA Technical Reports Server (NTRS)
Kim, Moon S.; Daughtry, C. S. T.; Chappelle, E. W.; Mcmurtrey, J. E.; Walthall, C. L.
1994-01-01
Most remote sensing estimations of vegetation variables such as Leaf Area Index (LAI), Absorbed Photosynthetically Active Radiation (APAR), and phytomass are made using broad band sensors with a bandwidth of approximately 100 nm. However, high resolution spectrometers are available and have not been fully exploited for the purpose of improving estimates of vegetation variables. A study directed to investigate the use of high spectral resolution spectroscopy for remote sensing estimates of APAR in vegetation canopies in the presence of nonphotosynthetic background materials such as soil and leaf litter is presented. A high spectral resolution method defined as the Chlorophyll Absorption Ratio Index (CARI) was developed for minimizing the effects of nonphotosynthetic materials in the remote estimates of APAR. CARI utilizes three bands at 550, 670, and 700 nm with bandwidth of 10 nm. Simulated canopy reflectance of a range of LAI were generated with the SAIL model using measurements of 42 different soil types as canopy background. CARI obtained from the simulated canopy reflectance was compared with the broad band vegetation indices (Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), and Simple Ratio (SR)). CARI reduced the effect of nonphotosynthetic background materials in the assessment of vegetation canopy APAR more effectively than broad band vegetation indices.
Practical energy response estimation of photon counting detectors for spectral X-ray imaging
NASA Astrophysics Data System (ADS)
Kang, Dong-Goo; Lee, Jongha; Sung, Younghun; Lee, SeongDeok
2010-04-01
Spectral X-ray imaging is a promising technique to drastically improve the diagnostic quality of radiography and computed tomography (CT), since it enables material decomposition and/or identification based on the energy dependency of material-specific X-ray attenuation. Unlike the charge-integration based X-ray detectors, photon counting X-ray detectors (PCXDs) can discriminate the energies of incident X-ray photons and thereby multi-energy images can be obtained in single exposure. However, the measured data are not accurate since the spectra of incident X-rays are distorted according to the energy response function (ERF) of a PCXD. Thus ERF should be properly estimated in advance for accurate spectral imaging. This paper presents a simple method for ERF estimation based on a polychromatic X-ray source that is widely used for medical imaging. The method consists of three steps: source spectra measurement, detector spectra reconstruction, and ERF inverse estimation. Real spectra of an X-ray tube are first measured at all kVs by using an X-ray spectrometer. The corresponding detector spectra are obtained by threshold scans. The ERF is then estimated by solving the inverse problem. Simulations are conducted to demonstrate the concept of the proposed method.
NASA Technical Reports Server (NTRS)
Howell, L. W.
2001-01-01
A simple power law model consisting of a single spectral index (alpha-1) is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at knee energy (E(sub k)) to a steeper spectral index alpha-2 > alpha-1 above E(sub k). The maximum likelihood procedure is developed for estimating these three spectral parameters of the broken power law energy spectrum from simulated detector responses. These estimates and their surrounding statistical uncertainty are being used to derive the requirements in energy resolution, calorimeter size, and energy response of a proposed sampling calorimeter for the Advanced Cosmic-ray Composition Experiment for the Space Station (ACCESS). This study thereby permits instrument developers to make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope.
NASA Astrophysics Data System (ADS)
Antunes, Jose; Borsoi, Laurent; Delaune, Xavier; Piteau, Philippe
2016-02-01
In this paper, we propose analytical and numerical straightforward approximate methods to estimate the unknown terms of incomplete spectral or correlation matrices, when the cross-spectra or cross-correlations available from multiple measurements do not cover all pairs of transducer locations. The proposed techniques may be applied whenever the available data includes the auto-spectra at all measurement locations, as well as selected cross-spectra which implicates all measurement locations. The suggested methods can also be used for checking the consistency between the spectral or correlation functions pertaining to measurement matrices, in cases of suspicious data. After presenting the proposed spectral estimation formulations, we discuss their merits and limitations. Then we illustrate their use on a realistic simulation of a multi-supported tube subjected to turbulence excitation from cross-flow. Finally, we show the effectiveness of the proposed techniques by extracting the modal responses of the simulated flow-excited tube, using the SOBI (Second Order Blind Identification) method, from an incomplete response matrix 1
NASA Technical Reports Server (NTRS)
Howell, L. W.; Rose, M. Franklin (Technical Monitor)
2000-01-01
A simple power law model consisting of a single spectral index alpha (sub 1), is believed to be an adequate description of the galactic cosmic ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at knee energy E(sub k) to a steeper spectral index alpha(sub 2) greater than alpha(sub 1) above E(sub k). The maximum likelihood procedure is developed for estimating these three spectral parameters of the broken power law energy spectrum from simulated detector responses. These estimates and their surrounding statistical uncertainty are being used to derive the requirements in energy resolution, calorimeter size, and energy response of a proposed sampling calorimeter for the Advanced Cosmic ray Composition Experiment for the Space Station (ACCESS). This study thereby permits instrument developers to make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope.
Spectral Doppler estimation utilizing 2-D spatial information and adaptive signal processing.
Ekroll, Ingvild K; Torp, Hans; Løvstakken, Lasse
2012-06-01
The trade-off between temporal and spectral resolution in conventional pulsed wave (PW) Doppler may limit duplex/triplex quality and the depiction of rapid flow events. It is therefore desirable to reduce the required observation window (OW) of the Doppler signal while preserving the frequency resolution. This work investigates how the required observation time can be reduced by adaptive spectral estimation utilizing 2-D spatial information obtained by parallel receive beamforming. Four adaptive estimation techniques were investigated, the power spectral Capon (PSC) method, the amplitude and phase estimation (APES) technique, multiple signal classification (MUSIC), and a projection-based version of the Capon technique. By averaging radially and laterally, the required covariance matrix could successfully be estimated without temporal averaging. Useful PW spectra of high resolution and contrast could be generated from ensembles corresponding to those used in color flow imaging (CFI; OW = 10). For a given OW, the frequency resolution could be increased compared with the Welch approach, in cases in which the transit time was higher or comparable to the observation time. In such cases, using short or long pulses with unfocused or focused transmit, an increase in temporal resolution of up to 4 to 6 times could be obtained in in vivo examples. It was further shown that by using adaptive signal processing, velocity spectra may be generated without high-pass filtering the Doppler signal. With the proposed approach, spectra retrospectively calculated from CFI may become useful for unfocused as well as focused imaging. This application may provide new clinical information by inspection of velocity spectra simultaneously from several spatial locations. PMID:22711413
Heasler, Patrick; Posse, Christian; Hylden, Jeff; Anderson, Kevin
2007-01-01
This paper presents a nonlinear Bayesian regression algorithm for detecting and estimating gas plume content from hyper-spectral data. Remote sensing data, by its very nature, is collected under less controlled conditions than laboratory data. As a result, the physics-based model that is used to describe the relationship between the observed remote-sensing spectra, and the terrestrial (or atmospheric) parameters that are estimated is typically littered with many unknown “nuisance” parameters. Bayesian methods are well-suited for this context as they automatically incorporate the uncertainties associated with all nuisance parameters into the error estimates of the parameters of interest. The nonlinear Bayesian regression methodology is illustrated on simulated data from a three-layer model for longwave infrared (LWIR) measurements from a passive instrument. The generated LWIR scenes contain plumes of varying intensities, and this allows estimation uncertainty and probability of detection to be quantified. The results show that this approach should permit more accurate estimation as well as a more reasonable description of estimate uncertainty. Specifically, the methodology produces a standard error that is more realistic than that produced by matched filter estimation.
NASA Technical Reports Server (NTRS)
Hardisky, M.; Klemas, V.
1984-01-01
Spectral radiance data were collected from the ground and from a low altitude aircraft in an attempt to gain some insight into the potential utility of actual Thematic Mapper data for biomass estimation in wetland plant communities. No attempt was made to distinguish individual plant species within brackish marsh plant associations. Rather, it was decided to lump plant species with similar canopy morphologies and then estimate from spectral radiance data the biomass of the group. The rationale for such an approach is that plants with a similar morphology will produce a similar reflecting or absorping surface (i.e., canopy) for incoming electromagnetic radiation. Variations in observed reflectance from different plant communities with a similar canopy morphology are more likely to be a result of biomass differences than a result of differences in canopy architecture. If the hypothesis that plants with a similar morphology exhibit similar reflectance characteristics is true, then biomass can be estimated based on a model for the dominant plant morphology within a plant association and the need for species discrimination has effectively been eliminated.
Modeling of uncertain spectra through stochastic autoregressive systems
NASA Astrophysics Data System (ADS)
Wang, Yiwei; Wang, X. Q.; Mignolet, Marc P.; Yang, Shuchi; Chen, P. C.
2016-03-01
The focus of this investigation is on the formulation and validation of a modeling strategy of the uncertainty that may exist on the specification of the power spectral density of scalar stationary processes and on the spectral matrices of vector ones. These processes may, for example, be forces on a structure originating from natural phenomena which are coarsely modeled (i.e., with epistemic uncertainty) or are specified by parameters unknown (i.e., with aleatoric uncertainty) in the application considered. The propagation of the uncertainty, e.g., to the response of the structure, may be carried out provided that a stochastic model of the uncertainty in the power spectral density/matrix is available from which admissible samples can be efficiently generated. Such a stochastic model will be developed here through an autoregressive-based parameterization of the specified baseline power spectral density/matrix and of its random samples. Autoregressive (AR) models are particularly well suited for this parametrization since their spectra are known to converge to a broad class of spectra (all non-pathological spectra) as the AR order increases. Note that the characterization of these models is not achieved directly in terms of their coefficients but rather in terms of their reflection coefficients which lie (or their eigenvalues in the vector process case) in the domain [0,1) as a necessary and sufficient condition for stability. Maximum entropy concepts are then employed to formulate the distribution of the reflection coefficients in both scalar and vector process case leading to a small set of hyperparameters of the uncertain model. Depending on the information available, these hyperparameters could either be varied in a parametric study format to assess the effects of uncertainty or could be identified, e.g., in a maximum likelihood format, from observed data. The validation and assessment of these concepts is finally achieved on several examples including the
Heasler, Patrick G.; Posse, Christian; Hylden, Jeff L.; Anderson, Kevin K.
2007-06-13
This paper presents a nonlinear Bayesian regression algorithm for the purpose of detecting and estimating gas plume content from hyper-spectral data. Remote sensing data, by its very nature, is collected under less controlled conditions than laboratory data. As a result, the physics-based model that is used to describe the relationship between the observed remotesensing spectra, and the terrestrial (or atmospheric) parameters that we desire to estimate, is typically littered with many unknown "nuisance" parameters (parameters that we are not interested in estimating, but also appear in the model). Bayesian methods are well-suited for this context as they automatically incorporate the uncertainties associated with all nuisance parameters into the error estimates of the parameters of interest. The nonlinear Bayesian regression methodology is illustrated on realistic simulated data from a three-layer model for longwave infrared (LWIR) measurements from a passive instrument. This shows that this approach should permit more accurate estimation as well as a more reasonable description of estimate uncertainty.
Remote estimation of phytoplankton size fractions using the spectral shape of light absorption.
Wang, Shengqiang; Ishizaka, Joji; Hirawake, Toru; Watanabe, Yuji; Zhu, Yuanli; Hayashi, Masataka; Yoo, Sinjae
2015-04-20
Phytoplankton size structure plays an important role in ocean biogeochemical processes. The light absorption spectra of phytoplankton provide a great potential for retrieving phytoplankton size structure because of the strong dependence on the packaging effect caused by phytoplankton cell size and on different pigment compositions related to phytoplankton taxonomy. In this study, we investigated the variability in light absorption spectra of phytoplankton in relation to the size structure. Based on this, a new approach was proposed for estimating phytoplankton size fractions. Our approach use the spectral shape of the normalized phytoplankton absorption coefficient (a(ph)(λ)) through principal component analysis (PCA). Values of a(ph)(λ) were normalized to remove biomass effects, and PCA was conducted to separate the spectral variance of normalized a(ph)(λ) into uncorrelated principal components (PCs). Spectral variations captured by the first four PC modes were used to build relationships with phytoplankton size fractions. The results showed that PCA had powerful ability to capture spectral variations in normalized a(ph)(λ), which were significantly related to phytoplankton size fractions. For both hyperspectral a(ph)(λ) and multiband a(ph)(λ), our approach is applicable. We evaluated our approach using wide in situ data collected from coastal waters and the global ocean, and the results demonstrated a good and robust performance in estimating phytoplankton size fractions in various regions. The model performance was further evaluated by a(ph)(λ) derived from in situ remote sensing reflectance (R(rs)(λ)) with a quasi-analytical algorithm. Using R(rs)(λ) only at six bands, accurate estimations of phytoplankton size fractions were obtained, with R(2) values of 0.85, 0.61, and 0.76, and root mean-square errors of 0.130, 0.126, and 0.112 for micro-, nano-, and picophytoplankton, respectively. Our approach provides practical basis for remote estimation of
An initial model for estimating soybean development stages from spectral data
NASA Technical Reports Server (NTRS)
Henderson, K. E.; Badhwar, G. D.
1982-01-01
A model, utilizing a direct relationship between remotely sensed spectral data and soybean development stage, has been proposed. The model is based upon transforming the spectral data in Landsat bands to greenness values over time and relating the area of this curve to soybean development stage. Soybean development stages were estimated from data acquired in 1978 from research plots at the Purdue University Agronomy Farm as well as Landsat data acquired over sample areas of the U.S. Corn Belt in 1978 and 1979. Analysis of spectral data from research plots revealed that the model works well with reasonable variation in planting date, row spacing, and soil background. The R-squared of calculated U.S. observed development stage exceeded 0.91 for all treatment variables. Using Landsat data the calculated U.S. observed development stage gave an R-squared of 0.89 in 1978 and 0.87 in 1979. No difference in the models performance could be detected between early and late planted fields, small and large fields, or high and low yielding fields.
Ebel, J.E.; Wald, D.J.
2003-01-01
We describe a new probabilistic method that uses observations of modified Mercalli intensity (MMI) from past earthquakes to make quantitative estimates of ground shaking parameters (i.e., peak ground acceleration, peak ground velocity, 5% damped spectral acceleration values, etc.). The method uses a Bayesian approach to make quantitative estimates of the probabilities of different levels of ground motions from intensity data given an earthquake of known location and magnitude. The method utilizes probability distributions from an intensity/ground motion data set along with a ground motion attenuation relation to estimate the ground motion from intensity. The ground motions with the highest probabilities are the ones most likely experienced at the site of the MMI observation. We test the method using MMI/ground motion data from California and published ground motion attenuation relations to estimate the ground motions for several earthquakes: 1999 Hector Mine, California (M7.1); 1988 Saguenay, Quebec (M5.9); and 1982 Gaza, New Hampshire (M4.4). In an example where the method is applied to a historic earthquake, we estimate that the peak ground accelerations associated with the 1727 (M???5.2) earthquake at Newbury, Massachusetts, ranged from 0.23 g at Newbury to 0.06 g at Boston.
NASA Astrophysics Data System (ADS)
Edwards, Matthew C.; Meyer, Renate; Christensen, Nelson
2015-09-01
The standard noise model in gravitational wave (GW) data analysis assumes detector noise is stationary and Gaussian distributed, with a known power spectral density (PSD) that is usually estimated using clean off-source data. Real GW data often depart from these assumptions, and misspecified parametric models of the PSD could result in misleading inferences. We propose a Bayesian semiparametric approach to improve this. We use a nonparametric Bernstein polynomial prior on the PSD, with weights attained via a Dirichlet process distribution, and update this using the Whittle likelihood. Posterior samples are obtained using a blocked Metropolis-within-Gibbs sampler. We simultaneously estimate the reconstruction parameters of a rotating core collapse supernova GW burst that has been embedded in simulated Advanced LIGO noise. We also discuss an approach to deal with nonstationary data by breaking longer data streams into smaller and locally stationary components.
An autoregressive repeatability animal model for test-day records in multiple lactations.
Carvalheira, J; Pollak, E J; Quaas, R L; Blake, R W
2002-08-01
Test-day (TD) models are becoming a standard for genetic evaluation of production traits in dairy cattle. Various approaches to model covariances between TD records include random regression, autoregressive repeatability, orthogonal polynomials, and models based on character processing. The applicability of these models is mainly associated with the number of parameters to estimate, incorporation of multiple lactations, and the accuracy of correlations generated by the cow's repeated expression of milking performance (TD yields) within and across lactations. We define and evaluate a multiple-lactation, autoregressive-repeatability model that disentangles environmental effects due to cow within and between lactations. Simulated records either included or ignored a long-term environmental effect between lactations. Our autoregressive TD animal model correctly detected presence and the absence of this effect and accurately recovered the assumed variance components and correlations underlying the data (10 parameters for three lactations). Estimates of variance components and autocorrelation coefficients were obtained using DFREML-simplex methodology. Given the value of this approach to reduce the size of residual variance components, autoregressive animal models are a preferable alternative to classical methods based on cumulative lactation yield to improve milk production in dairy cattle. PMID:12214997
Tropospheric Response to Estimated Spectrally Discriminated Solar Forcing Over the Past 500 Years
NASA Technical Reports Server (NTRS)
Rind, David; Hansen, James E. (Technical Monitor)
2000-01-01
The GISS Global Climate Middle Atmosphere Model (GCMAM) is used to investigate the effect of estimated solar irradiance changes on climate for the past 500 years. This model is employed to allow the impact of UV variations on the stratosphere to affect the troposphere via wave-mean flow interactions. Multiple experiments are done with only a total solar irradiance change (peaking at 0.2 percent from the Maunder Minimum to today); with estimated spectrally-varying irradiance changes (i.e., peak changes of 0.7 percent in the UV, 0.2 percent in the visible and near IR; and 0.07 percent in the IR greater than 1 micron); and the spectrally-varying changes in conjunction with model calculated ozone responses in the stratosphere. Results of the varying temperature patterns and radiation response will be discussed. Of interest is whether the different methods of forcing the solar-induced climate change produce different spatial surface temperature signatures, particularly ones that can be differentiated from greenhouse gas warming. In preliminary tests, spectrally-varying solar forcing with induced ozone changes for solar maximum minus solar minimum conditions results in a temperature signal that is primarily at high latitudes.The high latitude response arises due to solar/ozone-induced alterations in the stratospheric wind field that affect planetary wave propagation from the troposphere, and alter tropospheric advection patterns. In contrast, forcing by total solar irradiance changes produces significant response at low and subtropical latitudes as well, driven by water vapor and cloud feedbacks to the radiative perturbation.
NASA Technical Reports Server (NTRS)
Reginato, R. J.; Vedder, J. F.; Idso, S. B.; Jackson, R. D.; Blanchard, M. B.; Goettelman, R.
1977-01-01
For several days in March of 1975, reflected solar radiation measurements were obtained from smooth and rough surfaces of wet, drying, and continually dry Avondale loam at Phoenix, Arizona, with pyranometers located 50 cm above the ground surface and a multispectral scanner flown at a 300-m height. The simple summation of the different band radiances measured by the multispectral scanner proved equally as good as the pyranometer data for estimating surface soil water content if the multispectral scanner data were standardized with respect to the intensity of incoming solar radiation or the reflected radiance from a reference surface, such as the continually dry soil. Without this means of standardization, multispectral scanner data are most useful in a spectral band ratioing context. Our results indicated that, for the bands used, no significant information on soil water content could be obtained by band ratioing. Thus the variability in soil water content should insignificantly affect soil-type discrimination based on identification of type-specific spectral signatures. Therefore remote sensing, conducted in the 0.4- to 1.0-micron wavelength region of the solar spectrum, would seem to be much More suited to identifying crop and soil types than to estimating of soil water content.
NASA Technical Reports Server (NTRS)
Howell, Leonard W.
2002-01-01
The method of Maximum Likelihood (ML) is used to estimate the spectral parameters of an assumed broken power law energy spectrum from simulated detector responses. This methodology, which requires the complete specificity of all cosmic-ray detector design parameters, is shown to provide approximately unbiased, minimum variance, and normally distributed spectra information for events detected by an instrument having a wide range of commonly used detector response functions. The ML procedure, coupled with the simulated performance of a proposed space-based detector and its planned life cycle, has proved to be of significant value in the design phase of a new science instrument. The procedure helped make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope. This ML methodology is then generalized to estimate broken power law spectral parameters from real cosmic-ray data sets.
Geostatistical estimation of signal-to-noise ratios for spectral vegetation indices
NASA Astrophysics Data System (ADS)
Ji, Lei; Zhang, Li; Rover, Jennifer; Wylie, Bruce K.; Chen, Xuexia
2014-10-01
In the past 40 years, many spectral vegetation indices have been developed to quantify vegetation biophysical parameters. An ideal vegetation index should contain the maximum level of signal related to specific biophysical characteristics and the minimum level of noise such as background soil influences and atmospheric effects. However, accurate quantification of signal and noise in a vegetation index remains a challenge, because it requires a large number of field measurements or laboratory experiments. In this study, we applied a geostatistical method to estimate signal-to-noise ratio (S/N) for spectral vegetation indices. Based on the sample semivariogram of vegetation index images, we used the standardized noise to quantify the noise component of vegetation indices. In a case study in the grasslands and shrublands of the western United States, we demonstrated the geostatistical method for evaluating S/N for a series of soil-adjusted vegetation indices derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The soil-adjusted vegetation indices were found to have higher S/N values than the traditional normalized difference vegetation index (NDVI) and simple ratio (SR) in the sparsely vegetated areas. This study shows that the proposed geostatistical analysis can constitute an efficient technique for estimating signal and noise components in vegetation indices.
Geostatistical estimation of signal-to-noise ratios for spectral vegetation indices
Ji, Lei; Zhang, Li; Rover, Jennifer R.; Wylie, Bruce K.; Chen, Xuexia
2014-01-01
In the past 40 years, many spectral vegetation indices have been developed to quantify vegetation biophysical parameters. An ideal vegetation index should contain the maximum level of signal related to specific biophysical characteristics and the minimum level of noise such as background soil influences and atmospheric effects. However, accurate quantification of signal and noise in a vegetation index remains a challenge, because it requires a large number of field measurements or laboratory experiments. In this study, we applied a geostatistical method to estimate signal-to-noise ratio (S/N) for spectral vegetation indices. Based on the sample semivariogram of vegetation index images, we used the standardized noise to quantify the noise component of vegetation indices. In a case study in the grasslands and shrublands of the western United States, we demonstrated the geostatistical method for evaluating S/N for a series of soil-adjusted vegetation indices derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The soil-adjusted vegetation indices were found to have higher S/N values than the traditional normalized difference vegetation index (NDVI) and simple ratio (SR) in the sparsely vegetated areas. This study shows that the proposed geostatistical analysis can constitute an efficient technique for estimating signal and noise components in vegetation indices.
Estimating probability densities from short samples: A parametric maximum likelihood approach
NASA Astrophysics Data System (ADS)
Dudok de Wit, T.; Floriani, E.
1998-10-01
A parametric method similar to autoregressive spectral estimators is proposed to determine the probability density function (PDF) of a random set. The method proceeds by maximizing the likelihood of the PDF, yielding estimates that perform equally well in the tails as in the bulk of the distribution. It is therefore well suited for the analysis of short sets drawn from smooth PDF's and stands out by the simplicity of its computational scheme. Its advantages and limitations are discussed.
On the maximum-entropy/autoregressive modeling of time series
NASA Technical Reports Server (NTRS)
Chao, B. F.
1984-01-01
The autoregressive (AR) model of a random process is interpreted in the light of the Prony's relation which relates a complex conjugate pair of poles of the AR process in the z-plane (or the z domain) on the one hand, to the complex frequency of one complex harmonic function in the time domain on the other. Thus the AR model of a time series is one that models the time series as a linear combination of complex harmonic functions, which include pure sinusoids and real exponentials as special cases. An AR model is completely determined by its z-domain pole configuration. The maximum-entropy/autogressive (ME/AR) spectrum, defined on the unit circle of the z-plane (or the frequency domain), is nothing but a convenient, but ambiguous visual representation. It is asserted that the position and shape of a spectral peak is determined by the corresponding complex frequency, and the height of the spectral peak contains little information about the complex amplitude of the complex harmonic functions.
Estimating workload using EEG spectral power and ERPs in the n-back task
NASA Astrophysics Data System (ADS)
Brouwer, Anne-Marie; Hogervorst, Maarten A.; van Erp, Jan B. F.; Heffelaar, Tobias; Zimmerman, Patrick H.; Oostenveld, Robert
2012-08-01
Previous studies indicate that both electroencephalogram (EEG) spectral power (in particular the alpha and theta band) and event-related potentials (ERPs) (in particular the P300) can be used as a measure of mental work or memory load. We compare their ability to estimate workload level in a well-controlled task. In addition, we combine both types of measures in a single classification model to examine whether this results in higher classification accuracy than either one alone. Participants watched a sequence of visually presented letters and indicated whether or not the current letter was the same as the one (n instances) before. Workload was varied by varying n. We developed different classification models using ERP features, frequency power features or a combination (fusion). Training and testing of the models simulated an online workload estimation situation. All our ERP, power and fusion models provide classification accuracies between 80% and 90% when distinguishing between the highest and the lowest workload condition after 2 min. For 32 out of 35 participants, classification was significantly higher than chance level after 2.5 s (or one letter) as estimated by the fusion model. Differences between the models are rather small, though the fusion model performs better than the other models when only short data segments are available for estimating workload.
Eckhard, Timo; Valero, Eva M; Hernández-Andrés, Javier; Schnitzlein, Markus
2014-02-01
The performance of learning-based spectral estimation is greatly influenced by the set of training samples selected to create the reconstruction model. Training sample selection schemes can be categorized into global and local approaches. Most of the previously proposed global training schemes aim to reduce the number of training samples, or a selection of representative samples, to maintain the generality of the training dataset. This work relates to printed ink reflectance estimation for quality assessment in in-line print inspection. We propose what we believe is a novel global training scheme that models a large population of realistic printable ink reflectances. Based on this dataset, we used a recursive top-down algorithm to reject clusters of training samples that do not enhance the performance of a linear least-square regression (pseudoinverse-based estimation) process. A set of experiments with real camera response data of a 12-channel multispectral camera system illustrate the advantages of this selection scheme over some other state-of-the-art algorithms. For our data, our method of global training sample selection outperforms other methods in terms of estimation quality and, more importantly, can quickly handle large datasets. Furthermore, we show that reflectance modeling is a reasonable, convenient tool to generate large training sets for print inspection applications. PMID:24514188
Two-Microphone Noise Reduction Using Spatial Information-Based Spectral Amplitude Estimation
NASA Astrophysics Data System (ADS)
Li, Kai; Guo, Yanmeng; Fu, Qiang; Li, Junfeng; Yan, Yonghong
Traditional two-microphone noise reduction algorithms to deal with highly nonstationary directional noises generally use the direction of arrival or phase difference information. The performance of these algorithms deteriorate when diffuse noises coexist with nonstationary directional noises in realistic adverse environments. In this paper, we present a two-channel noise reduction algorithm using a spatial information-based speech estimator and a spatial-information-controlled soft-decision noise estimator to improve the noise reduction performance in realistic non-stationary noisy environments. A target presence probability estimator based on Bayes rules using both phase difference and magnitude squared coherence is proposed for soft-decision of noise estimation, so that they can share complementary advantages when both directional noises and diffuse noises are present. Performances of the proposed two-microphone noise reduction algorithm are evaluated by noise reduction, log-spectral distance (LSD) and word recognition rate (WRR) of a distant-talking ASR system in a real room's noisy environment. Experimental results show that the proposed algorithm achieves better noises suppression without further distorting the desired signal components over the comparative dual-channel noise reduction algorithms.
[Research on Oil Sands Spectral Characteristics and Oil Content by Remote Sensing Estimation].
You, Jin-feng; Xing, Li-xin; Pan, Jun; Shan, Xuan-long; Liang, Li-heng; Fan, Rui-xue
2015-04-01
Visible and near infrared spectroscopy is a proven technology to be widely used in identification and exploration of hydrocarbon energy sources with high spectral resolution for detail diagnostic absorption characteristics of hydrocarbon groups. The most prominent regions for hydrocarbon absorption bands are 1,740-1,780, 2,300-2,340 and 2,340-2,360 nm by the reflectance of oil sands samples. These spectral ranges are dominated by various C-H overlapping overtones and combination bands. Meanwhile, there is relatively weak even or no absorption characteristics in the region from 1,700 to 1,730 nm in the spectra of oil sands samples with low bitumen content. With the increase in oil content, in the spectral range of 1,700-1,730 nm the obvious hydrocarbon absorption begins to appear. The bitumen content is the critical parameter for oil sands reserves estimation. The absorption depth was used to depict the response intensity of the absorption bands controlled by first-order overtones and combinations of the various C-H stretching and bending fundamentals. According to the Pearson and partial correlation relationships of oil content and absorption depth dominated by hydrocarbon groups in 1,740-1,780, 2,300-2,340 and 2,340-2,360 nm wavelength range, the scheme of association mode was established between the intensity of spectral response and bitumen content, and then unary linear regression(ULR) and partial least squares regression (PLSR) methods were employed to model the equation between absorption depth attributed to various C-H bond and bitumen content. There were two calibration equations in which ULR method was employed to model the relationship between absorption depth near 2,350 nm region and bitumen content and PLSR method was developed to model the relationship between absorption depth of 1,758, 2,310, 2,350 nm regions and oil content. It turned out that the calibration models had good predictive ability and high robustness and they could provide the scientific
Estimation of glottal source features from the spectral envelope of the acoustic speech signal
NASA Astrophysics Data System (ADS)
Torres, Juan Felix
Speech communication encompasses diverse types of information, including phonetics, affective state, voice quality, and speaker identity. From a speech production standpoint, the acoustic speech signal can be mainly divided into glottal source and vocal tract components, which play distinct roles in rendering the various types of information it contains. Most deployed speech analysis systems, however, do not explicitly represent these two components as distinct entities, as their joint estimation from the acoustic speech signal becomes an ill-defined blind deconvolution problem. Nevertheless, because of the desire to understand glottal behavior and how it relates to perceived voice quality, there has been continued interest in explicitly estimating the glottal component of the speech signal. To this end, several inverse filtering (IF) algorithms have been proposed, but they are unreliable in practice because of the blind formulation of the separation problem. In an effort to develop a method that can bypass the challenging IF process, this thesis proposes a new glottal source information extraction method that relies on supervised machine learning to transform smoothed spectral representations of speech, which are already used in some of the most widely deployed and successful speech analysis applications, into a set of glottal source features. A transformation method based on Gaussian mixture regression (GMR) is presented and compared to current IF methods in terms of feature similarity, reliability, and speaker discrimination capability on a large speech corpus, and potential representations of the spectral envelope of speech are investigated for their ability represent glottal source variation in a predictable manner. The proposed system was found to produce glottal source features that reasonably matched their IF counterparts in many cases, while being less susceptible to spurious errors. The development of the proposed method entailed a study into the aspects
NASA Astrophysics Data System (ADS)
Smith-Boughner, Lindsay
Many Earth systems cannot be studied directly. One cannot measure the velocities of convecting fluid in the Earth's core but can measure the magnetic field generated by these motions on the surface. Examining how the magnetic field changes over long periods of time, using power spectral density estimation provides insight into the dynamics driving the system. The changes in the magnetic field can also be used to study Earth properties - variations in magnetic fields outside of Earth like the ring-current induce currents to flow in the Earth, generating magnetic fields. Estimating the transfer function between the external changes and the induced response characterizes the electromagnetic response of the Earth. From this response inferences can be made about the electrical conductivity of the Earth. However, these types of time series, and many others have long breaks in the record with no samples available and limit the analysis. Standard methods require interpolation or section averaging, with associated problems of introducing bias or reducing the frequency resolution. Extending the methods of Fodor and Stark (2000), who adapt a set of orthogonal multi-tapers to compensate for breaks in sampling- an algorithm and software package for applying these techniques is developed. Methods of empirically estimating the average transfer function of a set of tapers and confidence intervals are also tested. These methods are extended for cross-spectral, coherence and transfer function estimation in the presence of noise. With these methods, new analysis of a highly interrupted ocean sediment core from the Oligocene (Hartl et al., 1993) reveals a quasi-periodic signal in the calibrated paleointensity time series at 2.5 cpMy. The power in the magnetic field during this period appears to be dominated by reversal rate processes with less overall power than the early Oligocene. Previous analysis of the early Oligocene by Constable et al. (1998) detected a signal near 8 cp
NASA Astrophysics Data System (ADS)
Shi, C.; Wang, L.
2015-12-01
Spectral unmixing is the process of decomposing the measured spectrum of a mixed pixel into a set of pure spectral signatures called endmembers and their corresponding abundances indicating the fractional area coverage of each endmember present in the pixel. A substantial number of spectral unmixing studies rely on a spectral mixture model which assumes that spectral mixing only occurs within the extent of a pixel. However, due to adjacency effect, the spectral measurement of the pixel may be contaminated by spatial interactions from materials that are present in its spatial neighborhood. In this paper, a linear spatial spectral mixture model is developed to improve the accuracy of the estimated abundance of invasive saltcedar along the Forgotten River reach of the Rio Grande. A spatial weights matrix which specifies for each pixel the locations and the weights of its neighborhood set is used to summarize the spatial relationships among pixels in the Landsat data. A spatial lag operator, defined as a weighted average of the values at neighboring locations, is adopted as an expression of spectral contribution from nearby pixels and added to the classic linear mixture model. The fractional abundances are iteratively estimated using the alternating direction method of multipliers (ADMM) algorithm. With the incorporation of adjacency effect, RMSEs of the fractional cover of ground classes were reduced. The derived sub-pixel abundances of saltcedar are beneficial for ecological management.
Estimation of the spectral parameter kappa in the region of the Gulf of California, Mexico
NASA Astrophysics Data System (ADS)
Castro, Raúl R.; Ávila-Barrientos, Lenin
2015-10-01
We analyzed records from the Broadband Seismological Network of the Gulf of California (RESBAN) and from stations of the NARS-Baja array, operated by CICESE, Ensenada, Baja California, Mexico, to make estimates of the spectral decay parameter kappa ( κ). This attenuation parameter is important for evaluating the seismic risk and hazard of this region. Thirteen shallow earthquakes with focal depths less than 20 km and magnitudes between 5.1 and 6.6 were selected to calculate κ and the near-site attenuation κ 0. We used three different approaches to estimate κ 0: (a) with individual measurements of κ from vector modulus of three-component spectral amplitudes at different epicentral distances and extrapolating to zero distance to estimate κ 0, (b) with individual measurements using vertical component spectra, and (c) measuring from the high-frequency part of the site transfer function determined calculating the horizontal-to-vertical spectral ratio (HVSR) method. For most stations, the three methods give similar results. At short distances (50-60 km), κ takes values close to 0.04 s at NE76, the station located in the middle of the array. κ increases with distance taking an average value of up to 0.18 s for distances close to 500 km. κ 0 at most sites is close to 0.03 s, except for GUYB (Guaymas) that has a κ 0 = 0.05 s and NE83 (Navolato) with κ 0 = 0.065 s, both stations located in the continent, on the eastern side of the gulf, where the soils are less consolidated. Finally, we analyze if κ 0 correlates with magnitude and back azimuth, and we found that for most stations, κ 0 does not correlate with either one. However, station TOPB, located on basalt, shows a moderate correlation with magnitude, with κ 0 increasing with increasing M W in a short back-azimuth range. We also found that for station NE74, located on soft soil, κ 0 correlates with back azimuth, having lower values for azimuths near 120°.
ENSO Prediction using Vector Autoregressive Models
NASA Astrophysics Data System (ADS)
Chapman, D. R.; Cane, M. A.; Henderson, N.; Lee, D.; Chen, C.
2013-12-01
A recent comparison (Barnston et al, 2012 BAMS) shows the ENSO forecasting skill of dynamical models now exceeds that of statistical models, but the best statistical models are comparable to all but the very best dynamical models. In this comparison the leading statistical model is the one based on the Empirical Model Reduction (EMR) method. Here we report on experiments with multilevel Vector Autoregressive models using only sea surface temperatures (SSTs) as predictors. VAR(L) models generalizes Linear Inverse Models (LIM), which are a VAR(1) method, as well as multilevel univariate autoregressive models. Optimal forecast skill is achieved using 12 to 14 months of prior state information (i.e 12-14 levels), which allows SSTs alone to capture the effects of other variables such as heat content as well as seasonality. The use of multiple levels allows the model advancing one month at a time to perform at least as well for a 6 month forecast as a model constructed to explicitly forecast 6 months ahead. We infer that the multilevel model has fully captured the linear dynamics (cf. Penland and Magorian, 1993 J. Climate). Finally, while VAR(L) is equivalent to L-level EMR, we show in a 150 year cross validated assessment that we can increase forecast skill by improving on the EMR initialization procedure. The greatest benefit of this change is in allowing the prediction to make effective use of information over many more months.
A Recommended Procedure for Estimating the Cosmic Ray Spectral Parameter of a Simple Power Law
NASA Technical Reports Server (NTRS)
Howell, Leonard W.; Rose, M. Franklin (Technical Monitor)
2000-01-01
A simple power law model consisting of a single spectral index a(f(sub i)) is believed to be an adequate description of the galactic cosmic ray (GQ proton flux at energies below 1013 eV. Two procedures for estimating a(f(sub i)), referred as (1) the method of moments, and (2) maximum likelihood, are developed and their statistical performance compared. I concluded that the maximum likelihood procedure attains the most desirable statistical properties and is hence the recommended statistic estimation procedure for estimating a1. The maximum likelihood procedure is then generalized for application to a set of real cosmic ray data and thereby makes this approach applicable to existing cosmic ray data sets. Several other important results, such as the relationship between collecting power and detector energy resolution, as well as inclusion of a non-Gaussian detector response function, are presented. These results have many practical benefits in the design phase of a cosmic ray detector because they permit instrument developers to make important trade studies in design parameters as a function of one of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose practical limits to the design envelope.
NASA Astrophysics Data System (ADS)
Cui, Qian; Shi, Jiancheng; Xu, Yuanliu
2011-12-01
Water is the basic needs for human society, and the determining factor of stability of ecosystem as well. There are lots of lakes on Tibet Plateau, which will lead to flood and mudslide when the water expands sharply. At present, water area is extracted from TM or SPOT data for their high spatial resolution; however, their temporal resolution is insufficient. MODIS data have high temporal resolution and broad coverage. So it is valuable resource for detecting the change of water area. Because of its low spatial resolution, mixed-pixels are common. In this paper, four spectral libraries are built using MOD09A1 product, based on that, water body is extracted in sub-pixels utilizing Multiple Endmembers Spectral Mixture Analysis (MESMA) using MODIS daily reflectance data MOD09GA. The unmixed result is comparing with contemporaneous TM data and it is proved that this method has high accuracy.
Super-resolution spectral estimation in short-time non-contact vital sign measurement
NASA Astrophysics Data System (ADS)
Sun, Li; Li, Yusheng; Hong, Hong; Xi, Feng; Cai, Weidong; Zhu, Xiaohua
2015-04-01
Non-contact techniques for measuring vital signs attract great interest due to the benefits shown in medical monitoring, military application, etc. However, the presence of respiration harmonics caused by nonlinear phase modulation will result in performance degradation. Suffering from smearing and leakage problems, conventional discrete Fourier transform (DFT) based methods cannot distinguish the heartbeat component from closely located respiration harmonics in frequency domain, especially in short-time processing. In this paper, the theory of sparse reconstruction is merged with an extended harmonic model of vital signals, aiming at achieving a super-resolution spectral estimation of vital signals by additionally exploiting the inherent sparse prior information. Both simulated and experimental results show that the proposed algorithm has superior performance to DFT-based methods and the recently applied multiple signal classification algorithm, and the required processing window length has been shortened to 5.12 s.
Super-resolution spectral estimation in short-time non-contact vital sign measurement.
Sun, Li; Li, Yusheng; Hong, Hong; Xi, Feng; Cai, Weidong; Zhu, Xiaohua
2015-04-01
Non-contact techniques for measuring vital signs attract great interest due to the benefits shown in medical monitoring, military application, etc. However, the presence of respiration harmonics caused by nonlinear phase modulation will result in performance degradation. Suffering from smearing and leakage problems, conventional discrete Fourier transform (DFT) based methods cannot distinguish the heartbeat component from closely located respiration harmonics in frequency domain, especially in short-time processing. In this paper, the theory of sparse reconstruction is merged with an extended harmonic model of vital signals, aiming at achieving a super-resolution spectral estimation of vital signals by additionally exploiting the inherent sparse prior information. Both simulated and experimental results show that the proposed algorithm has superior performance to DFT-based methods and the recently applied multiple signal classification algorithm, and the required processing window length has been shortened to 5.12 s. PMID:25933881
A robust approach to estimate stellar angular diameters from photometry and spectral type
NASA Astrophysics Data System (ADS)
Chelli, A.; Bourges, L.; Duvert, G.; Lafrasse, S.; Mella, G.; Le Bouquin, J.-B.; Chesneau, O.
2014-07-01
Observing reference stars with a known diameter is almost the only possibility to calibrate optical interferometry observations. The JMMC Calibrator Workgroup develops methods to ascertain the angular diameter of stars since 2000 and provides this expertise in the SearchCal software and associated databases. We provide on a regularly basis the JSDC, a catalogue of such stars, and an open access to our server that dynamically finds calibrators near science objects by querying CDS hosted catalogs. Here we propose a novel approach in the estimation of angular stellar diameters based on observational quantities only. It bypasses the knowledge of the visual extinction and intrinsic colors, thanks to the use of absorption free pseudo-colors (AFC) and the spectral type number on the x-axis. This new methodology allows to compute the angular diameter of 443 703 stars with a relative precision of about 1%. This calibrator set will become after filtering the next JSDC release.
Spectral parameter estimation of CAT radar echoes in the presence of fading clutter
NASA Technical Reports Server (NTRS)
Sato, T.; Woodman, R. F.
1980-01-01
The analysis technique and a part of the results obtained from CAT radar echoes from higher troposphere and lower stratosphere are presented. First, the effect of processing distortion caused by the periodogram method using FFT algorithm on the slowly fading ground clutter echo is discussed. It is shown that an extremely narrow clutter spectrum can spill over the entire frequency range if the data are truncated at a tie sorter than their correlation time affecting largely the estimation of the CAT spectrum contribution, especially when the latter is a few tens of dB weaker than the former. A nonlinear least squares fitting procedure is used to parameterize the observed power spectrum in terms of CAT echo power, Doppler shift, spectral width, and the parameters which specify the shape of the clutter component.
NASA Technical Reports Server (NTRS)
Eren, K.
1980-01-01
The mathematical background in spectral analysis as applied to geodetic applications is summarized. The resolution (cut-off frequency) of the GEOS 3 altimeter data is examined by determining the shortest wavelength (corresponding to the cut-off frequency) recoverable. The data from some 18 profiles are used. The total power (variance) in the sea surface topography with respect to the reference ellipsoid as well as with respect to the GEM-9 surface is computed. A fast inversion algorithm for matrices of simple and block Toeplitz matrices and its application to least squares collocation is explained. This algorithm yields a considerable gain in computer time and storage in comparison with conventional least squares collocation. Frequency domain least squares collocation techniques are also introduced and applied to estimating gravity anomalies from GEOS 3 altimeter data. These techniques substantially reduce the computer time and requirements in storage associated with the conventional least squares collocation. Numerical examples given demonstrate the efficiency and speed of these techniques.
NASA Technical Reports Server (NTRS)
Bergstrom, Robert W.; Pilewskie, Peter; Schmid, Beat; Russell, Philip B.
2003-01-01
Using measurements of the spectral solar radiative flux and optical depth for 2 days (24 August and 6 September 2000) during the SAFARI 2000 intensive field experiment and a detailed radiative transfer model, we estimate the spectral single scattering albedo of the aerosol layer. The single scattering albedo is similar on the 2 days even though the optical depth for the aerosol layer was quite different. The aerosol single scattering albedo was between 0.85 and 0.90 at 350 nm, decreasing to 0.6 in the near infrared. The magnitude and decrease with wavelength of the single scattering albedo are consistent with the absorption properties of small black carbon particles. We estimate the uncertainty in the single scattering albedo due to the uncertainty in the measured fractional absorption and optical depths. The uncertainty in the single scattering albedo is significantly less on the high-optical-depth day (6 September) than on the low-optical-depth day (24 August). On the high-optical-depth day, the uncertainty in the single scattering albedo is 0.02 in the midvisible whereas on the low-optical-depth day the uncertainty is 0.08 in the midvisible. On both days, the uncertainty becomes larger in the near infrared. We compute the radiative effect of the aerosol by comparing calculations with and without the aerosol. The effect at the top of the atmosphere (TOA) is to cool the atmosphere by 13 W/sq m on 24 August and 17 W/sq m on 6 September. The effect on the downward flux at the surface is a reduction of 57 W/sq m on 24 August and 200 W/sq m on 6 September. The aerosol effect on the downward flux at the surface is in good agreement with the results reported from the Indian Ocean Experiment (INDOEX).
Rainfall Estimation over the Nile Basin using Multi-Spectral, Multi- Instrument Satellite Techniques
NASA Astrophysics Data System (ADS)
Habib, E.; Kuligowski, R.; Sazib, N.; Elshamy, M.; Amin, D.; Ahmed, M.
2012-04-01
Management of Egypt's Aswan High Dam is critical not only for flood control on the Nile but also for ensuring adequate water supplies for most of Egypt since rainfall is scarce over the vast majority of its land area. However, reservoir inflow is driven by rainfall over Sudan, Ethiopia, Uganda, and several other countries from which routine rain gauge data are sparse. Satellite- derived estimates of rainfall offer a much more detailed and timely set of data to form a basis for decisions on the operation of the dam. A single-channel infrared (IR) algorithm is currently in operational use at the Egyptian Nile Forecast Center (NFC). In this study, the authors report on the adaptation of a multi-spectral, multi-instrument satellite rainfall estimation algorithm (Self- Calibrating Multivariate Precipitation Retrieval, SCaMPR) for operational application by NFC over the Nile Basin. The algorithm uses a set of rainfall predictors that come from multi-spectral Infrared cloud top observations and self-calibrate them to a set of predictands that come from the more accurate, but less frequent, Microwave (MW) rain rate estimates. For application over the Nile Basin, the SCaMPR algorithm uses multiple satellite IR channels that have become recently available to NFC from the Spinning Enhanced Visible and Infrared Imager (SEVIRI). Microwave rain rates are acquired from multiple sources such as the Special Sensor Microwave/Imager (SSM/I), the Special Sensor Microwave Imager and Sounder (SSMIS), the Advanced Microwave Sounding Unit (AMSU), the Advanced Microwave Scanning Radiometer on EOS (AMSR-E), and the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm has two main steps: rain/no-rain separation using discriminant analysis, and rain rate estimation using stepwise linear regression. We test two modes of algorithm calibration: real- time calibration with continuous updates of coefficients with newly coming MW rain rates, and calibration using static
NASA Astrophysics Data System (ADS)
Krezhova, Dora D.; Yanev, Tony K.
2007-04-01
Results from a remote sensing study of the influence of stress factors on the leaf spectral reflectance of wheat and tomato plants contaminated by viruses and pea plants treated with herbicides are presented and discussed. The changes arising in the spectral reflectance characteristics of control and treated plants are estimated through statistical methods as well as through derivative analysis to determine specific reflectance features in the red edge region.
Spectral Feature Analysis for Quantitative Estimation of Cyanobacteria Chlorophyll-A
NASA Astrophysics Data System (ADS)
Lin, Yi; Ye, Zhanglin; Zhang, Yugan; Yu, Jie
2016-06-01
In recent years, lake eutrophication caused a large of Cyanobacteria bloom which not only brought serious ecological disaster but also restricted the sustainable development of regional economy in our country. Chlorophyll-a is a very important environmental factor to monitor water quality, especially for lake eutrophication. Remote sensed technique has been widely utilized in estimating the concentration of chlorophyll-a by different kind of vegetation indices and monitoring its distribution in lakes, rivers or along coastline. For each vegetation index, its quantitative estimation accuracy for different satellite data might change since there might be a discrepancy of spectral resolution and channel center between different satellites. The purpose this paper is to analyze the spectral feature of chlorophyll-a with hyperspectral data (totally 651 bands) and use the result to choose the optimal band combination for different satellites. The analysis method developed here in this study could be useful to recognize and monitor cyanobacteria bloom automatically and accrately. In our experiment, the reflectance (from 350nm to 1000nm) of wild cyanobacteria in different consistency (from 0 to 1362.11ug/L) and the corresponding chlorophyll-a concentration were measured simultaneously. Two kinds of hyperspectral vegetation indices were applied in this study: simple ratio (SR) and narrow band normalized difference vegetation index (NDVI), both of which consists of any two bands in the entire 651 narrow bands. Then multivariate statistical analysis was used to construct the linear, power and exponential models. After analyzing the correlation between chlorophyll-a and single band reflectance, SR, NDVI respetively, the optimal spectral index for quantitative estimation of cyanobacteria chlorophyll-a, as well corresponding central wavelength and band width were extracted. Results show that: Under the condition of water disturbance, SR and NDVI are both suitable for quantitative
Garde, Ainara; Karlen, Walter; Ansermino, J Mark; Dumont, Guy A
2014-01-01
The photoplethysmogram (PPG) obtained from pulse oximetry measures local variations of blood volume in tissues, reflecting the peripheral pulse modulated by heart activity, respiration and other physiological effects. We propose an algorithm based on the correntropy spectral density (CSD) as a novel way to estimate respiratory rate (RR) and heart rate (HR) from the PPG. Time-varying CSD, a technique particularly well-suited for modulated signal patterns, is applied to the PPG. The respiratory and cardiac frequency peaks detected at extended respiratory (8 to 60 breaths/min) and cardiac (30 to 180 beats/min) frequency bands provide RR and HR estimations. The CSD-based algorithm was tested against the Capnobase benchmark dataset, a dataset from 42 subjects containing PPG and capnometric signals and expert labeled reference RR and HR. The RR and HR estimation accuracy was assessed using the unnormalized root mean square (RMS) error. We investigated two window sizes (60 and 120 s) on the Capnobase calibration dataset to explore the time resolution of the CSD-based algorithm. A longer window decreases the RR error, for 120-s windows, the median RMS error (quartiles) obtained for RR was 0.95 (0.27, 6.20) breaths/min and for HR was 0.76 (0.34, 1.45) beats/min. Our experiments show that in addition to a high degree of accuracy and robustness, the CSD facilitates simultaneous and efficient estimation of RR and HR. Providing RR every minute, expands the functionality of pulse oximeters and provides additional diagnostic power to this non-invasive monitoring tool. PMID:24466088
Zhu, Zai-Chun; Chen, Lian-Qun; Zhang, Jin-Shui; Pan, Yao-Zhong; Zhu, Wen-Quan; Hu, Tan-Gao
2012-07-01
Crop yield estimation division is the basis of crop yield estimation; it provides an important scientific basis for estimation research and practice. In the paper, MODIS EVI time-series data during winter wheat growth period is selected as the division data; JiangSu province is study area; A division method combined of advanced spectral angle mapping(SVM) and K-means clustering is presented, and tested in winter wheat yield estimation by remote sensing. The results show that: division method of spectral angle clustering can take full advantage of crop growth process that is reflected by MODIS time series data, and can fully reflect region differences of winter wheat that is brought by climate difference. Compared with the traditional division method, yield estimation result based on division result of spectral angle clustering has higher R2 (0.702 6 than 0.624 8) and lower RMSE (343.34 than 381.34 kg x hm(-2)), reflecting the advantages of the new division method in the winter wheat yield estimation. The division method in the paper only use convenient obtaining time-series remote sensing data of low-resolution as division data, can divide winter wheat into similar and well characterized region, accuracy and stability of yield estimation model is also very good, which provides an efficient way for winter wheat estimation by remote sensing, and is conducive to winter wheat yield estimation. PMID:23016349
Kepler AutoRegressive Planet Search
NASA Astrophysics Data System (ADS)
Caceres, Gabriel Antonio; Feigelson, Eric
2016-01-01
The Kepler AutoRegressive Planet Search (KARPS) project uses statistical methodology associated with autoregressive (AR) processes to model Kepler lightcurves in order to improve exoplanet transit detection in systems with high stellar variability. We also introduce a planet-search algorithm to detect transits in time-series residuals after application of the AR models. One of the main obstacles in detecting faint planetary transits is the intrinsic stellar variability of the host star. The variability displayed by many stars may have autoregressive properties, wherein later flux values are correlated with previous ones in some manner. Our analysis procedure consisting of three steps: pre-processing of the data to remove discontinuities, gaps and outliers; AR-type model selection and fitting; and transit signal search of the residuals using a new Transit Comb Filter (TCF) that replaces traditional box-finding algorithms. The analysis procedures of the project are applied to a portion of the publicly available Kepler light curve data for the full 4-year mission duration. Tests of the methods have been made on a subset of Kepler Objects of Interest (KOI) systems, classified both as planetary `candidates' and `false positives' by the Kepler Team, as well as a random sample of unclassified systems. We find that the ARMA-type modeling successfully reduces the stellar variability, by a factor of 10 or more in active stars and by smaller factors in more quiescent stars. A typical quiescent Kepler star has an interquartile range (IQR) of ~10 e-/sec, which may improve slightly after modeling, while those with IQR ranging from 20 to 50 e-/sec, have improvements from 20% up to 70%. High activity stars (IQR exceeding 100) markedly improve. A periodogram based on the TCF is constructed to concentrate the signal of these periodic spikes. When a periodic transit is found, the model is displayed on a standard period-folded averaged light curve. Our findings to date on real
Power spectral density of velocity fluctuations estimated from phase Doppler data
NASA Astrophysics Data System (ADS)
Jedelsky, Jan; Lizal, Frantisek; Jicha, Miroslav
2012-04-01
Laser Doppler Anemometry (LDA) and its modifications such as PhaseDoppler Particle Anemometry (P/DPA) is point-wise method for optical nonintrusive measurement of particle velocity with high data rate. Conversion of the LDA velocity data from temporal to frequency domain - calculation of power spectral density (PSD) of velocity fluctuations, is a non trivial task due to nonequidistant data sampling in time. We briefly discuss possibilities for the PSD estimation and specify limitations caused by seeding density and other factors of the flow and LDA setup. Arbitrary results of LDA measurements are compared with corresponding Hot Wire Anemometry (HWA) data in the frequency domain. Slot correlation (SC) method implemented in software program Kern by Nobach (2006) is used for the PSD estimation. Influence of several input parameters on resulting PSDs is described. Optimum setup of the software for our data of particle-laden air flow in realistic human airway model is documented. Typical character of the flow is described using PSD plots of velocity fluctuations with comments on specific properties of the flow. Some recommendations for improvements of future experiments to acquire better PSD results are given.
NASA Astrophysics Data System (ADS)
Trizna, D. B.; McNeal, G. D.
1985-12-01
This work is the second paper in a series of studies of the application of spectral estimations techniques to Doppler processing of coherent radar signals. In this work, simulated high-frequency (HF) radar sea scatter time series are generated and processed by use of three different spectral estimation algorithms and the fast Fourier transform (FFT). The sea clutter is simulated by narrowband filtering a wideband Gaussian noise spectrum in the frequency domain, with filter widths appropriate to describe first-order Bragg lines and second-order continuum. Targets are introduced as sinusoids, stepped by 5 dB for eight different echo power values, and stepped in Doppler frequency for four different values relative to the clutter. These simulations identify problems that appear unique to Doppler processing of coherent radar data in the presence of broadband clutter, and are in distinction to the application of spectral estimation to processing in the spatial domain. In the latter case, the spectral contributions are generally narrowly confined in the angular power spectral estimate, and the aim is to separate these contributions in the presence of noise. The HF radar application is concerned with separation of weak targets in the presence of stronger clutter returns, which are relati vely broad compared to the target return. It appears that the Burg maximum entropy method allows the detection of targets in clutter under conditions which the FFT is incapable of detection with any degree of accuracy.
NASA Astrophysics Data System (ADS)
Ward, A. L.; Draper, K.; Hasan, N.
2010-12-01
Knowledge of spatially variable aquifer hydraulic and sorption parameters is a pre-requisite for an improved understanding of the transport and spreading of sorbing solutes and for the development of effective strategies for remediation. Local-scale estimates of these parameters are often derived from core measurements but are typically not representative of field values. Fields-scale estimates are typically derived from pump and tracer tests but often lack the spatial resolution necessary to deconvolve the effects of fine-scale heterogeneities. Geophysical methods have the potential to bridge this gap both in terms of coverage and resolution, provided meaningful petrophysical relationships can be developed. The objective of this study was to develop a petrophysical relationship between soil textural attributes and the gamma-energy response of natural sediments. Measurements from Hanford’s 300 Area show the best model to be a linear relationship between 232Th concentration and clay content (R2 = 94%). This relationship was used to generate a 3-D distribution of clay mass fraction based on borehole spectral gamma logs. The distribution of clay was then used to predict distributions of permeability, porosity, bubbling pressure, and the pore-size distribution index, all of which are required for predicting variably saturated flow, as well as the specific surface area and cation exchange capacity needed for reactive transport predictions. With this approach, it is possible to obtain reliable estimates of hydraulic properties in zones that could not be characterized by field or laboratory measurements. The spatial distribution of flow properties is consistent with lithologic transitions inferred from geologist’s logs. A preferential flow path, identified from solute and heat tracer experiments and attributed to an erosional incision in the low-permeability Ringold Formation, is also evident. The resulting distributions can be used as a starting model for the
NASA Astrophysics Data System (ADS)
Dybus, W.; Benoit, M. H.; Ebinger, C. J.
2011-12-01
The crustal thickness beneath much of the eastern half of the US is largely unconstrained. Though there have been several controlled source seismic surveys of the region, many of these studies suffer from rays that turn in the crust above the Moho, resulting in somewhat ambiguous crustal thickness values. Furthermore, the broadband seismic station coverage east of the Mississippi has been limited, and most of the region remains largely understudied. In this study, we estimated the depth to the Moho using both spectral analysis and inversion of Bouguer gravity anomalies. We systematically estimated depths to lithospheric density contrasts from radial power spectra of Bouguer gravity within 100 km X 100 km windows eastward from the Mississippi River to the Atlantic Coast, and northward from North Carolina to Maine. The slopes and slope breaks in the radial power spectra were computed using an automated algorithm. The slope values for each window were visually inspected and then used to estimate the depth to the Moho and other lithospheric density contrasts beneath each windowed region. Additionally, we performed a standard Oldenburg-Parker inversion for lithospheric density contrasts using various reference depths and density contrasts that are realistic for the different physiographic provinces in the Eastern US. Our preliminary results suggest that the gravity-derived Moho depths are similar to those found using seismic data, and that the crust is relatively thinner (~28-33 km) than expected in beneath the Piedmont region (~35-40 km). Given the relative paucity of seismic data in the eastern US, analysis of onshore gravity data is a valuable tool for interpolating between seismic stations.
NASA Astrophysics Data System (ADS)
Qu, Liqin; Civco, Daniel; Lei, Tingwu; Yang, Xiusheng
2014-06-01
The dynamic sediment distribution in large rivers with dams constructed has often been the focus of considerable attention because of their potential adverse environmental impacts. Sedimentation modeling and environmental assessment of man-made projects are often hindered by the lack of sediment measurements with spatial details. This study aimed to investigate the method used to estimate the suspended sediment concentrations (SSCs) from on-site spectral measurements. The study investigated the spectral signature of river water from the natural channel and Sanmenxia Reservoir on the Yellow River. A field spectral survey was conducted through on-site spectral measurements by using a spectroradimeter and SSC estimation by sampling. Reectance at 750 nm to 950 nm, with all correlation coefficient (r) between SSC and reectance > 0:7, seemed to be the appropriate range for SSC estimation. Simulated Landsat Enhanced Thematic Mapper Plus Band 4 (760 nm to 900 nm) was used to build the single band model for estimating SSC. The results confirmed that the exponential model based on the relationship between SSC and reectance (R2 = 0:92, root mean square error [RMSE]= 0:241 g=l) was better than the linear model between reectance and logarithm-transformed SSC (R2 = 0:90, RMSE = 0:310 g=l). We also applied the Spectral Mixing Algorithm (SMA) from the tank experiment to the on-site spectral measurements. The result showed that the SMA models performed as well as the single band exponential model (R2 = 0:86, RMSE = 0:280 g=l). However, the valid range for application was improved from 1:99 g=l to 347 g=l. This study could provide critical instructional assistance for estimating SSC directly from remote sensing data.
NASA Astrophysics Data System (ADS)
Hansen, V.
1984-05-01
The distribution of direct and scattered solar radiant energy in the UV (295-385-nm), blue (385-495-nm), green-orange (495-630-nm), red (630-695-nm), and IR (695-2800-nm) bands commonly used in precision spectral pyranometers is estimated for clear sky conditions as a function of solar height, using a plane-parallel atmosphere model and data on the seasonal variation of the UV component at latitude 59.7 deg N. Integrated daily fluxes are also calculated for selected days of the year, and the results are compared with experimental measurements in graphs and tables. The model is found to give reasonable agreement with the observations, but fails to account for a significant blue shift with increasing solar height at heights above 15 deg. The measured distribution for March is given as UV 4.2, blue 8.9, green-orange 21.8, red 12.4, and IR 52.9 percent; for July, the respective values are 4.6, 16.1, 18.0, 10.7, and 50.6 percent.
Applications methods of spectral ratios in the estimation of site effects: Case Damien (Haiti)
NASA Astrophysics Data System (ADS)
Jean, B. J.; ST Fleur, S.
2014-12-01
Measurements of H/V type were carried out on the Damien site with Tromino hardware an « all in one » station which includes both the sensor and the integrated digitizer. A total of 32 measurements of seismic noise have been completed on this site in order to see if lithological site effects are detectable with this H/V method. After checking the H/V curve reliability criteria (length of the window to be analyzed, the number of windows analyzed, standard deviation) and the criteria for clear peaks in H/V (conditions for the amplitude, conditions for stability) found in the SESAME project in 2004, the results of the H/V spectra obtained are generally very consistent and clearly indicate site effects with peak resonance frequencies between 3 and 14 Hz. The presence of these well defined frequency peaks in the H/V spectral ratio indicates that the ground motion can be amplified by geomorphological site effects. Comparative analyzes of these H/V measurements with Grilla and Geopsy software were made in this paper to estimate the amplification magnitude of these effects. Graphical comparisons between the Grilla and Geopsy H/V maps were completed in this study and allow us to identify typical areas and their associated fundamental resonance frequencies.
NASA Astrophysics Data System (ADS)
Chang, Gao-Wei; Jian, Hong-Da; Yeh, Zong-Mu; Cheng, Chin-Pao
2004-02-01
For estimating spectral responsivities of digital video cameras, a filter-based optical system is designed with sophisticated filter selections, in this paper. The filter consideration in the presence of noise is central to the optical systems design, since the spectral filters primarily prescribe the structure of the perturbed system. A theoretical basis is presented to confirm that sophisticated filter selections can make this system as insensitive to noise as possible. Also, we propose a filter selection method based on the orthogonal-triangular (QR) decomposition with column pivoting (QRCP). To investigate the noise effects, we assess the estimation errors between the actual and estimated spectral responsivities, with the different signal-to-noise ratio (SNR) levels of an eight-bit/channel camera. Simulation results indicate that the proposed method yields satisfactory estimation accuracy. That is, the filter-based optical system with the spectral filters selected from the QRCP-based method is much less sensitive to noise than those with other filters from different selections.
No-reference image sharpness assessment in autoregressive parameter space.
Gu, Ke; Zhai, Guangtao; Lin, Weisi; Yang, Xiaokang; Zhang, Wenjun
2015-10-01
In this paper, we propose a new no-reference (NR)/blind sharpness metric in the autoregressive (AR) parameter space. Our model is established via the analysis of AR model parameters, first calculating the energy- and contrast-differences in the locally estimated AR coefficients in a pointwise way, and then quantifying the image sharpness with percentile pooling to predict the overall score. In addition to the luminance domain, we further consider the inevitable effect of color information on visual perception to sharpness and thereby extend the above model to the widely used YIQ color space. Validation of our technique is conducted on the subsets with blurring artifacts from four large-scale image databases (LIVE, TID2008, CSIQ, and TID2013). Experimental results confirm the superiority and efficiency of our method over existing NR algorithms, the stateof-the-art blind sharpness/blurriness estimators, and classical full-reference quality evaluators. Furthermore, the proposed metric can be also extended to stereoscopic images based on binocular rivalry, and attains remarkably high performance on LIVE3D-I and LIVE3D-II databases. PMID:26054063
Image restoration using 2D autoregressive texture model and structure curve construction
NASA Astrophysics Data System (ADS)
Voronin, V. V.; Marchuk, V. I.; Petrosov, S. P.; Svirin, I.; Agaian, S.; Egiazarian, K.
2015-05-01
In this paper an image inpainting approach based on the construction of a composite curve for the restoration of the edges of objects in an image using the concepts of parametric and geometric continuity is presented. It is shown that this approach allows to restore the curved edges and provide more flexibility for curve design in damaged image by interpolating the boundaries of objects by cubic splines. After edge restoration stage, a texture restoration using 2D autoregressive texture model is carried out. The image intensity is locally modeled by a first spatial autoregressive model with support in a strongly causal prediction region on the plane. Model parameters are estimated by Yule-Walker method. Several examples considered in this paper show the effectiveness of the proposed approach for large objects removal as well as recovery of small regions on several test images.
NASA Astrophysics Data System (ADS)
Baumann, Sean M.; Keenan, Cameron; Marciniak, Michael A.; Perram, Glen P.
2014-10-01
A database of spectral and temperature-dependent emissivities was created for painted Al-alloy laser-damage-testing targets for the purpose of improving the uncertainty to which temperature on the front and back target surfaces may be estimated during laser-damage testing. Previous temperature estimates had been made by fitting an assumed gray-body radiance curve to the calibrated spectral radiance data collected from the back surface using a Telops Imaging Fourier Transform Spectrometer (IFTS). In this work, temperature-dependent spectral emissivity measurements of the samples were made from room temperature to 500 °C using a Surface Optics Corp. SOC-100 Hemispherical Directional Reflectometer (HDR) with Nicolet FTS. Of particular interest was a high-temperature matte-black enamel paint used to coat the rear surfaces of the Al-alloy samples. The paint had been assumed to have a spectrally flat and temperatureinvariant emissivity. However, the data collected using the HDR showed both spectral variation and temperature dependence. The uncertainty in back-surface temperature estimation during laser-damage testing made using the measured emissivities was improved from greater than +10 °C to less than +5 °C for IFTS pixels away from the laser burn-through hole, where temperatures never exceeded those used in the SOC-100 HDR measurements. At beam center, where temperatures exceeded those used in the SOC-100 HDR, uncertainty in temperature estimates grew beyond those made assuming gray-body emissivity. Accurate temperature estimations during laser-damage testing are useful in informing a predictive model for future high-energy-laser weapon applications.
Circular Conditional Autoregressive Modeling of Vector Fields.
Modlin, Danny; Fuentes, Montse; Reich, Brian
2012-02-01
As hurricanes approach landfall, there are several hazards for which coastal populations must be prepared. Damaging winds, torrential rains, and tornadoes play havoc with both the coast and inland areas; but, the biggest seaside menace to life and property is the storm surge. Wind fields are used as the primary forcing for the numerical forecasts of the coastal ocean response to hurricane force winds, such as the height of the storm surge and the degree of coastal flooding. Unfortunately, developments in deterministic modeling of these forcings have been hindered by computational expenses. In this paper, we present a multivariate spatial model for vector fields, that we apply to hurricane winds. We parameterize the wind vector at each site in polar coordinates and specify a circular conditional autoregressive (CCAR) model for the vector direction, and a spatial CAR model for speed. We apply our framework for vector fields to hurricane surface wind fields for Hurricane Floyd of 1999 and compare our CCAR model to prior methods that decompose wind speed and direction into its N-S and W-E cardinal components. PMID:24353452
Circular Conditional Autoregressive Modeling of Vector Fields*
Modlin, Danny; Fuentes, Montse; Reich, Brian
2013-01-01
As hurricanes approach landfall, there are several hazards for which coastal populations must be prepared. Damaging winds, torrential rains, and tornadoes play havoc with both the coast and inland areas; but, the biggest seaside menace to life and property is the storm surge. Wind fields are used as the primary forcing for the numerical forecasts of the coastal ocean response to hurricane force winds, such as the height of the storm surge and the degree of coastal flooding. Unfortunately, developments in deterministic modeling of these forcings have been hindered by computational expenses. In this paper, we present a multivariate spatial model for vector fields, that we apply to hurricane winds. We parameterize the wind vector at each site in polar coordinates and specify a circular conditional autoregressive (CCAR) model for the vector direction, and a spatial CAR model for speed. We apply our framework for vector fields to hurricane surface wind fields for Hurricane Floyd of 1999 and compare our CCAR model to prior methods that decompose wind speed and direction into its N-S and W-E cardinal components. PMID:24353452
NASA Technical Reports Server (NTRS)
Tittle, R. A.
1988-01-01
The primary purpose of many in-situ airborne light scattering experiments in natural waters is to spectrally characterize the subsurface fluorescent organics and estimate their relative concentrations. This is often done by shining a laser beam into the water and monitoring its subsurface return signal. To do this with the proper interpretation, depth must be taken into account. If one disregards depth dependence when taking such estimates, both their spectral characteristics and their concentrations estimates can be rather ambiguous. A simple airborne lidar configuration is used to detect the subsurface return signal from a particular depth and wavelength. Underwater scatterometer were employed to show that in-situ subsurface organics are very sensitive to depth, but they also require the use of slow moving boats to cover large sample areas. Also, their very entry into the water disturbs the sample it is measuring. The method described is superior and simplest to any employed thus far.
Stable continuous-time autoregressive process driven by stable subordinator
NASA Astrophysics Data System (ADS)
Wyłomańska, Agnieszka; Gajda, Janusz
2016-02-01
In this paper we examine the continuous-time autoregressive moving average process driven by α-stable Lévy motion delayed by inverse stable subordinator. This process can be applied to high-frequency data with visible jumps and so-called "trapping-events". Those properties are often visible in financial time series but also in amorphous semiconductors, technical data describing the rotational speed of a machine working under various load regimes or data related to indoor air quality. We concentrate on the main characteristics of the examined subordinated process expressed in the language of the measures of dependence which are main tools used in statistical investigation of real data. However, because the analyzed system is based on the α-stable distribution therefore we cannot consider here the correlation (or covariance) as a main measure which indicates at the dependence inside the process. In the paper we examine the codifference, the more general measure of dependence defined for wide class of processes. Moreover we present the simulation procedure of the considered system and indicate how to estimate its parameters. The theoretical results we illustrate by the simulated data analysis.
A fuzzy-autoregressive model of daily river flows
NASA Astrophysics Data System (ADS)
Greco, R.
2012-04-01
A model for the identification of daily river flows has been developed, consisting of the combination of an autoregressive model with a fuzzy inference system. The AR model is devoted to the identification of base flow, supposed to be described by linear laws. The fuzzy model identifies the surface runoff, by applying a small set of linguistic statements, deriving from the knowledge of the physical features of the non linear rainfall-runoff transformation, to the inflow entering the river basin. The model has been applied to the identification of the daily flow series of river Volturno at Cancello-Arnone (Southern Italy), with a drainage basin of around 5560 km2, observed between 1970 and 1974. The inflow was estimated on the basis of daily precipitations registered during the same years at six rain gauges located throughout the basin. The first two years were used for model training, the remaining three for the validation. The obtained results show that the proposed model provides good predictions of either low river flows or high floods, although the analysis of residuals, which do not turn out to be a white noise, indicates that the cause and effect relationship between rainfall and runoff is not completely identified by the model.
A fuzzy-autoregressive model of daily river flows
NASA Astrophysics Data System (ADS)
Greco, Roberto
2012-06-01
A model for the identification of daily river flows has been developed, consisting of the combination of an autoregressive model with a fuzzy inference system. The AR model is devoted to the identification of base flow, supposed to be described by linear laws. The fuzzy model identifies the surface runoff, by applying a small set of linguistic statements, deriving from the knowledge of the physical features of the nonlinear rainfall-runoff transformation, to the inflow entering the river basin. The model has been applied to the identification of the daily flow series of river Volturno at Cancello-Arnone (Southern Italy), with a drainage basin of around 5560 km2, observed between 1970 and 1974. The inflow was estimated on the basis of daily precipitations registered during the same years at six rain gauges located throughout the basin. The first two years were used for model training, the remaining three for the validation. The obtained results show that the proposed model provides good predictions of either low river flows or high floods, although the analysis of residuals, which do not turn out to be a white noise, indicates that the cause and effect relationship between rainfall and runoff is not completely identified by the model.
ANALYSIS OF ROLLING GROUP THERAPY DATA USING CONDITIONALLY AUTOREGRESSIVE PRIORS
Paddock, Susan M.; Hunter, Sarah B.; Watkins, Katherine E.; McCaffrey, Daniel F.
2010-01-01
Group therapy is a central treatment modality for behavioral health disorders such as alcohol and other drug use (AOD) and depression. Group therapy is often delivered under a rolling (or open) admissions policy, where new clients are continuously enrolled into a group as space permits. Rolling admissions policies result in a complex correlation structure among client outcomes. Despite the ubiquity of rolling admissions in practice, little guidance on the analysis of such data is available. We discuss the limitations of previously proposed approaches in the context of a study that delivered group cognitive behavioral therapy for depression to clients in residential substance abuse treatment. We improve upon previous rolling group analytic approaches by fully modeling the interrelatedness of client depressive symptom scores using a hierarchical Bayesian model that assumes a conditionally autoregressive prior for session-level random effects. We demonstrate improved performance using our method for estimating the variance of model parameters and the enhanced ability to learn about the complex correlation structure among participants in rolling therapy groups. Our approach broadly applies to any group therapy setting where groups have changing client composition. It will lead to more efficient analyses of client-level data and improve the group therapy research community’s ability to understand how the dynamics of rolling groups lead to client outcomes. PMID:21857889
ANALYSIS OF ROLLING GROUP THERAPY DATA USING CONDITIONALLY AUTOREGRESSIVE PRIORS.
Paddock, Susan M; Hunter, Sarah B; Watkins, Katherine E; McCaffrey, Daniel F
2011-06-01
Group therapy is a central treatment modality for behavioral health disorders such as alcohol and other drug use (AOD) and depression. Group therapy is often delivered under a rolling (or open) admissions policy, where new clients are continuously enrolled into a group as space permits. Rolling admissions policies result in a complex correlation structure among client outcomes. Despite the ubiquity of rolling admissions in practice, little guidance on the analysis of such data is available. We discuss the limitations of previously proposed approaches in the context of a study that delivered group cognitive behavioral therapy for depression to clients in residential substance abuse treatment. We improve upon previous rolling group analytic approaches by fully modeling the interrelatedness of client depressive symptom scores using a hierarchical Bayesian model that assumes a conditionally autoregressive prior for session-level random effects. We demonstrate improved performance using our method for estimating the variance of model parameters and the enhanced ability to learn about the complex correlation structure among participants in rolling therapy groups. Our approach broadly applies to any group therapy setting where groups have changing client composition. It will lead to more efficient analyses of client-level data and improve the group therapy research community's ability to understand how the dynamics of rolling groups lead to client outcomes. PMID:21857889
Dealing with Multiple Solutions in Structural Vector Autoregressive Models.
Beltz, Adriene M; Molenaar, Peter C M
2016-01-01
Structural vector autoregressive models (VARs) hold great potential for psychological science, particularly for time series data analysis. They capture the magnitude, direction of influence, and temporal (lagged and contemporaneous) nature of relations among variables. Unified structural equation modeling (uSEM) is an optimal structural VAR instantiation, according to large-scale simulation studies, and it is implemented within an SEM framework. However, little is known about the uniqueness of uSEM results. Thus, the goal of this study was to investigate whether multiple solutions result from uSEM analysis and, if so, to demonstrate ways to select an optimal solution. This was accomplished with two simulated data sets, an empirical data set concerning children's dyadic play, and modifications to the group iterative multiple model estimation (GIMME) program, which implements uSEMs with group- and individual-level relations in a data-driven manner. Results revealed multiple solutions when there were large contemporaneous relations among variables. Results also verified several ways to select the correct solution when the complete solution set was generated, such as the use of cross-validation, maximum standardized residuals, and information criteria. This work has immediate and direct implications for the analysis of time series data and for the inferences drawn from those data concerning human behavior. PMID:27093380
NASA Astrophysics Data System (ADS)
Dettmer, Jan; Molnar, Sheri; Steininger, Gavin; Dosso, Stan E.; Cassidy, John F.
2012-02-01
This paper applies a general trans-dimensional Bayesian inference methodology and hierarchical autoregressive data-error models to the inversion of microtremor array dispersion data for shear wave velocity (vs) structure. This approach accounts for the limited knowledge of the optimal earth model parametrization (e.g. the number of layers in the vs profile) and of the data-error statistics in the resulting vs parameter uncertainty estimates. The assumed earth model parametrization influences estimates of parameter values and uncertainties due to different parametrizations leading to different ranges of data predictions. The support of the data for a particular model is often non-unique and several parametrizations may be supported. A trans-dimensional formulation accounts for this non-uniqueness by including a model-indexing parameter as an unknown so that groups of models (identified by the indexing parameter) are considered in the results. The earth model is parametrized in terms of a partition model with interfaces given over a depth-range of interest. In this work, the number of interfaces (layers) in the partition model represents the trans-dimensional model indexing. In addition, serial data-error correlations are addressed by augmenting the geophysical forward model with a hierarchical autoregressive error model that can account for a wide range of error processes with a small number of parameters. Hence, the limited knowledge about the true statistical distribution of data errors is also accounted for in the earth model parameter estimates, resulting in more realistic uncertainties and parameter values. Hierarchical autoregressive error models do not rely on point estimates of the model vector to estimate data-error statistics, and have no requirement for computing the inverse or determinant of a data-error covariance matrix. This approach is particularly useful for trans-dimensional inverse problems, as point estimates may not be representative of the
Best estimation of spectrum profiles for diagnosing femoral prostheses loosening.
Díaz-Pérez, Francisco; García-Nieto, Evelyn; Ros, Antonio; Claramunt, Rafael
2014-02-01
For the past few years, some authors have proposed several vibration analysis techniques to detect the prosthetic femoral stem loosening, having found some differences in the frequency response between secure and loose stems. Classical methods like periodogram have been used in most studies for the spectral estimation, and their conclusions have been reached only by visual inspection. A new metric called Non-linear Logarithmic Weighted Distance (NLWD), based on log-spectral distance is presented. As its name suggests, the spectral power is weighted in order to highlight discriminatory patterns of the spectral profiles. A Generalized Discriminant Ratio (GDR) based on NLWD metric has been also defined. In this study, experiments on a cadaveric dried bone with two kinds of fixation, Loose Stem class (LS) and Secure Stem class (SS), have been analyzed. To select the most discriminating approach to spectral estimation, five well known algorithms (Welch's, Burg's Auto-Regressive (AR), Auto-Regressive Moving Average (ARMA), Multiple Signal Classification (MUSIC) and Thomson's Multi-taper (MTM)) have been compared by using GDR. Finally, the use of the MTM method is proposed for the analysis of bone-stem interface vibratory signals, since it yields the most discriminatory profiles. PMID:24332977
Estimation of site-dependent spectral decay parameter from seismic array data
NASA Astrophysics Data System (ADS)
Park, Seon Jeong; Lee, Jung Mo; Baag, Chang-Eob; Choi, Hoseon; Noh, Myunghyun
2016-04-01
The kappa (κ), attenuation of acceleration amplitude at high frequencies, is one of the most important parameters in ground motion evaluation and seismic hazard analysis at sites. κ simply indicates the high frequency decay of the acceleration spectrum in log-linear space. The decay trend can be considered as linear for frequencies higher than a specific frequency, fe which is starting point of the linear regression at the acceleration spectrum. The κ has been investigated using the data from seismic arrays in the south-eastern part of Korea in which nuclear facilities such as power plant and radiological waste depository are located. The seismic array consists of 20 seismic stations and it was operated from October in 2010 through March in 2013. A classical method by Anderson and Hough (1984) and a standard procedure recently suggested by Ktenidou et al. (2013) were applied for computation of κ. There have been just a few studies on spectral attenuation characteristics for Korean Peninsula so far and even those studies utilized small amount of earthquake events whose frequency range was lower than 25 Hz. In this study, the available frequency range is up to 60 Hz based on the sampling rate of 200 and instrument response. This allows us to use a large range of frequencies for κ computations. It is outstanding advantage that we couldn't obtain from earlier κ studies in Korea. In addition, we investigate the regional κ characteristics through calculating the κ using data of 20 seismic stations which are highly extensive seismic array. It allows us to find the more specific attenuation characteristics of high frequencies in study area. Distance and magnitude dependence of κ has also been investigated. Before calculating the κ, the corner frequency (f_c) has been checked so that the fe can lie to the right of fc to exclude source effects in the computation. Manually picked fe is generally in the range of 10 to 25 Hz. The resulting κR is 9.2e-06 and κ0 is 0
NASA Astrophysics Data System (ADS)
Hirose, Misa; Akaho, Rina; Maita, Chikashi; Sugawara, Mai; Tsumura, Norimichi
2016-06-01
In this paper, the spectral sensitivities of a mosaic five-band camera were optimized using a numerical skin phantom to perform the separation of chromophore densities, shading and surface reflection. To simulate the numerical skin phantom, the spectral reflectance of skin was first calculated by Monte Carlo simulation of photon migration for different concentrations of melanin, blood and oxygen saturation levels. The melanin and hemoglobin concentration distributions used in the numerical skin phantom were obtained from actual skin images by independent component analysis. The calculated components were assigned as concentration distributions. The spectral sensitivities of the camera were then optimized using a nonlinear technique to estimate the spectral reflectance for skin separation. In this optimization, the spectral sensitivities were assumed to be normally distributed, and the sensor arrangement was identical to that of a conventional mosaic five-band camera. Our findings demonstrated that spectral estimation could be significantly improved by optimizing the spectral sensitivities.
NASA Astrophysics Data System (ADS)
Hirose, Misa; Akaho, Rina; Maita, Chikashi; Sugawara, Mai; Tsumura, Norimichi
2016-02-01
In this paper, the spectral sensitivities of a mosaic five-band camera were optimized using a numerical skin phantom to perform the separation of chromophore densities, shading and surface reflection. To simulate the numerical skin phantom, the spectral reflectance of skin was first calculated by Monte Carlo simulation of photon migration for different concentrations of melanin, blood and oxygen saturation levels. The melanin and hemoglobin concentration distributions used in the numerical skin phantom were obtained from actual skin images by independent component analysis. The calculated components were assigned as concentration distributions. The spectral sensitivities of the camera were then optimized using a nonlinear technique to estimate the spectral reflectance for skin separation. In this optimization, the spectral sensitivities were assumed to be normally distributed, and the sensor arrangement was identical to that of a conventional mosaic five-band camera. Our findings demonstrated that spectral estimation could be significantly improved by optimizing the spectral sensitivities.
Tanikawa, Tomonori; Li, Wei; Kuchiki, Katsuyuki; Aoki, Teruo; Hori, Masahiro; Stamnes, Knut
2015-11-30
A new retrieval algorithm for estimation of snow grain size and impurity concentration from spectral radiation data is developed for remote sensing applications. A radiative transfer (RT) model for the coupled atmosphere-snow system is used as a forward model. This model simulates spectral radiant quantities for visible and near-infrared channels. The forward RT calculation is, however, the most time-consuming part of the forward-inverse modeling. Therefore, we replaced it with a neural network (NN) function for fast computation of radiances and Jacobians. The retrieval scheme is based on an optimal estimation method with a priori constraints. The NN function was also employed to obtain an accurate first guess in the retrieval scheme. Validation with simulation data shows that a combination of NN techniques and optimal estimation method can provide more accurate retrievals than by using only NN techniques. In addition, validation with in-situ measurements conducted by using ground-based spectral radiometer system shows that comparison between retrieved snow parameters with in-situ measurements is acceptable with satisfactory accuracy. The algorithm provides simultaneous, accurate and fast retrieval of the snow properties. The algorithm presented here is useful for airborne/satellite remote sensing. PMID:26698793
Xu, Haojie; Lu, Yunfeng; Zhu, Shanan
2014-01-01
It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The non-zero covariance of the model’s residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the “causal ordering” is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In the present study, we firstly investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in
NASA Astrophysics Data System (ADS)
Homolová, L.; Janoutová, R.; Malenovský, Z.
2016-06-01
In this study we evaluated various spectral inputs for retrieval of forest chlorophyll content (Cab) and leaf area index (LAI) from high spectral and spatial resolution airborne imaging spectroscopy data collected for two forest study sites in the Czech Republic (beech forest at Štítná nad Vláří and spruce forest at Bílý Kříž). The retrieval algorithm was based on a machine learning method - support vector regression (SVR). Performance of the four spectral inputs used to train SVR was evaluated: a) all available hyperspectral bands, b) continuum removal (CR) 645 - 710 nm, c) CR 705 - 780 nm, and d) CR 680 - 800 nm. Spectral inputs and corresponding SVR models were first assessed at the level of spectral databases simulated by combined leaf-canopy radiative transfer models PROSPECT and DART. At this stage, SVR models using all spectral inputs provided good performance (RMSE for Cab < 10 μg cm-2 and for LAI < 1.5), with consistently better performance for beech over spruce site. Since application of trained SVRs on airborne hyperspectral images of the spruce site produced unacceptably overestimated values, only the beech site results were analysed. The best performance for the Cab estimation was found for CR bands in range of 645 - 710 nm, whereas CR bands in range of 680 - 800 nm were the most suitable for LAI retrieval. The CR transformation reduced the across-track bidirectional reflectance effect present in airborne images due to large sensor field of view.
Technology Transfer Automated Retrieval System (TEKTRAN)
Vegetation indices derived from multispectral imagery are commonly used to extract crop growth and yield information. Spectral unmixing techniques provide an alternative approach to quantifying crop canopy abundance within each image pixel and have the potential for mapping crop yield variability. T...
Seasonal Patterns and Remote Spectral Estimation of Canopy Chemistry Across the Oregon Transect
NASA Technical Reports Server (NTRS)
Matson, Pamela; Johnson, Lee; Billow, Christine; Miller, John; Pu, Ruiliang
1994-01-01
We examined seasonal changes in canopy chemical concentrations and content in conifer forests growing along a climate gradient in western Oregon, as part of the Oregon Transect Ecosystem Research (OTTER) study. The chemical variables were related to seasonal patterns of growth and production. Statistical comparisons of chemical variables with data collected from two different airborne remote-sensing platforms were also carried out. Total nitrogen (N) concentrations in foliage varied significantly both seasonally and among sites; when expressed as content in the forest canopy, nitrogen varied to a much greater extent and was significantly related to aboveground net primary production (r = 0.99). Chlorophyll and free amino acid concentrations varied more strongly than did total N and may have reflected changes in physiological demands for N. Large variations in starch concentrations were measured from pre- to post-budbreak in all conifer sites. Examination of remote-sensing data from two different airborne instruments suggests the potential for remote measurement of some canopy chemicals. Multivariate analysis of high-resolution spectral data in the near infrared region indicated significant correlations between spectral signals and N concentration and canopy N content; the correlation with canopy N content was stronger and was probably associated in part with water absorption features of the forest canopy. The spectral bands that were significantly correlated with lignin concentration and content were similar to bands selected in the other laboratory and airborne studies; starch concentrations were not significantly related to spectral reflectance data. Strong relationships between the spectral position of specific reflectance features in the visible region and chlorophyll were also found.
NASA Technical Reports Server (NTRS)
Rutledge, Charles K.
1988-01-01
The validity of applying chi-square based confidence intervals to far-field acoustic flyover spectral estimates was investigated. Simulated data, using a Kendall series and experimental acoustic data from the NASA/McDonnell Douglas 500E acoustics test, were analyzed. Statistical significance tests to determine the equality of distributions of the simulated and experimental data relative to theoretical chi-square distributions were performed. Bias and uncertainty errors associated with the spectral estimates were easily identified from the data sets. A model relating the uncertainty and bias errors to the estimates resulted, which aided in determining the appropriateness of the chi-square distribution based confidence intervals. Such confidence intervals were appropriate for nontonally associated frequencies of the experimental data but were inappropriate for tonally associated estimate distributions. The appropriateness at the tonally associated frequencies was indicated by the presence of bias error and noncomformity of the distributions to the theoretical chi-square distribution. A technique for determining appropriate confidence intervals at the tonally associated frequencies was suggested.
NASA Astrophysics Data System (ADS)
Gusev, A. A.; Guseva, E. M.
2016-07-01
The parameters of S-wave attenuation (the total effect of absorption and scattering) near the Petropavlovsk (PET) station in Kamchatka were estimated by means of the spectral method through an original procedure. The spectral method typically analyzes the changes with distance of the shape of spectra of the acceleration records assuming that the acceleration spectrum at the earthquake source is flat. In reality, this assumption is violated: the source acceleration spectra often have a high-frequency cutoff (the source-controlled f max) which limits the spectral working bandwidth. Ignoring this phenomenon not only leads to a broad scatter of the individual estimates but also causes systematic errors in the form of overestimation of losses. In the approach applied in the present study, we primarily estimated the frequency of the mentioned high-frequency cutoff and then constructed the loss estimates only within the frequency range where the source spectrum is approximately flat. The shape of the source spectrum was preliminarily assessed by the approximate loss compensation technique. For this purpose, we used the tentative attenuation estimates which are close to the final ones. The difference in the logarithms of the spectral amplitudes at the edges of the working bandwidth is the input for calculating the attenuation. We used the digital accelerograms from the PET station, with 80 samples per second digitization rate, and based on them, we calculated the averaged spectrum of the S-waves as the root mean square along two horizontal components. Our analysis incorporates 384 spectra from the local earthquakes with M = 4-6.5 at the hypocentral distances ranging from 80 to 220 km. By applying the nonlinear least-square method, we found the following parameters of the loss model: the Q-factor Q 0 = 156 ± 33 at frequency f = 1 Hz for the distance interval r = 0-100 km; the exponent in the power-law relationship describing the growth of the Q-factor with frequency,
An autoregressive growth model for longitudinal item analysis.
Jeon, Minjeong; Rabe-Hesketh, Sophia
2016-09-01
A first-order autoregressive growth model is proposed for longitudinal binary item analysis where responses to the same items are conditionally dependent across time given the latent traits. Specifically, the item response probability for a given item at a given time depends on the latent trait as well as the response to the same item at the previous time, or the lagged response. An initial conditions problem arises because there is no lagged response at the initial time period. We handle this problem by adapting solutions proposed for dynamic models in panel data econometrics. Asymptotic and finite sample power for the autoregressive parameters are investigated. The consequences of ignoring local dependence and the initial conditions problem are also examined for data simulated from a first-order autoregressive growth model. The proposed methods are applied to longitudinal data on Korean students' self-esteem. PMID:26645083
Testing and modelling autoregressive conditional heteroskedasticity of streamflow processes
NASA Astrophysics Data System (ADS)
Wang, W.; van Gelder, P. H. A. J. M.; Vrijling, J. K.; Ma, J.
2005-01-01
Conventional streamflow models operate under the assumption of constant variance or season-dependent variances (e.g. ARMA (AutoRegressive Moving Average) models for deseasonalized streamflow series and PARMA (Periodic AutoRegressive Moving Average) models for seasonal streamflow series). However, with McLeod-Li test and Engle's Lagrange Multiplier test, clear evidences are found for the existence of autoregressive conditional heteroskedasticity (i.e. the ARCH (AutoRegressive Conditional Heteroskedasticity) effect), a nonlinear phenomenon of the variance behaviour, in the residual series from linear models fitted to daily and monthly streamflow processes of the upper Yellow River, China. It is shown that the major cause of the ARCH effect is the seasonal variation in variance of the residual series. However, while the seasonal variation in variance can fully explain the ARCH effect for monthly streamflow, it is only a partial explanation for daily flow. It is also shown that while the periodic autoregressive moving average model is adequate in modelling monthly flows, no model is adequate in modelling daily streamflow processes because none of the conventional time series models takes the seasonal variation in variance, as well as the ARCH effect in the residuals, into account. Therefore, an ARMA-GARCH (Generalized AutoRegressive Conditional Heteroskedasticity) error model is proposed to capture the ARCH effect present in daily streamflow series, as well as to preserve seasonal variation in variance in the residuals. The ARMA-GARCH error model combines an ARMA model for modelling the mean behaviour and a GARCH model for modelling the variance behaviour of the residuals from the ARMA model. Since the GARCH model is not followed widely in statistical hydrology, the work can be a useful addition in terms of statistical modelling of daily streamflow processes for the hydrological community.
NASA Technical Reports Server (NTRS)
Asrar, G.; Kanemasu, E. T.; Yoshida, M.
1985-01-01
The influence of management practices and solar illumination angle on the leaf area index (LAI) was estimated from measurements of wheat canopy reflectance evaluated by two methods, a regression formula and an indirect technique. The date of planting and the time of irrigation in relation to the stage of plant growth were found to have significant effects on the development of leaves in spring wheat. A reduction in soil moisture adversely affected both the duration and magnitude of the maximum LAI for late planting dates. In general, water stress during vegetative stages resulted in a reduction in maximum LAI, while water stress during the reproductive period shortened the duration of green LAI in spring wheat. Canopy geometry and solar angle also affected the spectral properties of the canopies, and hence the estimated LAI. Increase in solar zenith angles resulted in a general increase in estimated LAI obtained from both methods.
NASA Technical Reports Server (NTRS)
Giono, G.; Katsukawa, Y.; Ishikawa, R.; Narukage, N.; Bando, T.; Kano, R.; Suematsu, Y.; Winebarger, A.; Kobayashi, K.; Auchere, F.
2015-01-01
The Chromospheric Lyman-Alpha SpectroPolarimeter is a sounding rocket experiment design to measure for the first time the polarization signal of the Lyman-Alpha line (121.6nm), emitted in the solar upper-chromosphere and transition region. This instrument aims to detect the Hanle effect's signature hidden in the Ly-alpha polarization, as a tool to probe the chromospheric magnetic field. Hence, an unprecedented polarization accuracy is needed ((is) less than 10 (exp -3). Nevertheless, spatial and spectral resolutions are also crucial to observe chhromospheric feature such as spicules, and to have precise measurement of the Ly-alpha line core and wings. Hence, this poster will present how the telescope and the spectrograph were separately aligned, and their combined spatial and spectral resolutions.
An adaptive spectral estimation technique to detect cavitation in HIFU with high spatial resolution.
Hsieh, Chang-Yu; Probert Smith, Penny; Mayia, Fares; Ye, Guoliang
2011-07-01
In ultrasound-guided high-intensity focused ultrasound (HIFU) therapy, the changes observed on tissue are subtle during treatment; some ultrasound-guided HIFU protocols rely on the observation of significant brightness changes as the indicator of tissue lesions. The occurrence of a distinct hyperechogenic region ("bright-up") around the focus is often associated with acoustic cavitation resulting in microbubble formation, but it may indicate different physical events such as larger bubbles from boiling (known to alter acoustic impedance) or sometimes lesion formation. A reliable method to distinguish and spatially localize these causes within the tissue would assist the control of HIFU delivery, which is the subject of this paper. Spectral analysis of the radio frequency (RF) signal underlying the B-mode image provides more information on the physical cause, but the usual techniques that are methods on the Fourier transform require a long series for good spectral resolution and so they give poor spatial resolution. This paper introduces an active spectral cavitation detection method to attain high spatial resolution (0.15 × 0.15 mm per pixel) through a parametric statistical method (ARMA modeling) used on finite-length data sets, which enables local changes to be identified more easily. This technique uses the characteristics of the signal itself to optimize the model parameters and structure. Its performance is assessed using synthesized cavitation RF data, and it is then demonstrated in ex vivo bovine liver during and after HIFU exposure. The results suggest that good spatial and spectral resolution can be obtained by the design of suitable algorithms. In ultrasound-guided HIFU, the technique provides a useful supplement to B-mode analysis, with no additional time penalty in data acquisition. PMID:21684454
On estimating attenuation from the amplitude of the spectrally whitened ambient seismic field
NASA Astrophysics Data System (ADS)
Weemstra, Cornelis; Westra, Willem; Snieder, Roel; Boschi, Lapo
2014-06-01
Measuring attenuation on the basis of interferometric, receiver-receiver surface waves is a non-trivial task: the amplitude, more than the phase, of ensemble-averaged cross-correlations is strongly affected by non-uniformities in the ambient wavefield. In addition, ambient noise data are typically pre-processed in ways that affect the amplitude itself. Some authors have recently attempted to measure attenuation in receiver-receiver cross-correlations obtained after the usual pre-processing of seismic ambient-noise records, including, most notably, spectral whitening. Spectral whitening replaces the cross-spectrum with a unit amplitude spectrum. It is generally assumed that cross-terms have cancelled each other prior to spectral whitening. Cross-terms are peaks in the cross-correlation due to simultaneously acting noise sources, that is, spurious traveltime delays due to constructive interference of signal coming from different sources. Cancellation of these cross-terms is a requirement for the successful retrieval of interferometric receiver-receiver signal and results from ensemble averaging. In practice, ensemble averaging is replaced by integrating over sufficiently long time or averaging over several cross-correlation windows. Contrary to the general assumption, we show in this study that cross-terms are not required to cancel each other prior to spectral whitening, but may also cancel each other after the whitening procedure. Specifically, we derive an analytic approximation for the amplitude difference associated with the reversed order of cancellation and normalization. Our approximation shows that an amplitude decrease results from the reversed order. This decrease is predominantly non-linear at small receiver-receiver distances: at distances smaller than approximately two wavelengths, whitening prior to ensemble averaging causes a significantly stronger decay of the cross-spectrum.
NASA Astrophysics Data System (ADS)
Kirby, Jon F.
2014-09-01
The effective elastic thickness (Te) is a geometric measure of the flexural rigidity of the lithosphere, which describes the resistance to bending under the application of applied, vertical loads. As such, it is likely that its magnitude has a major role in governing the tectonic evolution of both continental and oceanic plates. Of the several ways to estimate Te, one has gained popularity in the 40 years since its development because it only requires gravity and topography data, both of which are now readily available and provide excellent coverage over the Earth and even the rocky planets and moons of the solar system. This method, the ‘inverse spectral method’, develops measures of the relationship between observed gravity and topography data in the spatial frequency (wavenumber) domain, namely the admittance and coherence. The observed measures are subsequently inverted against the predictions of thin, elastic plate models, giving estimates of Te and other lithospheric parameters. This article provides a review of inverse spectral methodology and the studies that have used it. It is not, however, concerned with the geological or geodynamic significance or interpretation of Te, nor does it discuss and compare Te results from different methods in different provinces. Since the three main aspects of the subject are thin elastic plate flexure, spectral analysis, and inversion methods, the article broadly follows developments in these. The review also covers synthetic plate modelling, and concludes with a summary of the controversy currently surrounding inverse spectral methods, whether or not the large Te values returned in cratonic regions are artefacts of the method, or genuine observations.
Pittman, Jeremy Joshua; Arnall, Daryl Brian; Interrante, Sindy M; Moffet, Corey A; Butler, Twain J
2015-01-01
Non-destructive biomass estimation of vegetation has been performed via remote sensing as well as physical measurements. An effective method for estimating biomass must have accuracy comparable to the accepted standard of destructive removal. Estimation or measurement of height is commonly employed to create a relationship between height and mass. This study examined several types of ground-based mobile sensing strategies for forage biomass estimation. Forage production experiments consisting of alfalfa (Medicago sativa L.), bermudagrass [Cynodon dactylon (L.) Pers.], and wheat (Triticum aestivum L.) were employed to examine sensor biomass estimation (laser, ultrasonic, and spectral) as compared to physical measurements (plate meter and meter stick) and the traditional harvest method (clipping). Predictive models were constructed via partial least squares regression and modeled estimates were compared to the physically measured biomass. Least significant difference separated mean estimates were examined to evaluate differences in the physical measurements and sensor estimates for canopy height and biomass. Differences between methods were minimal (average percent error of 11.2% for difference between predicted values versus machine and quadrat harvested biomass values (1.64 and 4.91 t·ha(-1), respectively), except at the lowest measured biomass (average percent error of 89% for harvester and quad harvested biomass < 0.79 t·ha(-1)) and greatest measured biomass (average percent error of 18% for harvester and quad harvested biomass >6.4 t·ha(-1)). These data suggest that using mobile sensor-based biomass estimation models could be an effective alternative to the traditional clipping method for rapid, accurate in-field biomass estimation. PMID:25635415
Pittman, Jeremy Joshua; Arnall, Daryl Brian; Interrante, Sindy M.; Moffet, Corey A.; Butler, Twain J.
2015-01-01
Non-destructive biomass estimation of vegetation has been performed via remote sensing as well as physical measurements. An effective method for estimating biomass must have accuracy comparable to the accepted standard of destructive removal. Estimation or measurement of height is commonly employed to create a relationship between height and mass. This study examined several types of ground-based mobile sensing strategies for forage biomass estimation. Forage production experiments consisting of alfalfa (Medicago sativa L.), bermudagrass [Cynodon dactylon (L.) Pers.], and wheat (Triticum aestivum L.) were employed to examine sensor biomass estimation (laser, ultrasonic, and spectral) as compared to physical measurements (plate meter and meter stick) and the traditional harvest method (clipping). Predictive models were constructed via partial least squares regression and modeled estimates were compared to the physically measured biomass. Least significant difference separated mean estimates were examined to evaluate differences in the physical measurements and sensor estimates for canopy height and biomass. Differences between methods were minimal (average percent error of 11.2% for difference between predicted values versus machine and quadrat harvested biomass values (1.64 and 4.91 t·ha−1, respectively), except at the lowest measured biomass (average percent error of 89% for harvester and quad harvested biomass < 0.79 t·ha−1) and greatest measured biomass (average percent error of 18% for harvester and quad harvested biomass >6.4 t·ha−1). These data suggest that using mobile sensor-based biomass estimation models could be an effective alternative to the traditional clipping method for rapid, accurate in-field biomass estimation. PMID:25635415
NASA Astrophysics Data System (ADS)
Aranha dos Santos, Valentin; Schmetterer, Leopold; Gröschl, Martin; Garhofer, Gerhard; Werkmeister, René M.
2016-03-01
Dry eye syndrome is a highly prevalent disease of the ocular surface characterized by an instability of the tear film. Traditional methods used for the evaluation of tear film stability are invasive or show limited repeatability. Here we propose a new noninvasive approach to measure tear film thickness using an efficient delay estimator and ultrahigh resolution spectral domain OCT. Silicon wafer phantoms with layers of known thickness and group index were used to validate the estimator-based thickness measurement. A theoretical analysis of the fundamental limit of the precision of the estimator is presented and the analytical expression of the Cramér-Rao lower bound (CRLB), which is the minimum variance that may be achieved by any unbiased estimator, is derived. The performance of the estimator against noise was investigated using simulations. We found that the proposed estimator reaches the CRLB associated with the OCT amplitude signal. The technique was applied in vivo in healthy subjects and dry eye patients. Series of tear film thickness maps were generated, allowing for the visualization of tear film dynamics. Our results show that the central tear film thickness precisely measured in vivo with a coefficient of variation of about 0.65% and that repeatable tear film dynamics can be observed. The presented method has the potential of being an alternative to breakup time measurements (BUT) and could be used in clinical setting to study patients with dry eye disease and monitor their treatments.
NASA Technical Reports Server (NTRS)
Vasquez, R. P.; Klein, J. D.; Barton, J. J.; Grunthaner, F. J.
1981-01-01
A comparison is made between maximum-entropy spectral estimation and traditional methods of deconvolution used in electron spectroscopy. The maximum-entropy method is found to have higher resolution-enhancement capabilities and, if the broadening function is known, can be used with no adjustable parameters with a high degree of reliability. The method and its use in practice are briefly described, and a criterion is given for choosing the optimal order for the prediction filter based on the prediction-error power sequence. The method is demonstrated on a test case and applied to X-ray photoelectron spectra.
Estimation of tissue optical parameters with hyperspectral imaging and spectral unmixing
NASA Astrophysics Data System (ADS)
Lu, Guolan; Qin, Xulei; Wang, Dongsheng; Chen, Zhuo G.; Fei, Baowei
2015-03-01
Early detection of oral cancer and its curable precursors can improve patient survival and quality of life. Hyperspectral imaging (HSI) holds the potential for noninvasive early detection of oral cancer. The quantification of tissue chromophores by spectral unmixing of hyperspectral images could provide insights for evaluating cancer progression. In this study, non-negative matrix factorization has been applied for decomposing hyperspectral images into physiologically meaningful chromophore concentration maps. The approach has been validated by computer-simulated hyperspectral images and in vivo tumor hyperspectral images from a head and neck cancer animal model.
Estimation of Tissue Optical Parameters with Hyperspectral Imaging and Spectral Unmixing
Lu, Guolan; Qin, Xulei; Wang, Dongsheng; Chen, Zhuo Georgia; Fei, Baowei
2015-01-01
Early detection of oral cancer and its curable precursors can improve patient survival and quality of life. Hyperspectral imaging (HSI) holds the potential for noninvasive early detection of oral cancer. The quantification of tissue chromophores by spectral unmixing of hyperspectral images could provide insights for evaluating cancer progression. In this study, non-negative matrix factorization has been applied for decomposing hyperspectral images into physiologically meaningful chromophore concentration maps. The approach has been validated by computer-simulated hyperspectral images and in vivo tumor hyperspectral images from a head and neck cancer animal model. PMID:26855467
NASA Technical Reports Server (NTRS)
Pearson, R. L.; Miller, L. D.; Tucker, C. J.
1976-01-01
A simple hand-held instrument has been designed and constructed to nondestructively estimate above-ground gramineous biomass using radiometric measurements. The prototype unit consists of a modified two-channel digital radiometer interfaced to a pocket calculator. A digital interface was constructed to join electronically and control the radiometer and calculator to enable the radiometer-calculator system to solve a linear conversion solution from radiometric units to estimated biomass. This instrument has been used to estimate radiometrically gramineous biomass in a more efficient fashion with a high degree of accuracy.
Estimations of Mo X-pinch plasma parameters on QiangGuang-1 facility by L-shell spectral analyses
Wu, Jian; Qiu, Aici; Li, Mo; Wang, Liangping; Wu, Gang; Ning, Guo; Qiu, Mengtong; Li, Xingwen
2013-08-15
Plasma parameters of molybdenum (Mo) X-pinches on the 1-MA QiangGuang-1 facility were estimated by L-shell spectral analysis. X-ray radiation from X-pinches had a pulsed width of 1 ns, and its spectra in 2–3 keV were measured with a time-integrated X-ray spectrometer. Relative intensities of spectral features were derived by correcting for the spectral sensitivity of the spectrometer. With an open source, atomic code FAC (flexible atomic code), ion structures, and various atomic radiative-collisional rates for O-, F-, Ne-, Na-, Mg-, and Al-like ionization stages were calculated, and synthetic spectra were constructed at given plasma parameters. By fitting the measured spectra with the modeled, Mo X-pinch plasmas on the QiangGuang-1 facility had an electron density of about 10{sup 21} cm{sup −3} and the electron temperature of about 1.2 keV.
On estimating frequency response function envelopes using the spectral element method and fuzzy sets
NASA Astrophysics Data System (ADS)
Nunes, R. F.; Klimke, A.; Arruda, J. R. F.
2006-04-01
The influence of uncertain input data on response spectra of dynamic structures is considered. Traditionally, frequency response analyses are based on finite or boundary element models of the objective structure. In the case of the mid-frequency range problem, however, a very fine mesh is required to correctly approximate the frequency response. This is particularly problematic in uncertainty modeling where the computational effort is usually increased significantly by the need for multiple runs (e.g. when conducting a Monte Carlo analysis) to achieve reliable results. In this paper, the spectral element method, combined with a fuzzy set-based uncertainty modeling approach, is presented as an appealing alternative, provided that the models are simple enough to yield a spectral element representation. To conduct the fuzzy analysis part, three different implementations of the extension principle of fuzzy arithmetic are applied and compared. The suitability of each method depends on the number of uncertain parameters, the problem characteristics, and the required accuracy of the results. The performance of the proposed approach is illustrated by two test problems, a simple coupled rod and a reinforced plate model. To verify the fuzzy-valued results, a Monte Carlo simulation has also been included.
Chakraborty, Somsubhra; Das, Bhabani S; Ali, Md Nasim; Li, Bin; Sarathjith, M C; Majumdar, K; Ray, D P
2014-03-01
The aim of this study was to investigate the feasibility of using visible near-infrared (VisNIR) diffuse reflectance spectroscopy (DRS) as an easy, inexpensive, and rapid method to predict compost enzymatic activity, which traditionally measured by fluorescein diacetate hydrolysis (FDA-HR) assay. Compost samples representative of five different compost facilities were scanned by DRS, and the raw reflectance spectra were preprocessed using seven spectral transformations for predicting compost FDA-HR with six multivariate algorithms. Although principal component analysis for all spectral pretreatments satisfactorily identified the clusters by compost types, it could not separate different FDA contents. Furthermore, the artificial neural network multilayer perceptron (residual prediction deviation=3.2, validation r(2)=0.91 and RMSE=13.38 μg g(-1) h(-1)) outperformed other multivariate models to capture the highly non-linear relationships between compost enzymatic activity and VisNIR reflectance spectra after Savitzky-Golay first derivative pretreatment. This work demonstrates the efficiency of VisNIR DRS for predicting compost enzymatic as well as microbial activity. PMID:24398221
Sadek, H.S.; Rashad, S.M.; Blank, H.R.
1984-01-01
If proper account is taken of the constraints of the method, it is capable of providing depth estimates to within an accuracy of about 10 percent under suitable circumstances. The estimates are unaffected by source magnetization and are relatively insensitive to assumptions as to source shape or distribution. The validity of the method is demonstrated by analyses of synthetic profiles and profiles recorded over Harrat Rahat, Saudi Arabia, and Diyur, Egypt, where source depths have been proved by drilling.
On the Estimation of T-Wave Alternans Using the Spectral Fast Fourier Transform Method
Armoundas, Antonis A; Mela, Theofanie; Merchant, Faisal M
2012-01-01
BACKGROUND T-wave alternans (TWA), has been associated with increased vulnerability to ventricular tachyarrhythmias and sudden cardiac death (SCD). However, both random (white) noise and (patho)physiologic processes (i.e. premature ventricular contractions [PVCs], heart and respiration rates) may hamper TWA estimation and therefore, lessen its clinical utility for risk stratification. OBJECTIVE To investigate the effect of random noise and certain (patho)physiologic processes on the estimation of TWA using the Fast Fourier Transform (FFT) method and to develop methods to overcome these potential sources of error. METHODS We used a combination of human electrocardiogram data and computer simulations to assess the effects of a PVC, random and colored noise on the accuracy of TWA estimation. RESULTS We quantitatively demonstrate that replacing a “bad” beat with an odd/even median beat is a more accurate approach than replacing it with the overall average or the overall median beat. We also show that phase resetting may have a significant effect on alternans estimation and that estimation of alternans using frequencies greater than 0.4922 cycles/beat in a 128-point FFT provides the most accurate approach for estimating the alternans when phase resetting is likely to occur. Additionally, our data demonstrate that the number of indeterminate TWA tests due to high levels of noise can be reduced when the alternans voltage exceeds a new higher threshold. Also, the amplitude of random noise has a significant effect on alternans estimation and should be considered to adjust the alternans voltage threshold for noise levels greater than 1.8 μV. Finally, we quantitatively demonstrate that colored noise may lead to a false positive or a false negative result. We propose methods to estimate the effect of these (patho)physiologic processes on the alternans estimation in order to determine whether a TWA test is likely to be a true positive or a true negative. CONCLUSION This
NASA Astrophysics Data System (ADS)
Molnar, S.; Dettmer, J.; Steininger, G.; Dosso, S. E.; Cassidy, J. F.
2013-12-01
This paper applies hierarchical, trans-dimensional Bayesian models for earth and residual-error parametrizations to the inversion of microtremor array dispersion data for shear-wave velocity (Vs) structure. The earth is parametrized in terms of flat-lying, homogeneous layers and residual errors are parametrized with a first-order autoregressive data-error model. The inversion accounts for the limited knowledge of the optimal earth and residual error model parametrization (e.g. the number of layers in the Vs profile) in the resulting Vs parameter uncertainty estimates. The assumed earth model parametrization influences estimates of parameter values and uncertainties due to different parametrizations leading to different ranges of data predictions. The support of the data for a particular model is often non-unique and several parametrizations may be supported. A trans-dimensional formulation accounts for this non-uniqueness by including a model-indexing parameter as an unknown so that groups of models (identified by the index) are considered in the results. In addition, serial residual-error correlations are addressed by augmenting the geophysical forward model with a hierarchical autoregressive error model that can account for a wide range of error processes with a small number of parameters. Hence, the limited knowledge about the true statistical distribution of data errors is also accounted for in the earth model parameter estimates, resulting in more realistic uncertainties and parameter values. Hierarchical autoregressive error models do not rely on point estimates of the model vector to estimate residual-error statistics, and have no requirement for computing the inverse or determinant of a covariance matrix. This approach is particularly useful for trans-dimensional inverse problems, as point estimates may not be representative of the state space that spans multiple subspaces of different dimensions. The autoregressive process is restricted to first order and
Kepler AutoRegressive Planet Search: Motivation & Methodology
NASA Astrophysics Data System (ADS)
Caceres, Gabriel; Feigelson, Eric; Jogesh Babu, G.; Bahamonde, Natalia; Bertin, Karine; Christen, Alejandra; Curé, Michel; Meza, Cristian
2015-08-01
The Kepler AutoRegressive Planet Search (KARPS) project uses statistical methodology associated with autoregressive (AR) processes to model Kepler lightcurves in order to improve exoplanet transit detection in systems with high stellar variability. We also introduce a planet-search algorithm to detect transits in time-series residuals after application of the AR models. One of the main obstacles in detecting faint planetary transits is the intrinsic stellar variability of the host star. The variability displayed by many stars may have autoregressive properties, wherein later flux values are correlated with previous ones in some manner. Auto-Regressive Moving-Average (ARMA) models, Generalized Auto-Regressive Conditional Heteroskedasticity (GARCH), and related models are flexible, phenomenological methods used with great success to model stochastic temporal behaviors in many fields of study, particularly econometrics. Powerful statistical methods are implemented in the public statistical software environment R and its many packages. Modeling involves maximum likelihood fitting, model selection, and residual analysis. These techniques provide a useful framework to model stellar variability and are used in KARPS with the objective of reducing stellar noise to enhance opportunities to find as-yet-undiscovered planets. Our analysis procedure consisting of three steps: pre-processing of the data to remove discontinuities, gaps and outliers; ARMA-type model selection and fitting; and transit signal search of the residuals using a new Transit Comb Filter (TCF) that replaces traditional box-finding algorithms. We apply the procedures to simulated Kepler-like time series with known stellar and planetary signals to evaluate the effectiveness of the KARPS procedures. The ARMA-type modeling is effective at reducing stellar noise, but also reduces and transforms the transit signal into ingress/egress spikes. A periodogram based on the TCF is constructed to concentrate the signal
The use of large-area spectral data in wheat yield estimation
NASA Technical Reports Server (NTRS)
Barnett, T. L.; Thompson, D. R.
1982-01-01
Large-area relations between satellite spectral data and end-of-season crop yield were investigated. Green Index Number (GIN) values from Landsat MSS data of sample segments throughout the U.S. Great Plains winter wheat belt in 1978 were correlated to county USDA-SRS reported yields. A linear relation between GIN and yield appeared to exist up to GIN values of 40 or 50, covering cases of severe to moderate stress. In a test on 1978 Texas winter wheat at the county level, GIN values for sample segments in the counties were used in conjunction with an agronomic-meteorological yield model. The combined fit explained significantly more of the observed yield variation at the county level than the agromet model alone.
Modal identification based on Gaussian continuous time autoregressive moving average model
NASA Astrophysics Data System (ADS)
Xiuli, Du; Fengquan, Wang
2010-09-01
A new time-domain modal identification method of the linear time-invariant system driven by the non-stationary Gaussian random force is presented in this paper. The proposed technique is based on the multivariate continuous time autoregressive moving average (CARMA) model. This method can identify physical parameters of a system from the response-only data. To do this, we first transform the structural dynamic equation into the CARMA model, and subsequently rewrite it in the state-space form. Second, we present the exact maximum likelihood estimators of parameters of the continuous time autoregressive (CAR) model by virtue of the Girsanov theorem, under the assumption that the uniformly modulated function is approximately equal to a constant matrix over a very short period of time. Then, based on the relation between the CAR model and the CARMA model, we present the exact maximum likelihood estimators of parameters of the CARMA model. Finally, the modal parameters are identified by the eigenvalue analysis method. Numerical results show that the method we introduced here not only has high precision and robustness, but also has very high computing efficiency. Therefore, it is suitable for real-time modal identification.
A fast, personal-computer based method of estimating molecular weights of organic compounds from low resolution mass I spectra has been thoroughly evaluated. he method is based on a rule-based pattern,recognition/expert system approach which uses empirical linear corrections whic...
Fu, Xiao-Ning; Wang, Jie; Yang, Lin
2013-01-01
It is a typical passive ranging technology that estimation of distance of an object is based on transmission characteristic of infrared radiation, it is also a hotspot in electro-optic countermeasures. Because of avoiding transmitting energy in the detection, this ranging technology will significantly enhance the penetration capability and infrared conceal capability of the missiles or unmanned aerial vehicles. With the current situation in existing passive ranging system, for overcoming the shortage in ranging an oncoming target object with small temperature difference from background, an improved distance estimation scheme was proposed. This article begins with introducing the concept of signal transfer function, makes clear the working curve of current algorithm, and points out that the estimated distance is not unique due to inherent nonlinearity of the working curve. A new distance calculation algorithm was obtained through nonlinear correction technique. It is a ranging formula by using sensing information at 3-5 and 8-12 microm combined with background temperature and field meteorological conditions. The authors' study has shown that the ranging error could be mainly kept around the level of 10% under the condition of the target and background apparent temperature difference equal to +/- 5 K, and the error in estimating background temperature is no more than +/- 15 K. PMID:23586223
Technology Transfer Automated Retrieval System (TEKTRAN)
Hyperspectral forage canopy absorbance was estimated on eight random plots in each of three 1.2 ha common bermudagrass pastures weekly over a period of 9 weeks from June through early August, 2005 using spectroradiometers measuring light reflectance from 410 nm to 1010 nm. Forage in each plot was ...
Fusion of spectral and electrochemical sensor data for estimating soil macronutrients
Technology Transfer Automated Retrieval System (TEKTRAN)
Rapid and efficient quantification of plant-available soil phosphorus (P) and potassium (K) is needed to support variable-rate fertilization strategies. Two methods that have been used for estimating these soil macronutrients are diffuse reflectance spectroscopy in visible and near-infrared (VNIR) w...
Combined Spectral Index to Improve Ground-Based Estimates of Nitrogen Status in Dryland Wheat
Technology Transfer Automated Retrieval System (TEKTRAN)
Recent studies have demonstrated the usefulness of the single ratio Normalized Difference Vegetation Index (NDVI) and ground-based remote sensing for estimating crop yield potential and basing in-season nitrogen (N) fertilizer application. The NDVI is positively related to crop N status and leaf ar...
Estimation of response-spectral values as functions of magnitude, distance, and site conditions
Joyner, W.B.; Boore, David M.
1982-01-01
We have developed empirical predictive equations for the horizontal pseudo-velocity response at 5-percent damping for 12 different periods from 0.1 to 4.0 s. Using a multiple linear-regression method similar to the one we used previously for peak horizontal acceleration and velocity, we analyzed response spectra period by period for 64 records of 12 shallow earthquakes in western North America, including the recent Coyote Lake and Imperial Valley, California, earthquakes. The resulting predictive equations show amplification of the response values at soil sites for periods greater than or equal to 0.5 s, with maximum amplification exceeding a factor of 2 at 1.5 s. For periods less than 0.5 s there is no statistically significant difference between rock sites and the soil sites represented in the data set. These results are consistent with those of several earlier studies. A particularly significant aspect of the predictive equations is that the response values at different periods are different functions of magnitude (confirming earlier results by McGuire and by Trifunac and Anderson). The slope of the least-squares straight line relating log response to moment magnitude ranges from 0.21 at a period of 0.1 s to greater than 0.5 at periods of 1 s and longer. This result indicates that the conventional practice of scaling a constant spectral shape by peak acceleration will not give accurate answers. The Newmark and Hall method of spectral scaling, using both peak acceleration and peak velocity, largely avoids this error. Comparison of our spectra with the Regulatory Guide 1.60 spectrum anchored at the same value at 0.1 s shows that the Regulatory Guide 1.60 spectrum is exceeded at soil sites for a magnitude of 7.5 at all distances for periods greater than about 0.5 s. Comparison of our spectra for soil sites with the corresponding ATC-3 curve of lateral design force coefficients for the highest seismic zone indicates that the ATC-3 curve is exceeded within about 5 km
NASA Astrophysics Data System (ADS)
Leite, Argentina; Paula Rocha, Ana; Eduarda Silva, Maria
2013-06-01
Heart Rate Variability (HRV) series exhibit long memory and time-varying conditional variance. This work considers the Fractionally Integrated AutoRegressive Moving Average (ARFIMA) models with Generalized AutoRegressive Conditional Heteroscedastic (GARCH) errors. ARFIMA-GARCH models may be used to capture and remove long memory and estimate the conditional volatility in 24 h HRV recordings. The ARFIMA-GARCH approach is applied to fifteen long term HRV series available at Physionet, leading to the discrimination among normal individuals, heart failure patients, and patients with atrial fibrillation.
NASA Astrophysics Data System (ADS)
Privette, Jeffrey L.; Eck, Thomas F.; Deering, Donald W.
1997-12-01
In recent years, many computationally efficient bidirectional reflectance models have been developed to account for angular effects in land remote sensing data, particularly those from the NOAA advanced very high resolution radiometer (AVHRR), polarization and directionality of the Earth's reflectances (POLDER), and the planned EOS moderate-resolution imaging spectrometer (MODIS) and multi-angle imaging spectroradiometer (MISR) sensors. In this study, we assessed the relative ability of 10 such models to predict commonly used remote sensing products (nadir reflectance and albedo). Specifically, we inverted each model with ground-based data from the portable apparatus for rapid acquisition of bidirectional observations of the land and atmosphere (PARABOLA) arranged in subsets representative of satellite sampling geometries. We used data from nine land cover types, ranging from soil to grassland (First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE)) to forest (Boreal Ecosystem-Atmosphere Study (BOREAS)). Retrieved parameters were used in forward model runs to estimate nadir reflectance and spectral albedo over a wide range of solar angles. We rank the models by the accuracy of the estimated products and find results to be strongly dependent on the view azimuth angle range of the inversion data, and less dependent on the spectral band and land cover type. Overall, the nonlinear model of Rahman et al. [993] and the linear kernel-driven RossThickLiSparse model [Wanner et al., 1995] were most accurate. The latter was at least 25 times faster to invert than the former. Interestingly, we found these two models were not able to match the various bidirectional reflectance distribution function (BRDF) shapes as well as other models, suggesting their superior performance lies in their ability to be more reliably inverted with sparse data sets. These results should be useful to those interested in the computationally fast normalization
NASA Astrophysics Data System (ADS)
Ghorbanidehno, H.; Kokkinaki, A.; Darve, E. F.; Kitanidis, P. K.
2014-12-01
The Kalman Filter has been widely used for dynamic monitoring in reservoir engineering, and has recently gained popularity in hydrogeologic applications. A common characteristic of such applications is that the physical processes of interest are greatly affected by preferential flow (e.g., contaminant spreading, CO2 leakage), which can only be delineated if the problem is finely discretized into a large number of unknowns. However, for problems with large numbers of unknowns (e.g., larger than 10,000), the Kalman Filter has prohibitively expensive computation and storage costs. The EnKF, which is typically used to reduce the cost of computing the covariance in such cases converges slowly to the best estimate, and for a reasonable number of realizations, the estimate may not be accurate, especially for strongly heterogeneous systems. We present the Spectral Kalman Filter, a new Kalman Filter implementation that has a dramatically reduced computational cost compared to the full Kalman Filter, with comparable or higher accuracy than the EnKF for the same computational cost. Our algorithm's computational efficiency is achieved by a recurrence that updates small cross-covariance matrices instead of large covariance matrices, in combination with a low-rank approximation of the noise covariance matrix. In addition, instead of computing the expensive Jacobian matrix, a matrix-free method is used to obtain sensitivities. Finally, the error of our method can be explicitly controlled by reducing the time between matrix updates. The frequency of these updates is controlled independently from the data assimilation steps. We demonstrate the performance of the Spectral Kalman Filter for the joint estimation of domain properties and state evolution by assimilation of quasi-continuous data during a hypothetical CO2 injection in a heterogeneous domain.
Estimating high mosquito-producing rice fields using spectral and spatial data
NASA Technical Reports Server (NTRS)
Wood, B. L.; Beck, L. R.; Washino, R. K.; Hibbard, K. A.; Salute, J. S.
1992-01-01
The cultivation of irrigated rice provides ideal larval habitat for a number of anopheline vectors of malaria throughout the world. Anopheles freeborni, a potential vector of human malaria, is associated with the nearly 240,000 hectares of irrigated rice grown annually in Northern and Central California; therefore, this species can serve as a model for the study of rice field anopheline population dynamics. Analysis of field data revealed that rice fields with early season canopy development, that are located near bloodmeal sources (i.e., pastures with livestock) were more likely to produce anopheline larvae than fields with less developed canopies located further from pastures. Remote sensing reflectance measurements of early-season canopy development and geographic information system (GIS) measurements of distanes between rice fields and pastures with livestock were combined to distinguish between high and low mosquito-producing rice fields. Using spectral and distance measures in either a discriminant or Bayesian analysis, the identification of high mosquito-producing fields was made with 85 percent accuracy nearly two months before anopheline larval populations peaked. Since omission errors were also minimized by these approaches, they could provide a new basis for directing abatement techniques for the control of malaria vectors.
[Estimation of Winter Wheat Biomass Using Visible Spectral and BP Based Artificial Neural Networks].
Cui, Ri-xian; Liu, Ya-dong; Fu, Jin-dong
2015-09-01
The objective of this study was to evaluate the feasibility of using color digital image analysis and back propagation (BP) based artificial neural networks (ANN) method to estimate above ground biomass at the canopy level of winter wheat field. Digital color images of winter wheat canopies grown under six levels of nitrogen treatments were taken with a digital camera for four times during the elongation stage and at the same time wheat plants were sampled to measure above ground biomass. Canopy cover (CC) and 10 color indices were calculated from winter wheat canopy images by using image analysis program (developed in Microsoft Visual Basic). Correlation analysis was carried out to identify the relationship between CC, 10 color indices and winter wheat above ground biomass. Stepwise multiple linear regression and BP based ANN methods were used to establish the models to estimate winter wheat above ground biomass. The results showed that CC, and two color indices had a significant cor- relation with above ground biomass. CC revealed the highest correlation with winter wheat above ground biomass. Stepwise multiple linear regression model constituting CC and color indices of NDI and b, and BP based ANN model with four variables (CC, g, b and NDI) for input was constructed to estimate winter wheat above ground biomass. The validation results indicate that the model using BP based ANN method has a better performance with higher R2 (0.903) and lower RMSE (61.706) and RRMSE (18.876) in comparation with the stepwise regression model. PMID:26669174
NASA Technical Reports Server (NTRS)
Bell, Thomas L.; Kundu, Prasun K.; Lau, William K. M. (Technical Monitor)
2002-01-01
Validation of satellite remote-sensing methods for estimating rainfall against rain-gauge data is attractive because of the direct nature of the rain-gauge measurements. Comparisons of satellite estimates to rain-gauge data are difficult, however, because of the extreme variability of rain and the fact that satellites view large areas over a short time while rain gauges monitor small areas continuously. In this paper, a statistical model of rainfall variability developed for studies of sampling error in averages of satellite data is used to examine the impact of spatial and temporal averaging of satellite and gauge data on intercomparison results. The model parameters were derived from radar observations of rain, but the model appears to capture many of the characteristics of rain-gauge data as well. The model predicts that many months of data from areas containing a few gauges are required to validate satellite estimates over the areas, and that the areas should be of the order of several hundred km in diameter. Over gauge arrays of sufficiently high density, the optimal areas and averaging times are reduced. The possibility of using time-weighted averages of gauge data is explored.
Speaker height estimation from speech: Fusing spectral regression and statistical acoustic models.
Hansen, John H L; Williams, Keri; Bořil, Hynek
2015-08-01
Estimating speaker height can assist in voice forensic analysis and provide additional side knowledge to benefit automatic speaker identification or acoustic model selection for automatic speech recognition. In this study, a statistical approach to height estimation that incorporates acoustic models within a non-uniform height bin width Gaussian mixture model structure as well as a formant analysis approach that employs linear regression on selected phones are presented. The accuracy and trade-offs of these systems are explored by examining the consistency of the results, location, and causes of error as well a combined fusion of the two systems using data from the TIMIT corpus. Open set testing is also presented using the Multi-session Audio Research Project corpus and publicly available YouTube audio to examine the effect of channel mismatch between training and testing data and provide a realistic open domain testing scenario. The proposed algorithms achieve a highly competitive performance to previously published literature. Although the different data partitioning in the literature and this study may prevent performance comparisons in absolute terms, the mean average error of 4.89 cm for males and 4.55 cm for females provided by the proposed algorithm on TIMIT utterances containing selected phones suggest a considerable estimation error decrease compared to past efforts. PMID:26328721
NASA Astrophysics Data System (ADS)
Yu, Kang; Lenz-Wiedemann, Victoria; Chen, Xinping; Bareth, Georg
2014-11-01
Monitoring in situ chlorophyll (Chl) content in agricultural crop leaves is of great importance for stress detection, nutritional state diagnosis, yield prediction and studying the mechanisms of plant and environment interaction. Numerous spectral indices have been developed for chlorophyll estimation from leaf- and canopy-level reflectance. However, in most cases, these indices are negatively affected by variations in canopy structure and soil background. The objective of this study was to develop spectral indices that can reduce the effects of varied canopy structure and growth stages for the estimation of leaf Chl. Hyperspectral reflectance data was obtained through simulation by a radiative transfer model, PROSAIL, and measurements from canopies of barley comprising different cultivars across growth stages using spectroradiometers. We applied a comprehensive band-optimization algorithm to explore five types of spectral indices: reflectance difference (RD), reflectance ratio (RR), normalized reflectance difference (NRD), difference of reflectance ratio (DRR) and ratio of reflectance difference (RRD). Indirectly using the multiple scatter correction (MSC) theory, we hypothesized that RRD can eliminate adverse effects of soil background, canopy structure and multiple scattering. Published indices and multivariate models such as optimum multiple band regression (OMBR), partial least squares regression (PLSR) and support vector machines for regression (SVR) were also employed. Results showed that the ratio of reflectance difference index (RRDI) optimized for simulated data significantly improved the correlation with Chl (R2 = 0.98, p < 0.0001) and was insensitive to LAI variations (1-8), compared to widely used indices such as MCARI/OSAVI (R2 = 0.64, p < 0.0001) and TCARI/OSAVI (R2 = 0.74, p < 0.0001). The RRDI optimized for barley explained 76% of the variation in Chl and outperformed multivariate models. However, the accuracy decreased when employing the indices
Satellite Estimation of Spectral Surface UV Irradiance. 2; Effect of Horizontally Homogeneous Clouds
NASA Technical Reports Server (NTRS)
Krothov, N.; Herman, J. R.; Bhartia, P. K.; Ahmad, Z.a; Fioletov, V.
1998-01-01
The local variability of UV irradiance at the Earth's surface is mostly caused by clouds in addition to the seasonal variability. Parametric representations of radiative transfer RT calculations are presented for the convenient solution of the transmission T of ultraviolet radiation through plane parallel clouds over a surface with reflectivity R(sub s). The calculations are intended for use with the Total Ozone Mapping Spectrometer (TOMS) measured radiances to obtain the calculated Lambert equivalent scene reflectivity R for scenes with and without clouds. The purpose is to extend the theoretical analysis of the estimation of UV irradiance from satellite data for a cloudy atmosphere. Results are presented for a range of cloud optical depths and solar zenith angles for the cases of clouds over a low reflectivity surface R(sub s) less than 0.1, over a snow or ice surface R(sub s) greater than 0.3, and for transmission through a non-conservative scattering cloud with single scattering albedo omega(sub 0) = 0.999. The key finding for conservative scattering is that the cloud-transmission function C(sub T), the ratio of cloudy-to clear-sky transmission, is roughly C(sub T) = 1 - R(sub c) with an error of less than 20% for nearly overhead sun and snow-free surfaces. For TOMS estimates of UV irradiance in the presence of both snow and clouds, independent information about snow albedo is needed for conservative cloud scattering. For non-conservative scattering with R(sub s) greater than 0.5 (snow) the satellite measured scene reflectance cannot be used to estimate surface irradiance. The cloud transmission function has been applied to the calculation of UV irradiance at the Earth's surface and compared with ground-based measurements.