problem in krige function
Roger Bivand wrote:
On Thu, 5 Apr 2007, temiz wrote:
following the topic GRASS and R kriging interpolation from grass6 tutorial
I came here without problems
But here:
zz <- krige(z ~ 1, locations = el1 , newdata = mask_SG , model = efitted)
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, :
var1 : data item in gstat object and newdata have different
coordinate reference systems
If: proj4string(el1) proj4string(mask_SG) differ, the function will stop you predicting from data declared to be in one projection to newdata points in a different projection. You can set them by: proj4string(mask_SG) <- CRS(proj4string(el1)) (untried). If this does not work, post the commands you used to create the objects. Roger
how can I solve this problem ?
regards
Thank you I regarded some modifications in http://grass.gdf-hannover.de/wiki/GRASS_6_Tutorial#GRASS_and_R_kriging_interpolation I encountered a new problem : > OK_pred <- krige(z ~ 1, el3, newdata=mask_SG, model=efitted) [using ordinary kriging] "memory.c", line 57: can't allocate memory in function m_get() Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : m_get Why does this problem arise ? I selected a great deal smaller area in R than I did in grass and carried out with spline tension without any problem. is the data too dense or grid size too small ? regards -- Ahmet Temiz
This message has been scanned for viruses and\ dangerous con...{{dropped}}