On Sat, 10 Jan 2009, Markus Neteler wrote:
On Sat, Jan 3, 2009 at 11:09 AM, Roger Bivand <[email protected]> wrote:
The current settings of readRAST6() will use the plugin if available (but
try to check that the regions match) when plugin=NULL. If you want to force
the use of r.out.gdal, set plugin=FALSE and useGDAL=TRUE (version 0.5-16).
Please let me know if this helps.
I have installed R (R version 2.8.1 (2008-12-22)) + extensions from scratch.
I get now this message (GRASS book p. 358):
elev <- readRAST6("elev_state_500m", mapset="PERMANENT", ignore.stderr=TRUE)
cols rows origin.northing origin.easting
FALSE FALSE FALSE TRUE
Warning message:
In readRAST6("elev_state_500m", mapset = "PERMANENT", :
set plugin=FALSE - raster/current window mismatch
or plugin=TRUE to override; continuing with plugin=FALSE
To me it is not clear from the message that I should (better) set
useGDAL=TRUE
Ideally, the interface would work as simple as this
elev <- readRAST6("elev_state_500m")
which would
* find the map in the current mapset search path (internally, use g.findfile),
in my tests I was in a mapset different from PERMANENT which leads to
elev <- readRAST6("elev_state_500m", ignore.stderr=TRUE)
Error in .local(.Object, ...) :
GDAL Error 4:
`/home/neteler/grassdata/nc_spm_07/user1/cellhd/elev_state_500m' does
not exist in the file system,
and is not recognised as a supported dataset name.
* use the current region (internally, use r.in.gdal and not the plugin
if GRASS >= 6.4)
Markus,
This is effectively what I've tried to now. However, the user output is:
library(spgrass6)
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.0, released 2008/12/04
GDAL_DATA: /home/rsb/lib/r_libs/rgdal/gdal
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
PROJ_LIB: /home/rsb/lib/r_libs/rgdal/proj
GRASS GIS interface loaded with GRASS version: 6.4.0RC1
and location: nc_spm_08
elev <- readRAST6("elev_state_500m")
cols rows origin.northing origin.easting
FALSE FALSE FALSE FALSE
ERROR 6: SetColorTable() only supported for Byte or UInt16 bands in TIFF
format.
WARNING: Input raster map contains cells with NULL-value (no-data). The
value 999 was used to represent no-data values in the input map.
You can specify nodata value by nodata parameter.
/home/rsb/topics/grassdata/nc_spm_08/rsb/.tmp/reclus2/elev_state_500m has
GDAL driver GTiff
and has 620 rows and 1630 columns
Warning message:
In readRAST6("elev_state_500m") :
set plugin=FALSE - raster/current window mismatch
or plugin=TRUE to override; continuing with plugin=FALSE
which I agree is unhelpful (and most of it cannot be avoided on Windows,
where stderr cannot be ignored). I can work on some of this, but the two
r.out.gdal messages are not on my side. Even if I set GRASS_VERBOSE to 0,
won't they still be displayed?
Roger
Successfully, I have used today:
elev <- readRAST6("elev_state_500m", mapset="PERMANENT",
ignore.stderr=TRUE, plugin=FALSE, useGDAL=TRUE)
which is rather long/complex.
Hope I am not asking too much - try to compensate with GRASS 6.4.x work :)
Markus
--
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: [email protected]
_______________________________________________
grass-stats mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-stats
