如何在R中的GetMap (RgoogleMaps)上叠加一个经度网格？

Question

我使用(RgoogleMaps)包检索了一张地图，使用以下Rcode在上绘制点：

library(RgoogleMaps)
lat = c(-30.3022,-30.5000,-33.48569)
lon = c(153.1189,151.6500,145.5316)
center = c(mean(lat), mean(lon))
zoom <- min(MaxZoom(range(lat), range(lon)))
mymap <- GetMap(center=center, zoom=zoom, maptype= "terrain", destfile = "MyTile1.png")

我还成功地在地图上绘制了我的3个点，使用：

NewMap <- PlotOnStaticMap(mymap, lat = c(-30.3022,-30.5000,-32.24300),  
      lon = c(153.1189,151.6500,148.6019), destfile = "MyTile1.png", cex=1.5,pch=20, 
      col=c('red', 'purple', 'green'), add=FALSE)       

现在我需要覆盖一个Lat-long值的网格，甚至是任何类型的网格。有什么想法吗？我做了一些相当广泛的研究，似乎(RgoogleMaps)没有一个简单的方法来做到这一点。

谢谢，D.A.S。

Answer 1

为此，您可以使用dismo::gmap。它以RasterLayer的形式返回谷歌地图，您可以使用它来叠加空间*对象(sp包)，例如SpatialPolygons或SpatialLines，以及栅格*对象(栅格对象)，或者只使用points或lines

library(dismo)
lat = c(-30.3022,-30.5000,-33.48569)
lon = c(153.1189,151.6500,145.5316)
xy <- cbind(lon, lat)

g <- gmap(xy, lonlat=TRUE, scale=2)
plot(g, interpolate=TRUE)
points(xy, col='red', pch=20)