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)