两个栅格之间的散点图，从第三个栅格给出颜色

Question

我正在绘制两个栅格数据，生成下面的图像。

我想用取自第三个栅格数据的变量（具有相同的 bbox、像素大小等）为图中的每个点着色。 来自 R 用户的任何想法？这个操作在绘制数据集中的数据时非常容易，但是我不知道栅格……

我在这里附上生成图像的代码（简化，我认为您不需要所有绘图参数，例如 abline、xlab 等）：

plot(mask(raster1, my_mask,maskvalue=0), #first raster, masked by my_mask
      mask(raster2, my_mask,maskvalue=0),  #second raster, masked by my_mask           
      col = alpha('black', 0.1), #the current color scheme     
      )
 raster3 #raster with categorical variable, 
         #that should give the colors to the points in the graph

非常感谢！

Answer 1

使用中定义的xyplot 方法 rasterVis 您可以使用 RasterStack 的层，就好像它们是 a 的列一样 data.frame。因此，它们可以是公式的组成部分或 groups 参数。

例如，

library(raster)
library(rasterVis)

f <- system.file("external/test.grd", package="raster")
r <- raster(f)
r2 <- r + 500 * init(r, rnorm)
## categorical variable
r3 <- cut(r, 3)

s <- stack(r, r2, r3)
names(s) <- c('r', 'r2', 'r3')

xyplot(r ~ r2, groups = r3, data = s,
       auto.key = list(space = 'right'),
       alpha = 1)