r - Levelplot 颜色键 - 范围和极值

Question

是否可以在 R 中创建一个像下面这样的颜色键？（这个来自软件Grid Analysis and Display System - Grads）。

有两个我无法在 R 中重现的功能：

1. 该序列是非线性的，但它显示为
2. 大于 200 的值为灰色/小于 0 的值为白色

我正在使用来自 rastervis 的 levelplot，它使用 lattice levelplot 绘制栅格：

require(raster)
require(rasterVis)

set.seed(200)

X = seq(-40,0,by=1)
Y = seq(-60,-40,by=1)
grid = expand.grid(list(X=X,Y=Y))
Z = rnorm(nrow(grid),mean=10,sd=100)

data = data.frame(grid,Z)
r = rasterFromXYZ(data)
mapTheme <- rasterTheme(region=c('#EEF7FA','#D6F8F7',"#BEDAFF",'#5DA4FF',"#0000FF","#D4F9E2","#00FF7F","#008000","#FFFF00","#FFD27F", "#FFB732"    ,"#EE7600",
                                 "#D53E4F","#FF6A6A"))  
my.at = c(0,1,5,10,15,20,25,30,40,50,75,100,150,200)
myColorkey <- list(at=my.at,
                   space="bottom",
                   labels=list(at=my.at))
p=levelplot(r, par.settings=mapTheme,at = my.at, colorkey=myColorkey,margin=F)
print(p)

结果：

如您所见，小于 0 和大于 200 的值都是白色的，我不知道如何将大于或小于某个值的值设置为特定颜色。Morover，我怎样才能使颜色键中连续粗标记之间的空间具有相同的大小，尽管间隔不一样？

Answer 1

这是非线性序列等大小间隔的解决方法：

library(raster)
library(rasterVis)

set.seed(200)
X = seq(-40,0,by=1)
Y = seq(-60,-40,by=1)
grid = expand.grid(list(X=X,Y=Y))
Z = rnorm(nrow(grid),mean=10,sd=100)

data = data.frame(grid,Z)
r = rasterFromXYZ(data)
mapTheme <- rasterTheme(region=c('#EEF7FA','#D6F8F7',"#BEDAFF",'#5DA4FF',"#0000FF","#D4F9E2","#00FF7F",
                                 "#008000","#FFFF00","#FFD27F", "#FFB732" ,"#EE7600", "#D53E4F","#FF6A6A"))  

my.at=c(0,1,5,10,15,20,25,30,40,50,75,100,150,200)
my.brks=seq(0, 200, by=15)

myColorkey <- list(at=my.brks, labels=list(at=my.brks, labels=my.at), space="bottom")
p=levelplot(r, par.settings=mapTheme, at=my.at, colorkey=myColorkey, margin=F)
print(p)

这可能是值小于 0 且大于 200 的解决方案：

library(raster)
library(rasterVis)

set.seed(200)
X = seq(-40,0,by=1)
Y = seq(-60,-40,by=1)
grid = expand.grid(list(X=X,Y=Y))
Z = rnorm(nrow(grid),mean=10,sd=100)

data = data.frame(grid,Z)
r = rasterFromXYZ(data)
mapTheme <- rasterTheme(region=c('white','#EEF7FA','#D6F8F7',"#BEDAFF",'#5DA4FF',"#0000FF","#D4F9E2","#00FF7F",
                                 "#008000","#FFFF00","#FFD27F", "#FFB732" ,"#EE7600", "#D53E4F","#FF6A6A", "gray"))  
max(values(r))
min(values(r))

my.at=c(min(values(r)), 0,1,5,10,15,20,25,30,40,50,75,100,150,200, max(values(r)))
my.brks=seq(0, 200, by=13)

myColorkey <- list(at=my.brks, labels=list(at=my.brks, labels=c(-276,0,1,5,10,15,20,25,30,40,50,75,100,150,200, 388)), space="bottom")
p=levelplot(r, par.settings=mapTheme, at=my.at, colorkey=myColorkey, margin=F)
print(p)

你的颜色没有从浅到深。您可以使用 RColorBrewer 包来解决此问题。

library(RColorBrewer)
reds = brewer.pal(5, "YlOrRd")
greens = brewer.pal(3, "Greens")
blues = brewer.pal(5, "Blues")
mapTheme <- rasterTheme(region=c('white', blues, greens, reds, "gray"))

Answer 2

这是一个非常有用的解决方法。虽然没有解决问题 1，但我发现对问题 2 有用的东西（为低于/高于颜色条范围限制的值添加三角形）可以通过添加以下内容来实现：

library(s2dverification)
data_array <- array(Z, dim = c(length(X), length(Y)))
PlotEquiMap(data_array, X, Y,bar_limits=c(0,200),col_inf='white',col_sup='gray')

带颜色条的光栅

Answer 3

另一个对 lattice 进行更新的解决方案：

library(raster)
library(rasterVis)

set.seed(200)
X = seq(-40,0,by=1)
Y = seq(-60,-40,by=1)
grid = expand.grid(list(X=X,Y=Y))
Z = rnorm(nrow(grid),mean=10,sd=100)
data = data.frame(grid,Z)
r = rasterFromXYZ(data)

levelplot(r, margin=F, at=c(-Inf, seq(0, 200, 20), Inf),
          colorkey = list(tri.lower = TRUE, tri.upper = TRUE))

只要在at定义中添加“-Inf”和“Inf”，就会激活将三角形添加到颜色条的选项。