plot() сдвигает растр вправо

[rhot]

Пытаюсь отрисовать одновременно 4 растра с разной цветовой палитрой. Функция image() работает отлично, но не умеет показывать условные обозначения. А вот plot() умеет, но 4-ый растр сдвигает почему-то вправо. См. приложение.

Код: Выделить всё

par(mfrow=c(2,1))
plot(r.lc.aspen,col=pal.a, legend=TRUE)
plot(r.lc.birch,col=pal.b,add=TRUE,legend=FALSE)
plot(r.lc.spruce,col=pal.s,add=TRUE, legend=FALSE)
plot(r.lc.pine,col=pal.p,add=TRUE, legend=FALSE)
title("plot()")
image(out[[1]], col=pal.s)
image(out[[2]], col=pal.p, add=TRUE)
image(out[[3]], col=pal.b, add=TRUE)
image(out[[4]], col=pal.a, add=TRUE)
title("image()")

Является ли это ествественным ограничением функции plot? Если использовать image(), есть ли решение проще, чем создавать свою легенду отдельно?

[bolotoved]

Ощущение, что у 4-го растра ругая проекция (зона UTM?), а plot при этом не может перепроецировать на лету.

[gamm]

не очень понятно, что у вас за данные, поскольку на самом деле plot() - это только общее имя, а вызывается что-то, соответствующее типу данных, как и, скорее всего, для image().

а чтобы выводить макеты со всеми чудесами, придуман sp::spplot(), у которого есть sp.layout, там можно любую композицию создать. Если пару часиков попотеть

P.S. Насколько я понял, нужно вывести цветовую легенду по отдельным значениям, так это делается из одного растра, зачем 4?

[rhot]

Ощущение, что у 4-го растра ругая проекция (зона UTM?), а plot при этом не может перепроецировать на лету.

Проекции идентичны. Пробовал в другом порядке отрисовывать, то же самое выходит - последний растр слетает.

не очень понятно, что у вас за данные

У меня четыре растра, которые содержат возраст лесных насаждений (ель, сосна, берёза и осина). Для каждой породы - своя легенда (цветвой градиент).

Код: Выделить всё

> class(r.lc.aspen)
[1] "RasterLayer"
attr(,"package")
[1] "raster"

[gamm]

не очень понятно, что у вас за данные

У меня четыре растра, которые содержат возраст лесных насаждений (ель, сосна, берёза и осина).

яснее не стало. Покажите вывод str() для одного из растров, а лучше для всех - сохраните в текстовый файл с помощью sink() и приложите.

Для каждой породы - своя легенда (цветвой градиент).

что мешает сделать из них один RGB растр, и его нарисовать? Или они устроены по разному?