<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>Hi Helmut and Margherita,</div>
<div> </div>
<div>Thanks for your mails! Here are the answers to your questions:</div>
<div> </div>
<div><strong>a) </strong><span style="line-height: 1.6em;"><strong>This is the result of r.info and g.region: </strong> </span></div>
<div> </div>
<div>r.info map@Elevation </div>
<div> +----------------------------------------------------------------------------+</div>
<div> | Map: map@Elevation Date: Sat Nov 15 16:35:08 2014 |</div>
<div> | Mapset: Elevation Login of Creator: Andrea |</div>
<div> | Location: ElevationData |</div>
<div> | DataBase: C:\Users\Andrea\Documents\grassdata |</div>
<div> | Title: ( map ) |</div>
<div> | Timestamp: none |</div>
<div> |----------------------------------------------------------------------------|</div>
<div> | |</div>
<div> | Type of Map: raster Number of Categories: 0 |</div>
<div> | Data Type: FCELL |</div>
<div> | Rows: 21612 |</div>
<div> | Columns: 10812 |</div>
<div> | Total Cells: 233668944 |</div>
<div> | Projection: Latitude-Longitude |</div>
<div> | N: 45:00:02N S: 42:59:58N Res: 0:00:00.333333 |</div>
<div> | E: 70:59:58W W: 72:00:02W Res: 0:00:00.333333 |</div>
<div> | Range of data: min = 18.40538 max = 1916.398 |</div>
<div> | |</div>
<div> | Data Description: |</div>
<div> | generated by r.patch |</div>
<div> | |</div>
<div> | Comments: |</div>
<div> | r.patch input="n44-72@Elevation,n45_72@Elevation" output="map" |</div>
<div> | |</div>
<div> +----------------------------------------------------------------------------+</div>
<div> </div>
<div>g.region -p </div>
<div>projection: 3 (Latitude-Longitude)</div>
<div>zone: 0</div>
<div>datum: nad83</div>
<div>ellipsoid: grs80</div>
<div>north: 45:00:02N</div>
<div>south: 42:59:58N</div>
<div>west: 72:00:02W</div>
<div>east: 70:59:58W</div>
<div>nsres: 0:00:01</div>
<div>ewres: 0:00:01</div>
<div>rows: 7204</div>
<div>cols: 3604</div>
<div>cells: 25963216</div>
<div> </div>
<div><strong>b) These are all my steps:</strong></div>
<div>++ open grass</div>
<div>++ load my two raster files</div>
<div>++ r.patch --overwrite input=n44_72@Elevation,n45_72@Elevation output=map</div>
<div>++ g.remove -f type=all name=original@Elevation,o_map_accumulation@Elevation,o_map_aspect@Elevation,o_map_aspect_mod@Elevation,o_map_average_hillslope@Elevation,o_map_basin@Elevation,o_map_dist2out@Elevation,o_map_drainage@Elevation,o_map_drainage_e@Elevation,o_map_hack@Elevation,o_map_height_average@Elevation,o_map_hillslope_distance@Elevation,o_map_horton@Elevation,o_map_mainchannel@Elevation,o_map_mainchannel_dim@Elevation,o_map_mainchannel_thin@Elevation,o_map_mask@Elevation,o_map_ord_1@Elevation,o_map_ord_1_thin@Elevation,o_map_r_outlet@Elevation,o_map_shreve@Elevation,o_map_slope@Elevation,o_map_strahler@Elevation,o_map_stream_e@Elevation,o_map_stream_e_thin@Elevation,r_elevation_crop@Elevation,o_map_basin@Elevation,o_map_mainchannel@Elevation,o_map_mainchannel_dim@Elevation,o_map_mainchannel_dim_point@Elevation,o_map_network@Elevation,o_map_ord_1@Elevation,o_map_outlet@Elevation,o_map_outlet_snap@Elevation,o_map_mainchannel_dim_thin@Elevation</div>
<div> </div>
<div>++ g.region rast=map@Elevation res=0:00:01</div>
<div>(also tried with g.region <strong>-a</strong> rast=map@Elevation res=0:00:01, and got the same result)</div>
<div> </div>
<div>++ r.basin map=map@Elevation prefix=o coordinates=-71.10394196,43.9865230801 threshold=19005 dir="C:\\Users\\Andrea\\Basins5"</div>
<div> </div>
<div><strong>c) Try to use the NC dataset</strong></div>
<div>I downloaded it, but did not any changes in the resolution and computational region, since the manual said that the maps were ready for use...</div>
<div>I chose some coordinates just by a visual inspection of the DEM and used just random threshold. I run:</div>
<div> </div>
<div>r.basin map=elev_ned_30m@PERMANENT prefix=o coordinates=640856.761198,215050.690725 threshold=400 dir="C:\\Users\\Andrea\\Basins6</div>
<div> </div>
<div>and get:</div>
<div>“<em>...</em></div>
<div><em>...</em></div>
<div><em>Tot. cells 543.0</em></div>
<div><em>===========================</em></div>
<div><em>Hypsometric | quantiles</em></div>
<div><em>===========================</em></div>
<div><em>106 | 0.025</em></div>
<div><em>105 | 0.05</em></div>
<div><em>103 | 0.1</em></div>
<div><em>100 | 0.25</em></div>
<div><em>95 | 0.5</em></div>
<div><em>87 | 0.75</em></div>
<div><em>89 | 0.7</em></div>
<div><em>81 | 0.9</em></div>
<div><em>76 | 0.975</em></div>
<div><em>Done!</em></div>
<div><em>------------------------------</em></div>
<div><em>------------------------------</em></div>
<div><em>Traceback (most recent call last):</em></div>
<div><em> File "C:\Users\Andrea\AppData\Roaming\GRASS7\addons/scri</em></div>
<div><em>pts/r.width.funct.py", line 130, in <module></em></div>
<div><em> sys.exit(main())</em></div>
<div><em> File "C:\Users\Andrea\AppData\Roaming\GRASS7\addons/scri</em></div>
<div><em>pts/r.width.funct.py", line 86, in main</em></div>
<div><em> prc[4,0] , prc[4,1] = findint(kl,0.5) , 0.5</em></div>
<div><em> File "C:\Users\Andrea\AppData\Roaming\GRASS7\addons/scri</em></div>
<div><em>pts/r.width.funct.py", line 123, in findint</em></div>
<div><em> z1 , z2 , f1 , f2 = kl[float(Xf[0])][0] ,</em></div>
<div><em>kl[float(Xf[0]-1)][0] , kl[float(Xf[0])][1] ,</em></div>
<div><em>kl[float(Xf[0]-1)][1]</em></div>
<div><em>TypeError: only length-1 arrays can be converted to Python</em></div>
<div><em>scalars</em></div>
<div><em>Tot. cells 543.0</em></div>
<div><em>Tot. area 441051.75</em></div>
<div><em>Max distance 1693.190486</em></div>
<div>
<div><br/>
<em>------------------------------</em></div>
<div><em>An ERROR occurred running r.basin<br/>
Please check for error messages above or try with another pairs of outlet coordinates"</em></div>
<div> </div>
<div> </div>
</div>
<div><strong>d) Results of r.info and g.region for the NC dataset</strong></div>
<div> </div>
<div>r.info elev_ned_30m@PERMANENT </div>
<div> +----------------------------------------------------------------------------+</div>
<div> | Map: elev_ned_30m@PERMANENT Date: Tue Nov 7 00:35:18 2006 |</div>
<div> | Mapset: PERMANENT Login of Creator: helena |</div>
<div> | Location: nc_spm_08_grass7 |</div>
<div> | DataBase: C:\Users\Andrea\Documents\grassdata\nc_spm_08_grass7 |</div>
<div> | Title: South-West Wake county: National Elevation Data 30m ( elev_ned30 |</div>
<div> | Timestamp: none |</div>
<div> |----------------------------------------------------------------------------|</div>
<div> | |</div>
<div> | Type of Map: raster Number of Categories: 255 |</div>
<div> | Data Type: FCELL |</div>
<div> | Rows: 450 |</div>
<div> | Columns: 500 |</div>
<div> | Total Cells: 225000 |</div>
<div> | Projection: Lambert Conformal Conic |</div>
<div> | N: 228500 S: 215000 Res: 30 |</div>
<div> | E: 645000 W: 630000 Res: 30 |</div>
<div> | Range of data: min = 55.1736 max = 156.3865 |</div>
<div> | |</div>
<div> | Data Description: |</div>
<div> | generated by r.proj |</div>
<div> | |</div>
<div> | Comments: |</div>
<div> | r.proj input="ned1arcsec" location="northcarolina_latlong" mapset="h\ |</div>
<div> | elena" output="elev_ned30m" method="cubic" resolution=30 |</div>
<div> | |</div>
<div> +----------------------------------------------------------------------------+</div>
<div> </div>
<div> </div>
<div>g.region -p </div>
<div>projection: 99 (Lambert Conformal Conic)</div>
<div>zone: 0</div>
<div>datum: nad83</div>
<div>ellipsoid: a=6378137 es=0.006694380022900787</div>
<div>north: 228513</div>
<div>south: 214975.5</div>
<div>west: 629992.5</div>
<div>east: 645012</div>
<div>nsres: 28.5</div>
<div>ewres: 28.5</div>
<div>rows: 475</div>
<div>cols: 527</div>
<div>cells: 250325</div>
<div>
<div> </div>
<div> </div>
<div><strong style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;">e) The DEMs I am using</strong></div>
<div><span style="line-height: 1.6em;">They are from the NED webpage (USGS), here is the link they sent me for downloading their files, but they are kind of large. If you want to download them I could select a catchment which only uses one DEM ...</span></div>
<div> </div>
<div>
<div><a href="https://3c.gmx.net/mail/client/dereferrer?redirectUrl=ftp%3A%2F%2Frockyftp.cr.usgs.gov%2Fvdelivery%2FDatasets%2FStaged%2FNED%2F13%2FArcGrid%2Fn45w072.zip">https://3c.gmx.net/mail/client/dereferrer?redirectUrl=ftp%3A%2F%2Frockyftp.cr.usgs.gov%2Fvdelivery%2FDatasets%2FStaged%2FNED%2F13%2FArcGrid%2Fn45w072.zip</a></div>
<div> </div>
<div><a href="https://3c.gmx.net/mail/client/dereferrer?redirectUrl=ftp%3A%2F%2Frockyftp.cr.usgs.gov%2Fvdelivery%2FDatasets%2FStaged%2FNED%2F13%2FArcGrid%2Fn44w072.zip">https://3c.gmx.net/mail/client/dereferrer?redirectUrl=ftp%3A%2F%2Frockyftp.cr.usgs.gov%2Fvdelivery%2FDatasets%2FStaged%2FNED%2F13%2FArcGrid%2Fn44w072.zip</a></div>
</div>
<div> </div>
<div>Hope this helps, </div>
<div>Regards, Andrea. </div>
<div> </div>
<div style="margin: 10.0px 5.0px 5.0px 10.0px;padding: 10.0px 0 10.0px 10.0px;border-left: 2.0px solid rgb(195,217,229);">
<div style="margin: 0 0 10.0px 0;"><b>Gesendet:</b> Dienstag, 18. November 2014 um 16:54 Uhr<br/>
<b>Von:</b> "Helmut Kudrnovsky" <hellik@web.de><br/>
<b>An:</b> grass-user@lists.osgeo.org<br/>
<b>Betreff:</b> Re: [GRASS-user] Problem with r.basin in grass7</div>
<div>Andrea Timmermann wrote<br/>
> Hi Margherita, <br/>
><br/>
> <br/>
><br/>
> Thanks for your answer. I am sorry I did not explain well what I<br/>
> meant. <br/>
><br/>
> <br/>
><br/>
> When I delete all files and start again, I get again stuck at that point.<br/>
><br/>
> <br/>
><br/>
> What I do is: delete all files, set region, run r.basin, get the error ...<br/>
> detele all files, set the region, run r.basin, get the error ....<br/>
> <br/>
><br/>
><br/>
> <br/>
><br/>
> I just tried it again ..<br/>
><br/>
> <br/>
><br/>
> Thanks, <br/>
><br/>
> Andrea<br/>
><br/>
><br/>
> Gesendet: Dienstag, 18. November 2014 um 09:38 Uhr<br/>
> Von: "Margherita Di Leo" <<br/>
<br/>
> diregola@<br/>
<br/>
> ><br/>
> An: "Andrea Timmermann" <<br/>
<br/>
> timmermann@<br/>
<br/>
> ><br/>
> Cc: "Helmut Kudrnovsky" <<br/>
<br/>
> hellik@<br/>
<br/>
> >, "GRASS user list" <<br/>
<br/>
> grass-user@.osgeo<br/>
<br/>
> ><br/>
> Betreff: Re: [GRASS-user] Problem with r.basin in grass7<br/>
><br/>
><br/>
> Hi Andrea,<br/>
> <br/>
> On Tue, Nov 18, 2014 at 3:55 AM, Andrea Timmermann <<br/>
<br/>
> timmermann@<br/>
<br/>
> > wrote:<br/>
><br/>
><br/>
><br/>
><br/>
> Hi Helmut, <br/>
><br/>
> <br/>
><br/>
> Thanks for your email. <br/>
><br/>
> I set again the area and resolution with g.region before running r.basin.<br/>
> I also deleted all files except my DEM. <br/>
><br/>
> <br/>
><br/>
> Now I get a different error ...<br/>
><br/>
> <br/>
><br/>
> When I run<br/>
><br/>
><br/>
> "r.basin map=map@Elevation prefix=o<br/>
> coordinates=-71.10394196,43.9865230801 threshold=19005<br/>
> dir=C:\\Users\\Andrea\\Basins"<br/>
><br/>
> <br/>
><br/>
> I get:<br/>
><br/>
> "...<br/>
><br/>
> Building topology for vector map <o_map_ord_1@Elevation>...<br/>
> Registering primitives...<br/>
> 10 primitives registered<br/>
> 558 vertices registered<br/>
> Building areas...<br/>
> 0 areas built<br/>
> 0 isles built<br/>
> Attaching islands...<br/>
> Attaching centroids...<br/>
> Number of nodes: 20<br/>
> Number of primitives: 10<br/>
> Number of points: 0<br/>
> Number of lines: 10<br/>
> Number of boundaries: 0<br/>
> Number of centroids: 0<br/>
> Number of areas: 0<br/>
> Number of isles: 0<br/>
> Updating attributes...<br/>
> r.to.vect complete.<br/>
> All in RAM calculation...<br/>
> ERROR: Region resolution and raster map <o_map_strahler> resolution<br/>
> differs. Run 'g.region rast=o_map_strahler' to set proper region<br/>
> resolution"<br/>
><br/>
> <br/>
><br/>
> I run g.region rast=o_map_strahler, but then I have to delete<br/>
> all files and start over again running r.basin.<br/>
><br/>
><br/>
><br/>
><br/>
><br/>
> <br/>
><br/>
> Yes if you change the region you have to delete the files created in the<br/>
> previous run (if they have same name) or allow overwrite.<br/>
><br/>
><br/>
> <br/>
> --<br/>
><br/>
><br/>
><br/>
> Best regards,<br/>
><br/>
> <br/>
><br/>
> Dr. Margherita DI LEO <br/>
><br/>
> Scientific / technical project officer<br/>
><br/>
> <br/>
><br/>
> European Commission - DG JRC <br/>
><br/>
> Institute for Environment and Sustainability (IES)<br/>
><br/>
> Via Fermi, 2749<br/>
><br/>
> I-21027 Ispra (VA) - Italy - TP 261<br/>
><br/>
> <br/>
><br/>
> Tel. +39 0332 78 3600 <br/>
<br/>
> margherita.di-leo@.europa<br/>
<br/>
><br/>
><br/>
> <br/>
><br/>
> Disclaimer: The views expressed are purely those of the writer and may not<br/>
> in any circumstance be regarded as stating an official position of the<br/>
> European Commission.<br/>
><br/>
><br/>
><br/>
><br/>
><br/>
><br/>
><br/>
><br/>
> _______________________________________________<br/>
> grass-user mailing list<br/>
<br/>
> grass-user@.osgeo<br/>
<br/>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br/>
<br/>
could you post also the exact commands of your steps?<br/>
e.g.<br/>
<br/>
g.region -a -p rast=yourdem align=yourdem<br/>
r.basin --verbose elevation=yourdem....<br/>
<br/>
g.region<br/>
<br/>
<br/>
<br/>
<br/>
-----<br/>
best regards<br/>
Helmut<br/>
--<br/>
View this message in context: <a href="http://osgeo-org.1560.x6.nabble.com/Problem-with-r-basin-in-grass7-tp5169155p5173590.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/Problem-with-r-basin-in-grass7-tp5169155p5173590.html</a><br/>
Sent from the Grass - Users mailing list archive at Nabble.com.<br/>
_______________________________________________<br/>
grass-user mailing list<br/>
grass-user@lists.osgeo.org<br/>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a></div>
</div>
</div></div></body></html>