Множество (rasterVis) levelplots

Question

Я довольно отчаянно пытаюсь настроить два уровня уровня одной растровой таблицы каждый на одном участке. Похоже, что rasterVis::levelplot не принимает параметр par(mfrow = c(...)) для разделения панели. Пример использования двух rasterstacks, которые должны быть расположены бок о бок в новом окне участка:

f <- system.file("external/test.grd", package="raster") 
r1 <- stack(raster(f),log(raster(f))) 
r2 <- stack(raster(f)-500,raster(f)+500) 

par(mfrow=c(2,2)) 
    levelplot(r1) 
    levelplot(r2) 

Здесь levelplot(r1) это затевается на полномасштабную окне, в то время как levelplot(r2), к сожалению, закрашивать levelplot(r1).

Я пытался играть вокруг, обернув вызов для levelplot с функцией print, которая принимает split, а также newpage = false варианты. К сожалению, я не понимаю, как правильно использовать split, поэтому я получаю только разочарование.

Я очень ценю вашу помощь, заранее спасибо

Andi

Answer 1

Для сетки -На графики, в том числе производства решетки (которая лежит в основе rasterVis "Построение графиков функций) сеткаExtra функция grid.arrange() делает +/- то же, что и par(mfcol=). Для базы R графика.

library(gridExtra) 
p1 <- levelplot(r1) 
p2 <- levelplot(r2) 
grid.arrange(p1, p2, ncol=2) 

Редактировать: Альтернативного решетка -специфическое решения использует split= аргумента print.trellis(), метод построения для решетчатых участков (Н.Т. Батист & Оскар Perpiñán). split= принимает вектор из четырех чисел. 3-й и 4-й элементы вектора дают количество столбцов и строк на дисплее, а его 1-й и 2-й элементы предоставляют столбцы и строки для объекта, который печатается.

library(gridExtra) 
p1 <- levelplot(r1) 
p2 <- levelplot(r2) 
print(p1, split=c(1,1,2,1), more=TRUE) 
print(p2, split=c(2,1,2,1))