def test_local_distances(self):
for reflat, reflon in [(0.0, 0.0), (10.0, 10.0), (90.0, 0.0),
(-90.0, 0.0), (0.0, 180.0), (0.0, -180.0),
(90.0, 180.0)]:
north, east = serialgrid(num.linspace(-10 * km, 10 * km, 21),
num.linspace(-10 * km, 10 * km, 21))
lat, lon = orthodrome.ne_to_latlon2(reflat, reflon, north, east)
north2, east2 = orthodrome.latlon_to_ne2(reflat, reflon, lat, lon)
dist1 = num.sqrt(north**2 + east**2)
dist2 = num.sqrt(north2**2 + east2**2)
dist3 = orthodrome.distance_accurate15nm(reflat, reflon, lat, lon)
assert num.all(num.abs(dist1 - dist2) < 0.0001)
assert num.all(num.abs(dist1 - dist3) < 0.0001)
def test_local_distances(self):
for reflat, reflon in [
(0.0, 0.0),
(10.0, 10.0),
(90.0, 0.0),
(-90.0, 0.0),
(0.0, 180.0),
(0.0, -180.0),
(90.0, 180.0) ]:
north, east = serialgrid(num.linspace(-10*km, 10*km, 21),
num.linspace(-10*km, 10*km, 21))
lat, lon = orthodrome.ne_to_latlon2(reflat, reflon, north, east)
north2, east2 = orthodrome.latlon_to_ne2(reflat, reflon, lat, lon)
dist1 = num.sqrt(north**2 + east**2)
dist2 = num.sqrt(north2**2 + east2**2)
dist3 = orthodrome.distance_accurate15nm(reflat, reflon, lat, lon)
assert num.all(num.abs(dist1-dist2) < 0.0001)
assert num.all(num.abs(dist1-dist3) < 0.0001)