Сегодня я искал дорожный шейп-файл для Австралии. Я закончил тем, что пошел в Geofabrik . Открытая карта улиц содержит данные в .osmформате. Но они не имеют в .shpАвстралии.

Если возможно, я хочу преобразовать osmфайлы в shapefilesформат R. Я искал вокруг, но пока не нашел решения. Есть ли способ к этому? Как вариант, можно ли прочитать файлы osm в R и преобразовать данные во фрейм данных, чтобы я мог рисовать фигуры, используя ggplot2?

Я получил следующий рисунок NZ, используя два шейп-файла с maptoolsи ggplot2. В идеале я хочу выпустить аналогичную карту Австралии.

Это не решение R, но Quantum GIS (QGIS) - отличный способ достичь того, чего вы хотите.

Вы можете просто загрузить файл .osm (инструмент Добавить вектор), щелкнуть его правой кнопкой мыши в содержании и сохранить как ESRI Shapefile.

QGIS может аварийно завершить работу с таким большим извлечением, поэтому во избежание этого вы можете использовать OSM Tools, такие как OverPass API, для загрузки только того, что вам нужно, с помощью ограничительных рамок.

Путепровод-Turbo Api также доступна для получения экстрактов, краткое руководство по это здесь!

Я запустил быстрый пример, основанный на тегах highway = primary и highway = primary_link (схему тегирования OSM Highway можно посмотреть здесь! ) С использованием мастера на Overpass-Turbo, и изображение ниже было результатом для Victoria.

Затем я экспортировал данные как GeoJSON, загрузил их в QGIS и затем сохранил результат в виде файла формы. (Второе изображение показывает линии и полисы, загруженные в QGIS)

Другой альтернативой является загрузка файла PBF или OSM для области из GeoFabrik и подмножество данных путем извлечения тегов шоссе = * с использованием Osmosis . Если вы хотите регулярно обновлять данные, рекомендуется использовать Osmosis. Если это одноразовая выдержка, Overpass, вероятно, будет проще, даже если вам придется делать это в меньших ограничивающих прямоугольниках из-за ограничений памяти. Вы бы просто применили одни и те же запросы Overpass для разных ограничивающих рамок.

@jazzurro, вы прекрасно можете сделать это с помощью R, просто посмотрите пакет osmar! Прочитайте документацию osmar (osmar.r-forge.r-project.org/RJpreprint.pdf). На страницах 11 с. Вы можете найти подробный пример извлечения дорог / автомагистралей по соответствующим тегам для munich.osm! После извлечения и извлечения данных из файла планеты для Австралии вы можете преобразовать их в любой формат!

Редактировать:

Поскольку некоторые комментаторы жаловались на отсутствие примеров, я выложу пример из документации. ИМХО, нет необходимости перепечатывать существующие примеры здесь, не так ли?

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")

Хорошо, вот правильный ответ:

• Убедитесь, что rgdal(версия> = 1.0.4) установлена

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

• Убедитесь, что gdal(версия> = 1.11.0) установлена

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

• Убедитесь, что gdalон скомпилирован с поддержкой Expat / OSM и SQLite:

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

• Убедитесь, что вы знаете, какой слой вы хотите сохранить как шейп-файл:

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

• Мы готовы пойти:

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

После того, как вы прочитаете файл readOGR, следуйте этим рекомендациям, чтобы узнать, как построить его ggplot2.

Обратите внимание, что вы также можете читать .osmфайлы в формате XML, просто убедитесь, что они не сжаты (то есть расширение .osmотсутствует .osm.bz2), но попробуйте использовать файл .osm.pbf, так как они намного меньше.

osm2shp.ru здесь вы можете скачать данные openstreetmap в формате shapefiles. Данные разделены по регионам: Северная и Южная Америка, Австралия и Океания, Африка, Европа и Азия. 61 слой для загрузки. Данные отфильтрованы по условиям «Особенности карты».