Hey,
I have had another idea since. Is it possible to join these points together
(by lines) and then created a polygon from them?
Thanks
On Wed, Dec 11, 2013 at 10:43 PM, Shane Carey <careys...@gmail.com> wrote:
> Hey,
>
> Thanks for this. I think I need some way of transposing the data onto a
> horizontal plane and then running the interpolation. Would this work maybe?
> I really dont know at this stage.
> Thanks for your help
>
>
> On Wednesday, December 11, 2013, Ben Bolker wrote:
>
>>
>> I got a little bit farther. The results are still very ugly, jagged
>> for linear interpolation and with a huge range of values for cubic
>> spline extrapolation. I still don't know about the crash (can't
>> reproduce it), but don't know how you expect to be able to extrapolate a
>> sensible 3-D surface from points measured along a single narrow 1-D
>> curve in the X-Y plane ... Can you propose *any* software tool/method
>> that can do a reasonable job of this?
>>
>> source("akimatmp2.R") ## BMB: contents from your data dump below,
>> ## assigned to "dat2"
>> dat <- dat2
>> nrow(dat) ## 1500
>>
>> ## BMB: this makes things a little prettier, but just cosmetic
>> dat <- transform(dat,X=X-min(X),Y=Y-min(Y))
>> library(rgl)
>> library(akima)
>> 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),
>> duplicate="mean")
>> # interp surface:
>>
>> ## BMB: this worked for me once I specified *lowercase* (x,y,z)
>> ## rather than *uppercase* (X,Y,Z) elements.
>>
>> with(akima.li,rgl.surface(x,y,z,color="green",alpha=c(0.5)))
>>
>> akima.ci0 <- 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,
>> duplicate="mean")
>> akima.ci0$z[!is.finite(akima.ci0$z)] <- NA
>> with(akima.ci0,image(x,y,z))
>> range(na.omit(c(akima.ci0$z))) ## c() to treat as vector
>> ## -250000 to 101000
>> range(dat$Z) ## 3500 to 4400
>> summary(na.omit(c(akima.ci0$z)))
>> points(dat$X,dat$Y)
>>
>> rgl.close()
>> rgl.spheres(dat$X,dat$Z , dat$Y,1,color="red")
>> rgl.bbox()
>> with(akima.ci0,rgl.surface(x,y,z,color="green",alpha=c(0.5)))
>>
>>
>> # 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,
>> duplicate="mean")
>> with(akima.si,rgl.surface(x,y,z,color="blue",alpha=c(0.5)))
>> with(akima.si,image(x,y,z))
>> filled.contour(akima.si$x,akima.si$y,akima.si$z)
>>
>>
>> On 13-12-11 11:44 AM, Shane Carey wrote:
>> > I want it as one block, if you know what I mean. Like the akima example.
>> >
>> > Thanks
>> >
>> >
>> > On Wed, Dec 11, 2013 at 4:41 PM, Shane Carey <careys...@gmail.com>
>> wrote:
>> >
>> >> Ok, here is the first 1500 points, but Im giving up hope on it to be
>> >> honest.
>> >> It seems to be going crazy, creating triangles everywhere. Do the
>> points
>> >> need to be evenly spaced?
>> >> Thanks for your help.
>> >>
>> >> structure(list(X = c(9816966.951, 9816963.08, 9816947.516, 9816939.51,
>> >> 9816924.005, 9816916.096, 9816901.984, 9816896.967, 9816892.928,
>> >> 9816890.743, 9816886.427, 9816884.006, 9816879.185, 9816876.468,
>> >> 9816871.027, 9816868.276, 9816863.404, 9816860.712, 9816855.409,
>> >> 9816852.487, 9816846.77, 9816843.635, 9816837.53, 9816834.196,
>> >> 9816827.733, 9816824.214, 9816817.419, 9816811.926, 9816797.531,
>> >> 9816789.588, 9816774.134, 9816765.716, 9816749.492, 9816740.717,
>> >> 9816723.946, 9816714.931, 9816697.836, 9816675.781, 9816606.59,
>> >> 9816578.983, 9816544.487, 9816531.179, 9816514.376, 9816505.891,
>> >> 9816489.52, 9816481.288, 9816465.377, 9816457.374, 9816441.836,
>> >> 9816434.026, 9816418.807, 9816411.161, 9816396.205, 9816385.052,
>> >> 9816355.719, 9816341.037, 9816312.075, 9816304.796, 9816304.784,
>> >> 9816916.962, 9816917.365, 9816918.982, 9816919.414, 9816919.519,
>> >> 9816919.545, 9816599.694, 9816615.198, 9816677.27, 9816697.31,
>> >> 9816714.671, 9816723.431, 9816740.464, 9816748.995, 9816765.474,
>> >> 9816773.661, 9816789.36, 9816797.089, 9816811.764, 9816817.212,
>> >> 9816824.113, 9816827.536, 9816834.1, 9816837.345, 9816843.545,
>> >> 9816846.596, 9816852.404, 9816855.248, 9816860.635, 9816863.257,
>> >> 9816868.198, 9816870.861, 9816876.39, 9816879.039, 9816883.937,
>> >> 9816886.296, 9816890.68, 9816892.809, 9816896.829, 9816901.555,
>> >> 9816915.878, 9816921.977, 9816931.382, 9816933.721, 9816619.581,
>> >> 9816631.77, 9816680.572, 9816697.262, 9816714.648, 9816723.385,
>> >> 9816740.441, 9816748.95, 9816765.452, 9816773.618, 9816789.339,
>> >> 9816797.049, 9816811.749, 9816817.193, 9816824.103, 9816827.518,
>> >> 9816834.092, 9816837.328, 9816843.537, 9816846.58, 9816852.396,
>> >> 9816855.233, 9816860.628, 9816863.244, 9816868.191, 9816870.846,
>> >> 9816876.383, 9816879.025, 9816883.93, 9816886.284, 9816890.674,
>> >> 9816892.798, 9816896.817, 9816901.516, 9816915.859, 9816921.464,
>> >> 9816929.009, 9816930.886, 9816639.468, 9816648.343, 9816683.875,
>> >> 9816697.215, 9816714.624, 9816723.338, 9816740.418, 9816748.905,
>> >> 9816765.43, 9816773.575, 9816789.318, 9816797.009, 9816811.734,
>> >> 9816817.175, 9816824.094, 9816827.5, 9816834.083, 9816837.311,
>> >> 9816843.528, 9816846.565, 9816852.389, 9816855.218, 9816860.621,
>> >> 9816863.23, 9816868.184, 9816870.831, 9816876.376, 9816879.012,
>> >> 9816883.924, 9816886.272, 9816890.669, 9816892.787, 9816896.805,
>> >> 9816901.477, 9816915.839, 9816920.952, 9816926.637, 9816928.05,
>> >> 9816659.354, 9816664.915, 9816687.177, 9816697.167, 9816714.6,
>> >> 9816723.291, 9816740.395, 9816748.859, 9816765.407, 9816773.532,
>> >> 9816789.298, 9816796.969, 9816811.719, 9816817.156, 9816824.085,
>> >> 9816827.482, 9816834.074, 9816837.294, 9816843.52, 9816846.549,
>> >> 9816852.381, 9816855.204, 9816860.614, 9816863.217, 9816868.176,
>> >> 9816870.816, 9816876.369, 9816878.999, 9816883.918, 9816886.26,
>> >> 9816890.663, 9816892.777, 9816896.792, 9816901.438, 9816915.819,
>> >> 9816920.439, 9816924.264, 9816925.215, 9816679.241, 9816681.487,
>> >> 9816690.48, 9816697.119, 9816714.577, 9816723.244, 9816740.372,
>> >> 9816748.814, 9816765.385, 9816773.489, 9816789.277, 9816796.929,
>> >> 9816811.705, 9816817.137, 9816824.076, 9816827.464, 9816834.065,
>> >> 9816837.277, 9816843.512, 9816846.533, 9816852.373, 9816855.189,
>> >> 9816860.607, 9816863.203, 9816868.169, 9816870.801, 9816876.362,
>> >> 9816878.985, 9816883.911, 9816886.248, 9816890.657, 9816892.766,
>> >> 9816896.78, 9816901.399, 9816915.799, 9816919.926, 9816921.892,
>> >> 9816922.38, 9816541.016, 9816548.667, 9816579.391, 9816605.021,
>> >> 9816675.293, 9816697.454, 9816714.742, 9816723.572, 9816740.533,
>> >> 9816749.131, 9816765.54, 9816773.79, 9816789.422, 9816797.21,
>> >> 9816811.808, 9816817.269, 9816824.14, 9816827.589, 9816834.126,
>> >> 9816837.395, 9816843.569, 9816846.643, 9816852.427, 9816855.292,
>>
>
>
> --
> Shane
>
--
Shane
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
