ArcGIS 类似于 R for Rasters 中的拉伸配色方案

Question

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

我不知道如何告诉 R 绘制所有值 above/below 阈值作为 max/min 颜色值（它不会为单元格着色）。

#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)

ArcGIS 类似于 R for Rasters 中的拉伸配色方案

ArcGIS like stretched color scheme in R for Rasters

plot

r

raster

arcgis