SpatialGridDataFrame
Hi, The problem is that your coordinates are not exactly on a grid. The very low suggested threshold indicates that this is only a marginal deviation, probably due to rounding of the coordinates. Try using the points2grid function directly, see ?points2grid. Specifically, take a look at the tolerance argument of points2grid. Paul
On 09/23/2011 08:01 AM, Eugenio Trumpy wrote:
Hi to all, i'm new of the list, and i have a poor experience unfortunatly in geostatistic. I'm working with grass and R in a linux debian pc. I'm using grass 6.4.1 with R 2.13.1. I have loaded spgrass6, gstat, lattice libraries to create some interpolated map by kriging method. I created the map with universal kriging, and now i'm tryng to write this map in grass with:
writeRAST6(UK_avampaese,"UK_avampaese",zcol="var1.pred")
but i have this error message: Errore: gridded(dataset) is not TRUE I cheked the dataset typology, UK_avampaese is a SpatialPointDataFrame, and writeRAST6 needs SpatialGridDataFrame. Which is the procedure to convert the a SpatialPointDataFrame to SpatialGridDataFrame? Is possible to make this changeover? I also tried:
gridded(UK_avampaese)=T
but i received this message: suggested tolerance minimum: 4.61319e-08 Errore in points2grid(points, tolerance, round) : dimension 2 : coordinate intervals are not constant Inoltre: Warning message: In points2grid(points, tolerance, round) : grid has empty column/rows in dimension 1 I would like to specify that the prediction grid used as newdata, is a SpatialPointDataFrame coming from grass point vector data, who have an irregular boundary and 5 km spaced. Any hits are very appreced, bast regards Eugenio [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Paul Hiemstra, Ph.D. Global Climate Division Royal Netherlands Meteorological Institute (KNMI) Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39 P.O. Box 201 | 3730 AE | De Bilt tel: +31 30 2206 494 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770