[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
CFLAGS="-g -O2 -D_FILE_OFFSET_BITS=64
-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
CXX='g++-3.3' CXXFLAGS='-g -O2 -D_LARGE_FILE_SOURCE
-D_FILE_OFFSET_BITS=64' CFLAGS='-g -O2 -D_LARGE_FILE_SOURCE
-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.
-Andy
More information about the grass-user
mailing list