Rich,<br><br><div class="gmail_quote">On Fri, Jun 22, 2012 at 6:37 PM, Rich Shepard <span dir="ltr"><<a href="mailto:rshepard@appl-ecosys.com" target="_blank">rshepard@appl-ecosys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  The r.basin wiki page tells me, "The prevalent orientation: in GRASS the<br>
aspect categories represent the number degrees of east and they increase<br>
counterclockwise: (90deg is North, 180 is West, 270 is South 360 is East).<br>
The aspect value 0 is used to indicate undefined aspect in flat areas with<br>
slope=0. We instead calculate the orientation as the number of degree from<br>
north, increasing counterclockwise."<br>
<br>
  I have 7 sub-basins for which r.basin presents prevalent orientation<br>
ranging from 0.1 to 1.5 degrees west of North. That isn't correct if this<br>
parameter is the overall orientation of the main channel in the basin. For<br>
example, sub-basin #32 has the main channel headwaters in the northwest and<br>
the outlet in the southeast with reported prevalent orientation of<br>
0.81195844521. If a compass was set on the map and degrees increased west of<br>
North, the main channel direction from headwaters to mouth is approximately<br>
225 degrees.<br>
<br>
  Ergo, I must be mis-interpreting prevalent orientation in the r.basin<br>
output and would appreciate having someone explain it to me.<br></blockquote><div> <br>The prevalent orientation is defined as the orientation of the directing vector. The directing vector is a segment which joins the basin's barycenter and the outlet.<br>
The orientation is given by the formula:<br>prevalent_orientation = math.atan(delta_y/delta_x)<br><br>where<br><br>delta_x = abs(float(basin_east) - east)<br>delta_y = abs(float(basin_north) - north)<br><br>being basin_east and basin_north the coordinates of barycenter and east and north those of the outlet.<br>
<br>HTH,<br>madi<br><br><br></div></div>-- <br>Dr. Margherita Di Leo<br>