2016-11-10 26 views
3

例えば、ラスタ画像をプロットしながら:リバース伝説

library(raster) 

r <- raster(nrow = 3, ncol = 3) 
values(r) <- 1:9 

plot(r, col = terrain.colors(255)) 

enter image description here

私は9(下)に伝説1(上)から、すなわち、昇順であることを取得できますか?

私はlegend.argsを考えましたが、正しい引数が見つかりませんでした。

答えて

0

私はちょっと試しましたが、私が解決策を見つけたと思うのは、それが最もエレガントな方法ではありませんが。

library(raster) 

r <- raster(nrow = 3, ncol = 3) 
values(r) <- 1:9 

par(mar = c(3, 3, 4, 3)) 
plot(r, col = terrain.colors(255),legend = FALSE, xlim = c(-200,200), 
    ylim = c(-200,200)) 

vz = matrix(1:100, nrow = 1) 
vx = c(185, 195) 
vy = seq(-10, 10, length.out = 100) 

par(new = TRUE, mar = c(3, 3, 4, 3)) 
plot(1, xlab = "", ylab = "", axes = FALSE, type = "n", 
    xlim = c(-200, 180), ylim = c(-20, 20)) 
image(vx, vy, vz, col = rev(terrain.colors(255)), axes = FALSE, 
    xlab = "", ylab = "", add = TRUE) 
polygon(c(185, 195, 195, 185), c(-10, -10, 10, 10)) 
axis(4, at = seq(-10, 10, length.out = 9), labels = 9:1, las = 1) 

enter image description here

とにかく、私は他のアイデアをいただければと思います!