[GRASS-user] Strange issue with v.buffer and latest 6.4

2010-08-17 Thread thedok78
Hello everybody, I have a strange issue with the v.buffer command: Here you can see 2 vectors: red is input and grey is the output of v.buffer: http://dl.dropbox.com/u/348430/buffer.png The upper left of the polygon it's not correct Here you can download the GML file (EPSG:4326):

[GRASS-user] Two questions about r.stream.*

2010-05-10 Thread thedok78
Hello everyone, before beginning to write script from scratch I'd like to know if I can use the r.stream.* for my tasks. What I need to this is: - For each subbasin in a certain area extract each subbasin-outlet; - Input a certain number of points and extract the subbasins using the points as

[GRASS-user] Re: r.mapcalc and nested if

2010-04-16 Thread thedok78
Thank you all for yours answers! I'll explain my problem: I have got a raster map with a lot of small area like this one: 112 112 112 122 112 112 112 1 1 122 122 1 1 1 122 112 122 122 112 112 The boundary is made cells with values from 100 to 200 instead the inner ones have

[GRASS-user] r.mapcalc and nested if

2010-04-14 Thread thedok78
Hello everybody, I need to run a complex r.mapcalc expression, but I can't understand how to use nested if. This is the expression: r.mapcalc

[GRASS-user] Re: r.mapcalc and nested if

2010-04-14 Thread thedok78
Sorry, this is the expression: r.mapcalc 'map2=if(map ==1 ,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))' Thanks Luca -- View this message in context:

[GRASS-user] Re: Parsing output with python

2010-04-08 Thread thedok78
Thank you very much, it worked fine! Ciao Luca -- View this message in context: http://n2.nabble.com/Parsing-output-with-python-tp4864622p4870862.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list

[GRASS-user] Output parsing in python

2010-04-07 Thread thedok78
Hello everybody, I am porting a bash script to python and for it I need to know the number of cells of a certain category. Right now I use this syntax: grass.parse_command('r.stats',flags='c',input='map'), but the output is something like this: {'112 525': None.. }. Is there a way to have

[GRASS-user] Parsing output with python

2010-04-07 Thread thedok78
Hello everybody, I am porting a bash script to python and for it I need to know the number of cells of a certain category. Right now I use this syntax: grass.parse_command('r.stats',flags='c',input='map'), but the output is something like this: {'112 525': None.. }. Is there a way to have

[GRASS-user] Re: r.out.gdal and INTERLEAVE=PIXEL

2010-02-23 Thread thedok78
Thanks for all your answer. This is the reply from the Gdal list: http://trac.osgeo.org/gdal/ticket/3433 My problem is how to export a raster to be used inside ArcGIS..the procedure on the wiki [1] uses INTERLEAVE=PIXEL . If I use a standard r.out.gdal syntax, the min shown is 2.22 *e^-308 and

[GRASS-user] r.out.gdal and INTERLEAVE=PIXEL

2010-02-22 Thread thedok78
Hello everybody, I need to export some raster data from GRASS to ArcView; as suggested on the wiki [1] I used the option INTERLEAVE=PIXEL with r.out.gdal: r.out.gdal in=elevation.10m format=GTiff type=Float32 createopt=INTERLEAVE=PIXEL out=output.tif The problem is that, using gdalinfo,