[GRASS-user] Calculate stream order from digitized hydrology network

filipesdias · May 14, 2012, 11:26am

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

Johannes_Radinger1 · May 14, 2012, 1:47pm

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don't know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

filipesdias · May 14, 2012, 1:53pm

Hi Johannes

Yes, its in vector format, but it could be converted to raster if necessary.

I installed v.strahler with g.extension, but loading it produced an error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at> wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

filipesdias · May 15, 2012, 9:41am

I found a similar question here:

http://gis.stackexchange.com/questions/22735/how-to-calculate-stream-order-for-vector-data-without-a-dem

On Mon, May 14, 2012 at 2:53 PM, Filipe Silva Dias <filipesdias@gmail.com> wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if necessary.

I installed v.strahler with g.extension, but loading it produced an error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at> wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

madi · May 15, 2012, 10:39am

Hi Filipe,

could you fill a bug report reporting this? Thanks

Best,
madi

On Mon, May 14, 2012 at 3:53 PM, Filipe Silva Dias <filipesdias@gmail.com> wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if necessary.

I installed v.strahler with g.extension, but loading it produced an error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at> wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

–
Ing. Margherita Di Leo, Ph.D.

filipesdias · May 15, 2012, 11:14am

Sure, is there a bug tracker for grass gis extensions…?

On Tue, May 15, 2012 at 11:39 AM, Margherita Di Leo <diregola@gmail.com> wrote:

Hi Filipe,

could you fill a bug report reporting this? Thanks

Best,
madi

On Mon, May 14, 2012 at 3:53 PM, Filipe Silva Dias <filipesdias@gmail.com> wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if necessary.

I installed v.strahler with g.extension, but loading it produced an error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at> wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

–
Ing. Margherita Di Leo, Ph.D.

Daniel_Victoria · May 15, 2012, 11:36am

There is a software from the people at Rio de Janeiro State University
that should do what you want. It's called HydroFlow

http://www.fgel.uerj.br/labgis/hydroflow/en/index_ingles.htm

I used it once and it's quite simple to run. Results are also very
good but you can't have loops in your drainadge.

Better yet, it's open source

Cheers
Daniel

On Tue, May 15, 2012 at 8:14 AM, Filipe Silva Dias
<filipesdias@gmail.com> wrote:

Sure, is there a bug tracker for grass gis extensions...?

On Tue, May 15, 2012 at 11:39 AM, Margherita Di Leo <diregola@gmail.com>
wrote:

Hi Filipe,

could you fill a bug report reporting this? Thanks

Best,
madi

On Mon, May 14, 2012 at 3:53 PM, Filipe Silva Dias <filipesdias@gmail.com>
wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if
necessary.

I installed v.strahler with g.extension, but loading it produced an
error. Is there another one?

Thanks

Traceback (most recent call last):
File "C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py", line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File "C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py", line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File "C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py", line 1732, in ParseCommand

blackList = _blackList)
File "C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py", line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
"C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py",
line 1283, in XML

return parser.close()
File
"C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py",
line 1636, in close

self._raiseerror(v)
File
"C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py",
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at>
wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don't know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

> Dear all,
>
> I have a digitized hydrological network and I need to calculate the
> stream order (Horton, Strahler etc).
>
> I found r.stream.order but it requires that the hydrological network
> comes from a DEM, which in my case produced an hydrological network that is
> considerably different than the digitized one (tried both with r.wateshed,
> r.terraflow and r.stream.extract).
>
> Is there a tool capable of doing this? Apparently ArcGIS can do this,
> but id prefer na open source solution.
>
> Thanks
>
> Best regards
>
> Filipe
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Ing. Margherita Di Leo, Ph.D.

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

filipesdias · May 15, 2012, 12:11pm

Thanks Daniel
Easy and simple and it worked perfectly =)
Best

On Tue, May 15, 2012 at 12:36 PM, Daniel Victoria <daniel.victoria@gmail.com> wrote:

There is a software from the people at Rio de Janeiro State University
that should do what you want. It’s called HydroFlow

http://www.fgel.uerj.br/labgis/hydroflow/en/index_ingles.htm

I used it once and it’s quite simple to run. Results are also very
good but you can’t have loops in your drainadge.

Better yet, it’s open source

Cheers
Daniel

On Tue, May 15, 2012 at 8:14 AM, Filipe Silva Dias

