Corinne Mailhes

Quality criteria benchmark for hyperspectral imagery

Hyperspectral data appear to be of a growing interest over the past few years. However, applications for hyperspectral data are still in their infancy as handling the significant size of the data presents a challenge for the user community. Efficient compression techniques are required, and lossy compression, specifically, will have a role to play, provided its impact on remote sensing applications remains insignificant. To assess the data quality, suitable distortion measures relevant to end-user applications are required. Quality criteria are also of a major interest for the conception and development of new sensors to define their requirements and specifications. This paper proposes a method to evaluate quality criteria in the context of hyperspectral images. The purpose is to provide quality criteria relevant to the impact of degradations on several classification applications. Different quality criteria are considered. Some are traditionally used in image and video coding and are adapted here to hyperspectral images. Others are specific to hyperspectral data. We also propose the adaptation of two advanced criteria in the presence of different simulated degradations on AVIRIS hyperspectral images. Finally, five criteria are selected to give an accurate representation of the nature and the level of the degradation affecting hyperspectral data.

Bayesian Estimation of Linear Mixtures Using the Normal Compositional Model. Application to Hyperspectral Imagery

This paper studies a new Bayesian unmixing algorithm for hyperspectral images. Each pixel of the image is modeled as a linear combination of so-called endmembers. These endmembers are supposed to be random in order to model uncertainties regarding their knowledge. More precisely, we model endmembers as Gaussian vectors whose means have been determined using an endmember extraction algorithm such as the famous N-finder (N-FINDR) or Vertex Component Analysis (VCA) algorithms. This paper proposes to estimate the mixture coefficients (referred to as abundances) using a Bayesian algorithm. Suitable priors are assigned to the abundances in order to satisfy positivity and additivity constraints whereas conjugate priors are chosen for the remaining parameters. A hybrid Gibbs sampler is then constructed to generate abundance and variance samples distributed according to the joint posterior of the abundances and noise variances. The performance of the proposed methodology is evaluated by comparison with other unmixing algorithms on synthetic and real images.

Hyperspectral Image Compression: Adapting SPIHT and EZW to Anisotropic 3-D Wavelet Coding

Hyperspectral images present some specific characteristics that should be used by an efficient compression system. In compression, wavelets have shown a good adaptability to a wide range of data, while being of reasonable complexity. Some wavelet-based compression algorithms have been successfully used for some hyperspectral space missions. This paper focuses on the optimization of a full wavelet compression system for hyperspectral images. Each step of the compression algorithm is studied and optimized. First, an algorithm to find the optimal 3-D wavelet decomposition in a rate-distortion sense is defined. Then, it is shown that a specific fixed decomposition has almost the same performance, while being more useful in terms of complexity issues. It is shown that this decomposition significantly improves the classical isotropic decomposition. One of the most useful properties of this fixed decomposition is that it allows the use of zero tree algorithms. Various tree structures, creating a relationship between coefficients, are compared. Two efficient compression methods based on zerotree coding (EZW and SPIHT) are adapted on this near-optimal decomposition with the best tree structure found. Performances are compared with the adaptation of JPEG 2000 for hyperspectral images on six different areas presenting different statistical properties.

P- and T-Wave Delineation in ECG Signals Using a Bayesian Approach and a Partially Collapsed Gibbs Sampler

Detection and delineation of P- and T-waves are important issues in the analysis and interpretation of electrocardiogram (ECG) signals. This paper addresses this problem by using Bayesian inference to represent a priori relationships among ECG wave components. Based on the recently introduced partially collapsed Gibbs sampler principle, the wave delineation and estimation are conducted simultaneously by using a Bayesian algorithm combined with a Markov chain Monte Carlo method. This method exploits the strong local dependency of ECG signals. The proposed strategy is evaluated on the annotated QT database and compared to other classical algorithms. An important feature of this paper is that it allows not only for the detection of P- and T-wave peaks and boundaries, but also for the accurate estimation of waveforms for each analysis window. This can be useful for some ECG analysis that require wave morphology information.

On the choice of an electromyogram data compression method

Electromyogram (EMG) data compression is of great importance within the framework of telemedicine. For example, there is an increasing demand in medicine to achieve patient healthcare directly from the office of the specialist. The aim of the research presented in this paper is to investigate several kinds of compression methods applied to EMG signals in order to find the method that is most well-suited to EMG data compression. Thus, in this paper, the application of several compression methods to EMG data is studied: predictive linear methods, transform methods and, more specifically, methods based on the wavelet transform. Each method is briefly discussed and experimental results are presented in terms of the signal-to-noise ratio and the compression ratio. The results show that methods based on the wavelet transform outperform the other compression methods.

Parameter Estimation for Peaky Altimetric Waveforms

Much attention has been recently devoted to the analysis of coastal altimetric waveforms. When approaching the coast, altimetric waveforms are sometimes corrupted by peaks caused by high reflective areas inside the illuminated land surfaces or by the modification of the sea state close to the shoreline. This paper introduces a new parametric model for these peaky altimetric waveforms. This model assumes that the received altimetric waveform is the sum of a Brown echo and an asymmetric Gaussian peak. The asymmetric Gaussian peak is parameterized by a location, an amplitude, a width, and an asymmetry coefficient. A maximum-likelihood estimator is studied to estimate the Brown plus peak model parameters. The Cramér-Rao lower bounds of the model parameters are then derived providing minimum variances for any unbiased estimator, i.e., a reference in terms of estimation error. The performance of the proposed model and the resulting estimation strategy are evaluated via many simulations conducted on synthetic and real data. Results obtained in this paper show that the proposed model can be used to retrack efficiently standard oceanic Brown echoes as well as coastal echoes corrupted by symmetric or asymmetric Gaussian peaks. Thus, the Brown with Gaussian peak model is useful for analyzing altimetric measurements closer to the coast.

