[GRASS5] [bug #4022] (grass) problem with r.shaded.relief
Request Tracker
grass-bugs at intevation.de
Fri Jan 20 03:45:18 EST 2006
this bug's URL: http://intevation.de/rt/webrt?serial_num=4022
-------------------------------------------------------------------------
Subject: problem with r.shaded.relief
Platform: Mac OSX
grass obtained from: Other (CDROM etc)
grass binary for platform: Downloaded precompiled Binaries
GRASS Version: 060407
problem with the last version of Grass6.1cvs of 7 january (http://wwwamb.bologna.enea.it/forgrass/index.htm), particularly with r.shade.relief :
"r.shaded.relief map=dempierre shadedmap=dempierreshaded50 altitude=50 azimuth=270 zmult=4 scale=1
Calculating shading, please stand by.
/Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/scripts/r.shaded.relief: line 186: 4287 Bus error r.mapcalc <<EOF
$ELEVSHADE = eval( \\
x=($zmult*$elev[-1,-1] + 2*$zmult*$elev[0,-1] + $zmult*$elev[1,-1] \\
-$zmult*$elev[-1,1] - 2*$zmult*$elev[0,1] - $zmult*$elev[1,1])/(8.*ewres()*$scale) , \\
y=($zmult*$elev[-1,-1] + 2*$zmult*$elev[-1,0] + $zmult*$elev[-1,1] \\
-$zmult*$elev[1,-1] - 2*$zmult*$elev[1,0] - $zmult*$elev[1,1])/(8.*nsres()*$scale) , \\
slope=90.-atan(sqrt(x*x + y*y)), \\
a=round(atan(x,y)), \\
a=if(isnull(a),1,a), \\
aspect=if(x!=0||y!=0,if(a,a,360.)), \\
cang = sin($alt)*sin(slope) + cos($alt)*cos(slope) * cos($az-aspect), \\
if(cang < 0.,0.,100.*cang), \\
if(isnull(cang), null(), 100.*cang))
EOF
ERROR: dempierreshaded50 - map not found
Shaded relief map created and named [dempierreshaded50]"
.
-------------------------------------------- Managed by Request Tracker
More information about the grass-dev
mailing list