[GRASS-dev] r.volume - error: reading raster at row x outside of region
Helmut Kudrnovsky
hellik at web.de
Sun Apr 24 13:06:28 PDT 2016
Helmut Kudrnovsky wrote
> Hi devs,
>
> I have 2 rasters:
>
> r.info map=sample_tallard_eu_dem_westernalps_slope
>
> +----------------------------------------------------------------------------+
> | Map: sample_tallard_eu_dem_western Date: Sun Apr 24 20:36:43 2016
> |
> | Mapset: westernalpsriver Login of Creator: bugs
> |
> | Location: eudem_epsg3035
> |
> | DataBase: /home/bugs/grassdata
> |
> | Title: slope in degrees ( sample_tallard_eu_dem_westernalps_slope )
> |
> | Timestamp: none
> |
>
> |----------------------------------------------------------------------------|
> |
> |
> | Type of Map: raster Number of Categories: 77
> |
> | Data Type: FCELL
> |
> | Rows: 4713
> |
> | Columns: 4808
> |
> | Total Cells: 22660104
> |
> | Projection: Lambert Azimuthal Equal Area
> |
> | N: 2460525 S: 2342700 Res: 25
> |
> | E: 4098800 W: 3978600 Res: 25
> |
> | Range of data: min = 0 max = 76.56607
> |
> |
> |
> | Data Source:
> |
> | raster elevation file r_elevation_crop
> |
> |
> |
> |
> |
> | Data Description:
> |
> | generated by r.slope.aspect
> |
> |
> |
> | Comments:
> |
> | slope map elev = r_elevation_crop
> |
> | zfactor = 1.00 format = degrees
> |
> | min_slope = 0.000000
> |
> |
> |
> | r.slope.aspect elevation="r_elevation_crop"
> slope="sample_tallard_eu\ |
> | _dem_westernalps_slope"
> aspect="sample_tallard_eu_dem_westernalps_as\ |
> | pect" format="degrees" precision="FCELL" zscale=1.0 min_slope=0.0
> |
> |
> |
>
> +----------------------------------------------------------------------------+
>
>
> r.info map=sample_tallard_eu_dem_westernalps_basin
>
> +----------------------------------------------------------------------------+
> | Map: sample_tallard_eu_dem_western Date: Sun Apr 24 20:36:03 2016
> |
> | Mapset: westernalpsriver Login of Creator: bugs
> |
> | Location: eudem_epsg3035
> |
> | DataBase: /home/bugs/grassdata
> |
> | Title: ( sample_tallard_eu_dem_westernalps_basin )
> |
> | Timestamp: none
> |
>
> |----------------------------------------------------------------------------|
> |
> |
> | Type of Map: raster Number of Categories: 0
> |
> | Data Type: CELL
> |
> | Rows: 4713
> |
> | Columns: 4808
> |
> | Total Cells: 22660104
> |
> | Projection: Lambert Azimuthal Equal Area
> |
> | N: 2460525 S: 2342700 Res: 25
> |
> | E: 4098800 W: 3978600 Res: 25
> |
> | Range of data: min = 1 max = 1
> |
> |
> |
> | Data Description:
> |
> | generated by r.stream.basins
> |
> |
> |
> | Comments:
> |
> | r.stream.basins
> direction="sample_tallard_eu_dem_westernalps_drainag\ |
> | e_e" points="sample_tallard_eu_dem_westernalps_outlet_snap"
> memory=3\ |
> | 00 basins="sample_tallard_eu_dem_westernalps_basin"
> |
> |
> |
>
> +----------------------------------------------------------------------------+
>
> when I do r.volume with these rasters, I get:
>
> r.volume --verbose input=sample_tallard_eu_dem_westernalps_slope
> clump=sample_tallard_eu_dem_westernalps_basin
> ERROR: Die Leseanfrage der Rasterkarte
> <sample_tallard_eu_dem_westernalps_basin@> für Zeile -82 ist
> außerhalb der Region.
> ERROR: Reading of raster map
> <sample_tallard_eu_dem_westernalps_basin@> for row -82 is outside
> the region.
>
> D3/3: row=4712 col=4805: zero or negs ignored
> D3/3: row=4712 col=4806: zero or negs ignored
> D3/3: row=4712 col=4807: zero or negs ignored
> 100%
> FEHLER: Die Leseanfrage der Rasterkarte
> <sample_tallard_eu_dem_westernalps_basin@> für Zeile -82 ist
> außerhalb der Region.
>
> the computational region:
>
> g.region -p -a raster=eu_dem_westernalps at westernalpsriver
> vector=v_region_durance_buffer5000 at westernalpsriver
> align=eu_dem_westernalps at westernalpsriver
> projection: 99 (Lambert Azimuthal Equal Area)
> zone: 0
> datum: etrs89
> ellipsoid: grs80
> north: 2460525
> south: 2342700
> west: 3978600
> east: 4098800
> nsres: 25
> ewres: 25
> rows: 4713
> cols: 4808
> cells: 22660104
>
> any idea what causes such an error?
>
> r.volume, tested with a lot of other rasters of different extents, worked
> so far ...
Maybe some memory issue?
-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/r-volume-error-reading-raster-at-row-x-outside-of-region-tp5262679p5262683.html
Sent from the Grass - Dev mailing list archive at Nabble.com.
More information about the grass-dev
mailing list