Spatiotemporal Continuous Impervious Surface Mapping by Fusion of Landsat Time Series Data and Google Earth Imagery

The monitoring of impervious surfaces in urban areas using remote sensing with fine spatial and temporal resolutions is crucial for monitoring urban development and environmental changes in urban areas. Spatiotemporal super-resolution mapping (STSRM) fuses fine-spatial-coarse-temporal remote sensing data with coarse-spatial-fine-temporal data, allowing for urban impervious surface mapping at both fine-spatial and fine-temporal resolutions. The STSRM involves two main steps: unmixing the coarse-spatial-fine-temporal remote sensing data to class fraction images, and downscaling the fraction images to sub-pixel land cover maps. Yet, challenges exist in each step when applying STSRM in mapping impervious surfaces. First, the impervious surfaces have high spectral variability (i.e., high intra-class and low inter-class variability), which impacts the accurate extraction of sub-pixel scale impervious surface fractions. Second, downscaling the fraction images to sub-pixel land cover maps is an ill-posed problem and would bring great uncertainty and error in the predictions. This paper proposed a new Spatiotemporal Continuous Impervious Surface Mapping (STCISM) method to deal with these challenges in fusing Landsat and Google Earth imagery. The STCISM used the Multiple Endmember Spectral Mixture Analysis and the Fisher Discriminant Analysis to minimize the within-class variability and maximize the between-class variability to reduce the spectral unmixing uncertainty. In addition, the STCISM adopted a new temporal consistency check model to incorporate temporal contextual information to reduce the uncertainty in the time-series impervious surface prediction maps. Unlike the traditional temporal consistency check model that assumed the impervious-to-pervious conversion is unlikely to happen, the new model allowed the bidirectional conversions between pervious and impervious surfaces. The temporal consistency check was used as a post-procession method to correct the errors in the prediction maps. The proposed STCISM method was used to predict time-series impervious surface maps at 5 m resolution of Google Earth image at the Landsat frequency. The results showed that the proposed STCISM outperformed the STSRM model without using the temporal consistency check and the STSRM model using the temporal consistency check based on the unidirectional pervious-to-impervious surface conversion rule.

Multiple Endmember Spectral Mixture Analysis (MESMA) Applied to the Study of Habitat Diversity in the Fine-Grained Landscapes of the Cantabrian Mountains

Heterogeneous and patchy landscapes where vegetation and abiotic factors vary at small spatial scale (fine-grained landscapes) represent a challenge for habitat diversity mapping using remote sensing imagery. In this context, techniques of spectral mixture analysis may have an advantage over traditional methods of land cover classification because they allow to decompose the spectral signature of a mixed pixel into several endmembers and their respective abundances. In this work, we present the application of Multiple Endmember Spectral Mixture Analysis (MESMA) to quantify habitat diversity and assess the compositional turnover at different spatial scales in the fine-grained landscapes of the Cantabrian Mountains (northwestern Iberian Peninsula). A Landsat-8 OLI scene and high-resolution orthophotographs (25 cm) were used to build a region-specific spectral library of the main types of habitats in this region (arboreal vegetation; shrubby vegetation; herbaceous vegetation; rocks–soil and water bodies). We optimized the spectral library with the Iterative Endmember Selection (IES) method and we applied MESMA to unmix the Landsat scene into five fraction images representing the five defined habitats (root mean square error, RMSE ≤ 0.025 in 99.45% of the pixels). The fraction images were validated by linear regressions using 250 reference plots from the orthophotographs and then used to calculate habitat diversity at the pixel (α-diversity: 30 × 30 m), landscape (γ-diversity: 1 × 1 km) and regional (ε-diversity: 110 × 33 km) scales and the compositional turnover (β- and δ-diversity) according to Simpson’s diversity index. Richness and evenness were also computed. Results showed that fraction images were highly related to reference data (R2 ≥ 0.73 and RMSE ≤ 0.18). In general, our findings indicated that habitat diversity was highly dependent on the spatial scale, with values for the Simpson index ranging from 0.20 ± 0.22 for α-diversity to 0.60 ± 0.09 for γ-diversity and 0.72 ± 0.11 for ε-diversity. Accordingly, we found β-diversity to be higher than δ-diversity. This work contributes to advance in the estimation of ecological diversity in complex landscapes, showing the potential of MESMA to quantify habitat diversity in a comprehensive way using Landsat imagery.

Mapping Burned Areas of Mato Grosso State Brazilian Amazon Using Multisensor Datasets

