[R-sig-Geo] ordinary kriging using 'krige.conv' function

djrwicks djrwicks at gmail.com
Thu Jun 21 13:57:28 CEST 2012


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.



More information about the R-sig-Geo mailing list