0

我正在使用rasterVis horizonplot功能,我想更改区域间隔;例如,我想显示. 中每 10 个纬度区域的平均值x axis。这是图书馆的一个例子。rasterVis

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

1 回答 1

2

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

horizonplot(SSTanom, digits = -1)

水平线图和数字

另一种解决方案是提供更多灵活性的dirXY论点。horizonplot例如:

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

水平线图和切割

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

于 2016-06-28T19:38:51.520 回答