使用R中的levelplot可視化負/正數據

Question

我需要您對如何可視化/映射跨越負值和正值的數據的見解。 我有6個柵格已經堆疊在一起（ s ）和隨后將被繪制經由levelplot在R.函數以下是在max和min的值對於每個raster 。

39.2887, 53.09207  (min, max) # r1
-32.4956, -27.25534  (min, max)# r2   
-14.37683, -11.37742  (min, max)# r3    
9.512934, 13.60197  (min, max)# r4    
-4.993901, -1.851784  (min, max)# r5    
-8.190711, -5.104764  (min, max)# r6

目前，我可以通過以下方式制作地圖：

library(raster)
library(rasterVis)
library(colorRamp)
    s <- stack(r1,r2,r3,r4,r5,r6)
    themes2 <- colorRampPalette(c("darkred", "red3", "orange", "yellow", "lightskyblue", "royalblue3", "darkblue"))(19)
    myat =unique(seq(floor(min(s)) ,ceiling(max(s)),length.out=20))
    myat=round(myat,digits = 0)#     
    #themes <- rasterTheme(region=rev(brewer.pal(11,'RdYlBu')))
    myColorkey <- list(at=myat,space = "right",labels=list(cex=1,at=myat))

    if (dev.cur() == 1) x11(width=18,height=18)

    levelplot(s, layout=c(3, 2), index.cond=list(c(1, 3, 5, 2, 4, 6)),col.regions=themes2,  
              margin=FALSE,xlab=NULL,at =unique(seq(floor(min(s)) ,ceiling(max(s)),length.out=20)),
              par.strip.text=list(cex=0),colorkey=myColorkey,scales=list(alternating=F))

注意：我需要所有柵格通用的單個colorkey。

問題：對於范圍較小的柵格，像元之間的變化難以傳達-地圖上的信息無法很好地可視化。 如何在如此大的范圍內使用levelplot可視化此類數據？

您可以生成6個柵格，其值在上面顯示的范圍內，以給我進一步的線索。

謝謝你的幫助！ 色鍵不需要從紅色到藍色。 像“默認彩虹調色板”之類的東西就可以。

“

對於所有這些柵格，我都需要一個共享的色鍵。 目前，像元范圍較小的柵格顯示出很小的變化。

Answer 1

我認為，您有兩種選擇：

1. 為每個圖形使用不同的鍵。
2. 重新縮放數據以創建通用比例。

因為只需要一個色鍵，所以應該重新縮放數據。 盡管此文檔是有關時間序列的，但它提供了有關此問題的良好建議。