更改區域間隔rasterVis horizo​​nplot

Question

我正在與rasterVishorizonplot功能，我想改變區間間隔;例如我想要顯示x axis中每10個緯度區域的平均值。這是an example從rasterVis庫。

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不能存儲字符）的標籤，所以每個面板的名稱沒有直接關係到該區域。