오늘 호주의 도로 모양 파일을 찾고있었습니다. 나는 Geofabrik에 가는 것을 끝내었다 . Open Street Map에는 데이터 .osm형식이 있습니다. 그러나 그들은 .shp호주에 하나도 없습니다 .

가능한 경우 R osm을 shapefiles사용하여 파일 을 변환하려고합니다. 검색했지만 아직 해결책을 찾지 못했습니다. 이 방법이 있습니까? 또는 osm 파일을 R로 읽고 데이터를 데이터 프레임으로 변환하여 ggplot2?를 사용하여 그림을 그릴 수 있습니까?

나는 두 개의 shape 파일을 사용하여 뉴질랜드의 다음 그래픽 달성 maptools과 ggplot2. 이상적으로는 비슷한 호주지도를 만들고 싶습니다.

이 솔루션은 R 솔루션이 아니지만 Quantum GIS (QGIS) 는 원하는 것을 달성 할 수있는 좋은 방법입니다.

.osm 파일 (벡터 추가 도구)을로드하고 목차에서 마우스 오른쪽 단추로 클릭 한 후 ESRI 셰이프 파일로 저장을 클릭하면됩니다.

QGIS는 이러한 큰 추출로 충돌 할 수 있으므로이를 피하기 위해 OverPass API 와 같은 OSM 도구를 사용하여 경계 상자를 사용하여 필요한 것만 다운로드 할 수 있습니다.

육교 터보 API는 또한 추출물을 얻을 수 있습니다, 그에 관한 짧은 튜토리얼은 여기에!

Overpass-Turbo의 마법사를 사용하여 highway = primary 및 highway = primary_link 태그 (OSM Highway Tagging Scheme은 여기를 참조하십시오 ! )를 기반으로 한 빠른 예를 실행했으며 아래 이미지는 Victoria의 결과입니다.

그런 다음 데이터를 GeoJSON으로 내보내고 QGIS에로드 한 다음 결과를 모양 파일로 저장했습니다. (두 번째 이미지는 QGIS에로드 된 선과 폴리를 보여줍니다)

다른 대안은 GeoFabrik에서 해당 영역에 대한 PBF 또는 OSM 파일을 다운로드하고 Osmosis를 사용하여 highway = * 태그를 추출하여 데이터를 서브 세트하는 것 입니다. yor 데이터를 정기적으로 업데이트하려면 Osmosis가 권장되는 방법입니다. 추출이 한 번이면 메모리 제한으로 인해 작은 경계 상자에서해야하더라도 오버 패스가 더 쉬울 것입니다. 동일한 오버 패스 쿼리를 적용하여 다른 경계 상자를 적용하면됩니다.

@ jazzurro, 당신은 완벽하게 R로 이것을 할 수 있습니다, 그냥 osmar 패키지를 찾으십시오! osmar 문서 (osmar.r-forge.r-project.org/RJpreprint.pdf)를 읽으십시오. 11 페이지의 munich.osm에 대한 태그로 도로 / 고속도로를 추출하는 자세한 예를 찾을 수 있습니다! 호주의 행성 파일에서 데이터를 가져와 추출한 후 원하는 형식으로 변환 할 수 있습니다!

편집하다:

일부 주석가가 부족한 사례에 대해 불평하면서 문서에서 예제를 게시합니다. IMHO 여기에 기존 예제를 다시 입력 할 필요가 없습니까?

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"

• gdalExpat / 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.

.osmXML 형식의 파일을 읽을 수도 있으며 압축되지 않았는지 확인하십시오 (예 : 확장자가 .osm아닌지 확인 .osm.bz2). .osm.pbf 파일은 훨씬 작으므로 사용하십시오.

osm2shp.ru 여기서 shapestreet 형식의 오픈 스트리트 맵 데이터를 다운로드 할 수 있습니다. 북미 및 남미, 호주 및 오세아니아, 아프리카, 유로파 및 아시아 지역으로 나눈 데이터 : 61 개의 레이어를 다운로드 할 수 있습니다. "지도 기능"조건으로 필터링 된 데이터.