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

[Paulo van Breugel]

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" <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 ?
>>

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 at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev