[GRASS-user] r.width.funct in correct version

Helmut Kudrnovsky hellik at web.de
Thu Feb 21 02:37:21 PST 2019


[keep the discussion always on the ML!]

>Gesendet: Donnerstag, 21. Februar 2019 um 11:25 Uhr
>
>Von: "Juan Lopez" <cubarro at gmail.com>
>An:
>Betreff: Re: [GRASS-user] r.width.funct in correct version
>
>Dear Helmut.
>This is the output of g.region command:
> 
>
>GRASS 7.6.0 (nc):/departamento/Test/nc_basic_spm > g.region -p
>projection: 99 (Lambert Conformal Conic)
>zone:       0
>datum:      nad83
>ellipsoid:  a=6378137 es=0.006694380022900787
>north:      228500
>south:      215000
>west:       630000
>east:       645000
>nsres:      10
>ewres:      10
>rows:       1350
>cols:       1500
>cells:      2025000

and what is the exact r.width.funct command you used?

see 

https://grass.osgeo.org/grass74/manuals/addons/r.width.funct.html

for an example:

###########
g.region raster=elevation -p

# Calculate flow direction
r.stream.extract elevation=elevation threshold=1000 \
direction=direction 

# Create outlet point 
echo "637304.924954,218617.100523" | v.in.ascii input=- sep=',' out=outlet

# Convert outlet point to raster
v.to.rast input=outlet type=point output=outlet use=cat

# Calculate distance to outlet map
r.stream.distance -o stream_rast=outlet \
direction=direction distance=dist2out

# Calculate width function
r.width.funct map=dist2out image=/tmp/my_basin
###########


kind regards
Helmut


 


More information about the grass-user mailing list