<div dir="ltr">

<span style="font-size:12.8px;text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Thanks all for your suggestions!</span><div style="font-size:12.8px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div style="font-size:12.8px;text-decoration-style:initial;text-decoration-color:initial">I created a comparison of methods [1] which were mentioned. But unfortunately, the results are still very different (at least from my point of view). What do you think about it? What approach should be chosen?</div><div style="font-size:12.8px;text-decoration-style:initial;text-decoration-color:initial"><br></div><div style="font-size:12.8px;text-decoration-style:initial;text-decoration-color:initial">[1] <a href="https://docs.google.com/document/d/15MyV-l2BH5agjN1T3LtQBotH7RJsIL_5K7exSZypTnQ/edit?usp=sharing" target="_blank" style="color:rgb(17,85,204)">https://docs.google.com/<wbr>document/d/15MyV-<wbr>l2BH5agjN1T3LtQBotH7RJsIL_<wbr>5K7exSZypTnQ/edit?usp=sharing</a></div>

<br></div><div class="gmail_extra"><br><div class="gmail_quote">2018-05-02 22:17 GMT+02:00 Markus Metz <span dir="ltr"><<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><span class=""><div><br><br>On Wed, May 2, 2018 at 8:29 PM, Anna Petrášová <<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>> wrote:<br>><br>> On Wed, May 2, 2018 at 2:21 PM, Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>> wrote:<br>> > On Wed, May 2, 2018 at 7:13 PM, Martin Landa <<a href="mailto:landa.martin@gmail.com" target="_blank">landa.martin@gmail.com</a>> wrote:<br>> >> Hi,<br>> >><br>> >> one of my colleague is dealing with flow accumulation calculation in<br>> >> GRASS, and comparing results with Esri ArcGIS (builtin Flow<br>> >> Accumulation [1]). Tool in ArcGIS to me knowledge uses D8, so we run<br>> >> r.watershed with -s flag and also -a flag to avoid negative values.<br>> >> Resultant values are very different.<br>> >><br>> >> * Esri ArcGIS (range: 0 - 20741)<br>> >> * r.watershed -s (range: -34068 - 5697)<br>> >> * r.watershed -sa (1 - 34068)<br>> ><br>> > I'd suggest to look at the univariate statistics first (r.univar -e)<br>> > and compare that instead of min/max only.<br>><br>> Also, r.watershed uses least-cost path, unlike the arcmap tool. Did<br>> you fill the depressions before running the arcmap tool?<br><br></div></span>If you want to get nearly identical results, use r.hydrodem to remove sinks from the DEM, then use this conditioned DEM as input for both ArcGIS and r.watershed. Note that with D8 there is some randomness about where surface flow will go to in case of several neighbors with the same elevation. In any case, if surface accumulation of ArcGIS and r.watershed differ, literature cited in the corresponding manuals should allow to decide which method is more appropriate for the current task.<br><br></div>Markus M<div><div class="h5"><br><div><div><br>><br>> Anna<br>><br>> ><br>> > markusN<br>> ><br>> >> Do you have idea why results are so different? Thanks for any pointers<br>> >> in advance! Martin<br>> >><br>> >> [1] <a href="http://desktop.arcgis.com/en/arcmap/10.5/tools/spatial-analyst-toolbox/flow-accumulation.htm" target="_blank">http://desktop.arcgis.com/en/<wbr>arcmap/10.5/tools/spatial-<wbr>analyst-toolbox/flow-<wbr>accumulation.htm</a><br>> >><br>> >> --<br>> >> Martin Landa<br>> >> <a href="http://geo.fsv.cvut.cz/gwiki/Landa" target="_blank">http://geo.fsv.cvut.cz/gwiki/<wbr>Landa</a><br>> >> <a href="http://gismentors.cz/mentors/landa" target="_blank">http://gismentors.cz/mentors/<wbr>landa</a><br>> >> ______________________________<wbr>_________________<br>> >> grass-user mailing list<br>> >> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>> >> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br>> ><br>> ><br>> ><br>> > --<br>> > Markus Neteler, PhD<br>> > <a href="http://www.mundialis.de" target="_blank">http://www.mundialis.de</a> - free data with free software<br>> > <a href="http://grass.osgeo.org" target="_blank">http://grass.osgeo.org</a><br>> > <a href="http://courses.neteler.org/blog" target="_blank">http://courses.neteler.org/<wbr>blog</a><br>> > ______________________________<wbr>_________________<br>> > grass-user mailing list<br>> > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br>> ______________________________<wbr>_________________<br>> grass-user mailing list<br>> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br><br></div></div></div></div></div>
</blockquote></div><br></div>