<filipesdias@gmail.com> wrote:

Sure, is there a bug tracker for grass gis extensions…?

On Tue, May 15, 2012 at 11:39 AM, Margherita Di Leo <diregola@gmail.com>
wrote:

Hi Filipe,

could you fill a bug report reporting this? Thanks

Best,
madi

On Mon, May 14, 2012 at 3:53 PM, Filipe Silva Dias <filipesdias@gmail.com>
wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if
necessary.

I installed v.strahler with g.extension, but loading it produced an
error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at>
wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the
stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network
comes from a DEM, which in my case produced an hydrological network that is
considerably different than the digitized one (tried both with r.wateshed,
r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this,
but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

–
Ing. Margherita Di Leo, Ph.D.

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

madi · May 15, 2012, 12:55pm

It’s the same bug tracker of main code: http://trac.osgeo.org/grass/
You can just indicate the keyword “add-ons”.
Thank you for collaboration.

madi

On Tue, May 15, 2012 at 1:14 PM, Filipe Silva Dias <filipesdias@gmail.com> wrote:

Sure, is there a bug tracker for grass gis extensions…?

On Tue, May 15, 2012 at 11:39 AM, Margherita Di Leo <diregola@gmail.com> wrote:

Hi Filipe,

could you fill a bug report reporting this? Thanks

Best,
madi

On Mon, May 14, 2012 at 3:53 PM, Filipe Silva Dias <filipesdias@gmail.com> wrote:

Hi Johannes

Yes, its in vector format, but it could be converted to raster if necessary.

I installed v.strahler with g.extension, but loading it produced an error. Is there another one?

Thanks

Traceback (most recent call last):
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\prompt.py”, line 1105, in OnKeyPressed

self.parent.RunCmd(cmd)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\goutput.py”, line 523, in RunCmd

task = menuform.GUI(show = None).ParseCommand(command)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\wxpython\gui_m
odules\menuform.py”, line 1732, in ParseCommand

blackList = _blackList)
File “C:\OSGeo4W\apps\grass\grass-6.4.2\etc\python\grass\s
cript\task.py”, line 460, in parse_interface

tree = etree.fromstring(get_interface_description(name))
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1283, in XML

return parser.close()
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1636, in close

self._raiseerror(v)
File
“C:\OSGeo4W\apps\Python27\lib\xml\etree\ElementTree.py”,
line 1488, in _raiseerror

raise err
xml.etree.ElementTree
.
ParseError
:
no element found: line 1, column 0

On Mon, May 14, 2012 at 2:47 PM, Johannes Radinger <jradinger@gmx.at> wrote:

Hi Filipe,

you mean your hydrological network is in vector format?

Maybe you can try the v.strahler add-on but I don’t know
how uptodate this module is:

http://grass.osgeo.org/wiki/GRASS_AddOns#v.strahler

/johannes

Am 14.05.2012 um 13:26 schrieb Filipe Silva Dias:

Dear all,

I have a digitized hydrological network and I need to calculate the stream order (Horton, Strahler etc).

I found r.stream.order but it requires that the hydrological network comes from a DEM, which in my case produced an hydrological network that is considerably different than the digitized one (tried both with r.wateshed, r.terraflow and r.stream.extract).

Is there a tool capable of doing this? Apparently ArcGIS can do this, but id prefer na open source solution.

Thanks

Best regards

Filipe

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

–
Ing. Margherita Di Leo, Ph.D.

–
Ing. Margherita Di Leo, Ph.D.

neteler · May 15, 2012, 9:15pm

On Tue, May 15, 2012 at 2:55 PM, Margherita Di Leo <diregola@gmail.com> wrote:

It's the same bug tracker of main code: http://trac.osgeo.org/grass/
You can just indicate the keyword "add-ons".
Thank you for collaboration.

I have now added a new "Addons" component in trac.

Markus

filipesdias · May 15, 2012, 9:19pm

Done

http://trac.osgeo.org/grass/ticket/1657

On Tue, May 15, 2012 at 10:15 PM, Markus Neteler <neteler@osgeo.org> wrote:

On Tue, May 15, 2012 at 2:55 PM, Margherita Di Leo <diregola@gmail.com> wrote:

It’s the same bug tracker of main code: http://trac.osgeo.org/grass/
You can just indicate the keyword “add-ons”.
Thank you for collaboration.

I have now added a new “Addons” component in trac.

Markus