[R-sig-Geo] Problem when loading GRASS 6 raster with large integer values

[Rainer M Krug]

Hi

when I try to import a raster layer from grass with large integer
values, I get wrong values in R (see below). When I change the
datatype in GRASS to double, it works.

Thanks,

Rainer



GRASS 6.3.0 (grass):~/Documents/Projects/AlienSpread/R > r.stats -c
input=tmpDispLocal
 100%
0 39798
11 831
58 78
105 38
242 1652
889 816
2205 294
2310 56
3315 1171
4935 116
7875 9
8085 67
11550 6
13010 299
24675 2
30135 425
47702 1
48510 981
103635 1255
165375 8
177797 160
177870 6
379995 324
606375 8
662970 292
1416345 445
2260125 41
* 50661
[Raster MASK present]


BUT:


>     seeds <- readRAST6(
+                        "tmpDispLocal",
+                        ignore.stderr=!Debug
+                        )
Creating BIL support files...
Header File =
/home/rkrug/Documents/Projects/AlienSpread/R/../grass/simulation/.tmp/ecolmod/tmpDispLocal.hdr
World File =
/home/rkrug/Documents/Projects/AlienSpread/R/../grass/simulation/.tmp/ecolmod/tmpDispLocal.wld
Exporting raster as integer values (bytes=2)
Using the current region settings...
north=6247000.000000
south=6195000.000000
east=270200.000000
west=251000.000000
r=520
c=192

> table(seeds[[1]])

-31233 -27437 -25447 -18811 -18738 -17834 -17026 -13221      0     11     58
     8   1255    445    160      6      1    981    324  39798    831     78
   105    242    889   2205   2310   3315   4935   7610   7875   8085  11550
    38   1652    816    294     56   1171    116    292      9     67      6
 13010  16551  24675  30135  31901
   299      8      2    425     41
>


Additional info:

GRASS 6.3.0 (grass):~/Documents/Projects/AlienSpread/R > r.info map=tmpDispLocal
 +----------------------------------------------------------------------------+
 | Layer:    tmpDispLocal                   Date: Tue Sep 16 16:45:35 2008    |
 | Mapset:   simulation                     Login of Creator: rkrug           |
 | Location: grass                                                            |
 | DataBase: /home/rkrug/Documents/Projects/AlienSpread/R/..                  |
 | Title:     ( tmpDispLocal )                                                |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 2260125         |
 |   Data Type:    CELL                                                       |
 |   Rows:         520                                                        |
 |   Columns:      192                                                        |
 |   Total Cells:  99840                                                      |
 |        Projection: UTM (zone -34)                                          |
 |            N:    6247000    S:    6195000   Res:   100                     |
 |            E:     270200    W:     251000   Res:   100                     |
 |   Range of data:    min = 0  max = 2260125                                 |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.mapcalc                                                  |
 |                                                                            |
 |   Comments:                                                                |
 |    round(tmpDisp * 0.300000)                                               |
 |                                                                            |
 +----------------------------------------------------------------------------+

[Raster MASK present]
GRASS 6.3.0 (grass):~/Documents/Projects/AlienSpread/R >

-- 
Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation
Biology, UCT), Dipl. Phys. (Germany)

Centre of Excellence for Invasion Biology
Faculty of Science
Natural Sciences Building
Private Bag X1
University of Stellenbosch
Matieland 7602
South Africa