def test_explicit_grid(self):
grid = TileGrid(res=[0.1, 0.05, 0.01])
eq_(grid.resolution(0), 0.1)
eq_(grid.resolution(1), 0.05)
eq_(grid.resolution(2), 0.01)
eq_(grid.closest_level(0.00001), 2)
def test_explicit_grid(self):
grid = TileGrid(res=[0.1, 0.05, 0.01])
eq_(grid.resolution(0), 0.1)
eq_(grid.resolution(1), 0.05)
eq_(grid.resolution(2), 0.01)
eq_(grid.closest_level(0.00001), 2)
class TestWGS84TileGrid(object):
def setup(self):
self.grid = TileGrid(is_geodetic=True)
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_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)])
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_(bbox, (0.0, 0.0, 180.0, 90.0))
eq_(list(tiles), [(2, 1, 2), (3, 1, 2)])
class TestWGS84TileGrid(object):
def setup(self):
self.grid = TileGrid(is_geodetic=True)
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_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)])
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_(bbox, (0.0, 0.0, 180.0, 90.0))
eq_(list(tiles), [(2, 1, 2), (3, 1, 2)])
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
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 TestFixedResolutionsTileGrid(TileGridTest):
def setup(self):
self.res = [1000.0, 500.0, 200.0, 100.0, 50.0, 20.0, 5.0]
bbox = (3250000, 5230000, 3930000, 6110000)
self.grid = TileGrid(SRS(31467), bbox=bbox, res=self.res)
def test_resolution(self):
for level, res in enumerate(self.res):
assert res == self.grid.resolution(level)
def test_closest_level(self):
assert self.grid.closest_level(2000) == 0
assert self.grid.closest_level(1000) == 0
assert self.grid.closest_level(950) == 0
assert self.grid.closest_level(210) == 2
def test_affected_tiles(self):
req_bbox = (3250000, 5230000, 3930000, 6110000)
self.grid.max_shrink_factor = 10
bbox, grid_size, tiles = \
self.grid.get_affected_tiles(req_bbox, (256, 256))
assert bbox == (req_bbox[0], req_bbox[1],
req_bbox[0]+1000*256*3, req_bbox[1]+1000*256*4)
assert grid_size == (3, 4)
tiles = list(tiles)
assert tiles == [(0, 3, 0), (1, 3, 0), (2, 3, 0),
(0, 2, 0), (1, 2, 0), (2, 2, 0),
(0, 1, 0), (1, 1, 0), (2, 1, 0),
(0, 0, 0), (1, 0, 0), (2, 0, 0),
]
def test_affected_tiles_2(self):
req_bbox = (3250000, 5230000, 3930000, 6110000)
self.grid.max_shrink_factor = 2.0
try:
bbox, grid_size, tiles = \
self.grid.get_affected_tiles(req_bbox, (256, 256))
except NoTiles:
pass
else:
assert False, 'got no exception'
def test_grid(self):
for level, grid_size in [(0, (3, 4)), (1, (6, 7)), (2, (14, 18))]:
yield self.check_grid, level, grid_size
def test_tile_bbox(self):
tile_bbox = self.grid.tile_bbox((0, 0, 0)) # w: 1000x256
assert tile_bbox == (3250000.0, 5230000.0, 3506000.0, 5486000.0)
tile_bbox = self.grid.tile_bbox((0, 0, 1)) # w: 500x256
assert tile_bbox == (3250000.0, 5230000.0, 3378000.0, 5358000.0)
tile_bbox = self.grid.tile_bbox((0, 0, 2)) # w: 200x256
assert tile_bbox == (3250000.0, 5230000.0, 3301200.0, 5281200.0)
class TestFixedResolutionsTileGrid(TileGridTest):
def setup(self):
self.res = [1000.0, 500.0, 200.0, 100.0, 50.0, 20.0, 5.0]
bbox = (3250000, 5230000, 3930000, 6110000)
self.grid = TileGrid(SRS(31467), bbox=bbox, res=self.res)
def test_resolution(self):
for level, res in enumerate(self.res):
assert res == self.grid.resolution(level)
def test_closest_level(self):
assert self.grid.closest_level(2000) == 0
assert self.grid.closest_level(1000) == 0
assert self.grid.closest_level(950) == 0
assert self.grid.closest_level(210) == 2
def test_affected_tiles(self):
req_bbox = (3250000, 5230000, 3930000, 6110000)
self.grid.max_shrink_factor = 10
bbox, grid_size, tiles = \
self.grid.get_affected_tiles(req_bbox, (256, 256))
assert bbox == (req_bbox[0], req_bbox[1],
req_bbox[0]+1000*256*3, req_bbox[1]+1000*256*4)
assert grid_size == (3, 4)
tiles = list(tiles)
assert tiles == [(0, 3, 0), (1, 3, 0), (2, 3, 0),
(0, 2, 0), (1, 2, 0), (2, 2, 0),
(0, 1, 0), (1, 1, 0), (2, 1, 0),
(0, 0, 0), (1, 0, 0), (2, 0, 0),
]
def test_affected_tiles_2(self):
req_bbox = (3250000, 5230000, 3930000, 6110000)
self.grid.max_shrink_factor = 2.0
try:
bbox, grid_size, tiles = \
self.grid.get_affected_tiles(req_bbox, (256, 256))
except NoTiles:
pass
else:
assert False, 'got no exception'
def test_grid(self):
for level, grid_size in [(0, (3, 4)), (1, (6, 7)), (2, (14, 18))]:
yield self.check_grid, level, grid_size
def test_tile_bbox(self):
tile_bbox = self.grid.tile_bbox((0, 0, 0)) # w: 1000x256
assert tile_bbox == (3250000.0, 5230000.0, 3506000.0, 5486000.0)
tile_bbox = self.grid.tile_bbox((0, 0, 1)) # w: 500x256
assert tile_bbox == (3250000.0, 5230000.0, 3378000.0, 5358000.0)
tile_bbox = self.grid.tile_bbox((0, 0, 2)) # w: 200x256
assert tile_bbox == (3250000.0, 5230000.0, 3301200.0, 5281200.0)
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 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]
class MapProxyConfiguration(object):
def __init__(self, srs, target_dir):
self.service_srs = srs
self.sources = {}
self.caches = {}
self.layers = []
self.yaml_file = os.path.join(target_dir, 'mapproxy.yaml')
self.grid = TileGrid(SRS(3857))
def _load_sources(self):
public_wmts = db.session.query(WMTS, ST_Transform(WMTS.view_coverage, 3857)).filter_by(is_public=True).group_by(WMTS).all()
for wmts, view_coverage in public_wmts:
self.sources['%s_source' % wmts.name] = {
'type': 'tile',
'url': wmts.url,
'grid': 'GoogleMapsCompatible',
'coverage': {
'srs': 'EPSG:3857',
'bbox': list(to_shape(view_coverage).bounds)
}
}
self.caches['%s_cache' % wmts.name] = {
'sources': ['%s_source' % wmts.name],
'grids': ['GoogleMapsCompatible'],
'disable_storage': True
}
self.layers.append({
'name': '%s_layer' % wmts.name,
'title': wmts.title,
'sources': ['%s_cache' % wmts.name],
'min_res': self.grid.resolution(wmts.view_level_start),
# increase max_res to allow a bit of oversampling
'max_res': self.grid.resolution(wmts.view_level_end) / 2,
})
def _write_mapproxy_yaml(self):
grids = {
'GoogleMapsCompatible': {
'base': 'GLOBAL_MERCATOR',
'srs': 'EPSG:3857',
'num_levels': 20,
'origin': 'nw'
}
}
services = {
'demo': None,
'wms': {
'srs': self.service_srs,
'md': {'title': _('external_wms_title')}
},
}
globals_ = {
'image': {
'paletted': True,
},
'cache': {
'meta_size': [8, 8],
'meta_buffer': 50,
},
'http': {
'client_timeout': 120,
},
}
config = {}
if globals_: config['globals'] = globals_
if grids: config['grids'] = grids
if self.layers:
config['layers'] = self.layers
config['services'] = services
if self.caches: config['caches'] = self.caches
if self.sources: config['sources'] = self.sources
# safe_dump does not output !!python/unicode, etc.
mapproxy_yaml = yaml.safe_dump(config, default_flow_style=False)
f = open(self.yaml_file, 'w')
f.write(mapproxy_yaml)
f.close()
def test_sqrt_grid(self):
grid = TileGrid(is_geodetic=True, res='sqrt2', tile_size=(360, 180))
eq_(grid.resolution(0), 1.0)
assert_almost_equal(grid.resolution(2), 0.5)
assert_almost_equal(grid.resolution(4), 0.25)
def test_factor_grid(self):
grid = TileGrid(is_geodetic=True, res=1/0.75, tile_size=(360, 180))
eq_(grid.resolution(0), 1.0)
eq_(grid.resolution(1), 0.75)
eq_(grid.resolution(2), 0.75*0.75)
class MapProxyConfiguration(object):
def __init__(self, db_session, srs, target_dir, couch_url, couch_port):
self.db_session = db_session
self.couch_url = couch_url
self.couch_port = couch_port
self.service_srs = srs
self.sources = {}
self.caches = {}
self.layers = []
self.yaml_file = os.path.join(target_dir, 'mapproxy.yaml')
self.grid = TileGrid(SRS(3857))
def _load_sources(self):
local_sources = self.db_session.query(LocalWMTSSource).all()
for local_source in local_sources:
wmts_source = local_source.wmts_source
self.sources[wmts_source.name + '_source'] = {
'type': 'wms',
'req': {
'url': 'http://dummy.example.org/service?'
},
'seed_only': True,
'coverage': {
'srs':
'EPSG:3857',
'bbox':
list(
coverage_from_geojson(
wmts_source.download_coverage).bbox)
}
}
self.caches[wmts_source.name + '_cache'] = {
'sources': [wmts_source.name + '_source'],
'grids': [wmts_source.matrix_set],
'cache': {
'type': 'couchdb',
'url': '%s:%s' % (self.couch_url, self.couch_port),
'db_name': wmts_source.layer,
'tile_metadata': {
'tile_col': '{{x}}',
'tile_row': '{{y}}',
'tile_level': '{{z}}',
'created_ts': '{{timestamp}}',
'created': '{{utc_iso}}',
'center': '{{wgs_tile_centroid}}'
}
}
}
self.layers.append({
'name':
wmts_source.name + '_layer',
'title':
wmts_source.title,
'sources': [wmts_source.name + '_cache'],
'min_res':
self.grid.resolution(local_source.download_level_start),
# increase max_res to allow a bit of oversampling
'max_res':
self.grid.resolution(local_source.download_level_end) / 2,
})
def _write_mapproxy_yaml(self):
grids = {
'GoogleMapsCompatible': {
'base': 'GLOBAL_MERCATOR',
'srs': 'EPSG:3857',
'num_levels': 19,
'origin': 'nw'
}
}
services = {
'demo': None,
'wmts': None,
'wms': {
'srs': self.service_srs,
'md': {
'title': 'Geobox WMS'
}
},
}
globals_ = {
'image': {
'paletted': True,
},
'cache': {
'meta_size': [8, 8],
'meta_buffer': 50,
},
'http': {
'client_timeout': 120,
},
}
config = {}
if globals_: config['globals'] = globals_
if grids: config['grids'] = grids
if self.layers:
config['layers'] = self.layers
config['services'] = services
if self.caches: config['caches'] = self.caches
if self.sources: config['sources'] = self.sources
# safe_dump does not output !!python/unicode, etc.
mapproxy_yaml = yaml.safe_dump(config, default_flow_style=False)
f = open(self.yaml_file, 'w')
f.write(mapproxy_yaml)
f.close()
log.info('Mapproxy configuration written to %s', self.yaml_file)
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_factor_grid(self):
grid = TileGrid(is_geodetic=True, res=1 / 0.75, tile_size=(360, 180))
eq_(grid.resolution(0), 1.0)
eq_(grid.resolution(1), 0.75)
eq_(grid.resolution(2), 0.75 * 0.75)
def test_sqrt_grid(self):
grid = TileGrid(is_geodetic=True, res='sqrt2', tile_size=(360, 180))
eq_(grid.resolution(0), 1.0)
assert_almost_equal(grid.resolution(2), 0.5)
assert_almost_equal(grid.resolution(4), 0.25)
class MapProxyConfiguration(object):
def __init__(self, db_session, srs, target_dir, couchdb_url, template_dir):
self.db_session = db_session
self.couchdb_url = couchdb_url
self.service_srs = srs
self.template_dir = template_dir
self.sources = {}
self.caches = {}
self.layers = []
self.yaml_file = os.path.join(target_dir, 'mapproxy.yaml')
self.grid = TileGrid(SRS(3857))
def _load_sources(self):
local_sources = self.db_session.query(LocalWMTSSource).all()
for local_source in local_sources:
wmts_source = local_source.wmts_source
self.sources[wmts_source.name + '_source'] = {
'type': 'wms',
'req': {
'url': 'http://dummy.example.org/service?'
},
'seed_only': True,
'coverage': {
'srs': 'EPSG:3857',
'bbox': list(coverage_from_geojson(wmts_source.download_coverage).bbox)
}
}
self.caches[wmts_source.name + '_cache'] = {
'sources': [wmts_source.name + '_source'],
'grids': [wmts_source.matrix_set],
'cache': {
'type': 'couchdb',
'url': '%s' % self.couchdb_url,
'db_name': wmts_source.name,
'tile_metadata': {
'tile_col': '{{x}}',
'tile_row': '{{y}}',
'tile_level': '{{z}}',
'created_ts': '{{timestamp}}',
'created': '{{utc_iso}}',
'center': '{{wgs_tile_centroid}}'
}
}
}
self.layers.append({
'name': wmts_source.name + '_layer',
'title': wmts_source.title,
'sources': [wmts_source.name + '_cache'],
'min_res': self.grid.resolution(local_source.download_level_start),
# increase max_res to allow a bit of oversampling
'max_res': self.grid.resolution(local_source.download_level_end) / 2,
})
def _write_mapproxy_yaml(self):
grids = {
'GoogleMapsCompatible': {
'base': 'GLOBAL_MERCATOR',
'srs': 'EPSG:3857',
'num_levels': 19,
'origin': 'nw'
}
}
services = {
'demo': None,
'wmts': None,
'wms': {
'srs': self.service_srs,
'md': {'title': 'Geobox WMS'}
},
}
globals_ = {
'image': {
'paletted': True,
},
'cache': {
'meta_size': [8, 8],
'meta_buffer': 50,
},
'http': {
'client_timeout': 120,
},
}
if self.template_dir:
globals_['template_dir'] = self.template_dir
config = {}
if globals_: config['globals'] = globals_
if grids: config['grids'] = grids
if self.layers:
config['layers'] = self.layers
config['services'] = services
if self.caches: config['caches'] = self.caches
if self.sources: config['sources'] = self.sources
# safe_dump does not output !!python/unicode, etc.
mapproxy_yaml = yaml.safe_dump(config, default_flow_style=False)
f = open(self.yaml_file, 'w')
f.write(mapproxy_yaml)
f.close()
log.info('Mapproxy configuration written to %s', self.yaml_file)
class MapProxyConfiguration(object):
def __init__(self, srs, target_dir):
self.service_srs = srs
self.sources = {}
self.caches = {}
self.layers = []
self.yaml_file = os.path.join(target_dir, 'mapproxy.yaml')
self.grid = TileGrid(SRS(3857))
def _load_sources(self):
public_wmts = db.session.query(
WMTS,
ST_Transform(WMTS.view_coverage,
3857)).filter_by(is_public=True).group_by(WMTS).all()
for wmts, view_coverage in public_wmts:
self.sources['%s_source' % wmts.name] = {
'type': 'tile',
'url': wmts.url,
'grid': 'GoogleMapsCompatible',
'coverage': {
'srs': 'EPSG:3857',
'bbox': list(to_shape(view_coverage).bounds)
}
}
self.caches['%s_cache' % wmts.name] = {
'sources': ['%s_source' % wmts.name],
'grids': ['GoogleMapsCompatible'],
'disable_storage': True
}
self.layers.append({
'name':
'%s_layer' % wmts.name,
'title':
wmts.title,
'sources': ['%s_cache' % wmts.name],
'min_res':
self.grid.resolution(wmts.view_level_start),
# increase max_res to allow a bit of oversampling
'max_res':
self.grid.resolution(wmts.view_level_end) / 2,
})
def _write_mapproxy_yaml(self):
grids = {
'GoogleMapsCompatible': {
'base': 'GLOBAL_MERCATOR',
'srs': 'EPSG:3857',
'num_levels': 20,
'origin': 'nw'
}
}
services = {
'demo': None,
'wms': {
'srs': self.service_srs,
'md': {
'title': _('external_wms_title')
}
},
}
globals_ = {
'image': {
'paletted': True,
},
'cache': {
'meta_size': [8, 8],
'meta_buffer': 50,
},
'http': {
'client_timeout': 120,
},
}
config = {}
if globals_: config['globals'] = globals_
if grids: config['grids'] = grids
if self.layers:
config['layers'] = self.layers
config['services'] = services
if self.caches: config['caches'] = self.caches
if self.sources: config['sources'] = self.sources
# safe_dump does not output !!python/unicode, etc.
mapproxy_yaml = yaml.safe_dump(config, default_flow_style=False)
f = open(self.yaml_file, 'w')
f.write(mapproxy_yaml)
f.close()
