Esempio n. 1
0
    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')