def testClipy(self, module): # MCCABE 15
self.subtitle(module)
LatLon = module.LatLon
# <https://www.CS.Helsinki.Fi/group/goa/viewing/leikkaus/example1.html>
ll, ur = LatLon(0, 0), LatLon(10, 10)
ps = LatLon(3, -5), LatLon(5, 5), LatLon(5, 5), LatLon(9, 15)
t = tuple(clipCS4(ps, ll, ur, closed=True, inull=True))
p1, p2, i, j = t[0] # closing edge
self.test('clipCS4.p1', p1, '07.5°N, 010.0°E')
self.test('clipCS4.p2', p2, '04.5°N, 000.0°E')
self.test('clipCS4.i', i, 3)
self.test('clipCS4.j', j, 0)
p1, p2, i, j = t[2] # only null edge
self.test('clipCS4.p1', p1, '05.0°N, 005.0°E')
self.test('clipCS4.p2', p2, '05.0°N, 005.0°E')
self.test('clipCS4.i', i, 1)
self.test('clipCS4.j', j, 2)
t = tuple(clipLB6(ps, ll, ur, closed=True, inull=True))
p1, p2, i, fi, fj, j = t[0] # closing edge
self.test('clipLB6.p1', p1, '07.5°N, 010.0°E')
self.test('clipLB6.p2', p2, '04.5°N, 000.0°E')
self.test('clipLB6.i', i, 3)
self.test('clipLB6.fi', fi, 3.25, prec=2)
self.test('clipLB6.fi', fi.fractional(ps, LatLon=LatLon), p1)
self.test('clipLB6.fj', fj, 3.75, prec=2)
self.test('clipLB6.fj', fj.fractional(ps, LatLon=LatLon), p2)
self.test('clipLB6.j', j, 0)
self.test('clipLB6.fin', fi.fin, 4)
p1, p2, i, fi, fj, j = t[2] # only null edge
self.test('clipLB6.p1', p1, '05.0°N, 005.0°E')
self.test('clipLB6.p2', p2, '05.0°N, 005.0°E')
self.test('clipLB6.i', i, 1)
self.test('clipLB6.fi', fi, 1.00, prec=2)
self.test('clipLB6.fi', fi.fractional(ps, LatLon=LatLon), p1)
self.test('clipLB6.fj', fj, 2.00, prec=2)
self.test('clipLB6.fj', fj.fractional(ps, LatLon=LatLon), p2)
self.test('clipLB6.j', j, 2)
self.test('clipLB6.fin', fi.fin, 4)
ll, ur = LatLon(60, 70), LatLon(70, 130)
ps = LatLon(20, 30), LatLon(80, 170)
for p1, p2, i, j in clipCS4(ps, ll, ur):
self.test('clipCS4.p1', p1, '60.0°N, 123.333333°E')
self.test('clipCS4.p2', p2, '62.857143°N, 130.0°E')
self.test('clipCS4.i', i, 0)
self.test('clipCS4.j', j, 1)
for p1, p2, i, fi, fj, j in clipLB6(ps, ll, ur):
self.test('clipLB6.p1', p1, '60.0°N, 123.333333°E')
self.test('clipLB6.p2', p2, '62.857143°N, 130.0°E')
self.test('clipLB6.i', i, 0)
self.test('clipLB6.fi', fi, 0.666667, prec=6)
self.test('clipLB6.fi', fi.fractional(ps, LatLon=LatLon), p1)
self.test('clipLB6.fj', fj, 0.714286, prec=6)
self.test('clipLB6.fj', fj.fractional(ps, LatLon=LatLon), p2)
self.test('clipLB6.j', j, 1)
self.test('clipLB6.fin', fi.fin, 0)
ll, ur = LatLon(15, 15), LatLon(20, 20)
ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
for p1, p2, i, j in clipCS4(ps, ll, ur, closed=True, inull=False):
self.test('clipCS4.p1', p1, '17.5°N, 020.0°E')
self.test('clipCS4.p2', p2, '16.25°N, 015.0°E')
self.test('clipCS4.i', i, 2) # closing edge
self.test('clipCS4.j', j, 0)
for p1, p2, i, fi, fj, j in clipLB6(ps,
ll,
ur,
closed=True,
inull=False):
self.test('clipLB6.p1', p1, '17.5°N, 020.0°E')
self.test('clipLB6.p2', p2, '16.25°N, 015.0°E')
self.test('clipLB6.i', i, 2) # closing edge
self.test('clipLB6.fi', fi, 2.500, prec=3) # closing edge
self.test('clipLB6.fi', fi.fractional(ps, LatLon=LatLon), p1)
self.test('clipLB6.fj', fj, 2.750, prec=3) # closing edge
self.test('clipLB6.fj', fj.fractional(ps, LatLon=LatLon), p2)
self.test('clipLB6.j', j, 0) # closing edge
self.test('clipLB6.fin', fi.fin, 3)
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
sh = tuple(clipSH(ps, (ll, ur)))
self.testClipSH_('clipSH1.', sh, [
LatLon(20, 20),
LatLon(17.5, 20),
LatLon(16.25, 15),
LatLon(20, 15)
])
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(30, 10), LatLon(30, 30), LatLon(10, 20)
for r in ('', 'reversed.'):
sh = tuple(clipSH(ps, cs, closed=False))
self.testClipSH_('clipSH2.' + r, sh, [
LatLon(18.571, 24.286),
LatLon(16.667, 16.667),
LatLon(20, 15),
LatLon(25, 20),
LatLon(22, 26)
])
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
# cs = LatLon(30, 10), LatLon(30, 30), LatLon(10, 20)
sh = clipSH3(ps, cs, closed=True)
for r3, x3 in zip(
sh,
_lles3((LatLon(16.667, 16.667), True),
(LatLon(20, 15), False), (LatLon(25, 20), True),
(LatLon(22, 26), True),
(LatLon(18.571, 24.286), False))):
self.testClipSH_('clipSH3.' + r, r3[:2], x3[:2])
self.test('clipSH3.edge.' + r, r3[2], x3[2])
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(30, -10), LatLon(30, -30), LatLon(10, -20) # all outside
for r in ('', '.reversed'):
sh = tuple(clipSH(ps, cs, closed=True))
self.test('clipSH.allout' + r, sh, ())
sh = tuple(clipSH3(ps, cs, closed=True))
self.test('clipSH3.allout' + r, sh, ())
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(15, 10), LatLon(30,
10), LatLon(30,
30), LatLon(15,
30) # all inside
for r in ('', 'reversed.'):
sh = tuple(clipSH(ps, cs, closed=True))
self.testClipSH_('clipSH.allin.' + r, sh,
(ps[-1], ) + ps) # close ps
sh = tuple(clipSH3(ps, cs, closed=False))
for r3, x3 in zip(sh, _llsTrue3(ps)):
self.testClipSH_('clipSH3.allin.' + r, r3[:2], x3[:2])
self.test('clipSH3.edge.' + r, r3[2], x3[2])
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(10, 10), LatLon(20, 20), LatLon(10,
20), LatLon(20,
10) # warped
for r in ('', 'reversed.'):
try:
# use list to force exception, see
# <https://RickardLindberg.me/writing/bitten-by-python-generators>
sh = list(clipSH(ps, cs))
t = ClipError.__name__
except ClipError as x:
t = sh = str(x) # .split(':')[0]
self.test('clipSH.warped' + r, sh, t)
try:
# use list to force exception, see
# <https://RickardLindberg.me/writing/bitten-by-python-generators>
sh = list(clipSH3(ps, cs))
t = ClipError.__name__
except ClipError as x:
t = sh = str(x) # .split(':')[0]
self.test('clipSH3.warped' + r, sh, t)
cs = tuple(reversed(cs))
def testClipy(self, LatLon): # MCCABE 13
ll, ur = LatLon(60, 70), LatLon(70, 130)
ps = LatLon(20, 30), LatLon(80, 170)
for p1, p2, i in clipCS3(ps, ll, ur):
self.test('clipCS3.p1', p1, '60.0°N, 123.333333°E')
self.test('clipCS3.p2', p2, '62.857143°N, 130.0°E')
self.test('clipCS3.i', i, 1)
ll, ur = LatLon(15, 15), LatLon(20, 20)
ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
for p1, p2, i in clipCS3(ps, ll, ur, closed=True, inull=False):
self.test('clipCS3.p1', p1, '17.5°N, 020.0°E')
self.test('clipCS3.p2', p2, '16.25°N, 015.0°E')
self.test('clipCS3.i', i, 0) # closing edge
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
sh = tuple(clipSH(ps, (ll, ur)))
self.testClipSH_('clipSH1.', sh, [LatLon(20, 20), LatLon(17.5, 20), LatLon(16.25, 15), LatLon(20, 15)])
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(30, 10), LatLon(30, 30), LatLon(10, 20)
for r in ('', 'reversed.'):
sh = tuple(clipSH(ps, cs, closed=False))
self.testClipSH_('clipSH2.' + r, sh, [LatLon(18.571, 24.286), LatLon(16.667, 16.667), LatLon(20, 15), LatLon(25, 20), LatLon(22, 26)])
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
# cs = LatLon(30, 10), LatLon(30, 30), LatLon(10, 20)
sh = clipSH3(ps, cs, closed=True)
for r3, x3 in zip(sh, _lles3((LatLon(16.667, 16.667), True),
(LatLon(20, 15), False),
(LatLon(25, 20), True),
(LatLon(22, 26), True),
(LatLon(18.571, 24.286), False))):
self.testClipSH_('clipSH3.' + r, r3[:2], x3[:2])
self.test('clipSH3.edge.' + r, r3[2], x3[2])
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(30, -10), LatLon(30, -30), LatLon(10, -20) # all outside
for r in ('', '.reversed'):
sh = tuple(clipSH(ps, cs, closed=True))
self.test('clipSH.allout' + r, sh, ())
sh = tuple(clipSH3(ps, cs, closed=True))
self.test('clipSH3.allout' + r, sh, ())
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(15, 10), LatLon(30, 10), LatLon(30, 30), LatLon(15, 30) # all inside
for r in ('', 'reversed.'):
sh = tuple(clipSH(ps, cs, closed=True))
self.testClipSH_('clipSH.allin.' + r, sh, (ps[-1],) + ps) # close ps
sh = tuple(clipSH3(ps, cs, closed=False))
for r3, x3 in zip(sh, _llsTrue3(ps)):
self.testClipSH_('clipSH3.allin.' + r, r3[:2], x3[:2])
self.test('clipSH3.edge.' + r, r3[2], x3[2])
cs = tuple(reversed(cs))
# ps = LatLon(15, 10), LatLon(25, 20), LatLon(20, 30)
cs = LatLon(10, 10), LatLon(20, 20), LatLon(10, 20), LatLon(20, 10) # warped
for r in ('', 'reversed.'):
try:
# use list to force exception, see
# <https://RickardLindberg.me/writing/bitten-by-python-generators>
sh = list(clipSH(ps, cs))
t = ClipError.__name__
except ClipError as x:
t = sh = str(x) # .split(':')[0]
self.test('clipSH.warped' + r, sh, t)
try:
# use list to force exception, see
# <https://RickardLindberg.me/writing/bitten-by-python-generators>
sh = list(clipSH3(ps, cs))
t = ClipError.__name__
except ClipError as x:
t = sh = str(x) # .split(':')[0]
self.test('clipSH3.warped' + r, sh, t)
cs = tuple(reversed(cs))