Hello Marcello. It's good to _finally_ see you here :)
I tried your script and it worked for me, but I had to change the last line from sys.exit(main()) to: main() could you try that? best Carlos On Fri, May 28, 2010 at 08:23, Marcello Gorini <[email protected]> wrote: > > > Nikos Alexandris wrote: >> >> Marcello Gorini wrote: >>> Thanks for the try. I tested it again with other rasters, but I got the >>> same error. Could it possibly be something related to the 64bits >>> architecture? >> >> I don't think so. I work with Kubuntu Lucid 64 bit :-) >> >>> Or could GRASS be sending the outcome somewhere else? >> >> Negative. You are running the script from within a grass-session (as you >> wrote >> in the first post). Let's try to take this step by step. What is the >> result of >> each of the following commands: >> >> ls -l myscript.py >> g.gisenv >> g.region -p >> r.info map1 >> r.info map2 >> >>> Well, I guess that eventually I will go around this problem or even bump >>> into the solution, but the most important thing to say is that I >>> appreciated the effort and that it was a very good reception to this >>> list. >> >> Bah... it was only me trying to not forget the basics with python scripts. >> You will get better and more precise replies from other advanced users and >> developers. >> >>> Thank you very much Nikos. And thank you also for the explanations with >>> respect to the right way to answer posts. I will do my best to make my >>> posts as clean and explanatory as possible. >> >> You are warned, you might get yourself into trouble ;-p >> Cheers, Nikos >> _______________________________________________ >> grass-user mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/grass-user >> >> > > OK, thank you, so let's go step by step: > ls -l myscript.py > result: > -rwxr-xr-x 1 marcello marcello 1227 2010-05-28 08:13 myscript.py > > g.gisenv > result: > GISDBASE=/home/marcello/grassdata > LOCATION_NAME=Ocean_floor > MAPSET=level4 > MONITOR=x0 > GRASS_GUI=tcltk > > g.region -p > result: > projection: 99 (Mercator) > zone: 0 > datum: wgs84 > ellipsoid: wgs84 > north: 3947900 > south: 3633400 > west: -7366700 > east: -7022600 > nsres: 3700 > ewres: 3700 > rows: 85 > cols: 93 > cells: 7905 > > r.info map1 > result: > +----------------------------------------------------------------------------+ > | Layer: map1 Date: Thu May 27 18:48:30 2010 > | > | Mapset: level4 Login of Creator: marcello > | > | Location: Ocean_floor > | > | DataBase: /home/marcello/grassdata > | > | Title: ( topmod ) > | > | Timestamp: none > | > > |----------------------------------------------------------------------------| > | > | > | Type of Map: raster Number of Categories: 255 > | > | Data Type: FCELL > | > | Rows: 1177 > | > | Columns: 1388 > | > | Total Cells: 1633676 > | > | Projection: Mercator > | > | N: 6358450 S: 2003550 Res: 3700 > | > | E: -3920150 W: -9055750 Res: 3700 > | > | Range of data: min = 0.000000 max = 10000.000000 > | > | > | > | Data Description: > | > | generated by r.mapcalc > | > | > | > | Comments: > | > | if(topmod > 0, 0, -topmod) > | > | > | > > +----------------------------------------------------------------------------+ > > r.info map2 > result: > +----------------------------------------------------------------------------+ > | Layer: map2 Date: Thu May 27 18:49:29 2010 > | > | Mapset: level4 Login of Creator: marcello > | > | Location: Ocean_floor > | > | DataBase: /home/marcello/grassdata > | > | Title: ( topmod2 ) > | > | Timestamp: none > | > > |----------------------------------------------------------------------------| > | > | > | Type of Map: raster Number of Categories: 255 > | > | Data Type: FCELL > | > | Rows: 1177 > | > | Columns: 1388 > | > | Total Cells: 1633676 > | > | Projection: Mercator > | > | N: 6358450 S: 2003550 Res: 3700 > | > | E: -3920150 W: -9055750 Res: 3700 > | > | Range of data: min = 1000.000000 max = 1000.000000 > | > | > | > | Data Description: > | > | generated by r.mapcalc > | > | > | > | Comments: > | > | if(topmod > -2, 1000, topmod) > | > | > | > > +----------------------------------------------------------------------------+ > > Now just to emphasize :) > Result of "r.info outmap": > ERROR: Raster map <outmap> not found > > Can you make anything out of it? > > Best regards, > Marcello. > -- > View this message in context: > http://osgeo-org.1803224.n2.nabble.com/Trying-python-scripts-tp5109105p5112658.html > Sent from the Grass - Users mailing list archive at Nabble.com. > _______________________________________________ > grass-user mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user > -- Prof. Carlos Henrique Grohmann - Geologist D.Sc. Institute of Geosciences - Univ. of São Paulo, Brazil http://www.igc.usp.br/pessoais/guano http://lattes.cnpq.br/5846052449613692 Linux User #89721 ________________ Can’t stop the signal. _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
