Skip to content

Instantly share code, notes, and snippets.

@andrewharvey
Created January 7, 2020 07:05
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save andrewharvey/f30beadf25c91c59abedd95ee0c5aaf6 to your computer and use it in GitHub Desktop.
Save andrewharvey/f30beadf25c91c59abedd95ee0c5aaf6 to your computer and use it in GitHub Desktop.
Compare a road network to OpenStreetMap
wget http://download.openstreetmap.fr/extracts/oceania/australia/new_south_wales.osm.pbf
osmium tags-filter new_south_wales.osm.pbf w/highway -o nsw-roads.osm.pbf
osmium export -o nsw-roads.geojson nsw-roads.osm.pbf
ogr2ogr -f GPKG nsw-roads.gpkg nsw-roads.geojson
ogr2ogr -dialect OGRSQL -sql "SELECT geom FROM \"nsw-roads\" WHERE OGR_GEOMETRY='LineString'" nsw-roads.shp nsw-roads.gpkg

# in QGIS
# Vector > Geoprocessing Tools > Buffer
# buffer distance 0.0005

# Vector > Geoprocessing Tools > Difference
# Input Layer: AMS_Roads
# Overlay Layer: Buffered
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment