6

我正在使用 RasterVis 和 levelplot 制作一些栅格的格状图。我目前对大多数事情都很好,但我想将每个面板的标题从文件名更改为选定的字符串(文件名很复杂而且很长,我只想使用一年,例如“2004”)。

查看levelplot 页面,这表明 levelplot 会根据参数“useRaster”寻找一些设置,它要么转到panel.levelplot 要么 panel.levelplot.raster,但我正在努力使用后面的这些功能。

非常感谢任何帮助,这里有一些示例代码;

require(rasterVis)

layers <- c(1:4)
s2 <- stack()

for (i in layers) {
  r <- raster(nrows=100, ncols=100,ext)
  r[] <- sample(seq(from = 1, to = 6, by = 1), size = 10000, replace = TRUE)
  rasc <- ratify(r)
  rat <- levels(rasc)[[1]]
  rat$legend <- c("A","B","C","D","E","F")
  levels(rasc) <- rat
  s2 <- stack(s2, rasc)
}

levelplot(s2, col.regions=rev(terrain.colors(6)),main = "example")

在上面例如,我希望“layer.1.1”为“2004”,依此类推到 2007

4

1 回答 1

8
require(rasterVis)

layers <- c(1:4)
s2 <- stack()

for (i in layers) {
  r <- raster(nrows=100, ncols=100)
  r[] <- sample(seq(from = 1, to = 6, by = 1), size = 10000, replace = TRUE)
  rasc <- ratify(r)
  rat <- levels(rasc)[[1]]
  rat$legend <- c("A","B","C","D","E","F")
  levels(rasc) <- rat
  s2 <- stack(s2, rasc)
}
levelplot(s2, col.regions=rev(terrain.colors(6)),main = "example", names.attr=2004:2007)

在此处输入图像描述

p.strip <- list(cex=1.5, lines=1, col="blue", fontfamily='Serif')

levelplot(s2, col.regions=rev(terrain.colors(6)), main = "example",
          names.attr=2004:2007, par.strip.text=p.strip)

在此处输入图像描述

于 2016-01-07T10:24:03.717 回答