Я шукав дорожню форму для Австралії сьогодні. Я врешті-решт поїхав до Геофабрику . Відкрита карта вулиць має дані у .osmформаті. Але їх немає в .shpАвстралії.

Якщо можливо, я хочу перетворити osmфайли на shapefilesR. Я шукав, але ще не знайшов рішення. Чи є до цього спосіб? Як варіант, чи можна зчитувати файли osm в R і перетворювати дані в кадр даних, щоб я міг малювати фігури, використовуючи ggplot2?

Я домігся наступної графіки NZ, використовуючи два файли з maptoolsі ggplot2. В ідеалі я хочу створити подібну карту Австралії.

Це не рішення R, але квантовий ГІС (QGIS) - прекрасний спосіб досягти того, що ви хочете.

Ви можете просто завантажити.

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

Шляхопровід-Turbo Api також доступна для отримання екстрактів, інструкція по це тут!

Я наводив короткий приклад на основі тегів шосе = первинне та шосе = первинне посилання (схему маркування шосе OSM можна переглянути тут! ) За допомогою Майстра на Overpass-Turbo, і зображення нижче було результатом для Вікторії.

Потім я експортував дані як GeoJSON, завантажив їх у QGIS, а потім зберег результат як файл форми. (На другому зображенні показані лінії та поліси, завантажені в QGIS)

Інша альтернатива - завантажити файл PBF або OSM для області з GeoFabrik та підмножити дані, витягнувши теги шосе = * за допомогою осмосу . Якщо ви хочете регулярно оновлювати свої дані, рекомендованим способом продовжуватиметься осмос. Якщо це одноразовий витяг, Overpass, мабуть, буде простішим, хоча вам доведеться робити це в менших обмежувальних полях через обмеження пам'яті. Ви просто застосуєте одні й ті ж запити Overpass для різних обмежувальних полів.

@jazzurro, ви прекрасно можете це зробити з R, просто шукайте 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

• Переконайтеся, що ви знаєте, який шар ви хочете зберегти як shapefile:

  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. Дані розділені за регіонами: Північна та Південна Америка, Австралія та Океанія, Африка, Європа та Азія. Дані відфільтровані за умовами "Особливості карти".