如何在R中使用levelplot添加框層？

Question

我想在情節中添加一些特殊范圍的框。

例

     gh <- raster()
     gh[] <- 1:ncell(gh)
     SP <- spsample(Spatial(bbox=bbox(gh)), 10, type="random")

然后畫他們

  levelplot(gh, col.regions = rev(terrain.colors(255)), cuts=254, margin=FALSE) +
  layer(sp.points(SP, col = "red"))

這將繪制一個帶有多個十字的地圖，但是我需要繪制一個具有一定范圍的框：

     extent(gh) = extent(c(xmn=-180,xmx=180,ymn=-90,ymx=90))
      e6 <- extent( 2  , 8 , 45   , 51  )

我想在情節上加上e6並將數字2放在方框內。

Answer 1

將Extent對象轉換為SpatialPolygons ，並使用以下coordinates提取其質心：

library("raster")
library("sp")
library("rasterVis")

gh <- raster()
gh[] <- 1:ncell(gh)
SP <- spsample(Spatial(bbox=bbox(gh)), 10, type="random")

e6 <- extent( 2, 8, 45, 51)
e6pol <- as(e6, 'SpatialPolygons')
centroid <- coordinates(e6pol)

levelplot(gh, col.regions = rev(terrain.colors(255)), cuts=254, margin=FALSE) +
    layer({sp.points(SP, col = "red")
           sp.polygons(e6pol)
           panel.text(centroid[,1], centroid[,2], '2')
           })