On 03-07-15 11:04, Moritz Lennert wrote:
On 02/07/15 19:26, Markus Neteler wrote:
On Jul 2, 2015 5:18 PM, "Moritz Lennert" <mlenn...@club.worldonline.be
<mailto:mlenn...@club.worldonline.be>> wrote:
>
> Hello,
>
> When I run a script that loops over a long series of point data sets
and then does a series of raster calculations based on these data sets,
I sometimes get the following error:
>
> ERROR: Error reading raster data for row 239 of <MASK>
>
> Can someone explain what this means and how to debug this ?
Just to note that this does not seem to be an isolated problem:
* http://lists.osgeo.org/pipermail/grass-dev/2014-September/070584.html
* http://lists.osgeo.org/pipermail/grass-dev/2015-May/074937.html
Maybe some off_t issue?
Please post the region settings and related info...
> g.version -reb
GRASS 7.1.svn (2015)
./configure --prefix=/usr/lib --sysconfdir=/etc
--sharedstatedir=/var --enable-socket --enable-shared
--enable-largefile --with-postgres --with-mysql --with-pthread
--with-cxx --with-x --with-gdal --with-freetype --with-motif
--with-readline --with-nls --with-odbc --with-sqlite
--with-freetype-includes=/usr/include/freetype2
--with-tcltk-includes=/usr/include/tcl
--with-postgres-includes=/usr/include/postgresql
--with-mysql-includes=/usr/include/mysql
--with-proj-share=/usr/share/proj --with-python=/usr/bin/python-config
--with-cairo --with-geos --with-blas --with-lapack
--with-liblas=/usr/bin/liblas-config
libgis Revision: 64732
libgis Date: 2015-02-25 01:54:05 +0100 (mer 25 fév 2015)
PROJ.4: 4.9.1
GDAL/OGR: 1.10.1
GEOS: 3.4.2
SQLite: 3.8.10.2
> g.region -p
projection: 99 (Lambert Azimuthal Equal Area)
zone: 0
datum: etrs89
ellipsoid: grs80
north: 3177000
south: 2937000
west: 3790000
east: 4071000
nsres: 1000
ewres: 1000
rows: 240
cols: 281
cells: 67440
> r.info -g MASKnorth=3177000
south=2937000
east=4071000
west=3790000
nsres=1000
ewres=1000
rows=240
cols=281
cells=67440
datatype=CELL
ncats=1
Here's the error:
ERREUR: Error reading raster data for row 239 of <MASK>
Traceback (most recent call last):
File "../../calulate_huff.py", line 346, in <module>
quiet=True)
File
"/data/home/mlennert/SRC/GRASS/grass_trunk/dist.x86_64-unknown-linux-gnu/etc/python/grass/script/core.py",
line 376, in run_command
return handle_errors(returncode, returncode, args, kwargs)
File
"/data/home/mlennert/SRC/GRASS/grass_trunk/dist.x86_64-unknown-linux-gnu/etc/python/grass/script/core.py",
line 312, in handle_errors
returncode=returncode)
grass.exceptions.CalledModuleError: Module run None ['r.mapcalc',
'--o', '--q', 'expression=temp_prob = float(firm_rate_364596) /
float(sum_rates)'] ended with error
Process ended with non-zero return code 1. See errors in the (error)
output.
And here the info for the two maps in the r.mapcalc call:
> r.info -g firm_rate_364596
north=3177000
south=2937000
east=4071000
west=3790000
nsres=1000
ewres=1000
rows=240
cols=281
cells=67440
datatype=DCELL
ncats=0
> r.info -g sum_rates
north=3177000
south=2937000
east=4071000
west=3790000
nsres=1000
ewres=1000
rows=240
cols=281
cells=67440
datatype=DCELL
ncats=0
But running the same script on the same maps, the error always appears
at a different stage (other firm_rate_* map).
Moritz
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev