如何在 R 中为不同的 spplots 使用相同的颜色间隔

Question

我有不同日子的温度数据插值图，比如：

我想并排显示不同日子的地图，但希望两者在传奇中具有相同的颜色渐变和间隔。 最高温度值为：29.8 分钟：-13.4

我确实尝试使用来自How to use the same color scale for points in different maps in R? 的答案，但这样做确实出错了：

maxim <- 30
minim <- -14
cuts <- 2*(-7:15)
p1 <- spplot(idw_ger, "var1.pred", main="IDW January, 10", cuts=cuts, col.regions=bpy.colors(30)) +
latticeExtra::layer(sp.polygons(ger, fill = NA, col = 'blue')) +
latticeExtra::layer(sp.points(tdata,
                            fill = NA, col = 'black'))

seq.default(zrng1, zrng[2], length.out = cut + 2) 中的错误： 长度必须是非负数 另外：警告信息： 在 seq.default(zrng1, zrng[2], length.out = cut + 2) 中： 'length.out' 参数使用的第一个元素

我做错了什么，我该如何让它发挥作用？

Answer 1

鉴于只能为 SpatialPointsDataFrame 指定精确的cuts，我认为您需要为spplot 使用替代方案。这是一个使用image 和image.plot 作为图例的示例（改编自https://stat.ethz.ch/pipermail/r-sig-geo/2007-June/002143.html）。通过设置数据的上限和下限（zlim），您可以创建两个相同的色带。只需将zinc.idw 和copper.idw 替换为您自己的数据即可。

加载数据：

library(sp)
library(gstat)
library(fields)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
coordinates(meuse.grid) = ~x+y
gridded(meuse.grid) = TRUE

创建两个 IDW 网格：

zinc.idw = krige(zinc~1, meuse, meuse.grid)
copper.idw = krige(copper~1, meuse, meuse.grid)

设置情节：

par(mfrow=c(1,2),oma=c(2,2,2,2))
cols = bpy.colors(30)
zmin = 0
zmax = 1000

绘制地图和测量位置：

image(zinc.idw, zlim=c(zmin,zmax),col=cols)
box()
mtext(side=3,text='zinc')
points(meuse, pch = 3)

添加图例：

image.plot( zlim=c(zmin,zmax), nlevel=5,legend.only=TRUE, horizontal=F, col=cols)

其他铜网也一样；

image(copper.idw, zlim=c(zmin,zmax),col=cols)
box()
mtext(side=3,text='copper')
points(meuse, pch = 3)
image.plot( zlim=c(zmin,zmax), nlevel=5,legend.only=TRUE, horizontal=F, col=cols)