我想把SpatialPolygonsDataFrame画成一个半透明的主要物体(右边有图例),但我想要绘制更多的光栅(山坡)作为背景--只是为了制作更好的地图。我需要这样的东西:
spplot(polygons, sp.layout = list(list("raster", myRaster)))
但是看看?spplot
,在sp.layout
中指定光栅似乎是不可能的。我不能将光栅指定为主对象,因为主对象是polygons
dataFrame,它的值比例尺要在右侧的图例中绘制。
怎么可能在斯派特中画一个额外的光栅呢?
发布于 2014-11-21 05:59:59
这里有一种方法。也许有一个更整洁的方法来实现它,而不需要绘制两次多边形对象,尽管.
library(sp)
library(rasterVis)
r <- raster(nrow=18, ncol=36)
r[] <- runif(ncell(r)) * 10
r[r > 8] <- NA
pol <- rasterToPolygons(r, function(x) x > 6)
spplot(pol) + levelplot(r) + spplot(pol)
或者另一种选择:
library(latticeExtra)
spplot(pol) + spplot(r) + spplot(pol)
编辑
根据@OscarPerpián的评论,一个更好的方法是:
spplot(pol) + as.layer(spplot(r), under = TRUE)
https://stackoverflow.com/questions/27062768
复制相似问题