Geografic-Information-System/python/grass/temporal/metadata.py

"""
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 __future__ import print_function
from .base import SQLDatabaseInterface

###############################################################################


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)

    def set_id(self, ident):
        """Convenient method to set the unique identifier (primary key)"""
        self.ident = ident
        self.D["id"] = ident

    def set_datatype(self, datatype):
        """Set the datatype"""
        self.D["datatype"] = datatype

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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)

    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()))


###############################################################################


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
        >>> 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,
        ident=None,
        datatype=None,
        cols=None,
        rows=None,
        number_of_cells=None,
        nsres=None,
        ewres=None,
        min=None,
        max=None,
        band_reference=None,
    ):

        RasterMetadataBase.__init__(
            self,
            "raster_metadata",
            ident,
            datatype,
            cols,
            rows,
            number_of_cells,
            nsres,
            ewres,
            min,
            max,
        )

        self.set_band_reference(band_reference)

    def set_band_reference(self, band_reference):
        """Set the band reference identifier"""
        self.D["band_reference"] = band_reference

    def get_band_reference(self):
        """Get the band reference identifier
        :return: None if not found"""
        if "band_reference" in self.D:
            return self.D["band_reference"]
        else:
            return None

    band_reference = property(fget=get_band_reference, fset=set_band_reference)

    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
        """
        super()._print_info_body(shell)
        # band reference section (raster specific only)
        if shell:
            print("band_reference=" + str(self.get_band_reference()))
        else:
            print(" | Band reference:............. " + str(self.get_band_reference()))


###############################################################################


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)

    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

    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

    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

    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)

    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
        """
        super()._print_info_body(shell)
        if shell:
            print("depths=" + str(self.get_depths()))
            print("tbres=" + str(self.get_tbres()))
        else:
            print(" | Number of depths:........... " + str(self.get_depths()))
            print(" | Top-Bottom resolution:...... " + str(self.get_tbres()))


###############################################################################


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)

    def set_id(self, ident):
        """Convenient method to set the unique identifier (primary key)"""
        self.ident = ident
        self.D["id"] = ident

    def set_3d_info(self, is_3d):
        """Set True if the vector map is three dimensional"""
        self.D["is_3d"] = is_3d

    def set_number_of_points(self, number_of_points):
        """Set the number of points of the vector map"""
        self.D["points"] = number_of_points

    def set_number_of_lines(self, number_of_lines):
        """Set the number of lines of the vector map"""
        self.D["lines"] = number_of_lines

    def set_number_of_boundaries(self, number_of_boundaries):
        """Set the number of boundaries of the vector map"""
        self.D["boundaries"] = number_of_boundaries

    def set_number_of_centroids(self, number_of_centroids):
        """Set the number of centroids of the vector map"""
        self.D["centroids"] = number_of_centroids

    def set_number_of_faces(self, number_of_faces):
        """Set the number of faces of the vector map"""
        self.D["faces"] = number_of_faces

    def set_number_of_kernels(self, number_of_kernels):
        """Set the number of kernels of the vector map"""
        self.D["kernels"] = number_of_kernels

    def set_number_of_primitives(self, number_of_primitives):
        """Set the number of primitives of the vector map"""
        self.D["primitives"] = number_of_primitives

    def set_number_of_nodes(self, number_of_nodes):
        """Set the number of nodes of the vector map"""
        self.D["nodes"] = number_of_nodes

    def set_number_of_areas(self, number_of_areas):
        """Set the number of areas of the vector map"""
        self.D["areas"] = number_of_areas

    def set_number_of_islands(self, number_of_islands):
        """Set the number of islands of the vector map"""
        self.D["islands"] = number_of_islands

    def set_number_of_holes(self, number_of_holes):
        """Set the number of holes of the vector map"""
        self.D["holes"] = number_of_holes

    def set_number_of_volumes(self, number_of_volumes):
        """Set the number of volumes of the vector map"""
        self.D["volumes"] = number_of_volumes

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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)

    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("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()))
        else:
            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())
            )


###############################################################################


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

    def set_id(self, ident):
        """Convenient method to set the unique identifier (primary key)"""
        self.ident = ident
        self.D["id"] = ident

    def set_title(self, title):
        """Set the title"""
        self.D["title"] = title

    def set_description(self, description):
        """Set the number of cols"""
        self.D["description"] = description

    def set_command(self, command):
        """Set the number of cols"""
        self.D["command"] = command

    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

    def get_title(self):
        """Get the title
        :return: None if not found"""
        if "title" in self.D:
            return self.D["title"]
        else:
            return None

    def get_description(self):
        """Get description
        :return: None if not found"""
        if "description" in self.D:
            return self.D["description"]
        else:
            return None

    def get_command(self):
        """Get command
        :return: None if not found"""
        if "command" in self.D:
            return self.D["command"]
        else:
            return None

    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)

    def print_info(self):
        """Print information about this class in human readable style"""
        self._print_info_head(shell=False)
        self._print_info_body(shell=False)
        self._print_info_tail(shell=False)

    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)
        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 ----------------------------------+"
            )

    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))

    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)


###############################################################################


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

    def set_aggregation_type(self, aggregation_type):
        """Set the aggregation type of the dataset (mean, min, max, ...)"""
        self.D["aggregation_type"] = aggregation_type

    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

    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

    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

    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

    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

    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

    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

    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

    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()))


###############################################################################


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 registered bands:. 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_bands=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
        )

        self.D["number_of_bands"] = None

        self.set_raster_register(raster_register)

    def set_raster_register(self, raster_register):
        """Set the raster map register table name"""
        self.D["raster_register"] = raster_register

    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

    def get_number_of_bands(self):
        """Get the number of registered bands
        :return: None if not found
        """
        if "number_of_bands" in self.D:
            return self.D["number_of_bands"]
        else:
            return None

    raster_register = property(fget=get_raster_register, fset=set_raster_register)
    number_of_bands = property(fget=get_number_of_bands)

    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_bands=" + str(self.get_number_of_bands()))
        else:
            print(" | Number of registered bands:. " + str(self.get_number_of_bands()))


###############################################################################


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
        nsres_min=None
        nsres_max=None
        ewres_min=None
        ewres_max=None
        min_min=None
        min_max=None
        max_min=None
        max_max=None
        tbres_min=None
        tbres_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

    def set_raster3d_register(self, raster3d_register):
        """Set the raster map register table name"""
        self.D["raster3d_register"] = raster3d_register

    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

    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

    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)

    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)


###############################################################################


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

    def set_vector_register(self, vector_register):
        """Set the vector map register table name"""
        self.D["vector_register"] = vector_register

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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

    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)

    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()