<div dir="ltr">Hi,<br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 14, 2016 at 2:33 PM, Priscilla_I <span dir="ltr"><<a href="mailto:isacco.silvia@gmail.com" target="_blank">isacco.silvia@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote><div><br></div><div>I think the problem lies in the "overwrite", which is not handled by r.basin. Try to remove all your output map, then set the region on the DEM and then run r.basin, and let us know.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div><br></div><div>and these are my input raster info</div><div><br></div><div><div>+-----------------------------<wbr>------------------------------<wbr>-----------------+</div><div> | Map:      Ghiandone_DEM@Silvia           Date: Fri Sep 02 12:32:06 2016    |</div><span class=""><div> | Mapset:   Silvia                         Login of Creator: Silvia          |</div><div> | Location: Piemonte                                                         |</div><div> | DataBase: C:\Users\Silvia\Desktop\Grass                                    |</div></span><div> | Title:     ( Ghiandone_DEM )                                               |</div><span class=""><div> | Timestamp: none                                                            |</div><div> |----------------------------<wbr>------------------------------<wbr>------------------|</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></span><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\<wbr>Desktop\Grass\Piemonte\Silvia\<wbr>Ghian\   |</div><div> |    done_DEM.tif" output="Ghiandone_DEM" memory=300 offset=0                |</div><div> |                                                                            |</div><div> +----------------------------<wbr>------------------------------<wbr>------------------+</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"><div><div class="h5">2016-10-14 2:45 GMT+02:00 SBL [via OSGeo.org] <span dir="ltr"><<a href="http:///user/SendEmail.jtp?type=node&node=5290740&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span>:<br></div></div><blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote">

        






<div class="m_-7075550699890084018m_388836769193234384WordSection1"><div><div class="h5"><div><div class="m_-7075550699890084018h5">
<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/grass7<wbr>0/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_strah<wbr>ler@Silvia align=out_Ghiandone_DEM_strahl<wbr>er@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_strah<wbr>ler@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_-7075550699890084018m_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_-7075550699890084018m_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_-7075550699890084018m_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_-7075550699890084018m_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></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=""><span><br>
<b>Sent:</b> 13. oktober 2016 18:18<br>
</span></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=""><span><br>
<b>Subject:</b> Re: [GRASS-user] Problem in running r.basin<u></u><u></u></span></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><div class="m_-7075550699890084018h5">
<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">
<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_strah<wbr>ler' 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_strah<wbr>ler@Silvia<br>
align=out_Ghiandone_DEM_strahl<wbr>er@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_DEM<wbr>_stream_e" direction="o\<br>
|<br>
 |    ut_Ghiandone_DEM_drainage_e" strahler="out_Ghiandone_DEM_st<wbr>rahler" h\<br>
|<br>
 |    orton="out_Ghiandone_DEM_horto<wbr>n" shreve="out_Ghiandone_DEM_shre<wbr>ve" 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></div>
</div>


<br><span>______________________________<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/mailman<wbr>/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>
                <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" rel="nofollow" link="external" target="_blank">http://osgeo-org.1560.x6.nabbl<wbr>e.com/Problem-in-running-r-<wbr>basin-tp5290587p5290654.html</a>
        </div><div class="m_-7075550699890084018HOEnZb"><div class="m_-7075550699890084018h5">
        <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 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" link="external" target="_blank">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" target="_blank">Re: Problem in running r.basin</a><div class="HOEnZb"><div class="h5"><br>
Sent from the <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html" target="_blank">Grass - Users mailing list archive</a> at Nabble.com.<br></div></div><br>______________________________<wbr>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">http://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><font color="#666666">Margherita Di Leo</font></div></div></div></div></div></div></div></div></div>
</div></div>