This data sets include yearly maps of land cover classification for the state of Mato Grosso, Brasil, from 2001 to 2016, based on MODIS image time series at 250 meter spatial resolution. Ground samples consisting of 2,115 time series with known labels are used as training data for a support vector machine classifier. The classes include natural and human-transformed land areas, discriminating among different agricultural crops in state of Mato Grosso, Brazil's agricultural frontier. The results provide spatially explicit estimates of productivity increases in agriculture as well as the trade-offs between crop and pasture expansion. Quality assessment using a 5-fold cross-validation of the training samples indicates an overall accuracy of 93% and the following user's and producer's accuracy for the land cover classes:
Cerrado: UA - 99% PA - 98%
Fallow_Cotton UA - 100% PA - 100%
Forest UA - 99% PA - 98%
Pasture UA - 95% PA - 96%
Soy-Corn UA- 87% PA - 97%
Soy-Cotton UA - 99% PA - 94%
Soy-Fallow UA - 100% PA - 100%
Soy-Millet UA- 84% PA - 84%
Soy-Sunflower UA - 85% PA - 85%
The correlation coefficients between the agricultural areas classified by our method and the estimates by IBGE (Brazil's Census Bureau) for the harvests from 2005 to 2016, were equal to 0.98. At state level the soybean, cotton, corn and sunflower areas had a correlation equal 0.98, 0.73, 0.96 and 0.80.
The following data sets are provided:
(a) The classified maps in compressed TIFF format (one per year) at MODIS resolution.
(b) A QGIS style file for displaying the data in the QGIS software
(c) An RDS file (R compressed format) with the training data set (2,115 ground samples).
The software used to produce the analysis is available as open source on https://github.com/e-sensing.
Note: The TIFF raster files use the Sinusoidal Projection, which is the same cartographical projection used by the input MODIS images. When opening the TIFF raster maps in QGIS, to ensure correct navigation please use the Sinusoidal Projection, by selecting in QGIS projection menu, the following option:
"Generated CRS (+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs)"