[GRASS-SVN] r58553 - grass/trunk/raster/r.horizon
svn_grass at osgeo.org
svn_grass at osgeo.org
Sun Dec 29 16:10:47 PST 2013
Author: neteler
Date: 2013-12-29 16:10:47 -0800 (Sun, 29 Dec 2013)
New Revision: 58553
Added:
grass/trunk/raster/r.horizon/rhorizon_shaded_dem_point.png
grass/trunk/raster/r.horizon/rhorizon_singlepoint_plot.png
grass/trunk/raster/r.horizon/rhorizon_testplot.sh
Modified:
grass/trunk/raster/r.horizon/r.horizon.html
Log:
r.horizon manual: screenshots added; rhorizon_testplot.sh added
Modified: grass/trunk/raster/r.horizon/r.horizon.html
===================================================================
--- grass/trunk/raster/r.horizon/r.horizon.html 2013-12-29 23:43:08 UTC (rev 58552)
+++ grass/trunk/raster/r.horizon/r.horizon.html 2013-12-30 00:10:47 UTC (rev 58553)
@@ -43,7 +43,7 @@
<p>The <i>horizon_step</i> parameter gives the angle step (in degrees)
between successive azimuthal directions for the calculation of the
horizon. Thus, a value of 5 for the <i>horizon_step</i> will give a total of
-360/5=72 directions (72 raster maps if used in the raster mode).
+360/5=72 directions (72 raster maps if used in the raster map mode).
<p>The <i>direction</i> parameter gives the initial direction of the
first output. This parameter acts as an direction angle offset. For
@@ -95,7 +95,7 @@
<p>The <i>output</i> parameter allows to save the resulting horizon
angles in a comma separated ASCII file (single point mode only). If
-you use <b>r.horizon</b> in the raster mode this option will be ignored.
+you use <b>r.horizon</b> in the raster map mode this option will be ignored.
<p>At the moment the elevation and maximum distance must be measured in meters,
even if you use geographical coordinates (longitude/latitude). If your
@@ -124,8 +124,10 @@
<h2>EXAMPLES</h2>
-The examples are intended for the North Carolina sample dataset:
-<p>
+The examples are intended for the North Carolina sample dataset.
+
+<h3>Single point mode</h3>
+
Single point mode (output of horizon angles CCW from East):
<div class="code"><pre>
# determine horizon angle in 225 degree direction:
@@ -136,9 +138,22 @@
# save result in CSV file
r.horizon elev_in=elevation direction=90 horizon_step=5 bufferzone=200 \
coordinate=638871.6,223384.4 maxdistance=5000 output=horizon.csv
+
+# test point near high way intersection
+g.region n=223540 s=220820 w=634650 e=638780 res=10 -p
+r.horizon elev_in=elevation direction=0 horizon_step=5 bufferzone=200 \
+ coordinate=636483.54,222176.25 maxdistance=5000 -d output=horizon.csv
</pre></div>
+<img src="rhorizon_shaded_dem_point.png"><br>
+Test point near high way intersection (North Carolina sample dataset)
+<p>
+<img src="rhorizon_singlepoint_plot.png"><br>
+Horizon angles for test point (CCW from East)
+<p>
+<h3>Raster map mode</h3>
+
Raster map mode (output maps "horangle*" become input for <em>r.sun</em>):
<div class="code"><pre>
# we put a bufferzone of 10% of maxdistance around the study area
Added: grass/trunk/raster/r.horizon/rhorizon_shaded_dem_point.png
===================================================================
(Binary files differ)
Property changes on: grass/trunk/raster/r.horizon/rhorizon_shaded_dem_point.png
___________________________________________________________________
Added: svn:mime-type
+ image/png
Added: grass/trunk/raster/r.horizon/rhorizon_singlepoint_plot.png
===================================================================
(Binary files differ)
Property changes on: grass/trunk/raster/r.horizon/rhorizon_singlepoint_plot.png
___________________________________________________________________
Added: svn:mime-type
+ image/png
Added: grass/trunk/raster/r.horizon/rhorizon_testplot.sh
===================================================================
--- grass/trunk/raster/r.horizon/rhorizon_testplot.sh (rev 0)
+++ grass/trunk/raster/r.horizon/rhorizon_testplot.sh 2013-12-30 00:10:47 UTC (rev 58553)
@@ -0,0 +1,22 @@
+# Horizon angles CCW from East
+
+# Horizon angle output in degree
+
+# NC dataset
+DEM=elevation
+# test point near high way intersection
+coords=636483.54,222176.25
+
+g.region n=223540 s=220820 w=634650 e=638780 res=10 -p
+
+# direction=0 is East, direction=90 is North
+r.horizon elev_in=$DEM direction=0 horizon_step=5 bufferzone=200 \
+ coordinate=$coords maxdistance=5000 -d output=horizon.csv --o
+
+echo 'set datafile separator ","
+
+set ylabel "Horizon angle"
+set xlabel "Angle (CCW from East)"
+
+plot "horizon.csv" using 1:2 with lines' > gnuplot.txt
+gnuplot -persist gnuplot.txt
Property changes on: grass/trunk/raster/r.horizon/rhorizon_testplot.sh
___________________________________________________________________
Added: svn:executable
+ *
Added: svn:mime-type
+ text/x-sh
Added: svn:eol-style
+ native
More information about the grass-commit
mailing list