[R-sig-Geo] SpatialGridDataFrame

[Paul Hiemstra]

 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