[GRASS-user] Problem with slope modeling .

Radomir dod.mort at gmail.com
Wed May 2 08:14:35 EDT 2007


Dnia 01-05-2007, wto o godzinie 21:23 +1200, Hamish napisał(a):
> Radomir wrote:
> > The r.surf.contour methods gives much worse results then RST (see the
> > Elevation graph).
> 
> r.surf.contour will only work with interger values (so far).
> can you try something like:
> 
> r.mapcalc "map10000= int(map * 10000)"
> r.surf.contour in=map10000 out=map10000_rsc
> r.mapcalc "map_rsc = map10000_rsc / 10000.0"
> 
> then you can produce a smoother result in spite of r.surf.contour's
> limitations.
> 
I did as you suggested:
r.mapcalc "izolinie_rsc_izolinie_10k = int(izolinie_rsc_izolinie
*10000)"
r.surf.contour input=izolinie_rsc_izolinie_10k output=wysokosci_rsc_10k
r.mapcalc "wysokosci_rsc_male = wysokosci_rsc_10k / 10000.0"
then:
r.slope.aspect
and again I have those "waves" between each isoline....
See picture:
http://mort.no-ip.org/grass/indexgrass2.html

hmmm. Maybe there is something wrong with my whole project??
Basics about the project:
projection: 0 (x,y)
zone:       0
north:      444979.06151534
south:      438995.17651068
west:       597964.6049651
east:       604980.31251938
nsres:      1.26992466
ewres:      1.27004119
rows:       4712
cols:       5524
cells:      26029088
Elevation isolines are every 2.5m. Original map scale 1:10000.  

> 
> > It was really interesting to analyse the output vector map. It seems
> > that when the distance between two isolines is "large" GRASS "puts the
> > next isoline" very close to the previous and this is how the "waves"
> > comes from.
> 
> Could you specify the module name, not just "GRASS 'puts the next
> isoline'"? GRASS is the whole, there are many different methods within
> its modules, not just one, so it's hard to know which one you mean to
> refer to.
> 
> 
> Hamish

-- 
Radomir 

----
The Future Is Now ;)




More information about the grass-user mailing list