Ich habe heute nach einem Road Shapefile für Australien gesucht. Am Ende ging ich zur Geofabrik . Open Street Map enthält Daten im .osmFormat. Aber sie haben keinen .shpfür Australien.

Wenn möglich, möchte ich die osmDateien shapefilesmit R konvertieren . Ich habe gesucht, aber noch keine Lösungen gefunden. Gibt es einen Weg dorthin? Ist es alternativ möglich, Osm-Dateien in R zu lesen und die Daten in Datenrahmen umzuwandeln, damit ich mit Zahlen zeichnen kann ggplot2?

Ich habe die folgende Grafik von NZ mit zwei Shapefiles mit maptoolsund erstellt ggplot2. Idealerweise möchte ich eine ähnliche Karte von Australien produzieren.

Dies ist keine R-Lösung, aber Quantum GIS (QGIS) ist ein großartiger Weg, um das zu erreichen, was Sie wollen.

Sie können einfach die .osm-Datei laden (Vektor-Tool hinzufügen), im Inhaltsverzeichnis mit der rechten Maustaste darauf klicken und als ESRI-Shapefile speichern.

QGIS kann mit einem so großen Auszug abstürzen. Um dies zu vermeiden, können Sie OSM-Tools wie die OverPass-API verwenden , um nur das herunterzuladen, was Sie mithilfe von Begrenzungsrahmen benötigen.

Die OverPass-Turbo-API ist auch verfügbar, um Auszüge zu erhalten. Ein kurzes Tutorial dazu finden Sie hier!

Ich habe ein kurzes Beispiel basierend auf den Tags highway = primary und highway = primary_link (das OSM-Highway-Tagging-Schema finden Sie hier! ) Mit dem Assistenten auf Overpass-Turbo ausgeführt. Das folgende Bild war das Ergebnis für Victoria.

Ich exportierte die Daten dann als GeoJSON, lud sie in QGIS und speicherte das Ergebnis als Formdatei. (Das zweite Bild zeigt die Linien und Polys, die in QGIS geladen wurden.)

Die andere Alternative besteht darin, die PBF- oder OSM-Datei für das Gebiet von GeoFabrik herunterzuladen und die Daten durch Extrahieren der highway = * -Tags mit Osmosis zu unterteilen . Wenn Sie Ihre Daten regelmäßig aktualisieren möchten, ist Osmosis die empfohlene Vorgehensweise. Wenn es sich um einen einmaligen Auszug handelt, ist die Überführung wahrscheinlich einfacher, obwohl Sie dies aufgrund von Speicherbeschränkungen in kleineren Begrenzungsrahmen tun müssen. Sie würden einfach dieselben Überführungsabfragen anwenden, um unterschiedliche Begrenzungsrahmen zu erstellen.

@jazzurro, du kannst das perfekt mit R machen, schau einfach nach dem osmar-Paket! Lesen Sie die osmar-Dokumentation (osmar.r-forge.r-project.org/RJpreprint.pdf). Auf den Seiten 11 S. finden Sie ein ausführliches Beispiel zum Extrahieren von Straßen / Autobahnen mit den entsprechenden Tags für munich.osm! Nachdem Sie die Daten aus einer Planetendatei für Australien gezogen und extrahiert haben, können Sie sie in ein beliebiges Format konvertieren!

Bearbeiten:

Da sich einige Kommentatoren über fehlende Beispiele beschwerten, werde ich ein Beispiel aus den Dokumenten veröffentlichen. IMHO wäre es nicht nötig, existierende Beispiele hier noch einmal einzutippen, oder?

library(maptools)
library(osmar)
url <- "http://osmar.r-forge.r-project.org/"
file <- "muenchen.osm.gz"
download.file(sprintf("%s%s", url, file), file)
unzip("gzip -d muenchen.osm.gz") # gzip is linux only, on windows I unzipped this manually with 7zip!

src <- osmsource_osmosis(file = "muenchen.osm")
muc_bbox <- center_bbox(11.575278, 48.137222, 3000, 3000)
muc <- get_osm(muc_bbox, src)
muc
summary(muc)

hw_ids <- find(muc, way(tags(k == "highway")))
hw_ids <- find_down(muc, way(hw_ids))
hw <- subset(muc, ids = hw_ids)

plot(muc)
plot_ways(hw, add = TRUE, col = "green")

# convert to spatial object (SpatialLinesDataFrame)
# and save to whatever format you like..
hw_line <- as_sp(hw, "lines")

OK, hier kommt die richtige Antwort:

• Stellen Sie sicher, dass rgdal(Version> = 1.0.4) installiert ist

  install.packages('rgdal')
  packageVersion('rgdal')
  
  [1] ‘1.0.4’

• gdalStellen Sie sicher, dass (Version> = 1.11.0) installiert ist

  library(rgdal)
  getGDALVersionInfo()
  
  [1] "GDAL 1.11.2, released 2015/02/10"

• gdalStellen Sie sicher, dass mit Expat / OSM- und SQLite-Unterstützung kompiliert wurde:

  c('SQLite', 'OSM') %in% ogrDrivers()$name
  
  [1] TRUE TRUE

• Stellen Sie sicher, dass Sie wissen, welche Ebene Sie als Shapefile speichern möchten:

  ogrListLayers('filename.osm.pbf')
  
  [1] "points" "lines" "multilinestrings" "multipolygons"
  attr(,"driver")
  [1] "OSM"
  attr(,"nlayers")
  [1] 4

• Wir sind bereit zu gehen:

  osm <- readOGR('filename.osm.pbf', 'lines')
  writeOGR(osm, 'myshapedir', 'mylayer', driver = 'ESRI Shapefile')

Wenn Sie die Datei über gelesen haben readOGR, befolgen Sie diese Richtlinien , um zu erfahren, wie Sie sie plotten ggplot2.

Beachten Sie, dass Sie auch .osmDateien im XML-Format lesen können. Achten Sie jedoch darauf, dass diese nicht komprimiert sind (dh die Erweiterung ist .osmnicht vorhanden .osm.bz2). Versuchen Sie jedoch, die Datei .osm.pbf zu verwenden, da sie viel kleiner sind.

osm2shp.ru hier können sie openstreetmap daten im shapefiles format herunterladen. Daten unterteilt nach Regionen: Nord- und Südamerika, Australien und Ozeanien, Afrika, Europa und Asien.61 Ebenen zum Herunterladen. Daten gefiltert nach "Map Features" -Bedingungen.