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@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@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@hydrology accumulation=fac@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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user