如何在水平图栅格中标记位置？

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控制大小。