def setup(self):
# box from 10 10 to 80 80 with small spike/corner to -10 60 (upper left)
self.geom = shapely.wkt.loads(
"POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
self.coverage1 = coverage(self.geom, SRS(4326))
self.coverage2 = coverage([100, 0, 120, 20], SRS(4326))
self.coverage = MultiCoverage([self.coverage1, self.coverage2])
def test_eq(self):
g1 = shapely.wkt.loads("POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
g2 = shapely.wkt.loads("POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
assert coverage(g1, SRS(4326)) == coverage(g2, SRS(4326))
assert coverage(g1, SRS(4326)) != coverage(g2, SRS(31467))
g3 = shapely.wkt.loads("POLYGON((10.0 10, 10 50.0, -10 60, 10 80, 80 80, 80 10, 10 10))")
assert coverage(g1, SRS(4326)) == coverage(g3, SRS(4326))
g4 = shapely.wkt.loads("POLYGON((10 10, 10.1 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
assert coverage(g1, SRS(4326)) != coverage(g4, SRS(4326))
def test_eq(self):
g1 = shapely.wkt.loads("POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
g2 = shapely.wkt.loads("POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
assert MultiCoverage([coverage(g1, SRS(4326))]) == MultiCoverage([coverage(g2, SRS(4326))])
assert MultiCoverage([coverage(g1, SRS(4326))]) != MultiCoverage([coverage(g2, SRS(31467))])
c = coverage([-10, 10, 80, 80], SRS(4326))
assert MultiCoverage([c, coverage(g1, SRS(4326))]) != MultiCoverage([c, coverage(g2, SRS(31467))])
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 test_eq(self):
assert coverage([-10, 10, 80, 80], SRS(4326)) == coverage([-10, 10, 80, 80], SRS(4326))
assert coverage([-10, 10, 80, 80], SRS(4326)) == coverage([-10, 10.0, 80.0, 80], SRS(4326))
assert coverage([-10, 10, 80, 80], SRS(4326)) != coverage([-10.1, 10.0, 80.0, 80], SRS(4326))
assert coverage([-10, 10, 80, 80], SRS(4326)) != coverage([-10, 10.0, 80.0, 80], SRS(31467))
def setup(self):
self.coverage = coverage([-10, 10, 80, 80], SRS(4326))
def setup(self):
# box from 10 10 to 80 80 with small spike/corner to -10 60 (upper left)
self.geom = shapely.wkt.loads(
"POLYGON((10 10, 10 50, -10 60, 10 80, 80 80, 80 10, 10 10))")
self.coverage = coverage(self.geom, SRS(4326))
