overpass is a packge with tools to work with the OpenStreetMap (OSM) Overpass API. To explore simple Overpass queries interactively, try overpass turbo.
Here's an RPub for overpass that I'll continually update as this goes (that will eventually be a vignette).
The Overpass API (or OSM3S) is a read-only API that serves up custom selected parts of the OSM map data. It acts as a database over the web: the client sends a query to the API and gets back the data set that corresponds to the query.
Unlike the main API, which is optimized for editing, Overpass API is optimized for data consumers that need a few elements within a glimpse or up to roughly 100 million elements in some minutes, both selected by search criteria like e.g. location, type of objects, tag properties, proximity, or combinations of them. It acts as a database backend for various services.
Overpass API has a powerful query language (language guide, language reference, an IDE) beyond XAPI, but also has a compatibility layer to allow a smooth transition from XAPI.
This package pairs nicely with nominatim.
The following functions are implemented:
overpass_query: Issue OSM Overpass Queryread_osm: Read an XML OSM Overpass response from pathopq: Begin building an Overpass queryadd_feature: Add a feature to an Overpass queryissue_query: Finalize and issue an Overpass queryavailable_features: List recognized features in OSM Overpassavailable_tags: List tags associated with a feature
- Version 0.1.1.9000 released - Bounding box tweaks
- Version 0.1.0.9000 released - trying out
opq/add_feature/issue_query - Version 0.0.0.9005 released - HUGE speed improvement; passes CRAN checks;
- Version 0.0.0.9000 released
devtools::install_github("hrbrmstr/overpass")library(overpass)
library(sp)
library(ggmap)
# current verison
packageVersion("overpass")
#> [1] '0.1.1.9000'# CSV example
osmcsv <- '[out:csv(::id,::type,"name")];
area[name="Bonn"]->.a;
( node(area.a)[railway=station];
way(area.a)[railway=station];
rel(area.a)[railway=station]; );
out;'
opq <- overpass_query(osmcsv)
read.table(text = opq, sep="\t", header=TRUE,
check.names=FALSE, stringsAsFactors=FALSE)
#> @id @type name
#> 1 26945519 node Bonn-Oberkassel
#> 2 1271017705 node Bonn-Beuel
#> 3 2428355974 node Bonn-Bad Godesberg
#> 4 2713060210 node Bonn Hauptbahnhof
#> 5 3400717493 node Bonn-Mehlem# just nodes
only_nodes <- '[out:xml];
node
["highway"="bus_stop"]
["shelter"]
["shelter"!~"no"]
(50.7,7.1,50.8,7.25);
out body;'
pts <- overpass_query(only_nodes)
plot(pts)
# ways & nodes
nodes_and_ways <- '[out:xml];
(node["amenity"="fire_station"]
(50.6,7.0,50.8,7.3);
way["amenity"="fire_station"]
(50.6,7.0,50.8,7.3);
rel["amenity"="fire_station"]
(50.6,7.0,50.8,7.3););
(._;>;);
out;'
wys <- overpass_query(nodes_and_ways)
plot(wys)
# xml version of the query
actual_ways <- '<osm-script output="xml">
<query type="way">
<bbox-query e="7.157" n="50.748" s="50.746" w="7.154"/>
</query>
<union>
<item/>
<recurse type="down"/>
</union>
<print/>
</osm-script>'
awy <- overpass_query(actual_ways)
plot(awy)
# more complex example from Robin: motorways surrounding London
# warning: may take a few minutes to run
from_robin <- '[out:xml][timeout:100];
(
node["highway"="motorway"](51.24,-0.61,51.73,0.41);
way["highway"="motorway"](51.24,-0.61,51.73,0.41);
relation["highway"="motorway"](51.24,-0.61,51.73,0.41);
);
out body;
>;
out skel qt;'
frb <- overpass_query(from_robin)
gg <- ggplot()
gg <- gg + geom_path(data=fortify(frb),
aes(x=long, y=lat, group=group),
color="black", size=0.25)
gg <- gg + coord_quickmap()
gg <- gg + ggthemes::theme_map()
gg
ggsave("README-from_robin-1.png")
library(overpass)
library(testthat)
date()
#> [1] "Wed Aug 12 17:06:32 2015"
test_dir("tests/")
#> testthat results ===========================================================
#> OK: 4 SKIPPED: 0 FAILED: 0
#>
#> DONEPlease note that this project is released with a Contributor Code of Conduct. By participating in this project you agree to abide by its terms.
React
Typescript
Django
Laravel
Facebook
Microsoft
Google
Alibaba
D3
Tencent