我今天在寻找澳大利亚的道路形状文件。我最终去了Geofabrik。Open Street Map具有.osm格式的数据。但是他们没有一个.shp适合澳大利亚。

如果可能的话，我想将osm文件转换为shapefiles使用R。我已经搜索过，但是还没有找到解决方案。有什么办法吗？另外，是否可以将osm文件读入R并将数据转换为数据帧，以便我可以使用绘制图形ggplot2？

我通过使用带有maptools和的两个shapefile实现了NZ的以下图形ggplot2。理想情况下，我想制作一张类似的澳大利亚地图。

这不是R解决方案，但是Quantum GIS（QGIS）是实现所需目标的好方法。

您可以简单地加载.osm文件（添加矢量工具），在目录中右键单击它，然后另存为ESRI Shapefile。

QGIS提取的内容如此庞大可能会崩溃，因此，为避免这种情况，您可以使用OSM工具（如OverPass API）使用边界框仅下载所需的内容。

该立交桥涡轮增压 API是提供获得提取物，对一个简短的教程是在这里！

我使用立交桥上的向导运行了一个基于Highway = primary和Highway = primary_link标签的快速示例（可在此处查看OSM高速公路标签方案！），下图是Victoria的结果。

然后，我将数据导出为GeoJSON，将其加载到QGIS中，然后将结果保存为形状文件。（第二张图显示了加载到QGIS中的线和多边形）

另一种选择是从GeoFabrik下载该区域的PBF或OSM文件，并通过使用Osmosis提取Highway = *标签来对数据进行子集化。如果您希望定期更新您的数据，则建议使用渗透法进行操作。如果它是一次性提取，则由于内存限制，即使必须在较小的边界框中执行此操作，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

• 确保知道要保存为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格式的文件，只需确保未压缩文件（即，扩展名为.osmnot .osm.bz2），但是请尝试使用.osm.pbf文件，因为它们要小得多。

osm2shp.ru在这里，您可以下载shapefile格式的openstreetmap数据。数据按地区划分：北美和南美，澳大利亚和大洋洲，非洲，欧洲和亚洲61层可供下载。根据“地图要素”条件过滤的数据。