def test_line_cx_selection(gp_line, rect): x0, y0, x1, y1 = rect for xslice in get_slices(x0, x1): for yslice in get_slices(y0, y1): expected = LineArray.from_geopandas(gp_line.cx[xslice, yslice]) result = LineArray.from_geopandas(gp_line).cx[xslice, yslice] assert all(expected == result)
def test_line_intersects_rect(gp_line, rect): sg_rect = sg.box(*rect) expected = gp_line.intersects(sg_rect) lines = LineArray.from_geopandas(gp_line) # Test LineArray.intersects_rect result = lines.intersects_bounds(rect) np.testing.assert_equal(result, expected) # Test LineArray.intersects_rect with inds inds = np.flipud(np.arange(0, len(lines))) result = lines.intersects_bounds(rect, inds) np.testing.assert_equal(result, np.flipud(expected)) # Test Line.intersects_rect result = np.array([line.intersects_bounds(rect) for line in lines]) np.testing.assert_equal(result, expected)