Search

In applications related to environmental monitoring and disaster management, multichannel synthetic aperture radar (SAR) data present a great potential, owing both to their insensitivity to atmospheric and Sun-illumination conditions and to the improved discrimination capability they may provide as compared with single-channel SAR. However, exploiting this potential requires accurate and automatic techniques to generate change maps from (multichannel) SAR images acquired over the same geographic region in different polarizations or at different frequencies at different times. In this paper, a contextual unsupervised change-detection technique (based on a data-fusion approach) is proposed for two-date multichannel SAR images. Each SAR channel is modeled as a distinct information source, and a Markovian approach to data fusion is adopted. A Markov random field model is introduced that combines together the information conveyed by each SAR channel and the spatial contextual information concerning the correlation among neighboring pixels and formulated by using 'energy functions.' In order to address the task of the estimation of the model parameters, the expectation--maximization algorithm is combined with the recently proposed 'method of log-cumulants.' The proposed technique was experimentally validated with semisimulated multipolarization and multifrequency data and with real SIR-C/XSAR images. Index Terms--Change detection, data fusion, expectation--maximization (EM), Markov random fields (MRFs), multichannel SAR, synthetic aperture radar (SAR).

Land surface temperature (LST) and sea surface temperature (SST) are important quantities for many environmental models. Remote sensing is a source of information for their estimation on both regional and global scales. Many algorithms have been devised to estimate LST and SST from satellite data, most of which require a priori information about the surface and the atmosphere. A recently proposed approach involves the use of support vector machines (SVMs). Based on satellite data and corresponding in situ measurements, they generate an approximation of the relation between them, which can subsequently be used to estimate unknown surface temperatures from additional satellite data. Such a strategy requires the user to set several internal parameters. In this paper, a method is proposed for automatically setting these parameters to quasi-optimal values in the sense of minimum estimation errors. This is achieved by minimizing a functional correlated to regression errors (i.e., the 'span-bound' upper bound on the leave-one-out (LOO) error) which can be computed by using only the training set, without need for a further validation set. In order to minimize this functional, Powell's algorithm is adopted, since it is applicable also to nondifferentiable functions. Experimental results yielded by the proposed method are similar in accuracy to those achieved by cross-validation and by a grid search for the parameter configuration which yields the best test-set accuracy. However, the proposed method gives a dramatic reduction in the computational time required, particularly when many training samples are available. Index Terms--Generalization error bounds, land surface temperature (LST), Powell's method, sea surface temperature (SST), supervised regression, support vector machines (SVMs).

This paper proposes a procedure to extract spectral channels of variable handwidths and spectral positions from the hyperspectral image in such a way as to optimize the accuracy for a specific classification problem. In particular each spectral channel (s-band') is obtained by averaging a group of contiguous channels of the hyperspectral image (h-bands'). Therefore, if one wants to define m s-bands, the problem can be formulated as the optimization of the related m starting and m ending h-bands. Toward this end, we propose to adopt, as an optimization criterion. an interclass distance computed on a training set and to generate a sequence of possible solutions by one of three possible search strategies. As the proposed formalization of the problem makes it analogous to a feature-selection problem, the proposed three strategies have been derived by modifying three feature-selection strategies, namely 1) the ''sequential forward selection.' 2) the ''steepest ascent.' and 3) the 'fast constrained search.' Experimental results on a well-known hyperspectral data set confirm the effectiveness of the approach, which yields better results than other widely used methods. The importance of this kind of procedure lies in feature reduction for hyperspectral image classification or in the ease-based design of the spectral hands of a programmable sensor. It represents a special case of feature extraction that is expected to be more powerful than feature selection. The kind of transformation used allows the interpretability of the new features (i.e., the spectral bands) to be saved. Index Terms--Feature extraction, feature reduction, hyperspectral images, remote-sensing image classification, spectral channels.

