epifanio
/
grassgis
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

testsuite rhemisphere.sh: update r.mapcalc call 

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@65691 15284696-431f-4ddb-bdfa-cd5b030d7da7
py3
neteler 2015-07-20 10:37:29 +00:00
parent ea1e61112f
commit 5d8d7ebaf4
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
testsuite/raster/rhemisphere.sh
View File

@ -36,13 +36,13 @@ Y="($RADIUS - row())"
r="sqrt($X^2 + $Y^2)"
#Mask out unwanted parts (check for <= ??):
r.mapcalc "$TMP=if($r<$RADIUS,$r,null())"
r.mapcalc "$TMP = if($r<$RADIUS,$r,null())"
ALPHA="acos ($TMP/$RADIUS)"
HEIGHT="$RADIUS * sin($ALPHA)"
r.mapcalc "hemisphere=$HEIGHT"
r.mapcalc "hemisphere = $HEIGHT"
cleanup
g.message "Generated raster map <hemisphere>"
#echo "Now generate aspect + slope on <hemisphere>"