def testUnregister(self):
self.subtitle(named, 'ing %s ' % ('unregister',))
from pygeodesy import Conics, Ellipsoids, RefFrames, Transforms
self.test('Conics', len(Conics), 8)
for n, c in tuple(Conics.items()):
c.unregister() # coverage _NamedEnum.unregister
self.test('Conics.' + n + '.unregister', getattr(Conics, n, None), None)
self.test('Conics', len(Conics), 0)
self.test('Datums', len(Datums), 18)
for n, d in tuple(Datums.items()):
Datums.unregister(d) # coverage _NamedEnum.unregister
self.test('Datums.unregister(%s)' % (n,), getattr(Datums, n, None), None)
self.test('Datums', len(Datums), 0)
self.test('Ellipsoids', len(Ellipsoids), 41)
for n, e in tuple(Ellipsoids.items()):
e.unregister() # coverage _NamedEnum.unregister
self.test('Ellipsoids.' + n + '.unregister', getattr(Ellipsoids, n, None), None)
self.test('Ellipsoids', len(Ellipsoids), 0)
self.test('RefFrames', len(RefFrames), 12)
for n, r in tuple(RefFrames.items()):
r.unregister() # coverage _NamedEnum.unregister
self.test('RefFrames.' + n + '.unregister', getattr(RefFrames, n, None), None)
self.test('RefFrames', len(RefFrames), 0)
self.test('Transforms', len(Transforms), 18)
for n, x in tuple(Transforms.items()):
x.unregister() # coverage _NamedEnumItem.unregister
self.test('Transforms.' + n + '.unregister', getattr(Transforms, n, None), None)
self.test('Transforms', len(Transforms), 0)
def testLcc(self, module):
self.subtitle(module, 'Lcc')
LatLon = module.LatLon
lb = Lcc(448251, 5411932.0001, name='lb1')
self.test('lb1', lb.toStr(4), '448251.0 5411932.0001')
self.test('lb1', lb.toStr(sep=', '), '448251, 5411932')
self.test('lb1', lb.toStr2(), '[E:448251, N:5411932]')
self.test('lb1', lb.conic.name2, 'WRF_Lb.WGS84')
self.test('lb1', lb.name, 'lb1')
self.test('lb1', lb.latlon, '(81.929348, -79.558697)')
self.test('lb1', lb.philam, '(1.429937, -1.388561)')
ll = LatLon(46.5, 3)
self.test('LatLon', ll, '46.5°N, 003.0°E')
self.test('LatLon', ll.toStr(form=F_DMS),
'46°30′00.0″N, 003°00′00.0″E')
lb = toLcc(ll, conic=Conics.Fr93Lb)
self.test('toLcc1', str(lb), '700000 6600000')
self.test('toLcc1', lb.toLatLon(LatLon), '46.5°N, 003.0°E')
lb = Lcc(1894410.9, 1564649.5, conic=Conics.Snyder, name='lb2')
self.test('lb2', lb, '1894411 1564650')
self.test('lb2', lb.conic.datum.ellipsoid.name, 'Clarke1866')
self.test('lb2', lb.name, 'lb2')
ll = lb.toLatLon(LatLon) # Clark1866
self.test('toLatLon2', ll.toStr(prec=6, form=F_D), '35.0°N, 075.0°W')
self.test('toLatLon2', ll.toStr(prec=4, form=F_DMS),
'35°00′00.0007″N, 074°59′59.9997″W')
self.test('toLatLon2', ll.datum.name, 'NAD27')
lb = toLcc(ll, conic=Conics.Snyder)
self.test('toLcc2', lb.toStr(prec=1), '1894410.9 1564649.5')
self.test('toLcc2', lb.toStr2(), '[E:1894411, N:1564649]')
self.test('toLcc2', lb.conic.name2, 'Snyder.NAD27')
for n, c in sorted(Conics.items()):
d = abs(c.par1 - c.par2)
if d > 0: # test corners of the conic
for ll in (LatLon(c.par1, c.lon0 - d, datum=c.datum),
LatLon(c.par1, c.lon0, datum=c.datum),
LatLon(c.par1, c.lon0 + d, datum=c.datum),
LatLon(c.par2, c.lon0 - d, datum=c.datum),
LatLon(c.par2, c.lon0, datum=c.datum),
LatLon(c.par2, c.lon0 + d, datum=c.datum)):
# self.test(n, ll, str(ll)) # PYCHOK expected
lb = toLcc(ll, conic=c)
# self.test(n, lb, '')
ll_ = lb.toLatLon(LatLon)
self.test(n, ll, str(ll_))
self.test(n, ll.datum.name, ll_.datum.name)
def testLcc(self, LatLon):
lb = Lcc(448251, 5411932.0001)
self.test('lb1', lb.toStr(4), '448251.0 5411932.0001')
self.test('lb1', lb.toStr(sep=', '), '448251, 5411932')
self.test('lb1', lb.conic.name2, 'WRF_Lb.WGS84')
ll = LatLon(46.5, 3)
self.test('LatLon', ll, '46.5°N, 003.0°E')
self.test('LatLon', ll.toStr(form=F_DMS),
'46°30′00.0″N, 003°00′00.0″E')
lb = toLcc(ll, conic=Conics.Fr93Lb)
self.test('toLcc1', str(lb), '700000 6600000')
self.test('toLcc1', lb.toLatLon(LatLon), '46.5°N, 003.0°E')
lb = Lcc(1894410.9, 1564649.5, conic=Snyder)
self.test('lb2', lb, '1894411 1564650')
self.test('lb2', lb.conic.datum.ellipsoid.name, 'Clarke1866')
ll = lb.toLatLon(LatLon) # Clark1866
self.test('toLatLon2', ll.toStr(prec=6, form=F_D), '35.0°N, 075.0°W')
self.test('toLatLon2', ll.toStr(prec=4, form=F_DMS),
'35°00′00.0007″N, 074°59′59.9997″W')
self.test('toLatLon2', ll.datum.name, 'NAD27')
lb = toLcc(ll, conic=Snyder)
self.test('toLcc2', lb.toStr(prec=1), '1894410.9 1564649.5')
self.test('toLcc2', lb.conic.name2, 'Snyder.NAD27')
for n, c in sorted(Conics.items()):
d = abs(c.par1 - c.par2)
if d > 0: # test corners of the conic
for ll in (LatLon(c.par1, c.lon0 - d, datum=c.datum),
LatLon(c.par1, c.lon0, datum=c.datum),
LatLon(c.par1, c.lon0 + d, datum=c.datum),
LatLon(c.par2, c.lon0 - d, datum=c.datum),
LatLon(c.par2, c.lon0, datum=c.datum),
LatLon(c.par2, c.lon0 + d, datum=c.datum)):
# self.test(n, ll, str(ll)) # PYCHOK expected
lb = toLcc(ll, conic=c)
# self.test(n, lb, '')
ll_ = lb.toLatLon(LatLon)
self.test(n, ll, str(ll_))
self.test(n, ll.datum.name, ll_.datum.name)
ll_ = lb.toLatLon(LatLon)
self.test(n, ll, str(ll_))
self.test(n, ll.datum.name, ll_.datum.name)
if __name__ == '__main__':
from pygeodesy import ellipsoidalNvector, ellipsoidalVincenty, lcc
# Snyder, pp 297 <http://pubs.er.USGS.gov/djvu/PP/PP_1395.pdf>
Conic(ellipsoidalVincenty.LatLon(23, -96, datum=Datums.NAD27),
33,
45,
E0=0,
N0=0,
name='Snyder')
t = Tests(__file__, __version__, lcc)
t.testLcc(ellipsoidalNvector)
t.testLcc(ellipsoidalVincenty)
t.testConic(ellipsoidalNvector, 'N')
t.testConic(ellipsoidalVincenty, 'V')
for n in ('', 'N', 'V'):
Conics.unregister('Snyder' + n)
t.results()
t.exit()