def setUp(self):
loader = nhpn.Loader() # Loading the data takes a few
# seconds. This statement can be moved
# to just below "import nhpn" if you
# don't modify nodes or links in your
# code.
self.nodes = loader.nodes()
self.links = loader.links()
def flat_earth(filename='path_flat.kml'):
"""Test long path"""
loader = nhpn.Loader()
nodes = loader.nodes()
links = loader.links()
source = dijkstra.node_by_name(nodes, 'CAMDEN', 'NJ')
destination = dijkstra.node_by_name(nodes, 'SAN DIEGO', 'CA')
ans = dijkstra.shortest_path(nodes, links, dijkstra.distance, source,
destination)
nhpn.Visualizer.toKML(ans, 'path_flat.kml')
def curved_earth(filename='path_curved.kml'):
"""Test long path"""
loader = nhpn.Loader()
nodes = loader.nodes()
links = loader.links()
def distance(node1, node2):
"""Returns the distance between node1 and node2, including the
Earth's curvature."""
A = node1.latitude * pi / 10**6 / 180
B = node1.longitude * pi / 10**6 / 180
C = node2.latitude * pi / 10**6 / 180
D = node2.longitude * pi / 10**6 / 180
return acos(sin(A) * sin(C) + cos(A) * cos(C) * cos(B - D))
source = dijkstra.node_by_name(nodes, 'CAMDEN', 'NJ')
destination = dijkstra.node_by_name(nodes, 'SAN DIEGO', 'CA')
ans = dijkstra.shortest_path(nodes, links, distance, source, destination)
nhpn.Visualizer.toKML(ans, 'path_curved.kml')