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)