我想在分类栅格图中添加一些箭头。以Base R(在未经批准的栅格上)为例:
r <- raster(nrow=10, ncol=10, crs='+proj=utm +zone=1')
r[] = 1
r[51:100] = 3
r[3:6, 1:5] = 5
plot(r)
arrows(100, 50,
0, 50,
length = 0.1,
lwd=2,
col="red"
)
在levelplot中,您必须首先批准栅格并添加适当的级别:
r <- raster(nrow=10, ncol=10, crs='+proj=utm +zone=1')
r[] = 1
r[51:100] = 3
r[3:6, 1:5] = 5
r <- ratify(r)
rat <- levels(r)[[1]]
rat$landcover <- c('Pine', 'Oak', 'Meadow')
rat$class <- c('A1', 'B2', 'C3')
levels(r) <- rat
levelplot(r)
如何将红色箭头(如示例1中所示)添加到第二段代码中显示的levelplot()
上?
谢谢。
答案 0 :(得分:1)
您可以使用layer
和panel.arrows
来解决它:
levelplot(r) +
layer(panel.arrows(100, 50, 0, 50,
col = 'white'))