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

[Jarek Jasiewicz]

On 06/27/2011 01:44 PM, jentjr wrote:
> 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.
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
it looks like the accum map has bigger extend than DEM. Try to cut both 
maps to the same extend. It should help