[R] 3-D interpretation

Sarah Goslee sarah.goslee at gmail.com
Tue Dec 10 16:48:37 CET 2013


Hi,

On Tue, Dec 10, 2013 at 10:42 AM, Shane Carey <careyshan at gmail.com> wrote:
> Hi, I have since solved that problem,

It would be useful to post the solution, for the archives, so that
others with the same problem can benefit from your insight.


> but now R keeps crashing. I have over 20,000 points. Can R handle that
> amount?

R can. Your computer may not be able to.

But again, we don't have enough information to answer you, since that
depends on OS, available RAM, etc.

Sarah

>
> Cheers
>
>
> On Tue, Dec 10, 2013 at 3:35 PM, 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.
>>
>> 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