Question

我有4个rasters，其值包含正负标度。对于4 rasters： 1）我想在colorkey的两端只使用一个颜色中断来表示+/- 20范围之外的所有值。这是因为我对较小值的光栅内和光栅间可变性更感兴趣。

2）如何指定colorkey的标签，使其看起来类似于：at=c(<=-20, -10, 0 , 10 , >=20)。这里最重要的是如何在colorkey标签

中包含不等号

可重复的示例和示例代码如下所示。该图像使用levelplot函数显示我的实际数据显示。

library(raster)
library (colorRamps)
set.seed(100)
ras <- raster(ncol=100, nrow=100)
ras1 <- setValues(ras, (1:ncell(ras))/100 + rnorm(ncell(ras)) - 50)

s=stack(ras1,ras1,ras1,ras1)


color_levels=14 #the number of colors to use
max_abolute_value=max(abs(c(cellStats(s, min), cellStats(s, max)))) #what is the maximum absolute value of raster?
color_sequence=unique(round(seq(-max_abolute_value,max_abolute_value,length.out=color_levels+1),0))

myColorkey <- list(at=color_sequence,space = "bottom", ## where the colors change 
                   labels=list(axis.line = list(col = NA),at=color_sequence,rot=0,cex=0.9,font=6,
                               fontface=1),height=1,width=1.4)

col1 <- colorRampPalette(c("darkred", "red3","red",
                           "gray96",  "lightskyblue", "royalblue3", 
                           "darkblue"))

levelplot(s,contour=F, layout=c(4, 3),  col.regions = col1,colorkey=myColorkey,margin=FALSE,xlab=NULL,ylab=NULL,par.strip.text=list(cex=0))

Answer 1

对于（i），您可以简单地reclassify栅格图层并将所有值设置为小于-20（或大于20），例如到-20.5（或20.5）。我想这很有意义，因为无论如何你对较小或较大的值都不感兴趣。至于（ii），我建议使用spplot并修改?levelplot中详细说明的颜色键（来自包网格）。请注意，可以使用Unicode字符轻松插入所需的标签（例如greater than or equal to）。

 ## reclassify data
m_rcl <- matrix(c(-100, -20, -20.5, 
                  20, 100, 20.5), 
                byrow = TRUE, ncol = 3)

s_rcl <- reclassify(s, m_rcl)

## colors
library(RColorBrewer)
cols <- brewer.pal(6, "RdBu")
cols <- rev(cols)
cols <- colorRampPalette(cols)

## display data
spplot(s_rcl, col.regions = cols(100), at = seq(-21, 21, 1),
       colorkey = list(space = "bottom", 
                       labels = list(at = seq(-20, 20, 10), 
                                     labels = c("\u2264 -20", 10, 0, 10, "\u2265 20"))))

Answer 2

按照@fdetsch的建议，这对我有用：

m_rcl <- matrix(c(-375, -100, -100.5, 
                  100, 484, 100.5), 
                byrow = TRUE, ncol = 3)

s_rcl <- reclassify(s, m_rcl)


levelplot(s_rcl,contour=F,margin=FALSE,xlab=NULL,ylab=NULL,par.strip.text=list(cex=0), scales = list(x=x.scale, y=y.scale),
          col.regions = col1, at = seq(-110, 110,20),layout=c(4, 3),index.cond=list(c( 1,2,3,4,5,6,7,8,9,10,11,12)),
          colorkey = list(space = "bottom", 
                          labels = list(at = seq(-100, 100, 20), rot=0,cex=0.9,font=6,fontface=1,
                                        labels = c("\u2264 -100", "-80",  "-60" , "-40"  ,"-20" ,  
                                                   "0" ,  "20" ,  "40" ,  "60" ,  "80" ,  "\u2265 100")),height=1,width=1.4))

非常感谢。

对颜色范围使用单一颜色，并在colorkey levelplot R

2 个答案: