Hi, everybody!
i created a imagem by kriging using geoR package. I imported points
from GRASS("zn", after converted to geodata "zn_geo"), the
border
"zn_border" and a raster mask. Then i interpolated the points by
kriging and created a raster image. Now, i need export this image back
to GRASS to use it in the module r.mapcalc. I can't do it. I tried use
writeRAST6, but it returns a message: "tentativa de obter um slot
"data" de um objeto cuja classe ("kriging") n?o ? um objeto
S4" or
"attempt to obtain a slot from an object whose class ("kriging")
is
not a S4 object". what does this mean?
Help me, please!
Thyb?rio Luna - Brazil
On Tue, Dec 9, 2008 at 6:03 PM, Thyb?rio Luna Freire <thyberio at gmail.com> wrote:> Hi, everybody! > i created a imagem by kriging using geoR package. I imported points > from GRASS("zn", after converted to geodata "zn_geo"), the border > "zn_border" and a raster mask. Then i interpolated the points by > kriging and created a raster image. Now, i need export this image back > to GRASS to use it in the module r.mapcalc. I can't do it. I tried use > writeRAST6, but it returns a message: "tentativa de obter um slot > "data" de um objeto cuja classe ("kriging") n?o ? um objeto S4" or > "attempt to obtain a slot from an object whose class ("kriging") is > not a S4 object". what does this mean? > Help me, please! > > Thyb?rio Luna - Brazil >Hi, Don't know much about geoR, but you probably need to coerce your kriging results into an sp-class object. Some ideas here: http://casoilresource.lawr.ucdavis.edu/drupal/node/438 specifically: library(spGRASS6) ?SpatialGridDataFrame Cheers, Dylan
Thyb?rio Luna Freire <thyberio <at> gmail.com> writes:> > Hi, everybody! > i created a imagem by kriging using geoR package. I imported points > from GRASS("zn", after converted to geodata "zn_geo"), the border > "zn_border" and a raster mask. Then i interpolated the points by > kriging and created a raster image. Now, i need export this image back > to GRASS to use it in the module r.mapcalc. I can't do it. I tried use > writeRAST6, but it returns a message: "tentativa de obter um slot > "data" de um objeto cuja classe ("kriging") n?o ? um objeto S4" or > "attempt to obtain a slot from an object whose class ("kriging") is > not a S4 object". what does this mean?It means what it says, as the help page to krige.conv() says. To make progress, you need to create a SpatialGridDataFrame matching the GRASS location (region, resolution). Use this to generate the prediction locations. The output object "predict" component contains the predictions. Add it to the SpatialGridDataFrame (or SpatialPixelsDataFrame), and proceed as writeRAST6() suggests in the helpfile: kc <- krige.conv(coords=coordinates(SPz2), data=SPz2$day16, locations=coordinates(SPDF_utm), krige=krige.control(obj.model=fvg)) SPDF_utm$geoR_pred <- kc$predict writeRAST6(SPDF_utm["geoR_pred"], "geoR_pred") or similar (untried). Note that either of the statgrass or R-sig-geo lists are more relevant than this one. Roger Bivand> Help me, please! > > Thyb?rio Luna - Brazil >