Added module for temporal sampling of raster space time datasets.

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@49057 15284696-431f-4ddb-bdfa-cd5b030d7da7

temporal/tr.sample/Makefile

@@ -0,0 +1,7 @@
+MODULE_TOPDIR = ../../
+
+PGM = tr.sample
+
+include $(MODULE_TOPDIR)/include/Make/Script.make
+
+default: script $(TEST_DST)

temporal/tr.sample/test.tr.sample.sh

@@ -0,0 +1,76 @@
+# This is a test to sample a space time raster dataset with a space time vector dataset
+
+# We need to set a specific region in the
+# @preprocess step of this test. We generate
+# raster with r.mapcalc and create a space time raster datasets
+# The region setting should work for UTM and LL test locations
+g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3
+
+r.mapcalc --o expr="prec_1 = rand(0, 550)"
+r.mapcalc --o expr="prec_2 = rand(0, 450)"
+r.mapcalc --o expr="prec_3 = rand(0, 320)"
+r.mapcalc --o expr="prec_4 = rand(0, 510)"
+r.mapcalc --o expr="prec_5 = rand(0, 300)"
+r.mapcalc --o expr="prec_6 = rand(0, 650)"
+
+v.random --o -z output=pnts1 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts2 n=20 zmin=0 zmax=100 column=height
+
+v.random --o -z output=pnts3 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts4 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts5 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts6 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts7 n=20 zmin=0 zmax=100 column=height
+v.random --o -z output=pnts8 n=20 zmin=0 zmax=100 column=height
+
+n1=`g.tempfile pid=1 -d` 
+n2=`g.tempfile pid=2 -d`
+n3=`g.tempfile pid=3 -d`
+
+cat > $n1 << EOF
+prec_1
+prec_2
+prec_3
+prec_4
+prec_5
+prec_6
+EOF
+
+cat > $n2 << EOF
+pnts1|2001-01-01|2001-03-01
+pnts2|2001-05-01|2001-07-01
+EOF
+
+cat > $n3 << EOF
+pnts3|2001-01-01|2001-01-05
+pnts4|2001-01-05|2001-01-10
+pnts5|2001-01-10|2001-01-15
+pnts6|2001-01-15|2001-01-20
+pnts7|2001-01-20|2001-01-25
+pnts8|2001-01-25|2001-01-30
+EOF
+
+
+t.create --o type=strds temporaltype=absolute output=precip_abs0 title="A test with raster input files" descr="A test with raster input files"
+t.create --o type=stvds temporaltype=absolute output=pnts_abs0 title="A test with vector input files" descr="A test with vector input files"
+t.create --o type=stvds temporaltype=absolute output=pnts_abs1 title="A test with vector input files" descr="A test with vector input files"
+
+tr.register -i input=precip_abs0 file=$n1 start="2001-01-01" increment="1 months"
+tv.register    input=pnts_abs0 file=$n2 start=file end=file
+tv.register    input=pnts_abs1 file=$n3 start=file end=file
+
+# The @test
+tr.sample input=precip_abs0 type=strds sample=precip_abs0 -h
+tr.sample input=precip_abs0 type=stvds sample=pnts_abs0 -h
+tr.sample input=precip_abs0 type=stvds sample=pnts_abs1 -h
+
+# Test with temporal point data
+tr.register    input=precip_abs0 file=$n1 start="2001-01-01" increment="1 months"
+tr.sample input=precip_abs0 type=strds sample=precip_abs0 -hs
+tr.sample input=precip_abs0 type=stvds sample=pnts_abs0 -hs
+tr.sample input=precip_abs0 type=stvds sample=pnts_abs1 -hs
+
+
+t.remove type=rast input=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
+t.remove type=strds input=precip_abs0
+t.remove type=stvds input=pnts_abs0,pnts_abs1

temporal/tr.sample/tr.sample.py

@@ -0,0 +1,91 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+############################################################################
+#
+# MODULE:	t.sample
+# AUTHOR(S):	Soeren Gebbert
+#               
+# PURPOSE:	 Sample the input space time raster dataset with a sample space time dataset and print the result to stdout
+# COPYRIGHT:	(C) 2011 by the GRASS Development Team
+#
+#		This program is free software under the GNU General Public
+#		License (version 2). Read the file COPYING that comes with GRASS
+#		for details.
+#
+#############################################################################
+
+#%module
+#% description:  Sample the input space time dataset with a sample space time dataset and print the result to stdout
+#% keywords: dataset
+#% keywords: spacetime
+#% keywords: raster
+#% keywords: sample
+#%end
+
+#%option
+#% key: input
+#% type: string
+#% description: Name of a space time raster dataset
+#% required: yes
+#% multiple: no
+#%end
+
+#%option
+#% key: sample
+#% type: string
+#% description: Name of the sample space time raster dataset 
+#% required: yes
+#% multiple: no
+#%end
+
+#%option
+#% key: type
+#% type: string
+#% description: Type of the sample space time dataset, default is space time raster dataset (strds)
+#% required: no
+#% options: strds, str3ds, stvds
+#% answer: strds
+#%end
+
+
+#%option
+#% key: fs
+#% type: string
+#% description: The field separator character between the columns, default is tabular " | ". Do not use "," as this char is reserved to list several map ids in a sample granule
+#% required: no
+#%end
+
+#%flag
+#% key: s
+#% description: Consider only the start time of maps in input space time dataset while sampling
+#%end
+
+
+#%flag
+#% key: h
+#% description: Print column names 
+#%end
+
+import grass.script as grass
+import grass.temporal as tgis
+
+############################################################################
+
+def main():
+
+    # Get the options
+    input = options["input"]
+    sampler = options["sample"]
+    sampler_type = options["type"]
+    separator = options["fs"]
+    header = flags["h"]
+    use_simple = flags["s"]
+
+    # Make sure the temporal database exists
+    tgis.create_temporal_database()
+
+    tgis.sample_stds_by_stds_topology("strds", sampler_type, input, sampler, header, separator, use_simple)
+
+if __name__ == "__main__":
+    options, flags = grass.parser()
+    main()