Olá, pessoal. Eu baixei os dados .tif de “uso e cobertura da terra” que possui apenas 1 layer. Mas eu gostaria de ter acesso apenas aos dados de urbanização. Como eu faço?
Olá @nascimento.flucas ,
Os dados de cobertura e uso da terra do MapBiomas possuem apenas um layer/banda, no qual o valor de cada pixel representa uma classe de cobertura e uso da terra, conforme a legenda (https://mapbiomas-br-site.s3.amazonaws.com/downloads/_EN__Códigos_da_legenda_Coleção_7.pdf). No caso da classe de áreas urbanizadas o ID é 24. Se precisa apenas desta classe, sugiro reclassificar o mapa mantendo apenas a classe de interesse (função reclassify no QGIS ou ArcGIS ou a função reclassify do pacote raster no R).
Vou Postar o código aqui. Pode ser ser útil!!
===========================================================
march 2023
load packages
library(raster)
load data “uso e cobertura da terra”
mapbiomas ← raster(“data/brasil_coverage_2021.tif”)
mapbiomas
matrix for reclassify function
m ← c(-1, 23.9, 0,
23.91, 24.1, 1,
24.11, 70, 0)
rclmat ← matrix(m, ncol=3, byrow=TRUE)
reclassify (demora p caramba!)
urban_raster_mapbiomas_br ← reclassify(mapbiomas, rclmat)
save in GTiff (demora p caramba tbm!)
writeRaster(urban_raster_mapbiomas_br, “urban_raster_mapbiomas_br.tif”,
format = “GTiff”,
datatype = “INT1U”,
options = c(“COMPRESS=DEFLATE”, “TFW=YES”),
overwrite = T)
===========================================================
1 curtida