http://r-sig-geo.2731867.n2.nabble.com/file/n7580281/Moran_XY.csv Moran_XY.csv List readers, I am trying to perform an ordinary kriging with the above data using the 'krige.conv' function in the geoR package. I am however, encountering problems with successfully completing this. I have provided the data and code I have been running and would welcome some suggestions. data<-read.csv("Moran_XY.csv") sb<-as.geodata(data) # Defining a prediction grid grid<-expand.grid(seq(300000,360000,l=60),seq(3613000,3650000,l=60)) # predicting by ordinary kriging vario <- variog(sb) vf <- variofit(vario) kc <- krige.conv(sb, loc=grid, krige=krige.control(obj.m=vf)) #generating image image(kc) Thanks in advance, Dan. -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/ordinary-kriging-using-krige-conv-function-tp7580281.html Sent from the R-sig-geo mailing list archive at Nabble.com.
ordinary kriging using 'krige.conv' function
2 messages · djrwicks, Giuseppe Calamita
6 days later
I never worked with geoR package for interpolation pourposes. Have you already checked the gstat package ? Further, do you have any idea of which model would be the best for fitting the data you're working on (exponential, spherical, gaussian, etc ?). Give me more information so I can help you better. Cheers Giuseppe Calamita -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/ordinary-kriging-using-krige-conv-function-tp7580281p7580349.html Sent from the R-sig-geo mailing list archive at Nabble.com.