A Semi-Analytical Model for Delay/Doppler Altimetry and Its Estimation Algorithm

The concept of delay/Doppler (DD) altimetry (DDA) has been under study since the mid-1990s, aiming at reducing the measurement noise and increasing the along-track resolution in comparison with the conventional pulse-limited altimetry. This paper introduces a new model for the mean backscattered power waveform acquired by a radar altimeter operating in synthetic aperture radar mode, as well as an associated least squares (LS) estimation algorithm. As in conventional altimetry (CA), the mean power can be expressed as the convolution of three terms: the flat surface impulse response (FSIR), the probability density function of the heights of the specular scatterers, and the time/frequency point target response of the radar. An important contribution of this paper is to derive an analytical formula for the FSIR associated with DDA. This analytical formula is obtained for a circular antenna pattern, no mispointing, no vertical speed effect, and a uniform scattering. The double convolution defining the mean echo power can then be computed numerically, resulting in a 2-D semi-analytical model called the DD map (DDM). This DDM depends on three altimetric parameters: the epoch, the sea surface wave height, and the amplitude. A multi-look model is obtained by summing all the reflected echoes from the same along-track surface location of interest after applying appropriate delay compensation (range migration) to align the DDM on the same reference. The second contribution of this paper concerns the estimation of the parameters associated with the multi-look semi-analytical model. An LS approach is investigated by means of the Levenberg-Marquardt algorithm. Simulations conducted on simulated altimetric waveforms allow the performance of the proposed estimation algorithm to be appreciated. The analysis of Cryosat-2 waveforms shows an improvement in parameter estimation when compared to the CA.

Time-Frequency Tracking of Spectral Structures Estimated by a Data-Driven Method

The installation of a condition monitoring system (CMS) aims to reduce the operating costs of the monitored system by applying a predictive maintenance strategy. However, a system-driven configuration of the CMS requires the knowledge of the system kinematics and could induce a lot of false alarms because of predefined thresholds. The purpose of this paper is to propose a complete data-driven method to automatically generate system health indicators without any a priori on the monitored system or the acquired signals. This method is composed of two steps. First, every acquired signal is analyzed: the spectral peaks are detected and then grouped in a more complex structure as harmonic series or modulation sidebands. Then, a time-frequency tracking operation is applied on all available signals: the spectral peaks and the spectral structures are tracked over time and grouped in trajectories, which will be used to generate the system health indicators. The proposed method is tested on real-world signals coming from a wind turbine test rig. The detection of a harmonic series and a modulation sideband reports the birth of a fault on the main bearing inner ring. The evolution of the fault severity is characterized by three automatically generated health indicators and is confirmed by experts.

Biomonitoring of the Genotoxic Potential of Draining Water from Dredged Sediments, Using the Comet and Micronucleus Tests on Amphibian (Xenopus Laevis) Larvae And Bacterial Assays (Mutatox® and Ames Tests)

Management of contaminated dredged sediments is a matter of great human concern. The present investigation evaluates the genotoxic potential of aqueous extracts of five sediments from French channels (draining water from dredged sediments), using larvae of the frog Xenopus laevis. Two genotoxic endpoints were analyzed in larvae: clastogenic and/or aneugenic effects (micronucleus induction after 12 d of exposure) and DNA-strand breaking potency (comet assay after 1 and 12 d of exposure) in the circulating blood. Additionally, in vitro bacterial assays (Microtox and Ames tests) were carried out and the results were compared with those obtained with larvae. Physicochemical analyses were also taken into account. Analytical analyses highlighted in the five draining waters a heavy load of contaminants such as metals and hydrocarbons. The results obtained with the micronucleus test established the genotoxicity of three draining waters. The comet assay showed that all 5 draining waters were genotoxic after 1 d of exposure. Although 3 of them were still genotoxic after 12 d of exposure, DNA damage globally decreased between d 1 and 12. The comet assay can be considered as a genotoxicity-screening tool. Data indicate that both tests should be used in conjunction in Xenopus. Bacterial tests (Ames) revealed genotoxicity for only one draining water. The results confirm the relevance of the amphibian model and the need to resort to bioassays in vivo such as the Xenopus micronucleus and comet assays for evaluation of the ecotoxicological impact, an essential complement to the physicochemical data.

Including Antenna Mispointing in a Semi-Analytical Model for Delay/Doppler Altimetry

Delay/Doppler altimetry (DDA) aims at reducing the measurement noise and increasing the along-track resolution in comparison with conventional pulse-limited altimetry. In a previous paper, we have proposed a semi-analytical model for DDA, which considers some simplifications as the absence of mispointing antenna. This paper first proposes a new analytical expression for the flat surface impulse response (FSIR), considering antenna mispointing angles, a circular antenna pattern, no vertical speed effect, and uniform scattering. The 2-D delay/Doppler map is then obtained by a numerical computation of the convolution between the proposed analytical function, the probability density function of the heights of the specular scatterers, and the time/frequency point target response of the radar. The approximations used to obtain the semi-analytical model are analyzed, and the associated errors are quantified by analytical bounds for these errors. The second contribution of this paper concerns the estimation of the parameters associated with the multilook semi-analytical model. Two estimation strategies based on the least squares procedure are proposed. The proposed model and algorithms are validated on both synthetic and real waveforms. The obtained results are very promising and show the accuracy of this generalized model with respect to the previous model assuming zero antenna mispointing.