osm-admin-lookup

The project aims to provide performant lookups of administrative hierarchies for coordinates using OSM data and fast lookup trees.

Prepare

Download OSM bundle.

curl https://download.geofabrik.de/europe/germany/berlin-latest.osm.pbf -o berlin-latest.osm.pbf

Extract administrative boundaries from OSM to speed up processing.

brew install osmium-tool
osmium tags-filter berlin-latest.osm.pbf \
  r/boundary=administrative \
  -o berlin-boundaries.pbf

Build

cargo build --release

Create RTree

./target/release/build-rtree --bin rtree.bin --pbf berlin-boundaries.pbf

Locate point

List boundaries.

./target/release/locate -b rtree.bin -l 13.4,52.5
boundary: Berlin, level: 4
boundary: Kreuzberg, level: 10
boundary: Friedrichshain-Kreuzberg, level: 9

Compile geojson file with boundaries.

./target/release/locate -b rtree.bin -l 13.4,52.5 -g boundaries.geojson
cat boundaries.geojson | pbcopy
# paste in geojson.io or similar

Benchmark

The benchmark requires a pre-built rtree (w/ build-rtree) and a CSV file with locations (columns: id, lng, lat).

./target/release/bench -- single \
  --bin brandenburg-rtree.bin \
  --locs 4000_locs.csv \
  -m 16
took 648.188979ms for 4000 requests
took 672.047414ms for 4000 requests
took 648.485565ms for 4000 requests

Web Service

The web service requires a pre-built rtree (w/ build-rtree) and accepts coordinates in longitude,latitude format.

cargo run --release --bin admin-lookup -- --bin rtree.bin

export LOC=13.425979614257812,52.53919655252312
curl -s "localhost:8080/locate?loc=$LOC | jq .
{
  "boundaries": [
    {
      "level": 4,
      "name": "Berlin"
    },
    {
      "level": 10,
      "name": "Prenzlauer Berg"
    },
    {
      "level": 9,
      "name": "Pankow"
    }
  ]
}