Hi all,
I would like to consider measurement uncertainty for kriging. Have you done this before?
Is there any option in the kriging function to add this uncertainty?
Thanks a lot,
Best,
Alicia
library(RGeostats)
# Variogram model (exponential)
model = model.create(2,range=20)
# Targets grid
grid = db.create(nx=c(100,100))
# Neighborhood (unique)
neigh = neigh.create(type=0)
# Input dataset
dat = db.create(x1=c(20,30),x2=c(50,50),z1=c(-1,1))
# Simple kriging without measurement error
res1 = kriging(dat,grid,model,neigh,uc="",mean=0)
plot(res1,title="Kriging without VEM",pos.legend=1)
plot(dat,add=T,pos.legend=7)
# Add measurement error on the second point only
datm = db.add(dat,error=c(0,100),loctype="v")
# Simple kriging with measurement error
res2 = kriging(datm,grid,model,neigh,uc="")
plot(res2,title="Kriging with VEM",pos.legend=1)
plot(datm,add=T,name.prop="error",pos.legend=7)
Users browsing this forum: No registered users and 2 guests