如何在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控制尺寸。