[R] 3-D interpretation

Sarah Goslee sarah.goslee at gmail.com
Tue Dec 10 16:41:21 CET 2013


Hi Shane,

On Tue, Dec 10, 2013 at 10:35 AM, Shane Carey <careyshan at gmail.com> wrote:
> 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.

This is a plain-text email list, so no highlighting in red is available.

More importantly, you should take a look at:

http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example

and provide a simple example with data that demonstrates your problem.

Sarah

>
> 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)
>
>
>

-- 
Sarah Goslee
http://www.functionaldiversity.org



More information about the R-help mailing list