Re: [GRASS-user] Re: Building topologically correct stream network
Dylan Beaudette wrote: On Tuesday 21 July 2009, Jarek Jasiewicz wrote: BTW I'm about to finish regular grass modue for calculating Strahler Stream order and Shreeve Stream Magnitude. I hope I publish it for testing in next few days. Any help appreciated. I would be happy to help test! Cheers, Dylan +1 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=visualizzaHomePageid=106 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Building topologically correct stream network
Hi Stephen, I actually am fighting against almost the same problem, so let's do it together :) I found that r.strahler.sh works well to delineate streams and correct topology, (I use a treshold of 1km^2, it works well but don't ask me why :) ), but (in my little experience) does not work to order streams, so I must run v.strahler after r.strahler.sh. I know it's a very small hint but hope this helps. Best margherita I have a low relief area in the costal plain of georgia that I am trying to produce a topologically correct stream network. There are many breaks in the flow accumulation grid and I would like to delineate the streams and correct the topology after r.thin and r.to.vect have been preformed so that I can then use v.strahler to order the streams. Thanks for any help and I will provide and data that would be helpful. thanks, -- 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 -- 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=visualizzaHomePageid=106 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Building topologically correct stream network
Hi maybe this will help: http://gis.vsb.cz/GIS_Ostrava/GIS_Ova_2009/sbornik/Lists/Papers/060.pdf it require another gis program but gives great results! Jarek Margherita Di Leo pisze: Hi Stephen, I actually am fighting against almost the same problem, so let's do it together :) I found that r.strahler.sh works well to delineate streams and correct topology, (I use a treshold of 1km^2, it works well but don't ask me why :) ), but (in my little experience) does not work to order streams, so I must run v.strahler after r.strahler.sh. I know it's a very small hint but hope this helps. Best margherita I have a low relief area in the costal plain of georgia that I am trying to produce a topologically correct stream network. There are many breaks in the flow accumulation grid and I would like to delineate the streams and correct the topology after r.thin and r.to.vect have been preformed so that I can then use v.strahler to order the streams. Thanks for any help and I will provide and data that would be helpful. thanks, -- 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
Re: [GRASS-user] Re: Building topologically correct stream network
Hi Jarek Jarek Jasiewicz wrote: Hi maybe this will help: http://gis.vsb.cz/GIS_Ostrava/GIS_Ova_2009/sbornik/Lists/Papers/060.pdf What were the limitations of the GRASS raster-to-vector conversion tools you encountered? It's not clear to me from Fig. 3, but I would like to understand the problem, and I hope you can explain it in a bit more detail. Thanks a lot, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Building topologically correct stream network
Markus Metz pisze: Hi Jarek Jarek Jasiewicz wrote: Hi maybe this will help: http://gis.vsb.cz/GIS_Ostrava/GIS_Ova_2009/sbornik/Lists/Papers/060.pdf What were the limitations of the GRASS raster-to-vector conversion tools you encountered? that problems are shown on fig 4) 1) streams directions are not correct (are random) because we have no information about its directions in raster file 2) there are problems if streamlines toutch each other like It's not clear to me from Fig. 3, but I would like to understand the problem, and I hope you can explain it in a bit more detail. In fact I use script to convert Terrain Analysis System data to Grass topological data. I attached it here (note it is not GRASS script with full error handling). Sorry for polish comment but I have not time to translate it. If you need more information how it works contact me off list. Jarek BTW I'm about to finish regular grass modue for calculating Strahler Stream order and Shreeve Stream Magnitude. I hope I publish it for testing in next few days. Any help appreciated. Thanks a lot, Markus M #!/bin/bash zlewnia=`eval g.gisenv get=MAPSET` cieki='cieki_'$zlewnia vertex='vertex_'$zlewnia vector='index_vector' links='index' #unikalne numery strumieni class='class' #klasy punktów/odcinków (arc/nodes) strahler='strahler' #plik z strahler stream order shrevee='sreeve' #plik z strahler stream order main='main' #hierarchia zlewni zlewnia_elem='elem' z1='w1' z2='w2' z3='w3' z4='w4' z5='w5' z6='w6' #parametry morfometryczne strumieni elev='elev' elev2main='elev2main' slope='slope' planc='planc' profc='profc' maxic='maxic' diff5='diff5' diff9='diff9' diff15='diff15' ### Koniec definicji ### db.connect driver=pg database='rzeki_mgr' g.mremove -f rast=tmp* vect=tmp* 2/dev/null #usuÅ tempy jeżeli istniejÄ g.remove -f vect=$cieki 2/dev/null #usuÅ wynikowy plik jeÅżeli istnieje echo 'DROP table tmp_streams'|db.execute 2/dev/null #usuÅ tabelÄ tmp_streams jeżeli istnieje echo 'DROP table tmp_streams_param'|db.execute 2/dev/null #usuÅ tabelÄ tmp_streams jeżeli istnieje #echo 'DROP table tmp_katy'|db.execute 2/dev/null #usuÅ tabelÄ tmp_streams jeżeli istnieje echo 'DROP table tmp_join'|db.execute 2/dev/null #usuÅ tabelÄ tmp_streams jeżeli istnieje okreÅlenie regionu analizy g.region rast=$class #ustaw region wzglÄdem pliku class ### procedura wygenerowania punktów z mapy # najpierw sprawdzamy, czy plik ma kategorie jeżeli nie koÅczymy tp=`eval v.category -g input=$vector type=line option=report`/dev/null if [ -z $tp ] then echo wektor nie ma kategorii. PrzewidujÄ problemy, wiÄc koÅczÄ exit 1 fi v.to.points -n input=$vector output=tmp_point type=point,line llayer=1 dmax=100 v.db.addcol map=tmp_point layer=2 'columns=rzad INTEGER, magn INTEGER, link INTEGER, is_main INTEGER, zlewnia_elem INTEGER, z1 INTEGER, z2 INTEGER, z3 INTEGER, z4 INTEGER, z5 INTEGER, z6 INTEGER' #probkowanie rastrow v.what.rast vector=tmp_point raster=$links layer=2 column=link v.what.rast vector=tmp_point raster=$strahler layer=2 column=rzad v.what.rast vector=tmp_point raster=$shrevee layer=2 column=magn v.what.rast vector=tmp_point raster=$main layer=2 column=is_main echo 'zlewnia_elem' v.what.rast vector=tmp_point raster=$zlewnia_elem layer=2 column=zlewnia_elem v.what.rast vector=tmp_point raster=$z1 layer=2 column=z1 v.what.rast vector=tmp_point raster=$z2 layer=2 column=z2 v.what.rast vector=tmp_point raster=$z3 layer=2 column=z3 v.what.rast vector=tmp_point raster=$z4 layer=2 column=z4 v.what.rast vector=tmp_point raster=$z5 layer=2 column=z5 v.what.rast vector=tmp_point raster=$z6 layer=2 column=z6 #tworzenie tabeli zawierajÄ cej uÅrednione parametry morfometryczne strumieni v.to.points -i input=$vector output=$vertex type=line dmax=20 v.db.addcol map=$vertex layer=2 'columns=elev REAL, elev2main REAL, elev2stream REAL, slope REAL, planc REAL, profc REAL, maxic REAL, diff5 REAL, diff9 REAL, diff15 REAL' #probkowanie dużej tabeli. Tabela vertex zostaje zachowana ze wzglÄdu na przydatnoÅÄ do innych celów v.what.rast vector=$vertex raster=$elev layer=2 column=elev v.what.rast vector=$vertex raster=$elev2main layer=2 column=elev2main v.what.rast vector=$vertex raster=$slope layer=2 column=slope v.what.rast vector=$vertex raster=$planc layer=2 column=planc v.what.rast vector=$vertex raster=$profc layer=2 column=profc v.what.rast vector=$vertex raster=$maxic layer=2 column=maxic v.what.rast vector=$vertex raster=$diff5 layer=2 column=diff5 v.what.rast vector=$vertex raster=$diff9 layer=2 column=diff9 v.what.rast vector=$vertex raster=$diff15 layer=2 column=diff15 echo 'CREATE TABLE tmp_streams_param AS SELECT lcat AS lcat, avg(elev) AS rzedna, avg(elev2main) AS elev2mean, avg(slope) AS slope, avg(planc) AS planc, avg(profc) AS profc, avg(maxic) AS maxic, avg(diff5) AS diff5, avg(diff9) AS diff9, avg(diff15) AS diff15 FROM '$vertex'_2 GROUP BY
Re: [GRASS-user] Re: Building topologically correct stream network
Jarek Jasiewicz pisze: Markus Metz pisze: Hi Jarek Jarek Jasiewicz wrote: Hi maybe this will help: http://gis.vsb.cz/GIS_Ostrava/GIS_Ova_2009/sbornik/Lists/Papers/060.pdf What were the limitations of the GRASS raster-to-vector conversion tools you encountered? that problems are shown on fig 4) 1) streams directions are not correct (are random) because we have no information about its directions in raster file 2) there are problems if streamlines toutch each other like Opps... there is the problem in the paper: fig 3 shoud be fig 4 and fig 4 shoud be 3 sorry Jarek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Building topologically correct stream network
On Tuesday 21 July 2009, Jarek Jasiewicz wrote: BTW I'm about to finish regular grass modue for calculating Strahler Stream order and Shreeve Stream Magnitude. I hope I publish it for testing in next few days. Any help appreciated. I would be happy to help test! Cheers, Dylan -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University of California at Davis 530.754.7341 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Building topologically correct stream network
I just don't know how. Tell me what to do - I am new to grass, and I am a intermediate in compiling software. Also, I am running mac os x if that makes a difference, and I can also test on a debian ppc linux laptop. thanks Stephen Sefick On Tue, Jul 21, 2009 at 2:32 PM, Dylan Beaudettedebeaude...@ucdavis.edu wrote: On Tuesday 21 July 2009, Jarek Jasiewicz wrote: BTW I'm about to finish regular grass modue for calculating Strahler Stream order and Shreeve Stream Magnitude. I hope I publish it for testing in next few days. Any help appreciated. I would be happy to help test! Cheers, Dylan -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University of California at Davis 530.754.7341 ___ 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