On Mon, Apr 8, 2013 at 11:44 AM,  <[email protected]> wrote:
> I discovered the following in r.watershed/ram/slope_len.c where the slope
> length is calculated given two cell positions in the raster.
>
> if (r == dr)
>     res = window.ns_res;
> else if (c == dc)
>     res = window.ew_res;
> else
>     res = diag;
>
> Here we are dealing with a cell given by coordinates (r,c) which flows into
> another cell (dr,dc). So the code reads to me that if the two cells are on
> the same row, add to the slope the map's North-South resolution. Should it
> not be West-East since the slope is horizontal rather than vertical and vice
> versa for the columns?

yes

>
> This seems like a bug to me and the code is taken from the latest snapshot.

This bug has been present ever since (GRASS 4.0). Fixed in trunk r55666.

Markus M
_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev

Reply via email to