如何在 levelplot 栅格中标记位置？

Question

我有一个光栅文件，我想显示我在这张地图中使用的站点的位置。不知道你有没有这方面的帮助：

 library(rasterVis)
 r <- raster(nrow=10, ncol=10)
 r[] = 1
 r[51:100] = 3
 r[3:6, 1:5] = 5
 r <- ratify(r)
 rat <- levels(r)[[1]]
 rat$landcover <- c('Pine', 'Oak', 'Meadow')
 rat$class <- c('A1', 'B2', 'C3')
 levels(r) <- rat
 levelplot(r, col.regions=c('palegreen', 'midnightblue', 'indianred1'))

我想显示 lat = 50 和 long=100 的地方。只是在地图上的这个点上画一个圆圈或正方形。

Answer 1

以下代码将在 lat=50，lon=100 处添加一个点：

levelplot(r, col.regions=c('palegreen', 'midnightblue', 'indianred1')) +
  layer(panel.points(100,50, pch=21, cex=3, colour='black', fill='red'))

pch 设置标记类型（参见 ?pch 了解各种标记类型），cex 控制大小。

如何在 levelplot 栅格中标记位置？

How to mark a place in a levelplot raster?

r

raster