Here are the examples of the python api numpy.testing.assert_arr_almost_eq taken from open source projects. By voting up you can indicate which examples are most useful and appropriate.
10 Examples
3
Example 1
Project: cartopy Source File: test_robinson.py
def test_transform_point():
# this way has always worked
result = _CRS_ROB.transform_point(35.0, 70.0, _CRS_PC)
assert_arr_almost_eq(result, (2376187.27182751, 7275317.81573085), _TOL)
# this always did something, but result has altered
result = _CRS_ROB.transform_point(_NAN, 70.0, _CRS_PC)
assert_true(np.all(np.isnan(result)))
# this used to crash + is now fixed
result = _CRS_ROB.transform_point(35.0, _NAN, _CRS_PC)
assert_true(np.all(np.isnan(result)))
3
Example 2
Project: cartopy Source File: test_crs.py
def test_osni(self):
osni = ccrs.OSNI()
ll = ccrs.Geodetic()
# results obtained by nearby.org.uk.
lat, lon = np.array([54.5622169298669, -5.54159863617957],
dtype=np.double)
east, north = np.array([359000, 371000], dtype=np.double)
assert_arr_almost_eq(osni.transform_point(lon, lat, ll),
np.array([east, north]),
-1)
assert_arr_almost_eq(ll.transform_point(east, north, osni),
np.array([lon, lat]),
3)
3
Example 3
Project: cartopy Source File: test_crs.py
def test_transform_points_nD(self):
rlons = np.array([[350., 352., 354.], [350., 352., 354.]])
rlats = np.array([[-5., -0., 1.], [-4., -1., 0.]])
src_proj = ccrs.RotatedGeodetic(pole_longitude=178.0,
pole_latitude=38.0)
target_proj = ccrs.Geodetic()
res = target_proj.transform_points(x=rlons, y=rlats,
src_crs=src_proj)
unrotated_lon = res[..., 0]
unrotated_lat = res[..., 1]
# Solutions derived by proj4 direct.
solx = np.array([[-16.42176094, -14.85892262, -11.90627520],
[-16.71055023, -14.58434624, -11.68799988]])
soly = np.array([[46.00724251, 51.29188893, 52.59101488],
[46.98728486, 50.30706042, 51.60004528]])
assert_arr_almost_eq(unrotated_lon, solx)
assert_arr_almost_eq(unrotated_lat, soly)
3
Example 4
Project: cartopy Source File: test_crs.py
def test_globe(self):
# Ensure the globe affects output.
rugby_globe = ccrs.Globe(semimajor_axis=9000000,
semiminor_axis=1000000)
footy_globe = ccrs.Globe(semimajor_axis=1000000,
semiminor_axis=1000000)
rugby_moll = ccrs.Mollweide(globe=rugby_globe)
footy_moll = ccrs.Mollweide(globe=footy_globe)
rugby_pt = rugby_moll.transform_point(10, 10, ccrs.Geodetic())
footy_pt = footy_moll.transform_point(10, 10, ccrs.Geodetic())
assert_arr_almost_eq(rugby_pt, (1400915, 1741319), decimal=0)
assert_arr_almost_eq(footy_pt, (155657, 193479), decimal=0)
3
Example 5
Project: cartopy Source File: test_crs.py
def test_project_point(self):
point = sgeom.Point([0, 45])
multi_point = sgeom.MultiPoint([point, sgeom.Point([180, 45])])
pc = ccrs.PlateCarree()
pc_rotated = ccrs.PlateCarree(central_longitude=180)
result = pc_rotated.project_geometry(point, pc)
assert_arr_almost_eq(result.xy, [[-180.], [45.]])
result = pc_rotated.project_geometry(multi_point, pc)
self.assertIsInstance(result, sgeom.MultiPoint)
self.assertEqual(len(result), 2)
assert_arr_almost_eq(result[0].xy, [[-180.], [45.]])
assert_arr_almost_eq(result[1].xy, [[0], [45.]])
0
Example 6
def test_transform_points():
# these always worked
result = _CRS_ROB.transform_points(_CRS_PC,
np.array([35.0]),
np.array([70.0]))
assert_arr_almost_eq(result,
[[2376187.27182751, 7275317.81573085, 0]], _TOL)
result = _CRS_ROB.transform_points(_CRS_PC,
np.array([35.0]),
np.array([70.0]),
np.array([0.0]))
assert_arr_almost_eq(result,
[[2376187.27182751, 7275317.81573085, 0]], _TOL)
# this always did something, but result has altered
result = _CRS_ROB.transform_points(_CRS_PC,
np.array([_NAN]),
np.array([70.0]))
assert_true(np.all(np.isnan(result)))
# this used to crash + is now fixed
result = _CRS_ROB.transform_points(_CRS_PC,
np.array([35.0]),
np.array([_NAN]))
assert_true(np.all(np.isnan(result)))
# multipoint case
x = np.array([10.0, 21.0, 0.0, 77.7, _NAN, 0.0])
y = np.array([10.0, _NAN, 10.0, 77.7, 55.5, 0.0])
z = np.array([10.0, 0.0, 0.0, _NAN, 55.5, 0.0])
expect_result = np.array(
[[9.40422591e+05, 1.06952091e+06, 1.00000000e+01],
[11.1, 11.2, 11.3],
[0.0, 1069520.91213902, 0.0],
[22.1, 22.2, 22.3],
[33.1, 33.2, 33.3],
[0.0, 0.0, 0.0]])
result = _CRS_ROB.transform_points(_CRS_PC, x, y, z)
assert_equal(result.shape, (6, 3))
assert_true(np.all(np.isnan(result[[1, 3, 4], :])))
result[[1, 3, 4], :] = expect_result[[1, 3, 4], :]
assert_false(np.any(np.isnan(result)))
assert_true(np.allclose(result, expect_result))
0
Example 7
Project: cartopy Source File: test_crs.py
def _check_osgb(self, osgb):
ll = ccrs.Geodetic()
# results obtained by streetmap.co.uk.
lat, lon = np.array([50.462023, -3.478831], dtype=np.double)
east, north = np.array([295131, 63511], dtype=np.double)
# note the handling of precision here...
assert_arr_almost_eq(np.array(osgb.transform_point(lon, lat, ll)),
np.array([east, north]),
1)
assert_arr_almost_eq(ll.transform_point(east, north, osgb),
[lon, lat],
2)
r_lon, r_lat = ll.transform_point(east, north, osgb)
r_inverted = np.array(osgb.transform_point(r_lon, r_lat, ll))
assert_arr_almost_eq(r_inverted, [east, north], 3)
r_east, r_north = osgb.transform_point(lon, lat, ll)
r_inverted = np.array(ll.transform_point(r_east, r_north, osgb))
assert_arr_almost_eq(r_inverted, [lon, lat])
0
Example 8
Project: cartopy Source File: test_crs.py
def test_transform_points_1D(self):
rlons = np.array([350., 352., 354., 356.])
rlats = np.array([-5., -0., 5., 10.])
src_proj = ccrs.RotatedGeodetic(pole_longitude=178.0,
pole_latitude=38.0)
target_proj = ccrs.Geodetic()
res = target_proj.transform_points(x=rlons, y=rlats,
src_crs=src_proj)
unrotated_lon = res[..., 0]
unrotated_lat = res[..., 1]
# Solutions derived by proj4 direct.
solx = np.array([-16.42176094, -14.85892262,
-12.88946157, -10.35078336])
soly = np.array([46.00724251, 51.29188893,
56.55031485, 61.77015703])
assert_arr_almost_eq(unrotated_lon, solx)
assert_arr_almost_eq(unrotated_lat, soly)
0
Example 9
Project: cartopy Source File: test_crs.py
def test_transform_points_xyz(self):
# Test geodetic transforms when using z value
rx = np.array([2574.32516e3])
ry = np.array([837.562e3])
rz = np.array([5761.325e3])
src_proj = ccrs.Geocentric()
target_proj = ccrs.Geodetic()
res = target_proj.transform_points(x=rx, y=ry, z=rz,
src_crs=src_proj)
glat = res[..., 0]
glon = res[..., 1]
galt = res[..., 2]
# Solution generated by pyproj
solx = np.array([18.0224043189])
soly = np.array([64.9796515089])
solz = np.array([5048.03893734])
assert_arr_almost_eq(glat, solx)
assert_arr_almost_eq(glon, soly)
assert_arr_almost_eq(galt, solz)
0
Example 10
Project: cartopy Source File: test_crs.py
def test_utm(self):
utm30n = ccrs.UTM(30)
ll = ccrs.Geodetic()
lat, lon = np.array([51.5, -3.0], dtype=np.double)
east, north = np.array([500000, 5705429.2], dtype=np.double)
assert_arr_almost_eq(utm30n.transform_point(lon, lat, ll),
[east, north],
decimal=1)
assert_arr_almost_eq(ll.transform_point(east, north, utm30n),
[lon, lat],
decimal=1)
utm38s = ccrs.UTM(38, southern_hemisphere=True)
lat, lon = np.array([-18.92, 47.5], dtype=np.double)
east, north = np.array([763316.7, 7906160.8], dtype=np.double)
assert_arr_almost_eq(utm38s.transform_point(lon, lat, ll),
[east, north],
decimal=1)
assert_arr_almost_eq(ll.transform_point(east, north, utm38s),
[lon, lat],
decimal=1)