[GRASS-SVN] r66607 - grass/branches/releasebranch_7_0/raster/r.regression.line
svn_grass at osgeo.org
svn_grass at osgeo.org
Mon Oct 26 06:42:58 PDT 2015
Author: neteler
Date: 2015-10-26 06:42:58 -0700 (Mon, 26 Oct 2015)
New Revision: 66607
Modified:
grass/branches/releasebranch_7_0/raster/r.regression.line/r.regression.line.html
Log:
r.regression.line manual: update example to NC dataset; HTML cosmetics
Modified: grass/branches/releasebranch_7_0/raster/r.regression.line/r.regression.line.html
===================================================================
--- grass/branches/releasebranch_7_0/raster/r.regression.line/r.regression.line.html 2015-10-26 13:42:45 UTC (rev 66606)
+++ grass/branches/releasebranch_7_0/raster/r.regression.line/r.regression.line.html 2015-10-26 13:42:58 UTC (rev 66607)
@@ -1,8 +1,18 @@
<h2>DESCRIPTION</h2>
-<em>r.regression.line</em> Calculates linear regression from two raster maps,
-according to the formula y = a + b*x, where x and y represent raster maps.
-Optionally saves regression coefficients to an ASCII file.
+<em>r.regression.line</em> calculates a linear regression from two
+raster maps, according to the formula
+<div class="code"><pre>
+y = a + b*x
+</pre></div>
+where
+<div class="code"><pre>
+x
+y
+</pre></div>
+represent the input raster maps.
+<p>
+Optionally, it saves regression coefficients as a ASCII file.
The result includes the following coefficients:
offset/intercept (a) and gain/slope (b), correlation coefficient (R),
number of elements (N), means (medX, medY), standard deviations
@@ -15,29 +25,43 @@
<h2>EXAMPLE</h2>
-Comparison of the old and the new DEM in Spearfish:
+Comparison of two DEMs (SRTM and NED, both at 30m resolution),
+provided in the North Carolina sample dataset:
+
<div class="code"><pre>
-g.region raster=elevation.10m -p
-r.regression.line map1=elevation.dem map2=elevation.10m
+g.region raster=elev_srtm_30m -p
+r.regression.line mapx=elev_ned_30m mapy=elev_srtm_30m
+ y = a + b*x
+ a (Offset): -1.659279
+ b (Gain): 1.043968
+ R (sumXY - sumX*sumY/N): 0.894038
+ N (Number of elements): 225000
+ F (F-test significance): 896093.366283
+ meanX (Mean of map1): 110.307571
+ sdX (Standard deviation of map1): 20.311998
+ meanY (Mean of map2): 113.498292
+ sdY (Standard deviation of map2): 23.718307
</pre></div>
<p>
Using the script style flag AND <em>eval</em> to make results
available in the shell:
<div class="code"><pre>
-g.region raster=elevation.10m -p
-eval `r.regression.line -g map1=elevation.dem map2=elevation.10m`
+g.region raster=elev_srtm_30m -p
+eval `r.regression.line -g mapx=elev_ned_30m mapy=elev_srtm_30m`
+
+# print result stored in respective variables
echo $a
-479.615
+-1.659279
echo $b
-0.645631
+1.043968
echo $R
-0.804441
+0.894038
</pre></div>
-<h2>AUTHOR</h2>
+<h2>AUTHORS</h2>
Dr. Agustin Lobo - alobo at ija.csic.es<br>
Updated to GRASS 5.7 Michael Barton, Arizona State University<br>
More information about the grass-commit
mailing list