[GRASS-user] Re: announcement of strahler stream order analysis

GRASS GRASS user
1

Hi Jarek,
I downloaded and compiled r.stream without problems but i'm afraid it does not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps...
*** glibc detected *** r.stream: corrupted double-linked list: 0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 ---p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 ---p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 ---p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 ---p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r--p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 ---p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 ---p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 ---p 00008000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 ---p 00053000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r--p 00000000 08:03 6521179 /usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r--p 00000000 08:03 6521188 /usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r--s 00000000 08:03 6496774 /usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one dataset.

To compile it copy r.streams directory to your grassXX../raster/ cd
r.stream run make and make install

What module calculate:
- strahler stream order
- shreeve stream magnitude
- class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream network.

Analysis of r.stream shall be run on same regions setings as r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check if
region and rasters are identical. It is planned to add in nearest feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons (or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I don't
know C.
General code is based on GRASS r.topindex module (by Keith Beven and Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya, modified)
Shreeve stream magnitude and class of topological network are my own ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL: http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106

2

It works fine on my machine. Mac OS X, compilied addon and compilied
GRASS from scratch. The strahler output also looks sensable. I do
have a question, though. Say I have my own delineated streams with
another method, but in raster format- will this still work? And is
there an easy way to a vector out of a raster stream network?

Stephen Sefick

On Tue, Jul 28, 2009 at 4:32 AM, Margherita Di Leo<diregola@gmail.com> wrote:

Hi Jarek,
I downloaded and compiled r.stream without problems but i'm afraid it does
not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler
shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps...
*** glibc detected *** r.stream: corrupted double-linked list:
0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 ---p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 ---p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 ---p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 ---p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r--p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 ---p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 ---p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 ---p 00008000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 ---p 00053000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r--p 00000000 08:03 6521179
/usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r--p 00000000 08:03 6521188
/usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r--s 00000000 08:03 6496774
/usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is
regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code
you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one
dataset.

To compile it copy r.streams directory to your grassXX../raster/ cd
r.stream run make and make install

What module calculate:
- strahler stream order
- shreeve stream magnitude
- class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream
network.

Analysis of r.stream shall be run on same regions setings as
r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check
if
region and rasters are identical. It is planned to add in nearest
feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and
have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more
than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton
analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is
directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons
(or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I
don't
know C.
General code is based on GRASS r.topindex module (by Keith Beven and
Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya,
modified)
Shreeve stream magnitude and class of topological network are my own
ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL:
http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106

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

--
Stephen Sefick

Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.

                -K. Mullis

3

I tried it just right now to see the strahler ordering and… it works greatly!
well done Jarek…:slight_smile:
I think I’m going to read the code because I’ve got to understand how did you solved the problem of cleaning topology…
congratulations! great work!!!

Annalisa

2009/7/28 stephen sefick <ssefick@gmail.com>

It works fine on my machine. Mac OS X, compilied addon and compilied
GRASS from scratch. The strahler output also looks sensable. I do
have a question, though. Say I have my own delineated streams with
another method, but in raster format- will this still work? And is
there an easy way to a vector out of a raster stream network?

Stephen Sefick

On Tue, Jul 28, 2009 at 4:32 AM, Margherita Di Leo<diregola@gmail.com> wrote:

Hi Jarek,
I downloaded and compiled r.stream without problems but i’m afraid it does
not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler
shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps…
*** glibc detected *** r.stream: corrupted double-linked list:
0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 —p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 —p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 —p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 —p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r–p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 —p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 —p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 —p 00008000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 —p 00053000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r–p 00000000 08:03 6521179
/usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r–p 00000000 08:03 6521188
/usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r–s 00000000 08:03 6496774
/usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is
regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code
you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one
dataset.

To compile it copy r.streams directory to your grassXX…/raster/ cd
r.stream run make and make install

What module calculate:

  • strahler stream order
  • shreeve stream magnitude
  • class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream
network.

Analysis of r.stream shall be run on same regions setings as
r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check
if
region and rasters are identical. It is planned to add in nearest
feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and
have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more
than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton
analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is
directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons
(or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I
don’t
know C.
General code is based on GRASS r.topindex module (by Keith Beven and
Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya,
modified)
Shreeve stream magnitude and class of topological network are my own
ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute

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

–
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL:
http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106

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

–
Stephen Sefick

Let’s not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.

-K. Mullis

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

4

2009/7/28 Annalisa Minelli <annagrass6@gmail.com>

I tried it just right now to see the strahler ordering and… it works greatly!
well done Jarek…:slight_smile:

I think I’m going to read the code because I’ve got to understand how did you solved the problem of cleaning topology…

oops… forgive me… you didn’t use v.strahler so you don’t pass through vector format…
good idea!

congratulations! great work!!!

Annalisa

2009/7/28 stephen sefick <ssefick@gmail.com>

It works fine on my machine. Mac OS X, compilied addon and compilied
GRASS from scratch. The strahler output also looks sensable. I do
have a question, though. Say I have my own delineated streams with
another method, but in raster format- will this still work? And is
there an easy way to a vector out of a raster stream network?

Stephen Sefick

On Tue, Jul 28, 2009 at 4:32 AM, Margherita Di Leo<diregola@gmail.com> wrote:

Hi Jarek,
I downloaded and compiled r.stream without problems but i’m afraid it does
not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler
shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps…
*** glibc detected *** r.stream: corrupted double-linked list:
0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752
/usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 —p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 —p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 —p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 —p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r–p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 —p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 —p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242
/usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 —p 00008000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193
/usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 —p 00053000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236
/usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r–p 00000000 08:03 6521179
/usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r–p 00000000 08:03 6521188
/usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r–s 00000000 08:03 6496774
/usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is
regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code
you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one
dataset.

To compile it copy r.streams directory to your grassXX…/raster/ cd
r.stream run make and make install

What module calculate:

  • strahler stream order
  • shreeve stream magnitude
  • class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream
network.

Analysis of r.stream shall be run on same regions setings as
r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check
if
region and rasters are identical. It is planned to add in nearest
feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and
have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more
than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton
analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is
directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons
(or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I
don’t
know C.
General code is based on GRASS r.topindex module (by Keith Beven and
Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya,
modified)
Shreeve stream magnitude and class of topological network are my own
ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute

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

–
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL:
http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106

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

–
Stephen Sefick

Let’s not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.

-K. Mullis

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

5

Hi Jarek, Hi list

I downloaded and compiled the latest version, commented lines 87 and 88 in io.c file

/* G_free(rDirs);*/
/* G_free(rStreams);*/

it does not abort now, but it does not go after this:

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 strahler=strahler shreeve=shreeve class=class
Reading maps...
Finding nodes...
Creating springs index...

I tried first in a big region, then in the spearfish. I thought .. can it be a 64 bit problem? My system is Ubuntu 8.04:

uname -a
Linux marghe 2.6.24-24-generic #1 SMP Fri Jul 24 22:15:50 UTC 2009 x86_64 GNU/Linux

Did anyone try r.stream on a 64 bit architecture?

Any suggestion or hint appreciate

Cheers

margherita

Jarek Jasiewicz wrote:

thanks,
Please download and compile newer wersion
it's look like alocation memory error
I have some reports that program work but I will try to find potential buc
Jarek

Margherita Di Leo pisze:

Hi Jarek,
I downloaded and compiled r.stream without problems but i'm afraid it does not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps...
*** glibc detected *** r.stream: corrupted double-linked list: 0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 ---p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 ---p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 ---p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 ---p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r--p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 ---p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 ---p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 ---p 00008000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 ---p 00053000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r--p 00000000 08:03 6521179 /usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r--p 00000000 08:03 6521188 /usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r--s 00000000 08:03 6496774 /usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one dataset.

To compile it copy r.streams directory to your grassXX../raster/ cd
r.stream run make and make install

What module calculate:
- strahler stream order
- shreeve stream magnitude
- class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream network.

Analysis of r.stream shall be run on same regions setings as r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check if
region and rasters are identical. It is planned to add in nearest feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons (or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I don't
know C.
General code is based on GRASS r.topindex module (by Keith Beven and Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya, modified)
Shreeve stream magnitude and class of topological network are my own ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL: http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106

6

How much RAM you have?

try to calculate ony one output map
if p[rogram do not responce after 3 munutes brak it.

Jarek

Margherita Di Leo pisze:

Hi Jarek, Hi list

I downloaded and compiled the latest version, commented lines 87 and 88 in io.c file

/* G_free(rDirs);*/
/* G_free(rStreams);*/

it does not abort now, but it does not go after this:

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 strahler=strahler shreeve=shreeve class=class
Reading maps...
Finding nodes...
Creating springs index...

I tried first in a big region, then in the spearfish. I thought .. can it be a 64 bit problem? My system is Ubuntu 8.04:

uname -a
Linux marghe 2.6.24-24-generic #1 SMP Fri Jul 24 22:15:50 UTC 2009 x86_64 GNU/Linux

Did anyone try r.stream on a 64 bit architecture?

Any suggestion or hint appreciate

Cheers

margherita

Jarek Jasiewicz wrote:

thanks,
Please download and compile newer wersion
it's look like alocation memory error
I have some reports that program work but I will try to find potential buc
Jarek

Margherita Di Leo pisze:

Hi Jarek,
I downloaded and compiled r.stream without problems but i'm afraid it does not work:

r.stream dir=drainage@ex_20 stream=basin1@ex_20 strahler=strahler shreeve=shreve class=class

dirs, 3
streams, 4
Reading maps...
*** glibc detected *** r.stream: corrupted double-linked list: 0x000000000067c2c0 ***
======= Backtrace: =========
/lib/libc.so.6[0x7f77dfb86322]
/lib/libc.so.6(cfree+0x8c)[0x7f77dfb89c1c]
r.stream(createMaps+0x353)[0x401df3]
r.stream(main+0x1c1)[0x402011]
/lib/libc.so.6(__libc_start_main+0xf4)[0x7f77dfb301c4]
r.stream[0x401619]
======= Memory map: ========
00400000-00404000 r-xp 00000000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00603000-00604000 rw-p 00003000 08:03 5955752 /usr/local/grass-6.5.svn/bin/r.stream
00604000-41e9f000 rw-p 00604000 00:00 0 [heap]
7f77d8000000-7f77d8021000 rw-p 7f77d8000000 00:00 0
7f77d8021000-7f77dc000000 ---p 7f77d8021000 00:00 0
7f77df700000-7f77df70d000 r-xp 00000000 08:03 4071544 /lib/libgcc_s.so.1
7f77df70d000-7f77df90d000 ---p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90d000-7f77df90e000 rw-p 0000d000 08:03 4071544 /lib/libgcc_s.so.1
7f77df90e000-7f77df910000 r-xp 00000000 08:03 6677083 /lib/libdl-2.7.so
7f77df910000-7f77dfb10000 ---p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb10000-7f77dfb12000 rw-p 00002000 08:03 6677083 /lib/libdl-2.7.so
7f77dfb12000-7f77dfc6a000 r-xp 00000000 08:03 6677080 /lib/libc-2.7.so
7f77dfc6a000-7f77dfe6a000 ---p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6a000-7f77dfe6d000 r--p 00158000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6d000-7f77dfe6f000 rw-p 0015b000 08:03 6677080 /lib/libc-2.7.so
7f77dfe6f000-7f77dfe74000 rw-p 7f77dfe6f000 00:00 0
7f77dfe74000-7f77dfef4000 r-xp 00000000 08:03 6677084 /lib/libm-2.7.so
7f77dfef4000-7f77e00f3000 ---p 00080000 08:03 6677084 /lib/libm-2.7.so
7f77e00f3000-7f77e00f5000 rw-p 0007f000 08:03 6677084 /lib/libm-2.7.so
7f77e00f5000-7f77e010b000 r-xp 00000000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e010b000-7f77e030b000 ---p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030b000-7f77e030c000 rw-p 00016000 08:03 6490242 /usr/lib/libz.so.1.2.3.3
7f77e030c000-7f77e0314000 r-xp 00000000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0314000-7f77e0513000 ---p 00008000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0513000-7f77e0514000 rw-p 00007000 08:03 984193 /usr/local/grass-6.5.svn/lib/libgrass_datetime.6.5.svn.so
7f77e0514000-7f77e0567000 r-xp 00000000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0567000-7f77e0766000 ---p 00053000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0766000-7f77e0769000 rw-p 00052000 08:03 984236 /usr/local/grass-6.5.svn/lib/libgrass_gis.6.5.svn.so
7f77e0769000-7f77e0770000 rw-p 7f77e0769000 00:00 0
7f77e0770000-7f77e078d000 r-xp 00000000 08:03 6676548 /lib/ld-2.7.so
7f77e0926000-7f77e0965000 r--p 00000000 08:03 6521179 /usr/lib/locale/en_GB.utf8/LC_CTYPE
7f77e0965000-7f77e0969000 rw-p 7f77e0965000 00:00 0
7f77e0982000-7f77e0983000 r--p 00000000 08:03 6521188 /usr/lib/locale/en_GB.utf8/LC_MESSAGES/SYS_LC_MESSAGES
7f77e0983000-7f77e098a000 r--s 00000000 08:03 6496774 /usr/lib/gconv/gconv-modules.cache
7f77e098a000-7f77e098d000 rw-p 7f77e098a000 00:00 0
7f77e098d000-7f77e098f000 rw-p 0001d000 08:03 6676548 /lib/ld-2.7.so
7fffe8979000-7fffe898e000 rw-p 7ffffffea000 00:00 0 [stack]
7fffe89fe000-7fffe8a00000 r-xp 7fffe89fe000 00:00 0 [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall]
Aborted

Hi list!

I just finish first version of new grass module r.stream. Module is regular
Grass module, written in C an has no addational dependencies. It was
compiled with GRASS65 devel downloaded in last two week. Source code you
can find here:

http://heretic.livenet.pl/heretic/r.stream.tar.gz

Module compiles and seems to work. Analysis was tested only on one dataset.

To compile it copy r.streams directory to your grassXX../raster/ cd
r.stream run make and make install

What module calculate:
- strahler stream order
- shreeve stream magnitude
- class of topological elements of drainage networks

Module use r.watershed outputs: SDF flow direction and SDF stream network.

Analysis of r.stream shall be run on same regions setings as r.watershed. If
no, run g.region rast=your_dir_file first. For now module do not check if
region and rasters are identical. It is planned to add in nearest feature.
It also do not check if dirs and streams are SDF result. Threre are no
addational description. It will be added after testing.

It works only on cells data. Floating points data are not allowed (and have
no sense).

Module is quite fast. 3000 x3000 cells all three analysis runs no more than
30 s. on AMD 3800 2GB, Ubuntu 8.04

Module is part of comprehensive project which will add full Horton analysis
and some addational terrrains analysis to Grass. More details in the
feature. Any cooperation is welcome. The next module I plan is directional
vectorisation of stream network.

For now module is available only form web address pointed below. After
testing and removing critical errors it will be added to grass add-ons (or
maybe to main branch?).

Any coments, improvements and critics is very welcome.

More advanced Grass developers: Please check my code. Its funny, but I don't
know C.
General code is based on GRASS r.topindex module (by Keith Beven and Huidae
Cho, strongly modified)
Strahler stream order algorithm from Saga GIS (by Victor Olaya, modified)
Shreeve stream magnitude and class of topological network are my own ideas

Test and enjoy.

Jarek Jasiewicz
Adam Mickiewicz University, Poznan
Geoecology and Geoinformation Institute
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

7

Jarek Jasiewicz wrote:

How much RAM you have?

4 GB

try to calculate ony one output map
if p[rogram do not responce after 3 munutes brak it.

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 strahler=strahler
Reading maps...
Finding nodes...
Creating springs index...

- too much time -

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 shreeve=shreeve
Reading maps...
Finding nodes...
Creating springs index...

- too much time -

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 class=class
Reading maps...
Finding nodes...
Creating springs index...
Writing maps...
100%
class OK

margherita

Jarek

Margherita Di Leo pisze:

Hi Jarek, Hi list

I downloaded and compiled the latest version, commented lines 87 and 88 in io.c file

/* G_free(rDirs);*/
/* G_free(rStreams);*/

it does not abort now, but it does not go after this:

GRASS 6.5.svn (spearfish60):~ > r.stream dir=drainage_SFD@user1 stream=stream_SFD@user1 strahler=strahler shreeve=shreeve class=class
Reading maps...
Finding nodes...
Creating springs index...

I tried first in a big region, then in the spearfish. I thought .. can it be a 64 bit problem? My system is Ubuntu 8.04:

uname -a
Linux marghe 2.6.24-24-generic #1 SMP Fri Jul 24 22:15:50 UTC 2009 x86_64 GNU/Linux

Did anyone try r.stream on a 64 bit architecture?

Any suggestion or hint appreciate

Cheers

margherita

--
Eng. Margherita Di Leo
Ph.D. Student Methods and Technologies for Environmental Monitoring
Department of Environmental Engineering and Physics (DIFA)

University of Basilicata Campus Macchia Romana
85100 - Potenza Italy

Office: +39-0971205363
Fax: +39-0971205160
E-mail: dileomargherita AT gmail DOT com
Skype: dileomargherita
URL: http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePage&id=106