Get MapBiomas Peru land use / land cover raster
Source:R/get_mapbiomas_peru_lulc.R
get_mapbiomas_peru_lulc.RdLazily reads a single-year Land Use and Land Cover (LULC) classification
raster from the MapBiomas Peru collection, hosted as a Cloud Optimized
GeoTIFF (COG) on Google Cloud Storage. Only the bytes required for the
requested extent are downloaded (via GDAL's /vsicurl/ driver).
Optionally crops and masks the raster to an area of interest.
Arguments
- year
Integer. Year of the classification (e.g.
2024).- crop_to
Optional. An
sf/sfcobject,SpatVector, orSpatExtentdefining the area of interest. IfNULL(default), the full raster for Peru is returned.- collection
Integer. MapBiomas Peru collection number. Currently available collections are
1,2, and3. Default is3(the latest).
Examples
# \donttest{
library(geoidep)
# Download departmental boundaries
lima <- get_departaments("LIMA")
#>
|
| | 0%
|
|= | 1%
|
|== | 2%
|
|== | 3%
|
|=== | 4%
|
|===== | 7%
|
|=========== | 16%
|
|==================== | 28%
|
|====================== | 32%
|
|======================= | 33%
|
|================================ | 46%
|
|======================================== | 58%
|
|============================================= | 65%
|
|====================================================== | 77%
|
|=============================================================== | 89%
|
|======================================================================| 100%
# Download and crop MapBiomas Peru LULC for 2024
lulc_2024 <- get_mapbiomas_peru_lulc(year = 2024, crop_to = lima)
lulc_2024
#> class : SpatRaster
#> size : 11315, 8827, 1 (nrow, ncol, nlyr)
#> resolution : 0.0002694946, 0.0002694946 (x, y)
#> extent : -77.88636, -75.50753, -13.32354, -10.27421 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84 (EPSG:4326)
#> source(s) : memory
#> varname : peru_collection3_integration_v1-classification_2024
#> name : classification_2024
#> min value : 0
#> max value : 72
#> class : SpatRaster
#> dimensions : 4521, 5103, 1 (nrow, ncol, nlyr)
#> resolution : 0.00025, 0.00025 (x, y)
#> extent : -77.45, -75.7, -13.5, -10.4 (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84
#> source : peru_collection3_integration_v1-classification_2024.tif
#> name : classification_2024
# }