GET TO KNOW THE STEPS OF THE MAPBIOMAS BOLIVIA METHODOLOGY
Here we detail step by step the methodology of MapBiomas Bolivia. For each class treated in the map, there are particularities and specific characteristics that can be consulted in detail in the . ATBD (Algorithm Theory Base Document).
Access full methodology - ATBD
Satellite images from Landsat are processed, with a resolution of 30 meters, freely available on the Google Earth Engine platform for a temporal series of more than 39 years. Around 5,577 tiles were constructed across the national boundary, each containing tens of millions of pixels in total. These pixels are the work units of MapBiomas. Images may contain clouds, haze, and other conditions that may affect their quality. To produce a clean image, the clear pixels are selected from the images available for the selected period. For each of these pixels, metrics are extracted that explain the behavior of the pixel in that year. This is done with each of the satellite's 7 spectral bands, as well as with the calculated spectral fractions and indices. For example, for Band 1, the median of the values of the band in the period, the maximum value, the minimum value in the year and the amplitude of variation are collected. In the end, each pixel for a year carries up to 156 layers of information.
For each year, a mosaic covering the entire country is created, representing the behavior of each pixel through 156 metrics or layers of information. This set of mosaics is saved as a data collection within the Google Earth Engine platform. These mosaics are used in two main ways. First, as a source of parameters for the algorithm to classify the images (see next step). This mosaic is also used to generate the RGB composition that allows visualizing the background image on the MapBiomas platform. This composition is also used to collect training samples and evaluate accuracy through visual interpretation.
From the image mosaics, the land cover and land use map is developed (forest, non-forest natural, agriculture, pasture, non-vegetated area, water bodies, etc.). To do this, analysts at MapBiomas Bolivia use an automatic classifier called "random forest," which runs on Google's cloud processors. This system is based on machine learning: for each class to be classified, the algorithm is "trained" with samples of the targets to be classified. These samples are obtained through reference maps, generation of stable class maps from previous series of MapBiomas Bolivia, and direct collection of complementary samples through visual interpretation of Landsat images.
The spatial filter aims to increase the spatial consistency of the data by eliminating isolated or edge pixels. Rules of neighborhood are defined that can lead to a change in the classification of pixels. For example, a pixel that has fewer than two of its nine neighboring pixels in the same class will be reclassified to the predominant class in the neighborhood. The filter is applied to all classes and years of the collection.
Para reducir las inconsistencias temporales, en particular cambios en la cobertura o el uso que son imposibles o no están permitidos (por ejemplo: Bosque natural > No bosque > Bosque natural) y para corregir fallas debido al desbordamiento de nubes o falta de datos, se aplica un filtro temporal. Cada bioma, tema o región puede tener reglas de filtrado temporales específicas. En total en la Colección 2.0 se aplicaron 3 reglas. El filtro temporal se aplica a cada píxel analizando todos los años de la Colección (por ejemplo, la Colección 2 tiene 39 años).
In this step, the maps of each class are integrated into a single map that represents the land coverage and use of the entire territory for each year. Prevalence rules are applied: in this way, if the same pixel is classified into two different class maps, it is possible to define which one it belongs to in the final map. Prevalence rules may vary depending on the peculiarities of biomes, themes or regions. The integration is performed for each year in the series and generates an integrated map for each year, usually saved as a single ASSET with the number of annual layers for the analyzed period. The embedded map goes through one more spatial filter step to clean up any edges and loose pixels as a result of the embedding process.
To understand changes in land cover and land use, maps with class transitions between different pairs of selected years are produced. This allows for the analysis of changes in the territory and answers questions such as how much of the woody class converted to pasture from one year to another. Transition maps are generated pixel by pixel and, once completed, they are also subjected to a spatial filter to eliminate isolated or edge pixels. Transition matrices are constructed from these maps for different territorial divisions available on the MapBiomas Bolivia platform.