<div dir="ltr"><div><div><br><br>On Wed, May 2, 2018 at 8:29 PM, Anna Petrášová <<a href="mailto:kratochanna@gmail.com">kratochanna@gmail.com</a>> wrote:<br>><br>> On Wed, May 2, 2018 at 2:21 PM, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br>> > On Wed, May 2, 2018 at 7:13 PM, Martin Landa <<a href="mailto:landa.martin@gmail.com">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>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<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">http://desktop.arcgis.com/en/arcmap/10.5/tools/spatial-analyst-toolbox/flow-accumulation.htm</a><br>> >><br>> >> --<br>> >> Martin Landa<br>> >> <a href="http://geo.fsv.cvut.cz/gwiki/Landa">http://geo.fsv.cvut.cz/gwiki/Landa</a><br>> >> <a href="http://gismentors.cz/mentors/landa">http://gismentors.cz/mentors/landa</a><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>> > --<br>> > Markus Neteler, PhD<br>> > <a href="http://www.mundialis.de">http://www.mundialis.de</a> - free data with free software<br>> > <a href="http://grass.osgeo.org">http://grass.osgeo.org</a><br>> > <a href="http://courses.neteler.org/blog">http://courses.neteler.org/blog</a><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>> 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></div></div></div>