In the context of remote-sensing image analysis, Markov random field (MRF) models play an important role, due to their ability to integrate the use of contextual information associated with the image data in the analysis process, through the definition of suitable energy functions. However, especially when dealing with MRF models for supervised classification, the estimation of the parameters characterizing an MRF model is still an open issue, which is typically addressed by using time-expensive interactive 'trial-and-error' procedures. In this paper, by focusing on a broad class of MRF models, in which the energy functions can be expressed as linear combinations of distinct contributions (e.g., representing different typologies of contextual interactions), an automatic supervised procedure for the optimization of the weight parameters of the combinations is proposed. Formulating the parameter-estimation problem as the solution of a system of linear inequalities, the problem is solved by extending to the present context the Ho-Kashyap algorithm, which is typically applied by the pattern-recognition community in computing linear discriminant functions for binary classification. The method is validated experimentally on three different (both single-date and multitemporal) datasets, which are endowed with distinct MRF models formalizing the spatial contextual information associated with a single remote-sensing image and/or the spatio-temporal contextual information associated with an image sequence. Index Terms--Automatic parameter setting, Ho-Kashyap (HK) algorithm, iterated conditional mode (ICM), Markov random field (MRF), remote-sensing image classification.

The availability of synthetic aperture radar (SAR) data offers great potential for environmental monitoring due to the insensitiveness of SAR imagery to atmospheric and sunlight-illumination conditions. In addition, the short revisit time provided by future SAR-based missions will allow a huge amount of multitemporal SAR data to become systematically available for monitoring applications. In this paper, the problem of detecting the changes that occurred on the ground by analyzing SAR imagery is addressed by a completely unsupervised approach, i.e., by developing an automatic thresholding technique. The image-ratioing approach to SAR change detection is adopted, and the Kittler and Illingworth minimum-error thresholding algorithm is generalized to take into account the non-Gaussian distribution of the amplitude values of SAR images. In particular, a SAR-specific parametric modeling approach for the ratio image is proposed and integrated into the thresholding process. Experimental results, which confirm the accuracy of the method for real X-band SAR and spaceborne imaging radar C-band images, are presented. Index Terms--Change detection, Kittler and Illingworth method, method of log-cumulants (MoLC), parametric density estimation, synthetic aperture radar (SAR), thresholding.

A distinct parametric estimation technique for synthetic aperture radar (SAR) amplitude data that assumes a Generalized Gaussian (GG) model for complex SAR backscattered signals is described. The tests show that the method-of-log cumulants (MoLC) estimates are acceptable, consistent, and the amplitude probability density function (PDF) is better than other proposed models.

In remotely sensed data analysis, a crucial problem is represented by the need to develop accurate models for the statistics of the pixel intensities. This paper deals with the problem of probability density function (pdf) estimation in the context of synthetic aperture radar (SAR) amplitude data analysis. Several theoretical and heuristic models for the pdfs of SAR data have been proposed in the literature, which have been proved to be effective for different land-cover typologies, thus making the choice of a single optimal parametric pdf a hard task, especially when dealing with heterogeneous SAR data. In this paper, an innovative estimation algorithm is described, which faces such a problem by adopting a finite mixture model for the amplitude pdf, with mixture components belonging to a given dictionary of SAR-specific pdfs. The proposed method automatically integrates the procedures of selection of the optimal model for each component, of parameter estimation, and of optimization of the number of components by combining the stochastic expectation-maximization iterative methodology with the recently developed 'method-of-log-cumulants' for parametric pdf estimation in the case of nonnegative random variables. Experimental results on several real SAR images are reported, showing that the proposed method accurately models the statistics of SAR amplitude data. Index Terms--Finite mixture models (FMMs), parametric estimation, probability density function (pdf) estimation, stochastic expectation maximization (SEM), synthetic aperture radar (SAR) images.

A general problem of supervised remotely sensed image classification assumes prior knowledge to be available for all the thematic classes that are present in the considered dataset. However, the ground-truth map representing that prior knowledge usually does not really describe all the land-cover typologies iii the image, and the generation of a complete training set often represents a time-consuming, difficult and expensive task. This problem affects the performances of supervised classifiers, which erroneously assign each sample drawn from an unknown class to one of the known classes. In the present paper, a classification strategy is described that allows the identification of samples drawn from unknown classes through the application of a suitable Bayesian decision rule. The proposed approach is based on support vector machines (SVMs) for the estimation of probability density functions and on a recursive procedure to generate prior probability estimates for known and unknown classes. In the experiments, both a synthetic dataset and two real datasets were used. Index Terms--Probability density function (pdf) estimation, supervised classification, support vector machines (SVMs), unknown classes.

