[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