ArcGIS就像R for Rasters中的拉伸颜色方案一样

Question

我正在尝试复制ArcGIS用于在R中显示栅格的拉伸颜色方案。我想使用将颜色标度的最小值/最大值定义为平均值的方法（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)