我想改变海洋的颜色(形状文件边界之外)。我可以剪辑光栅和改变背景颜色,但在这里,我想用添加的形状文件。
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'))

发布于 2016-06-22 05:58:51
首先,用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'))

https://stackoverflow.com/questions/37951456
复制相似问题