[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