Изменение интервала зоны в растровом горизонте горизонтального экрана

Question

Я работаю с rasterVishorizonplot функцией, и я хочу изменить интервал зоны; например, я хочу показать среднее значение для каждых 10 широтных зон в x axis. Это 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) 

Другим решением является dirXY аргумент horizonplot, что дает большую гибкость. Например:

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

К сожалению, cut не сохраняет метку интервала (поскольку Raster не может хранить символы), так что название каждой панели не имеет прямого отношения к зоне.