Depois de estudar por um tempo como criar um gráfico a partir de um objeto raster, obtive quase tudo o que precisava, exceto para adicionar um título à legenda do spplot.
O que eu tentei:
- usar
mtext("XXX (m)", side=4)
, - para adicionar um
sp.text
fora da grade, - usar
legend.args = list(title="XXX (m)")
.
A imagem abaixo identifica o recurso que eu adicionei manualmente ao meu gráfico e no qual gostaria de incluir o uso de R.
Aqui está um exemplo do que recebi até agora:
require(raster)
require(sp)
require(lattice)
north <- list("SpatialPolygonsRescale", layout.north.arrow(type=1), offset = c(0.95,0.85), scale=0.1)
scale <- list("SpatialPolygonsRescale", layout.scale.bar(),
offset = c(0.55, 0.03), scale = 0.4, fill = c("transparent","black"))
txt1 <- list("sp.text", c(0.55, 0.08), "0")
txt2 <- list("sp.text", c(0.75, 0.08), "0.2")
txt3 <- list("sp.text", c(0.95, 0.08), "0.4")
raster_layout <- list(north,scale,txt1,txt2,txt3)
cuts <- c(110,120,130,140,150,160,170,180,190,200)
spplot(raster(volcano), scales = list(draw=T), at=cuts, col.regions = palette(gray(seq(0,0.9,len=9))), sp.layout= raster_layout, add=T)
Obrigado pela ajuda.