[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