Exemplo 1: Mapa de cobertura do solo (Mapbiomas)

#Explorando dados geoespaciais em R

#Este exemplo mostra como:
#1) ler um arquivo raster 
#2) processar dados em variaveis raster
#3) gerar mapa de variaveis raster

#IMPORTANTE: use Session no menu acima para definir o diretorio de trabalho (set working directory) se necessario.

#carrega bibliotecas
library(raster)
library(rasterVis)

#le o arquivo de dados de biomas em 1985 (mapbiomas.org)
#arquivo contem 31 tipos de biomas

bio1985=raster("biomas_1985_rs.tif")

#metadados (informacoes sobre o arquivo de dados)
bio1985

## class      : RasterLayer 
## dimensions : 4694, 4708, 22099352  (nrow, ncol, ncell)
## resolution : 0.0083334, 0.0083334  (x, y)
## extent     : -74.00941, -34.77576, -33.8047, 5.312277  (xmin, xmax, ymin, ymax)
## crs        : +proj=longlat +datum=WGS84 +no_defs 
## source     : biomas_1985_rs.tif 
## names      : biomas_1985_rs 
## values     : 0, 33  (min, max)

#Arquivo contem 33 tipos de cobertura de solo
#Aqui a funcao reclassify eh usada para simplificar os dados
#para 3 caterogias: natural (cat 0), alterado (cat 1) e agua (cat 2)
#e salva resultado como desf1985
m = c(-Inf,0,NA, 0.5,13,0, 14,21,1, 22,Inf,2)
rcm = matrix(m, ncol=3, byrow=T)
desf1985=reclassify(bio1985,rcm)

#As proximas linhas convertem a variable desf2015 para dados 
#categoricos e insere nomes as novas categorias
desf1985 = ratify(desf1985)
rat = levels(desf1985)[[1]]
rat$landcover = c("Natural","Alterado","Agua")
levels(desf1985) = rat

#Escolhe as cores para o grafico
#mais cores: http://www.stat.columbia.edu/~tzheng/files/Rcolor.pdf
cc=c("darkolivegreen4","bisque3","blue")

#Gera mapa da variable desf1985
p1 = levelplot(desf1985,margin=F, colorkey=T,
               par.settings=custom.theme(region=cc),
               scales = list(cex=1, tck=-0.5),
               xlab="",ylab="",main="Cobertura do solo em 1985")
#Plota o mapa
print(p1)