小编mas*_*don的帖子

我有一个包含科隆所有街道的shapefile(SpatialLinesDataFrame),可以从这里下载.我将此@data与来自外部源的数据合并.我如何绘制这些街道(如果可能在谷歌地图上使用ggmaps),以便每条街道都有不同的颜色(或厚度),具体取决于它的个别价值？

到目前为止,我已经这样做了:

shapefile <- readOGR(shapfile, "Strasse", stringsAsFactors=FALSE, 
encoding="latin-9")
shp <- spTransform(shapefile, CRS("+proj=longlat +datum=WGS84"))

此时我在shp @ data数据框中添加了另一列,其中包含每条街道的特定值.然后我强化了shapefile,因此可以使用ggplot绘制:

shp$id <- rownames(shp@data)
shp.df <- as.data.frame(shp)
data_fort <- fortify(shp, region = "id")
data_merged <- join(data_fort, shp.df, by="id")

当我使用geom_lines时,这些线看起来不太好并且不容易识别:

ggplot(data_merged, aes(x=long, y=lat,
                      group=group,
                      colour=values)) +
 geom_line()

在这里,我看到可以转换shapefile,以便可以使用geom_segement(或者在这种情况下是修改后的函数"geom_segment2"),但随后会丢失我的街道特定值.