def get_dist(n1, n2, cheat=False):
if cheat:
return haversine(n1, n2) * 1250
data = LoadOsm("foot")
node1 = data.findNode(*n1)
node2 = data.findNode(*n2)
router = Router(data)
result, route = router.doRoute(node1, node2)
return 1000 * calk_route_dist(route, data)
def routeToGpx(lat1,
lon1,
lat2,
lon2,
transport,
description="",
style="track"):
"""Format a route (as list of nodes) into a GPX file"""
data = LoadOsm(transport)
node1 = data.findNode(lat1, lon1)
node2 = data.findNode(lat2, lon2)
router = Router(data)
result, route = router.doRoute(node1, node2)
if result != 'success':
return
output = ''
output = output + "<?xml version='1.0'?>\n"
output = output + "<gpx version='1.1' creator='pyroute' xmlns='http://www.topografix.com/GPX/1/1' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xsi:schemaLocation='http://www.topografix.com/GPX/1/1 http://www.topografix.com/GPX/1/1/gpx.xsd'>\n"
if (style == 'track'):
output = output + " <trk>\n"
output = output + " <name>%s</name>\n" % description
output = output + " <trkseg>\n"
count = 0
for i in route:
node = data.rnodes[i]
output = output + " <trkpt lat='%f' lon='%f'>\n" % ( \
node[0],
node[1])
output = output + " </trkpt>\n"
count = count + 1
output = output + " </trkseg>\n </trk>\n</gpx>\n"
elif (style == 'route'):
output = output + " <rte>\n"
output = output + " <name>%s</name>\n" % description
count = 0
for i in route:
node = data.rnodes[i]
output = output + " <rtept lat='%f' lon='%f'>\n" % ( \
node[0],
node[1])
output = output + " <name>%d</name>\n" % count
output = output + " </rtept>\n"
count = count + 1
output = output + " </rte>\n</gpx>\n"
return (output)
def routeToCSV(lat1, lon1, lat2, lon2, transport):
"""Format a route (as list of nodes)"""
data = LoadOsm(transport)
node1 = data.findNode(lat1, lon1)
node2 = data.findNode(lat2, lon2)
router = Router(data)
result, route = router.doRoute(node1, node2)
if result != 'success':
return ("Fail")
output = ''
for i in route:
node = data.rnodes[i]
output = output + "%d,%f,%f\n" % ( \
i,
node[0],
node[1])
return (output)
52.282673935069106,
104.28143367544139,
)
n2 = (
52.282374539713466,
104.28054854646648,
)
n4 = (
52.281596431286715,
104.27944884077036,
)
print(get_dist(n1, n2))
print(haversine(n1, n2) * 1250)
node1 = data.findNode(*n1) #52.282673935069106, 104.28143367544139
node2 = data.findNode(*n4) #52.282374539713466, 104.28054854646648
# print(node1)
# print(node2)
router = Router(data)
result, route = router.doRoute(node1, node2)
if result == 'success':
# list the nodes
print("ROUTE")
print("---------------------------")
print(route)
print("---------------------------")
# print("Points")
def route_geojson(input_f, output_f, mode='foot', local_planet=None):
osmdata = LoadOsm(mode)
if local_planet != None:
osmdata.getArea = lambda lat, lon: None
osmdata.api = None
print('loading osm data (this may take a while)...')
osmdata.loadOsm(local_planet)
print('starting router...')
router = Router(osmdata)
print('processing shapes...')
# First load up the shapes
layer = geojson.load(input_f)
non_linestring = 0
not_two_points = 0
unsuccessful = 0
successful = 0
very_long = 0
first = True
output_f.write('{"crs": {"type": "name", "properties": {"name": "urn:ogc:def:crs:OGC:1.3:CRS84"}}, "type": "FeatureCollection", "features": [\n')
for feature in layer.features:
if feature.geometry.type != 'LineString':
# Not a LineString, skip!
non_linestring += 1
continue
geometry = list(feature.geometry.coordinates)
if len(geometry) != 2:
# LineString with other than two points, skip!
not_two_points += 1
continue
if pythagoras(*geometry[0] + geometry[1]) > 1.0:
very_long += 1
continue
# Now find a route. Data has x,y coordinates, but function is y,x, so
# reverse the parameters.
start = osmdata.findNode(*geometry[0][::-1])
end = osmdata.findNode(*geometry[1][::-1])
result, route = router.doRoute(start, end)
if result != 'success':
unsuccessful += 1
continue
routed_geometry = []
for node_id in route:
node = osmdata.rnodes[node_id]
routed_geometry.append((node[1], node[0]))
new_feature = geojson.Feature(
geometry=geojson.LineString(coordinates=routed_geometry),
properties=feature.properties,
id=feature.id,
)
if not first:
output_f.write(',\n')
first = False
geojson.dump(new_feature, output_f)
output_f.flush()
successful += 1
output_f.write('\n]}\n')
output_f.close()
print('%d LineStrings routed. Errors: %d non-linestring(s), %d linestring(s) with !=2 points, %d very long, %d unsuccessful routings' % (successful, non_linestring, not_two_points, very_long, unsuccessful))
#!/bin/env python3
from pyroutelib2.route import Router
from pyroutelib2.loadOsm import LoadOsm
osm = LoadOsm("foot")
#home = osm.findNode(54.355561, 10.131721)
home = osm.findNode(54.357069, 10.130382)
#finance = osm.findNode(54.344639, 10.132518)
finance = osm.findNode(54.310372, 10.131111)
print(home)
print(finance)
router = Router(osm)
result, route = router.doRoute(home, finance)
if result == 'success':
for i in route:
node = osm.rnodes[i]
print("%d: %f,%f" % (i,node[0],node[1]))
else:
print('Failed to get route: ' + result)
# Try to insert, keeping the queue ordered by decreasing worst-case distance
count = 0
for test in self.queue:
if test['maxdistance'] > queueItem['maxdistance']:
self.queue.insert(count,queueItem)
break
count = count + 1
else:
self.queue.append(queueItem)
if __name__ == "__main__":
# Test suite - do a little bit of easy routing in birmingham
data = LoadOsm("cycle")
node1 = data.findNode(52.552394,-1.818763)
node2 = data.findNode(52.563368,-1.818291)
print(node1)
print(node2)
router = Router(data)
result, route = router.doRoute(node1, node2)
if result == 'success':
# list the nodes
print(route)
# list the lat/long
for i in route:
node = data.rnodes[i]
print("%d: %f,%f" % (i,node[0],node[1]))