Markov random fields (MRFs) provide a useful and theoretically well-established tool for integrating temporal contextual information into the classification process. In particular, when dealing with a sequence of temporal images, the usual MRF-based approach consists in adopting a 'cascade' scheme, i.e., in propagating the temporal information from the current image to the next one of the sequence. The simplicity of the cascade scheme makes it attractive; on the other hand, it does not fully exploit the temporal information available in a sequence of temporal images. In this paper, a 'mutual' MRF approach is proposed that aims at improving both the accuracy and the reliability of the classification process by means of a better exploitation of the temporal information. It involves carrying out a bidirectional exchange of the temporal information between the defined single-time MRF models of consecutive images. A difficult issue related to MRFs is the determination of the MRF model parameters that weight the energy terms related to the available information sources. To solve this problem, we propose a simple and fast method based on the concept of 'minimum perturbation' and implemented with the pseudoinverse technique for the minimization of the sum of squared errors. Experimental results on a multitemporal dataset made up of two multisensor (Landsat Thematic Mapper and European Remote Sensing 1 synthetic aperture radar) images are reported. The results obtained by the proposed 'mutual' approach show a clear improvement in terms of classification accuracy over those yielded by a reference MRF-based classifier. The presented method to automatically estimate the MRF parameters yielded significant results that make it an attractive alternative to the usual trial-and-error search procedure. Index Terms--Artificial neural networks, data fusion, Markov random fields (MRFs), Markov random field (MRF) parameter estimation, multitemporal image classification, spatio-temporal context.

In order to apply a statistical approach to the classification of multisource remote-sensing data, one of the main problems to face lies in the estimation of probability distribution functions. This problem arises out of the difficulty of defining a common statistical model for such heterogeneous data. A possible solution is to adopt nonparametric approaches, which rely on the availability of training samples without any assumption about the related statistical distributions. The purpose of this paper is to investigate the suitability of the concept of dependence trees for the integration of multisource information through estimation of probability distributions. First, this concept, introduced by Chow and Liu, is used to provide an approximation of a probability distribution defined in an N-dimensional space by a product of N--1 probability distributions defined in two-dimensional (2-D) spaces; this approximation corresponds, in terms of graph theoretical interpretation, to a tree of dependence. For each land cover class, a dependence tree is generated by minimizing an appropriate closeness measure. Then, a nonparametric estimation of the second-order probability distributions is carried out through the Parzen window approach, based on the implementation of 2-D Gaussian kernels. In this way, it is possible to reduce the complexity of the estimation, while capturing a significant part of the interdependence among variables. A comparison with other multisource data fusion methods, namely, the multilayer perceptron (MLP) method, the k-nearest neighbor (k-NN) method, and a Bayesian hierarchical classifier (BHC), is made. Experimental results obtained on multisensor [airborne thematic mapper (ATM) and synthetic aperture radar (SAR)] and multisource (experimental synthetic aperture radar (E-SAR) and a textural feature) data sets show that the proposed fusion method based on dependence trees is able to provide a classification accuracy similar to those of the other methods considered, but with the advantage of a reduced computational load. Index Terms--Dependence trees, discrete probability distribution approximation, multisource fusion, Parzen window approach, remote sensing image classification.

We propose a supervised nonparametric technique, based on the 'compound classification rule' for minimum error, to detect land-cover transitions between two remote-sensing images acquired at different times. Thanks to a simplifying hypothesis, the compound classification rule is transformed into a form easier to compute. In the obtained rule, an important role is played by the probabilities of transitions, which take into account the temporal dependence between two images. In order to avoid requiring that training sets be representative of all possible types of transitions, we propose an iterative algorithm which allows the probabilities of transitions to be estimated directly from the images under investigation. Experimental results on two Thematic Mapper images confirm that the proposed algorithm may provide remarkably better detection accuracy than the 'Post-Classification Comparison' algorithm, which is based on the separate classifications of the two images.

The problem of extending the Jeffreys-Matusita distance to multiclass cases for feature-selection purposes is addressed and a solution equivalent to the Bhattacharyya bound is presented. This extension is compared with the widely used weighted average Jeffreys-Matusita distance both by examining the respective formulae and by experimenting on an optical remote-sensing data set.

This paper proposes the application of structured neural networks to classification of multisensor remote-sensing images. The purpose of our approach is to allow the interpretation of the 'network behavior,' as it can be utilized by photointerpreters for the validation of the neural classifier. In addition, our approach gives a criterion for defining the network architecture, so avoiding the classical trial-and-error process. First of all, the architecture of structured multilayer feedforward networks is tailored to a multisensor classification problem. Then, such networks are trained to solve the problem by the error backpropagation algorithm. Finally, they are transformed into equivalent networks to obtain a simplified representation. The resulting equivalent networks may be interpreted as a hierarchical arrangement of 'committees' that accomplish the classification task by checking on a set of explicit constraints on input data. Experimental results on a multisensor (optical and SAR) data set are described in terms of both classification accuracy and network interpretation. Comparisons with fully connected neural networks and with the k-nearest neighbor classifier are also made.

