/usr/share/sumo/tools/osmGet.py is in sumo-tools 0.25.0+dfsg1-2.
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 | #!/usr/bin/env python
"""
@file osmGet.py
@author Daniel Krajzewicz
@author Jakob Erdmann
@author Michael Behrisch
@date 2009-08-01
@version $Id: osmGet.py 19535 2015-12-05 13:47:18Z behrisch $
Retrieves an area from OpenStreetMap.
SUMO, Simulation of Urban MObility; see http://sumo.dlr.de/
Copyright (C) 2009-2015 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.
"""
import os
import sys
import httplib
import StringIO
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), "w")
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()
|