[GRASS-user] r.shaded.relief strangeness

2008-03-17 Thread Dylan Beaudette
Hi,


Not sure if this is very critical, but I have noticed a problem with 
r.shaded.relief. If the region is set to 10m and then a shaded rel. map is 
produced from a 30m input data, a corrupt map is produced. 


example with spearfish data:

# set region, and force 10m resolution
g.region rast=elevation.10m -p

# compute shaded rel. map from 30m res data:
# see image [1]
r.shaded.relief map=elevation.dem shade=shade

# compute shaded rel. map from 10m res data:
# see image [2]
r.shaded.relief map=elevation.10m shade=shade --o


1. http://169.237.35.250/~dylan/temp/bad_shade.png
2. http://169.237.35.250/~dylan/temp/good_shade.png


This was done with today's SVN. 


-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


RE: [GRASS-user] r.shaded.relief strangeness

2008-03-17 Thread Glynn Clements

Patton, Eric wrote:

 Not sure if this is very critical, but I have noticed a problem with 
 r.shaded.relief. If the region is set to 10m and then a shaded rel. map is 
 produced from a 30m input data, a corrupt map is produced. 
 
 example with spearfish data:
 
 # set region, and force 10m resolution
 g.region rast=elevation.10m -p
 
 I assume you meant g.region rast=elevation.10m res=10 -ap ?
  
 # compute shaded rel. map from 30m res data:
 # see image [1]
 r.shaded.relief map=elevation.dem shade=shade
 
 I can't reproduce this. I can produce a Moire-patterned shaded-relief map, 
 since I didn't resample the 30m data to 10m:

 Try r.resample first?

r.resample won't help. That uses libgis' built-in nearest-neighbour
resampling, so you will get exactly the same result as with using the
lower resolution map.

Using r.resamp.interp with method=bilinear or method=bicubic should
eliminate the problems with generating slope and aspect maps.

-- 
Glynn Clements [EMAIL PROTECTED]
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user