r - ArcGIS 类似于 R for Rasters 中的拉伸(stretch)配色方案

Question

我正在尝试复制 ArcGIS 用于在 R 中显示栅格的拉伸(stretch)配色方案。我想使用将色标的最小/最大值定义为 mean(raster_values) +/- 2*sd( raster_values)，然后超出此范围的任何颜色都将根据不等式设置为最小值或最大值的颜色。

我不知道如何告诉 R 将高于/低于阈值的所有值绘制为最大/最小颜色值(它只是不会为单元格着色)。

#Load libraries
  library(raster)
  library(RColorBrewer)

#Color scale
  bluescale = colorRampPalette(brewer.pal(9,"PuBu"))

#Create raster and plot it 
  x = raster(matrix(rnorm(625,20,5),25,25))
  plot(x,col=bluescale(50),breaks=seq(mean(values(x))-2*sd(values(x)),mean(values(x))+2*sd(values(x)),length.out=51))

它不会绘制所有值。我怎样才能做到这一点？

Answer 1

延长最后一组的休息时间如何

br <- seq(mean(values(x))-2*sd(values(x)),
    mean(values(x))+2*sd(values(x)),length.out=51)
br[1] <- min(values(x));
br[length(br)] <- max(values(x))
plot(x,col=bluescale(50),breaks=br)