我正在尝试复制 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))
它不会绘制所有值。我怎样才能做到这一点?
最佳答案
延长最后一组的休息时间如何
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)
关于r - ArcGIS 类似于 R for Rasters 中的拉伸(stretch)配色方案,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/27788184/