with the second proposed method
In the second proposed method, the case-specific problem of jointly classifying a multimission, multifrequency (radar X band, radar C band and optical VNIR) and multiresolution series of images acquired by COSMO-SkyMed, RADARSAT-2 and Pléiades is addressed. All of these sensors support multiple spatial resolutions, up to 0.5 m for Pléiades, approximately 1 m for COSMO-SkyMed and 1 m × 2 m for RADARSAT-2.
Let us propose a time series of three images, collected by COSMO-SkyMed, RADARSAT-2 and Pléiades, in the same area, at times close enough to assume land cover stability. A case-specific hierarchical model, again using two quad-trees in cascade, is introduced to generate a supervised classification map at the finest among the spatial resolutions of the input images in the series. Unlike in the case of the hierarchical model of the first proposed approach, in which distinct quad-trees are related to different sensors, in this second approach, for each of the two considered radar sensors, the input image is included in a separate quad-tree according to its spatial resolution (see Figure 1.5). These resolutions are generally expected to be coarser than the finest resolution that is observed through Pléiades. Hence, the input Pléiades data are inserted in the finest resolution layers of both aforementioned quad-trees (see Figure 1.5). The rationale is to both benefit from all input multisensor imagery and to map land cover at the finest resolution available. Empty layers of each of the two quad-trees are filled in with wavelet transforms computed from the Pléiades image included in the leaves layer (see Figure 1.5). The Markovian formulation defined in the previous section, the recursive equations [1.5] and the remarks on the power-of-2 relation among the resolutions of the input images in the series are also valid in this case.
Figure 1.5. Cascaded quad-trees associated with the second proposed method. In this example, we assume that the input time series includes Pléiades data at a 0.5 m resolution, COSMO-SkyMed spotlight data at a 1 m resolution and RADARSAT-2 data sampled on a 2 m pixel lattice. Accordingly, the Pléiades image is included in the leaves layer of both quad-trees. The COSMO-SkyMed and RADARSAT-2 images are inserted in the intermediate layers of separate quad-trees, according to their resolutions. The empty layers of both quad-trees are filled in with wavelet transforms of the Pléiades imagery
1.2.4. Multisensor hierarchical MPM inference
Equation [1.5] allows us to calculate the posterior marginal
Specifically, under appropriate assumptions, in Hedhli et al. (2016) we proved that, for all
where
1.2.4.1. Initializing on the first quad-tree
First, in order to initialize the process, classification is performed using only the data included in the first quad-tree through a classical MPM on a single quad-tree. We refer the reader to Laferté et al. (2000) for more detail. Here, we only briefly recall that the algorithm is initialized by choosing the prior distribution on the root of this first quad-tree. In order to favor spatial regularity, this prior is selected here according to a Potts MRF model (Li 2009; Kato and Zerubia 2012). Details can be found in Hedhli et al. (2016). From the perspective of the classification of the input series of multisensor images, the key point is that, after this initialization stage, the posterior marginal
1.2.4.2. First top-down pass
In the first top-down pass, the second quad-tree is swept downward from the root to the leaves to calculate the prior P(cs) recursively. This prior is initialized in each site
[1.7]
This formulation encourages parent and children sites to share the same class label, although it does not deterministically enforce this condition. It also implies a model for the statistical relations between labels in consecutive layers. Here, the parent–child transition probability
where
