Carlos Grohmann-2 wrote: > > 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 <gor...@gmail.com> 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 >>> grass-user@lists.osgeo.org >>> 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 >> grass-user@lists.osgeo.org >> 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 > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > >
Hey Carlos, how are you my friend? I used to understand only 10% of the posts, but I raised it to 20%, so I took the courage to come in :) Thanks for the hint, but I tried it and it still doesn't work. That is a very weird problem! Cheers, Marcello. -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Trying-python-scripts-tp5109105p5113324.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