[GRASS-SVN] r33997 - grass-addons/raster/r.viewshed
svn_grass at osgeo.org
svn_grass at osgeo.org
Fri Oct 24 04:11:46 EDT 2008
Author: neteler
Date: 2008-10-24 04:11:46 -0400 (Fri, 24 Oct 2008)
New Revision: 33997
Modified:
grass-addons/raster/r.viewshed/testscript.sh
Log:
line breaks fixed
Modified: grass-addons/raster/r.viewshed/testscript.sh
===================================================================
--- grass-addons/raster/r.viewshed/testscript.sh 2008-10-24 01:22:20 UTC (rev 33996)
+++ grass-addons/raster/r.viewshed/testscript.sh 2008-10-24 08:11:46 UTC (rev 33997)
@@ -5,14 +5,12 @@
# create first hemisphere
g.region n=1000 s=0 w=0 e=1000 -p res=1
-r.mapcalc 'disk.15031=if(sqrt((col() - 500)^2 + (500 -
-row())^2)<500,sqrt((col() - 500)^2 + (500 - row())^2),null())'
+r.mapcalc 'disk.15031=if(sqrt((col() - 500)^2 + (500 - row())^2)<500,sqrt((col() - 500)^2 + (500 - row())^2),null())'
r.mapcalc 'hemisphere1=500 * sin(acos (disk.15031/500))'
# create second hemisphere
g.region n=500 s=0 w=0 e=500 -p res=1
-r.mapcalc 'disk.14947=if(sqrt((col() - 500)^2 + (500 -
-row())^2)<500,sqrt((col() - 500)^2 + (500 - row())^2),null())'
+r.mapcalc 'disk.14947=if(sqrt((col() - 500)^2 + (500 - row())^2)<500,sqrt((col() - 500)^2 + (500 - row())^2),null())'
r.mapcalc 'hemisphere2=500 * sin(acos (disk.14947/500))'
g.remove --q rast=disk.14947,disk.15031
# merge both
@@ -22,8 +20,7 @@
d.rast hemisphere
# run r.viewshed
-r.viewshed hemisphere out=hemisphere_viewshed coord=250,250
-max=1000000 obs=100 mem=2000 --o
+r.viewshed hemisphere out=hemisphere_viewshed coord=250,250 max=1000000 obs=100 mem=2000 --o
r.shaded.relief hemisphere --o
d.his h=hemisphere_viewshed i=hemisphere.shade
@@ -33,3 +30,4 @@
d.his h=hemisphere_los i=hemisphere.shade
nviz hemisphere col=hemisphere_viewshed
+
More information about the grass-commit
mailing list