在rasterVis levelplot

Question

我想改变海洋颜色（在shapefile边界之外）。我可以剪辑光栅并更改背景颜色，但是我想在这里添加shapefile。

library(raster)
library(rasterVis)
library(maps)
library(maptools)
library(mapdata)

r <- raster(nrow=361, ncol=576, ymn=-90, ymx=90) 
values(r) <- 1:ncell(r)
data(wrld_simpl, package = "maptools")
levelplot(r)+ layer(sp.polygons(wrld_simpl, lwd=0.1, col='gray'))

Answer 1

首先，使用SpatialPolygons对象屏蔽Raster。细胞未被覆盖 由它设置为NA。

land <- mask(r, wrld_simpl)

现在，更改背景颜色（用于NA细胞）：

catTheme <- rasterTheme(panel.background = list(col='lightskyblue1'))

最后，显示结果：

levelplot(land, par.settings = catTheme) +
    layer(sp.polygons(wrld_simpl,
                      lwd=0.1, col='gray'))