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

Question

我正在使用 rasterVis horizonplot 函数，我想更改区域间隔；例如，我想显示 x axis 中每 10 个纬度区域的平均值。这是来自 rasterVis 图书馆的 an example。

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

Answer 1

区域由 zonal 函数确定，可通过 digits 参数（整数，要保留的位数）间接访问。使用这种方法，您无法完全控制间隔数。例如：

horizonplot(SSTanom, digits = -1)

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

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

遗憾的是，cut没有保留区间的标签（因为一个Raster不能存储字符），所以每个面板的名称与区域没有直接关系。

Change the zone interval in rasterVis horizonplot