[GRASSLIST:7976] error reading large rasters in 5.4/6.0

Andrew Danner adanner at cs.duke.edu
Thu Aug 18 18:01:05 EDT 2005

I'm having problems working with some large raster files in GRASS, even
though I have enabled large file support. 

I built GRASS from source

my 5.4 config

-D_LARGE_FILE_SOURCE" ../grass-5.4.0/configure
--with-tcltk-includes=/usr/include/tcl8.3 --without-postgres
--with-mysql --with-mysql-includes=/usr/include/mysql --with-cxx

my 6.0 config

-D_FILE_OFFSET_BITS=64' ./configure
--with-tcltk-includes=/usr/include/tcl/ --without-postgres --with-mysql
--with-mysql-includes=/usr/include/mysql/ --with-cxx

I get the following problems:

1) "r.info test" displays the wrong number of total cells. e.g., 

>  |   Type of Map:  cell 
>  |   Data Type:    DCELL         
>  |   Rows:         50590     
>  |   Columns:      71840 
>  |   Total Cells:  -660581696         
>  |        Projection: Lambert Conformal Conic (zone 0)      
>  |            N:     967200    S:     461300   Res:    10   
>  |            E:    2638600    W:    1920200   Res:    10
2) "d.rast test" complains. 
WARNING: error reading compressed map [test] in mapset [local], row 0

sometimes the image displays fine even after the warning, on another
data set I get a green bar at the top and if I change the region,
sometimes I don't get the message at all. For example, if I set my
region to be

g.region rast=test

I get the error using d.rast, but if I expand the northing by 10 feet
(one cell), I get no error. 

Even though the image tends to display fine and look correct, the
problem causes some programs to quit, e.g. r.mapcalc. 

r.mapcalc test2=test 
ERROR: error reading compressed map [Neuse10] in mapset [local], row 0

The problem seems to be on rasters whose fcell file is bigger than 4GB. 

Comments/ideas/fixes appreciated.


More information about the grass-user mailing list