在r中的一个栅格中显示3个变量

Question

我发现可以使用二元等值线来同时显示两个栅格的值。例如： https://flowingdata.com/2015/03/13/bivariate-choropleth-how-to/

我想知道 r 是否有一个包来显示三个栅格的值？

在这种情况下，比例尺可以是一个 rgb 立方体，比如

http://www.dig.cs.gc.cuny.edu/manuals/Gimp2/Grokking-the-GIMP-v1.0/node50.html,

三变量栅格的比例是三个栅格中每一个的红色、绿色、蓝色的组合？

我一直在尝试生成立方体（结合所有 RGB 值）以将其用作比例。但不幸的是，文件变得非常大：

require(reshape2)
library(rgl)
Cols<-data.frame()

for (blue in 0:255)
{
rg <- outer(0:255, 0:255, rgb, blue=blue, maxColorValue = 255)
rg <-data.frame(blue,melt(rg))
colnames(rg)<-c("blue","red", "green","col")
Cols<-rbind(Cols,rg)
}

plot3d(Cols$blue, Cols$red, Cols$green, col=Cols$col)

Answer 1

我认为raster 包的plotRGB() 功能可以满足您的需求。请参阅此处以获取具有三个栅格层的模拟堆栈的示例：

set.seed(1002)

library(raster)

r1 <- raster(matrix(rnorm(10000), 100, 100))
r2 <- raster(matrix(rnorm(10000), 100, 100))
r3 <- raster(matrix(rnorm(10000), 100, 100))

s <- stack(r1, r2, r3)

plotRGB(s, stretch='hist')