<div dir="ltr"><div>Hi all,</div><div><br></div><div>r.fill.dir with the bellow artificial DEM leaves part of the directions unresolved and this area is marked in areas, although the area is filled in the output. Here is the code to reproduce:<br></div><div><br>g.region n=220750 s=220000 w=638300 e=639000 nsres=1 ewres=1 rows=750 cols=700 -p<br>r.mapcalc "hills = 500 * sin(row()) + 500 * sin(col())"<br>r.fill.dir input=hills output=filled direction=directions areas=areas<br><br></div>In r71809, I have committed a fix which transforms these unresolved directions to NULLs because in the previous implementation (current 7.4svn and older) internal representation of angles was leaking outside of the module.<br><div><br></div><div>I have 2 questions:</div><div><br></div><div>Should I backport r71809 to 7.4?<br></div><div><br></div><div>Should we pursue this issue further to resolve why these edge cases are not resolved, although they are filled? Does anybody has suggestions at this point? (Paper with the method is linked in the doc.)<br></div><div><br></div><div></div><div>Thank you,</div><div>Vaclav<br></div><div><br><a href="https://trac.osgeo.org/grass/changeset/71809">https://trac.osgeo.org/grass/changeset/71809</a></div><div><a href="https://grass.osgeo.org/grass74/manuals/r.fill.dir.html">https://grass.osgeo.org/grass74/manuals/r.fill.dir.html</a><br></div></div>