Quantifying forest fires remain a challenging task for the implementation of public policies aimed to mitigate climate change. In this paper, we propose a new method to provide an annual burned area map of Mato Grosso State located in the Brazilian Amazon region, taking advantage of the high spatial and temporal resolution sensors. The method consists of generating the vegetation, soil, and shade fraction images by applying the Linear Spectral Mixing Model (LSMM) to the Landsat-8 OLI (Operational Land Imager), PROBA-V (Project for On-Board Autonomy–Vegetation), and Suomi NPP-VIIRS (National Polar-Orbiting Partnership-Visible Infrared Imaging Radiometer Suite) datasets. The shade fraction images highlight the burned areas, in which values are represented by low reflectance of ground targets, and the mapping was performed using an unsupervised classifier. Burned areas were evaluated in terms of land use and land cover classes over the Amazon, Cerrado and Pantanal biomes in the Mato Grosso State. Our results showed that most of the burned areas occurred in non-forested areas (66.57%) and old deforestation (21.54%). However, burned areas over forestlands (11.03%), causing forest degradation, reached more than double compared with burned areas identified in consolidated croplands (5.32%). The results obtained were validated using the Sentinel-2 data and compared with active fire data and existing global burned areas products, such as the MODIS (Moderate Resolution Imaging Spectroradiometer product) MCD64A1 and MCD45A1, and Fire CCI (ESA Climate Change Initiative) products. Although there is a good visual agreement among the analyzed products, the areas estimated were quite different. Our results presented correlation of 51% with Sentinel-2 and agreement of r2 = 0.31, r2 = 0.29, and r2 = 0.43 with MCD64A1, MCD45A1, and Fire CCI products, respectively. However, considering the active fire data, it was achieved the better performance between active fire presence and burn mapping (92%). The proposed method provided a general perspective about the patterns of fire in various biomes of Mato Grosso State, Brazil, that are important for the environmental studies, specially related to fire severity, regeneration, and greenhouse gas emissions.

Maximum Fraction Images Derived from Year-Based Project for On-Board Autonomy-Vegetation (PROBA-V) Data for the Rapid Assessment of Land Use and Land Cover Areas in Mato Grosso State, Brazil

This paper presents a new approach for rapidly assessing the extent of land use and land cover (LULC) areas in Mato Grosso state, Brazil. The novel idea is the use of an annual time series of fraction images derived from the linear spectral mixing model (LSMM) instead of original bands. The LSMM was applied to the Project for On-Board Autonomy-Vegetation (PROBA-V) 100-m data composites from 2015 (~73 scenes/year, cloud-free images, in theory), generating vegetation, soil, and shade fraction images. These fraction images highlight the LULC components inside the pixels. The other new idea is to reduce these time series to only six single bands representing the maximum and standard deviation values of these fraction images in an annual composite, reducing the volume of data to classify the main LULC classes. The whole image classification process was conducted in the Google Earth Engine platform using the pixel-based random forest algorithm. A set of 622 samples of each LULC class was collected by visual inspection of PROBA-V and Landsat-8 Operational Land Imager (OLI) images and divided into training and validation datasets. The performance of the method was evaluated by the overall accuracy and confusion matrix. The overall accuracy was 92.4%, with the lowest misclassification found for cropland and forestland (<9% error). The same validation data set showed 88% agreement with the LULC map made available by the Landsat-based MapBiomas project. This proposed method has the potential to be used operationally to accurately map the main LULC areas and to rapidly use the PROBA-V dataset at regional or national levels.

Vegetation Fraction Images Derived from PROBA-V Data for Rapid Assessment of Annual Croplands in Brazil

This paper presents a new method for rapid assessment of the extent of annual croplands in Brazil. The proposed method applies a linear spectral mixing model (LSMM) to PROBA-V time series images to derive vegetation, soil, and shade fraction images for regional analysis. We used S10-TOC (10 days synthesis, 1 km spatial resolution, and top-of-canopy) products for Brazil and S5-TOC (five days synthesis, 100 m spatial resolution, and top-of-canopy) products for Mato Grosso State (Brazilian Legal Amazon). Using the time series of the vegetation fraction images of the whole year (2015 in this case), only one mosaic composed with maximum values of vegetation fraction was generated, allowing detecting and mapping semi-automatically the areas occupied by annual crops during the year. The results (100 m spatial resolution map) for the Mato Grosso State were compared with existing global datasets (Finer Resolution Observation and Monitoring—Global Land Cover (FROM-GLC) and Global Food Security—Support Analyses Data (GFSAD30)). Visually those maps present a good agreement, but the area estimated are not comparable since the agricultural class definition are different for those maps. In addition, we found 11.8 million ha of agricultural areas in the entire Brazilian territory. The area estimation for the Mato Grosso State was 3.4 million ha for 1 km dataset and 5.3 million ha for 100 m dataset. This difference is due to the spatial resolution of the PROBA-V datasets used. A coefficient of determination of 0.82 was found between PROBA-V 100 m and Landsat-8 OLI area estimations for the Mato Grosso State. Therefore, the proposed method is suitable for detecting and mapping annual croplands distribution operationally using PROBA-V datasets for regional analysis.