[GRASS-user] error: accumulation map does not match with elevation map!

Tue Jan 29 00:26:01 PST 2013

Hi,

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.

HTH
madi

On Tue, Jan 29, 2013 at 8:45 AM, Giuliano Urgeghe <giuliano.u at gmail.com>wrote:

> I think I have solved the error accumulation map does not match with the
> elevation map installing 6.4.3
> but now I get this error:
> Tue Jan 29 08:32:15 2013)
> r.basin map=dem at PERMANENT prefix=LuLitarroni61 easting=1509351.753173
> northing=4553749.816433 threshold=100
> SECTION 1a (of 4): Initiating Memory.
> SECTION 1b (of 4): Determining Offmap Flow.
> SECTION 2: A * Search.
> SECTION 3: Accumulating Surface Flow with SFD.
> SECTION 4: Closing Maps.
> WARNING: Writing out only positive flow accumulation values.
> WARNING: Cells with a likely underestimate for flow accumulation can no
> longer be identified.
> Load input maps and get start points...
> Set edge points...
> A* Search...
> Extract streams...
> Thin stream segments...
> Write raster maps
> Reading maps...
> Calculate basins using coordinates...
> LuLitarroni61_dem_basin Done!
> Delineation of basin done
> Rename raster <tmp> to <LuLitarroni61_dem_accumulation>
> Rename raster <tmp> to <LuLitarroni61_dem_drainage>
> Rename raster <tmp> to <LuLitarroni61_dem_drainage_e>
> Rename raster <tmp> to <LuLitarroni61_dem_stream_e>
> File LuLitarroni61_dem_stream_e -- 2426 rows X 3003 columns
> ERROR: (null): Unable to find bounding box for lines
> ERROR: Raster map <LuLitarroni61_dem_stream_e_thin> not found
> Aspect raster map <LuLitarroni61_dem_aspect> complete
> Slope raster map <LuLitarroni61_dem_slope> complete
> Extracting areas...
> Building topology for vector map <LuLitarroni61_dem_basin>...
> Registering primitives...
> 37 primitives registered
> 189 vertices registered
> Building areas...
> 13 areas built
> 1 isles built
> Attaching islands...
> Attaching centroids...
> Number of nodes: 25
> Number of primitives: 37
> Number of points: 0
> Number of lines: 0
> Number of boundaries: 24
> Number of centroids: 13
> Number of areas: 13
> Number of isles: 1
> r.to.vect complete.
> Reading areas...
> Updating database...
> 1 categories read from vector map (layer 1)
> 1 records selected from table (layer 1)
> 1 categories read from vector map exist in selection from table
> 1 records updated/inserted (layer 1)
> Reading areas...
> 1 categories read from vector map (layer 1)
> 0 records updated/inserted (layer 1)
> Reading areas...
> Updating database...
> 1 categories read from vector map (layer 1)
> 1 records selected from table (layer 1)
> 1 categories read from vector map exist in selection from table
> 1 records updated/inserted (layer 1)
> Reading areas...
> 1 categories read from vector map (layer 1)
> 0 records updated/inserted (layer 1)
> Creating LuLitarroni61_dem_hack
> Reading maps...
> ERROR: No streams found
> Scanning input for column types...
> Maximum input row length: 33
> Maximum number of columns: 3
> Minimum number of columns: 3
> Column: 1  type: double
> Column: 2  type: double
> Column: 3  type: integer
> Importing points...
> Populating table...
> Building topology for vector map <LuLitarroni61_dem_outlet>...
> Registering primitives...
> 1 primitives registered
> 1 vertices registered
> Building areas...
> 0 areas built
> 0 isles built
> Attaching islands...
> Attaching centroids...
> Number of nodes: 1
> Number of primitives: 1
> Number of points: 1
> Number of lines: 0
> Number of boundaries: 0
> Number of centroids: 0
> Number of areas: 0
> Number of isles: 0
> v.in.ascii complete.
> Loading data...
> Reading features...
> Writing raster map...
> Converted points/lines: 1 of 1
> v.to.rast complete.
> Reading maps...
> Finding nodes...
> Calculate downstream distance
> LuLitarroni61_dem_dist2out Done
> ##################################
> Tot. cells 13.0
> ===========================
> Ipsometric | quantiles
> ===========================
> 10 | 0.025
> 10 | 0.05
> 10 | 0.1
> 10 | 0.25
> 10 | 0.5
> 7 | 0.75
> 9 | 0.7
> 4 | 0.9
> 4 | 0.975
> Done!
> ##################################
> ##################################
> Tot. cells 13.0
> Tot. area 1300.0
> Max distance 169.372871
> ===========================
> Whidth Function | quantiles
> ===========================
> -131 | 0.05
> -113 | 0.15
> -86 | 0.3
> -67 | 0.4
> -49 | 0.5
> -31 | 0.6
> -12 | 0.7
> 15 | 0.85
> 33 | 0.95
> Done!
> ##################################
> Reading maps...
> Finding nodes...
> Calculate downstream distance
> LuLitarroni61_dem_hillslope_distance Done
> r.stats:
> WARNING: Cats for raster map <r_elevation_crop at aglientu> are either
> missing
> or have no explicit labels. Using nsteps=255.
> r.recode:
> r.recode complete. Raster map <LuLitarroni61_dem_height_average> created.
> Dev note: Adapted sites library used for vector points. (module should be
> updated to GRASS 6 vector library)
> WARNING: Site/vector attributes ignored.
> Complete ...
> Invalid map <LuLitarroni61_dem_hack>
> Parse error
> ERROR: An error occurred while running r.mapcalc
>
> ##################################
>
> An error occurred with the parameters calculation.
> Please note that outlet coordinates must belong to the river network.
> You might want to run r.stream.extract and choose coordinates matching with
> the extracted stream map.
> Please report to the authors any other problem not related with coordinates
> outlet.
> (Tue Jan 29 08:45:27 2013) Command finished (13 min 12 sec)
>
>
>
> --
> View this message in context:
> http://osgeo-org.1560.n6.nabble.com/error-accumulation-map-does-not-match-with-elevation-map-tp5029128p5030492.html
> Sent from the Grass - Users mailing list archive at Nabble.com.
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>

-- 
Margherita DI LEO
Postdoctoral Researcher

European Commission - DG JRC
Institute for Environment and Sustainability (IES). Unit H03 – FRC
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo at jrc.ec.europa.eu

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20130129/497e6952/attachment-0001.html>