#openstreetmap #osm #topology #osm-pbf

app osm-lump-ways

Group OpenStreetMap lines based on graph topology and OSM tags

31 releases (breaking)

new 0.25.0 Apr 25, 2024
0.23.0 Mar 19, 2024
0.18.0 Nov 30, 2023
0.3.0 Jul 23, 2023
0.1.0 Mar 31, 2023

#8 in Geospatial

Download history 5/week @ 2024-01-25 116/week @ 2024-02-15 110/week @ 2024-02-22 218/week @ 2024-02-29 167/week @ 2024-03-07 229/week @ 2024-03-14 33/week @ 2024-03-21 16/week @ 2024-03-28 104/week @ 2024-04-04 43/week @ 2024-04-11

233 downloads per month

MIT/Apache and AGPL-3.0

250KB
5K SLoC

osm-lump-ways group OSM ways based on topology & shared tags

Crates.io Number of Downloads Crates.io Latest Version

Answer questions about OSM data like:

  • “What's the longest Main Street?”
  • “How far can I drive on unpaved roads in this region?”
  • “How long is the M1 motorway?”
  • “Are these rivers connected?”
  • “What's the river drainage basins?”

2 similar programmes are included: osm-lump-ways, which ignores the direction of the OSM way, and osm-lump-ways-down, which uses direction of OSM ways to produce data, incl. QA suitable files. Both share similarities.

Background

OSM linear features (eg roads, rivers, walls) are stored as way object. The OSM tagging model often requires one feature to be mapped as many different ways. osm-lump-ways will assemble them all together.

Filtering OSM Data

There are 2 ways to select which OSM ways will be used. All relations are currently ignored.

Tag Filter Rules

  • key / ∃key way has this tag
  • ~key_regex / ~key_regex There is a key, which matches this regex.
  • ∄key way does not has this key
  • key=value way has this key and this value
  • key≠value way either doesn't have this key, or if it does, it's not equal to value
  • key∈value1,value2, way has this key and the value is one of these
  • key∉value1,value2, way either doesn't have this key, or if it does, it's not one of these values
  • key~regex way has this key and the value matches this regex.
  • F1F2F3 logical OR of the other tag filters F1, F2, …
  • F1F2F3 logical AND of the other tag filters F1, F2, …

The popular regex crate is used for matching. Regexes, and string comparison, are case sensitive. Add (?i) at start of regex to switch to case insensitive (e.g. name~(?i).* street) Regexes match the whole value, name~[Ss]treet will match Street or street, but not Main Street North nor Main Street. Use name~.*[Ss]treet.* to match all.

Simple Tag Filtering

The -f/--tag-filter can be specified one or more times, and an OSM object is only included if it matches all defined filter's, i.e. a logical AND of all filters.

  • -f highway: Only ways with a highway tag are included
  • -f highway -f name: Only ways with a highway and name tag are included.
  • -f highway -f ∄name: Only ways with a highway and without a name tag are included.

More complex Tag Filtering Functions

The -F/--tag-filter-func takes a single ordered list of tag filters (separated by ;), and includes (with →T), or excludes (with →F), the OSM object based on the first filter function which matches. A bare T or F applies to all.

Example: We want to include all waterways. But not waterway=canal. But we want a waterway=canal iff it also has a lock=yes tag.

-F "waterway=canal∧lock=yes→T; waterway=canal→F; waterway→T; F

If the argument to -F/--tag-filter-func starts with @, the rest is a filename containing the tag filter func code. e.g. -F @myrules.txt.

Comments start with # and continue to the end of the line.

Output format

If a filename ends with .geojson, a GeoJSON file (RFC 7946 will be created. For .geojsons, a GeoJSON Text Sequences (RFC 8142), aka GeoJSONSeq, file.

osm-lump-ways

Usage

Generate river drainage basins

osm-lump-ways -i path/to/region-latest.osm.pbf -o region-rivers.geojson -f waterway=river

To group based on the river's name:

osm-lump-ways -i path/to/region-latest.osm.pbf -o region-rivers.geojson -f waterway=river -g name

To find long streets and assemble them into connected (Multi)LineStrings:

osm-lump-ways -i path/to/region-latest.osm.pbf -o long-streets.geojson -f highway -g name

Installation

cargo install osm-lump-ways

Full Options

Run with --help to see all options.

Frames

Here, a “frame” of a grouping is a shortest path through 2 points in the grouped together ways. This can be useful for waterways to find where a big group is connected.

--output-frames FILENAME.geojsons will save these features to a file. GeoJSONSeq output format only.

Examples of usage

osm-lump-ways-down

It reads & groups an OSM PBF file, like osm-lump-ways, but it uses the direction of the OSM way, to produce

The output filename must contain %s, which will be replaced

Output files

loops

Cycles in the network. Each is a strongly connected component.

upstreams

Each way segment (a 2 point LineString) with the upstream

upstreams-points

Above, but just the first point. Each feature a Point.

ends

Points where waterways end.

Ends membership of tagged ways

With the --ends-membership TAGFILTER arg every end will have a boolean property called is_in:TAGFILTER which is only true if this end node is a member of a way with this tag filter. This argument can be applied many times. An addditional property is_in_count is an integer of how many of these properties are true.

e.g. --ends-membership natural=coastline will cause each end point to have a JSON property is_in:natural=coastline which is true iff this node is also a member of a way with the natural=coastline tag, false otherwise.

Loop removal

After the loops are detected, all the edges (way segments) in the loops are contracted together, producing a new graph which is loop-free.

TODOs

This software isn't finished, here's what I'd like to add. Feel free to send a patch.

  • All tags need to be specified in advance to join on. Perhaps add something to match all possible tags? (inspired by this q). (“Group by all tags the same” might do it)

External Mentions

TBC

Copyright & Licence

Copyright 2023, 2024, MIT/Apache2.0. Source code is on Github (osm-lump-ways).

The output data file(s) are a Derived Database of the OpenStreetMap database, and hence under the ODbL 1.0 licence, the same as the OpenStreetMap copyright, and contains data © OpenStreetMap contributors

Dependencies

~53MB
~825K SLoC