更改rasterVis水平图中的区域间隔

Question

我正在使用rasterVis horizonplot函数，我想改变区域间隔；例如，我想显示x axis中每10个纬度区域的平均值。这是来自一个例子库的rasterVis。

library(raster)
library(rasterVis)
horizonplot(SSTanom, col.regions = rev(brewer.pal(n = 10, 'RdBu')))

Answer 1

区域由zonal函数确定，该函数可以通过digits参数间接访问(整数，要维护的数字数)。使用这种方法，您无法完全控制间隔的数量。例如：

horizonplot(SSTanom, digits = -1)

​

​

另一个解决方案是dirXY的horizonplot参数，它提供了更多的灵活性。例如：

horizonplot(SSTanom, dirXY = cut(y, 10))

​

​

不幸的是，cut没有保留间隔的标签(因为Raster不能存储字符)，因此每个面板的名称与区域没有直接关系。