<div dir="ltr"><div><div><br><br>On Fri, Nov 17, 2017 at 7:52 AM, Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> wrote:<br>><br>> Hello<br>><br>> On 11/16/2017 11:51 PM, Giuseppe Amatulli wrote:<br>><br>><br>> Hi Markus  M.<br>><br>> I was testing the r.stream.extract for 2 tiff <br>><br>> 1) 80040 x 72870 = 5,832,514,800 i got the stream results - no error <br>><br>> 2) 80490 x 142690 = 11,485,118,100 i got the following error <br>><br>> A* Search...<br>> 0..2../var/spool/slurmd/job6514787/slurm_script: line 80: 28925 Bus error               r.stream.extract elevation=elv accumulation=upa threshold=0.5 depression=dep direction=dir stream_raster=stream memory=45000 --o --verbose<br>><br>><br>> I think it was already mentioned in this thread: the threshold=0.5 is certainly wrong. The threshold value is *number of pixels* for minimum basin size. Usually it would be in the thousands. With your region size of 11 billion you probably want a threshold of tens of thousands.<br><br></div>As Giuseppe explained previously, accumulation has been rescaled to square kilometers, therefore the threshold is ok (still a bit small, but not nonsense).<br><br></div>Markus M<br><div><div>><br>> Do you think that is something with the slurm ram limitation or is something with the r.extract.stream? <br>> If the the stream output have a number of stream segments larger than 2,147,483,647 what is happen? <br>> Do I get an error or all the value larger than   2,147,483,647 are just rounded to 2,147,483,647 ? <br>><br>><br>> Moreover, if use the stream obtain from option 1) and I use the stream as input for the r.stream.basins  <br>> I got the following error<br>><br>> reading raster map <stream>...<br>> 0..3..6..9..12..15..18..21..24..27..30..33..36..39..42..45../var/spool/slurmd/job6514788/slurm_script: line 82: 17687 Bus error               /gpfs/home/fas/sbsc/ga254/.grass7/addons/bin/r.stream.basins -l stream_rast=stream direction=dir<br>><br>> is this something that need to be fixed in r.stream.basins, or should i think that is due to other problems <br>><br>> Thank you <br>> Giuseppe <br>><br>><br>><br>><br>> On 10 November 2017 at 11:21, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br>>><br>>> On Wed, Nov 1, 2017 at 10:12 PM, Markus Metz<br>>> <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> wrote:<br>>> > On Wed, Nov 1, 2017 at 7:15 PM, Giuseppe Amatulli<br>>> > <<a href="mailto:giuseppe.amatulli@gmail.com">giuseppe.amatulli@gmail.com</a>> wrote:<br>>> >><br>>> >> Thanks Markus!!<br>>> >> I will test and I will let you know how it works.<br>>> ><br>>> > Your feedback is very helpful!<br>>> >><br>>> >> I have few  more questions<br>>> >> 1) now how much is the upper limit matrix cell number that<br>>> >> r.stream.extract can handle?<br>>> ><br>>> > About 1.15e+18 cells.<br>>> ><br>>> > Another limitation is the number of detected stream segments. This must not<br>>> > be larger than 2,147,483,647 streams,<br>>> ...<br>>><br>>> (Added as a note to<br>>>  <a href="https://grasswiki.osgeo.org/wiki/GRASS_GIS_Performance#Some_benchmarks">https://grasswiki.osgeo.org/wiki/GRASS_GIS_Performance#Some_benchmarks</a><br>>> )<br>>><br>>> best,<br>>> markusN<br>><br>><br>><br>><br>> --<br>> Giuseppe Amatulli, Ph.D.<br>><br>> Research scientist at<br>> Yale School of Forestry & Environmental Studies<br>> Yale Center for Research Computing<br>> Center for Science and Social Science Information<br>> New Haven, 06511<br>> Teaching: <a href="http://spatial-ecology.net">http://spatial-ecology.net</a><br>> Work:  <a href="https://environment.yale.edu/profile/giuseppe-amatulli/">https://environment.yale.edu/profile/giuseppe-amatulli/</a><br>><br>><br>> _______________________________________________<br>> grass-user mailing list<br>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>><br>><br>> -- <br>> Micha Silver<br>> Ben Gurion Univ.<br>> Sde Boker, Remote Sensing Lab<br>> cell: +972-523-665918<br><br></div></div></div>