1

次の Rcode を使用して (RgoogleMaps) パッケージを使用して、ポイントをプロットするマップを取得しました。

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)

ここで、緯度経度値のグリッド、または任意のタイプのグリッドをオーバーレイする必要があります。何か案は？かなり広範な調査を行ったところ、(RgoogleMaps) にはこれを行う簡単な方法がないようです。

ありがとう、DAS

4

1 に答える 1

1

これには dismo::gmap を使用できます。pointsこれは Google マップを RasterLayer として返し、これを使用して Spatial* オブジェクト (sp パッケージ)、たとえば SpatialPolygons または SpatialLines、および Raster* オブジェクト (ラスター オブジェクト) をオーバーレイするか、または単に使用することができます。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)
于 2015-11-23T09:16:56.640 に答える