繁体   English   中英

顶部栅格背景上的空间对象的图例

[英]legend of spatial objects on top raster background

有很好的示例如何在栅格图层上覆盖矢量/栅格地图,例如在R中将矢量地图绘制在栅格背景之上,在R 的背景GM图层上覆盖栅格图层 但是合并的图形中没有图例。 在示例中, rproblv具有图例,但migmaplv + rproblv没有。

有没有一种简单的方法可以合并两个levelplotspplot对象,并在结果中保留一个或两个图例?

R中将栅格图层覆盖在R的背景GM图层上

library(utils); library(colorRamps); library(rgdal);
library(raster); library(dismo); library(rasterVis)

download.file(
  "https://dl.dropboxusercontent.com/u/3180464/rprob520.tif", 
  "rprob520.tif", method = "curl")

rprob <- raster("rprob520.tif")
rprobGM <- projectRaster(from = rprob, crs = CRS("+init=epsg:3857"))

migmap <- gmap(x = rprob, type = "hybrid", zoom = 5)

migmaplv <- levelplot(
  migmap, maxpixels = ncell(migmap), 
  col.regions = migmap@legend@colortable, 
  at = 0:255, panel = panel.levelplot.raster, interpolate = TRUE, 
  colorkey = FALSE, margin = FALSE)

rproblv <- levelplot(
  rprobGM, margin = FALSE, contour = TRUE, 
  par.settings = rasterTheme(region = matlab.like(n = 10)), 
  alpha.regions = 0.35, at = (0:10)/10, 
  main = "p(dep > 1kg/m2 per eruption event)")

print(migmaplv + rproblv)

在此处输入图片说明

as.layer函数与under = TRUE

rproblv + as.layer(migmaplv, under = TRUE)

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM