class TestWGS83TileGridUL(object):
def setup(self):
self.grid = TileGrid(4326, bbox=(-180, -90, 180, 90), origin='ul')
def test_resolution(self):
assert_almost_equal(self.grid.resolution(0), 1.40625)
assert_almost_equal(self.grid.resolution(1), 1.40625 / 2)
def test_bbox(self):
eq_(self.grid.bbox, (-180.0, -90.0, 180.0, 90.0))
def test_tile_bbox(self):
eq_(self.grid.tile_bbox((0, 0, 0)), (-180.0, -270.0, 180.0, 90.0))
eq_(self.grid.tile_bbox((0, 0, 0), limit=True),
(-180.0, -90.0, 180.0, 90.0))
eq_(self.grid.tile_bbox((0, 0, 1)), (-180.0, -90.0, 0.0, 90.0))
def test_tile(self):
eq_(self.grid.tile(-170, -80, 0), (0, 0, 0))
eq_(self.grid.tile(-170, -80, 1), (0, 0, 1))
eq_(self.grid.tile(-170, -80, 2), (0, 1, 2))
def test_grid_size(self):
eq_(self.grid.grid_sizes[0], (1, 1))
eq_(self.grid.grid_sizes[1], (2, 1))
eq_(self.grid.grid_sizes[2], (4, 2))
def test_affected_tiles(self):
bbox, grid, tiles = self.grid.get_affected_tiles((-180, -90, 180, 90),
(512, 256))
eq_(bbox, (-180.0, -90.0, 180.0, 90.0))
eq_(grid, (2, 1))
eq_(list(tiles), [(0, 0, 1), (1, 0, 1)])
bbox, grid, tiles = self.grid.get_affected_tiles((-180, -90, 0, 90),
(512, 512))
eq_(bbox, (-180.0, -90.0, 0.0, 90.0))
eq_(grid, (2, 2))
eq_(list(tiles), [(0, 0, 2), (1, 0, 2), (0, 1, 2), (1, 1, 2)])
def test_affected_level_tiles(self):
bbox, grid, tiles = self.grid.get_affected_level_tiles(
(-180, -90, 180, 90), 1)
eq_(grid, (2, 1))
eq_(bbox, (-180.0, -90.0, 180.0, 90.0))
eq_(list(tiles), [(0, 0, 1), (1, 0, 1)])
bbox, grid, tiles = self.grid.get_affected_level_tiles((0, 0, 180, 90),
2)
eq_(grid, (2, 1))
eq_(list(tiles), [(2, 0, 2), (3, 0, 2)])
eq_(bbox, (0.0, 0.0, 180.0, 90.0))
bbox, grid, tiles = self.grid.get_affected_level_tiles(
(0, -90, 180, 90), 2)
eq_(grid, (2, 2))
eq_(list(tiles), [(2, 0, 2), (3, 0, 2), (2, 1, 2), (3, 1, 2)])
eq_(bbox, (0.0, -90.0, 180.0, 90.0))
class TestGeodeticTileGrid(TileGridTest):
def setup(self):
self.grid = TileGrid(is_geodetic=True, )
def test_auto_resolution(self):
grid = TileGrid(is_geodetic=True,
bbox=(-10, 30, 10, 40),
tile_size=(20, 20))
tile_bbox = grid.tile_bbox((0, 0, 0))
assert tile_bbox == (-10, 30, 10, 50)
assert grid.resolution(0) == 1.0
def test_grid(self):
for level, grid_size in [(0, (1, 1)), (1, (2, 1)), (2, (4, 2))]:
yield self.check_grid, level, grid_size
def test_adjacent_tile_bbox(self):
grid = TileGrid(is_geodetic=True,
bbox=(-10, 30, 10, 40),
tile_size=(20, 20))
t1 = grid.tile_bbox((0, 0, 2))
t2 = grid.tile_bbox((1, 0, 2))
t3 = grid.tile_bbox((0, 1, 2))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[2] == t3[2]
assert t1[3] == t3[1]
def test_w_resolution(self):
res = [1, 0.5, 0.2]
grid = TileGrid(is_geodetic=True,
bbox=(-10, 30, 10, 40),
tile_size=(20, 20),
res=res)
assert grid.grid_sizes[0] == (1, 1)
assert grid.grid_sizes[1] == (2, 1)
assert grid.grid_sizes[2] == (5, 3)
def test_tile(self):
assert self.grid.tile(-180, -90, 0) == (0, 0, 0)
assert self.grid.tile(-180, -90, 1) == (0, 0, 1)
assert self.grid.tile(-180, -90, 2) == (0, 0, 2)
assert self.grid.tile(180 - 0.001, 90 - 0.001, 0) == (0, 0, 0)
assert self.grid.tile(10, 50, 1) == (1, 0, 1)
def test_affected_tiles(self):
bbox, grid_size, tiles = \
self.grid.get_affected_tiles((-45,-45,45,45), (512,512))
assert self.grid.grid_sizes[3] == (8, 4)
assert bbox == (-45.0, -45.0, 45.0, 45.0)
assert grid_size == (2, 2)
tiles = list(tiles)
assert tiles == [(3, 2, 3), (4, 2, 3), (3, 1, 3), (4, 1, 3)]
class TestGKTileGridUL(TileGridTest):
"""
Custom grid with ul origin.
"""
def setup(self):
self.grid = TileGrid(SRS(31467),
bbox=(3300000, 5300000, 3900000, 6000000),
origin='ul',
res=[1500, 1000, 500, 300, 150, 100])
def test_bbox(self):
assert self.grid.bbox == (3300000, 5300000, 3900000, 6000000)
def test_tile_bbox(self):
eq_(self.grid.tile_bbox((0, 0, 0)),
(3300000.0, 5616000.0, 3684000.0, 6000000.0))
eq_(self.grid.tile_bbox((1, 0, 0)),
(3684000.0, 5616000.0, 4068000.0, 6000000.0))
eq_(self.grid.tile_bbox((1, 1, 0)),
(3684000.0, 5232000.0, 4068000.0, 5616000.0))
def test_tile(self):
x, y = 3310000, 5990000
eq_(self.grid.tile(x, y, 0), (0, 0, 0))
eq_(self.grid.tile(x, y, 1), (0, 0, 1))
eq_(self.grid.tile(x, y, 2), (0, 0, 2))
x, y = 3890000, 5310000
eq_(self.grid.tile(x, y, 0), (1, 1, 0))
eq_(self.grid.tile(x, y, 1), (2, 2, 1))
eq_(self.grid.tile(x, y, 2), (4, 5, 2))
def test_grids(self):
assert_grid_size(self.grid, 0, (2, 2))
assert_grid_size(self.grid, 1, (3, 3))
assert_grid_size(self.grid, 2, (5, 6))
def test_closest_level(self):
assert self.grid.closest_level(1500) == 0
assert self.grid.closest_level(1000) == 1
assert self.grid.closest_level(900) == 1
assert self.grid.closest_level(600) == 2
def test_adjacent_tile_bbox(self):
t1 = self.grid.tile_bbox((0, 0, 1))
t2 = self.grid.tile_bbox((1, 0, 1))
t3 = self.grid.tile_bbox((0, 1, 1))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[1] == t3[3]
class TestGKTileGridUL(TileGridTest):
"""
Custom grid with ul origin.
"""
def setup(self):
self.grid = TileGrid(SRS(31467),
bbox=(3300000, 5300000, 3900000, 6000000), origin='ul',
res=[1500, 1000, 500, 300, 150, 100])
def test_bbox(self):
assert self.grid.bbox == (3300000, 5300000, 3900000, 6000000)
def test_tile_bbox(self):
eq_(self.grid.tile_bbox((0, 0, 0)),
(3300000.0, 5616000.0, 3684000.0, 6000000.0))
eq_(self.grid.tile_bbox((1, 0, 0)),
(3684000.0, 5616000.0, 4068000.0, 6000000.0))
eq_(self.grid.tile_bbox((1, 1, 0)),
(3684000.0, 5232000.0, 4068000.0, 5616000.0))
def test_tile(self):
x, y = 3310000, 5990000
eq_(self.grid.tile(x, y, 0), (0, 0, 0))
eq_(self.grid.tile(x, y, 1), (0, 0, 1))
eq_(self.grid.tile(x, y, 2), (0, 0, 2))
x, y = 3890000, 5310000
eq_(self.grid.tile(x, y, 0), (1, 1, 0))
eq_(self.grid.tile(x, y, 1), (2, 2, 1))
eq_(self.grid.tile(x, y, 2), (4, 5, 2))
def test_grids(self):
assert_grid_size(self.grid, 0, (2, 2))
assert_grid_size(self.grid, 1, (3, 3))
assert_grid_size(self.grid, 2, (5, 6))
def test_closest_level(self):
assert self.grid.closest_level(1500) == 0
assert self.grid.closest_level(1000) == 1
assert self.grid.closest_level(900) == 1
assert self.grid.closest_level(600) == 2
def test_adjacent_tile_bbox(self):
t1 = self.grid.tile_bbox((0, 0, 1))
t2 = self.grid.tile_bbox((1, 0, 1))
t3 = self.grid.tile_bbox((0, 1, 1))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[1] == t3[3]
class TestWGS83TileGridUL(object):
def setup(self):
self.grid = TileGrid(4326, bbox=(-180, -90, 180, 90), origin='ul')
def test_resolution(self):
assert_almost_equal(self.grid.resolution(0), 1.40625)
assert_almost_equal(self.grid.resolution(1), 1.40625/2)
def test_bbox(self):
eq_(self.grid.bbox, (-180.0, -90.0, 180.0, 90.0))
def test_tile_bbox(self):
eq_(self.grid.tile_bbox((0, 0, 0)), (-180.0, -270.0, 180.0, 90.0))
eq_(self.grid.tile_bbox((0, 0, 0), limit=True), (-180.0, -90.0, 180.0, 90.0))
eq_(self.grid.tile_bbox((0, 0, 1)), (-180.0, -90.0, 0.0, 90.0))
def test_tile(self):
eq_(self.grid.tile(-170, -80, 0), (0, 0, 0))
eq_(self.grid.tile(-170, -80, 1), (0, 0, 1))
eq_(self.grid.tile(-170, -80, 2), (0, 1, 2))
def test_grid_size(self):
eq_(self.grid.grid_sizes[0], (1, 1))
eq_(self.grid.grid_sizes[1], (2, 1))
eq_(self.grid.grid_sizes[2], (4, 2))
def test_affected_tiles(self):
bbox, grid, tiles = self.grid.get_affected_tiles((-180,-90,180,90), (512,256))
eq_(bbox, (-180.0, -90.0, 180.0, 90.0))
eq_(grid, (2, 1))
eq_(list(tiles), [(0, 0, 1), (1, 0, 1)])
bbox, grid, tiles = self.grid.get_affected_tiles((-180,-90,0,90), (512, 512))
eq_(bbox, (-180.0, -90.0, 0.0, 90.0))
eq_(grid, (2, 2))
eq_(list(tiles), [(0, 0, 2), (1, 0, 2), (0, 1, 2), (1, 1, 2)])
def test_affected_level_tiles(self):
bbox, grid, tiles = self.grid.get_affected_level_tiles((-180,-90,180,90), 1)
eq_(grid, (2, 1))
eq_(bbox, (-180.0, -90.0, 180.0, 90.0))
eq_(list(tiles), [(0, 0, 1), (1, 0, 1)])
bbox, grid, tiles = self.grid.get_affected_level_tiles((0,0,180,90), 2)
eq_(grid, (2, 1))
eq_(list(tiles), [(2, 0, 2), (3, 0, 2)])
eq_(bbox, (0.0, 0.0, 180.0, 90.0))
bbox, grid, tiles = self.grid.get_affected_level_tiles((0,-90,180,90), 2)
eq_(grid, (2, 2))
eq_(list(tiles), [(2, 0, 2), (3, 0, 2), (2, 1, 2), (3, 1, 2)])
eq_(bbox, (0.0, -90.0, 180.0, 90.0))
class TestGeodeticTileGrid(TileGridTest):
def setup(self):
self.grid = TileGrid(is_geodetic=True, )
def test_auto_resolution(self):
grid = TileGrid(is_geodetic=True, bbox=(-10, 30, 10, 40), tile_size=(20, 20))
tile_bbox = grid.tile_bbox((0, 0, 0))
assert tile_bbox == (-10, 30, 10, 50)
assert grid.resolution(0) == 1.0
def test_grid(self):
for level, grid_size in [(0, (1, 1)), (1, (2, 1)), (2, (4, 2))]:
yield self.check_grid, level, grid_size
def test_adjacent_tile_bbox(self):
grid = TileGrid(is_geodetic=True, bbox=(-10, 30, 10, 40), tile_size=(20, 20))
t1 = grid.tile_bbox((0, 0, 2))
t2 = grid.tile_bbox((1, 0, 2))
t3 = grid.tile_bbox((0, 1, 2))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[2] == t3[2]
assert t1[3] == t3[1]
def test_w_resolution(self):
res = [1, 0.5, 0.2]
grid = TileGrid(is_geodetic=True, bbox=(-10, 30, 10, 40), tile_size=(20, 20), res=res)
assert grid.grid_sizes[0] == (1, 1)
assert grid.grid_sizes[1] == (2, 1)
assert grid.grid_sizes[2] == (5, 3)
def test_tile(self):
assert self.grid.tile(-180, -90, 0) == (0, 0, 0)
assert self.grid.tile(-180, -90, 1) == (0, 0, 1)
assert self.grid.tile(-180, -90, 2) == (0, 0, 2)
assert self.grid.tile(180-0.001, 90-0.001, 0) == (0, 0, 0)
assert self.grid.tile(10, 50, 1) == (1, 0, 1)
def test_affected_tiles(self):
bbox, grid_size, tiles = \
self.grid.get_affected_tiles((-45,-45,45,45), (512,512))
assert self.grid.grid_sizes[3] == (8, 4)
assert bbox == (-45.0, -45.0, 45.0, 45.0)
assert grid_size == (2, 2)
tiles = list(tiles)
assert tiles == [(3, 2, 3), (4, 2, 3), (3, 1, 3), (4, 1, 3)]
class TestGKTileGrid(TileGridTest):
def setup(self):
self.grid = TileGrid(SRS(31467),
bbox=(3250000, 5230000, 3930000, 6110000))
def test_bbox(self):
assert self.grid.bbox == (3250000, 5230000, 3930000, 6110000)
def test_resolution(self):
res = self.grid.resolution(0)
width = self.grid.bbox[2] - self.grid.bbox[0]
height = self.grid.bbox[3] - self.grid.bbox[1]
assert height == 880000.0 and width == 680000.0
assert res == 880000.0 / 256
def test_tile_bbox(self):
tile_bbox = self.grid.tile_bbox((0, 0, 0))
assert tile_bbox == (3250000.0, 5230000.0, 4130000.0, 6110000.0)
def test_tile(self):
x, y = 3450000, 5890000
assert [self.grid.tile(x, y, level) for level in range(5)] == \
[(0, 0, 0), (0, 1, 1), (0, 3, 2), (1, 6, 3), (3, 12, 4)]
def test_grids(self):
for level, grid_size in [(0, (1, 1)), (1, (2, 2)), (2, (4, 4)),
(3, (7, 8))]:
yield self.check_grid, level, grid_size
def test_closest_level(self):
assert self.grid.closest_level(880000.0 / 256) == 0
assert self.grid.closest_level(600000.0 / 256) == 1
assert self.grid.closest_level(440000.0 / 256) == 1
assert self.grid.closest_level(420000.0 / 256) == 1
def test_adjacent_tile_bbox(self):
t1 = self.grid.tile_bbox((0, 0, 1))
t2 = self.grid.tile_bbox((1, 0, 1))
t3 = self.grid.tile_bbox((0, 1, 1))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[3] == t3[1]
class TestGKTileGrid(TileGridTest):
def setup(self):
self.grid = TileGrid(SRS(31467), bbox=(3250000, 5230000, 3930000, 6110000))
def test_bbox(self):
assert self.grid.bbox == (3250000, 5230000, 3930000, 6110000)
def test_resolution(self):
res = self.grid.resolution(0)
width = self.grid.bbox[2] - self.grid.bbox[0]
height = self.grid.bbox[3] - self.grid.bbox[1]
assert height == 880000.0 and width == 680000.0
assert res == 880000.0/256
def test_tile_bbox(self):
tile_bbox = self.grid.tile_bbox((0, 0, 0))
assert tile_bbox == (3250000.0, 5230000.0, 4130000.0, 6110000.0)
def test_tile(self):
x, y = 3450000, 5890000
assert [self.grid.tile(x, y, level) for level in range(5)] == \
[(0, 0, 0), (0, 1, 1), (0, 3, 2), (1, 6, 3), (3, 12, 4)]
def test_grids(self):
for level, grid_size in [(0, (1, 1)), (1, (2, 2)), (2, (4, 4)), (3, (7, 8))]:
yield self.check_grid, level, grid_size
def test_closest_level(self):
assert self.grid.closest_level(880000.0/256) == 0
assert self.grid.closest_level(600000.0/256) == 1
assert self.grid.closest_level(440000.0/256) == 1
assert self.grid.closest_level(420000.0/256) == 1
def test_adjacent_tile_bbox(self):
t1 = self.grid.tile_bbox((0, 0, 1))
t2 = self.grid.tile_bbox((1, 0, 1))
t3 = self.grid.tile_bbox((0, 1, 1))
assert t1[1] == t2[1]
assert t1[3] == t2[3]
assert t1[2] == t2[0]
assert t1[0] == t3[0]
assert t1[3] == t3[1]
def test_tile_out_of_grid_bounds(self):
grid = TileGrid(is_geodetic=True)
eq_(grid.tile(-180.01, 50, 1), (-1, 0, 1))
def test_tile_out_of_grid_bounds(self):
grid = TileGrid(is_geodetic=True)
eq_(grid.tile(-180.01, 50, 1), (-1, 0, 1))
