[GRASS-SVN] r72554 - grass/branches/releasebranch_7_2/imagery/i.vi

[svn_grass at osgeo.org]

Author: veroandreo
Date: 2018-03-24 04:36:15 -0700 (Sat, 24 Mar 2018)
New Revision: 72554

Modified:
   grass/branches/releasebranch_7_2/imagery/i.vi/i.vi.html
   grass/branches/releasebranch_7_2/imagery/i.vi/msavi2.c
Log:
i.vi: corrected MSAVI2 equation according to https://doi.org/10.1016/0034-4257(94)90134-1, fix #3385 (trunk, r72549)

Modified: grass/branches/releasebranch_7_2/imagery/i.vi/i.vi.html
===================================================================
--- grass/branches/releasebranch_7_2/imagery/i.vi/i.vi.html	2018-03-24 11:33:59 UTC (rev 72553)
+++ grass/branches/releasebranch_7_2/imagery/i.vi/i.vi.html	2018-03-24 11:36:15 UTC (rev 72554)
@@ -179,7 +179,7 @@
 <div class="code"><pre>
 msavi2( redchan, nirchan )
 
-MSAVI2 = (1/2)*(2(NIR+1)-sqrt((2*NIR+1)^2-8(NIR-red)))
+MSAVI2 = (1/2)*(2*NIR+1-sqrt((2*NIR+1)^2-8*(NIR-red)))
 </pre></div>
 
 <p>

Modified: grass/branches/releasebranch_7_2/imagery/i.vi/msavi2.c
===================================================================
--- grass/branches/releasebranch_7_2/imagery/i.vi/msavi2.c	2018-03-24 11:33:59 UTC (rev 72553)
+++ grass/branches/releasebranch_7_2/imagery/i.vi/msavi2.c	2018-03-24 11:36:15 UTC (rev 72554)
@@ -15,7 +15,7 @@
     }
     else 
     {
-	result=0.5*(2*(nirchan+1)-sqrt((2*nirchan+1)*(2*nirchan+1)-8*(nirchan-redchan)));
+	result=0.5*(2*nirchan+1-sqrt((2*nirchan+1)*(2*nirchan+1)-8*(nirchan-redchan)));
     }
     return result;
 }