1

我需要您对如何可视化/映射跨越负值和正值的数据的见解。我有 6 个已堆叠在一起的栅格 ( s)，随后将通过levelplotR 中的函数进行绘制。下面是每个栅格的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))

注意：我需要一个所有栅格通用的颜色键。

问题：对于范围较小的栅格，单元格之间的变化很难传达 - 地图上的信息不能很好地可视化。如何使用levelplot如此大的整体范围来可视化此类数据？

您可以生成 6 个值在上述范围内的栅格，以便为我提供更多线索。

谢谢你的帮助！colorkey 不需要从红到蓝。像“默认彩虹调色板”这样的东西就可以了。

在此处输入图像描述"

对于所有这些栅格，我需要一个共享的颜色键。目前，像元范围较小的栅格显示的变化非常小。

4

1 回答 1

0

在我看来，你有两个选择：

  1. 为每个图形使用不同的键。
  2. 重新缩放数据以创建通用比例。

因为您需要一个颜色键，所以您应该重新调整数据。尽管本文档是关于时间序列的，但它对这个问题给出了很好的建议。

于 2015-02-28T06:58:56.053 回答