"""@package grass.temporal
Temporal raster algebra
(C) 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: Thomas Leppelt and Soeren Gebbert
.. code-block:: python
>>> p = TemporalRasterAlgebraLexer()
>>> p.build()
>>> p.debug = True
>>> expression = 'R = A {+,equal,l} B '
>>> p.test(expression)
R = A {+,equal,l} B
LexToken(NAME,'R',1,0)
LexToken(EQUALS,'=',1,2)
LexToken(NAME,'A',1,4)
LexToken(T_ARITH2_OPERATOR,'{+,equal,l}',1,6)
LexToken(NAME,'B',1,18)
>>> expression = 'R = A {*,equal|during,r} B '
>>> p.test(expression)
R = A {*,equal|during,r} B
LexToken(NAME,'R',1,0)
LexToken(EQUALS,'=',1,2)
LexToken(NAME,'A',1,4)
LexToken(T_ARITH1_OPERATOR,'{*,equal|during,r}',1,6)
LexToken(NAME,'B',1,25)
>>> expression = 'R = A {+,equal|during} B '
>>> p.test(expression)
R = A {+,equal|during} B
LexToken(NAME,'R',1,0)
LexToken(EQUALS,'=',1,2)
LexToken(NAME,'A',1,4)
LexToken(T_ARITH2_OPERATOR,'{+,equal|during}',1,6)
LexToken(NAME,'B',1,23)
"""
import grass.pygrass.modules as pymod
from temporal_operator import *
from temporal_algebra import *
##############################################################################
[docs]class TemporalRasterAlgebraLexer(TemporalAlgebraLexer):
"""Lexical analyzer for the GRASS GIS temporal algebra"""
def __init__(self):
TemporalAlgebraLexer.__init__(self)
# Supported r.mapcalc functions.
mapcalc_functions = {
'exp' : 'EXP',
'log' : 'LOG',
'sqrt' : 'SQRT',
'abs' : 'ABS',
'cos' : 'COS',
'acos' : 'ACOS',
'sin' : 'SIN',
'asin' : 'ASIN',
'tan' : 'TAN',
'double' : 'DOUBLE',
'float' : 'FLOATEXP',
'int' : 'INTEXP',
'isnull' : 'ISNULL',
'isntnull': 'ISNTNULL',
'null' : 'NULL',
'exist' : 'EXIST',
}
# Functions that defines single maps with time stamp and without temporal extent.
map_functions = {'map' : 'MAP'}
# This is the list of token names.
raster_tokens = (
'MOD',
'DIV',
'MULT',
'ADD',
'SUB',
'T_ARITH1_OPERATOR',
'T_ARITH2_OPERATOR',
'L_SPAREN',
'R_SPAREN',
)
# Build the token list
tokens = TemporalAlgebraLexer.tokens \
+ raster_tokens \
+ tuple(mapcalc_functions.values()) \
+ tuple(map_functions.values())
# Regular expression rules for simple tokens
t_MOD = r'[\%]'
t_DIV = r'[\/]'
t_MULT = r'[\*]'
t_ADD = r'[\+]'
t_SUB = r'[-]'
t_T_ARITH1_OPERATOR = r'\{[\%\*\/][,]?[a-zA-Z\| ]*([,])?([lrudi]|left|right|union|disjoint|intersect)?\}'
t_T_ARITH2_OPERATOR = r'\{[+-][,]?[a-zA-Z\| ]*([,])?([lrudi]|left|right|union|disjoint|intersect)?\}'
t_L_SPAREN = r'\['
t_R_SPAREN = r'\]'
# Parse symbols
[docs] def temporal_symbol(self, t):
# Check for reserved words
if t.value in TemporalRasterAlgebraLexer.time_functions.keys():
t.type = TemporalRasterAlgebraLexer.time_functions.get(t.value)
elif t.value in TemporalRasterAlgebraLexer.datetime_functions.keys():
t.type = TemporalRasterAlgebraLexer.datetime_functions.get(t.value)
elif t.value in TemporalRasterAlgebraLexer.conditional_functions.keys():
t.type = TemporalRasterAlgebraLexer.conditional_functions.get(t.value)
elif t.value in TemporalRasterAlgebraLexer.mapcalc_functions.keys():
t.type = TemporalRasterAlgebraLexer.mapcalc_functions.get(t.value)
elif t.value in TemporalRasterAlgebraLexer.map_functions.keys():
t.type = TemporalRasterAlgebraLexer.map_functions.get(t.value)
else:
t.type = 'NAME'
return t
##############################################################################
[docs]class TemporalRasterBaseAlgebraParser(TemporalAlgebraParser):
"""The temporal algebra class"""
# Get the tokens from the lexer class
tokens = TemporalRasterAlgebraLexer.tokens
# Setting equal precedence level for select and hash operations.
precedence = (
('left', 'T_SELECT_OPERATOR', 'T_SELECT', 'T_NOT_SELECT'), # 1
('left', 'ADD', 'SUB', 'T_ARITH2_OPERATOR', 'T_HASH_OPERATOR', 'HASH'), #2
('left', 'AND', 'OR', 'T_COMP_OPERATOR', 'MOD', 'DIV', 'MULT',
'T_ARITH1_OPERATOR'))
def __init__(self, pid=None, run = True, debug = False, spatial = False, \
nprocs = 1, register_null = False):
TemporalAlgebraParser.__init__(self, pid, run, debug, spatial)
self.nprocs = nprocs
self.register_null = register_null
[docs] def check_null(self, t):
try:
int(t)
return t
except ValueError:
return "null()"
######################### Temporal functions ##############################
[docs] def get_temporal_topo_list(self, maplistA, maplistB = None, topolist = ["EQUAL"],
assign_val = False, count_map = False, compare_bool = False,
compare_cmd = False, compop = None, aggregate = None,
new = False, convert = False, operator_cmd = False):
"""Build temporal topology for two space time data sets, copy map objects
for given relation into map list.
:param maplistA: List of maps.
:param maplistB: List of maps.
:param topolist: List of strings of temporal relations.
:param assign_val: Boolean for assigning a boolean map value based on
the map_values from the compared map list by
topological relationships.
:param count_map: Boolean if the number of topological related maps
should be returned.
:param compare_bool: Boolean for comparing boolean map values based on
related map list and compariosn operator.
:param compare_cmd: Boolean for comparing command list values based on
related map list and compariosn operator.
:param compop: Comparison operator, && or ||.
:param aggregate: Aggregation operator for relation map list, & or |.
:param new: Boolean if new temporary maps should be created.
:param convert: Boolean if conditional values should be converted to
r.mapcalc command strings.
:param operator_cmd: Boolean for aggregate arithmetic operators implicitly
in command list values based on related map lists.
:return: List of maps from maplistA that fulfil the topological relationships
to maplistB specified in topolist.
"""
topologylist = ["EQUAL", "FOLLOWS", "PRECEDES", "OVERLAPS", "OVERLAPPED", \
"DURING", "STARTS", "FINISHES", "CONTAINS", "STARTED", \
"FINISHED"]
complementdict = {"EQUAL": "EQUAL", "FOLLOWS" : "PRECEDES",
"PRECEDES" : "FOLLOWS", "OVERLAPS" : "OVERLAPPED",
"OVERLAPPED" : "OVERLAPS", "DURING" : "CONTAINS",
"CONTAINS" : "DURING", "STARTS" : "STARTED",
"STARTED" : "STARTS", "FINISHES" : "FINISHED",
"FINISHED" : "FINISHES"}
resultdict = {}
# Check if given temporal relation are valid.
for topo in topolist:
if topo.upper() not in topologylist:
raise SyntaxError("Unpermitted temporal relation name '" + topo + "'")
# Create temporal topology for maplistA to maplistB.
tb = SpatioTemporalTopologyBuilder()
# Dictionary with different spatial variables used for topology builder.
spatialdict = {'strds' : '2D', 'stvds' : '2D', 'str3ds' : '3D'}
# Build spatial temporal topology
if self.spatial:
tb.build(maplistA, maplistB, spatial = spatialdict[self.stdstype])
else:
tb.build(maplistA, maplistB)
# Iterate through maps in maplistA and search for relationships given
# in topolist.
for map_i in maplistA:
tbrelations = map_i.get_temporal_relations()
# Check for boolean parameters for further calculations.
if assign_val:
self.assign_bool_value(map_i, tbrelations, topolist)
elif compare_bool:
self.compare_bool_value(map_i, tbrelations, compop, aggregate, topolist)
elif compare_cmd:
self.compare_cmd_value(map_i, tbrelations, compop, aggregate, topolist, convert)
elif operator_cmd:
self.operator_cmd_value(map_i, tbrelations, compop, topolist)
for topo in topolist:
if topo.upper() in tbrelations.keys():
if count_map:
relationmaplist = tbrelations[topo.upper()]
gvar = GlobalTemporalVar()
gvar.td = len(relationmaplist)
if "map_value" in dir(map_i):
map_i.map_value.append(gvar)
else:
map_i.map_value = gvar
# Use unique identifier, since map names may be equal
resultdict[map_i.uid] = map_i
resultlist = resultdict.values()
# Sort list of maps chronological.
resultlist = sorted(resultlist, key = AbstractDatasetComparisonKeyStartTime)
return(resultlist)
[docs] def build_command_string(self, map_i, relmap, operator = None, cmd_type = None):
"""This function build the r.mapcalc command string for conditionals,
spatial variable combinations and boolean comparisons.
For Example: 'if(a1 == 1, b1, c2)' or 'exist(a1) && sin(b1)'
:param map_i: map object with temporal extent and built relations.
:param relmap: map object with defined temporal relation to map_i.
:param operator: String representing operator between two spatial variables
(&&,||,+,-,*,/).
:param cmd_type: map object with defined temporal relation to map_i:
condition, conclusion or operator.
:return: the resulting command string for conditionals or spatial variable
combinations
"""
def sub_cmdstring(map_i):
"""This function search for command string in a map object and
return substitute string (contained commandstring or map name)"""
if "cmd_list" in dir(map_i):
map_sub = map_i.cmd_list
elif "map_value" in dir(map_i) and len(map_i.map_value) > 0 and map_i.map_value[0].get_type() == "timediff":
map_sub = map_i.map_value[0].get_type_value()[0]
else:
try:
map_sub = map_i.get_id()
except:
map_sub = map_i
return(map_sub)
# Check for type of operation, conditional or spatial variable combination
# and Create r.mapcalc expression string for the operation.
cmdstring = ""
if cmd_type == 'condition':
conditionsub = sub_cmdstring(map_i)
conclusionsub = sub_cmdstring(relmap)
cmdstring = "if(%s, %s)" %(conditionsub, conclusionsub)
elif cmd_type == 'conclusion':
thensub = sub_cmdstring(map_i)
elsesub = sub_cmdstring(relmap)
cmdstring = "%s, %s" %(thensub, elsesub)
elif cmd_type == 'operator':
leftsub = sub_cmdstring(map_i)
rightsub = sub_cmdstring(relmap)
if operator == None:
self.msgr.fatal("Error: Can't build command string for map %s, operator is missing" \
%(map_i.get_map_id()))
cmdstring = "(%s %s %s)" %(leftsub, operator, rightsub)
return(cmdstring)
[docs] def compare_cmd_value(self, map_i, tbrelations, compop, aggregate,
topolist = ["EQUAL"], convert = False):
""" Function to evaluate two map lists with boolean values by boolean
comparison operator. Extended temporal algebra version with command
list builder for temporal raster algebra.
:param map_i: Map object with temporal extent.
:param tbrelations: List of temporal relation to map_i.
:param topolist: List of strings for given temporal relations.
:param compop: Comparison operator, && or ||.
:param aggregate: Aggregation operator for relation map list, & or |.
:param convert: Boolean if conditional values should be converted to
r.mapcalc command strings.
:return: Map object with conditional value that has been evaluated by
comparison operators.
"""
# Build comandlist list with elements from related maps and given relation operator.
if convert and "condition_value" in dir(map_i):
if map_i.condition_value != []:
cmdstring = str(int(map_i.condition_value[0]))
map_i.cmd_list = cmdstring
if "cmd_list" in dir(map_i):
leftcmd = map_i.cmd_list
cmd_value_list = [leftcmd]
count = 0
for topo in topolist:
if topo.upper() in tbrelations.keys():
relationmaplist = tbrelations[topo.upper()]
if count == 0 and "cmd_list" in dir(map_i):
cmd_value_list.append(compop)
cmd_value_list.append('(')
for relationmap in relationmaplist:
if convert and "condition_value" in dir(relationmap):
if relationmap.condition_value != []:
cmdstring = str(int(relationmap.condition_value[0]))
relationmap.cmd_list = cmdstring
if "cmd_list" in dir(relationmap):
if count > 0:
cmd_value_list.append(aggregate + aggregate)
cmd_value_list.append(relationmap.cmd_list)
count = count + 1
if count > 0:
cmd_value_list.append(')')
cmd_value_str = ''.join(map(str, cmd_value_list))
# Add command list to result map.
map_i.cmd_list = cmd_value_str
return(cmd_value_str)
[docs] def operator_cmd_value(self, map_i, tbrelations, operator, topolist = ["EQUAL"]):
""" Function to evaluate two map lists by given arithmetic operator.
:param map_i: Map object with temporal extent.
:param tbrelations: List of temporal relation to map_i.
:param topolist: List of strings for given temporal relations.
:param operator: Arithmetic operator, +-*/%.
:return: Map object with command list with operators that has been
evaluated by implicit aggregration.
"""
# Build comandlist list with elements from related maps and given relation operator.
leftcmd = map_i
cmdstring = ""
for topo in topolist:
if topo.upper() in tbrelations.keys():
relationmaplist = tbrelations[topo.upper()]
for relationmap in relationmaplist:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(leftcmd,
relationmap,
operator = operator,
cmd_type = "operator")
leftcmd = cmdstring
# Add command list to result map.
map_i.cmd_list = cmdstring
return(cmdstring)
[docs] def set_temporal_extent_list(self, maplist, topolist = ["EQUAL"], temporal = 'l' ,
cmd_bool = False, cmd_type = None, operator = None):
""" Change temporal extent of map list based on temporal relations to
other map list and given temporal operator.
:param maplist: List of map objects for which relations has been build
correctely.
:param topolist: List of strings of temporal relations.
:param temporal: The temporal operator specifying the temporal
extent operation (intersection, union, disjoint
union, right reference, left reference).
:param cmd_bool: Boolean if command string should be merged for related maps.
:param cmd_type: map object with defined temporal relation to map_i:
condition, conclusion or operator.
:param operator: String defining the type of operator.
:return: Map list with specified temporal extent and optional command string.
"""
resultdict = {}
for map_i in maplist:
# Loop over temporal related maps and create overlay modules.
tbrelations = map_i.get_temporal_relations()
# Generate an intermediate map for the result map list.
map_new = self.generate_new_map(base_map=map_i, bool_op = 'and',
copy = True, rename = True)
# Combine temporal and spatial extents of intermediate map with related maps.
for topo in topolist:
if topo in tbrelations.keys():
for map_j in (tbrelations[topo]):
if temporal == 'r':
# Generate an intermediate map for the result map list.
map_new = self.generate_new_map(base_map=map_i, bool_op = 'and',
copy = True, rename = True)
# Create overlayed map extent.
returncode = self.overlay_map_extent(map_new, map_j, 'and', \
temp_op = temporal)
# Stop the loop if no temporal or spatial relationship exist.
if returncode == 0:
break
# Append map to result map list.
elif returncode == 1:
# print(map_new.cmd_list)
# resultlist.append(map_new)
if cmd_bool:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(map_i,
map_j,
operator = operator,
cmd_type = cmd_type)
# Conditional append of module command.
map_new.cmd_list = cmdstring
# Write map object to result dictionary.
resultdict[map_new.uid] = map_new
if returncode == 0:
break
# Append map to result map list.
#if returncode == 1:
# resultlist.append(map_new)
# Get sorted map objects as values from result dictionoary.
resultlist = resultdict.values()
resultlist = sorted(resultlist, key = AbstractDatasetComparisonKeyStartTime)
return(resultlist)
[docs] def build_condition_cmd_list(self, iflist, thenlist, elselist = None,
condition_topolist = ["EQUAL"],
conclusion_topolist = ["EQUAL"],
temporal = 'l', null = False):
"""This function build the r.mapcalc command strings for spatial conditionals.
For Example: 'if(a1 == 1, b1, c2)'
:param iflist: Map list with temporal extents and command list.
:param thenlist: Map list with temporal extents and command list or numeric string.
:param elselist: Map list with temporal extents and command list or numeric string.
:param condition_topolist: List of strings for given temporal relations between
conditions and conclusions.
:param conclusion_topolist: List of strings for given temporal relations between
conditions (then and else).
:param temporal: The temporal operator specifying the temporal
extent operation (intersection, union, disjoint
union, right reference, left reference).
:param null: Boolean if null map support should be activated.
:return: map list with resulting command string for given condition type.
"""
resultlist = []
# First merge conclusion command maplists or strings.
# Check if alternative conclusion map list is given.
if all([isinstance(thenlist, list), isinstance(elselist, list)]):
# Build conclusion command map list.
conclusiontopolist = self.get_temporal_topo_list(thenlist, elselist, conclusion_topolist)
conclusionlist = self.set_temporal_extent_list(conclusiontopolist,
topolist = conclusion_topolist,
temporal = temporal ,
cmd_bool = True,
cmd_type = "conclusion")
# Check if any conclusion is a numeric statements.
elif any([isinstance(thenlist, str), isinstance(elselist, str)]):
conclusionlist = []
# Check if only alternative conclusion is a numeric statements.
if all([isinstance(thenlist, list), isinstance(elselist, str)]):
listinput = thenlist
numinput = elselist
for map_i in listinput:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string( map_i,
numinput,
cmd_type = 'conclusion')
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
conclusionlist.append(map_i)
# Check if only direct conclusion is a numeric statements.
elif all([isinstance(thenlist, str), isinstance(elselist, list)]):
listinput = elselist
numinput = thenlist
for map_i in listinput:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(numinput,
map_i,
cmd_type = 'conclusion')
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
conclusionlist.append(map_i)
elif all([isinstance(thenlist, str), isinstance(elselist, str)]):
conclusionlist = thenlist + ',' + elselist
else:
# The direct conclusion is used.
conclusionlist = thenlist
# Use the conclusion map or string to merge it with the condition and
# return maplist.
if isinstance(conclusionlist, str):
resultlist = []
for map_i in iflist:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(map_i,
conclusionlist,
cmd_type = 'condition')
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
return(resultlist)
elif isinstance(conclusionlist, list):
# Build result command map list between conditions and conclusions.
conditiontopolist = self.get_temporal_topo_list(iflist, conclusionlist,
topolist = condition_topolist)
resultlist = self.set_temporal_extent_list(conditiontopolist,
topolist = condition_topolist,
temporal = 'r' ,
cmd_bool = True,
cmd_type = "condition")
return(resultlist)
###########################################################################
[docs] def p_statement_assign(self, t):
# The expression should always return a list of maps.
"""
statement : stds EQUALS expr
"""
if self.run:
# Create the process queue for parallel mapcalc processing
process_queue = pymod.ParallelModuleQueue(int(self.nprocs))
if isinstance(t[3], list):
num = len(t[3])
count = 0
returncode = 0
register_list = []
for i in range(num):
# Check if resultmap names exist in GRASS database.
map_name = self.basename + "_" + str(i) + "@" + self.mapset
if self.stdstype == "strds":
new_map = RasterDataset(map_name)
else:
new_map = Raster3DDataset(map_name)
if new_map.map_exists() and self.overwrite == False:
self.msgr.fatal("Error maps with basename %s exist. Use --o flag to overwrite existing file" \
%(map_name))
map_test_list = []
for map_i in t[3]:
newident = self.basename + "_" + str(count)
if "cmd_list" in dir(map_i):
# Build r.mapcalc module and execute expression.
# Change map name to given basename.
# Create deepcopy of r.mapcalc module.
map_test = map_i.get_new_instance(newident + "@" + self.mapset)
map_test.set_temporal_extent(map_i.get_temporal_extent())
map_test.set_spatial_extent(map_i.get_spatial_extent())
map_test_list.append(map_test)
m = copy.deepcopy(self.m_mapcalc)
m_expression = newident + "=" + map_i.cmd_list
m.inputs["expression"].value = str(m_expression)
m.flags["overwrite"].value = self.overwrite
process_queue.put(m)
elif map_i.map_exists():
# Copy map if it exists
map_test = map_i.get_new_instance(newident + "@" + self.mapset)
map_test.set_temporal_extent(map_i.get_temporal_extent())
map_test.set_spatial_extent(map_i.get_spatial_extent())
map_test_list.append(map_test)
m = copy.deepcopy(self.m_mapcalc)
m_expression = newident + "=" + map_i.get_map_id()
m.inputs["expression"].value = str(m_expression)
m.flags["overwrite"].value = self.overwrite
print m.get_bash()
process_queue.put(m)
else:
self.msgr.error(_("Error computing map <%s>"%(map_i.get_id()) ))
count += 1
process_queue.wait()
for map_i in map_test_list:
register_list.append(map_i)
# Open connection to temporal database.
dbif, connect = init_dbif(self.dbif)
# Create result space time dataset.
resultstds = open_new_stds(t[1], self.stdstype, \
'absolute', t[1], t[1], \
'mean', self.dbif, \
overwrite = self.overwrite)
for map_i in register_list:
# Get meta data from grass database.
map_i.load()
# Do not register empty maps if not required
# In case of a null map continue, do not register null maps
if map_i.metadata.get_min() is None and \
map_i.metadata.get_max() is None:
if not self.register_null:
self.removable_maps[map_i.get_name()] = map_i
continue
if map_i.is_in_db(dbif) and self.overwrite:
# Update map in temporal database.
map_i.update_all(dbif)
elif map_i.is_in_db(dbif) and self.overwrite == False:
# Raise error if map exists and no overwrite flag is given.
self.msgr.fatal("Error raster map %s exist in temporal database. Use overwrite flag. : \n%s" \
%(map_i.get_map_id(), cmd.popen.stderr))
else:
# Insert map into temporal database.
map_i.insert(dbif)
# Register map in result space time dataset.
success = resultstds.register_map(map_i, dbif)
resultstds.update_from_registered_maps(dbif)
dbif.close()
t[0] = register_list
self.remove_maps()
[docs] def p_expr_spmap_function(self, t):
# Add a single map.
# Only the spatial extent of the map is evaluated.
# Temporal extent is not existing.
# Examples:
# R = map(A)
"""
mapexpr : MAP LPAREN stds RPAREN
"""
if self.run:
# Check input map.
input = t[3]
if not isinstance(input, list):
# Check for mapset in given stds input.
if input.find("@") >= 0:
id_input = input
else:
id_input = input + "@" + self.mapset
# Create empty map dataset.
map_i = dataset_factory(self.maptype, id_input)
# Check for occurrence of space time dataset.
if map_i.map_exists() == False:
raise FatalError(_("%s map <%s> not found in GRASS spatial database") %
(map_i.get_type(), id_input))
else:
# Select dataset entry from database.
map_i.select(dbif=self.dbif)
# Create command list for map object.
cmdstring = "(%s)" %(map_i.get_map_id())
map_i.cmd_list = cmdstring
# Return map object.
t[0] = cmdstring
else:
t[0] = "map(" + t[3] + ")"
if self.debug:
print "map(" + t[3] + ")"
[docs] def p_arith1_operation(self, t):
"""
expr : stds MOD stds
| expr MOD stds
| stds MOD expr
| expr MOD expr
| stds DIV stds
| expr DIV stds
| stds DIV expr
| expr DIV expr
| stds MULT stds
| expr MULT stds
| stds MULT expr
| expr MULT expr
| stds MOD t_td_var
| expr MOD t_td_var
| stds DIV t_td_var
| expr DIV t_td_var
| stds MULT t_td_var
| expr MULT t_td_var
"""
# Check input stds.
maplistA = self.check_stds(t[1])
maplistB = self.check_stds(t[3])
topolist = self.get_temporal_topo_list(maplistA, maplistB)
if self.run:
resultlist = []
for map_i in topolist:
# Generate an intermediate map for the result map list.
map_new = self.generate_new_map(base_map=map_i, bool_op = 'and', copy = True)
# Loop over temporal related maps and create overlay modules.
tbrelations = map_i.get_temporal_relations()
count = 0
for map_j in (tbrelations['EQUAL']):
# Create overlayed map extent.
returncode = self.overlay_map_extent(map_new, map_j, 'and', \
temp_op = 'l')
# Stop the loop if no temporal or spatial relationship exist.
if returncode == 0:
break
if count == 0:
# Set map name.
name = map_new.get_id()
else:
# Generate an intermediate map
name = self.generate_map_name()
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(map_i, map_j,
operator = t[2],
cmd_type = "operator")
# Conditional append of module command.
map_new.cmd_list = cmdstring
count += 1
# Append map to result map list.
if returncode == 1:
resultlist.append(map_new)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith1_operation_numeric1(self, t):
"""
expr : stds MOD number
| expr MOD number
| stds DIV number
| expr DIV number
| stds MULT number
| expr MULT number
| stds MOD numberstr
| expr MOD numberstr
| stds DIV numberstr
| expr DIV numberstr
| stds MULT numberstr
| expr MULT numberstr
| stds MOD mapexpr
| expr MOD mapexpr
| stds DIV mapexpr
| expr DIV mapexpr
| stds MULT mapexpr
| expr MULT mapexpr
"""
# Check input stds.
maplist = self.check_stds(t[1])
if self.run:
resultlist = []
for map_i in maplist:
mapinput = map_i.get_id()
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "(%s %s %s)" %(map_i.cmd_list, t[2], t[3])
else:
cmdstring = "(%s %s %s)" %(mapinput, t[2], t[3])
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith1_operation_numeric2(self, t):
"""
expr : number MOD stds
| number MOD expr
| number DIV stds
| number DIV expr
| number MULT stds
| number MULT expr
| numberstr MOD stds
| numberstr MOD expr
| numberstr DIV stds
| numberstr DIV expr
| numberstr MULT stds
| numberstr MULT expr
| mapexpr MOD stds
| mapexpr MOD expr
| mapexpr DIV stds
| mapexpr DIV expr
| mapexpr MULT stds
| mapexpr MULT expr
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
mapinput = map_i.get_id()
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "(%s %s %s)" %(t[1], t[2], map_i.cmd_list)
else:
cmdstring = "(%s %s %s)" %(t[1], t[2], mapinput)
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith2_operation(self, t):
"""
expr : stds ADD stds
| expr ADD stds
| stds ADD expr
| expr ADD expr
| stds SUB stds
| expr SUB stds
| stds SUB expr
| expr SUB expr
| stds ADD t_td_var
| expr ADD t_td_var
| expr SUB t_td_var
| stds SUB t_td_var
"""
# Check input stds.
maplistA = self.check_stds(t[1])
maplistB = self.check_stds(t[3])
topolist = self.get_temporal_topo_list(maplistA, maplistB)
if self.run:
resultlist = []
for map_i in topolist:
# Generate an intermediate map for the result map list.
map_new = self.generate_new_map(base_map=map_i, bool_op = 'and', copy = True)
# Loop over temporal related maps and create overlay modules.
tbrelations = map_i.get_temporal_relations()
count = 0
for map_j in (tbrelations['EQUAL']):
# Create overlayed map extent.
returncode = self.overlay_map_extent(map_new, map_j, 'and', \
temp_op = 'l')
# Stop the loop if no temporal or spatial relationship exist.
if returncode == 0:
break
if count == 0:
# Set map name.
name = map_new.get_id()
else:
# Generate an intermediate map
name = self.generate_map_name()
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(map_i, map_j,
operator = t[2],
cmd_type = "operator")
# Conditional append of module command.
map_new.cmd_list = cmdstring
count += 1
# Append map to result map list.
if returncode == 1:
resultlist.append(map_new)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith2_operation_numeric1(self, t):
"""
expr : stds ADD number
| expr ADD number
| stds SUB number
| expr SUB number
| stds ADD numberstr
| expr ADD numberstr
| stds SUB numberstr
| expr SUB numberstr
| stds ADD mapexpr
| expr ADD mapexpr
| stds SUB mapexpr
| expr SUB mapexpr
"""
# Check input stds.
maplist = self.check_stds(t[1])
if self.run:
resultlist = []
for map_i in maplist:
mapinput = map_i.get_id()
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "(%s %s %s)" %(map_i.cmd_list, t[2], t[3])
else:
cmdstring = "(%s %s %s)" %(mapinput, t[2], t[3])
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith2_operation_numeric2(self, t):
"""
expr : number ADD stds
| number ADD expr
| number SUB stds
| number SUB expr
| numberstr ADD stds
| numberstr ADD expr
| numberstr SUB stds
| numberstr SUB expr
| mapexpr ADD stds
| mapexpr ADD expr
| mapexpr SUB stds
| mapexpr SUB expr
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
mapinput = map_i.get_id()
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "(%s %s %s)" %(t[1], t[2], map_i.cmd_list)
else:
cmdstring = "(%s %s %s)" %(t[1], t[2], mapinput)
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith1_operation_relation(self, t):
"""
expr : stds T_ARITH1_OPERATOR stds
| expr T_ARITH1_OPERATOR stds
| stds T_ARITH1_OPERATOR expr
| expr T_ARITH1_OPERATOR expr
| stds T_ARITH1_OPERATOR t_td_var
| expr T_ARITH1_OPERATOR t_td_var
"""
if self.run:
# Check input stds.
maplistA = self.check_stds(t[1])
maplistB = self.check_stds(t[3])
relations, temporal, function, aggregate = self.eval_toperator(t[2], optype = 'raster')
# Build conditional values based on topological relationships.
complist = self.get_temporal_topo_list(maplistA, maplistB, topolist = relations,
operator_cmd = True, compop = function)
# Set temporal extent based on topological relationships.
resultlist = self.set_temporal_extent_list(complist, topolist = relations,
temporal = temporal)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith2_operation_relation(self, t):
"""
expr : stds T_ARITH2_OPERATOR stds
| expr T_ARITH2_OPERATOR stds
| stds T_ARITH2_OPERATOR expr
| expr T_ARITH2_OPERATOR expr
| stds T_ARITH2_OPERATOR t_td_var
| expr T_ARITH2_OPERATOR t_td_var
"""
if self.run:
# Check input stds.
maplistA = self.check_stds(t[1])
maplistB = self.check_stds(t[3])
relations, temporal, function, aggregate = self.eval_toperator(t[2], optype = 'raster')
# Build conditional values based on topological relationships.
complist = self.get_temporal_topo_list(maplistA, maplistB, topolist = relations,
operator_cmd = True, compop = function)
# Set temporal extent based on topological relationships.
resultlist = self.set_temporal_extent_list(complist, topolist = relations,
temporal = temporal)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_arith_operation_numeric_string(self, t):
"""
numberstr : number ADD number
| number SUB number
| number DIV number
| number MULT number
| number MOD number
"""
numstring = "(%s %s %s)" %(t[1], t[2], t[3])
t[0] = numstring
if self.debug:
print(numstring)
[docs] def p_mapcalc_function(self, t):
# Supported mapcalc functions.
"""
mapcalc_arith : ABS
| LOG
| SQRT
| EXP
| COS
| ACOS
| SIN
| ASIN
| TAN
| DOUBLE
| FLOATEXP
| INTEXP
"""
t[0] = t[1]
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_mapcalc_operation1(self, t):
# Examples:
# sin(A)
# log(B)
"""
expr : mapcalc_arith LPAREN stds RPAREN
| mapcalc_arith LPAREN expr RPAREN
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "%s(%s)" %(t[1].lower(), map_i.cmd_list)
else:
cmdstring = "%s(%s)" %(t[1].lower(), map_i.get_id())
# Set new command list for map.
map_i.cmd_list = cmdstring
# Append map with updated command list to result list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_mapexpr_operation(self, t):
# Examples:
# sin(map(a))
"""
mapexpr : mapcalc_arith LPAREN mapexpr RPAREN
"""
# Check input stds.
mapstring = t[3]
if self.run:
cmdstring = "%s(%s)" %(t[1].lower(), mapstring)
t[0] = cmdstring
if self.debug:
print mapstring
[docs] def p_s_var_expr_1(self, t):
# Examples:
# isnull(A)
"""
s_var_expr : ISNULL LPAREN stds RPAREN
| ISNULL LPAREN expr RPAREN
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "%s(%s)" %(t[1].lower(), map_i.cmd_list)
else:
cmdstring = "%s(%s)" %(t[1].lower(), map_i.get_id())
# Set new command list for map.
map_i.cmd_list = cmdstring
# Append map with updated command list to result list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_var_expr_2(self, t):
# Examples:
# isntnull(A)
"""
s_var_expr : ISNTNULL LPAREN stds RPAREN
| ISNTNULL LPAREN expr RPAREN
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "!isnull(%s)" %(map_i.cmd_list)
else:
cmdstring = "!isnull(%s)" %(map_i.get_id())
# Set new command list for map.
map_i.cmd_list = cmdstring
# Append map with updated command list to result list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_var_expr_3(self, t):
# Examples:
# A <= 2
"""
s_var_expr : stds comp_op number
| expr comp_op number
"""
# Check input stds.
maplist = self.check_stds(t[1])
if self.run:
resultlist = []
for map_i in maplist:
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "%s %s %s" %(map_i.cmd_list, t[2], t[3])
else:
cmdstring = "%s %s %s" %(map_i.get_id(), t[2], t[3])
# Set new command list for map.
map_i.cmd_list = cmdstring
# Append map with updated command list to result list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_var_expr_4(self, t):
# Examples:
# exist(B)
"""
s_var_expr : EXIST LPAREN stds RPAREN
| EXIST LPAREN expr RPAREN
"""
# Check input stds.
maplist = self.check_stds(t[3])
if self.run:
resultlist = []
for map_i in maplist:
# Create r.mapcalc expression string for the operation.
if "cmd_list" in dir(map_i):
cmdstring = "%s" %(map_i.cmd_list)
else:
cmdstring = "%s" %(map_i.get_id())
# Set new command list for map.
map_i.cmd_list = cmdstring
# Append map with updated command list to result list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_var_expr_comp(self, t):
# Examples:
# A <= 2 || B == 10
# A < 3 && A > 1
"""
s_var_expr : s_var_expr AND AND s_var_expr
| s_var_expr OR OR s_var_expr
"""
if self.run:
# Check input stds.
s_var_exprA = self.check_stds(t[1])
s_var_exprB = self.check_stds(t[4])
relations = ["EQUAL"]
temporal = "l"
function = t[2] + t[3]
aggregate = t[2]
# Build conditional values based on topological relationships.
complist = self.get_temporal_topo_list(s_var_exprA, s_var_exprB, topolist = relations,
compare_cmd = True, compop = function, aggregate = aggregate)
# Set temporal extent based on topological relationships.
resultlist = self.set_temporal_extent_list(complist, topolist = relations,
temporal = temporal)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_var_expr_comp_op(self, t):
# Examples:
# A <= 2 {||} B == 10
# A < 3 {&&, equal} A > 1
"""
s_var_expr : s_var_expr T_COMP_OPERATOR s_var_expr
"""
if self.run:
# Check input stds.
s_var_exprA = self.check_stds(t[1])
s_var_exprB = self.check_stds(t[3])
# Evaluate temporal comparison operator.
relations, temporal, function, aggregate = self.eval_toperator(t[2], optype = 'boolean')
# Build conditional values based on topological relationships.
complist = self.get_temporal_topo_list(s_var_exprA, s_var_exprB, topolist = relations,
compare_cmd = True, compop = function, aggregate = aggregate)
# Set temporal extent based on topological relationships.
resultlist = self.set_temporal_extent_list(complist, topolist = relations,
temporal = temporal)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_expr_condition_if(self, t):
# Examples:
# if(s_var_expr, B)
# if(A == 1, B)
"""
expr : IF LPAREN s_var_expr COMMA stds RPAREN
| IF LPAREN s_var_expr COMMA expr RPAREN
| IF LPAREN ts_var_expr COMMA stds RPAREN
| IF LPAREN ts_var_expr COMMA expr RPAREN
"""
ifmaplist = self.check_stds(t[3])
thenmaplist = self.check_stds(t[5])
resultlist = self.build_condition_cmd_list(ifmaplist, thenmaplist,
elselist = None,
condition_topolist = ["EQUAL"],
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_numeric_condition_if(self, t):
# Examples:
# if(s_var_expr, 1)
# if(A == 5, 10)
"""
expr : IF LPAREN s_var_expr COMMA number RPAREN
| IF LPAREN s_var_expr COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN ts_var_expr COMMA number RPAREN
| IF LPAREN ts_var_expr COMMA NULL LPAREN RPAREN RPAREN
"""
ifmaplist = self.check_stds(t[3])
resultlist = []
# Select input for r.mapcalc expression based on length of PLY object.
if len(t) == 7:
numinput = str(t[5])
elif len(t) == 9:
numinput = str(t[5] + t[6] + t[7])
# Iterate over condition map list.
for map_i in ifmaplist:
# Create r.mapcalc expression string for the operation.
cmdstring = self.build_command_string(map_i, numinput, cmd_type = 'condition')
# Conditional append of module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_expr_condition_if_relation(self, t):
# Examples:
# if({equal||during}, s_var_expr, A)
"""
expr : IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA expr RPAREN
"""
relations, temporal, function, aggregation = self.eval_toperator(t[3],
optype = 'relation')
ifmaplist = self.check_stds(t[5])
thenmaplist = self.check_stds(t[7])
resultlist = self.build_condition_cmd_list(ifmaplist, thenmaplist,
elselist = None,
condition_topolist = relations,
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_expr_condition_elif(self, t):
# Examples:
# if(s_var_expr, A, B)
"""
expr : IF LPAREN s_var_expr COMMA stds COMMA stds RPAREN
| IF LPAREN s_var_expr COMMA stds COMMA expr RPAREN
| IF LPAREN s_var_expr COMMA expr COMMA stds RPAREN
| IF LPAREN s_var_expr COMMA expr COMMA expr RPAREN
| IF LPAREN ts_var_expr COMMA stds COMMA stds RPAREN
| IF LPAREN ts_var_expr COMMA stds COMMA expr RPAREN
| IF LPAREN ts_var_expr COMMA expr COMMA stds RPAREN
| IF LPAREN ts_var_expr COMMA expr COMMA expr RPAREN
"""
# Check map list inputs.
ifmaplist = self.check_stds(t[3])
thenmaplist = self.check_stds(t[5])
elsemaplist = self.check_stds(t[7])
# Create conditional command map list.
resultlist = self.build_condition_cmd_list(ifmaplist, thenmaplist,
elselist = elsemaplist,
condition_topolist = ["EQUAL"],
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_numeric_condition_elif(self, t):
# Examples:
# if(s_var_expr, 1, 2)
# if(A == 5, 10, 0)
"""
expr : IF LPAREN s_var_expr COMMA number COMMA number RPAREN
| IF LPAREN s_var_expr COMMA NULL LPAREN RPAREN COMMA number RPAREN
| IF LPAREN s_var_expr COMMA number COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN s_var_expr COMMA NULL LPAREN RPAREN COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN ts_var_expr COMMA number COMMA number RPAREN
| IF LPAREN ts_var_expr COMMA NULL LPAREN RPAREN COMMA number RPAREN
| IF LPAREN ts_var_expr COMMA number COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN ts_var_expr COMMA NULL LPAREN RPAREN COMMA NULL LPAREN RPAREN RPAREN
"""
ifmaplist = self.check_stds(t[3])
resultlist = []
# Select input for r.mapcalc expression based on length of PLY object.
if len(t) == 9:
numthen = t[5]
numelse = t[7]
elif len(t) == 11 and t[6] == '(':
numthen = t[5] + t[6] + t[7]
numelse = t[9]
elif len(t) == 11 and t[6] == ',':
numthen = t[5]
numelse = t[7] + t[8] + t[9]
elif len(t) == 13:
numthen = t[5] + t[6] + t[7]
numelse = t[9] + t[10] + t[11]
numthen = str(numthen)
numelse = str(numelse)
print(numthen + " " +numelse )
# Create conditional command map list.
resultlist = self.build_condition_cmd_list(ifmaplist, numthen, numelse,
condition_topolist = ["EQUAL"],
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_numeric_expr_condition_elif(self, t):
# Examples:
# if(s_var_expr, 1, A)
# if(A == 5 && C > 5, A, null())
"""
expr : IF LPAREN s_var_expr COMMA number COMMA stds RPAREN
| IF LPAREN s_var_expr COMMA NULL LPAREN RPAREN COMMA stds RPAREN
| IF LPAREN s_var_expr COMMA number COMMA expr RPAREN
| IF LPAREN s_var_expr COMMA NULL LPAREN RPAREN COMMA expr RPAREN
| IF LPAREN s_var_expr COMMA stds COMMA number RPAREN
| IF LPAREN s_var_expr COMMA stds COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN s_var_expr COMMA expr COMMA number RPAREN
| IF LPAREN s_var_expr COMMA expr COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN ts_var_expr COMMA number COMMA stds RPAREN
| IF LPAREN ts_var_expr COMMA NULL LPAREN RPAREN COMMA stds RPAREN
| IF LPAREN ts_var_expr COMMA number COMMA expr RPAREN
| IF LPAREN ts_var_expr COMMA NULL LPAREN RPAREN COMMA expr RPAREN
| IF LPAREN ts_var_expr COMMA stds COMMA number RPAREN
| IF LPAREN ts_var_expr COMMA stds COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN ts_var_expr COMMA expr COMMA number RPAREN
| IF LPAREN ts_var_expr COMMA expr COMMA NULL LPAREN RPAREN RPAREN
"""
ifmaplist = self.check_stds(t[3])
# Select input for r.mapcalc expression based on length of PLY object.
if len(t) == 9:
if isinstance(t[5], int):
theninput = str(t[5])
elseinput = self.check_stds(t[7])
elif isinstance(t[7], int):
theninput = self.check_stds(t[5])
elseinput = str(t[7])
elif len(t) == 11:
if t[5] == 'null':
theninput = str(t[5] + t[6] + t[7])
elseinput = self.check_stds(t[9])
elif t[7] == 'null':
theninput = self.check_stds(t[5])
elseinput = str(t[7] + t[8] + t[9])
# Create conditional command map list.
resultlist = self.build_condition_cmd_list(ifmaplist, theninput, elseinput,
condition_topolist = ["EQUAL"],
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_numeric_expr_condition_elif_relation(self, t):
# Examples:
# if({during},s_var_expr, 1, A)
# if({during}, A == 5, A, null())
"""
expr : IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA number COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA NULL LPAREN RPAREN COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA number COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA NULL LPAREN RPAREN COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA stds COMMA number RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA stds COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA expr COMMA number RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA expr COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA number COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA NULL LPAREN RPAREN COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA number COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA NULL LPAREN RPAREN COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA stds COMMA number RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA stds COMMA NULL LPAREN RPAREN RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA expr COMMA number RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA expr COMMA NULL LPAREN RPAREN RPAREN
"""
relations, temporal, function, aggregation = self.eval_toperator(t[3], optype = 'relation')
ifmaplist = self.check_stds(t[5])
# Select input for r.mapcalc expression based on length of PLY object.
if len(t) == 11:
if isinstance(t[7], int):
theninput = str(t[7])
elseinput = self.check_stds(t[9])
elif isinstance(t[9], int):
theninput = self.check_stds(t[7])
elseinput = str(t[9])
elif len(t) == 13:
if isinstance(t[7], int):
theninput = str(t[7] + t[8] + t[9])
elseinput = self.check_stds(t[11])
elif isinstance(t[9], int):
theninput = self.check_stds(t[7])
elseinput = str(t[9] + t[10] + t[11])
# Create conditional command map list.
resultlist = self.build_condition_cmd_list(ifmaplist, theninput, elseinput,
condition_topolist = relations,
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_s_expr_condition_elif_relation(self, t):
# Examples:
# if({equal||during}, s_var_expr, A, B)
"""
expr : IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA stds COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA stds COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA expr COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA s_var_expr COMMA expr COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA stds COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA stds COMMA expr RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA expr COMMA stds RPAREN
| IF LPAREN T_REL_OPERATOR COMMA ts_var_expr COMMA expr COMMA expr RPAREN
"""
relations, temporal, function, aggregation = self.eval_toperator(t[3], optype = 'relation')
ifmaplist = self.check_stds(t[5])
thenmaplist = self.check_stds(t[7])
elsemaplist = self.check_stds(t[9])
# Create conditional command map list.
resultlist = self.build_condition_cmd_list(ifmaplist, thenmaplist, elsemaplist,
condition_topolist = relations,
conclusion_topolist = ["EQUAL"],
temporal = 'r', null = False)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
[docs] def p_ts_var_expr1(self, t):
# Combination of spatial and temporal conditional expressions.
# Examples:
# A <= 2 || start_date <= 2013-01-01
# end_date > 2013-01-15 && A > 10
# IMPORTANT: Only the intersection of map lists in conditionals are
# exported.
"""
ts_var_expr : s_var_expr AND AND t_var_expr
| t_var_expr AND AND s_var_expr
| t_var_expr OR OR s_var_expr
| s_var_expr OR OR t_var_expr
| ts_var_expr AND AND s_var_expr
| ts_var_expr AND AND t_var_expr
| ts_var_expr OR OR s_var_expr
| ts_var_expr OR OR t_var_expr
| s_var_expr AND AND ts_var_expr
| t_var_expr AND AND ts_var_expr
| s_var_expr OR OR ts_var_expr
| t_var_expr OR OR ts_var_expr
"""
if self.run:
# Check input stds.
s_var_exprA = self.check_stds(t[1])
s_var_exprB = self.check_stds(t[4])
relations = ["EQUAL"]
temporal = "l"
function = t[2] + t[3]
aggregate = t[2]
# Build conditional values based on topological relationships.
complist = self.get_temporal_topo_list(s_var_exprA, s_var_exprB,
topolist = relations, compare_cmd = True,
compop = function, aggregate = aggregate,
convert = True)
# Set temporal extent based on topological relationships.
resultlist = self.set_temporal_extent_list(complist, topolist = relations,
temporal = temporal)
t[0] = resultlist
[docs] def p_hash_operation(self, t):
# Calculate the number of maps within an interval of another map from a
# second space time dataset.
# A # B
# A {equal,r#} B
"""
expr : t_hash_var
"""
# Check input stds.
maplist = self.check_stds(t[1])
if self.run:
resultlist = []
for map_i in maplist:
for obj in map_i.map_value:
if isinstance(obj, GlobalTemporalVar):
n_maps = obj.td
mapinput = map_i.get_id()
# Create r.mapcalc expression string for the operation.
cmdstring = "(%s)" %(n_maps)
# Append module command.
map_i.cmd_list = cmdstring
# Append map to result map list.
resultlist.append(map_i)
t[0] = resultlist
if self.debug:
for map in resultlist:
print map.cmd_list
###############################################################################
if __name__ == "__main__":
import doctest
doctest.testmod()
Help Index | Topics Index | Keywords Index | Full Index
© 2003-2018 GRASS Development Team, GRASS GIS 7.0.7svn Reference Manual