1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
|
"""
Sampling functions for space time datasets
Usage:
.. code-block:: python
import grass.temporal as tgis
tgis.register_maps_in_space_time_dataset(type, name, maps)
(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 .factory import *
def sample_stds_by_stds_topology(intype, sampletype, inputs, sampler, header,
separator, method, spatial=False,
print_only=True):
"""Sample the input space time datasets with a sample
space time dataset, return the created map matrix and optionally
print the result to stdout
In case multiple maps are located in the current granule,
the map names are separated by comma.
In case a layer is present, the names map ids are extended
in this form: "name:layer@mapset"
Attention: Do not use the comma as separator for printing
:param intype: Type of the input space time dataset (strds, stvds or
str3ds)
:param sampletype: Type of the sample space time datasets (strds,
stvds or str3ds)
:param inputs: Name or comma separated names of space time datasets or
a list of map names
:param sampler: Name of a space time dataset used for temporal sampling
:param header: Set True to print column names
:param separator: The field separator character between the columns
:param method: The method to be used for temporal sampling
(start,during,contain,overlap,equal) as comma separated
string or as a list of methods
:param spatial: Perform spatial overlapping check
:param print_only: If set True (default) then the result of the
sampling will be printed to stdout, if set to False
the resulting map matrix will be returned.
:return: The map matrix or None if nothing found
"""
mapset = get_current_mapset()
msgr = get_tgis_message_interface()
# Make a method list
if not issubclass(type(method), type([])):
method = method.split(",")
# Split the inputs
if not issubclass(type(inputs), type([])):
inputs = inputs.split(",")
sts = []
for input in inputs:
if input.find("@") >= 0:
id = input
else:
id = input + "@" + mapset
st = dataset_factory(intype, id)
sts.append(st)
if sampler.find("@") >= 0:
sid = sampler
else:
sid = sampler + "@" + mapset
sst = dataset_factory(sampletype, sid)
dbif = SQLDatabaseInterfaceConnection()
dbif.connect()
for st in sts:
if st.is_in_db(dbif) is False:
msgr.fatal(_("Dataset <%s> not found in temporal database")
% (st.get_id()))
st.select(dbif)
if sst.is_in_db(dbif) is False:
msgr.fatal(_("Dataset <%s> not found in temporal database") % (sid))
sst.select(dbif)
if separator is None or separator == "" or separator.find(",") >= 0:
separator = " | "
mapmatrizes = []
for st in sts:
mapmatrix = st.sample_by_dataset(sst, method, spatial, dbif)
if mapmatrix and len(mapmatrix) > 0:
mapmatrizes.append(mapmatrix)
if len(mapmatrizes) > 0:
# Simply return the map matrix
if not print_only:
dbif.close()
return mapmatrizes
if header:
string = ""
string += "%s%s" % (sst.get_id(), separator)
for st in sts:
string += "%s%s" % (st.get_id(), separator)
string += "%s%s" % ("start_time", separator)
string += "%s%s" % ("end_time", separator)
string += "%s%s" % ("interval_length", separator)
string += "%s" % ("distance_from_begin")
print(string)
first_time, dummy = mapmatrizes[0][0]["granule"].get_temporal_extent_as_tuple()
for i in range(len(mapmatrizes[0])):
mapname_list = []
for mapmatrix in mapmatrizes:
mapnames = ""
count = 0
entry = mapmatrix[i]
for sample in entry["samples"]:
if count == 0:
mapnames += str(sample.get_id())
else:
mapnames += ",%s" % str(sample.get_id())
count += 1
mapname_list.append(mapnames)
entry = mapmatrizes[0][i]
map = entry["granule"]
start, end = map.get_temporal_extent_as_tuple()
if end:
delta = end - start
else:
delta = None
delta_first = start - first_time
if map.is_time_absolute():
if end:
delta = time_delta_to_relative_time(delta)
delta_first = time_delta_to_relative_time(delta_first)
string = ""
string += "%s%s" % (map.get_id(), separator)
for mapnames in mapname_list:
string += "%s%s" % (mapnames, separator)
string += "%s%s" % (start, separator)
string += "%s%s" % (end, separator)
string += "%s%s" % (delta, separator)
string += "%s" % (delta_first)
print(string)
dbif.close()
if len(mapmatrizes) > 0:
return mapmatrizes
return None
|
