[GRASS-dev] what is the meaning of: "Error reading raster data for row 239 of <MASK>"

[Moritz Lennert]

On 02/07/15 19:26, Markus Neteler wrote:
> On Jul 2, 2015 5:18 PM, "Moritz Lennert" <mlennert at club.worldonline.be
> <mailto:mlennert at 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 ?
>
> 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