Oops, I guess I quoted the text in the wrong way. Sorry for that.
Marcello Gorini wrote: > > Hello Nikos, > > Marcello Gorini. > > Nikos Alexandris wrote: >> >> Hi Marcello, >> >> Marcello Gorini: >>> This is kind of my first e-mail ever to ANY mailing list, so please >>> forgive >>> me if I don't express myself in the best way. I will improve eventually. >>> I am a total beginner in python and GRASS (and any other programming >>> language), so I am trying to run the first examples found in >>> http://grass.osgeo.org/wiki/GRASS_and_Python. >>> I managed to run the shaded_relief one OK, but the one with respect to >>> mapcalc, I couldn't. >>> I copied a script found in >>> http://www.mail-archive.com/grass-...@lists.osgeo.org/msg14669.html >>> (transcribed below) and tried to run it: >>> >>> #!/usr/bin/env python >>> >>> # g.parser demo script for python programing >>> >>> #%module >>> #% description: g.parser test script (python) TUN >>> #%end >>> #%flag >>> #% key: f >>> #% description: A flag >>> #%end >>> #%flag >>> #% key: n >>> #% description: Another flag >>> #%end >>> #%option >>> #% key: raster1 >>> #% type: string >>> #% gisprompt: old,cell,raster >>> #% description: Raster input map >>> #% required : yes >>> #%end >>> #%option >>> #% key: raster2 >>> #% type: string >>> #% gisprompt: old,cell,raster >>> #% description: Raster2 input map >>> #% required : yes >>> #%end >>> #%option >>> #% key: output >>> #% type: string >>> #% gisprompt: new,cell,raster >>> #% description: Output >>> #% required : yes >>> #%end >>> >>> import os >>> import sys >>> >>> import grass.script as grass >>> >>> def main(): >>> output = options['output'] >>> raster1 = options['raster1'] >>> raster2 = options['raster2'] >>> grass.mapcalc("$out = $rast1 + $rast2", out = output, rast1 = >>> raster1, >>> rast2 = raster2) >>> if flags['f']: >>> print "Flag -f set" >>> else: >>> print "Flag -f not set" >>> >>> if flags['n']: >>> print "Flag -n foi definida mas nao devia" >>> else: >>> print "Flag -n not set- GOOD" >>> print options >>> print flags >>> return 0 >>> >>> if __name__ == "__main__": >>> options, flags = grass.parser() >>> sys.exit(main()) >>> >>> >>> It seems to run OK and “100%” is displayed in the GRASS prompt, but when >>> I >>> tried to display the map, I found out that it wasn't created at all. >> >> If there are no errors and the script ran fine, maybe you just need to >> check >> the region settings? >> >>> Just for you to know, to run the script I simply create it (say >>> “myscript.py”), give permissions and run it in a GRASS session with the >>> necessary options: > ./myscript.py output=outmap raster1=map1 >>> raster2=map2 >>> I am using Ubuntu 9.04 64bits and GRASS 6.4.ORC5. >> >> For example: >> >> # check region extent, etc. >> g.region -p >> >> # set region extent and resolution to the map of your interest >> g.region rast=outmap -pa >> >> # give color >> r.colors rast=outmap color=rainbow >> >> # try to display >> d.mon x0 && d.rast outmap >> >> # is it there? >> >> Regards, Nikos >> >> >> No, it is not there. To make sure I set the region to the input map: >> >> g.region rast=map1 -pa >> d.mon x0 && d.rast map1 >> >> OK, it is there. Then I ran the script: >> >> ./myscript.py output=outmap raster1=map1 raster2=map2 >> >> After processing, I tried: >> >> d.rast outmap >> >> But then: >> >> ERROR: Raster map <outmap> not found >> >> What could be wrong? >> >> Best regards, >> >> Marcello. >> >> P.S.: By the way, is this the right way to quote answers? >> >> _______________________________________________ >> grass-user mailing list >> grass-user@lists.osgeo.org >> http://lists.osgeo.org/mailman/listinfo/grass-user >> >> > > -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Trying-python-scripts-tp5109105p5109444.html Sent from the Grass - Users mailing list archive at Nabble.com. _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user