[GRASS-user] Longest flow path

2010-08-03 Thread David Townshend
Hi all Can anyone suggest a way for me to determine the longest flow path in a catchment? I have several hundred to do, so I'd like a way to automate it if possible. Regards David ___ grass-user mailing list grass-user@lists.osgeo.org

Re: [GRASS-user] Re: Installing r.area

2010-08-03 Thread Jarek Jasiewicz
Hamish pisze: fyi I get what it does now, and just fixed a small bug in it which was stopping it from reading CELL maps. Thanks, Hamish, I was offline for over the week so I cannot fix it myself... Hamish If someone is interested in result of r.stats and r.area are the same (as I

Re: [GRASS-user] v.db.connect -p

2010-08-03 Thread Markus Metz
Markus Neteler wrote: Patricio Toledo wrote: Dear GRASS Developers This is a petition: Is it possible to add a option to v.db.connect (like g.region -p) so the output will be somelike:    $ v.db.connect -p map=test layer=2    driver=dbf    database=/the/database    table=thetable    

[GRASS-user] Re: Installing r.area

2010-08-03 Thread LeeDaniel
Alright... I tried the following: l...@pc19384:~ svn checkout https://svn.osgeo.org/grass/grass-addons/raster/r.area/ Ar.area/main.c Ar.area/description.html Ar.area/Makefile Ausgecheckt, Revision 42981. l...@pc19384:~ cd r.area/ l...@pc19384:~/r.area make MODULE_TOPDIR=/opt/grass/

Re: [GRASS-user] Longest flow path

2010-08-03 Thread Jarek Jasiewicz
David Townshend pisze: Hi all Can anyone suggest a way for me to determine the longest flow path in a catchment? I have several hundred to do, so I'd like a way to automate it if possible. Regards David ___ grass-user mailing list

Re: [GRASS-user] Angle between segments of a polyline

2010-08-03 Thread Jarek Jasiewicz
Rainer M Krug pisze: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi a friend of mine wants (needs?) to determine the angle between segments in an ESRI polyline (line vector feature). As there does not seem to be an easy way in ArcGIS 9.2, I thought about askeng here: Is there way of

[GRASS-user] Script from GRASS 6.4 to 7 - g.ask

2010-08-03 Thread Rebecca Bennett
Hello all, I am trying to run a script in GRASS 7 (which was kindly written for me by Micha Silva for GRASS 6.4) which requires the g.ask module that I believe has been removed (http://trac.osgeo.org/grass/wiki/Grass7/NewFeatures). With what should I replace this this command? Thanks in

[GRASS-user] minimun (summed) cost under a polygon

2010-08-03 Thread eklein
Hi, I'm trying to find the location of a polygon over a cost-surface layer (raster) that represents the minimum possible cost when you add the pixels (cost) under it. The cost-surface has an irregular shape and the size of the polygon (in this case a square) is about 1/1000 of the size

Re: [GRASS-user] Longest flow path

2010-08-03 Thread Rich Shepard
On Tue, 3 Aug 2010, Jarek Jasiewicz wrote: r.watershed elevation=elev drainage=dirs streams=stream treshold=some value r.stream.order streams=streams dirs=dirs hack=longest the stream with order 1 will be the longest streams of the catchment Jarek, Please correct me if I'm wrong, but I

Re: [GRASS-user] Longest flow path

2010-08-03 Thread stephen sefick
To add on to Rich's comments. I will gladly suggest things that would be helpful to applied stream ecologists that could be used by a great many folks. I will help in any way. I am finding that I need to write a fair amount of shell script to extract the things that I am interested in. I may

[GRASS-user] Retrieve pixel position

2010-08-03 Thread António Rocha
Greetings I have retriecve a Pixel coordinates but instead of lat lon I want row/col. As far as I can see, r.what only retrieves Lat/lon not row/col. How can I do this? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5337 (20100803

Re: [GRASS-user] Script from GRASS 6.4 to 7 - g.ask

2010-08-03 Thread Markus Neteler
On Tue, Aug 3, 2010 at 2:57 PM, Rebecca Bennett rabenn...@ymail.com wrote: Hello all, I am trying to run a script in GRASS 7 (which was kindly written for me by Micha Silva for GRASS 6.4) which requires the g.ask module that I believe has been removed

Re: [GRASS-user] Re: Installing r.area

2010-08-03 Thread Markus Neteler
On Tue, Aug 3, 2010 at 1:37 PM, LeeDaniel lee.daniel.1...@gmail.com wrote: Alright... I tried the following: l...@pc19384:~ svn checkout https://svn.osgeo.org/grass/grass-addons/raster/r.area/ A    r.area/main.c A    r.area/description.html A    r.area/Makefile Ausgecheckt, Revision 42981.

Re: [GRASS-user] Script from GRASS 6.4 to 7 - g.ask

2010-08-03 Thread Rebecca Bennett
Thank you Markus! From: Markus Neteler nete...@osgeo.org To: Rebecca Bennett rabenn...@ymail.com Cc: GRASS user list grass-user@lists.osgeo.org Sent: Tue, 3 August, 2010 16:36:43 Subject: Re: [GRASS-user] Script from GRASS 6.4 to 7 - g.ask On Tue, Aug 3, 2010

Re: [GRASS-user] Longest flow path

2010-08-03 Thread stephen sefick
I realized this later... I didn't mean to come off as a jerk if I did. I am in everyone's debt who has written or otherwise contributed to all free software projects. Thank you to all of those named and un-named. Stephen On Tue, Aug 3, 2010 at 8:32 AM, stephen sefick ssef...@gmail.com wrote:

Re: [GRASS-user] Longest flow path

2010-08-03 Thread Rich Shepard
On Tue, 3 Aug 2010, Jarek Jasiewicz wrote: yes this is correct answer. I suggest hack ordering not horton or strahler. Hack ordering ascribe first order to the longest stream in the catchement, second order to its tributaries etc. see more in documentation. Jarek, Thanks for the

[GRASS-user] Error while running script that uses g.mapsets

2010-08-03 Thread Helena Herrera
Hi I have a Python Script (for GRASS6.4) where I run this: grass.run_command(g.mapset, mapset = t_mapset) (to change active mapset) But I get all this in the command output: Erasing monitors... PNG: GRASS_TRUECOLOR status: TRUE PNG: collecting to file: map.png, GRASS_WIDTH=768, GRASS_HEIGHT=501

[GRASS-user] Warning from v.digit- Unable to read

2010-08-03 Thread Luisa Peña
I'm running v.digit in Windows and, after I draw my boundaries in TclTk's digitizer, when I close I get this warning in v.digit command output tab: Unable to read vector map But my vectorial is created without problems. Why is tbhis printed? And how can I avoid it? Thanks Luisa

[GRASS-user] Error while creating vectorial file in WinGRASS

2010-08-03 Thread Kim Besson
I have installed WinGRASS-6.4.SVN-r42868-1-Setup.exe in my WinXP machine and, when I select Vector/Develop Vector Map/Create new vector map, I place output name and check Create Attibute table I get an error windows stating: - Error in command execution db.execute Execution failed: 'db.execute

Re: [GRASS-user] Retrieve pixel position

2010-08-03 Thread Nikos Alexandris
(forgot to cc to the list) António: Thanls for the tip. But isn't there any other more interactive way (like r.what?) I thought you have a points vector map and want to retrieve row,col for each point. Maybe d.what.rast -c is what you are looking for then? Cheers, Nikos

[GRASS-user] Warning from v.digit- Unable to read

2010-08-03 Thread Helmut Kudrnovsky
Hi Luisa, I'm running v.digit in Windows and, after I draw my boundaries in TclTk's digitizer, when I close I get this warning in v.digit command output tab: Unable to read vector map But my vectorial is created without problems. Why is tbhis printed? And how can I avoid it? Thanks Luisa just

[GRASS-user] Announcement AGU Session Earth and Space Science Informatics IN10: Open Source Remote Sensing for Environmental Mapping and Analysis

2010-08-03 Thread Karl Benedict
Dear Colleagues, We are pleased to announce the following session for the Fall 2010 American Geophysical Union Meeting December 13-17: Session Earth and Space Science Informatics IN10: Open Source Remote Sensing for Environmental Mapping and Analysis (See below for session description.)

Re: [GRASS-user] v.db.connect -p

2010-08-03 Thread Markus Neteler
On Tue, Aug 3, 2010 at 12:21 PM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Markus Neteler wrote: Patricio Toledo wrote: Dear GRASS Developers This is a petition: Is it possible to add a option to v.db.connect (like g.region -p) so the output will be somelike:    $

[GRASS-user] Distance-weighted focal functions?

2010-08-03 Thread Dylan Beaudette
Hi, After hacking together some R + r.mapcalc code, I started to wonder about how hard it would be to add distance-weighting to r.neighbors... Specifically, the focal distance-weighted mean difference. An example of what I was trying to accomplish can be found here:

[GRASS-user] Feed and manipulate the result of v.distance -pa in python

2010-08-03 Thread Nikos Alexandris
Hi! I need to convert the following bash-one-liner: v.distance --q -pa from=ref_points to=coi_points upload=cat column=anycol | cut -d| -f2 | sort -nu | tail -1 into python. My unsolved problem is that v.distance -pa returns for example thousands of lines in the following format: ... 7839|16

Re: [GRASS-user] v.db.connect -p

2010-08-03 Thread Hamish
Markus N wrote: I have locally backported it to 6.4, too. But pls review the patch first. a number of scripts rely and expect the old style formatting. (have to think about users' personal scripts too, not just our own) it's ok for new flags but please don't touch the formatting of old flags