def get_poly_position(self):
position = plane.Position()
if self.caom_version >= 23:
v0 = shape.Vertex(0.0, 0.0, shape.SegmentType.MOVE)
v1 = shape.Vertex(1.0, 2.0, shape.SegmentType.LINE)
v2 = shape.Vertex(2.0, 3.0, shape.SegmentType.LINE)
v3 = shape.Vertex(3.0, 4.0, shape.SegmentType.LINE)
v4 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
vl = [v0, v1, v2, v3, v4]
samples = shape.MultiPolygon(vertices=vl)
p1 = shape.Point(0.0, 0.0)
p2 = shape.Point(1.0, 2.0)
p3 = shape.Point(2.0, 3.0)
p4 = shape.Point(3.0, 4.0)
p = [p1, p2, p3, p4]
polygon = shape.Polygon(points=p, samples=samples)
position.bounds = polygon
position.dimension = wcs.Dimension2D(10, 20)
position.resolution = 0.5
position.sample_size = 1.1
position.time_dependent = False
return position
def get_position(self):
position = plane.Position()
v1 = shape.Vertex(1.0, 2.0, shape.SegmentType.LINE)
v2 = shape.Vertex(2.0, 3.0, shape.SegmentType.LINE)
v3 = shape.Vertex(3.0, 4.0, shape.SegmentType.LINE)
v4 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
vl = [v1, v2, v3, v4]
polygon = shape.Polygon()
polygon.vertices = vl
position.bounds = polygon
position.dimension = wcs.Dimension2D(10, 20)
position.resolution = 0.5
position.sample_size = 1.1
position.time_dependent = False
return position
def test_open_polygon():
p1 = shape.Point(-117.246094, 52.942018)
p2 = shape.Point(-101.601563, 56.535258)
p3 = shape.Point(-97.382813, 44.809122)
p4 = shape.Point(-111.445313, 37.405074)
# SphericalPolygon requires p1 == p5 for a closed polygon
p5 = shape.Point(-117.246094, 52.942018)
too_few_points = [p1, p2]
min_closed_points = [p1, p2, p3]
closed_points = [p1, p2, p3, p4, p5]
counter_clockwise_points = [p4, p3, p2, p1]
# should detect that the polygons is not clockwise
with pytest.raises(AssertionError) as ex:
validate_polygon(shape.Polygon(counter_clockwise_points))
assert('not in clockwise direction' in str(ex.value))
# should detect that polygon is requires a minimum of 4 points
with pytest.raises(AssertionError) as ex:
validate_polygon(shape.Polygon(too_few_points))
assert('invalid polygon: 2 points' in str(ex.value))
# polygon default constructor
validate_polygon(shape.Polygon())
# should detect that polygon is closed
validate_polygon(shape.Polygon(min_closed_points))
validate_polygon(shape.Polygon(closed_points))
# should detect that multipolygon is not closed
v0 = shape.Vertex(-126.210938, 67.991108, shape.SegmentType.MOVE)
v1 = shape.Vertex(-108.984375, 70.480896, shape.SegmentType.LINE)
v2 = shape.Vertex(-98.789063, 66.912834, shape.SegmentType.LINE)
v3 = shape.Vertex(-75.234375, 60.217991, shape.SegmentType.LINE)
v4 = shape.Vertex(-87.890625, 52.241256, shape.SegmentType.LINE)
v5 = shape.Vertex(-110.742188, 54.136696, shape.SegmentType.LINE)
v6 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
v7 = shape.Vertex(24.609375, 62.895218, shape.SegmentType.MOVE)
v8 = shape.Vertex(43.593750, 67.322924, shape.SegmentType.LINE)
v9 = shape.Vertex(55.898438, 62.734601, shape.SegmentType.LINE)
v10 = shape.Vertex(46.757813, 56.145550, shape.SegmentType.LINE)
v11 = shape.Vertex(26.015625, 55.354135, shape.SegmentType.LINE)
v12 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
closed_vertices = [
v0, v1, v2, v3, v4, v5, v6, v7, v8, v9, v10, v11, v12]
# should detect that multipolygon is closed
validate_multipolygon(shape.MultiPolygon(closed_vertices))
def test_multipoly_self_intersect():
# should detect self segment intersection of the multipolygon not
# near a Pole
v1 = shape.Vertex(-115.488281, 45.867063, shape.SegmentType.MOVE)
v2 = shape.Vertex(-91.230469, 36.075742, shape.SegmentType.LINE)
v3 = shape.Vertex(-95.800781, 54.807017, shape.SegmentType.LINE)
v4 = shape.Vertex(-108.457031, 39.951859, shape.SegmentType.LINE)
v5 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
points_with_self_intersecting_segments = [v1, v2, v3, v4, v5]
with pytest.raises(AssertionError) as ex:
validate_multipolygon(
shape.MultiPolygon(points_with_self_intersecting_segments))
assert('self intersecting' in str(ex.value))
# should detect self segment intersection of the multipolygon near
# the South Pole, with the Pole outside the multipolygon
v1 = shape.Vertex(0.6128286003, -89.8967940441, shape.SegmentType.MOVE)
v2 = shape.Vertex(
210.6391743183, -89.9073892376, shape.SegmentType.LINE)
v3 = shape.Vertex(
90.6405151921, -89.8972874698, shape.SegmentType.LINE)
v4 = shape.Vertex(270.6114701911, -89.90689353, shape.SegmentType.LINE)
v5 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
points_with_self_intersecting_segments = [v1, v2, v3, v4, v5]
with pytest.raises(AssertionError) as ex:
validate_multipolygon(
shape.MultiPolygon(points_with_self_intersecting_segments))
assert('self intersecting' in str(ex.value))
# should detect self segment intersection of the multipolygon near the
# South Pole, with the Pole inside the multipolygon
v1 = shape.Vertex(0.6128286003, -89.8967940441, shape.SegmentType.MOVE)
v2 = shape.Vertex(
130.6391743183, -89.9073892376, shape.SegmentType.LINE)
v3 = shape.Vertex(
90.6405151921, -89.8972874698, shape.SegmentType.LINE)
v4 = shape.Vertex(270.6114701911, -89.90689353, shape.SegmentType.LINE)
v5 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
points_with_self_intersecting_segments = [v1, v2, v3, v4, v5]
with pytest.raises(AssertionError) as ex:
validate_multipolygon(
shape.MultiPolygon(points_with_self_intersecting_segments))
assert('self intersecting' in str(ex.value))
# should detect self segment intersection of the multipolygon which
# intersects with meridian = 0
v1 = shape.Vertex(-7.910156, 13.293411, shape.SegmentType.MOVE)
v2 = shape.Vertex(4.042969, 7.068185, shape.SegmentType.LINE)
v3 = shape.Vertex(4.746094, 18.030975, shape.SegmentType.LINE)
v4 = shape.Vertex(-6.855469, 6.369894, shape.SegmentType.LINE)
v5 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
points_with_self_intersecting_segments = [v1, v2, v3, v4, v5]
with pytest.raises(AssertionError) as ex:
validate_multipolygon(
shape.MultiPolygon(points_with_self_intersecting_segments))
assert('self intersecting' in str(ex.value))
def test_open_multipolygon():
# should detect that multipolygon is not closed
v0 = shape.Vertex(-126.210938, 67.991108, shape.SegmentType.MOVE)
v1 = shape.Vertex(-108.984375, 70.480896, shape.SegmentType.LINE)
v2 = shape.Vertex(-98.789063, 66.912834, shape.SegmentType.LINE)
v3 = shape.Vertex(-75.234375, 60.217991, shape.SegmentType.LINE)
v4 = shape.Vertex(-87.890625, 52.241256, shape.SegmentType.LINE)
v5 = shape.Vertex(-110.742188, 54.136696, shape.SegmentType.LINE)
v6 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
v7 = shape.Vertex(24.609375, 62.895218, shape.SegmentType.MOVE)
v8 = shape.Vertex(43.593750, 67.322924, shape.SegmentType.LINE)
v9 = shape.Vertex(55.898438, 62.734601, shape.SegmentType.LINE)
v10 = shape.Vertex(46.757813, 56.145550, shape.SegmentType.LINE)
v11 = shape.Vertex(26.015625, 55.354135, shape.SegmentType.LINE)
v12 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
no_vertices = []
too_few_vertices = [v0, v1, v6]
two_moves_vertices = [v0, v1, v7, v2, v3, v4, v5, v6]
no_move_vertices = [v1, v2, v3, v4, v5, v6]
two_closes_vertices = [
v0, v1, v2, v3, v4, v5, v7, v8, v9, v10, v11, v12]
no_close_vertices = [v0, v1, v2, v3, v4, v5]
min_closed_vertices = [v0, v1, v2, v6]
closed_vertices = [
v0, v1, v2, v3, v4, v5, v6, v7, v8, v9, v10, v11, v12]
rv0 = shape.Vertex(26.015625, 55.354135, shape.SegmentType.MOVE)
rv1 = shape.Vertex(46.757813, 56.145550, shape.SegmentType.LINE)
rv2 = shape.Vertex(55.898438, 62.734601, shape.SegmentType.LINE)
rv3 = shape.Vertex(43.593750, 67.322924, shape.SegmentType.LINE)
rv4 = shape.Vertex(24.609375, 62.895218, shape.SegmentType.LINE)
rv5 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
rv6 = shape.Vertex(-110.742188, 54.136696, shape.SegmentType.MOVE)
rv7 = shape.Vertex(-87.890625, 52.241256, shape.SegmentType.LINE)
rv8 = shape.Vertex(-75.234375, 60.217991, shape.SegmentType.LINE)
rv9 = shape.Vertex(-98.789063, 66.912834, shape.SegmentType.LINE)
rv10 = shape.Vertex(-108.984375, 70.480896, shape.SegmentType.LINE)
rv11 = shape.Vertex(-126.210938, 67.991108, shape.SegmentType.LINE)
rv12 = shape.Vertex(0.0, 0.0, shape.SegmentType.CLOSE)
counter_clockwise_vertices = [
rv0, rv1, rv2, rv3, rv4, rv5, rv6, rv7, rv8, rv9, rv10, rv11, rv12]
# should detect that the polygons is not clockwise
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(counter_clockwise_vertices))
assert('not in clockwise direction' in str(ex.value))
# should detect that there are not enough number of vertices to
# produce a multipolygon
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(no_vertices))
assert('invalid polygon: 0 vertices' in str(ex.value))
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(too_few_vertices))
assert('invalid polygon: 3 vertices' in str(ex.value))
# no close between two 'MOVE'
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(two_moves_vertices))
assert(
'invalid polygon: MOVE vertex when loop open' in str(ex.value))
# no 'MOVE' before a 'CLOSE'
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(no_move_vertices))
assert(
'invalid polygon: first vertex is not a MOVE' in str(ex.value))
# no 'MOVE' between two 'CLOSE'
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(two_closes_vertices))
assert(
'invalid polygon: MOVE vertex when loop open' in str(ex.value))
# no 'CLOSE' after a 'MOVE'
with pytest.raises(AssertionError) as ex:
validate_multipolygon(shape.MultiPolygon(no_close_vertices))
assert(
'invalid polygon: last vertex is not a CLOSE' in str(ex.value))
# multipolygon default constructor -> too few vertices
with pytest.raises(AssertionError):
validate_multipolygon(shape.MultiPolygon(None))
# should detect that multipolygon is closed
validate_multipolygon(shape.MultiPolygon(min_closed_vertices))
# should detect that multipolygon is closed
validate_multipolygon(shape.MultiPolygon(closed_vertices))
# instantiated multipolygon should contain the same vertices
p = shape.MultiPolygon(vertices=closed_vertices)
validate_multipolygon(p)
actual_vertices = p.vertices
assert(actual_vertices[0].cval1 == closed_vertices[0].cval1)
assert(actual_vertices[0].cval2 == closed_vertices[0].cval2)
assert(actual_vertices[0].type == shape.SegmentType.MOVE)
assert(actual_vertices[1].cval1 == closed_vertices[1].cval1)
assert(actual_vertices[1].cval2 == closed_vertices[1].cval2)
assert(actual_vertices[1].type == shape.SegmentType.LINE)
assert(actual_vertices[2].cval1 == closed_vertices[2].cval1)
assert(actual_vertices[2].cval2 == closed_vertices[2].cval2)
assert(actual_vertices[2].type == shape.SegmentType.LINE)
assert(actual_vertices[3].cval1 == closed_vertices[3].cval1)
assert(actual_vertices[3].cval2 == closed_vertices[3].cval2)
assert(actual_vertices[3].type == shape.SegmentType.LINE)
assert(actual_vertices[4].cval1 == closed_vertices[4].cval1)
assert(actual_vertices[4].cval2 == closed_vertices[4].cval2)
assert(actual_vertices[4].type == shape.SegmentType.LINE)
assert(actual_vertices[5].cval1 == closed_vertices[5].cval1)
assert(actual_vertices[5].cval2 == closed_vertices[5].cval2)
assert(actual_vertices[5].type == shape.SegmentType.LINE)
assert(actual_vertices[6].cval1 == closed_vertices[6].cval1)
assert(actual_vertices[6].cval2 == closed_vertices[6].cval2)
assert(actual_vertices[6].type == shape.SegmentType.CLOSE)
assert(actual_vertices[7].cval1 == closed_vertices[7].cval1)
assert(actual_vertices[7].cval2 == closed_vertices[7].cval2)
assert(actual_vertices[7].type == shape.SegmentType.MOVE)
assert(actual_vertices[8].cval1 == closed_vertices[8].cval1)
assert(actual_vertices[8].cval2 == closed_vertices[8].cval2)
assert(actual_vertices[8].type == shape.SegmentType.LINE)
assert(actual_vertices[9].cval1 == closed_vertices[9].cval1)
assert(actual_vertices[9].cval2 == closed_vertices[9].cval2)
assert(actual_vertices[9].type == shape.SegmentType.LINE)
assert(actual_vertices[10].cval1 == closed_vertices[10].cval1)
assert(actual_vertices[10].cval2 == closed_vertices[10].cval2)
assert(actual_vertices[10].type == shape.SegmentType.LINE)
assert(actual_vertices[11].cval1 == closed_vertices[11].cval1)
assert(actual_vertices[11].cval2 == closed_vertices[11].cval2)
assert(actual_vertices[11].type == shape.SegmentType.LINE)
assert(actual_vertices[12].cval1 == closed_vertices[12].cval1)
assert(actual_vertices[12].cval2 == closed_vertices[12].cval2)
assert(actual_vertices[12].type == shape.SegmentType.CLOSE)
