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}} _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo