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)