[GRASS-SVN] r66237 - grass-addons/grass7/raster/r.biodiversity

[svn_grass at osgeo.org]

Author: pvanbosgeo
Date: 2015-09-15 13:14:34 -0700 (Tue, 15 Sep 2015)
New Revision: 66237

Modified:
   grass-addons/grass7/raster/r.biodiversity/r.biodiversity.py
Log:
Correction in Renyi formula, to better deal with 0 values

Modified: grass-addons/grass7/raster/r.biodiversity/r.biodiversity.py
===================================================================
--- grass-addons/grass7/raster/r.biodiversity/r.biodiversity.py	2015-09-15 20:10:09 UTC (rev 66236)
+++ grass-addons/grass7/raster/r.biodiversity/r.biodiversity.py	2015-09-15 20:14:34 UTC (rev 66237)
@@ -205,7 +205,7 @@
         tmpl = []
         
         if Q[n] == 1:      
-            # When alpha = 1
+            # If alpha = 1
             for i in xrange(len(IN)):
                 tmpi = tmpname('shi' + str(i) + "_")
                 tmpl.append(tmpi)
@@ -226,7 +226,7 @@
             for i in xrange(len(IN)):
                 tmpi = tmpname('reni')
                 tmpl.append(tmpi)       
-                grass.mapcalc("$tmpi = pow($inl/$tmpt,$alpha)", 
+                grass.mapcalc("$tmpi = if($inl==0,0,pow($inl/$tmpt,$alpha))", 
                               tmpi=tmpi,
                               tmpt=tmpt,
                               inl=IN[i],