module.py 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492
  1. # -*- coding: utf-8 -*-
  2. """
  3. Created on Tue Apr 2 18:41:27 2013
  4. @author: pietro
  5. @code
  6. >>> import grass.pygrass.modules as pymod
  7. >>> import copy
  8. >>> region = pymod.Module("g.region")
  9. >>> region.flags["p"].value = True
  10. >>> region.flags["u"].value = True
  11. >>> region.flags["3"].value = True
  12. >>> region.get_bash()
  13. 'g.region -p -u -3'
  14. >>> new_region = copy.deepcopy(region)
  15. >>> new_region.inputs["res"].value = "10"
  16. >>> new_region.get_bash()
  17. 'g.region res=10 -p -3'
  18. >>> neighbors = pymod.Module("r.neighbors")
  19. >>> neighbors.inputs["input"].value = "mapA"
  20. >>> neighbors.outputs["output"].value = "mapB"
  21. >>> neighbors.inputs["size"].value = 5
  22. >>> neighbors.get_bash()
  23. 'r.neighbors input=mapA method=average size=5 quantile=0.5 output=mapB'
  24. >>> new_neighbors1 = copy.deepcopy(neighbors)
  25. >>> new_neighbors1.inputs["input"].value = "mapD"
  26. >>> new_neighbors1.inputs["size"].value = 3
  27. >>> new_neighbors1.get_bash()
  28. 'r.neighbors input=mapD method=average size=3 quantile=0.5 output=mapB'
  29. >>> new_neighbors2 = copy.deepcopy(neighbors)
  30. >>> new_neighbors2(input="mapD", size=3, run_=False)
  31. >>> new_neighbors2.get_bash()
  32. 'r.neighbors input=mapD method=average size=3 quantile=0.5 output=mapB'
  33. >>> neighbors = pymod.Module("r.neighbors")
  34. >>> neighbors.get_bash()
  35. 'r.neighbors method=average size=3 quantile=0.5'
  36. >>> new_neighbors3 = copy.deepcopy(neighbors)
  37. >>> new_neighbors3(input="mapA", size=3, output="mapB", run_=False)
  38. >>> new_neighbors3.get_bash()
  39. 'r.neighbors input=mapA method=average size=3 quantile=0.5 output=mapB'
  40. @endcode
  41. """
  42. from __future__ import (nested_scopes, generators, division, absolute_import,
  43. with_statement, print_function, unicode_literals)
  44. import sys
  45. if sys.version_info[0] == 2:
  46. from itertools import izip_longest as zip_longest
  47. else:
  48. from itertools import zip_longest
  49. from xml.etree.ElementTree import fromstring
  50. import time
  51. from grass.script.core import Popen, PIPE
  52. from grass.pygrass.errors import GrassError, ParameterError
  53. from .parameter import Parameter
  54. from .flag import Flag
  55. from .typedict import TypeDict
  56. from .read import GETFROMTAG, DOC
  57. class ParallelModuleQueue(object):
  58. """!This class is designed to run an arbitrary number of pygrass Module
  59. processes in parallel.
  60. Objects of type grass.pygrass.modules.Module can be put into the
  61. queue using put() method. When the queue is full with the maximum
  62. number of parallel processes it will wait for all processes to finish,
  63. sets the stdout and stderr of the Module object and removes it
  64. from the queue when its finished.
  65. This class will raise a GrassError in case a Module process exits
  66. with a return code other than 0.
  67. Usage:
  68. @code
  69. >>> import copy
  70. >>> import grass.pygrass.modules as pymod
  71. >>> mapcalc_list = []
  72. >>> mapcalc = pymod.Module("r.mapcalc",
  73. ... overwrite=True,
  74. ... run_=False,
  75. ... finish_=False)
  76. >>> queue = pymod.ParallelModuleQueue(max_num_procs=3)
  77. >>> for i in xrange(5):
  78. ... new_mapcalc = copy.deepcopy(mapcalc)
  79. ... mapcalc_list.append(new_mapcalc)
  80. ... new_mapcalc(expression="test_pygrass_%i = %i"%(i, i))
  81. ... queue.put(new_mapcalc)
  82. >>> queue.wait()
  83. >>> for mapcalc in mapcalc_list:
  84. ... print(mapcalc.popen.returncode)
  85. 0
  86. 0
  87. 0
  88. 0
  89. 0
  90. @endcode
  91. """
  92. def __init__(self, max_num_procs=1):
  93. """!Constructor
  94. @param max_num_procs The maximum number of Module processes that
  95. can be run in parallel
  96. """
  97. self._num_procs = int(max_num_procs)
  98. self._list = int(max_num_procs) * [None]
  99. self._proc_count = 0
  100. def put(self, module):
  101. """!Put the next Module object in the queue
  102. To run the Module objects in parallel the run_ and finish_ options
  103. of the Module must be set to False.
  104. @param module A preconfigured Module object with run_ and finish_
  105. set to False
  106. """
  107. self._list[self._proc_count] = module
  108. self._list[self._proc_count].run()
  109. self._proc_count += 1
  110. if self._proc_count == self._num_procs:
  111. self.wait()
  112. def get(self, num):
  113. """!Get a Module object from the queue
  114. @param num The number of the object in queue
  115. @return The Module object or None if num is not in the queue
  116. """
  117. if num < self._num_procs:
  118. return self._list[num]
  119. return None
  120. def get_num_run_procs(self):
  121. """!Get the number of Module processes that are in the queue running
  122. or finished
  123. @return The maximum number fo Module processes running/finished in
  124. the queue
  125. """
  126. return len(self._list)
  127. def get_max_num_procs(self):
  128. """!Return the maximum number of parallel Module processes
  129. """
  130. return self._num_procs
  131. def set_max_num_procs(self, max_num_procs):
  132. """!Set the maximum number of Module processes that should run
  133. in parallel
  134. """
  135. self._num_procs = int(max_num_procs)
  136. self.wait()
  137. def wait(self):
  138. """!Wait for all Module processes that are in the list to finish
  139. and set the modules stdout and stderr output options
  140. """
  141. for proc in self._list:
  142. if proc:
  143. stdout, stderr = proc.popen.communicate(input=proc.stdin)
  144. proc.outputs['stdout'].value = stdout if stdout else ''
  145. proc.outputs['stderr'].value = stderr if stderr else ''
  146. if proc.popen.returncode != 0:
  147. GrassError(("Error running module %s") % (proc.name))
  148. self._list = self._num_procs * [None]
  149. self._proc_count = 0
  150. class Module(object):
  151. """
  152. Python allow developers to not specify all the arguments and
  153. keyword arguments of a method or function.
  154. ::
  155. def f(*args):
  156. for arg in args:
  157. print arg
  158. therefore if we call the function like: ::
  159. >>> f('grass', 'gis', 'modules')
  160. grass
  161. gis
  162. modules
  163. or we can define a new list: ::
  164. >>> words = ['grass', 'gis', 'modules']
  165. >>> f(*words)
  166. grass
  167. gis
  168. modules
  169. we can do the same with keyword arguments, rewrite the above function: ::
  170. def f(*args, **kargs):
  171. for arg in args:
  172. print arg
  173. for key, value in kargs.items():
  174. print "%s = %r" % (key, value)
  175. now we can use the new function, with: ::
  176. >>> f('grass', 'gis', 'modules', os = 'linux', language = 'python')
  177. grass
  178. gis
  179. modules
  180. os = 'linux'
  181. language = 'python'
  182. or, as before we can, define a dictionary and give the dictionary to
  183. the function, like: ::
  184. >>> keywords = {'os' : 'linux', 'language' : 'python'}
  185. >>> f(*words, **keywords)
  186. grass
  187. gis
  188. modules
  189. os = 'linux'
  190. language = 'python'
  191. In the Module class we heavily use this language feature to pass arguments
  192. and keyword arguments to the grass module.
  193. """
  194. def __init__(self, cmd, *args, **kargs):
  195. self.name = cmd
  196. try:
  197. # call the command with --interface-description
  198. get_cmd_xml = Popen([cmd, "--interface-description"], stdout=PIPE)
  199. except OSError as e:
  200. print("OSError error({0}): {1}".format(e.errno, e.strerror))
  201. str_err = "Error running: `%s --interface-description`."
  202. raise GrassError(str_err % self.name)
  203. # get the xml of the module
  204. self.xml = get_cmd_xml.communicate()[0]
  205. # transform and parse the xml into an Element class:
  206. # http://docs.python.org/library/xml.etree.elementtree.html
  207. tree = fromstring(self.xml)
  208. for e in tree:
  209. if e.tag not in ('parameter', 'flag'):
  210. self.__setattr__(e.tag, GETFROMTAG[e.tag](e))
  211. #
  212. # extract parameters from the xml
  213. #
  214. self.params_list = [Parameter(p) for p in tree.findall("parameter")]
  215. self.inputs = TypeDict(Parameter)
  216. self.outputs = TypeDict(Parameter)
  217. self.required = []
  218. # Insert parameters into input/output and required
  219. for par in self.params_list:
  220. if par.input:
  221. self.inputs[par.name] = par
  222. else:
  223. self.outputs[par.name] = par
  224. if par.required:
  225. self.required.append(par.name)
  226. #
  227. # extract flags from the xml
  228. #
  229. flags_list = [Flag(f) for f in tree.findall("flag")]
  230. self.flags = TypeDict(Flag)
  231. for flag in flags_list:
  232. self.flags[flag.name] = flag
  233. #
  234. # Add new attributes to the class
  235. #
  236. self.run_ = True
  237. self.finish_ = True
  238. self.env_ = None
  239. self.stdin_ = None
  240. self.stdin = None
  241. self.stdout_ = None
  242. self.stderr_ = None
  243. diz = {'name': 'stdin', 'required': False,
  244. 'multiple': False, 'type': 'all',
  245. 'value': None}
  246. self.inputs['stdin'] = Parameter(diz=diz)
  247. diz['name'] = 'stdout'
  248. self.outputs['stdout'] = Parameter(diz=diz)
  249. diz['name'] = 'stderr'
  250. self.outputs['stderr'] = Parameter(diz=diz)
  251. self.popen = None
  252. self.time = None
  253. if args or kargs:
  254. self.__call__(*args, **kargs)
  255. self.__call__.__func__.__doc__ = self.__doc__
  256. def __call__(self, *args, **kargs):
  257. if not args and not kargs:
  258. self.run()
  259. return
  260. #
  261. # check for extra kargs, set attribute and remove from dictionary
  262. #
  263. if 'flags' in kargs:
  264. for flg in kargs['flags']:
  265. self.flags[flg].value = True
  266. del(kargs['flags'])
  267. if 'run_' in kargs:
  268. self.run_ = kargs['run_']
  269. del(kargs['run_'])
  270. if 'stdin_' in kargs:
  271. self.inputs['stdin'].value = kargs['stdin_']
  272. del(kargs['stdin_'])
  273. if 'stdout_' in kargs:
  274. self.outputs['stdout'].value = kargs['stdout_']
  275. del(kargs['stdout_'])
  276. if 'stderr_' in kargs:
  277. self.outputs['stderr'].value = kargs['stderr_']
  278. del(kargs['stderr_'])
  279. if 'env_' in kargs:
  280. self.env_ = kargs['env_']
  281. del(kargs['env_'])
  282. if 'finish_' in kargs:
  283. self.finish_ = kargs['finish_']
  284. del(kargs['finish_'])
  285. #
  286. # check args
  287. #
  288. for param, arg in zip(self.params_list, args):
  289. param.value = arg
  290. for key, val in kargs.items():
  291. if key in self.inputs:
  292. self.inputs[key].value = val
  293. elif key in self.outputs:
  294. self.outputs[key].value = val
  295. elif key in self.flags:
  296. # we need to add this, because some parameters (overwrite,
  297. # verbose and quiet) work like parameters
  298. self.flags[key].value = val
  299. else:
  300. raise ParameterError('%s is not a valid parameter.' % key)
  301. #
  302. # check if execute
  303. #
  304. if self.run_:
  305. #
  306. # check reqire parameters
  307. #
  308. for k in self.required:
  309. if ((k in self.inputs and self.inputs[k].value is None) or
  310. (k in self.outputs and self.outputs[k].value is None)):
  311. msg = "Required parameter <%s> not set."
  312. raise ParameterError(msg % key)
  313. return self.run()
  314. def get_bash(self):
  315. return ' '.join(self.make_cmd())
  316. def get_python(self):
  317. prefix = self.name.split('.')[0]
  318. name = '_'.join(self.name.split('.')[1:])
  319. params = ', '.join([par.get_python() for par in self.params_list
  320. if par.get_python() != ''])
  321. flags = ''.join([flg.get_python()
  322. for flg in self.flags.values()
  323. if not flg.special and flg.get_python() != ''])
  324. special = ', '.join([flg.get_python()
  325. for flg in self.flags.values()
  326. if flg.special and flg.get_python() != ''])
  327. # pre name par flg special
  328. if flags and special:
  329. return "%s.%s(%s, flags=%r, %s)" % (prefix, name, params,
  330. flags, special)
  331. elif flags:
  332. return "%s.%s(%s, flags=%r)" % (prefix, name, params, flags)
  333. elif special:
  334. return "%s.%s(%s, %s)" % (prefix, name, params, special)
  335. else:
  336. return "%s.%s(%s)" % (prefix, name, params)
  337. def __str__(self):
  338. """!Return the command string that can be executed in a shell
  339. """
  340. return ' '.join(self.make_cmd())
  341. def __repr__(self):
  342. return "Module(%r)" % self.name
  343. @property
  344. def __doc__(self):
  345. """{cmd_name}({cmd_params})
  346. """
  347. head = DOC['head'].format(cmd_name=self.name,
  348. cmd_params=('\n' + # go to a new line
  349. # give space under the function name
  350. (' ' * (len(self.name) + 1))).join([', '.join(
  351. # transform each parameter in string
  352. [str(param) for param in line if param is not None])
  353. # make a list of parameters with only 3 param per line
  354. for line in zip_longest(*[iter(self.params_list)] * 3)]),)
  355. params = '\n'.join([par.__doc__ for par in self.params_list])
  356. flags = self.flags.__doc__
  357. return '\n'.join([head, params, DOC['flag_head'], flags, DOC['foot']])
  358. def get_dict(self):
  359. """!Return a dictionary that includes the name, all valid
  360. inputs, outputs and flags
  361. """
  362. dic = {}
  363. dic['name'] = self.name
  364. dic['inputs'] = [(k, v.value) for k, v in self.inputs.items()
  365. if v.value]
  366. dic['outputs'] = [(k, v.value) for k, v in self.outputs.items()
  367. if v.value]
  368. dic['flags'] = [flg for flg in self.flags if self.flags[flg].value]
  369. return dic
  370. def make_cmd(self):
  371. """!Create the command string that can be executed in a shell
  372. @return The command string
  373. """
  374. skip = ['stdin', 'stdout', 'stderr']
  375. args = [self.name, ]
  376. for key in self.inputs:
  377. if key not in skip and self.inputs[key].value:
  378. args.append(self.inputs[key].get_bash())
  379. for key in self.outputs:
  380. if key not in skip and self.outputs[key].value:
  381. args.append(self.outputs[key].get_bash())
  382. for flg in self.flags:
  383. if self.flags[flg].value:
  384. args.append(str(self.flags[flg]))
  385. return args
  386. def run(self, node=None):
  387. """!Run the module
  388. This function will wait for the process to terminate
  389. in case finish_==True and sets up stdout and stderr.
  390. If finish_==False this function will return after starting
  391. the process. Use self.popen.communicate() of self.popen.wait()
  392. to wait for the process termination. The handling
  393. of stdout and stderr must then be done outside of this
  394. function.
  395. """
  396. if self.inputs['stdin'].value:
  397. self.stdin = self.inputs['stdin'].value
  398. self.stdin_ = PIPE
  399. if self.outputs['stdout'].value:
  400. self.stdout_ = self.outputs['stdout'].value
  401. if self.outputs['stderr'].value:
  402. self.stderr_ = self.outputs['stderr'].value
  403. cmd = self.make_cmd()
  404. start = time.time()
  405. self.popen = Popen(cmd,
  406. stdin=self.stdin_,
  407. stdout=self.stdout_,
  408. stderr=self.stderr_,
  409. env=self.env_)
  410. if self.finish_:
  411. stdout, stderr = self.popen.communicate(input=self.stdin)
  412. self.outputs['stdout'].value = stdout if stdout else ''
  413. self.outputs['stderr'].value = stderr if stderr else ''
  414. self.time = time.time() - start
  415. return self
  416. ###############################################################################
  417. if __name__ == "__main__":
  418. import doctest
  419. doctest.testmod()