繪制柵格堆棧時使R中的顏色條變暗？

Question

我有一個 rasterstack，我試圖通過 spplot 在 R 中繪制：

rasterstack
class      : RasterStack 
dimensions : 2803, 5303, 14864309, 3  (nrow, ncol, ncell, nlayers)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : 60.85, 105.0417, 15.95833, 39.31667  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
names      :   layer.1,   layer.2,   layer.3 
min values : -1.819715, -1.800774, -1.835778 
max values :    6.5533,    6.5533,    6.5533 

我想要一個從紅色到藍色的色階，中間是白色。 我通過這段代碼實現了色標：

devtools::source_gist('306e4b7e69c87b1826db')
plot2<-diverge0(rasterstack, ramp='RdBu')

然而，當繪制空間圖時，顯示的顏色非常淺，在視覺上不容易看到。 我想制作一個更暗且視覺可見的色階。

可重現的柵格：

library(raster)
r1 <- r2 <- r3 <- raster(nrow=10, ncol=10)
values(r1) <- runif(ncell(r1))
values(r2) <- runif(ncell(r2))
values(r3) <- runif(ncell(r3))
s <- stack(r1, r2, r3)

請注意，可重現的柵格沒有相同的柵格堆棧最小值和最大值，如果有人可以指導我的實際rasterstack堆棧最小值和最大值，我將不勝感激。 （其次，我希望色階從紅白變為藍）

Answer 1

這是一個使用rasterVis和ggplot2來獲得所需色階的解決方案。 scale_fill_gradientn()和facet_wrap()確保您的比例適合堆棧中的所有柵格。

library(raster)
#> Lade nötiges Paket: sp
r1 <- r2 <- r3 <- raster(nrow=10, ncol=10)
values(r1) <- runif(ncell(r1),-.2,.2)        # I adjusted your stack min max
values(r2) <- runif(ncell(r2),-.3,.3)
values(r3) <- runif(ncell(r3),-.5,.5)
s <- stack(r1, r2, r3)

library(ggplot2)
library(rasterVis)
#> Lade nötiges Paket: lattice
#> Lade nötiges Paket: latticeExtra
#> Lade nötiges Paket: RColorBrewer
#> 
#> Attache Paket: 'latticeExtra'
#> The following object is masked from 'package:ggplot2':
#> 
#>     layer
cols <- c('#ff0000', '#ffffff', '#002bff')   # red, white, blue
gplot(s) + 
  geom_tile(aes(fill=value)) +
  scale_fill_gradientn(colors = cols, aesthetics = "fill") +
  coord_equal() +
  facet_wrap(~variable)

^{由reprex 包（v0.3.0）於 2020 年 3 月 28 日創建}

希望有幫助！

Answer 2

您的色階（-1 到 7）實際上是正確的，例如 -1 的暗度對應於 +1 的暗度。 這是不同規模的一個重要方面。 如果要在中心保持零，則必須將其擴展為常見的絕對極值（在您的情況下為 -7 到 7）。 經過一番挖掘，我發現了這個更優雅的解決方案。 零現在位於刻度的中心。

library(raster)
#> Lade nötiges Paket: sp
r1 <- r2 <- r3 <- raster(nrow=10, ncol=10)
values(r1) <- runif(ncell(r1),-5,5)        # I adjusted your stack min max to [-1;7]
values(r2) <- runif(ncell(r2),-1,4)
values(r3) <- runif(ncell(r3),0,7)
s <- stack(r1, r2, r3)
print(cellStats(s, range))
#>        layer.1    layer.2   layer.3
#> [1,] -4.999763 -0.9014062 0.2718563
#> [2,]  4.891400  3.9860373 6.9711417

cols <- c('#ff0000', '#ffffff', '#002bff')   # red, white, blue
c.pal <- colorRampPalette(cols)              # interpolate between the 3 colors

s.max <- max(abs(cellStats(s,stat = unique)))  # get absolute max of stack

# make 100 even breaks with 0 in center
library(lattice)
breaks <- do.breaks(c(-s.max,s.max), 100)

spplot(s, col.regions = c.pal, at = breaks,
       colorkey = list(col = c.pal, at = breaks))

^{由reprex 包（v0.3.0）於 2020 年 3 月 28 日創建}

希望這可以澄清事情！