Hoy estaba buscando un camino para Australia. Terminé yendo a Geofabrik . Open Street Map tiene datos en .osmformato. Pero no tienen uno .shppara Australia.

Si es posible, quiero convertir los osmarchivos para shapefilesusar R. He buscado, pero aún no he encontrado soluciones. ¿Hay alguna forma de hacerlo? Alternativamente, ¿es posible leer archivos osm en R y convertir los datos en un marco de datos para poder dibujar figuras usando ggplot2?

He logrado el siguiente gráfico de NZ usando dos archivos de forma con maptoolsy ggplot2. Idealmente, quiero producir un mapa similar de Australia.

Esta no es una solución R, pero Quantum GIS (QGIS) es una excelente manera de lograr lo que desea.

Simplemente puede cargar el archivo .osm (herramienta Agregar vector), hacer clic con el botón derecho en la Tabla de contenido y Guardar como archivo de forma ESRI.

QGIS puede bloquearse con un extracto tan grande, por lo que para evitar esto, puede usar las herramientas OSM como la API OverPass para descargar solo lo que necesita usando cuadros delimitadores.

El OverPass-Turbo Api también está disponible para obtener extractos, ¡un breve tutorial sobre eso está aquí!

Ejecuté un ejemplo rápido basado en las etiquetas highway = primary y highway = primary_link (¡El esquema de etiquetado de autopistas OSM se puede ver aquí! ) Usando el asistente en Overpass-Turbo y la imagen a continuación fue el resultado para Victoria.

Luego exporté los datos como GeoJSON, los cargué en QGIS y luego guardé el resultado como un archivo de forma. (La segunda imagen muestra las líneas y polys cargados en QGIS)

La otra alternativa es descargar el archivo PBF u OSM para el área de GeoFabrik y subconjunto de los datos extrayendo las etiquetas highway = * usando Osmosis . Si desea actualizar sus datos regularmente, entonces Osmosis sería la forma recomendada de proceder. Si se trata de un extracto único, el Overpass probablemente sería más fácil, aunque tenga que hacerlo en cuadros delimitadores más pequeños debido a las limitaciones de memoria. Simplemente aplicaría las mismas consultas de Overpass en diferentes cuadros delimitadores.

@jazzurro, puedes hacer esto perfectamente con R, ¡solo busca el paquete osmar! Lea la documentación de osmar (osmar.r-forge.r-project.org/RJpreprint.pdf). ¡En las páginas 11 págs. Puede encontrar un ejemplo detallado para extraer carreteras / autopistas mediante las etiquetas correspondientes para munich.osm! ¡Después de extraer y extraer los datos de un archivo de planeta para Australia, puede convertir a cualquier formato que desee!

Editar:

Como algunos comentaristas se quejaban de la falta de ejemplos, publicaré un ejemplo de los documentos. En mi humilde opinión, no sería necesario volver a escribir los ejemplos existentes aquí, ¿verdad?

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, aquí viene la respuesta correcta:

• Asegúrese de que rgdal(versión> = 1.0.4) esté instalado

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

• Asegúrese de que gdal(versión> = 1.11.0) esté instalado

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

• Asegúrese de que gdalesté compilado con soporte para Expat / OSM y SQLite:

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

• Asegúrese de saber qué capa le gustaría guardar como shapefile:

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

• Estamos listos para irnos:

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

Una vez que haya leído el archivo readOGR, siga estas pautas para aprender a trazarlo ggplot2.

Tenga en cuenta que también puede leer .osmarchivos en formato XML, solo asegúrese de que no estén comprimidos (es decir, la extensión .osmno lo está .osm.bz2). Pero intente usar el archivo .osm.pbf ya que son mucho más pequeños.

osm2shp.ru aquí puede descargar datos de openstreetmap en formato shapefiles. Datos divididos por regiones: Norteamérica y Sudamérica, Australia y Oceanía, África, Europa y Asia.61 capas para descargar. Datos filtrados por las condiciones de "Características del mapa".