Читать книгу Change Detection and Image Time Series Analysis 2 - Группа авторов - Страница 12

Let us consider a time series composed of K images, acquired over the same area on K acquisition dates, by up to K optical and SAR different sensors. Each image in the series is generally composed of multiple features (i.e. it is vector-valued), possibly corresponding to distinct spectral bands or radar polarizations. Specifically, indicates the feature vector of pixel (p, q) in the k-th image in the series In general, each sensor may operate at a distinct spatial resolution; hence, a multisensor and multiresolution time series is being considered.

The acquisition times of the images in the series are assumed to be close enough so that no significant changes occur in the land cover of the observed area. In particular, we assume that no abrupt changes (e.g. due to natural disasters such as floods or earthquakes) occur within the overall time span of the series. This assumption makes it possible to use the whole time series to classify the land cover in the scene, by using the benefit of the complementary properties of the images acquired by different sensors and at different spatial resolutions. Furthermore, this assumption may be especially relevant when the temporal dynamic of the ground scene per se is an indicator of land cover membership, such as in the case of forested (e.g. deciduous vs. evergreen) or agricultural areas. We denote as the land cover classes in the scene and as their set. We operate in a supervised framework; hence, we assume that training samples are available for all of these classes.

The overall formulation introduced in Hedhli et al. (2016) to address multitemporal fusion in the case of single-sensor imagery, and based on multiple quad-trees in cascade, is generalized here to take benefit from the images acquired by different sensors and from their mutual synergy. The multiscale topology of the quad-trees and of hierarchical MRFs defined on quad-trees intrinsically allows multiresolution and multisensor data to be naturally fused in the land cover mapping process.

In this framework, two specific algorithms are defined. In the first one, the k-th image in the series is assigned to a separate quad-tree based on its own spatial resolution. A hierarchical MRF is defined on this quad-tree topology, and inference on the resulting probabilistic graphical model is addressed using the Bayesian marginal posterior mode (MPM) criterion (Kato and Zerubia 2012). In the second proposed algorithm, the focus is on a specific case of multimission, multifrequency and multiresolution time series: multifrequency X-band COSMO-SkyMed and C-band RADARSAT-2 SAR images are used alongside optical visible and near-infrared (VNIR) Pléiades data. This scenario is of special current interest, both because of the potential of exploiting the synergy among these missions and especially in view of the recent COSMO-SkyMed Second Generation and RADARSAT Constellation programs. In the case of the second method, different quad-trees are also used, but both optical and SAR data are associated with each quad-tree in order to benefit from the finest resolution available from the considered sensors. Both approaches exploit the potential of hierarchical probabilistic graphical models (Kato and Zerubia 2012) to address challenging problems of multimodal classification of an image time series.