Temporal extension for GRASS GIS 7: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
Line 1: Line 1:
= This page has been moved to [http://trac.osgeo.org/grass/wiki/Grass7/G3dLib New trac based wiki] =
= Adding the time dimension to GRASS GIS raster, vector and voxel data types for field based applications using a time stamp approach =
= Adding the time dimension to GRASS GIS raster, vector and voxel data types for field based applications using a time stamp approach =



Revision as of 05:25, 30 June 2011

This page has been moved to New trac based wiki

Adding the time dimension to GRASS GIS raster, vector and voxel data types for field based applications using a time stamp approach

Time stamps in SQL tables

Raster, vector and voxel maps in GRASS GIS will have temporal attributes by default:

  • Creation time
  • Modification time
  • An absolute or relative time stamp
    • absolute:
      • start date
      • end date
    • relative
      • a relative time stamp

To store the time stamps and to enable temporal queries a SQL DBMS, capable of temporal datatypes, will be used. As default a Sqlite3 database will be used to store map specific informations in database files. Additionally Postgresql should be supported too.

In case of SQlite a GRASS SQL database file is created for each mapset in a location. Views are used to access the raster, vector or raster3d map entries of different mapset databases at once.

The time stamps as well as attributes and map specific meta data of all created maps in GRASS GIS will be stored in SQL tables in the GRASS SQL databases. The entries in the SQL tables will be created or updated automatically while creation, modification or removement of maps. For each GRASS datatype (raster, raster3d and vector) four SQL tables will be added to GRASS SQL databases at mapset creation:

  • The base table has an equal layout for raster, vector and raster3d maps. This table stores:
    • The id which is of type string. As id the raster, raster3d or vector map names should be used, because these must be unique in a mapset and can be used as primary key
    • The temporal type specifies if the time stamp is relative or absolute
    • The creation date is of type datetime
    • The modification date is of type datetime


  • The absolute time table to store the absolute time stamps for each map. The layout is equal for raster, raster3d and vector maps. This table stores:
    • The id as primary key
    • The start time stamp is of type datetime
    • The end time stampis of type datetime
    • The time zoneis of type integer
  • The relative time table to store the relative time stamps for each map. The layout is equal for raster, raster3d and vector maps. This table stores:
    • The id as primary key
    • The relative time stamp is of type datetime

A map can only have an absolute or relative time stamp.

  • The metadata table is specific for each grass data type. It stores the spatial and temporal extent as well as map specific attributes (creator, range, cell type, features, registered spacetime datasets, ...)

In long term these kind of SQL tables may replace the text file based storage of map metadata, database links and time stamps.

New temporal GRASS datatypes

To manage temporal raster, raster3d and vector objects, three new datatypes will be added to GRASS GIS: Spacetime raster, vector and voxel datasets. Additionally thematic spatio-temporal datasets will be added to handle thematic related spacetime datasets of different type.

Maps can be registered and unregistered in spacetime datasets. Maps are registered by id (name). All available spatio-temporal data is availabe in SQL tables. Cross referencing must be implemented (maps storing the spacetime datasets in which they are registered and vis versa) Spacetime datasets have absolute or relative time type, these types cant be mixed up but converted.

A temporal resolution can be set for each spacetime dataset (1s - decades) defining the granularity of the datasets (like the resolution settings for raster/voxel maps).

The new spacetime datasets are handled like raster or vector maps. They are registered in mapset specific tables with creation and modification date and absolute/relative time stamp.

Raster spacetime datasets

The raster spacetime datasets store raster map specific attributes and metadata:

  • The spatial extent (e, w, n, s)
  • The temporal extent (start and end time or duration)
  • The data range (min and max)
  • The raster datatype (CELL, FCELL, DCELL, MIXED)

All these table entries are computed or updated via trigger from all registered raster maps when maps are added or removed.

Table layout:

The base spacetime tables are equal to the raster/vector and raster3d one.

  • Spacetime metadata raster table
    • id the name of the spacetime raster table
    • temporal type relative or absolute
    • Raster register table name the name of the table storing the ids (names) of all registered raster maps
    • spatial extent (4 columns: e, w, s,n) of all registered raster maps
    • temporal extent of all registered raster maps
    • data range of all registered raster maps
    • ...

Raster3d spacetime datasets

The raster3d spacetime tables are equal to the raster spacetime tablse, except that the spatial extent is computed based on east, west, south, north, top and bottom.

Vector spacetime datasets

New temporal modules

New temporal modules must be implemented for registration and un-registration and extraction of raster, voxel and vector maps in süacetime datasets:

  • tr.register, tr.unregister, tr.extract
  • t3.register, t3.unregister, t3.extract
  • tv.register, tv.unregister, tv.extract
  • Spacetime modules for resampling, interpolation, time stamp snapping and extraction must be implemented.

Basic Spatio-temporal algorithms

Extraction is based on where SQL statements or specific time stamps. Interpolated maps are created or existing map ids will be returned. The extraction tool will provide different types of extraction algorithms Nearest neighbour Interpolation (only raster and voxel maps) (linear, quadratic, ...) or first valid time stamp.


tr(3).resample -> will create a new spacetime raster (voxel) datasets with a different spatio-temporal resolution, several methods will be available:

  • averaging (mean, median, geometric mean, ...)
  • summation, multiplication
  • max and min functions

Implementation details

Access to spacetime tables must be implemented as GRASS GIS library:

  • Registration and unregistration of maps into spacetime datasets
  • Update/creation of base table entries, so the SQL tables can be modified in low level GRASS GIS library functions (creation, deletion of maps)
  • Selection of maps from spacetime datasets