def acreate_vector(self,
path,
geometry,
fields=(),
layer=None,
driver='ESRI Shapefile',
options=(),
sr=None):
"""Create a vector file anonymously in this DataSource. Only metadata are kept in memory.
See DataSource.create_vector
Example
-------
>>> lines = ds.acreate_vector('/path/to.shp', 'linestring')
>>> file_proj4 = lines.proj4_stored
"""
# Parameter checking ***************************************************
path = str(path)
geometry = conv.str_of_wkbgeom(conv.wkbgeom_of_str(geometry))
fields = _tools.normalize_fields_defn(fields)
if layer is None:
layer = '.'.join(ntpath.basename(path).split('.')[:-1])
else:
layer = str(layer)
driver = str(driver)
options = [str(arg) for arg in options]
if sr is not None:
sr = osr.GetUserInputAsWKT(sr)
# Construction dispatch ************************************************
if driver.lower() == 'memory':
# TODO: Check not concurrent
allocator = lambda: BackGDALFileVector.create_file(
'', geometry, fields, layer, 'Memory', options, sr)
prox = GDALMemoryVector(self, allocator, options)
elif True:
allocator = lambda: BackGDALFileVector.create_file(
path, geometry, fields, layer, driver, options, sr)
prox = GDALFileVector(self, allocator, options, 'w')
else:
pass
# DataSource Registering ***********************************************
self._register([], prox)
return prox
def aopen_vector(self,
path,
layer=None,
driver='ESRI Shapefile',
options=(),
mode='r'):
"""Open a vector file anonymously in this DataSource. Only metadata are kept in memory.
See DataSource.open_vector
Example
-------
>>> trees = ds.aopen_vector('/path/to.shp')
>>> features_bounds = trees.bounds
"""
path = str(path)
if layer is None:
layer = 0
elif isinstance(layer, numbers.Integral):
layer = int(layer)
else:
layer = str(layer)
driver = str(driver)
options = [str(arg) for arg in options]
_ = conv.of_of_mode(mode)
# Construction dispatch ************************************************
if driver.lower() == 'memory': # pragma: no cover
raise ValueError("Can't open a MEMORY vector, user create_vector")
elif True:
allocator = lambda: BackGDALFileVector.open_file(
path, layer, driver, options, mode)
prox = GDALFileVector(self, allocator, options, mode)
else:
pass
# DataSource Registering ***********************************************
self._register([], prox)
return prox
def create_vector(self,
key,
path,
geometry,
fields=(),
layer=None,
driver='ESRI Shapefile',
options=(),
sr=None):
"""Create a vector file and register it under `key` in this DataSource. Only metadata are
kept in memory.
Parameters
----------
key: hashable (like a string)
File identifier within DataSource
path: string
geometry: string
name of a wkb geometry type
http://www.gdal.org/ogr__core_8h.html#a800236a0d460ef66e687b7b65610f12a
(see example below)
fields: sequence of dict
Attributes of fields, one dict per field. (see `Field attributes` below)
layer: None or string
driver: string
ogr driver to use when opening the file
http://www.gdal.org/ogr_formats.html
options: sequence of str
options for ogr
sr: string or None
Spatial reference of the new file
if None: don't set a spatial reference
if string:
if path: Use same projection as file at `path`
if textual spatial reference:
http://gdal.org/java/org/gdal/osr/SpatialReference.html#SetFromUserInput-java.lang.String-
Returns
-------
one of {GDALFileVector, GDALMemoryVector} depending on the `driver` parameter
Field attributes
----------------
Attributes:
'name': string
'type': string (see `Field type` below)
'precision': int
'width': int
'nullable': bool
'default': same as `type`
An attribute missing or None is kept to default value.
Field types
-----------
Binary key: 'binary', bytes, np.bytes_, aliases of np.bytes_
Date key: 'date'
DateTime key: 'datetime', datetime.datetime, np.datetime64, aliases of np.datetime64
Time key: 'time'
Integer key: 'integer' np.int32, aliases of np.int32
Integer64 key: 'integer64', int, np.int64, aliases of np.int64
Real key: 'real', float, np.float64, aliases of np.float64
String key: 'string', str, np.str_, aliases of np.str_
Integer64List key: 'integer64list', 'int list'
IntegerList key: 'integerlist'
RealList key: 'reallist', 'float list'
StringList key: 'stringlist', 'str list'
Example
-------
>>> ds.create_vector('lines', '/path/to.shp', 'linestring')
>>> geometry_type = ds.lines.type
>>> fields = [
{'name': 'name', 'type': str},
{'name': 'count', 'type': 'int32'},
{'name': 'area', 'type': np.float64, 'width': 5, precision: 18},
{'name': 'when', 'type': np.datetime64},
]
>>> ds.create_vector('zones', '/path/to.shp', 'polygon', fields)
>>> field0_type = ds.zones.fields[0]['type']
"""
# Parameter checking ***************************************************
self._validate_key(key)
path = str(path)
geometry = conv.str_of_wkbgeom(conv.wkbgeom_of_str(geometry))
fields = _tools.normalize_fields_defn(fields)
if layer is None:
layer = '.'.join(ntpath.basename(path).split('.')[:-1])
else:
layer = str(layer)
driver = str(driver)
options = [str(arg) for arg in options]
if sr is not None:
sr = osr.GetUserInputAsWKT(sr)
# Construction dispatch ************************************************
if driver.lower() == 'memory':
# TODO: Check not concurrent
allocator = lambda: BackGDALFileVector.create_file(
'', geometry, fields, layer, 'Memory', options, sr)
prox = GDALMemoryVector(self, allocator, options)
elif True:
allocator = lambda: BackGDALFileVector.create_file(
path, geometry, fields, layer, driver, options, sr)
prox = GDALFileVector(self, allocator, options, 'w')
else:
pass
# DataSource Registering ***********************************************
self._register([key], prox)
return prox
def open_vector(self,
key,
path,
layer=None,
driver='ESRI Shapefile',
options=(),
mode='r'):
"""Open a vector file in this DataSource under `key`. Only metadata are kept in memory.
Parameters
----------
key: hashable (like a string)
File identifier within DataSource
path: string
layer: None or int or string
driver: string
ogr driver to use when opening the file
http://www.gdal.org/ogr_formats.html
options: sequence of str
options for ogr
mode: one of {'r', 'w'}
Returns
-------
GDALFileVector
Example
-------
>>> ds.open_vector('trees', '/path/to.shp')
>>> feature_count = len(ds.trees)
>>> ds.open_vector('roofs', '/path/to.json', driver='GeoJSON', mode='w')
>>> fields_list = ds.roofs.fields
"""
# Parameter checking ***************************************************
self._validate_key(key)
path = str(path)
if layer is None:
layer = 0
elif isinstance(layer, numbers.Integral):
layer = int(layer)
else:
layer = str(layer)
driver = str(driver)
options = [str(arg) for arg in options]
_ = conv.of_of_mode(mode)
# Construction dispatch ************************************************
if driver.lower() == 'memory': # pragma: no cover
raise ValueError("Can't open a MEMORY vector, user create_vector")
elif True:
allocator = lambda: BackGDALFileVector.open_file(
path, layer, driver, options, mode)
prox = GDALFileVector(self, allocator, options, mode)
else:
pass
# DataSource Registering ***********************************************
self._register([key], prox)
return prox