7

事实上,这个问题由两个针对相同行为的问题组成。

  1. 如何将文本(因每个面板而异)添加到面板区域的固定位置?我知道panel.textlatticeExtra::layer 解决方案,但它使用绘图区域坐标添加文本。例如,我想在每个面板的右下角添加文本,即使它们的比例不同。

  2. 如何在 levelplot 面板区域中添加文本?这里解释的方法要求 levelplot 有一个plot_01.legend.top.vp区域来添加我没有的文本,并且该trellis对象是之前绘制的。此外,我想ylab在下图中的左侧添加文本。我ylab在这里用来说明行的含义,但我需要第二个表示 y 轴值的 ylab。我为这个问题找到了另一个 问题,但它不起作用。

样本图

上面的图是由raster::stack对象和rasterVis::levelplot方法创建的。即使我更喜欢优雅的解决方案,我也同意一个肮脏的解决方案。此外,尽管有上述问题,我对其他使用levelplot.

4

1 回答 1

2

目前正在R-sig-Geo上讨论一个非常相似的问题,看看我在那里提供的解决方案。这是相应的示例代码,可让您使用trellis.focus(..., clip.off = TRUE)from lattice在格状图的面板区域内部或外部添加自定义文本注释。

library(rasterVis)
library(grid)

## sample data
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
s <- stack(r, r+500, r-500, r+200)

p <- levelplot(s, layout = c(2, 2), names.att = rep("", 4), 
               scales = list(y = list(rot = 90)))

## labels
cls <- c("col1", "col2")
rws <- c("row1", "row2")

png("~/rasterVis.png", width = 14, height = 16, units = "cm", res = 300L)
grid.newpage()
print(p, newpage = FALSE)

## loop over panels to be labelled (ie 1:3)
panels  = trellis.currentLayout()
for (i in 1:3) {

  # focus on current panel of interest and disable clipping
  ids <- which(panels == i, arr.ind = TRUE)
  trellis.focus("panel", ids[2], ids[1], clip.off = TRUE)

  # add labels
  if (i %in% c(1, 3)) {
    if (i == 1) {
      grid.text(cls[1], x = .5, y = 1.1)            # add 'col1'
      grid.text(rws[1], x = -.35, y = .5, rot = 90) # add 'row1'
    } else {
      grid.text(rws[2], x = -.35, y = .5, rot = 90) # add 'row2'
    }
  } else {
    grid.text(cls[2], x = .5, y = 1.1)              # add 'col2'
  }

  trellis.unfocus()
}

dev.off()

光栅可见

您可以在此处找到更多信息:

于 2017-07-21T08:34:29.370 回答