def load_coverage(conf):
if "ogr_datasource" in conf:
require_geom_support()
srs = conf["ogr_srs"]
datasource = conf["ogr_datasource"]
if not re.match(r"^\w{2,}:", datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource)
where = conf.get("ogr_where", None)
geom = load_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif "polygons" in conf:
require_geom_support()
srs = conf["polygons_srs"]
geom = load_polygons(abspath(conf["polygons"]))
bbox, geom = build_multipolygon(geom, simplify=True)
else:
srs = conf["bbox_srs"]
bbox = conf["bbox"]
if isinstance(bbox, basestring):
bbox = map(float, bbox.split(","))
geom = None
return coverage(geom or bbox, SRS(srs))
def load_coverage(conf):
if 'ogr_datasource' in conf:
require_geom_support()
srs = conf['ogr_srs']
datasource = conf['ogr_datasource']
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource)
where = conf.get('ogr_where', None)
geom = load_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'polygons' in conf:
require_geom_support()
srs = conf['polygons_srs']
geom = load_polygons(abspath(conf['polygons']))
bbox, geom = build_multipolygon(geom, simplify=True)
else:
srs = conf['bbox_srs']
bbox = conf['bbox']
if isinstance(bbox, basestring):
bbox = map(float, bbox.split(','))
geom = None
return coverage(geom or bbox, SRS(srs))
def test_wkt(self):
with TempFile() as fname:
with open(fname, 'wb') as f:
f.write(VALID_POLYGON1)
f.write(b"\n")
f.write(VALID_POLYGON1)
geoms = load_datasource(fname)
eq_(len(geoms), 2)
def test_wkt(self):
with TempFile() as fname:
with open(fname, 'wb') as f:
f.write(VALID_POLYGON1)
f.write(b"\n")
f.write(VALID_POLYGON1)
geoms = load_datasource(fname)
eq_(len(geoms), 2)
def test_geojson(self):
with TempFile() as fname:
with open(fname, 'wb') as f:
f.write(b'''{"type": "FeatureCollection", "features": [
{"type": "Feature", "geometry": {"type": "Polygon", "coordinates": [[[0, 0], [10, 0], [10, 10], [0, 0]]]} },
{"type": "Feature", "geometry": {"type": "MultiPolygon", "coordinates": [[[[0, 0], [10, 0], [10, 10], [0, 0]]], [[[0, 0], [10, 0], [10, 10], [0, 0]]], [[[0, 0], [10, 0], [10, 10], [0, 0]]]]} },
{"type": "Feature", "geometry": {"type": "Point", "coordinates": [0, 0]} }
]}''')
geoms = load_datasource(fname)
assert len(geoms) == 4
def test_geojson(self):
with TempFile() as fname:
with open(fname, 'wb') as f:
f.write(b'''{"type": "FeatureCollection", "features": [
{"type": "Feature", "geometry": {"type": "Polygon", "coordinates": [[[0, 0], [10, 0], [10, 10], [0, 0]]]} },
{"type": "Feature", "geometry": {"type": "MultiPolygon", "coordinates": [[[[0, 0], [10, 0], [10, 10], [0, 0]]], [[[0, 0], [10, 0], [10, 10], [0, 0]]], [[[0, 0], [10, 0], [10, 10], [0, 0]]]]} },
{"type": "Feature", "geometry": {"type": "Point", "coordinates": [0, 0]} }
]}''')
geoms = load_datasource(fname)
eq_(len(geoms), 4)
def load_coverage(conf, base_path=None):
if 'ogr_datasource' in conf:
require_geom_support()
srs = conf['ogr_srs']
datasource = conf['ogr_datasource']
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('ogr_where', None)
geom = load_ogr_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'polygons' in conf:
require_geom_support()
srs = conf['polygons_srs']
geom = load_polygons(abspath(conf['polygons'], base_path=base_path))
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'bbox' in conf:
srs = conf.get('bbox_srs') or conf['srs']
bbox = conf['bbox']
if isinstance(bbox, string_type):
bbox = [float(x) for x in bbox.split(',')]
geom = None
elif 'datasource' in conf:
require_geom_support()
datasource = conf['datasource']
srs = conf['srs']
if isinstance(datasource, (list, tuple)):
bbox = datasource
geom = None
elif bbox_string_re.match(datasource):
bbox = [float(x) for x in datasource.split(',')]
geom = None
else:
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('where', None)
geom = load_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
else:
return None
return coverage(geom or bbox, SRS(srs))
def load_coverage(conf):
if 'ogr_datasource' in conf:
require_geom_support()
srs = conf['ogr_srs']
datasource = conf['ogr_datasource']
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource)
where = conf.get('ogr_where', None)
bbox, geom = load_datasource(datasource, where)
elif 'polygons' in conf:
require_geom_support()
srs = conf['polygons_srs']
bbox, geom = load_polygons(abspath(conf['polygons']))
else:
srs = conf['bbox_srs']
bbox = conf['bbox']
if isinstance(bbox, basestring):
bbox = map(float, bbox.split(','))
geom = None
return coverage(geom or bbox, SRS(srs))
def load_coverage(conf, base_path=None):
clip = False
if 'clip' in conf:
clip = conf['clip']
if 'union' in conf:
parts = []
for cov in conf['union']:
parts.append(load_coverage(cov))
return union_coverage(parts, clip=clip)
elif 'intersection' in conf:
parts = []
for cov in conf['intersection']:
parts.append(load_coverage(cov))
return intersection_coverage(parts, clip=clip)
elif 'difference' in conf:
parts = []
for cov in conf['difference']:
parts.append(load_coverage(cov))
return diff_coverage(parts, clip=clip)
elif 'ogr_datasource' in conf:
require_geom_support()
srs = conf['ogr_srs']
datasource = conf['ogr_datasource']
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('ogr_where', None)
geom = load_ogr_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'polygons' in conf:
require_geom_support()
srs = conf['polygons_srs']
geom = load_polygons(abspath(conf['polygons'], base_path=base_path))
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'bbox' in conf:
srs = conf.get('bbox_srs') or conf['srs']
bbox = conf['bbox']
if isinstance(bbox, string_type):
bbox = [float(x) for x in bbox.split(',')]
geom = None
elif 'datasource' in conf:
require_geom_support()
datasource = conf['datasource']
srs = conf['srs']
if isinstance(datasource, (list, tuple)):
bbox = datasource
geom = None
elif bbox_string_re.match(datasource):
bbox = [float(x) for x in datasource.split(',')]
geom = None
else:
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('where', None)
geom = load_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'expire_tiles' in conf:
require_geom_support()
filename = abspath(conf['expire_tiles'])
geom = load_expire_tiles(filename)
_, geom = build_multipolygon(geom, simplify=False)
return coverage(geom, SRS(3857))
else:
return None
return coverage(geom or bbox, SRS(srs), clip=clip)
def test_shp(self):
polygon_file = os.path.join(os.path.dirname(__file__), 'polygons',
'polygons.shp')
geoms = load_datasource(polygon_file)
eq_(len(geoms), 3)
def load_coverage(conf, base_path=None):
clip = False
if 'clip' in conf:
clip = conf['clip']
if 'union' in conf:
parts = []
for cov in conf['union']:
parts.append(load_coverage(cov))
return union_coverage(parts, clip=clip)
elif 'intersection' in conf:
parts = []
for cov in conf['intersection']:
parts.append(load_coverage(cov))
return intersection_coverage(parts, clip=clip)
elif 'difference' in conf:
parts = []
for cov in conf['difference']:
parts.append(load_coverage(cov))
return diff_coverage(parts, clip=clip)
elif 'ogr_datasource' in conf:
require_geom_support()
srs = conf['ogr_srs']
datasource = conf['ogr_datasource']
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('ogr_where', None)
geom = load_ogr_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'polygons' in conf:
require_geom_support()
srs = conf['polygons_srs']
geom = load_polygons(abspath(conf['polygons'], base_path=base_path))
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'bbox' in conf:
srs = conf.get('bbox_srs') or conf['srs']
bbox = conf['bbox']
if isinstance(bbox, string_type):
bbox = [float(x) for x in bbox.split(',')]
geom = None
elif 'datasource' in conf:
require_geom_support()
datasource = conf['datasource']
srs = conf['srs']
if isinstance(datasource, (list, tuple)):
bbox = datasource
geom = None
elif bbox_string_re.match(datasource):
bbox = [float(x) for x in datasource.split(',')]
geom = None
else:
if not re.match(r'^\w{2,}:', datasource):
# looks like a file and not PG:, MYSQL:, etc
# make absolute path
datasource = abspath(datasource, base_path=base_path)
where = conf.get('where', None)
geom = load_datasource(datasource, where)
bbox, geom = build_multipolygon(geom, simplify=True)
elif 'expire_tiles' in conf:
require_geom_support()
filename = abspath(conf['expire_tiles'])
geom = load_expire_tiles(filename)
_, geom = build_multipolygon(geom, simplify=False)
return coverage(geom, SRS(3857))
else:
return None
return coverage(geom or bbox, SRS(srs), clip=clip)
def test_shp(self):
polygon_file = os.path.join(os.path.dirname(__file__), 'polygons', 'polygons.shp')
geoms = load_datasource(polygon_file)
eq_(len(geoms), 3)
