如何绘制SpatialLinesDataFrame特征映射在谷歌地图

Question

我想要绘制在谷歌地图的线路层。

数据

> dput(map) 
new("SpatialLinesDataFrame" 
     , data = structure(list(att = c(463643, 2291491, 315237340, 10348934, 
309845150, 674351, 58057, 55962, 302861, 1405635)), .Names = "att",  row.names = c(NA, 
10L), class = "data.frame") 
    , lines = list(<S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>, 
    <S4 object of class structure("Lines", package = "sp")>) 
    , bbox = structure(c(50.497608475813, 26.1186426230732, 50.6164182652142, 
26.2649832975207), .Dim = c(2L, 2L), .Dimnames = list(c("x", 
"y"), c("min", "max"))) 
    , proj4string = new("CRS" 
    , projargs = "+proj=longlat +ellps=WGS84 +towgs84=0,0,0,-0,-0,-0,0 +no_defs" 
) 
) 

方法

library(rgdal) 
library(ggmap) 

gmap <- get_map(location=rowMeans(bbox(segMap)), zoom = 11) # get Google map to use as background  

变我

plot(map, col = map$att, lwd = 1.5) 
plot(gMap) 

变II

plot(map, col = map$att, lwd = 1.5) 
ggmap(gMap) 

问题

背景地图绘制在特征图的顶部，而不是象背景，作为结果的特征图是不可见的。 为了澄清，无论是剧情要求（plot()和ggmap）做工精细独立。 谢谢

Answer 1

像这样的事情？

所有的

library(raster) # for getData(...), also loads sp 
library(ggmap)  # for get_map, also loads ggplot2 

map <- getData("GADM",country="GBR",level=2) # SpatialPolygonsDataFrame of UK 
map <- map[map$NAME_2=="London",]    # extract London 
# this builds a spatialLinesDataFrame object - random walks around London 
# you have this already... 
set.seed(1) # for reproducible example 
get.coords <- function() { 
    do.call(cbind,lapply(rowMeans(bbox(map)), 
           function(x)cumsum(sample(0.01*(-1:1),50,replace=TRUE))+x)) 
} 
route <- SpatialLines(lapply(1:3,function(i)Lines(list(Line(get.coords())),ID=i))) 
route <- SpatialLinesDataFrame(route, data.frame(att=c("A","B","C"))) 

# you would start here... 
gg.df <- do.call(rbind,lapply(route$att, function(x)data.frame(att=x,coordinates(route[route$att==x,])))) 

ggmap(get_map(location=rowMeans(bbox(map)),zoom=12)) + 
    geom_path(data=gg.df, aes(x,y, color=att), size=2)+ 
    geom_point(data=gg.df, aes(x=x[1], y=y[1]),color="black",size=5) 

首先，感谢试图包含您的数据。通常情况下，dput(...)是要走的路，但空间对象只是没用。您确实需要上传文件（或文件，通常）在某处并发布链接。

第二，这不是一个简单的问题。这ggmap(...)产生可以通过增加使用+附加功能来增强一个ggplot对象，而且它也确实ggplot提供了一个非常强大的功能fortify(...)转换SpatialPolygonsDataFrames成ggplot可用的东西这是真的。但fortify(...)不适用于对象，所以我们必须明确地执行此操作。

因此，这段代码的大部分只是创建一个带有三行的SpatialLinesDataFrame，每行在属性表中都有对应的属性att。这应该很好地模仿你的例子。由于这是一个空间对象，所以在sp,rgeos等上的所有方法都可以对其进行操作。

现在，ggplot需要一个data.frame，其中的列对应于在调用aes(...)时映射的美学，因此这里x和y代表点的位置，颜色代表线的颜色。我们创建一个使用：

gg.df <- do.call(rbind,lapply(route$att, function(x)data.frame(att=x,coordinates(route[route$att==x,])))) 

，提取用于att每个值到data.frames列表的坐标，并结合这些一起行方式。

已经这样做了，这是一个简单的事情，添加行使用

ggmap(get_map(location=rowMeans(bbox(map)),zoom=12)) + 
    geom_path(data=gg.df, aes(x,y, color=att), size=2)+ 
    geom_point(data=gg.df, aes(x=x[1], y=y[1]),color="black",size=5) 

需要注意的是，（一）你需要使用geom_path(...)，因为geom_line(...)将整理x轴，和（b）你需要告诉ggplot有一个特定于层的数据集（gg.df）。