<div dir="ltr">Hi, <br><div><span style="color:rgb(0,0,0);font-family:verdana;font-size:12px">Thanks for your answers, </span>here the command line I run</div><div>r.basin map=Ghiandone_DEM@Silvia prefix=out coordinates=368450.0,4954025.0 dir=Silvia threshold=1000 --overwrite<br></div><div><br></div><div>and these are my input raster info</div><div><br></div><div><div>+----------------------------------------------------------------------------+</div><div> | Map: Ghiandone_DEM@Silvia Date: Fri Sep 02 12:32:06 2016 |</div><div> | Mapset: Silvia Login of Creator: Silvia |</div><div> | Location: Piemonte |</div><div> | DataBase: C:\Users\Silvia\Desktop\Grass |</div><div> | Title: ( Ghiandone_DEM ) |</div><div> | Timestamp: none |</div><div> |----------------------------------------------------------------------------|</div><div> | |</div><div> | Type of Map: raster Number of Categories: 0 |</div><div> | Data Type: CELL |</div><div> | Rows: 888 |</div><div> | Columns: 1047 |</div><div> | Total Cells: 929736 |</div><div> | Projection: UTM (zone 32) |</div><div> | N: 4965768.6398801 S: 4943886.74090327 Res: 24.64177813 |</div><div> | E: 379279.66825865 W: 353479.72655961 Res: 24.64177813 |</div><div> | Range of data: min = 0 max = 2550 |</div><div> | |</div><div> | Data Description: |</div><div> | generato da r.in.gdal |</div><div> | |</div><div> | Comments: |</div><div> | r.in.gdal input="C:\Users\Silvia\Desktop\Grass\Piemonte\Silvia\Ghian\ |</div><div> | done_DEM.tif" output="Ghiandone_DEM" memory=300 offset=0 |</div><div> | |</div><div> +----------------------------------------------------------------------------+</div><div>(Fri Oct 14 10:11:48 2016) Comando terminato (0 sec) </div><div><br></div><div>I tried with Stefan's command line but unfortunately I didn't solve my problem. </div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-10-14 2:45 GMT+02:00 SBL [via OSGeo.org] <span dir="ltr"><<a href="/user/SendEmail.jtp?type=node&node=5290740&i=0" target="_top" rel="nofollow" link="external">[hidden email]</a>></span>:<br><blockquote style='border-left:2px solid #CCCCCC;padding:0 1em' class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="m_388836769193234384WordSection1"><div><div class="h5">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Hi Priscilla,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">I guess the problem is that you run g.region both with the –a flag and the align option.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Se:
<a href="https://grass.osgeo.org/grass70/manuals/g.region.html" rel="nofollow" link="external" target="_blank">https://grass.osgeo.org/<wbr>grass70/manuals/g.region.html</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Please try:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">g.region -p raster=out_Ghiandone_DEM_<wbr>strahler@Silvia align=out_Ghiandone_DEM_<wbr>strahler@Silvia<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">or even just:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">g.region -p raster=out_Ghiandone_DEM_<wbr>strahler@Silvia<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">(both commands should give the same result).<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">I have seen quite some confusion around the –a flag (e.g. it has been used in Processing in QGIS in a unsuitable way for
quite some time).<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Maybe the manual should be improved by:<u></u><u></u></span></p>
<p class="m_388836769193234384MsoListParagraph"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><span>1)<span style="font:7.0pt "Times New Roman"">
</span></span></span><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">altering the flag description and
<u></u><u></u></span></p>
<p class="m_388836769193234384MsoListParagraph"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><span>2)<span style="font:7.0pt "Times New Roman"">
</span></span></span><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">adding examples for aligning the region to
<u></u><u></u></span></p>
<p class="m_388836769193234384MsoListParagraph" style="margin-left:72.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><span>a.<span style="font:7.0pt "Times New Roman"">
</span></span></span><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">another raster map<u></u><u></u></span></p>
<p class="m_388836769193234384MsoListParagraph" style="margin-left:72.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><span>b.<span style="font:7.0pt "Times New Roman"">
</span></span></span><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">a user given resolution / map resolution<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Flag description could be e.g.:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">“Adjusted to be even multiples of the resolution in raster or res option” instead of the current:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">“Align region to resolution (default = align to bounds, works only for 2D resolution)”<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">BTW, is there any case where it makes sense to specify both –a flag and align option? If not, both could be made mutually
exclusive in order to avoid such confusion?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Finally, two possible examples for the manual of g.region, which might help to make that a bit more clear:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"># Set region extent to a vector map and align resolution to an existing raster map
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">g.region –p vector=soils align=elev_state_500m<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"># Define region to align resolution to be an even multiple of the resolution specified by the user
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">g.region –pa vector=soils res=100<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Kind regards,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Stefan<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
</div></div><p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> grass-user [mailto:<a href="http:///user/SendEmail.jtp?type=node&node=5290654&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>]
<b>On Behalf Of </b>Margherita Di Leo<span class=""><br>
<b>Sent:</b> 13. oktober 2016 18:18<br>
</span><b>To:</b> Priscilla_I <<a href="http:///user/SendEmail.jtp?type=node&node=5290654&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
<b>Cc:</b> GRASS user list <<a href="http:///user/SendEmail.jtp?type=node&node=5290654&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>><span class=""><br>
<b>Subject:</b> Re: [GRASS-user] Problem in running r.basin<u></u><u></u></span></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p><div><div class="h5">
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Thu, Oct 13, 2016 at 4:51 PM, Priscilla_I <<a href="http:///user/SendEmail.jtp?type=node&node=5290654&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a>> wrote:<u></u><u></u></p>
<blockquote style='border-left:2px solid #CCCCCC;padding:0 1em' style="border-left:2px solid #cccccc;padding:0 1em">
<p class="MsoNormal">Hi,<br>
I'm a new grass user, running "r.basin" module for Grass 7.0.4 on Windows<br>
10, I get this error message "ERROR: Region resolution and raster map<br>
<out_Ghiandone_DEM_strahler> resolution differs. Run 'g.region<br>
raster=out_Ghiandone_DEM_<wbr>strahler' to set proper region resolution."<br>
I'm using a map in projected coordinates, and I align region with the raster<br>
<br>
g.region -p -a raster=out_Ghiandone_DEM_<wbr>strahler@Silvia<br>
align=out_Ghiandone_DEM_<wbr>strahler@Silvia<br>
<br>
projection: 1 (UTM)<br>
zone: 32<br>
datum: wgs84<br>
ellipsoid: wgs84<br>
north: 4965768.6398801<br>
south: 4943886.74090327<br>
west: 353479.72655961<br>
east: 379279.66825865<br>
nsres: 24.64177813<br>
ewres: 24.64177813<br>
rows: 888<br>
cols: 1047<br>
cells: 929736<br>
<br>
<br>
raster info are<br>
<br>
Map: out_Ghiandone_DEM_strahler@Si Date: Thu Oct 13 11:31:23 2016 |<br>
| Mapset: Silvia Login of Creator: Silvia<br>
|<br>
| Location: Piemonte<br>
|<br>
| DataBase: C:\Users\Silvia\Desktop\Grass<br>
|<br>
| Title: ( out_Ghiandone_DEM_strahler )<br>
|<br>
| Timestamp: none<br>
|<br>
<br>
|-----------------------------<wbr>------------------------------<wbr>-----------------|<br>
|<br>
|<br>
| Type of Map: raster Number of Categories: 0<br>
|<br>
| Data Type: CELL<br>
|<br>
| Rows: 888<br>
|<br>
| Columns: 1047<br>
|<br>
| Total Cells: 929736<br>
|<br>
| Projection: UTM (zone 32)<br>
|<br>
| N: 4965768.6398801 S: 4943886.74090327 Res: 24.64177813<br>
|<br>
| E: 379279.66825865 W: 353479.72655961 Res: 24.64177813<br>
|<br>
| Range of data: min = 1 max = 4<br>
|<br>
|<br>
|<br>
| Data Description:<br>
|<br>
| generato da r.stream.order<br>
|<br>
|<br>
|<br>
| Comments:<br>
|<br>
| r.stream.order stream_rast="out_Ghiandone_<wbr>DEM_stream_e" direction="o\<br>
|<br>
| ut_Ghiandone_DEM_drainage_e" strahler="out_Ghiandone_DEM_<wbr>strahler" h\<br>
|<br>
| orton="out_Ghiandone_DEM_<wbr>horton" shreve="out_Ghiandone_DEM_<wbr>shreve" h\<br>
|<br>
| ack="out_Ghiandone_DEM_hack" memory=300<br>
|<br>
|<br>
|<br>
<br>
+-----------------------------<wbr>------------------------------<wbr>-----------------+<br>
<br>
Any ideas? Thanks in advance!<u></u><u></u></p>
</blockquote>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">What is the exact r.basin command line that you run? Can you post also
<a href="http://r.info" rel="nofollow" link="external" target="_blank">r.info</a> of the DEM that you use as input? <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">-- <u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#666666">Margherita Di Leo</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
</div>
<br><span class="">______________________________<wbr>_________________
<br>grass-user mailing list
<br><a href="http:///user/SendEmail.jtp?type=node&node=5290654&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="http://lists.osgeo.org/mailman/listinfo/grass-user" rel="nofollow" link="external" target="_blank">http://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a>
<br>
<br>
<hr noshade size="1" color="#cccccc">
</span><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><span class="">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</span><a href="http://osgeo-org.1560.x6.nabble.com/Problem-in-running-r-basin-tp5290587p5290654.html" target="_blank" rel="nofollow" link="external">http://osgeo-org.1560.x6.<wbr>nabble.com/Problem-in-running-<wbr>r-basin-tp5290587p5290654.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from Problem in running r.basin, <a href="" target="_blank" rel="nofollow" link="external">click here</a>.<br>
<a href="http://osgeo-org.1560.x6.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank" link="external">NAML</a>
</div></div></div></blockquote></div><br></div></div>
<br/><hr align="left" width="300" />
View this message in context: <a href="http://osgeo-org.1560.x6.nabble.com/Problem-in-running-r-basin-tp5290587p5290740.html">Re: Problem in running r.basin</a><br/>
Sent from the <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html">Grass - Users mailing list archive</a> at Nabble.com.<br/>