def test_helpers(self):
p = gpstk.cartesian(100, 200, 300)
self.assertEqual(gpstk.Position.Cartesian, p.getCoordinateSystem())
p = gpstk.spherical(45, 60, 100000, model=gpstk.PZ90Ellipsoid())
self.assertEqual(gpstk.Position.Spherical, p.getCoordinateSystem())
p = gpstk.geodetic(frame=gpstk.ReferenceFrame('WGS84'))
self.assertEqual(gpstk.Position.Geodetic, p.getCoordinateSystem())
p = gpstk.geocentric(latitude=60, radius=10000)
self.assertEqual(gpstk.Position.Geocentric, p.getCoordinateSystem())
def test_functions(self):
system = gpstk.Position.Cartesian
ell = gpstk.PZ90Ellipsoid()
frame = gpstk.ReferenceFrame(gpstk.ReferenceFrame.PZ90)
p = gpstk.Position(10000.0, 150000.0, 200000.0, system, ell, frame)
q = gpstk.Position(20000.0, 160000.0, 190000.0, system, ell, frame)
self.assertAlmostEqual(1.32756277187, q.elevation(p))
self.assertAlmostEqual(86.18592516570916, p.getPhi())
self.assertAlmostEqual(57.5141089193572, p.geodeticLatitude())
self.assertAlmostEqual(10000.0, p.X())
def test(self):
data = gpstk.Xvt()
data.x = gpstk.Triple(1000.0, 2000.0, 1500.0)
data.v = gpstk.Triple(50.0, 25.0, -500.0)
data.clkbias = 0.0001
data.clkdrift = 0.05
data.relcorr = 0.83
data.frame = gpstk.ReferenceFrame(gpstk.ReferenceFrame.WGS84)
self.assertAlmostEqual(0.0001, data.getClockBias())
expected = 1.446445072869704e-11
self.assertAlmostEqual(expected, data.computeRelativityCorrection())
expected = ('x:(1000, 2000, 1500), v:(50, 25, -500), clk bias:0.0001, '
'clk drift:0.05, relcorr:1.44645e-11')
self.assertEqual(expected, str(data))
def test_unknown(self):
r = gpstk.ReferenceFrame()
self.assertEqual('Unknown', str(r))
def test_constant_input(self):
r = gpstk.ReferenceFrame(gpstk.ReferenceFrame.WGS84)
self.assertEqual('WGS84', str(r))
def test_string_input(self):
r = gpstk.ReferenceFrame('PZ90')
self.assertEqual('PZ90', str(r))
