[GRASS-user] Re: r.stream.extract with external flow accumulation maps

jentjr jentjr at gmail.com
Mon Jun 27 13:44:56 EDT 2011


Thank you for the replies. I set the region to the flow accumulation map, and
I still received the error though. Here is the region output for both the
elevation and flow accumulation maps

Elevation:
 
GRASS 6.4.1 (NAD83_GRS80_LCC):~ > g.region -p rast=elev_cm at hydrology
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      276735
south:      182265
west:       111075
east:       398265
nsres:      30
ewres:      30
rows:       3149
cols:       9573
cells:      30145377

Flow accumulation:

GRASS 6.4.1 (NAD83_GRS80_LCC):~ > g.region -p rast=fac at hydrology
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      276705
south:      182025
west:       110445
east:       399105
nsres:      30
ewres:      30
rows:       3156
cols:       9622
cells:      30367032


I also tried setting a mask 

r.mapcalc 'new_fac=if(fac==null(),null(),1)'
r.mask input=new_fac

GRASS 6.4.1 (NAD83_GRS80_LCC):~ > r.stream.extract
elevation=elev_cm at hydrology accumulation=fac at hydrology stream_rast=rstream5K
stream_vect=rstream5K threshold=5000 direction=rstream5K_fdr

ERROR: Accumulation map does not match elevation map!
[Raster MASK present]




--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-stream-extract-with-external-flow-accumulation-maps-tp6518518p6521518.html
Sent from the Grass - Users mailing list archive at Nabble.com.


More information about the grass-user mailing list