<div dir="ltr"><br><br>On Tue, Nov 20, 2018 at 9:50 AM Ken Mankoff <<a href="mailto:mankoff@gmail.com">mankoff@gmail.com</a>> wrote:<br>><br>> Hi GRASS list,<br>><br>> I'm seeking some advice about how to determine the eventual edge-cell destination of every pixel in a DEM. I have an algorithm that works on rasters up to 0.5 billion pixels but am now trying to do it on a region about 8x as large and it is failing.<br>><br>> I've spent a week exploring r.watershed and r.terraflow and can't figure out how to make it work.<br>><br><div>> What has worked for the 0.5 B pixel region is running r.watershed and then for every negative pixel, I run r.water.outlet and find the upstream contributors. But I can't run r.watershed on the entire large region. I can run it quickly on sub-regions, but some very large interior basins cause this method to fail.</div><div><br></div><div>How does it fail? What is the error output? You did use the -m flag for r.watershed, right?<br></div>><br>> I've run r.terraflow, but the direction raster from that has a massive number of internal gaps and I don't see how to get it to work.<br><div><br></div><div>You don't need r.watershed or r.terraflow for this purpose., You can run r.stream.extract with a DEM as the only input, then you get the outlets and can feed them to r.stream.basins.<br></div><div><br></div><div></div><div>></div>> Does anyone have any advice how to process large rasters with r.watershed, or r.fill.dir, or some other method that can trace all pixels to their edge destination?<br><div><br></div><div>r.watershed should process very large rasters, granted that there is enough free disk space.</div><div><br></div><div>Markus M</div><div><br></div><div>></div>> Thanks,<br>><br>>   -k.<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></div>