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 __future__ import print_function
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_datatype(self, datatype): """Set the datatype""" self.D["datatype"] = datatype
[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) 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 >>> 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, semantic_label=None, ): RasterMetadataBase.__init__( self, "raster_metadata", ident, datatype, cols, rows, number_of_cells, nsres, ewres, min, max, ) if get_tgis_db_version_from_metadata() > 2: self.set_semantic_label(semantic_label)
[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) 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) # semantic label section (raster specific only) if shell: print("semantic_label=" + str(self.get_semantic_label())) else: 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) 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()))
###############################################################################
[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) 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()) )
###############################################################################
[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_title(self, title): """Set the title""" self.D["title"] = title
[docs] def set_description(self, description): """Set the number of cols""" self.D["description"] = description
[docs] def set_command(self, command): """Set the number of cols""" self.D["command"] = command
[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_head(shell=False) 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_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))
[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 lables 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) 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 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
[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) 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) 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()