Source code for grass.temporal.metadata
"""
Metadata classes for map layer and space time datasets
Usage:
.. code-block:: python
>>> import grass.temporal as tgis
>>> tgis.init()
>>> meta = tgis.RasterMetadata()
>>> meta = tgis.Raster3DMetadata()
>>> meta = tgis.VectorMetadata()
>>> meta = tgis.STRDSMetadata()
>>> meta = tgis.STR3DSMetadata()
>>> meta = tgis.STVDSMetadata()
(C) 2012-2013 by the GRASS Development Team
This program is free software under the GNU General Public
License (>=v2). Read the file COPYING that comes with GRASS
for details.
:authors: Soeren Gebbert
"""
from .base import SQLDatabaseInterface
from .core import SQLDatabaseInterfaceConnection, get_tgis_db_version_from_metadata
###############################################################################
[docs]class RasterMetadataBase(SQLDatabaseInterface):
"""This is the metadata base class for time stamped raster and raster3d maps
Usage:
.. code-block:: python
>>> init()
>>> meta = RasterMetadataBase(
... table="metadata",
... ident="soil@PERMANENT",
... datatype="CELL",
... cols=100,
... rows=100,
... number_of_cells=10000,
... nsres=0.1,
... ewres=0.1,
... min=0,
... max=100,
... )
>>> meta.datatype
'CELL'
>>> meta.cols
100
>>> meta.rows
100
>>> meta.number_of_cells
10000
>>> meta.nsres
0.1
>>> meta.ewres
0.1
>>> meta.min
0.0
>>> meta.max
100.0
>>> meta.print_info()
| Datatype:................... CELL
| Number of columns:.......... 100
| Number of rows:............. 100
| Number of cells:............ 10000
| North-South resolution:..... 0.1
| East-west resolution:....... 0.1
| Minimum value:.............. 0.0
| Maximum value:.............. 100.0
>>> meta.print_shell_info()
datatype=CELL
cols=100
rows=100
number_of_cells=10000
nsres=0.1
ewres=0.1
min=0.0
max=100.0
"""
def __init__(
self,
table=None,
ident=None,
datatype=None,
cols=None,
rows=None,
number_of_cells=None,
nsres=None,
ewres=None,
min=None,
max=None,
):
SQLDatabaseInterface.__init__(self, table, ident)
self.set_id(ident)
self.set_datatype(datatype)
self.set_cols(cols)
self.set_rows(rows)
self.set_number_of_cells(number_of_cells)
self.set_nsres(nsres)
self.set_ewres(ewres)
self.set_min(min)
self.set_max(max)
[docs] def set_id(self, ident):
"""Convenient method to set the unique identifier (primary key)"""
self.ident = ident
self.D["id"] = ident
[docs] def set_cols(self, cols):
"""Set the number of cols"""
if cols is not None:
self.D["cols"] = int(cols)
else:
self.D["cols"] = None
[docs] def set_rows(self, rows):
"""Set the number of rows"""
if rows is not None:
self.D["rows"] = int(rows)
else:
self.D["rows"] = None
[docs] def set_number_of_cells(self, number_of_cells):
"""Set the number of cells"""
if number_of_cells is not None:
self.D["number_of_cells"] = int(number_of_cells)
else:
self.D["number_of_cells"] = None
[docs] def set_nsres(self, nsres):
"""Set the north-south resolution"""
if nsres is not None:
self.D["nsres"] = float(nsres)
else:
self.D["nsres"] = None
[docs] def set_ewres(self, ewres):
"""Set the east-west resolution"""
if ewres is not None:
self.D["ewres"] = float(ewres)
else:
self.D["ewres"] = None
[docs] def set_min(self, min):
"""Set the minimum raster value"""
if min is not None:
self.D["min"] = float(min)
else:
self.D["min"] = None
[docs] def set_max(self, max):
"""Set the maximum raster value"""
if max is not None:
self.D["max"] = float(max)
else:
self.D["max"] = None
[docs] def get_id(self):
"""Convenient method to get the unique identifier (primary key)
:return: None if not found
"""
if "id" in self.D:
return self.D["id"]
else:
return None
[docs] def get_datatype(self):
"""Get the map type
:return: None if not found"""
if "datatype" in self.D:
return self.D["datatype"]
else:
return None
[docs] def get_cols(self):
"""Get number of cols
:return: None if not found"""
if "cols" in self.D:
return self.D["cols"]
else:
return None
[docs] def get_rows(self):
"""Get number of rows
:return: None if not found"""
if "rows" in self.D:
return self.D["rows"]
else:
return None
[docs] def get_number_of_cells(self):
"""Get number of cells
:return: None if not found"""
if "number_of_cells" in self.D:
return self.D["number_of_cells"]
else:
return None
[docs] def get_nsres(self):
"""Get the north-south resolution
:return: None if not found"""
if "nsres" in self.D:
return self.D["nsres"]
else:
return None
[docs] def get_ewres(self):
"""Get east-west resolution
:return: None if not found"""
if "ewres" in self.D:
return self.D["ewres"]
else:
return None
[docs] def get_min(self):
"""Get the minimum cell value
:return: None if not found"""
if "min" in self.D:
return self.D["min"]
else:
return None
[docs] def get_max(self):
"""Get the maximum cell value
:return: None if not found"""
if "max" in self.D:
return self.D["max"]
else:
return None
# Properties
datatype = property(fget=get_datatype, fset=set_datatype)
cols = property(fget=get_cols, fset=set_cols)
rows = property(fget=get_rows, fset=set_rows)
number_of_cells = property(fget=get_number_of_cells, fset=set_number_of_cells)
nsres = property(fget=get_nsres, fset=set_nsres)
ewres = property(fget=get_ewres, fset=set_ewres)
min = property(fget=get_min, fset=set_min)
max = property(fget=get_max, fset=set_max)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
self._print_info_body(shell=False)
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_body(shell=True)
def _print_info_head(self, shell=False):
"""Print information about this class (head part).
No header printed in shell style mode.
:param bool shell: True for human readable style otherwise shell style
"""
if not shell:
print(
" +-------------------- Metadata information ----------------------------------+" # noqa: E501
)
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("datatype=" + str(self.get_datatype()))
print("cols=" + str(self.get_cols()))
print("rows=" + str(self.get_rows()))
print("number_of_cells=" + str(self.get_number_of_cells()))
print("nsres=" + str(self.get_nsres()))
print("ewres=" + str(self.get_ewres()))
print("min=" + str(self.get_min()))
print("max=" + str(self.get_max()))
else:
print(" | Datatype:................... " + str(self.get_datatype()))
print(" | Number of columns:.......... " + str(self.get_cols()))
print(" | Number of rows:............. " + str(self.get_rows()))
print(" | Number of cells:............ " + str(self.get_number_of_cells()))
print(" | North-South resolution:..... " + str(self.get_nsres()))
print(" | East-west resolution:....... " + str(self.get_ewres()))
print(" | Minimum value:.............. " + str(self.get_min()))
print(" | Maximum value:.............. " + str(self.get_max()))
###############################################################################
[docs]class RasterMetadata(RasterMetadataBase):
"""This is the raster metadata class
This class is the interface to the raster_metadata table in the
temporal database that stores the metadata of all registered raster maps.
The metadata includes the datatype, number of cols, rows and cells and
the north-south and east west resolution of the map. Additionally the
minimum and maximum valuesare stored.
Usage:
.. code-block:: python
>>> init()
>>> meta = RasterMetadata(
... ident="soil@PERMANENT",
... datatype="CELL",
... cols=100,
... rows=100,
... number_of_cells=10000,
... nsres=0.1,
... ewres=0.1,
... min=0,
... max=100,
... )
>>> meta.datatype
'CELL'
>>> meta.cols
100
>>> meta.rows
100
>>> meta.number_of_cells
10000
>>> meta.nsres
0.1
>>> meta.ewres
0.1
>>> meta.min
0.0
>>> meta.max
100.0
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| Datatype:................... CELL
| Number of columns:.......... 100
| Number of rows:............. 100
| Number of cells:............ 10000
| North-South resolution:..... 0.1
| East-west resolution:....... 0.1
| Minimum value:.............. 0.0
| Maximum value:.............. 100.0
| Semantic label:............. None
>>> meta.print_shell_info()
datatype=CELL
cols=100
rows=100
number_of_cells=10000
nsres=0.1
ewres=0.1
min=0.0
max=100.0
semantic_label=None
"""
def __init__(
self,
ident=None,
datatype=None,
cols=None,
rows=None,
number_of_cells=None,
nsres=None,
ewres=None,
min=None,
max=None,
semantic_label=None,
):
RasterMetadataBase.__init__(
self,
"raster_metadata",
ident,
datatype,
cols,
rows,
number_of_cells,
nsres,
ewres,
min,
max,
)
[docs] def set_semantic_label(self, semantic_label):
"""Set the semantic label identifier"""
self.D["semantic_label"] = semantic_label
[docs] def get_semantic_label(self):
"""Get the semantic label identifier
:return: None if not found"""
if "semantic_label" in self.D:
return self.D["semantic_label"]
else:
return None
semantic_label = property(fget=get_semantic_label, fset=set_semantic_label)
[docs] def print_info(self):
"""Print information about this class."""
self._print_info_head(shell=False)
self._print_info_body(shell=False)
# semantic label section (raster specific only)
print(" | Semantic label:............. " + str(self.get_semantic_label()))
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_head(shell=True)
self._print_info_body(shell=True)
# semantic label section (raster specific only)
print("semantic_label=" + str(self.get_semantic_label()))
###############################################################################
[docs]class Raster3DMetadata(RasterMetadataBase):
"""This is the raster3d metadata class
This class is the interface to the raster3d_metadata table in the
temporal database that stores the metadata of all registered
3D raster maps.
The metadata includes all raster metadata variables and additional
the number of depths, the top-bottom resolution and the space time 3D
raster dataset register table is stored.
Usage:
.. code-block:: python
>>> init()
>>> meta = Raster3DMetadata(
... ident="soil@PERMANENT",
... datatype="FCELL",
... cols=100,
... rows=100,
... depths=100,
... number_of_cells=1000000,
... nsres=0.1,
... ewres=0.1,
... tbres=0.1,
... min=0,
... max=100,
... )
>>> meta.datatype
'FCELL'
>>> meta.cols
100
>>> meta.rows
100
>>> meta.depths
100
>>> meta.number_of_cells
1000000
>>> meta.nsres
0.1
>>> meta.ewres
0.1
>>> meta.tbres
0.1
>>> meta.min
0.0
>>> meta.max
100.0
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| Datatype:................... FCELL
| Number of columns:.......... 100
| Number of rows:............. 100
| Number of cells:............ 1000000
| North-South resolution:..... 0.1
| East-west resolution:....... 0.1
| Minimum value:.............. 0.0
| Maximum value:.............. 100.0
| Number of depths:........... 100
| Top-Bottom resolution:...... 0.1
>>> meta.print_shell_info()
datatype=FCELL
cols=100
rows=100
number_of_cells=1000000
nsres=0.1
ewres=0.1
min=0.0
max=100.0
depths=100
tbres=0.1
"""
def __init__(
self,
ident=None,
datatype=None,
cols=None,
rows=None,
depths=None,
number_of_cells=None,
nsres=None,
ewres=None,
tbres=None,
min=None,
max=None,
):
RasterMetadataBase.__init__(
self,
"raster3d_metadata",
ident,
datatype,
cols,
rows,
number_of_cells,
nsres,
ewres,
min,
max,
)
self.set_tbres(tbres)
self.set_depths(depths)
[docs] def set_depths(self, depths):
"""Set the number of depths"""
if depths is not None:
self.D["depths"] = int(depths)
else:
self.D["depths"] = None
[docs] def set_tbres(self, tbres):
"""Set the top-bottom resolution"""
if tbres is not None:
self.D["tbres"] = float(tbres)
else:
self.D["tbres"] = None
[docs] def get_depths(self):
"""Get number of depths
:return: None if not found"""
if "depths" in self.D:
return self.D["depths"]
else:
return None
[docs] def get_tbres(self):
"""Get top-bottom resolution
:return: None if not found"""
if "tbres" in self.D:
return self.D["tbres"]
else:
return None
depths = property(fget=get_depths, fset=set_depths)
tbres = property(fget=get_tbres, fset=set_tbres)
[docs] def print_info(self):
"""Print information about this class."""
self._print_info_head(shell=False)
self._print_info_body(shell=False)
print(" | Number of depths:........... " + str(self.get_depths()))
print(" | Top-Bottom resolution:...... " + str(self.get_tbres()))
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_head(shell=True)
self._print_info_body(shell=True)
print("depths=" + str(self.get_depths()))
print("tbres=" + str(self.get_tbres()))
###############################################################################
[docs]class VectorMetadata(SQLDatabaseInterface):
"""This is the vector metadata class
This class is the interface to the vector_metadata table in the
temporal database that stores the metadata of all registered
vector maps.
Usage:
.. code-block:: python
>>> init()
>>> meta = VectorMetadata(
... ident="lidar@PERMANENT",
... is_3d=True,
... number_of_points=1,
... number_of_lines=2,
... number_of_boundaries=3,
... number_of_centroids=4,
... number_of_faces=5,
... number_of_kernels=6,
... number_of_primitives=7,
... number_of_nodes=8,
... number_of_areas=9,
... number_of_islands=10,
... number_of_holes=11,
... number_of_volumes=12,
... )
>>> meta.id
'lidar@PERMANENT'
>>> meta.is_3d
True
>>> meta.number_of_points
1
>>> meta.number_of_lines
2
>>> meta.number_of_boundaries
3
>>> meta.number_of_centroids
4
>>> meta.number_of_faces
5
>>> meta.number_of_kernels
6
>>> meta.number_of_primitives
7
>>> meta.number_of_nodes
8
>>> meta.number_of_areas
9
>>> meta.number_of_islands
10
>>> meta.number_of_holes
11
>>> meta.number_of_volumes
12
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| Is map 3d .................. True
| Number of points ........... 1
| Number of lines ............ 2
| Number of boundaries ....... 3
| Number of centroids ........ 4
| Number of faces ............ 5
| Number of kernels .......... 6
| Number of primitives ....... 7
| Number of nodes ............ 8
| Number of areas ............ 9
| Number of islands .......... 10
| Number of holes ............ 11
| Number of volumes .......... 12
>>> meta.print_shell_info()
is_3d=True
points=1
lines=2
boundaries=3
centroids=4
faces=5
kernels=6
primitives=7
nodes=8
areas=9
islands=10
holes=11
volumes=12
"""
def __init__(
self,
ident=None,
is_3d=False,
number_of_points=None,
number_of_lines=None,
number_of_boundaries=None,
number_of_centroids=None,
number_of_faces=None,
number_of_kernels=None,
number_of_primitives=None,
number_of_nodes=None,
number_of_areas=None,
number_of_islands=None,
number_of_holes=None,
number_of_volumes=None,
):
SQLDatabaseInterface.__init__(self, "vector_metadata", ident)
self.set_id(ident)
self.set_3d_info(is_3d)
self.set_number_of_points(number_of_points)
self.set_number_of_lines(number_of_lines)
self.set_number_of_boundaries(number_of_boundaries)
self.set_number_of_centroids(number_of_centroids)
self.set_number_of_faces(number_of_faces)
self.set_number_of_kernels(number_of_kernels)
self.set_number_of_primitives(number_of_primitives)
self.set_number_of_nodes(number_of_nodes)
self.set_number_of_areas(number_of_areas)
self.set_number_of_islands(number_of_islands)
self.set_number_of_holes(number_of_holes)
self.set_number_of_volumes(number_of_volumes)
[docs] def set_id(self, ident):
"""Convenient method to set the unique identifier (primary key)"""
self.ident = ident
self.D["id"] = ident
[docs] def set_3d_info(self, is_3d):
"""Set True if the vector map is three dimensional"""
self.D["is_3d"] = is_3d
[docs] def set_number_of_points(self, number_of_points):
"""Set the number of points of the vector map"""
self.D["points"] = number_of_points
[docs] def set_number_of_lines(self, number_of_lines):
"""Set the number of lines of the vector map"""
self.D["lines"] = number_of_lines
[docs] def set_number_of_boundaries(self, number_of_boundaries):
"""Set the number of boundaries of the vector map"""
self.D["boundaries"] = number_of_boundaries
[docs] def set_number_of_centroids(self, number_of_centroids):
"""Set the number of centroids of the vector map"""
self.D["centroids"] = number_of_centroids
[docs] def set_number_of_faces(self, number_of_faces):
"""Set the number of faces of the vector map"""
self.D["faces"] = number_of_faces
[docs] def set_number_of_kernels(self, number_of_kernels):
"""Set the number of kernels of the vector map"""
self.D["kernels"] = number_of_kernels
[docs] def set_number_of_primitives(self, number_of_primitives):
"""Set the number of primitives of the vector map"""
self.D["primitives"] = number_of_primitives
[docs] def set_number_of_nodes(self, number_of_nodes):
"""Set the number of nodes of the vector map"""
self.D["nodes"] = number_of_nodes
[docs] def set_number_of_areas(self, number_of_areas):
"""Set the number of areas of the vector map"""
self.D["areas"] = number_of_areas
[docs] def set_number_of_islands(self, number_of_islands):
"""Set the number of islands of the vector map"""
self.D["islands"] = number_of_islands
[docs] def set_number_of_holes(self, number_of_holes):
"""Set the number of holes of the vector map"""
self.D["holes"] = number_of_holes
[docs] def set_number_of_volumes(self, number_of_volumes):
"""Set the number of volumes of the vector map"""
self.D["volumes"] = number_of_volumes
[docs] def get_id(self):
"""Convenient method to get the unique identifier (primary key)
:return: None if not found
"""
if "id" in self.D:
return self.D["id"]
else:
return None
[docs] def get_3d_info(self):
"""Return True if the map is three dimensional,
False if not and None if not info was found"""
if "is_3d" in self.D:
return self.D["is_3d"]
else:
return None
[docs] def get_number_of_points(self):
"""Get the number of points of the vector map
:return: None if not found"""
if "points" in self.D:
return self.D["points"]
else:
return None
[docs] def get_number_of_lines(self):
"""Get the number of lines of the vector map
:return: None if not found"""
if "lines" in self.D:
return self.D["lines"]
else:
return None
[docs] def get_number_of_boundaries(self):
"""Get the number of boundaries of the vector map
:return: None if not found"""
if "boundaries" in self.D:
return self.D["boundaries"]
else:
return None
[docs] def get_number_of_centroids(self):
"""Get the number of centroids of the vector map
:return: None if not found"""
if "centroids" in self.D:
return self.D["centroids"]
else:
return None
[docs] def get_number_of_faces(self):
"""Get the number of faces of the vector map
:return: None if not found"""
if "faces" in self.D:
return self.D["faces"]
else:
return None
[docs] def get_number_of_kernels(self):
"""Get the number of kernels of the vector map
:return: None if not found"""
if "kernels" in self.D:
return self.D["kernels"]
else:
return None
[docs] def get_number_of_primitives(self):
"""Get the number of primitives of the vector map
:return: None if not found"""
if "primitives" in self.D:
return self.D["primitives"]
else:
return None
[docs] def get_number_of_nodes(self):
"""Get the number of nodes of the vector map
:return: None if not found"""
if "nodes" in self.D:
return self.D["nodes"]
else:
return None
[docs] def get_number_of_areas(self):
"""Get the number of areas of the vector map
:return: None if not found"""
if "areas" in self.D:
return self.D["areas"]
else:
return None
[docs] def get_number_of_islands(self):
"""Get the number of islands of the vector map
:return: None if not found"""
if "islands" in self.D:
return self.D["islands"]
else:
return None
[docs] def get_number_of_holes(self):
"""Get the number of holes of the vector map
:return: None if not found"""
if "holes" in self.D:
return self.D["holes"]
else:
return None
[docs] def get_number_of_volumes(self):
"""Get the number of volumes of the vector map
:return: None if not found"""
if "volumes" in self.D:
return self.D["volumes"]
else:
return None
# Set the properties
id = property(fget=get_id, fset=set_id)
is_3d = property(fget=get_3d_info, fset=set_3d_info)
number_of_points = property(fget=get_number_of_points, fset=set_number_of_points)
number_of_lines = property(fget=get_number_of_lines, fset=set_number_of_lines)
number_of_boundaries = property(
fget=get_number_of_boundaries, fset=set_number_of_boundaries
)
number_of_centroids = property(
fget=get_number_of_centroids, fset=set_number_of_centroids
)
number_of_faces = property(fget=get_number_of_faces, fset=set_number_of_faces)
number_of_kernels = property(fget=get_number_of_kernels, fset=set_number_of_kernels)
number_of_primitives = property(
fget=get_number_of_primitives, fset=set_number_of_primitives
)
number_of_nodes = property(fget=get_number_of_nodes, fset=set_number_of_nodes)
number_of_areas = property(fget=get_number_of_areas, fset=set_number_of_areas)
number_of_islands = property(fget=get_number_of_islands, fset=set_number_of_islands)
number_of_holes = property(fget=get_number_of_holes, fset=set_number_of_holes)
number_of_volumes = property(fget=get_number_of_volumes, fset=set_number_of_volumes)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
print(
" +-------------------- Metadata information ----------------------------------+" # noqa: E501
)
print(" | Is map 3d .................. " + str(self.get_3d_info()))
print(" | Number of points ........... " + str(self.get_number_of_points()))
print(" | Number of lines ............ " + str(self.get_number_of_lines()))
print(" | Number of boundaries ....... " + str(self.get_number_of_boundaries()))
print(" | Number of centroids ........ " + str(self.get_number_of_centroids()))
print(" | Number of faces ............ " + str(self.get_number_of_faces()))
print(" | Number of kernels .......... " + str(self.get_number_of_kernels()))
print(" | Number of primitives ....... " + str(self.get_number_of_primitives()))
print(" | Number of nodes ............ " + str(self.get_number_of_nodes()))
print(" | Number of areas ............ " + str(self.get_number_of_areas()))
print(" | Number of islands .......... " + str(self.get_number_of_islands()))
print(" | Number of holes ............ " + str(self.get_number_of_holes()))
print(" | Number of volumes .......... " + str(self.get_number_of_volumes()))
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
print("is_3d=" + str(self.get_3d_info()))
print("points=" + str(self.get_number_of_points()))
print("lines=" + str(self.get_number_of_lines()))
print("boundaries=" + str(self.get_number_of_boundaries()))
print("centroids=" + str(self.get_number_of_centroids()))
print("faces=" + str(self.get_number_of_faces()))
print("kernels=" + str(self.get_number_of_kernels()))
print("primitives=" + str(self.get_number_of_primitives()))
print("nodes=" + str(self.get_number_of_nodes()))
print("areas=" + str(self.get_number_of_areas()))
print("islands=" + str(self.get_number_of_islands()))
print("holes=" + str(self.get_number_of_holes()))
print("volumes=" + str(self.get_number_of_volumes()))
###############################################################################
[docs]class STDSMetadataBase(SQLDatabaseInterface):
"""This is the space time dataset metadata base class for
strds, stvds and str3ds datasets
setting/getting the id, the title and the description
Usage:
.. code-block:: python
>>> init()
>>> meta = STDSMetadataBase(
... ident="soils@PERMANENT", title="Soils", description="Soils 1950 - 2010"
... )
>>> meta.id
'soils@PERMANENT'
>>> meta.title
'Soils'
>>> meta.description
'Soils 1950 - 2010'
>>> meta.number_of_maps
>>> meta.print_info()
| Number of registered maps:.. None
|
| Title:
| Soils
| Description:
| Soils 1950 - 2010
| Command history:
>>> meta.print_shell_info()
number_of_maps=None
"""
def __init__(
self, table=None, ident=None, title=None, description=None, command=None
):
SQLDatabaseInterface.__init__(self, table, ident)
self.set_id(ident)
self.set_title(title)
self.set_description(description)
self.set_command(command)
# No setter for this
self.D["number_of_maps"] = None
[docs] def set_id(self, ident):
"""Convenient method to set the unique identifier (primary key)"""
self.ident = ident
self.D["id"] = ident
[docs] def set_description(self, description):
"""Set the number of cols"""
self.D["description"] = description
[docs] def get_id(self):
"""Convenient method to get the unique identifier (primary key)
:return: None if not found
"""
if "id" in self.D:
return self.D["id"]
else:
return None
[docs] def get_title(self):
"""Get the title
:return: None if not found"""
if "title" in self.D:
return self.D["title"]
else:
return None
[docs] def get_description(self):
"""Get description
:return: None if not found"""
if "description" in self.D:
return self.D["description"]
else:
return None
[docs] def get_command(self):
"""Get command
:return: None if not found"""
if "command" in self.D:
return self.D["command"]
else:
return None
[docs] def get_number_of_maps(self):
"""Get the number of registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "number_of_maps" in self.D:
return self.D["number_of_maps"]
else:
return None
id = property(fget=get_id, fset=set_id)
title = property(fget=get_title, fset=set_title)
description = property(fget=get_description, fset=set_description)
number_of_maps = property(fget=get_number_of_maps)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
self._print_info_body(shell=False)
self._print_info_tail(shell=False)
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_body(shell=True)
self._print_info_tail(shell=True)
def _print_info_head(self, shell=False):
"""Print information about this class (head part).
No header printed in shell style mode.
:param bool shell: True for human readable style otherwise shell style
"""
if not shell:
print(
" +-------------------- Metadata information ----------------------------------+" # noqa: E501
)
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
def _print_info_tail(self, shell=False):
"""Print information about this class (tail part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("number_of_maps=" + str(self.get_number_of_maps()))
else:
print(" | Number of registered maps:.. " + str(self.get_number_of_maps()))
print(" |")
print(" | Title:")
print(" | " + str(self.get_title()))
print(" | Description:")
print(" | " + str(self.get_description()))
print(" | Command history:")
command = self.get_command()
if command:
for token in command.split("\n"):
print(" | " + str(token))
[docs] def print_history(self):
"""Print history information about this class in human readable
shell style
"""
# 0123456789012345678901234567890
print("# Title:")
print("# " + str(self.get_title()))
print("# Description:")
print("# " + str(self.get_description()))
print("# Command history:")
command = self.get_command()
if command:
tokens = command.split("\n")
print_list = []
for token in tokens:
token = str(token).rstrip().lstrip()
if len(token) > 1:
print_list.append(token)
count = 0
for token in print_list:
count += 1
if len(token) > 1:
if token[0] == "#":
print(token)
elif count < len(print_list):
print(token + " \\")
else:
print(token)
###############################################################################
[docs]class STDSRasterMetadataBase(STDSMetadataBase):
"""This is the space time dataset metadata base
class for strds and str3ds datasets
Most of the metadata values are set by SQL scripts in the database when
new maps are added. Therefor only some set- an many
get-functions are available.
Usage:
.. code-block:: python
>>> init()
>>> meta = STDSRasterMetadataBase(
... ident="soils@PERMANENT", title="Soils", description="Soils 1950 - 2010"
... )
>>> meta.id
'soils@PERMANENT'
>>> meta.title
'Soils'
>>> meta.description
'Soils 1950 - 2010'
>>> meta.number_of_maps
>>> meta.min_max
>>> meta.max_max
>>> meta.min_min
>>> meta.max_min
>>> meta.nsres_min
>>> meta.nsres_max
>>> meta.ewres_min
>>> meta.ewres_max
>>> meta.print_info()
| North-South resolution min:. None
| North-South resolution max:. None
| East-west resolution min:... None
| East-west resolution max:... None
| Minimum value min:.......... None
| Minimum value max:.......... None
| Maximum value min:.......... None
| Maximum value max:.......... None
| Aggregation type:........... None
| Number of registered maps:.. None
|
| Title:
| Soils
| Description:
| Soils 1950 - 2010
| Command history:
>>> meta.print_shell_info()
nsres_min=None
nsres_max=None
ewres_min=None
ewres_max=None
min_min=None
min_max=None
max_min=None
max_max=None
aggregation_type=None
number_of_maps=None
"""
def __init__(
self,
table=None,
ident=None,
title=None,
description=None,
aggregation_type=None,
):
STDSMetadataBase.__init__(self, table, ident, title, description)
# Initialize the dict to select all values from the db
self.D["min_max"] = None
self.D["max_max"] = None
self.D["min_min"] = None
self.D["max_min"] = None
self.D["nsres_min"] = None
self.D["nsres_max"] = None
self.D["ewres_min"] = None
self.D["ewres_max"] = None
self.D["aggregation_type"] = aggregation_type
[docs] def set_aggregation_type(self, aggregation_type):
"""Set the aggregation type of the dataset (mean, min, max, ...)"""
self.D["aggregation_type"] = aggregation_type
[docs] def get_aggregation_type(self):
"""Get the aggregation type of the dataset (mean, min, max, ...)
:return: None if not found
"""
if "aggregation_type" in self.D:
return self.D["aggregation_type"]
else:
return None
[docs] def get_max_min(self):
"""Get the minimal maximum of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "max_min" in self.D:
return self.D["max_min"]
else:
return None
[docs] def get_min_min(self):
"""Get the minimal minimum of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "min_min" in self.D:
return self.D["min_min"]
else:
return None
[docs] def get_max_max(self):
"""Get the maximal maximum of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "max_max" in self.D:
return self.D["max_max"]
else:
return None
[docs] def get_min_max(self):
"""Get the maximal minimum of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "min_max" in self.D:
return self.D["min_max"]
else:
return None
[docs] def get_nsres_min(self):
"""Get the minimal north-south resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "nsres_min" in self.D:
return self.D["nsres_min"]
else:
return None
[docs] def get_nsres_max(self):
"""Get the maximal north-south resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "nsres_max" in self.D:
return self.D["nsres_max"]
else:
return None
[docs] def get_ewres_min(self):
"""Get the minimal east-west resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "ewres_min" in self.D:
return self.D["ewres_min"]
else:
return None
[docs] def get_ewres_max(self):
"""Get the maximal east-west resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "ewres_max" in self.D:
return self.D["ewres_max"]
else:
return None
nsres_min = property(fget=get_nsres_min)
nsres_max = property(fget=get_nsres_max)
ewres_min = property(fget=get_ewres_min)
ewres_max = property(fget=get_ewres_max)
min_min = property(fget=get_min_min)
min_max = property(fget=get_min_max)
max_min = property(fget=get_max_min)
max_max = property(fget=get_max_max)
aggregation_type = property(fset=set_aggregation_type, fget=get_aggregation_type)
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("nsres_min=" + str(self.get_nsres_min()))
print("nsres_max=" + str(self.get_nsres_max()))
print("ewres_min=" + str(self.get_ewres_min()))
print("ewres_max=" + str(self.get_ewres_max()))
print("min_min=" + str(self.get_min_min()))
print("min_max=" + str(self.get_min_max()))
print("max_min=" + str(self.get_max_min()))
print("max_max=" + str(self.get_max_max()))
print("aggregation_type=" + str(self.get_aggregation_type()))
else:
print(" | North-South resolution min:. " + str(self.get_nsres_min()))
print(" | North-South resolution max:. " + str(self.get_nsres_max()))
print(" | East-west resolution min:... " + str(self.get_ewres_min()))
print(" | East-west resolution max:... " + str(self.get_ewres_max()))
print(" | Minimum value min:.......... " + str(self.get_min_min()))
print(" | Minimum value max:.......... " + str(self.get_min_max()))
print(" | Maximum value min:.......... " + str(self.get_max_min()))
print(" | Maximum value max:.......... " + str(self.get_max_max()))
print(" | Aggregation type:........... " + str(self.get_aggregation_type()))
###############################################################################
[docs]class STRDSMetadata(STDSRasterMetadataBase):
"""This is the raster metadata class
This class is the interface to the strds_metadata table in the
temporal database that stores the metadata of all registered
space time raster datasets
Most of the metadata values are set by SQL scripts in the database when
new raster maps are added. Therefor only some set- an many
get-functions are available.
Usage:
.. code-block:: python
>>> init()
>>> meta = STRDSMetadata(
... ident="soils@PERMANENT", title="Soils", description="Soils 1950 - 2010"
... )
>>> meta.id
'soils@PERMANENT'
>>> meta.title
'Soils'
>>> meta.description
'Soils 1950 - 2010'
>>> meta.number_of_maps
>>> meta.min_max
>>> meta.max_max
>>> meta.min_min
>>> meta.max_min
>>> meta.nsres_min
>>> meta.nsres_max
>>> meta.ewres_min
>>> meta.ewres_max
>>> meta.raster_register
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| Raster register table:...... None
| North-South resolution min:. None
| North-South resolution max:. None
| East-west resolution min:... None
| East-west resolution max:... None
| Minimum value min:.......... None
| Minimum value max:.......... None
| Maximum value min:.......... None
| Maximum value max:.......... None
| Aggregation type:........... None
| Number of semantic labels:.. None
| Semantic labels:............ None
| Number of registered maps:.. None
|
| Title:
| Soils
| Description:
| Soils 1950 - 2010
| Command history:
>>> meta.print_shell_info()
raster_register=None
nsres_min=None
nsres_max=None
ewres_min=None
ewres_max=None
min_min=None
min_max=None
max_min=None
max_max=None
aggregation_type=None
number_of_semantic_labels=None
semantic_labels=None
number_of_maps=None
"""
def __init__(self, ident=None, raster_register=None, title=None, description=None):
STDSRasterMetadataBase.__init__(
self, "strds_metadata", ident, title, description
)
if get_tgis_db_version_from_metadata() > 2:
self.D["number_of_semantic_labels"] = None
self.set_raster_register(raster_register)
[docs] def set_raster_register(self, raster_register):
"""Set the raster map register table name"""
self.D["raster_register"] = raster_register
[docs] def get_raster_register(self):
"""Get the raster map register table name
:return: None if not found"""
if "raster_register" in self.D:
return self.D["raster_register"]
else:
return None
[docs] def get_number_of_semantic_labels(self):
"""Get the number of registered semantic labels
:return: None if not found
"""
if "number_of_semantic_labels" in self.D:
return self.D["number_of_semantic_labels"]
else:
return None
[docs] def get_semantic_labels(self):
"""Get the distinct semantic labels of registered maps
The distinct semantic labels are not stored in the metadata table
and fetched on-the-fly
:return: None if not found
"""
if get_tgis_db_version_from_metadata() <= 2:
# band names supported from TGIS DB version 3
return None
sql = "SELECT distinct semantic_label FROM %s WHERE %s.id " % (
"raster_metadata",
"raster_metadata",
)
sql += "IN (SELECT id FROM %s)" % (str(self.get_raster_register()))
dbif = SQLDatabaseInterfaceConnection()
dbif.connect()
dbif.execute(sql, mapset=self.mapset)
rows = dbif.fetchall(mapset=self.mapset)
dbif.close()
if rows:
string = ""
count = 0
for row in rows:
if row["semantic_label"]:
if count == 0:
string += row["semantic_label"]
else:
string += ",%s" % row["semantic_label"]
count += 1
if count > 0:
return string
else:
return None
else:
return None
raster_register = property(fget=get_raster_register, fset=set_raster_register)
number_of_semantic_labels = property(fget=get_number_of_semantic_labels)
semantic_labels = property(fget=get_semantic_labels)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
self._print_info_head(shell=False)
super().print_info()
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_head(shell=True)
super().print_shell_info()
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("raster_register=" + str(self.get_raster_register()))
else:
print(" | Raster register table:...... " + str(self.get_raster_register()))
super()._print_info_body(shell)
if shell:
print(
"number_of_semantic_labels=" + str(self.get_number_of_semantic_labels())
)
print("semantic_labels=" + str(self.get_semantic_labels()))
else:
print(
" | Number of semantic labels:.. "
+ str(self.get_number_of_semantic_labels())
)
print(" | Semantic labels:............ " + str(self.get_semantic_labels()))
###############################################################################
[docs]class STR3DSMetadata(STDSRasterMetadataBase):
"""This is the space time 3D raster metadata class
This class is the interface to the str3ds_metadata table in the
temporal database that stores the metadata of all registered
space time 3D raster datasets
Most of the metadata values are set by SQL scripts in the database when
new 3D raster maps are added. Therefor only some set- an many
get-functions are available.
Usage:
.. code-block:: python
>>> init()
>>> meta = STR3DSMetadata(
... ident="soils@PERMANENT", title="Soils", description="Soils 1950 - 2010"
... )
>>> meta.id
'soils@PERMANENT'
>>> meta.title
'Soils'
>>> meta.description
'Soils 1950 - 2010'
>>> meta.number_of_maps
>>> meta.min_max
>>> meta.max_max
>>> meta.min_min
>>> meta.max_min
>>> meta.nsres_min
>>> meta.nsres_max
>>> meta.ewres_min
>>> meta.ewres_max
>>> meta.tbres_min
>>> meta.tbres_max
>>> meta.raster3d_register
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| 3D raster register table:... None
| Top-bottom resolution min:.. None
| Top-bottom resolution max:.. None
| North-South resolution min:. None
| North-South resolution max:. None
| East-west resolution min:... None
| East-west resolution max:... None
| Minimum value min:.......... None
| Minimum value max:.......... None
| Maximum value min:.......... None
| Maximum value max:.......... None
| Aggregation type:........... None
| Number of registered maps:.. None
|
| Title:
| Soils
| Description:
| Soils 1950 - 2010
| Command history:
>>> meta.print_shell_info()
raster3d_register=None
tbres_min=None
tbres_max=None
nsres_min=None
nsres_max=None
ewres_min=None
ewres_max=None
min_min=None
min_max=None
max_min=None
max_max=None
aggregation_type=None
number_of_maps=None
"""
def __init__(
self, ident=None, raster3d_register=None, title=None, description=None
):
STDSRasterMetadataBase.__init__(
self, "str3ds_metadata", ident, title, description
)
self.set_raster3d_register(raster3d_register)
self.D["tbres_min"] = None
self.D["tbres_max"] = None
[docs] def set_raster3d_register(self, raster3d_register):
"""Set the raster map register table name"""
self.D["raster3d_register"] = raster3d_register
[docs] def get_raster3d_register(self):
"""Get the raster3d map register table name
:return: None if not found"""
if "raster3d_register" in self.D:
return self.D["raster3d_register"]
else:
return None
[docs] def get_tbres_min(self):
"""Get the minimal top-bottom resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "tbres_min" in self.D:
return self.D["tbres_min"]
else:
return None
[docs] def get_tbres_max(self):
"""Get the maximal top-bottom resolution of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "tbres_max" in self.D:
return self.D["tbres_max"]
else:
return None
raster3d_register = property(fget=get_raster3d_register, fset=set_raster3d_register)
tbres_min = property(fget=get_tbres_min)
tbres_max = property(fget=get_tbres_max)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
self._print_info_head(shell=False)
super().print_info()
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_head(shell=True)
super().print_shell_info()
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("raster3d_register=" + str(self.get_raster3d_register()))
print("tbres_min=" + str(self.get_tbres_min()))
print("tbres_max=" + str(self.get_tbres_max()))
else:
print(
" | 3D raster register table:... " + str(self.get_raster3d_register())
)
print(" | Top-bottom resolution min:.. " + str(self.get_ewres_min()))
print(" | Top-bottom resolution max:.. " + str(self.get_ewres_max()))
super()._print_info_body(shell)
###############################################################################
[docs]class STVDSMetadata(STDSMetadataBase):
"""This is the space time vector dataset metadata class
This class is the interface to the stvds_metadata table in the
temporal database that stores the metadata of all registered
space time vector datasets
Most of the metadata values are set by SQL scripts in the database when
new vector maps are added. Therefor only some set- an many get-functions
are available.
Usage:
.. code-block:: python
>>> init()
>>> meta = STVDSMetadata(
... ident="lidars@PERMANENT",
... title="LIDARS",
... description="LIDARS 2008 - 2010",
... )
>>> meta.id
'lidars@PERMANENT'
>>> meta.title
'LIDARS'
>>> meta.description
'LIDARS 2008 - 2010'
>>> meta.number_of_maps
>>> meta.number_of_points
>>> meta.number_of_lines
>>> meta.number_of_boundaries
>>> meta.number_of_centroids
>>> meta.number_of_faces
>>> meta.number_of_kernels
>>> meta.number_of_primitives
>>> meta.number_of_nodes
>>> meta.number_of_areas
>>> meta.number_of_islands
>>> meta.number_of_holes
>>> meta.number_of_volumes
>>> meta.print_info()
+-------------------- Metadata information ----------------------------------+
| Vector register table:...... None
| Number of points ........... None
| Number of lines ............ None
| Number of boundaries ....... None
| Number of centroids ........ None
| Number of faces ............ None
| Number of kernels .......... None
| Number of primitives ....... None
| Number of nodes ............ None
| Number of areas ............ None
| Number of islands .......... None
| Number of holes ............ None
| Number of volumes .......... None
| Number of registered maps:.. None
|
| Title:
| LIDARS
| Description:
| LIDARS 2008 - 2010
| Command history:
>>> meta.print_shell_info()
vector_register=None
points=None
lines=None
boundaries=None
centroids=None
faces=None
kernels=None
primitives=None
nodes=None
areas=None
islands=None
holes=None
volumes=None
number_of_maps=None
"""
def __init__(self, ident=None, vector_register=None, title=None, description=None):
STDSMetadataBase.__init__(self, "stvds_metadata", ident, title, description)
self.set_vector_register(vector_register)
self.D["points"] = None
self.D["lines"] = None
self.D["boundaries"] = None
self.D["centroids"] = None
self.D["faces"] = None
self.D["kernels"] = None
self.D["primitives"] = None
self.D["nodes"] = None
self.D["areas"] = None
self.D["islands"] = None
self.D["holes"] = None
self.D["volumes"] = None
[docs] def set_vector_register(self, vector_register):
"""Set the vector map register table name"""
self.D["vector_register"] = vector_register
[docs] def get_vector_register(self):
"""Get the vector map register table name
:return: None if not found"""
if "vector_register" in self.D:
return self.D["vector_register"]
else:
return None
[docs] def get_number_of_points(self):
"""Get the number of points of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "points" in self.D:
return self.D["points"]
else:
return None
[docs] def get_number_of_lines(self):
"""Get the number of lines of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "lines" in self.D:
return self.D["lines"]
else:
return None
[docs] def get_number_of_boundaries(self):
"""Get the number of boundaries of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "boundaries" in self.D:
return self.D["boundaries"]
else:
return None
[docs] def get_number_of_centroids(self):
"""Get the number of centroids of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "centroids" in self.D:
return self.D["centroids"]
else:
return None
[docs] def get_number_of_faces(self):
"""Get the number of faces of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "faces" in self.D:
return self.D["faces"]
else:
return None
[docs] def get_number_of_kernels(self):
"""Get the number of kernels of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "kernels" in self.D:
return self.D["kernels"]
else:
return None
[docs] def get_number_of_primitives(self):
"""Get the number of primitives of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "primitives" in self.D:
return self.D["primitives"]
else:
return None
[docs] def get_number_of_nodes(self):
"""Get the number of nodes of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "nodes" in self.D:
return self.D["nodes"]
else:
return None
[docs] def get_number_of_areas(self):
"""Get the number of areas of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "areas" in self.D:
return self.D["areas"]
else:
return None
[docs] def get_number_of_islands(self):
"""Get the number of islands of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "islands" in self.D:
return self.D["islands"]
else:
return None
[docs] def get_number_of_holes(self):
"""Get the number of holes of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "holes" in self.D:
return self.D["holes"]
else:
return None
[docs] def get_number_of_volumes(self):
"""Get the number of volumes of all registered maps,
this value is set in the database
automatically via SQL, so no setter exists
:return: None if not found"""
if "volumes" in self.D:
return self.D["volumes"]
else:
return None
# Set the properties
vector_register = property(fget=get_vector_register, fset=set_vector_register)
number_of_points = property(fget=get_number_of_points)
number_of_lines = property(fget=get_number_of_lines)
number_of_boundaries = property(fget=get_number_of_boundaries)
number_of_centroids = property(fget=get_number_of_centroids)
number_of_faces = property(fget=get_number_of_faces)
number_of_kernels = property(fget=get_number_of_kernels)
number_of_primitives = property(fget=get_number_of_primitives)
number_of_nodes = property(fget=get_number_of_nodes)
number_of_areas = property(fget=get_number_of_areas)
number_of_islands = property(fget=get_number_of_islands)
number_of_holes = property(fget=get_number_of_holes)
number_of_volumes = property(fget=get_number_of_volumes)
[docs] def print_info(self):
"""Print information about this class in human readable style"""
self._print_info_head(shell=False)
super().print_info()
[docs] def print_shell_info(self):
"""Print information about this class in shell style"""
self._print_info_head(shell=True)
super().print_shell_info()
def _print_info_body(self, shell=False):
"""Print information about this class (body part).
:param bool shell: True for human readable style otherwise shell style
"""
if shell:
print("vector_register=" + str(self.get_vector_register()))
print("points=" + str(self.get_number_of_points()))
print("lines=" + str(self.get_number_of_lines()))
print("boundaries=" + str(self.get_number_of_boundaries()))
print("centroids=" + str(self.get_number_of_centroids()))
print("faces=" + str(self.get_number_of_faces()))
print("kernels=" + str(self.get_number_of_kernels()))
print("primitives=" + str(self.get_number_of_primitives()))
print("nodes=" + str(self.get_number_of_nodes()))
print("areas=" + str(self.get_number_of_areas()))
print("islands=" + str(self.get_number_of_islands()))
print("holes=" + str(self.get_number_of_holes()))
print("volumes=" + str(self.get_number_of_volumes()))
else:
print(" | Vector register table:...... " + str(self.get_vector_register()))
print(" | Number of points ........... " + str(self.number_of_points))
print(" | Number of lines ............ " + str(self.number_of_lines))
print(" | Number of boundaries ....... " + str(self.number_of_boundaries))
print(" | Number of centroids ........ " + str(self.number_of_centroids))
print(" | Number of faces ............ " + str(self.number_of_faces))
print(" | Number of kernels .......... " + str(self.number_of_kernels))
print(" | Number of primitives ....... " + str(self.number_of_primitives))
print(" | Number of nodes ............ " + str(self.number_of_nodes))
print(" | Number of areas ............ " + str(self.number_of_areas))
print(" | Number of islands .......... " + str(self.number_of_islands))
print(" | Number of holes ............ " + str(self.number_of_holes))
print(" | Number of volumes .......... " + str(self.number_of_volumes))
###############################################################################
if __name__ == "__main__":
import doctest
doctest.testmod()