[GRASS-SVN] r60531 - grass-addons/grass7/raster/r.basin

svn_grass at osgeo.org svn_grass at osgeo.org
Wed May 28 01:36:11 PDT 2014


Author: hellik
Date: 2014-05-28 01:36:11 -0700 (Wed, 28 May 2014)
New Revision: 60531

Modified:
   grass-addons/grass7/raster/r.basin/r.basin.py
Log:
r.basin: flat pixels remains flat pixels

Modified: grass-addons/grass7/raster/r.basin/r.basin.py
===================================================================
--- grass-addons/grass7/raster/r.basin/r.basin.py	2014-05-28 08:30:38 UTC (rev 60530)
+++ grass-addons/grass7/raster/r.basin/r.basin.py	2014-05-28 08:36:11 UTC (rev 60531)
@@ -378,7 +378,7 @@
         # counterclockwise: 90deg is North, 180 is West, 270 is South 360 is East. 
         # The aspect value 0 is used to indicate undefined aspect in flat areas with slope=0.
         # We calculate the number of degree from north, increasing counterclockwise.
-        grass.mapcalc("$r_aspect_mod = if($r_aspect > 90, 450 - $r_aspect, 90 - $r_aspect)",
+        grass.mapcalc("$r_aspect_mod = if($r_aspect == 0, 0, if($r_aspect > 90, 450 - $r_aspect, 90 - $r_aspect))",
                   r_aspect = r_aspect,
                   r_aspect_mod = r_aspect_mod)
         grass.message("r.mapcalc done")



More information about the grass-commit mailing list