def test_intersection(self):
arc1 = Arc(SCoordinate(0, 0),
SCoordinate(np.deg2rad(10), np.deg2rad(10)))
arc2 = Arc(SCoordinate(0, np.deg2rad(10)),
SCoordinate(np.deg2rad(10), 0))
lon, lat = arc1.intersection(arc2)
self.assertTrue(np.allclose(np.rad2deg(lon), 5))
self.assertEqual(
np.rad2deg(lat).round(7), round(5.0575148968282093, 7))
arc1 = Arc(SCoordinate(0, 0),
SCoordinate(np.deg2rad(10), np.deg2rad(10)))
self.assertTrue(arc1.intersection(arc1) is None)
arc1 = Arc(
SCoordinate(np.deg2rad(24.341215776575297),
np.deg2rad(44.987819588259327)),
SCoordinate(np.deg2rad(18.842727517611817),
np.deg2rad(46.512483610284178)))
arc2 = Arc(
SCoordinate(np.deg2rad(20.165961750361905),
np.deg2rad(46.177305385810541)),
SCoordinate(np.deg2rad(20.253297585831707),
np.deg2rad(50.935830837274324)))
inter = SCoordinate(np.deg2rad(20.165957021925202),
np.deg2rad(46.177022633103398))
self.assertEqual(arc1.intersection(arc2), inter)
arc1 = Arc(
SCoordinate(np.deg2rad(-2.4982818108326734),
np.deg2rad(48.596644847869655)),
SCoordinate(np.deg2rad(-2.9571441235622835),
np.deg2rad(49.165688435261394)))
arc2 = Arc(
SCoordinate(np.deg2rad(-3.4976667413531688),
np.deg2rad(48.562704872921373)),
SCoordinate(np.deg2rad(-5.893976312685715),
np.deg2rad(48.445795283217116)))
self.assertTrue(arc1.intersection(arc2) is None)
def test_intersecting_arcs_crossing_antimeridian(self):
import copy
arc1 = Arc(SCoordinate(*np.deg2rad((-180.0, -90.0))),
SCoordinate(*np.deg2rad((-180.0, 90.0))))
arc2 = Arc(SCoordinate(*np.deg2rad((-171.03, -76.75))),
SCoordinate(*np.deg2rad((170.34, -71.36))))
arc1_orig = copy.deepcopy(arc1)
arc2_orig = copy.deepcopy(arc2)
point = arc1.intersection(arc2)
# Assert original arcs are unaffected
assert arc1_orig.end.lon == arc1.end.lon
assert arc2_orig.end.lon == arc2.end.lon
# Assert intersection result
assert point == SCoordinate(*np.deg2rad((-180.0, -74.7884716)))
def test_disjoint_arcs_crossing_antimeridian(self):
import copy
arc1 = Arc(SCoordinate(*np.deg2rad((143.76, 0))),
SCoordinate(*np.deg2rad((143.95, 7.33))))
arc2 = Arc(SCoordinate(*np.deg2rad((170.34, 71.36))),
SCoordinate(*np.deg2rad((-171.03, 76.75))))
arc1_orig = copy.deepcopy(arc1)
arc2_orig = copy.deepcopy(arc2)
point = arc1.intersection(arc2)
# Assert original arcs are unaffected
assert arc1_orig.end.lon == arc1.end.lon
assert arc2_orig.end.lon == arc2.end.lon
# Assert disjoint arcs returns None
assert isinstance(point, type(None))
