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

【问题标题】：ArcGIS like stretched color scheme in R for RastersArcGIS 类似于 R for Rasters 中的拉伸配色方案
【发布时间】：2015-03-03 12:02:26
【问题描述】：

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

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

【问题讨论】：

标签： r plot arcgis raster

【解决方案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)

【讨论】：