<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>Hi ! </div>
<div> </div>
<div>I downloaded the 5-Nov Grass and also updated all the addons extensions I use.</div>
<div> </div>
<div>When I run:</div>
<div>"<em>r.basin map=map@Elevation prefix=o coordinates=-71.10394196,43.9865230801 threshold=19005 dir="C:\\Users\\Andrea\\Basins</em>"</div>
<div>
<div> </div>
<div>I get:</div>
<div>"<em>...</em></div>
<div>
<div><em>Building topology for vector map <o_map_outlet@Elevation>...<br/>
Registering primitives...<br/>
One primitive registered<br/>
One vertex registered<br/>
Building areas...<br/>
0 areas built<br/>
0 isles built<br/>
Attaching islands...<br/>
Attaching centroids...<br/>
Number of nodes: 0<br/>
Number of primitives: 1<br/>
Number of points: 1<br/>
Number of lines: 0<br/>
Number of boundaries: 0<br/>
Number of centroids: 0<br/>
Number of areas: 0<br/>
Number of isles: 0</em></div>
<div><em>------------------------------</em></div>
<div><em>An ERROR occurred running r.basin<br/>
Please try with another pairs of outlet coordinates</em>"</div>
<div> </div>
<div>I tried to check the return code of the module as mentioned by Markus:</div>
<div><span style="line-height: 19.2000007629395px;">> This is not exactly true. You will be told where in the source code</span></div>
<div><span style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;">> memory allocation failed. This can be somewhere in a library, then you</span><br style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;"/>
<span style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;">> don't know which module actually failed. You need to check the return</span><br style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;"/>
<span style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;">> code of the module, as of r62566 automatically done by</span><br style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;"/>
<span style="font-family: Verdana; font-size: 12px; line-height: 19.2000007629395px;">> grass.run_command in trunk.</span></div>
<div> </div>
<div>but I could not find the file. Can you tell me where I can find it?</div>
<div> </div>
<div> </div>
<div>I have now another question ...</div>
<div>I was thinking about the note on the wiki page which said that:"... <span style="color: rgb(37, 37, 37); font-family: sans-serif; font-size: 13.8181819915771px; line-height: 16.2909088134766px;">perimeter is affected by a fractal problem and the number is mostly meaningless taken by itself. Hence, perimeter should never be mixed with anything from another dataset or even the same data set at a different cell resolution; it is only good as a percentage comparison to its neighbors in the same dataset.</span></div>
<div> </div>
<div>And about the known issues at the same page: "<span style="color: rgb(37, 37, 37); font-family: sans-serif; font-size: 13.8181819915771px; line-height: 16.2909088134766px;">r.basin hasn't been designed for working in lat/long coordinates. This means that if you are working in lat/long coordinates, you need to reproject your map first in order to apply the tool".</span></div>
<div> </div>
<div><span style="color: rgb(37, 37, 37); font-family: sans-serif; font-size: 13.8181819915771px; line-height: 16.2909088134766px;">The DEMs I am using have a resolution of 0.33 degrees and when I reproject them (</span>projection: 3 (Latitude-Longitude); zone: 0; datum: nad83; ellipsoid: grs80), the areas covered by the DEMs files vary depending on the longitude: the surface covered by one file is larger near the equator than closer to the poles. So I think this means that the cell resolution of my reprojected map is larger near the poles than at the equator? If so, it would only make sense to compare the morphometric parameters for catchments at a similar latitude I think... Am I right? I am a bit confused with this ...</div>
<div> </div>
<div>Thanks for your answers!</div>
<div>Regards, Andrea.</div>
<div> </div>
<div> </div>
</div>
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b> Mittwoch, 05. November 2014 um 09:55 Uhr<br/>
<b>Von:</b> "Markus Metz" <markus.metz.giswork@gmail.com><br/>
<b>An:</b> "Markus Neteler" <neteler@osgeo.org><br/>
<b>Cc:</b> "Andrea Timmermann" <timmermann@gmx.at>, "GRASS user list" <grass-user@lists.osgeo.org><br/>
<b>Betreff:</b> Re: [GRASS-user] Problem with r.basin in grass7</div>
<div name="quoted-content">On Tue, Oct 28, 2014 at 5:13 PM, Markus Neteler <neteler@osgeo.org> wrote:<br/>
> Hi Andrea,<br/>
><br/>
> On Fri, Oct 24, 2014 at 3:27 PM, Andrea Timmermann <timmermann@gmx.at> wrote:<br/>
>> Hi Markus,<br/>
>><br/>
>> Thanks for your answer! I have now however another problem ...<br/>
>> When running:<br/>
>><br/>
>> r.basin map=map@Elevation prefix=o coordinates=-71.10394196,43.9865230801<br/>
>> threshold=19005 dir=C:\\Users\\Andrea\\Basins<br/>
>><br/>
>> I get:<br/>
>> " ....<br/>
>> v.to.rast complete.<br/>
>> All in RAM calculation...<br/>
>> Reading raster map <o_map_drainage_e>...<br/>
>> Calculating basins using vector point map...<br/>
>> Current region rows: 21612, cols: 21612<br/>
>> ERROR: G_calloc: unable to allocate 21612 * 4 bytes of memory at io.c:41<br/>
><br/>
> So, starting with GRASS 7 revision r62444 (see startup screen or<br/>
> g.version -g) you will be told *which* module actually failed,<br/>
<br/>
This is not exactly true. You will be told where in the source code<br/>
memory allocation failed. This can be somewhere in a library, then you<br/>
don't know which module actually failed. You need to check the return<br/>
code of the module, as of r62566 automatically done by<br/>
grass.run_command in trunk.<br/>
<br/>
Markus M<br/>
<br/>
> so the<br/>
> G_calloc error will be more useful. To get that, you will need to<br/>
> update your GRASS GIS 7 installation.<br/>
><br/>
> However, a few thoughts before:<br/>
><br/>
> Current region rows: 21612, cols: 21612<br/>
> i.e. 467078544 cells at 24 bytes/cell (AFAIK). That is:<br/>
><br/>
>> 467078544 * 24 / 1024 / 1024/ 1024<br/>
> [1] 10.44002<br/>
><br/>
> Leaving a bit for the operating system etc, you may exceed 12GB RAM if<br/>
> I got the numbers right.<br/>
><br/>
> Are you sure that you want to use 467 million cells? If not, check<br/>
> "g.region -p", perhaps the computational region resolution is too<br/>
> high?<br/>
><br/>
> best<br/>
> Markus<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>