On Thu, 16 Oct 2008, Jarek Jasiewicz wrote:
> upps...
> I've just updated and....
to what? Please include sessionInfo() output and exact GRASS version, and
exact GDAL version (what you see when you load rgdal) and whether this is
the old (deprecated) or the new plugin. There are lots of things that
could go wrong, and this looks like a plugin built against a different
version of GRASS.
This problem (r.out.bin always using 16-bit integers for CELL rasters) was
discussed in a thread on R-sig-geo recently:
https://stat.ethz.ch/pipermail/r-sig-geo/2008-September/004184.htmland resolved off-list by implementing the plugin. The CVS rgdal on
sourceforge and spgrass6 on sourceforge also have a useGDAL= argument to
use a GTiff for the intermediate file, so r.out.gdal from GRASS >= 6.3.
See if gdalinfo on your GRASS raster also crashes - that would suggest a
plugin that is not synchronised with your GRASS and/or GDAL.
Hope this helps,
Roger
>
>> r=readRAST6("sca_mroga",plugin=TRUE)
>
> *** caught segfault ***
> address (nil), cause 'unknown'
>
> Traceback:
> 1: .Call("RGDAL_OpenDataset", as.character(filename), TRUE, PACKAGE =
> "rgdal")
> 2: .local(.Object, ...)
> 3: initialize(value, ...)
> 4: initialize(value, ...)
> 5: new("GDALReadOnlyDataset", filename)
> 6: GDAL.open(fname)
> 7: readGDAL(fname, silent = ignore.stderr)
> 8: readRAST6("sca_mroga", plugin = TRUE)
>
> Possible actions:
> 1: abort (with core dump, if enabled)
> 2: normal R exit
> 3: exit R without saving workspace
> 4: exit R saving workspace
>
> and I have very big problem....
> ...
> Markus Neteler pisze:
>> On Thu, Oct 16, 2008 at 9:02 PM, Jarek Jasiewicz <
[hidden email]>
>> wrote:
>>
>> > Hi:
>> >
>> > data from Grass: (integer)
>> >
>> > Range of data: min = 30 max = 31589370
>> > (there are more or less 14000 categories ranges as above)
>> >
>> > imported to R:
>> >
>> > summary(r@data$sca_mroga)
>> > Min. 1st Qu. Median Mean 3rd Qu. Max.
>> > -32770.0 60.0 150.0 532.1 330.0 32760.0
>> >
>> > I think the problems lies in 16bit proxy .bin file during export, and we
>> > have overflow
>> >
>>
>> Jarek,
>>
>> if your GRASS-R interface isn't too old, please try
>> to add plugin=TRUE as parameter.
>> It should enable the new method of reading in data.
>>
>> Please report if it (not) works.
>>
>> Best,
>> Markus
>>
> _______________________________________________
> grass-stats mailing list
>
[hidden email]
>
http://lists.osgeo.org/mailman/listinfo/grass-stats>
>
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail:
[hidden email]
_______________________________________________
grass-stats mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-statsRoger Bivand
Economic Geography Section
Department of Economics
Norwegian School of Economics and Business Administration
Helleveien 30
N-5045 Bergen, Norway
