def testDms(self): # dms module tests self.test('parseDMS', parseDMS( '0.0°'), '0.0') self.test('parseDMS', parseDMS( '0°'), '0.0') self.test('parseDMS', parseDMS('''000°00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00'00.0"'''), '0.0') self.test('parseDMS', parseDMS('''000° 00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00 ' 00.0"'''), '0.0') x = parse3llh('000° 00′ 05.31″W, 51° 28′ 40.12″ N') x = ', '.join('%.6f' % a for a in x) # XXX fStr self.test('parse3llh', x, '51.477811, -0.001475, 0.000000') for a, x in (((), '''45°45'45.36"'''), ((F_D, None), '45.7626°'), ((F_DM, None), "45°45.756'"), ((F_DMS, None), '''45°45'45.36"'''), ((F_D, 6), '45.7626°'), ((F_DM, -4), "45°45.7560'"), ((F_DMS, 2), '''45°45'45.36"''')): self.test('toDMS', toDMS(45.76260, *a), x) for a, x in (((1,), 'N'), ((0,), 'N'), ((-1,), 'N'), ((359,), 'N'), ((24,), 'NNE'), ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((226,), 'SW'), ((226, 1), 'W'), ((226, 2), 'SW'), ((226, 3), 'SW'), ((237,), 'WSW'), ((237, 1), 'W'), ((237, 2), 'SW'), ((237, 3), 'WSW')): self.test('compassPoint', compassPoint(*a), x)
def testEcef(self, Ecef): self.test(Ecef.__name__, '...', '...', nl=1) Karney = Ecef is EcefKarney Sudano = Ecef is EcefSudano g = Ecef(Datums.WGS84, name='Test') self.test('name', g.name, 'Test') t = g.toStr2() self.test('toStr', t, g.classname, known=True) t = Ecef(g.a, g.f, name=g.name) # coverage self.test('a, f', t, t.classname, known=True) self.testCopy(g) # <https://GeographicLib.SourceForge.io/html/classGeographicLib_1_1Geocentric.html> t = g.forward(27.99, 86.93, 8820) # Mt Everest self.test('forward', fstr(t[3:6], prec=2), '27.99, 86.93, 8820.0') self.test('forward', fstr(t[0:3], prec=1), '302271.4, 5635928.4, 2979666.1') self.test('name', t.name, 'Test') t = g.reverse(302271.4, 5635928.4, 2979666.1) self.test('reverse', fstr(t[0:3], prec=1), '302271.4, 5635928.4, 2979666.1') self.test('reverse', fstr(t[3:6], prec=2), '27.99, 86.93, 8820.01', known=Sudano and _known(t, 27.99, 8820)) self.test('case', t.C, 2 if Karney else (5 if Sudano else 1)) self.test('iteration', t.iteration, t.iteration) self.test('name', t.name, 'Test') # <https://GeographicLib.SourceForge.io/html/classGeographicLib_1_1Geocentric.html> t = g.reverse(302e3, 5636e3, 2980e3) self.test('reverse', fstr(t[0:3], prec=1), '302000.0, 5636000.0, 2980000.0') self.test('reverse', fstr(t[3:6], prec=2), '27.99, 86.93, 9027.03', known=Sudano and _known(t, 27.99, 9027)) self.test('case', t.C, 2 if Karney else (5 if Sudano else 1)) self.test('iteration', t.iteration, t.iteration) t = g.forward(27.99, 86.93, 8820.0) self.test('forward', fstr(t[3:6], prec=2), '27.99, 86.93, 8820.0') self.test('forward', fstr(t[0:3], prec=2), '302271.43, 5635928.37, 2979666.13') # <https://GeographicLib.SourceForge.io/html/CartConvert.1.html> t = g.forward(33.3, 44.4, 6000) self.test('forward', fstr(t[3:6], prec=2), '33.3, 44.4, 6000.0') self.test('forward', fstr(t[0:3], prec=2), '3816209.6, 3737108.55, 3485109.57') t = g.reverse(3816209.6, 3737108.55, 3485109.57) self.test('reverse', fstr(t[0:3], prec=2), '3816209.6, 3737108.55, 3485109.57') self.test('reverse', fstr(t[3:6], prec=3), '33.3, 44.4, 5999.996', known=Sudano and _known(t, 33.3, 5999)) self.test('case', t.C, 2 if Karney else (5 if Sudano else 1)) self.test('iteration', t.iteration, t.iteration) # <https://GeographicLib.SourceForge.io/html/CartConvert.1.html> t = g.reverse(30000, 30000, 0) self.test('reverse', fstr(t[0:3], prec=1), '30000.0, 30000.0, 0.0') self.test('reverse', fstr(t[3:6], prec=3), '6.483, 45.0, -6335709.726', known=not Karney) self.test('case', t.C, 3 if Karney else (0 if Sudano else 1)) self.test('iteration', t.iteration, t.iteration) t = g.forward(6.483, 45.0, -6335709.726) self.test('forward', fstr(t[3:6], prec=3), '6.483, 45.0, -6335709.726') self.test('forward', fstr(t[0:3], prec=1), '30000.0, 30000.0, -0.0', known=True) # Rey-Jer You <https://www.ResearchGate.net/publication/240359424> for i, (x, y, z, h) in enumerate( ((-2259148.993, 3912960.837, 4488055.516, 1e3), (-2259502.546, 3913573.210, 4488762.622, 2e3), (-2259856.100, 3914185.582, 4489469.729, 3e3), (-2260209.653, 3914797.955, 4490176.836, 4e3), (-2262330.973, 3918472.189, 4494419.477, 1e4), (-2265866.507, 3924595.914, 4501490.544, 2e4), (-2294150.778, 3973585.709, 4558059.087, 1e5), (-2541638.152, 4402246.414, 5053033.834, 8e5), (-2612348.830, 4524720.901, 5194455.190, 1e6))): i = '-%d' % (i + 1, ) r = '45.0, 120.0, %.1f' % (h, ) # Zero and First order columns t = g.reverse(x, y, z) k = Sudano and _known(t, 45, h) self.test('reverse' + i, fstr(t[3:6], prec=3), r, known=k) f = g.forward(t.lat, t.lon, t.height) self.test('forward' + i, fstr(f[0:3], prec=1), fstr((x, y, z), prec=1), known=k) self.test('xyzh' + i, fstr(f.xyzh, prec=1), fstr((x, y, z, h), prec=1), known=k) f = f.phi, f.lam t = radians(t.lat), radians(t.lon) self.test('philam' + i, fstr(f, prec=4), fstr(t, prec=4)) # <https://www.researchgate.net/publication/3709199_An_exact_conversion_from_an_Earth-centered_coordinate_system_to_latitude_longitude_and_altitude> t = g.reverse(4588301.55696757, 0, 4558059.086984613) self.test('sudano', fstr(t[3:6], prec=3), '45.0, 0.0, 100000.0', known=Sudano) # <https://www.OrdnanceSurvey.co.UK/documents/resources/guide-coordinate-systems-great-britain.pdf> pp 47 g = Ecef(Ellipsoids.GRS80, name='OS-UK') self.test('name', g.name, 'OS-UK') t = g.forward( parse3llh('''53°36′43.1653"N, 001°39′51.9920"W, 299.800''')) self.test('forward', fstr(t[3:6], prec=8), '53.61199036, -1.66444222, 299.8') self.test('forward', fstr(t[0:3], prec=2), '3790644.9, -110149.21, 5111482.97') t = g.reverse(3790644.9, -110149.21, 5111482.97) self.test('reverse', fstr(t[0:3], prec=2), '3790644.9, -110149.21, 5111482.97') self.test('reverse', fstr(t[3:5], prec=8), '53.61199036, -1.66444223', known=Sudano) self.test('reverse.lat', latDMS(t.lat, prec=4), '53°36′43.1653″N', known=Sudano) self.test('reverse.lon', lonDMS(t.lon, prec=4), '001°39′51.992″W') self.test('reverse.height', fstr(t.height, prec=-3), '299.800', known=Sudano) self.test('case', t.C, 2 if Karney else (6 if Sudano else 1)) self.test('iteration', t.iteration, t.iteration) # coverage self.test(EcefError.__name__, g.reverse(0, 0, 0), '(0.0, 0.0, ...)', known=True) try: self.test(EcefError.__name__, g.forward(None, None, None), EcefError.__name__) except Exception as x: self.test(EcefError.__name__, str(x).split(':')[0], 'lat (None), lon (None) ...', known=True) try: self.test(Ecef.__name__, Ecef(None), EcefError.__name__) except Exception as x: self.test(Ecef.__name__, str(x), Ecef.__name__, known=True)
def testDms(self): # MCCABE 14 for t in ('0.0°', '0°', '''000°00'00"''', '''000°00'00.0"''', '''000° 00'00"''', '''000°00 ' 00.0"''', '''000° 00' 00.0'''): self.test("parseDMS('%s')" % (t, ), parseDMS(t), '0.0') t = '000°-00′-00.0"' self.test("parseDMS('%s')" % (t, ), parseDMS(t, sep='-'), '0.0') _X = '' for t, X, x in ((1, '1.0', _X), (12, '12.0', _X), (123, '123.0', _X), (1234, '12.567', '1234.0'), (12345, '123.75', '12345.0'), (123456, '12.582', '123456.0'), (1234567, '123.769', '1234567.0'), (12345678, '1234.955', '12345678.0'), (.1, 0.1, _X), (1.2, 1.2, _X), (12.3, 12.3, _X), (123.4, 123.4, _X), (1234.5, 12.575, 1234.5), (12345.6, 123.76, 12345.6), (123456.7, 12.582, 123456.7), ('1N', '1.0', _X), ('12S', '-12.0', _X), ('012.3W', '-12.3', _X), ('123E', '123.0', _X), ('1234N', '12.567', '1234.0'), ('12345E', '123.75', '12345.0'), ('1234.5S', '-12.575', '-1234.5'), ('12345.6E', '123.76', '12345.6'), ('123456.7S', '-12.582', '-123456.7'), ('1234567.8W', '-123.769', '-1234567.8'), ('12345678E', '12345678.0', _X)): self.test('parseDDDMMSS(%r)' % (t, ), fstr(parseDDDMMSS(t), prec=3), X) self.test('parseDMS(%r)' % (t, ), fstr(parseDMS(t), prec=3), x or X) t = 345.0 self.test('parseDDDMMSS(%r, NS)' % (t, ), fstr(parseDDDMMSS(t, suffix='NS'), prec=3), '3.75') self.test('parseDDDMMSS(%r, EW)' % (t, ), fstr(parseDDDMMSS(t, suffix='EW'), prec=3), '345.0') t = 5430.0 self.test('parseDDDMMSS(%r, NS)' % (t, ), fstr(parseDDDMMSS(t, suffix='NS'), prec=3), '54.5') self.test('parseDDDMMSS(%r, EW)' % (t, ), fstr(parseDDDMMSS(t, suffix='EW'), prec=3), '54.5') for t, x in (('12E', '12.0'), ('012.3S', '-12.3'), ('123N', '123.0'), ('1234E', '1234.0'), ('12345N', '12345.0'), ('1234.5W', '-1234.5'), ('123456E', '123456.0'), ('1234567S', '-1234567.0')): try: self.test('parseDDDMMSS(%r)' % (t, ), parseDDDMMSS(t), ParseError.__name__) except ParseError as X: self.test('parseDDDMMSS(%r)' % (t, ), repr(X), repr(X)) self.test('parseDMS(%r)' % (t, ), fstr(parseDMS(t), prec=5), x) r = rangerrors(True) try: self.test('parseDMS', parseDMS(181, clip=180), RangeError.__name__) except RangeError as x: self.test('parseDMS', str(x), str(x)) rangerrors(False) try: self.test('parseDMS', parseDMS(-91, clip=90), RangeError.__name__) except RangeError as x: self.test('parseDMS', str(x), str(x)) rangerrors(r) x = parse3llh('000° 00′ 05.31″W, 51° 28′ 40.12″ N') self.test('parse3llh', fstr(x, prec=6), '51.477811, -0.001475, 0.0') t = 'toDMS(%s)' % ('', ) self.test(t, toDMS(45.99999, F_DM, prec=1), '46°00.0′') # not 45°60.0′ self.test(t, toDMS(45.99999, F_DM, prec=2), '46°00.0′') self.test(t, toDMS(45.9999, F_DM, prec=2), '45°59.99′') self.test(t, toDMS(45.99999, F_DM, prec=3), '45°59.999′') self.test(t, toDMS(45.99999, F_DMS, prec=1), '46°00′00.0″') self.test(t, toDMS(45.99999, F_DMS, prec=2), '45°59′59.96″') self.test(t, toDMS(45.99999, F_DMS, prec=3), '45°59′59.964″') self.test(t, toDMS(45.76260), '45°45′45.36″') for F, p, x in ((F_D, None, '45.7626°'), (F_DM, None, "45°45.756'"), (F_DMS, None, '''45°45'45.36"'''), (F_DEG, None, '45.7626'), (F_MIN, None, '4545.756'), (F_SEC, None, '454545.36'), (F_RAD, None, '0.79871'), (F_D, 6, '45.7626°'), (F_DM, -4, "45°45.7560'"), (F_DMS, 2, '''45°45'45.36"'''), (F_DEG, -6, '45.762600'), (F_MIN, -5, '4545.75600'), (F_SEC, -3, '454545.360'), (F_RAD, -6, '0.798708'), (F_D_, 6, '45.7626°'), (F_DM_, -4, "45°45.7560'"), (F_DMS_, 2, '''45°45'45.36"'''), (F_DEG_, -6, '45.762600'), (F_MIN_, -5, '4545.75600'), (F_SEC_, -3, '454545.360'), (F_RAD_, -6, '0.798708')): t = 'toDMS(%s)' % (F, ) self.test(t, toDMS(45.76260, F, prec=p), x) self.test(t, toDMS(-45.76260, F, prec=p), '-' + x) for F, p, x in ((F_D__, 6, '45.7626°'), (F_DM__, -4, "45°45.7560'"), (F_DMS__, 2, '''45°45'45.36"'''), (F_DEG__, -6, '45.762600'), (F_MIN__, -5, '4545.75600'), (F_SEC__, -3, '454545.360'), (F_RAD__, -6, '0.798708')): t = 'toDMS(%s)' % (F, ) self.test(t, toDMS(45.76260, F, prec=p), '+' + x) self.test(t, toDMS(-45.76260, F, prec=p), '-' + x) # <https://GitHub.com/chrisveness/geodesy/blob/master/test/dms-tests.js> for a, x in ( ((1, ), 'N'), ((0, ), 'N'), ((-1, ), 'N'), ((359, ), 'N'), ((24, ), 'NNE'), ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((226, ), 'SW'), ((226, 1), 'W'), ((226, 2), 'SW'), ((226, 3), 'SW'), ((237, ), 'WSW'), ((237, 1), 'W'), ((237, 2), 'SW'), ((237, 3), 'WSW'), # Martin Schultz ((11.25, ), 'NNE'), ((11.249, ), 'N'), ((-11.25, ), 'N'), ((348.749, ), 'NNW'), ((45, 1), 'E'), ((44.99, 1), 'N'), ((45, 2), 'NE'), ((44.99, 2), 'NE'), ((45, 3), 'NE'), ((44.99, 3), 'NE'), ((45, 4), 'NE'), ((44.99, 4), 'NE'), # XXX ((22.5, 1), 'N'), ((22.49, 1), 'N'), ((22.5, 2), 'NE'), ((22.49, 2), 'N'), ((22.5, 3), 'NNE'), ((22.49, 3), 'NNE'), ((22.5, 4), 'NNE'), ((22.49, 4), 'NNE'), # XXX ((11.25, 1), 'N'), ((11.249, 1), 'N'), ((11.25, 2), 'N'), ((11.249, 2), 'N'), ((11.25, 3), 'NNE'), ((11.249, 3), 'N'), ((11.25, 4), 'NbE'), ((11.249, 4), 'NbE'), # XXX # examples ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((24, ), 'NNE'), ((18, 3), 'NNE'), ((11, 4), 'NbE'), ((30, 4), 'NEbN')): self.test('compassPoint%s' % (a, ), compassPoint(*a), x) for a, x in enumerate( ('NbE', 'NEbN', 'NEbE', 'EbN', 'EbS', 'SEbE', 'SEbS', 'SbE', 'SbW', 'SWbS', 'SWbW', 'WbS', 'WbN', 'NWbW', 'NWbN', 'NbW')): a = 11.25 + a * 22.5 self.test('compassPoint(%s)' % (a, ), compassPoint(a, 4), x) t = degDMS(1.0101, prec=5, s_D='', pos='+') # coverage self.test('_DEG', t, '+1.0101') t = degDMS(0.0101, prec=5, s_S='', pos='+') self.test('_MIN', t, "+0.606'") t = degDMS(0.0101, prec=5, s_M='', pos='+') self.test('_SEC', t, '+36.36"')
def testDms(self): # dms module tests self.test('parseDMS', parseDMS('0.0°'), '0.0') self.test('parseDMS', parseDMS('0°'), '0.0') self.test('parseDMS', parseDMS('''000°00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00'00.0"'''), '0.0') self.test('parseDMS', parseDMS('''000° 00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00 ' 00.0"'''), '0.0') r = rangerrors(True) try: self.test('parseDMS', parseDMS(181, clip=180), 'ValueError') except ValueError as x: self.test('parseDMS', str(x), '181.0 beyond 180 degrees') rangerrors(False) try: self.test('parseDMS', parseDMS(-91, clip=90), -90) except ValueError as x: self.test('parseDMS', str(x), '-90') rangerrors(r) x = parse3llh('000° 00′ 05.31″W, 51° 28′ 40.12″ N') x = ', '.join('%.6f' % a for a in x) # XXX fStr self.test('parse3llh', x, '51.477811, -0.001475, 0.000000') for a, x in (((), '''45°45'45.36"'''), ((F_D, None), '45.7626°'), ((F_DM, None), "45°45.756'"), ((F_DMS, None), '''45°45'45.36"'''), ((F_DEG, None), '45.7626'), ((F_MIN, None), '4545.756'), ((F_SEC, None), '454545.36'), ((F_RAD, None), '0.79871'), ((F_D, 6), '45.7626°'), ((F_DM, -4), "45°45.7560'"), ((F_DMS, 2), '''45°45'45.36"'''), ((F_DEG, -6), '45.762600'), ((F_MIN, -5), '4545.75600'), ((F_SEC, -3), '454545.360'), ((F_RAD, -6), '0.798708')): t = 'toDMS(%s)' % (a[:1] or '') self.test(t, toDMS(45.76260, *a), x) self.test('compassAngle0', compassAngle(0, 0, 0, 0), 0.0) self.test('compassAngle1', compassAngle(0, 0, 10, 10), 45.0) self.test('compassAngle2', compassAngle(0, 0, 0, 10), 90.0) self.test('compassAngle3', compassAngle(0, 0, -1, 0), 180.0) self.test('compassAngle4', compassAngle(0, 0, -1, -1), 225.0) self.test('compassAngle5', compassAngle(0, 0, 0, -1), 270.0) self.test('compassAngle6', compassAngle(0, 0, 1, -1), 315.0) self.test('compassAngle7', compassAngle(0, 0, 99, -1), 359.4, fmt='%.1f') for a, x in (((1, ), 'N'), ((0, ), 'N'), ((-1, ), 'N'), ((359, ), 'N'), ((24, ), 'NNE'), ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((226, ), 'SW'), ((226, 1), 'W'), ((226, 2), 'SW'), ((226, 3), 'SW'), ((237, ), 'WSW'), ((237, 1), 'W'), ((237, 2), 'SW'), ((237, 3), 'WSW')): self.test('compassPoint', compassPoint(*a), x) for a, x in enumerate( ('NbE', 'NEbN', 'NEbE', 'EbN', 'EbS', 'SEbE', 'SEbS', 'SbE', 'SbW', 'SWbS', 'SWbW', 'WbS', 'WbN', 'NWbW', 'NWbN', 'NbW')): a = 11.25 + a * 22.5 self.test('compassPoint', compassPoint(a, 4), x) a = compassAngle(_LHR.lat, _LHR.lon, _FRA.lat, _FRA.lon) self.test('compassAngle', a, 100.016848, fmt='%.6f') b = _LHR.initialBearingTo(_FRA) self.test('initialBearingTo', b, 102.432182, fmt='%.6f') d = equirectangular(_LHR.lat, _LHR.lon, _FRA.lat, _FRA.lon) self.test('equirectangular', m2km(d), 592.185, fmt='%.3f') d = _LHR.distanceTo(_FRA) self.test('distanceTo', m2km(d), 591.831, fmt='%.3f')
def testDms(self): # dms module tests self.test('parseDMS', parseDMS( '0.0°'), '0.0') self.test('parseDMS', parseDMS( '0°'), '0.0') self.test('parseDMS', parseDMS('''000°00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00'00.0"'''), '0.0') self.test('parseDMS', parseDMS('''000° 00'00"'''), '0.0') self.test('parseDMS', parseDMS('''000°00 ' 00.0"'''), '0.0') r = rangerrors(True) try: self.test('parseDMS', parseDMS(181, clip=180), 'ValueError') except ValueError as x: self.test('parseDMS', str(x), '181.0 beyond 180 degrees') rangerrors(False) try: self.test('parseDMS', parseDMS(-91, clip=90), -90) except ValueError as x: self.test('parseDMS', str(x), '-90') rangerrors(r) x = parse3llh('000° 00′ 05.31″W, 51° 28′ 40.12″ N') x = ', '.join('%.6f' % a for a in x) # XXX fStr self.test('parse3llh', x, '51.477811, -0.001475, 0.000000') for a, x in (((), '''45°45'45.36"'''), ((F_D, None), '45.7626°'), ((F_DM, None), "45°45.756'"), ((F_DMS, None), '''45°45'45.36"'''), ((F_DEG, None), '45.7626'), ((F_MIN, None), '4545.756'), ((F_SEC, None), '454545.36'), ((F_RAD, None), '0.79871'), ((F_D, 6), '45.7626°'), ((F_DM, -4), "45°45.7560'"), ((F_DMS, 2), '''45°45'45.36"'''), ((F_DEG, -6), '45.762600'), ((F_MIN, -5), '4545.75600'), ((F_SEC, -3), '454545.360'), ((F_RAD, -6), '0.798708')): t = 'toDMS(%s)' % (a[:1] or '') self.test(t, toDMS(45.76260, *a), x) # <https://GitHub.com/chrisveness/geodesy/blob/master/test/dms-tests.js> for a, x in (((1,), 'N'), ((0,), 'N'), ((-1,), 'N'), ((359,), 'N'), ((24,), 'NNE'), ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((226,), 'SW'), ((226, 1), 'W'), ((226, 2), 'SW'), ((226, 3), 'SW'), ((237,), 'WSW'), ((237, 1), 'W'), ((237, 2), 'SW'), ((237, 3), 'WSW'), # Martin Schultz ((11.25,), 'NNE'), ((11.249,), 'N'), ((-11.25,), 'N'), ((348.749,), 'NNW'), ((45, 1), 'E'), ((44.99, 1), 'N'), ((45, 2), 'NE'), ((44.99, 2), 'NE'), ((45, 3), 'NE'), ((44.99, 3), 'NE'), ((45, 4), 'NE'), ((44.99, 4), 'NE'), # XXX ((22.5, 1), 'N'), ((22.49, 1), 'N'), ((22.5, 2), 'NE'), ((22.49, 2), 'N'), ((22.5, 3), 'NNE'), ((22.49, 3), 'NNE'), ((22.5, 4), 'NNE'), ((22.49, 4), 'NNE'), # XXX ((11.25, 1), 'N'), ((11.249, 1), 'N'), ((11.25, 2), 'N'), ((11.249, 2), 'N'), ((11.25, 3), 'NNE'), ((11.249, 3), 'N'), ((11.25, 4), 'NbE'), ((11.249, 4), 'NbE'), # XXX # examples ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((24,), 'NNE'), ((18, 3), 'NNE'), ((11, 4), 'NbE'), ((30, 4), 'NEbN')): self.test('compassPoint%s' % (a,), compassPoint(*a), x) for a, x in enumerate(('NbE', 'NEbN', 'NEbE', 'EbN', 'EbS', 'SEbE', 'SEbS', 'SbE', 'SbW', 'SWbS', 'SWbW', 'WbS', 'WbN', 'NWbW', 'NWbN', 'NbW')): a = 11.25 + a * 22.5 self.test('compassPoint(%s)' % (a,), compassPoint(a, 4), x)
def testEcef(self, Ecef): self.test(Ecef.__name__, '...', '...', nl=1) Karney = Ecef is EcefKarney g = Ecef(Datums.WGS84, name='Test') self.test('name', g.name, 'Test') t = g.toStr2() self.test('toStr', t, g.classname, known=True) # <https://GeographicLib.SourceForge.io/html/classGeographicLib_1_1Geocentric.html> t = g.forward(27.99, 86.93, 8820) # Mt Everest self.test('forward', fStr(t[3:6], prec=2), '27.99, 86.93, 8820.0') self.test('forward', fStr(t[0:3], prec=1), '302271.4, 5635928.4, 2979666.1') self.test('name', t.name, 'Test') t = g.reverse(302271.4, 5635928.4, 2979666.1) self.test('reverse', fStr(t[0:3], prec=1), '302271.4, 5635928.4, 2979666.1') self.test('reverse', fStr(t[3:6], prec=2), '27.99, 86.93, 8820.01') self.test('case', t.C, 2 if Karney else 1) self.test('name', t.name, 'Test') # <https://GeographicLib.SourceForge.io/html/classGeographicLib_1_1Geocentric.html> t = g.reverse(302e3, 5636e3, 2980e3) self.test('reverse', fStr(t[0:3], prec=1), '302000.0, 5636000.0, 2980000.0') self.test('reverse', fStr(t[3:6], prec=2), '27.99, 86.93, 9027.03') # 8820.0 self.test('case', t.C, 2 if Karney else 1) t = g.forward(27.99, 86.93, 8820.0) self.test('forward', fStr(t[3:6], prec=2), '27.99, 86.93, 8820.0') self.test('forward', fStr(t[0:3], prec=2), '302271.43, 5635928.37, 2979666.13') # <https://GeographicLib.SourceForge.io/html/CartConvert.1.html> t = g.forward(33.3, 44.4, 6000) self.test('forward', fStr(t[3:6], prec=2), '33.3, 44.4, 6000.0') self.test('forward', fStr(t[0:3], prec=2), '3816209.6, 3737108.55, 3485109.57') t = g.reverse(3816209.6, 3737108.55, 3485109.57) self.test('reverse', fStr(t[0:3], prec=2), '3816209.6, 3737108.55, 3485109.57') self.test('reverse', fStr(t[3:6], prec=3), '33.3, 44.4, 5999.996') # 6000.0 self.test('case', t.C, 2 if Karney else 1) # <https://GeographicLib.SourceForge.io/html/CartConvert.1.html> t = g.reverse(30000, 30000, 0) self.test('reverse', fStr(t[0:3], prec=1), '30000.0, 30000.0, 0.0') self.test('reverse', fStr(t[3:6], prec=3), '6.483, 45.0, -6335709.726', known=not Karney) self.test('case', t.C, 3 if Karney else 1) t = g.forward(6.483, 45.0, -6335709.726) self.test('forward', fStr(t[3:6], prec=3), '6.483, 45.0, -6335709.726') self.test('forward', fStr(t[0:3], prec=1), '30000.0, 30000.0, -0.0', known=True) # Rey-Jer You <https://www.ResearchGate.net/publication/240359424> for i, (x, y, z, h) in enumerate( ((-2259148.993, 3912960.837, 4488055.516, 1e3), (-2259502.546, 3913573.210, 4488762.622, 2e3), (-2259856.100, 3914185.582, 4489469.729, 3e3), (-2260209.653, 3914797.955, 4490176.836, 4e3), (-2262330.973, 3918472.189, 4494419.477, 1e4), (-2265866.507, 3924595.914, 4501490.544, 2e4), (-2294150.778, 3973585.709, 4558059.087, 1e5), (-2541638.152, 4402246.414, 5053033.834, 8e5), (-2612348.830, 4524720.901, 5194455.190, 1e6))): i = '-%d' % (i + 1, ) r = '45.0, 120.0, %.1f' % (h, ) # Zero and First order columns t = g.reverse(x, y, z) self.test('reverse' + i, fStr(t[3:6], prec=3), r) t = g.forward(t.lat, t.lon, t.height) self.test('forward' + i, fStr(t[0:3], prec=1), fStr((x, y, z), prec=1)) # <https://www.OrdnanceSurvey.co.UK/documents/resources/guide-coordinate-systems-great-britain.pdf> pp 47 g = Ecef(Ellipsoids.GRS80, name='OS-UK') self.test('name', g.name, 'OS-UK') t = g.forward( parse3llh('''53°36′43.1653"N, 001°39′51.9920"W, 299.800''')) self.test('forward', fStr(t[3:6], prec=8), '53.61199036, -1.66444222, 299.8') self.test('forward', fStr(t[0:3], prec=2), '3790644.9, -110149.21, 5111482.97') t = g.reverse(3790644.9, -110149.21, 5111482.97) self.test('reverse', fStr(t[0:3], prec=2), '3790644.9, -110149.21, 5111482.97') self.test('reverse', fStr(t[3:5], prec=8), '53.61199036, -1.66444223') self.test('reverse.lat', latDMS(t.lat, prec=4), '53°36′43.1653″N') self.test('reverse.lon', lonDMS(t.lon, prec=4), '001°39′51.992″W') self.test('reverse.height', fStr(t.height, prec=-3), '299.800') self.test('case', t.C, 2 if Karney else 1)
def testDms(self): # dms module tests for i, t in enumerate( ('0.0°', '0°', '''000°00'00"''', '''000°00'00.0"''', '''000° 00'00"''', '''000°00 ' 00.0"''', '''000° 00' 00.0''')): self.test('parseDMS' + str(i + 1), parseDMS(t), '0.0') self.test('parseDMS' + str(i + 2), parseDMS('000°-00′-00.0"', sep='-'), '0.0') r = rangerrors(True) try: self.test('parseDMS', parseDMS(181, clip=180), 'ValueError') except ValueError as x: self.test('parseDMS', str(x), '181.0 beyond 180 degrees') rangerrors(False) try: self.test('parseDMS', parseDMS(-91, clip=90), -90) except ValueError as x: self.test('parseDMS', str(x), '-90') rangerrors(r) x = parse3llh('000° 00′ 05.31″W, 51° 28′ 40.12″ N') self.test('parse3llh', fStr(x, prec=6), '51.477811, -0.001475, 0.0') for a, x in (((), '''45°45'45.36"'''), ((F_D, None), '45.7626°'), ((F_DM, None), "45°45.756'"), ((F_DMS, None), '''45°45'45.36"'''), ((F_DEG, None), '45.7626'), ((F_MIN, None), '4545.756'), ((F_SEC, None), '454545.36'), ((F_RAD, None), '0.79871'), ((F_D, 6), '45.7626°'), ((F_DM, -4), "45°45.7560'"), ((F_DMS, 2), '''45°45'45.36"'''), ((F_DEG, -6), '45.762600'), ((F_MIN, -5), '4545.75600'), ((F_SEC, -3), '454545.360'), ((F_RAD, -6), '0.798708')): t = 'toDMS(%s)' % (a[:1] or '') self.test(t, toDMS(45.76260, *a), x) for a, x in (((F_D_, 6), '45.7626°'), ((F_DM_, -4), "45°45.7560'"), ((F_DMS_, 2), '''45°45'45.36"'''), ((F_DEG_, -6), '45.762600'), ((F_MIN_, -5), '4545.75600'), ((F_SEC_, -3), '454545.360'), ((F_RAD_, -6), '0.798708')): t = 'toDMS(%s)' % (a[:1] or '') self.test(t, toDMS(45.76260, *a), x) self.test(t, toDMS(-45.76260, *a), '-' + x) for a, x in (((F_D__, 6), '45.7626°'), ((F_DM__, -4), "45°45.7560'"), ((F_DMS__, 2), '''45°45'45.36"'''), ((F_DEG__, -6), '45.762600'), ((F_MIN__, -5), '4545.75600'), ((F_SEC__, -3), '454545.360'), ((F_RAD__, -6), '0.798708')): t = 'toDMS(%s)' % (a[:1] or '') self.test(t, toDMS(45.76260, *a), '+' + x) self.test(t, toDMS(-45.76260, *a), '-' + x) # <https://GitHub.com/chrisveness/geodesy/blob/master/test/dms-tests.js> for a, x in ( ((1, ), 'N'), ((0, ), 'N'), ((-1, ), 'N'), ((359, ), 'N'), ((24, ), 'NNE'), ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((226, ), 'SW'), ((226, 1), 'W'), ((226, 2), 'SW'), ((226, 3), 'SW'), ((237, ), 'WSW'), ((237, 1), 'W'), ((237, 2), 'SW'), ((237, 3), 'WSW'), # Martin Schultz ((11.25, ), 'NNE'), ((11.249, ), 'N'), ((-11.25, ), 'N'), ((348.749, ), 'NNW'), ((45, 1), 'E'), ((44.99, 1), 'N'), ((45, 2), 'NE'), ((44.99, 2), 'NE'), ((45, 3), 'NE'), ((44.99, 3), 'NE'), ((45, 4), 'NE'), ((44.99, 4), 'NE'), # XXX ((22.5, 1), 'N'), ((22.49, 1), 'N'), ((22.5, 2), 'NE'), ((22.49, 2), 'N'), ((22.5, 3), 'NNE'), ((22.49, 3), 'NNE'), ((22.5, 4), 'NNE'), ((22.49, 4), 'NNE'), # XXX ((11.25, 1), 'N'), ((11.249, 1), 'N'), ((11.25, 2), 'N'), ((11.249, 2), 'N'), ((11.25, 3), 'NNE'), ((11.249, 3), 'N'), ((11.25, 4), 'NbE'), ((11.249, 4), 'NbE'), # XXX # examples ((24, 1), 'N'), ((24, 2), 'NE'), ((24, 3), 'NNE'), ((24, ), 'NNE'), ((18, 3), 'NNE'), ((11, 4), 'NbE'), ((30, 4), 'NEbN')): self.test('compassPoint%s' % (a, ), compassPoint(*a), x) for a, x in enumerate( ('NbE', 'NEbN', 'NEbE', 'EbN', 'EbS', 'SEbE', 'SEbS', 'SbE', 'SbW', 'SWbS', 'SWbW', 'WbS', 'WbN', 'NWbW', 'NWbN', 'NbW')): a = 11.25 + a * 22.5 self.test('compassPoint(%s)' % (a, ), compassPoint(a, 4), x) t = degDMS(1.0101, prec=5, s_D='', pos='+') # coverage self.test('_DEG', t, '+1.0101') t = degDMS(0.0101, prec=5, s_S='', pos='+') self.test('_MIN', t, "+0.606'") t = degDMS(0.0101, prec=5, s_M='', pos='+') self.test('_SEC', t, '+36.36"')