Hi,

im trying to create a 3-D interpretation of a geological fault using
the akima package. But my code fails below highlighted in red:
Does anyone have any ideas why this is.

Thanks

dat<-read.delim("D:\\fault.txt",header=T,sep=",")
library(rgl)
# data
rgl.spheres(dat$X,dat$Z , dat$Y,1,color="red")
rgl.bbox()
# bivariate linear interpolation
# interp:
akima.li <- interp(dat$X, dat$Y, dat$Z,
                   xo=seq(min(dat$X), max(dat$X), length = 100),
                   yo=seq(min(dat$Y), max(dat$Y), length = 100))

This is my error:
Error in interp.old(x, y, z, xo = xo, yo = yo, ncp = 0, extrap = extrap,  :
duplicate data points: need to set 'duplicate = ..'

# interp surface:
rgl.surface(akima.li$X,akima.li$Y,akima.li$Z,color="green",alpha=c(0.5))
# interpp:
akima.p <- interpp(dat$X, dat$Y, dat$Z,
                   runif(200,min(dat$X),max(dat$X)),
                   runif(200,min(dat$Y),max(dat$Y)))
# interpp points:
rgl.points(akima.p$X,akima.p$Z , akima.p$Y,size=4,color="yellow")
# bivariate spline interpolation
# data
rgl.spheres(dat$X,dat$Z ,dat$Y,0.5,color="red")
rgl.bbox()
# bivariate cubic spline interpolation
# interp:
akima.si <- interp(dat$X, dat$Y, dat$Z,
                   xo=seq(min(dat$X), max(dat$X), length = 100),
                   yo=seq(min(dat$Y), max(dat$Y), length = 100),
                   linear = FALSE, extrap = TRUE)






-- 
Shane

	[[alternative HTML version deleted]]

