[GRASS-user] Problem with r.basin in grass7

[Helmut Kudrnovsky]

> r.basin map=elev_ned_30m at PERMANENT <=

>My result of g.region -p is:
>projection: 99 (Lambert Conformal Conic)
>zone:       0
>datum:      nad83
>ellipsoid:  a=6378137 es=0.006694380022900787
>north:      228513
>south:      214975.5
>west:       629992.5
>east:       645012
>nsres:      28.5 <=
>ewres:      28.5 <=
>rows:       475
>cols:       527
>cells:      250325 

are you setting your region to match elev_ned_30m at PERMANENT? 

here: 

r.info map=elev_ned_30m at PERMANENT                                               

+----------------------------------------------------------------------------+
 | Map:      elev_ned_30m at PERMANENT         Date: Tue Nov  7 00:35:18 2006   
|
 | Mapset:   PERMANENT                      Login of Creator: helena         
|
 | Location: nc_spm_08_grass7                                                
|
 | DataBase: C:\grassdata                                                    
|
 | Title:    South-West Wake county: National Elevation Data 30m (
elev_ned30 |
 | Timestamp: none                                                           
|

|----------------------------------------------------------------------------|
 |                                                                           
|
 |   Type of Map:  raster               Number of Categories: 255            
|
 |   Data Type:    FCELL                                                     
|
 |   Rows:         450                                                       
|
 |   Columns:      500                                                       
|
 |   Total Cells:  225000                                                    
|
 |        Projection: Lambert Conformal Conic                                
|
 |            N:     228500    S:     215000   Res:    30 <=                    
|
 |            E:     645000    W:     630000   Res:    30 <=                 
|
 |   Range of data:    min = 55.1736  max = 156.3865                         
|
 |                                                                           
|
 |   Data Description:                                                       
|
 |    generated by r.proj                                                    
|
 |                                                                           
|
 |   Comments:                                                               
|
 |    r.proj input="ned1arcsec" location="northcarolina_latlong" mapset="h\  
|
 |    elena" output="elev_ned30m" method="cubic" resolution=30               
|
 |                                                                           
|

+----------------------------------------------------------------------------+

nsres: 28.5 vs 30?

please copy/paste this command:

g.region -a raster=elev_ned_30m at PERMANENT align=elev_ned_30m at PERMANENT

this means: matching region extent to raster DEM _and_ align cells to raster
DEM

here I get:

g.region -p -a raster=elev_ned_30m at PERMANENT align=elev_ned_30m at PERMANENT       
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  a=6378137 es=0.006694380022900787
north:      228500
south:      215000
west:       630000
east:       645000
nsres:      30 <=
ewres:      30 <=
rows:       450
cols:       500
cells:      225000

tested your command here on my side:

r.basin map=elev_ned_30m at PERMANENT prefix=b
coordinates=636654.791181,218824.126649 dir=C:\tmp\test threshold=40
[...]
(Mon Dec 15 09:13:17 2014) Befehl ausgeführt (4 Min 53 Sek)

r.basin finishes without any error.

>ERROR: Error while executing: 'DROP TABLE summary'
>Traceback (most recent call last): 
>[...]

no idea about this one ...




-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Problem-with-r-basin-in-grass7-tp5169155p5177858.html
Sent from the Grass - Users mailing list archive at Nabble.com.