[GRASS-user] map calculator in python loop
Moritz Lennert
mlennert at club.worldonline.be
Wed Dec 12 07:22:23 PST 2018
On 12/12/18 12:59, Alessandro Sebastiani wrote:
> dear all,
> i'm dealing this time with r.mapcalc. I have different rasters which i
> grouped using g.list.
> Now r_list contains all my raster. I want to multiply those rasters to
> another raster (called "costante") in order to obtain a new map, that i
> want to call as r followed by the letter "k"
> Here it is what i've done:
> import grass.script as gscript
> >>> r_list = gscript.read_command("g.list" , type = "rast"
> ).strip().split("\n")
> >>> c = " costante"
> >>> for r in r_list:
> >>> gscript.mapcalc ( " %s = %s * %s " % ((r+"k") , r , c))
Try
" %sk = %s * %s " % (r , r , c) "
Moritz
More information about the grass-user
mailing list