/usr/share/sumo/tools/osmGet.py is in sumo-tools 0.28.0+dfsg1-1.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 | #!/usr/bin/env python
"""
@file osmGet.py
@author Daniel Krajzewicz
@author Jakob Erdmann
@author Michael Behrisch
@date 2009-08-01
@version $Id: osmGet.py 20687 2016-05-10 11:27:00Z behrisch $
Retrieves an area from OpenStreetMap.
SUMO, Simulation of Urban MObility; see http://sumo.dlr.de/
Copyright (C) 2009-2016 DLR (http://www.dlr.de/) and contributors
This file is part of SUMO.
SUMO is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 3 of the License, or
(at your option) any later version.
"""
from __future__ import absolute_import
from __future__ import print_function
import os
import sys
try:
import httplib
except ImportError:
import http.client as httplib
import gzip
import optparse
from os import path
import sumolib
def readCompressed(conn, query, filename):
conn.request("POST", "/api/interpreter", """
<osm-script timeout="240" element-limit="1073741824">
<union>
%s
<recurse type="node-relation" into="rels"/>
<recurse type="node-way"/>
<recurse type="way-relation"/>
</union>
<union>
<item/>
<recurse type="way-node"/>
</union>
<print mode="body"/>
</osm-script>""" % query)
response = conn.getresponse()
print(response.status, response.reason)
if response.status == 200:
out = open(path.join(os.getcwd(), filename), "wb")
out.write(response.read())
out.close()
optParser = optparse.OptionParser()
optParser.add_option("-p", "--prefix", default="osm", help="for output file")
optParser.add_option(
"-b", "--bbox", help="bounding box to retrieve in geo coordinates south,west,north,east")
optParser.add_option("-o", "--oldapi", action="store_true",
default=False, help="use old API for retrieval")
optParser.add_option("-t", "--tiles", type="int",
default=1, help="number of tiles when using old api")
optParser.add_option(
"-d", "--output-dir", help="optional output directory (must already exist)")
optParser.add_option("-a", "--area", type="int", help="area id to retrieve")
optParser.add_option(
"-x", "--polygon", help="calculate bounding box from polygon data in file")
def get(args=None):
(options, args) = optParser.parse_args(args=args)
if not options.bbox and not options.area and not options.polygon:
optParser.error(
"At least one of 'bbox' and 'area' and 'polygon' has to be set.")
if options.oldapi and options.area:
optParser.error("Only the new API supports 'area'.")
if options.polygon:
west = 1e400
south = 1e400
east = -1e400
north = -1e400
for area in sumolib.output.parse_fast(options.polygon, 'poly', ['shape']):
coordList = [tuple(map(float, x.split(',')))
for x in area.shape.split()]
for point in coordList:
west = min(point[0], west)
south = min(point[1], south)
east = max(point[0], east)
north = max(point[1], north)
if options.bbox:
south, west, north, east = [float(v) for v in options.bbox.split(',')]
if south > north or west > east:
optParser.error("Invalid geocoordinates in bbox.")
if options.output_dir:
options.prefix = path.join(options.output_dir, options.prefix)
if options.oldapi:
num = options.tiles
b = west
conn = httplib.HTTPConnection("api.openstreetmap.org")
for i in range(num):
e = b + (east - west) / float(num)
req = "/api/0.6/map?bbox=%s,%s,%s,%s" % (b, south, e, north)
conn.request("GET", req)
r = conn.getresponse()
print(req, r.status, r.reason)
out = open(
path.join(os.getcwd(), "%s%s_%s.osm.xml" % (options.prefix, i, num)), "w")
out.write(r.read())
out.close()
b = e
conn.close()
else:
conn = httplib.HTTPConnection("www.overpass-api.de")
#conn = httplib.HTTPConnection("overpass.osm.rambler.ru")
if options.area:
if options.area < 3600000000:
options.area += 3600000000
readCompressed(conn, '<area-query ref="%s"/>' %
options.area, options.prefix + "_city.osm.xml")
if options.bbox or options.polygon:
if options.tiles == 1:
readCompressed(conn, '<bbox-query n="%s" s="%s" w="%s" e="%s"/>' %
(north, south, west, east), options.prefix + "_bbox.osm.xml")
else:
num = options.tiles
b = west
for i in range(num):
e = b + (east - west) / float(num)
readCompressed(conn, '<bbox-query n="%s" s="%s" w="%s" e="%s"/>' % (
north, south, b, e), "%s%s_%s.osm.xml" % (options.prefix, i, num))
b = e
conn.close()
if __name__ == "__main__":
get()
|