Hi,<div><br></div><div>r.thin fails, most likely because there's nothing to thin. It is supposed to thin the stream network extracted by r.stream.extract, please take a look at the map produced. Most likely it will be empty, because the coordinates of the outlet don't match to the stream.</div>
<div><br></div><div>HTH</div><div>madi<br><br><div class="gmail_quote">On Tue, Jan 29, 2013 at 8:45 AM, Giuliano Urgeghe <span dir="ltr"><<a href="mailto:giuliano.u@gmail.com" target="_blank">giuliano.u@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I think I have solved the error accumulation map does not match with the<br>
elevation map installing 6.4.3<br>
but now I get this error:<br>
Tue Jan 29 08:32:15 2013)<br>
r.basin map=dem@PERMANENT prefix=LuLitarroni61 easting=1509351.753173<br>
northing=4553749.816433 threshold=100<br>
SECTION 1a (of 4): Initiating Memory.<br>
SECTION 1b (of 4): Determining Offmap Flow.<br>
SECTION 2: A * Search.<br>
SECTION 3: Accumulating Surface Flow with SFD.<br>
SECTION 4: Closing Maps.<br>
WARNING: Writing out only positive flow accumulation values.<br>
WARNING: Cells with a likely underestimate for flow accumulation can no<br>
longer be identified.<br>
Load input maps and get start points...<br>
Set edge points...<br>
A* Search...<br>
Extract streams...<br>
Thin stream segments...<br>
Write raster maps<br>
Reading maps...<br>
Calculate basins using coordinates...<br>
LuLitarroni61_dem_basin Done!<br>
Delineation of basin done<br>
Rename raster <tmp> to <LuLitarroni61_dem_accumulation><br>
Rename raster <tmp> to <LuLitarroni61_dem_drainage><br>
Rename raster <tmp> to <LuLitarroni61_dem_drainage_e><br>
Rename raster <tmp> to <LuLitarroni61_dem_stream_e><br>
File LuLitarroni61_dem_stream_e -- 2426 rows X 3003 columns<br>
ERROR: (null): Unable to find bounding box for lines<br>
ERROR: Raster map <LuLitarroni61_dem_stream_e_thin> not found<br>
Aspect raster map <LuLitarroni61_dem_aspect> complete<br>
Slope raster map <LuLitarroni61_dem_slope> complete<br>
Extracting areas...<br>
Building topology for vector map <LuLitarroni61_dem_basin>...<br>
Registering primitives...<br>
37 primitives registered<br>
189 vertices registered<br>
Building areas...<br>
13 areas built<br>
1 isles built<br>
Attaching islands...<br>
Attaching centroids...<br>
Number of nodes: 25<br>
Number of primitives: 37<br>
Number of points: 0<br>
Number of lines: 0<br>
Number of boundaries: 24<br>
Number of centroids: 13<br>
Number of areas: 13<br>
Number of isles: 1<br>
r.to.vect complete.<br>
Reading areas...<br>
Updating database...<br>
1 categories read from vector map (layer 1)<br>
1 records selected from table (layer 1)<br>
1 categories read from vector map exist in selection from table<br>
1 records updated/inserted (layer 1)<br>
Reading areas...<br>
1 categories read from vector map (layer 1)<br>
0 records updated/inserted (layer 1)<br>
Reading areas...<br>
Updating database...<br>
1 categories read from vector map (layer 1)<br>
1 records selected from table (layer 1)<br>
1 categories read from vector map exist in selection from table<br>
1 records updated/inserted (layer 1)<br>
Reading areas...<br>
1 categories read from vector map (layer 1)<br>
0 records updated/inserted (layer 1)<br>
Creating LuLitarroni61_dem_hack<br>
Reading maps...<br>
ERROR: No streams found<br>
Scanning input for column types...<br>
Maximum input row length: 33<br>
Maximum number of columns: 3<br>
Minimum number of columns: 3<br>
Column: 1 type: double<br>
Column: 2 type: double<br>
Column: 3 type: integer<br>
Importing points...<br>
Populating table...<br>
Building topology for vector map <LuLitarroni61_dem_outlet>...<br>
Registering primitives...<br>
1 primitives registered<br>
1 vertices registered<br>
Building areas...<br>
0 areas built<br>
0 isles built<br>
Attaching islands...<br>
Attaching centroids...<br>
Number of nodes: 1<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<br>
v.in.ascii complete.<br>
Loading data...<br>
Reading features...<br>
Writing raster map...<br>
Converted points/lines: 1 of 1<br>
v.to.rast complete.<br>
Reading maps...<br>
Finding nodes...<br>
Calculate downstream distance<br>
LuLitarroni61_dem_dist2out Done<br>
##################################<br>
Tot. cells 13.0<br>
===========================<br>
Ipsometric | quantiles<br>
===========================<br>
10 | 0.025<br>
10 | 0.05<br>
10 | 0.1<br>
10 | 0.25<br>
10 | 0.5<br>
7 | 0.75<br>
9 | 0.7<br>
4 | 0.9<br>
4 | 0.975<br>
Done!<br>
##################################<br>
##################################<br>
Tot. cells 13.0<br>
Tot. area 1300.0<br>
Max distance 169.372871<br>
===========================<br>
Whidth Function | quantiles<br>
===========================<br>
-131 | 0.05<br>
-113 | 0.15<br>
-86 | 0.3<br>
-67 | 0.4<br>
-49 | 0.5<br>
-31 | 0.6<br>
-12 | 0.7<br>
15 | 0.85<br>
33 | 0.95<br>
Done!<br>
##################################<br>
Reading maps...<br>
Finding nodes...<br>
Calculate downstream distance<br>
LuLitarroni61_dem_hillslope_distance Done<br>
r.stats:<br>
WARNING: Cats for raster map <r_elevation_crop@aglientu> are either missing<br>
or have no explicit labels. Using nsteps=255.<br>
r.recode:<br>
r.recode complete. Raster map <LuLitarroni61_dem_height_average> created.<br>
Dev note: Adapted sites library used for vector points. (module should be<br>
updated to GRASS 6 vector library)<br>
WARNING: Site/vector attributes ignored.<br>
Complete ...<br>
Invalid map <LuLitarroni61_dem_hack><br>
Parse error<br>
ERROR: An error occurred while running r.mapcalc<br>
<br>
##################################<br>
<br>
An error occurred with the parameters calculation.<br>
Please note that outlet coordinates must belong to the river network.<br>
You might want to run r.stream.extract and choose coordinates matching with<br>
the extracted stream map.<br>
Please report to the authors any other problem not related with coordinates<br>
outlet.<br>
(Tue Jan 29 08:45:27 2013) Command finished (13 min 12 sec)<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://osgeo-org.1560.n6.nabble.com/error-accumulation-map-does-not-match-with-elevation-map-tp5029128p5030492.html" target="_blank">http://osgeo-org.1560.n6.nabble.com/error-accumulation-map-does-not-match-with-elevation-map-tp5029128p5030492.html</a><br>
Sent from the Grass - Users mailing list archive at Nabble.com.<br>
_______________________________________________<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" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div><font color="#666666">Margherita DI LEO </font></div><div><font color="#666666">Postdoctoral Researcher</font></div><div><font color="#666666"><br>
</font></div><div><font color="#666666">European Commission - DG JRC </font></div><div><font color="#666666">Institute for Environment and Sustainability (IES). Unit H03 – FRC</font></div><div><font color="#666666">Via Fermi, 2749</font></div>
<div><font color="#666666">I-21027 Ispra (VA) - Italy - TP 261</font></div><div><font color="#666666"> </font></div><div><font color="#666666">Tel. +39 0332 78 3600 </font></div><div><font color="#666666"><a href="mailto:margherita.di-leo@jrc.ec.europa.eu" target="_blank">margherita.di-leo@jrc.ec.europa.eu</a></font></div>
<div><font color="#666666"><br></font></div><div><font color="#666666">Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.</font></div>
</div>
</div>