ฉันกำลังมองหาไฟล์รูปร่างถนนสำหรับออสเตรเลียในวันนี้ ฉันสิ้นสุดขึ้นไปGeofabrik Open Street Map มีข้อมูลใน.osmรูปแบบ แต่พวกเขาไม่ได้มีหนึ่งใน.shpสำหรับออสเตรเลีย

ถ้าเป็นไปได้ฉันต้องการแปลงosmไฟล์เป็นshapefilesใช้ R. ฉันได้ค้นหาไปแล้ว แต่ยังไม่พบวิธีแก้ไข มีวิธีการนี้หรือไม่? อีกวิธีหนึ่งเป็นไปได้ไหมที่จะอ่านไฟล์ osm เป็น R และแปลงข้อมูลเป็น data frame เพื่อให้ฉันสามารถวาดรูปโดยใช้ggplot2?

ฉันได้ประสบความสำเร็จในภาพต่อไปนี้ของนิวซีแลนด์ใช้สอง shapefiles ด้วยและmaptools ggplot2เป็นการดีที่ฉันต้องการสร้างแผนที่ที่คล้ายกันของออสเตรเลีย

นี่ไม่ใช่โซลูชัน R แต่Quantum GIS (QGIS)เป็นวิธีที่ยอดเยี่ยมในการบรรลุสิ่งที่คุณต้องการ

คุณสามารถโหลดไฟล์. osm (เพิ่มเครื่องมือ Vector) คลิกขวาในสารบัญและบันทึกเป็นไฟล์ ESRI Shapefile

QGIS อาจมีปัญหากับสารสกัดขนาดใหญ่ดังนั้นเพื่อหลีกเลี่ยงปัญหานี้คุณสามารถใช้เครื่องมือ OSMเช่นOverPass APIเพื่อดาวน์โหลดเฉพาะสิ่งที่คุณต้องการโดยใช้กล่องขอบ

สะพานลอย-Turbo Api นอกจากนี้ยังมีที่จะได้รับสารสกัดจากกวดวิชาสั้น ๆ เกี่ยวกับที่อยู่ที่นี่!

ฉันวิ่งตัวอย่างรวดเร็วโดยใช้ทางหลวง = primary และ highway = primary_link tags (โครงการการติดแท็กทางหลวง OSM สามารถดูได้ที่นี่! ) โดยใช้ตัวช่วยสร้างบน Overpass-Turbo และภาพด้านล่างเป็นผลลัพธ์สำหรับ Victoria

ฉันส่งออกข้อมูลเป็น GeoJSON แล้วโหลดลงใน QGIS จากนั้นบันทึกผลลัพธ์เป็นไฟล์รูปร่าง (ภาพที่สองแสดงเส้นและ polys ที่โหลดเข้าสู่ QGIS)

ทางเลือกอื่น ๆ คือการดาวน์โหลดไฟล์ PBF หรือ OSM พื้นที่จาก GeoFabrik และย่อยข้อมูลโดยแยกทางหลวง = * แท็กโดยใช้Osmosis หากคุณต้องการอัปเดตข้อมูลเล็กน้อยเป็นประจำออสโมซิสจะเป็นวิธีที่แนะนำให้ดำเนินการต่อไป ถ้ามันเป็นหนึ่งแยกแยก, Overpass อาจจะง่ายขึ้นแม้ว่าคุณจะต้องทำมันในกล่องขนาดเล็กที่มีขอบเขตเนื่องจากหน่วยความจำ จำกัด คุณเพียงแค่ใช้แบบสอบถามแบบยกระดับเดียวกันทำกล่องขอบเขตที่แตกต่างกัน

@ 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ติดตั้ง (version> = 1.0.4)

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

• ตรวจสอบให้แน่ใจว่าgdalติดตั้ง (version> = 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 ข้อมูลถูกแบ่งตามภูมิภาค: อเมริกาเหนือและใต้, ออสเตรเลียและโอเชียเนีย, แอฟริกา, ยูโรปาและเอเชียชั้น 69 สำหรับการดาวน์โหลด ข้อมูลถูกกรองตามเงื่อนไข "คุณสมบัติของแผนที่"