def test_rad():
"""Test rad function"""
assert geo.rad(90) == approx(np.pi/2.)
assert geo.rad(4, 30) == approx(np.pi/40.)
assert geo.rad(1, 7, 30) == approx(np.pi/160.)
assert geo.rad(np.array([90, 45])) == approx(np.array((np.pi/2.,
np.pi/4.)))
def test_WGS84toOSGB36():
"""Test WGS84toOSGB36."""
lat_long_wgs = np.array([[geo.rad(52, 39, 28.71)], [geo.rad(1, 42,
57.79)]])
lat_long_os = np.array([[geo.rad(52, 39, 27.2531)],
[geo.rad(1, 43, 4.5177)]])
assert np.array(geo.WGS84toOSGB36(*lat_long_wgs,
True)) == approx(lat_long_os, rel=1.0e-5)
def test_lat_long_to_xyz():
"""Test lat_long_to_xyz function."""
latitude = np.array([geo.rad(52, 39, 27.2531)])
longitude = np.array([geo.rad(1, 43, 4.5177)])
xyz = np.array([[3874938.849], [116218.624], [5047168.208]])
datum = copy.copy(geo.osgb36)
datum.F_0 = 1.0
assert geo.lat_long_to_xyz(latitude, longitude, True,
datum=datum) == approx(xyz, rel=1.0e-8)
def test_get_easting_northing_from_lat_long():
"""Test lat_long_to_easting_northing."""
latitude = np.array([geo.rad(52, 39, 27.2531)])
longitude = np.array([geo.rad(1, 43, 4.5177)])
e = 651409.903
n = 313177.270
assert np.array(geo.get_easting_northing_from_lat_long(latitude,
longitude, radians=True)) \
== approx(np.array((e,n)),
rel=1.0e-5)
def test_xyz_to_lat_long():
"""Test xyz_to_lat_long function."""
x = np.array([3874938.849])
y = np.array([116218.624])
z = np.array([5047168.208])
latitude = np.array([geo.rad(52, 39, 27.2531)])
longitude = np.array([geo.rad(1, 43, 4.5177)])
assert np.array(geo.xyz_to_lat_long(x, y, z, True)) \
== approx(np.array((latitude, longitude)), rel=1.0e-5)
