ArcGIS like R的拉伸配色方案

Question

我試圖複製ArcGIS用來在R中顯示柵格的拉伸配色方案。我想使用定義顏色比例的最小值/最大值作爲平均值的方法（ raster_values）+/- 2 * sd（raster_values），然後這個範圍之外的任何東西都會被設置爲min或max的顏色，這取決於不等式。

我不知道如何告訴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)