def test_ogr_copy_layer_by_name(self):
ds = gdal_create_dataset('Memory', 'test', gdal_type=gdal.OF_VECTOR)
georef.ogr_copy_layer_by_name(self.ds, 'agger_merge', ds)
self.assertTrue(isinstance(ds.GetLayer(), ogr.Layer))
self.assertRaises(
ValueError,
lambda: georef.ogr_copy_layer_by_name(self.ds, 'agger_merge1', ds))
def test_ogr_copy_layer_by_name(self):
filename = util.get_wradlib_data_file("shapefiles/agger/" "agger_merge.shp")
src_ds, layer = wradlib.io.open_vector(filename)
ds = wradlib.io.gdal_create_dataset("Memory", "test", gdal_type=gdal.OF_VECTOR)
georef.ogr_copy_layer_by_name(src_ds, "agger_merge", ds)
assert isinstance(ds.GetLayer(), ogr.Layer)
with pytest.raises(ValueError):
georef.ogr_copy_layer_by_name(src_ds, "agger_merge1", ds)
def test_ogr_copy_layer_by_name(self):
ds = gdal_create_dataset('Memory', 'test',
gdal_type=gdal.OF_VECTOR)
georef.ogr_copy_layer_by_name(self.ds, 'agger_merge', ds)
self.assertTrue(isinstance(ds.GetLayer(), ogr.Layer))
self.assertRaises(ValueError,
lambda: georef.ogr_copy_layer_by_name(self.ds,
'agger_merge1',
ds))
def load_vector(self, filename):
""" Load source/target grid points/polygons into in-memory Shapefile
Parameters
----------
filename : string
path to vector file
"""
# get input file handles
ds_in, tmp = io.open_shape(filename)
# create all DataSources
self.src = DataSource(name='src')
self.src.ds = io.gdal_create_dataset('Memory', 'src',
gdal_type=gdal.OF_VECTOR)
self.trg = DataSource(name='trg')
self.trg.ds = io.gdal_create_dataset('Memory', 'trg',
gdal_type=gdal.OF_VECTOR)
self.dst = DataSource(name='dst')
self.dst.ds = io.gdal_create_dataset('Memory', 'dst',
gdal_type=gdal.OF_VECTOR)
# copy all layers
georef.ogr_copy_layer_by_name(ds_in, "src", self.src.ds)
georef.ogr_copy_layer_by_name(ds_in, "trg", self.trg.ds)
georef.ogr_copy_layer_by_name(ds_in, "dst", self.dst.ds)
# get spatial reference object
self._srs = self.src.ds.GetLayer().GetSpatialRef()
self.src._srs = self.src.ds.GetLayer().GetSpatialRef()
self.trg._srs = self.trg.ds.GetLayer().GetSpatialRef()
self.dst._srs = self.trg.ds.GetLayer().GetSpatialRef()
# flush everything
del ds_in
def test_ogr_copy_layer_by_name(self):
ds = wradlib.io.gdal_create_dataset("Memory", "test", gdal_type=gdal.OF_VECTOR)
georef.ogr_copy_layer_by_name(self.ds, "agger_merge", ds)
assert isinstance(ds.GetLayer(), ogr.Layer)
with pytest.raises(ValueError):
georef.ogr_copy_layer_by_name(self.ds, "agger_merge1", ds)