This study deals with a learning system for generation of descriptions of objects to be recognized in 2-D images. After proposing a framework for handling fuzzy and relational descriptions, we present the system obtained by making such a framework manage a well-known learning methodology. Satisfactory results and comparisons are reported.

A new suboptimal search strategy suitable for feature selection in very high-dimensional remote sensing images (e.g., those acquired by hyperspectral sensors) is proposed. Each solution of the feature selection problem is represented as a binary string that indicates which features are selected and which are disregarded. In turn, each binary string corresponds to a point of a multidimensional binary space. Given a criterion function to evaluate the effectiveness of a selected solution, the proposed strategy is based on the search for constrained local extremes of such a function in the above-defined binary space. In particular, two different algorithms are presented that explore the space of solutions in different ways. These algorithms are compared with the classical sequential forward selection and sequential forward floating selection suboptimal techniques, using hyperspectral remote sensing images (acquired by the airborne visible/infrared imaging spectrometer [AVIRIS] sensor) as a data set. Experimental results point out the effectiveness of both algorithms, which can be regarded as valid alternatives to classical methods, as they allow interesting tradeoffs between the qualities of selected feature subsets and computational cost. Index Terms--Feature selection, hyperspectral data , remote sensing, search algorithms.

A new preprocessing technique is presented in this paper to automatically highlight changes in multitemporal strongly heterogeneous remotely sensed images. The proposed technique is devoted to the case where the two acquisitions, before and after a given event, are significantly different, due, for instance, to different sensors, acquisition modalities, or climatic conditions. In a previous study, it was proven that the local statistics of the images acquired at the two dates could be used to extract a relevant change indicator. Nevertheless, this measure is valid when the two observations have been derived from similar acquisitions. When the acquisition modalities differ, local statistics tend to be too different from one image to the other one to be relevant in highlighting the ground evolution without mixing with the changes at ground. The technique proposed in this paper to overcome this limitation is based on the assumption that some dependence indeed exists between the two images in unchanged areas. This dependence is modeled by quantile regression applied according to the copula theory and used to perform an estimation of the local statistics that would have been observed if the acquisition conditions of the first image had been similar to the ones of the second image. The method yields an estimate of the local statistics of the first image through the point of view of the second one. Then, usual Kullback-Leibler-based comparisons of those statistics are applied to define a change measure, which may be analyzed (e.g., by thresholding) in order to detect changes. Experimental results are shown to validate the proposed method by using a pair of Synthetic Aperture Radar (SAR) images onboard European Remote Sensing (ERS) Satellite images and a pair of optical-SAR images (from the High Resolution Visible (HRV) sensor onboard Satellite Pour l'Observation de la Terre (SPOT) satellite and from ERS-SAR) acquired before and after a flood. [ABSTRACT FROM AUTHOR]

In the framework of synthetic aperture radar (SAR) systems, current satellite missions make it possible to acquire images at very high and multiple spatial resolutions with short revisit times. This scenario conveys a remarkable potential in applications to, for instance, environmental monitoring and natural disaster recovery. In this context, data fusion and change detection methodologies play major roles. This paper proposes an unsupervised change detection algorithm for the challenging case of multimodal SAR data collected by sensors operating at multiple spatial resolutions. The method is based on Markovian probabilistic graphical models, graph cuts, linear mixtures, generalized Gaussian distributions, Gram–Charlier approximations, maximum likelihood and minimum mean squared error estimation. It benefits from the SAR images acquired at multiple spatial resolutions and with possibly different modalities on the considered acquisition times to generate an output change map at the finest observed resolution. This is accomplished by modeling the statistics of the data at the various spatial scales through appropriate generalized Gaussian distributions and by iteratively estimating a set of virtual images that are defined on the pixel grid at the finest resolution and would be collected if all the sensors could work at that resolution. A Markov random field framework is adopted to address the detection problem by defining an appropriate multimodal energy function that is minimized using graph cuts. [ABSTRACT FROM AUTHOR]

Introduces a series of articles that aims to describe advances in processing and analysis techniques for hyperspectral data as of July 2001. Dense sampling of the spectral range of the electromagnetic spectrum; Suboptimal methods for feature extraction; Spatial resolution of satellite hyperspectral images.