ggplot map: "for" -Schleife fügt nur die letzte Kartenebene hinzu

Question

Eine ähnliche Frage wurde schon früher gestellt, aber zum Hinzufügen von Daten aus verschiedenen Spalten eines Datenrahmens zum Plot. Link to the question

In meinem Fall möchte ich neue Kartenebenen, wenn sie in der globalen Umgebung existieren, meiner endgültigen Karte mit einer "for" -Schleife hinzufügen.

Mit dem Code, wie es jetzt ist, wird nur die letzte Schicht, die Punkte, auf der Karte eingezeichnet, aber nicht die Leitungsschicht:

rm(list = ls()) 
library(OpenStreetMap) 
library(ggplot2) 
library(sp) 

### The base layer map 
map <- OpenStreetMap::openmap(c(55.8759, -4.2946), c(55.8638, -4.2776), type="osm") 
map <- openproj(map) 
myMap <- ggplot2::autoplot(map) 

### Create some custom map files 
# Some spatial point 
new_shapefile <- data.frame(long = -4.290, lat = c(55.868, 55.872)) 
sp::coordinates(new_shapefile) <- c("long", "lat") 
sp::proj4string(new_shapefile) <- sp::CRS("+init=epsg:4326") 

# Some spatial line 
x <- c(-4.290,-4.285) 
y <- c(55.868,55.868) 
new_line <- SpatialLines(list(Lines(Line(cbind(x,y)), ID="a"))) 
new_line = SpatialLinesDataFrame(new_line, data.frame(Z = c("Road"), row.names = c("a"))) 


addLayer <- function(){ 
    glob <- globalenv() 
    customFiles <- data.frame(ls(pattern = "^(?i)new", envir = glob)) 
    colnames(customFiles) <- "X" 
    customFiles$X <- as.character(customFiles$X) 
    for(i in 1:length(customFiles$X)){ 
    plot <- myMap 
    gg.data <- get(paste(customFiles$X[i])) 
    if(grepl("^SpatialPolygons", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_polygon(data = gg.data, aes(long, lat, group= group)) 
    } 
    if(grepl("^SpatialLines", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_path(data = gg.data, aes(long, lat, group= group)) 
    } 
    if(grepl("^SpatialPoints", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_point(data = data.frame(coordinates(gg.data)), aes(long, lat)) 
    } 
    } 
    print(plot) 
} 


### Run the function 
addLayer() 

Da ich nicht verwenden können, melt, um dies zu lösen, wie in der vorherigen Frage vorgeschlagen: Gibt es eine andere Möglichkeit, alle Schichten in der "for" -Schleife zu plotten?

Answer 1

Sie initiieren derzeit in jeder Iteration der Schleife Ihre Basiszuordnung, die plot ersetzt, die Sie in der vorherigen Iteration erstellt haben. Initiieren Sie die Basiszuordnung, bevor Sie die Schleife starten, um alle Ebenen derselben Karte hinzuzufügen.

addLayer <- function(){ 
    glob <- globalenv() 
    customFiles <- data.frame(ls(pattern = "^(?i)new", envir = glob)) 
    colnames(customFiles) <- "X" 
    customFiles$X <- as.character(customFiles$X) 
    plot <- myMap 
    for(i in 1:length(customFiles$X)){ 
    gg.data <- get(paste(customFiles$X[i])) 
    if(grepl("^SpatialPolygons", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_polygon(data = gg.data, aes(long, lat, group= group)) 
    } 
    if(grepl("^SpatialLines", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_path(data = gg.data, aes(long, lat, group= group)) 
    } 
    if(grepl("^SpatialPoints", class(get(paste(customFiles$X[i])))) == TRUE){ 
     plot <- plot + geom_point(data = data.frame(coordinates(gg.data)), aes(long, lat)) 
    } 
    } 
    print(plot) 
}