[R] 3-D interpretation

Ben Bolker bbolker at gmail.com
Wed Dec 11 19:35:00 CET 2013


    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 <careyshan at 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,
>> 9816860.656, 9816863.297, 9816868.219, 9816870.906, 9816876.411,
>> 9816879.078, 9816883.956, 9816886.331, 9816890.697, 9816892.842,
>> 9816896.867, 9816901.672, 9816915.938, 9816923.514, 9816938.5,
>> 9816942.226, 9816560.486, 9816564.893, 9816582.587, 9816604.825,
>> 9816675.232, 9816697.406, 9816714.719, 9816723.525, 9816740.51,
>> 9816749.085, 9816765.518, 9816773.747, 9816789.401, 9816797.17,
>> 9816811.793, 9816817.25, 9816824.131, 9816827.572, 9816834.118,
>> 9816837.378, 9816843.561, 9816846.628, 9816852.419, 9816855.277,
>> 9816860.649, 9816863.284, 9816868.212, 9816870.891, 9816876.404,
>> 9816879.065, 9816883.949, 9816886.32, 9816890.692, 9816892.831,
>> 9816896.854, 9816901.633, 9816915.918, 9816923.002, 9816936.127,
>> 9816939.391, 9816579.957, 9816581.118, 9816585.783, 9816604.629,
>> 9816675.171, 9816697.358, 9816714.695, 9816723.478, 9816740.487,
>> 9816749.04, 9816765.496, 9816773.704, 9816789.381, 9816797.13,
>> 9816811.778, 9816817.231, 9816824.122, 9816827.554, 9816834.109,
>> 9816837.361, 9816843.553, 9816846.612, 9816852.411, 9816855.262,
>> 9816860.642, 9816863.27, 9816868.205, 9816870.876, 9816876.397,
>> 9816879.052, 9816883.943, 9816886.308, 9816890.686, 9816892.82,
>> 9816896.842, 9816901.594, 9816915.898, 9816922.489, 9816933.755,
>> 9816936.556, 9816521.879, 9816524.089, 9816532.945, 9816543.835,
>> 9816578.476, 9816605.217, 9816675.354, 9816697.501, 9816714.766,
>> 9816723.618, 9816740.556, 9816749.176, 9816765.562, 9816773.833,
>> 9816789.443, 9816797.25, 9816811.823, 9816817.287, 9816824.149,
>> 9816827.607, 9816834.135, 9816837.412, 9816843.577, 9816846.659,
>> 9816852.434, 9816855.306, 9816860.663, 9816863.311, 9816868.226,
>> 9816870.922, 9816876.419, 9816879.092, 9816883.962, 9816886.343,
>> 9816890.703, 9816892.852, 9816896.879, 9816901.711, 9816915.958,
>> 9816923.725, 9816939.367, 9816943.604, 9816944.908, 9816945.232,
>> 9816503.043, 9816504.183, 9816508.745, 9816514.107, 9816530.972,
>> 9816543.928, 9816578.549, 9816605.413, 9816675.415, 9816697.549,
>> 9816714.79, 9816723.665, 9816740.579, 9816749.221, 9816765.584,
>> 9816773.876, 9816789.464, 9816797.29, 9816811.837, 9816817.306,
>> 9816824.158, 9816827.625, 9816834.144, 9816837.429, 9816843.586,
>> 9816846.675, 9816852.442, 9816855.321, 9816860.67, 9816863.324,
>> 9816868.233, 9816870.937, 9816876.426, 9816879.105, 9816883.968,
>> 9816886.355, 9816890.709, 9816892.863, 9816896.892, 9816901.75,
>> 9816915.977, 9816923.765, 9816939.387, 9816944.163, 9816947.504,
>> 9816948.335, 9816483.505, 9816483.785, 9816484.906, 9816489.303,
>> 9816505.781, 9816514.152, 9816531.007, 9816544.021, 9816578.621,
>> 9816605.61, 9816675.476, 9816697.597, 9816714.813, 9816723.712,
>> 9816740.602, 9816749.266, 9816765.606, 9816773.919, 9816789.485,
>> 9816797.33, 9816811.852, 9816817.325, 9816824.168, 9816827.643,
>> 9816834.152, 9816837.446, 9816843.594, 9816846.691, 9816852.449,
>> 9816855.336, 9816860.677, 9816863.337, 9816868.24, 9816870.952,
>> 9816876.433, 9816879.118, 9816883.975, 9816886.367, 9816890.714,
>> 9816892.874, 9816896.904, 9816901.789, 9816915.997, 9816923.805,
>> 9816939.407, 9816944.722, 9816950.1, 9816951.437, 9816440.173,
>> 9816443.676, 9816457.692, 9816465.209, 9816481.202, 9816489.346,
>> 9816505.803, 9816514.197, 9816531.041, 9816544.114, 9816578.693,
>> 9816605.806, 9816675.537, 9816697.645, 9816714.837, 9816723.759,
>> 9816740.625, 9816749.311, 9816765.628, 9816773.962, 9816789.505,
>> 9816797.37, 9816811.867, 9816817.344, 9816824.177, 9816827.661,
>> 9816834.161, 9816837.463, 9816843.602, 9816846.707, 9816852.457,
>> 9816855.35, 9816860.684, 9816863.351, 9816868.247, 9816870.967,
>> 9816876.44, 9816879.132, 9816883.981, 9816886.379, 9816890.72,
>> 9816892.885, 9816896.917, 9816901.828, 9816916.017, 9816923.845,
>> 9816939.428, 9816945.281, 9816952.696, 9816954.54, 9816398.533,
>> 9816401.252, 9816412.135, 9816418.687, 9816433.966, 9816441.714,
>> 9816457.312, 9816465.251, 9816481.224, 9816489.39, 9816505.825,
>> 9816514.241, 9816531.076, 9816544.207, 9816578.766, 9816606.002,
>> 9816675.598, 9816697.693, 9816714.86, 9816723.806, 9816740.648,
>> 9816749.356, 9816765.65, 9816774.005, 9816789.526, 9816797.41,
>> 9816811.882, 9816817.363, 9816824.186, 9816827.679, 9816834.17,
>> 9816837.48, 9816843.61, 9816846.722, 9816852.465, 9816855.365,
>> 9816860.691, 9816863.364, 9816868.255, 9816870.982, 9816876.447,
>> 9816879.145, 9816883.987, 9816886.391, 9816890.726, 9816892.896,
>> 9816896.929, 9816901.867, 9816916.037, 9816923.885, 9816939.448,
>> 9816945.84, 9816955.292, 9816957.643, 9816358.103, 9816363.81,
>> 9816386.643, 9816396.126, 9816411.121, 9816418.727, 9816433.986,
>> 9816441.755, 9816457.333, 9816465.293, 9816481.245, 9816489.433,
>> 9816505.847, 9816514.286, 9816531.11, 9816544.301, 9816578.838,
>> 9816606.198, 9816675.659, 9816697.74, 9816714.884, 9816723.852,
>> 9816740.671, 9816749.402, 9816765.672, 9816774.048, 9816789.547,
>> 9816797.45, 9816811.896, 9816817.382, 9816824.195, 9816827.697,
>> 9816834.178, 9816837.496, 9816843.618, 9816846.738, 9816852.472,
>> 9816855.379, 9816860.698, 9816863.378, 9816868.262, 9816870.997,
>> 9816876.454, 9816879.158, 9816883.993, 9816886.403, 9816890.731,
>> 9816892.907, 9816896.942, 9816901.906, 9816916.056, 9816923.925,
>> 9816939.469, 9816946.399, 9816957.888, 9816960.745, 9816318.255,
>> 9816323.251, 9816343.25, 9816355.642, 9816385.023, 9816396.165,
>> 9816411.141, 9816418.767, 9816434.006, 9816441.796, 9816457.353,
>> 9816465.335, 9816481.266, 9816489.477, 9816505.869, 9816514.331,
>> 9816531.145, 9816544.394, 9816578.91, 9816606.394, 9816675.72,
>> 9816697.788, 9816714.908, 9816723.899, 9816740.694, 9816749.447,
>> 9816765.694, 9816774.091, 9816789.568, 9816797.491, 9816811.911,
>> 9816817.4, 9816824.204, 9816827.715, 9816834.187, 9816837.513,
>> 9816843.626, 9816846.754, 9816852.48, 9816855.394, 9816860.705,
>> 9816863.391, 9816868.269, 9816871.012, 9816876.461, 9816879.172,
>> 9816884, 9816886.415, 9816890.737, 9816892.917, 9816896.954,
>> 9816901.945, 9816916.076, 9816923.965, 9816939.489, 9816946.958,
>> 9816960.484, 9816963.848, 9816969.359, 9816969.355, 9816969.279,
>> 9816969.038, 9816968.979, 9816969.359, 9816969.345, 9816969.078,
>> 9816968.233, 9816968.025, 9816969.359, 9816969.334, 9816968.878,
>> 9816967.428, 9816967.071, 9816969.359, 9816969.324, 9816968.677,
>> 9816966.622, 9816966.118, 9816969.359, 9816969.313, 9816968.476,
>> 9816965.817, 9816965.164, 9816969.359, 9816969.303, 9816968.275,
>> 9816965.012, 9816964.21, 9816969.359, 9816969.292, 9816968.074,
>> 9816964.207, 9816963.256, 9816969.359, 9816969.282, 9816967.874,
>> 9816963.401, 9816962.303, 9816969.359, 9816969.271, 9816967.673,
>> 9816962.596, 9816961.349, 9816969.359, 9816969.261, 9816967.472,
>> 9816961.791, 9816960.395, 9816969.359, 9816969.25, 9816967.271,
>> 9816960.985, 9816959.441, 9816969.359, 9816969.24, 9816967.07,
>> 9816960.18, 9816958.488, 9816969.359, 9816969.23, 9816966.87,
>> 9816959.375, 9816957.534, 9816969.359, 9816969.219, 9816966.669,
>> 9816958.57, 9816956.58, 9816969.359, 9816969.209, 9816966.468,
>> 9816957.764, 9816955.626, 9816969.359, 9816969.198, 9816966.267,
>> 9816956.959, 9816954.672, 9816969.359, 9816969.188, 9816966.066,
>> 9816956.154, 9816953.719, 9816969.359, 9816969.177, 9816965.865,
>> 9816955.349, 9816952.765, 9816969.359, 9816969.167, 9816965.665,
>> 9816954.543, 9816951.811, 9816969.359, 9816969.156, 9816965.464,
>> 9816953.738, 9816950.857, 9816969.359, 9816969.146, 9816965.263,
>> 9816952.933, 9816949.904, 9816969.359, 9816969.135, 9816965.062,
>> 9816952.128, 9816948.95, 9816969.359, 9816969.125, 9816964.861,
>> 9816951.322, 9816947.996, 9816969.359, 9816969.114, 9816964.661,
>> 9816950.517, 9816947.042, 9816969.359, 9816969.104, 9816964.46,
>> 9816949.712, 9816946.089, 9816969.359, 9816969.094, 9816964.259,
>> 9816948.907, 9816945.135, 9816969.359, 9816969.083, 9816964.058,
>> 9816948.101, 9816944.181, 9816969.359, 9816969.073, 9816963.861,
>> 9816947.311, 9816943.242, 9816943.232, 9816943.229, 9816969.359,
>> 9816969.062, 9816963.819, 9816947.291, 9816943.057, 9816942.495,
>> 9816942.357, 9816969.359, 9816969.052, 9816963.777, 9816947.27,
>> 9816942.872, 9816941.759, 9816941.486, 9816969.359, 9816969.041,
>> 9816963.735, 9816947.25, 9816942.686, 9816941.023, 9816940.614,
>> 9816969.359, 9816969.031, 9816963.693, 9816947.229, 9816942.501,
>> 9816940.286, 9816939.743, 9816969.359, 9816969.02, 9816963.651,
>> 9816947.209, 9816942.316, 9816939.55, 9816938.871, 9816969.359,
>> 9816969.01, 9816963.609, 9816947.189, 9816942.13, 9816938.814,
>> 9816938, 9816969.359, 9816968.999, 9816963.568, 9816947.168,
>> 9816941.945, 9816938.078, 9816937.128, 9816969.359, 9816968.989,
>> 9816963.526, 9816947.148, 9816941.76, 9816937.341, 9816936.257,
>> 9816969.359, 9816968.978, 9816963.484, 9816947.127, 9816941.575,
>> 9816936.605, 9816935.385, 9816969.359, 9816968.968, 9816963.442,
>> 9816947.107, 9816941.389, 9816935.869, 9816934.514, 9816969.359,
>> 9816968.958, 9816963.4, 9816947.086, 9816941.204, 9816935.133,
>> 9816933.642, 9816969.359, 9816968.947, 9816963.358, 9816947.066,
>> 9816941.019, 9816934.396, 9816932.771, 9816969.359, 9816968.937,
>> 9816963.317, 9816947.045, 9816940.833, 9816933.66, 9816931.899,
>> 9816969.359, 9816968.926, 9816963.275, 9816947.025, 9816940.648,
>> 9816932.924, 9816931.028, 9816969.359, 9816968.916, 9816963.233,
>> 9816947.004, 9816940.463, 9816932.187, 9816930.156, 9816969.359,
>> 9816968.905, 9816963.191, 9816946.984, 9816940.277, 9816931.451,
>> 9816929.285, 9816969.359, 9816968.895, 9816963.149, 9816946.963,
>> 9816940.092, 9816930.715, 9816928.413, 9816969.359, 9816968.884,
>> 9816963.107, 9816946.943, 9816939.907, 9816929.979, 9816927.542,
>> 9816969.359, 9816968.874, 9816963.066, 9816946.922, 9816939.721,
>> 9816929.242, 9816926.67, 9816969.359, 9816968.863, 9816963.024,
>> 9816946.902, 9816939.536, 9816928.506, 9816925.799, 9816969.359,
>> 9816968.853, 9816962.982, 9816946.881, 9816939.351, 9816927.77,
>> 9816924.927, 9816969.359, 9816968.842, 9816962.94, 9816946.861,
>> 9816939.165, 9816927.034, 9816924.056, 9816969.359, 9816968.832,
>> 9816962.898, 9816946.84, 9816938.98, 9816926.297, 9816923.184,
>> 9816969.359, 9816968.822, 9816962.856, 9816946.82, 9816938.795,
>> 9816925.561, 9816922.313, 9816969.359, 9816968.811, 9816962.814,
>> 9816946.799, 9816938.61, 9816924.825, 9816921.441, 9816969.359,
>> 9816968.801, 9816962.773, 9816946.779, 9816938.424, 9816924.089,
>> 9816920.57, 9816969.359, 9816968.79, 9816962.731, 9816946.759,
>> 9816938.239, 9816923.352, 9816919.698, 9816969.359, 9816968.78,
>> 9816962.689, 9816946.738, 9816938.147, 9816923.075, 9816915.758,
>> 9816901.367, 9816896.756, 9816892.751, 9816890.633, 9816886.232,
>> 9816883.882, 9816878.967, 9816876.326, 9816870.783, 9816868.135,
>> 9816863.186, 9816860.574, 9816855.172, 9816852.344, 9816846.518,
>> 9816843.488, 9816837.264, 9816834.046, 9816827.453, 9816824.063,
>> 9816817.127, 9816814.573, 9816811.312, 9816810.498, 9816969.359,
>> 9816968.769, 9816962.647, 9816946.718, 9816938.107, 9816923.055,
>> 9816915.719, 9816901.355, 9816896.745, 9816892.746, 9816890.621,
>> 9816886.226, 9816883.869, 9816878.96, 9816876.311, 9816870.776,
>> 9816868.121, 9816863.179, 9816860.56, 9816855.165, 9816852.329,
>> 9816846.509, 9816843.471, 9816837.255, 9816834.028, 9816827.444,
>> 9816824.044, 9816817.112, 9816811.653, 9816796.899, 9816791.307,
>> 9816783.741, 9816781.85, 9816969.359, 9816968.759, 9816962.605,
>> 9816946.697, 9816938.067, 9816923.035, 9816915.68, 9816901.342,
>> 9816896.734, 9816892.74, 9816890.61, 9816886.219, 9816883.856,
>> 9816878.953, 9816876.296, 9816870.769, 9816868.108, 9816863.172,
>> 9816860.545, 9816855.157, 9816852.313, 9816846.501, 9816843.454,
>> 9816837.246, 9816834.01, 9816827.434, 9816824.025, 9816817.098,
>> 9816811.612, 9816796.878, 9816789.209, 9816773.449, 9816766.432,
>> 9816754.204, 9816751.148, 9816969.359, 9816968.748, 9816962.563,
>> 9816946.677, 9816938.026, 9816923.016, 9816915.641, 9816901.33,
>> 9816896.723, 9816892.734, 9816890.598, 9816886.213, 9816883.842,
>> 9816878.946, 9816876.281, 9816870.762, 9816868.094, 9816863.165,
>> 9816860.53, 9816855.15, 9816852.297, 9816846.493, 9816843.437,
>> 9816837.237, 9816833.992, 9816827.425, 9816824.006, 9816817.083,
>> 9816811.572, 9816796.857, 9816789.166, 9816773.427, 9816765.303,
>> 9816748.764, 9816740.333, 9816723.262, 9816718.996, 9816969.359,
>> 9816968.738, 9816962.522, 9816946.656, 9816937.986, 9816922.996,
>> 9816915.602, 9816901.317, 9816896.712, 9816892.728, 9816890.586,
>> 9816886.207, 9816883.829, 9816878.939, 9816876.266, 9816870.755,
>> 9816868.081, 9816863.158, 9816860.516, 9816855.142, 9816852.281,
>> 9816846.485, 9816843.42, 9816837.229, 9816833.975, 9816827.416,
>> 9816823.987, 9816817.068, 9816811.532, 9816796.836, 9816789.123,
>> 9816773.405, 9816765.258, 9816748.741, 9816740.275, 9816723.183,
>> 9816714.514, 9816697.065, 9816675.046, 9816604.224, 9816576.996,
>> 9816543.212, 9816529.997, 9816513.594, 9816504.709, 9816488.747,
>> 9816480.024, 9816464.605, 9816456.114, 9816441.085, 9816432.766,
>> 9816418.074, 9816409.898, 9816395.185, 9816382.56, 9816354.343,
>> 9816338.456, 9816311.391, 9816304.784, 9816969.359, 9816968.727,
>> 9816962.48, 9816946.636, 9816937.946, 9816922.976, 9816915.563,
>> 9816901.305, 9816896.702, 9816892.723, 9816890.574, 9816886.201,
>> 9816883.816, 9816878.932, 9816876.251, 9816870.748, 9816868.068,
>> 9816863.151, 9816860.501, 9816855.134, 9816852.265, 9816846.477,
>> 9816843.403, 9816837.22, 9816833.957, 9816827.407, 9816823.969,
>> 9816817.053, 9816811.492, 9816796.815, 9816789.08, 9816773.383,
>> 9816765.213, 9816748.718, 9816740.228, 9816723.16, 9816714.466,
>> 9816697.004, 9816674.85, 9816604.151, 9816576.903, 9816543.177,
>> 9816529.952, 9816513.572, 9816504.666, 9816488.725, 9816479.982,
>> 9816464.584, 9816456.073, 9816441.065, 9816432.726, 9816418.054,
>> 9816409.858, 9816395.156, 9816382.483, 9816354.304, 9816338.38,
>> 9816311.372, 9816304.784, 9816969.359, 9816968.717, 9816962.438,
>> 9816946.615, 9816937.906, 9816922.956, 9816915.524, 9816901.292,
>> 9816896.691, 9816892.717, 9816890.562, 9816886.194, 9816883.802,
>> 9816878.925, 9816876.236, 9816870.741, 9816868.054, 9816863.144,
>> 9816860.486, 9816855.127, 9816852.25, 9816846.469, 9816843.386,
>> 9816837.211, 9816833.939, 9816827.398, 9816823.95, 9816817.039,
>> 9816811.452, 9816796.795, 9816789.037, 9816773.361, 9816765.168,
>> 9816748.695, 9816740.182, 9816723.136, 9816714.418, 9816696.943,
>> 9816674.654, 9816604.079, 9816576.81, 9816543.143, 9816529.907,
>> 9816513.549, 9816504.622, 9816488.704, 9816479.94, 9816464.564,
>> 9816456.032, 9816441.045, 9816432.686, 9816418.034, 9816409.819,
>> 9816395.127, 9816382.407, 9816354.266, 9816338.305, 9816311.353,
>> 9816304.784, 9816969.359, 9816968.706, 9816962.396, 9816946.595,
>> 9816937.866, 9816922.936, 9816915.485, 9816901.28, 9816896.68,
>> 9816892.711, 9816890.55, 9816886.188, 9816883.789, 9816878.918,
>> 9816876.221, 9816870.734, 9816868.041, 9816863.137, 9816860.472,
>> 9816855.119, 9816852.234, 9816846.46, 9816843.369, 9816837.203,
>> 9816833.921, 9816827.389, 9816823.931, 9816817.024, 9816811.412,
>> 9816796.774, 9816788.994, 9816773.339, 9816765.122, 9816748.672,
>> 9816740.135, 9816723.112, 9816714.371, 9816696.882, 9816674.458,
>> 9816604.007, 9816576.717, 9816543.109, 9816529.863, 9816513.527,
>> 9816504.579, 9816488.682, 9816479.898, 9816464.543, 9816455.991,
>> 9816441.025, 9816432.646, 9816418.015, 9816409.78, 9816395.098,
>> 9816382.33, 9816354.228, 9816338.229, 9816311.335, 9816304.784,
>> 9816969.359, 9816968.696, 9816962.354, 9816946.574, 9816937.826,
>> 9816922.917, 9816915.446, 9816901.267, 9816896.669, 9816892.706,
>> 9816890.538, 9816886.182, 9816883.776, 9816878.911, 9816876.206,
>> 9816870.727, 9816868.027, 9816863.13, 9816860.457, 9816855.112,
>> 9816852.218, 9816846.452, 9816843.353, 9816837.194, 9816833.903,
>> 9816827.379, 9816823.912, 9816817.009, 9816811.372, 9816796.753,
>> 9816788.951, 9816773.317, 9816765.077, 9816748.649, 9816740.088,
>> 9816723.089, 9816714.323, 9816696.821, 9816674.262, 9816603.934,
>> 9816576.623, 9816543.074, 9816529.818, 9816513.505, 9816504.535,
>> 9816488.661, 9816479.856, 9816464.522, 9816455.95, 9816441.005,
>> 9816432.606, 9816417.995, 9816409.741, 9816395.069, 9816382.253,
>> 9816354.19, 9816338.153, 9816311.316, 9816304.784, 9816969.359,
>> 9816968.686, 9816962.312, 9816946.554, 9816937.786, 9816922.897,
>> 9816915.407), Y = c(6685246.262, 6685259.661, 6685313.533, 6685342.653,
>> 6685402.003, 6685433.827, 6685494.079, 6685515.065, 6685530.452,
>> 6685538.336, 6685553.068, 6685560.845, 6685575.358, 6685583.002,
>> 6685597.245, 6685604.053, 6685615.247, 6685621.127, 6685632.091,
>> 6685637.868, 6685648.642, 6685654.322, 6685664.913, 6685670.499,
>> 6685680.917, 6685686.414, 6685696.667, 6685704.514, 6685724.452,
>> 6685734.976, 6685754.466, 6685764.715, 6685783.716, 6685793.726,
>> 6685812.306, 6685822.114, 6685840.341, 6685863.242, 6685934.549,
>> 6685963.196, 6685999.539, 6686013.753, 6686032.223, 6686041.737,
>> 6686060.477, 6686070.119, 6686089.207, 6686099.015, 6686118.482,
>> 6686128.449, 6686148.244, 6686158.343, 6686178.415, 6686193.73,
>> 6686234.484, 6686255.141, 6686296.412, 6686306.785, 6686306.803,
>> 6685430.056, 6685428.329, 6685421.421, 6685419.586, 6685419.183,
>> 6685419.082, 6685941.655, 6685925.678, 6685861.712, 6685840.901,
>> 6685822.396, 6685812.876, 6685794.014, 6685784.297, 6685765.009,
>> 6685755.061, 6685735.277, 6685725.063, 6685704.745, 6685696.98,
>> 6685686.572, 6685681.234, 6685670.659, 6685665.235, 6685654.484,
>> 6685648.969, 6685638.033, 6685632.423, 6685621.294, 6685615.585,
>> 6685604.246, 6685597.676, 6685583.22, 6685575.798, 6685561.067,
>> 6685553.515, 6685538.561, 6685530.905, 6685515.637, 6685495.914,
>> 6685434.708, 6685409.771, 6685373.751, 6685364.795, 6685921.162,
>> 6685908.601, 6685858.31, 6685840.952, 6685822.422, 6685812.928,
>> 6685794.04, 6685784.35, 6685765.036, 6685755.115, 6685735.305,
>> 6685725.118, 6685704.766, 6685697.008, 6685686.586, 6685681.263,
>> 6685670.673, 6685665.264, 6685654.499, 6685648.999, 6685638.048,
>> 6685632.454, 6685621.31, 6685615.616, 6685604.263, 6685597.715,
>> 6685583.24, 6685575.838, 6685561.087, 6685553.555, 6685538.581,
>> 6685530.946, 6685515.689, 6685496.08, 6685434.788, 6685411.734,
>> 6685382.838, 6685375.652, 6685900.669, 6685891.523, 6685854.907,
>> 6685841.003, 6685822.448, 6685812.98, 6685794.066, 6685784.403,
>> 6685765.063, 6685755.169, 6685735.332, 6685725.174, 6685704.787,
>> 6685697.036, 6685686.6, 6685681.292, 6685670.688, 6685665.293,
>> 6685654.513, 6685649.029, 6685638.063, 6685632.484, 6685621.325,
>> 6685615.646, 6685604.281, 6685597.755, 6685583.259, 6685575.878,
>> 6685561.107, 6685553.596, 6685538.602, 6685530.987, 6685515.741,
>> 6685496.247, 6685434.868, 6685413.697, 6685391.924, 6685386.51,
>> 6685880.176, 6685874.446, 6685851.504, 6685841.054, 6685822.473,
>> 6685813.031, 6685794.093, 6685784.456, 6685765.09, 6685755.224,
>> 6685735.36, 6685725.229, 6685704.808, 6685697.065, 6685686.615,
>> 6685681.321, 6685670.702, 6685665.323, 6685654.528, 6685649.058,
>> 6685638.078, 6685632.514, 6685621.34, 6685615.677, 6685604.298,
>> 6685597.794, 6685583.279, 6685575.918, 6685561.127, 6685553.637,
>> 6685538.622, 6685531.028, 6685515.793, 6685496.414, 6685434.948,
>> 6685415.66, 6685401.01, 6685397.367, 6685859.683, 6685857.368,
>> 6685848.101, 6685841.105, 6685822.499, 6685813.083, 6685794.119,
>> 6685784.509, 6685765.116, 6685755.278, 6685735.387, 6685725.285,
>> 6685704.829, 6685697.093, 6685686.629, 6685681.349, 6685670.717,
>> 6685665.352, 6685654.543, 6685649.088, 6685638.093, 6685632.544,
>> 6685621.355, 6685615.708, 6685604.316, 6685597.833, 6685583.299,
>> 6685575.958, 6685561.148, 6685553.677, 6685538.643, 6685531.069,
>> 6685515.845, 6685496.581, 6685435.028, 6685417.623, 6685410.096,
>> 6685408.225, 6686003.197, 6685995.133, 6685962.757, 6685936.166,
>> 6685863.748, 6685840.748, 6685822.319, 6685812.72, 6685793.936,
>> 6685784.139, 6685764.929, 6685754.899, 6685735.195, 6685724.896,
>> 6685704.682, 6685696.895, 6685686.529, 6685681.148, 6685670.615,
>> 6685665.147, 6685654.44, 6685648.88, 6685637.988, 6685632.333,
>> 6685621.249, 6685615.493, 6685604.193, 6685597.558, 6685583.16,
>> 6685575.678, 6685561.007, 6685553.393, 6685538.5, 6685530.781,
>> 6685515.481, 6685495.413, 6685434.468, 6685403.881, 6685346.492,
>> 6685332.222, 6685982.677, 6685978.033, 6685959.387, 6685936.368,
>> 6685863.811, 6685840.799, 6685822.345, 6685812.772, 6685793.962,
>> 6685784.192, 6685764.956, 6685754.953, 6685735.223, 6685724.952,
>> 6685704.703, 6685696.923, 6685686.543, 6685681.177, 6685670.63,
>> 6685665.176, 6685654.454, 6685648.91, 6685638.003, 6685632.363,
>> 6685621.264, 6685615.524, 6685604.211, 6685597.598, 6685583.18,
>> 6685575.718, 6685561.027, 6685553.434, 6685538.52, 6685530.822,
>> 6685515.533, 6685495.58, 6685434.548, 6685405.845, 6685355.579,
>> 6685343.08, 6685962.158, 6685960.933, 6685956.018, 6685936.57,
>> 6685863.874, 6685840.85, 6685822.371, 6685812.824, 6685793.988,
>> 6685784.245, 6685764.983, 6685755.007, 6685735.25, 6685725.007,
>> 6685704.724, 6685696.951, 6685686.557, 6685681.205, 6685670.644,
>> 6685665.206, 6685654.469, 6685648.939, 6685638.018, 6685632.393,
>> 6685621.279, 6685615.554, 6685604.228, 6685597.637, 6685583.2,
>> 6685575.758, 6685561.047, 6685553.474, 6685538.541, 6685530.864,
>> 6685515.585, 6685495.747, 6685434.628, 6685407.808, 6685364.665,
>> 6685353.937, 6686023.971, 6686021.541, 6686011.804, 6686000.226,
>> 6685963.722, 6685935.964, 6685863.685, 6685840.697, 6685822.293,
>> 6685812.669, 6685793.909, 6685784.086, 6685764.902, 6685754.845,
>> 6685735.168, 6685724.84, 6685704.661, 6685696.866, 6685686.514,
>> 6685681.119, 6685670.601, 6685665.118, 6685654.425, 6685648.85,
>> 6685637.973, 6685632.303, 6685621.233, 6685615.462, 6685604.176,
>> 6685597.519, 6685583.141, 6685575.638, 6685560.986, 6685553.352,
>> 6685538.479, 6685530.74, 6685515.429, 6685495.246, 6685434.388,
>> 6685403.077, 6685343.175, 6685327.077, 6685322.563, 6685321.44,
>> 6686044.994, 6686043.689, 6686038.466, 6686032.518, 6686013.974,
>> 6686000.128, 6685963.647, 6685935.762, 6685863.621, 6685840.646,
>> 6685822.268, 6685812.617, 6685793.883, 6685784.033, 6685764.876,
>> 6685754.791, 6685735.14, 6685724.785, 6685704.64, 6685696.838,
>> 6685686.5, 6685681.09, 6685670.586, 6685665.089, 6685654.41,
>> 6685648.82, 6685637.958, 6685632.272, 6685621.218, 6685615.431,
>> 6685604.158, 6685597.48, 6685583.121, 6685575.598, 6685560.966,
>> 6685553.312, 6685538.459, 6685530.699, 6685515.377, 6685495.08,
>> 6685434.308, 6685402.923, 6685343.1, 6685325.142, 6685313.577,
>> 6685310.7, 6686067.455, 6686067.12, 6686065.775, 6686060.726,
>> 6686041.86, 6686032.469, 6686013.937, 6686000.03, 6685963.571,
>> 6685935.559, 6685863.558, 6685840.595, 6685822.242, 6685812.565,
>> 6685793.857, 6685783.98, 6685764.849, 6685754.737, 6685735.113,
>> 6685724.729, 6685704.619, 6685696.809, 6685686.486, 6685681.061,
>> 6685670.572, 6685665.059, 6685654.395, 6685648.791, 6685637.943,
>> 6685632.242, 6685621.203, 6685615.401, 6685604.141, 6685597.441,
>> 6685583.101, 6685575.558, 6685560.946, 6685553.271, 6685538.438,
>> 6685530.658, 6685515.325, 6685494.913, 6685434.228, 6685402.77,
>> 6685343.026, 6685323.207, 6685304.591, 6685299.961, 6686120.566,
>> 6686116.177, 6686098.614, 6686089.408, 6686070.22, 6686060.676,
>> 6686041.836, 6686032.42, 6686013.9, 6685999.932, 6685963.496,
>> 6685935.357, 6685863.495, 6685840.545, 6685822.216, 6685812.513,
>> 6685793.831, 6685783.927, 6685764.822, 6685754.682, 6685735.086,
>> 6685724.674, 6685704.598, 6685696.781, 6685686.472, 6685681.032,
>> 6685670.557, 6685665.03, 6685654.381, 6685648.761, 6685637.928,
>> 6685632.212, 6685621.188, 6685615.37, 6685604.123, 6685597.401,
>> 6685583.081, 6685575.518, 6685560.926, 6685553.23, 6685538.418,
>> 6685530.617, 6685515.273, 6685494.746, 6685434.148, 6685402.617,
>> 6685342.951, 6685321.272, 6685295.605, 6685289.221, 6686175.29,
>> 6686171.64, 6686157.035, 6686148.4, 6686128.526, 6686118.635,
>> 6686099.091, 6686089.358, 6686070.194, 6686060.626, 6686041.811,
>> 6686032.371, 6686013.864, 6685999.833, 6685963.421, 6685935.155,
>> 6685863.431, 6685840.494, 6685822.191, 6685812.461, 6685793.805,
>> 6685783.875, 6685764.795, 6685754.628, 6685735.058, 6685724.618,
>> 6685704.577, 6685696.753, 6685686.457, 6685681.004, 6685670.543,
>> 6685665.001, 6685654.366, 6685648.731, 6685637.913, 6685632.182,
>> 6685621.172, 6685615.339, 6685604.106, 6685597.362, 6685583.061,
>> 6685575.478, 6685560.906, 6685553.19, 6685538.397, 6685530.576,
>> 6685515.221, 6685494.579, 6685434.068, 6685402.463, 6685342.877,
>> 6685319.338, 6685286.619, 6685278.481, 6686231.171, 6686223.242,
>> 6686191.517, 6686178.52, 6686158.396, 6686148.348, 6686128.5,
>> 6686118.584, 6686099.066, 6686089.307, 6686070.169, 6686060.576,
>> 6686041.786, 6686032.321, 6686013.827, 6685999.735, 6685963.346,
>> 6685934.953, 6685863.368, 6685840.443, 6685822.165, 6685812.41,
>> 6685793.779, 6685783.822, 6685764.769, 6685754.574, 6685735.031,
>> 6685724.563, 6685704.556, 6685696.724, 6685686.443, 6685680.975,
>> 6685670.528, 6685664.972, 6685654.351, 6685648.702, 6685637.898,
>> 6685632.151, 6685621.157, 6685615.308, 6685604.088, 6685597.323,
>> 6685583.042, 6685575.438, 6685560.886, 6685553.149, 6685538.377,
>> 6685530.535, 6685515.169, 6685494.412, 6685433.988, 6685402.31,
>> 6685342.802, 6685317.403, 6685277.633, 6685267.741, 6686287.605,
>> 6686280.486, 6686251.986, 6686234.591, 6686193.769, 6686178.467,
>> 6686158.369, 6686148.296, 6686128.475, 6686118.533, 6686099.04,
>> 6686089.257, 6686070.144, 6686060.526, 6686041.761, 6686032.272,
>> 6686013.79, 6685999.637, 6685963.271, 6685934.751, 6685863.305,
>> 6685840.392, 6685822.139, 6685812.358, 6685793.752, 6685783.769,
>> 6685764.742, 6685754.52, 6685735.003, 6685724.507, 6685704.535,
>> 6685696.696, 6685686.429, 6685680.946, 6685670.514, 6685664.942,
>> 6685654.336, 6685648.672, 6685637.883, 6685632.121, 6685621.142,
>> 6685615.278, 6685604.071, 6685597.284, 6685583.022, 6685575.398,
>> 6685560.865, 6685553.108, 6685538.357, 6685530.493, 6685515.117,
>> 6685494.246, 6685433.908, 6685402.156, 6685342.728, 6685315.468,
>> 6685268.647, 6685257.001, 6685237.925, 6685237.939, 6685238.202,
>> 6685239.036, 6685239.241, 6685237.925, 6685237.975, 6685238.897,
>> 6685241.824, 6685242.543, 6685237.925, 6685238.012, 6685239.592,
>> 6685244.611, 6685245.844, 6685237.925, 6685238.048, 6685240.287,
>> 6685247.398, 6685249.145, 6685237.925, 6685238.084, 6685240.982,
>> 6685250.186, 6685252.447, 6685237.925, 6685238.12, 6685241.677,
>> 6685252.973, 6685255.748, 6685237.925, 6685238.156, 6685242.372,
>> 6685255.761, 6685259.05, 6685237.925, 6685238.193, 6685243.067,
>> 6685258.548, 6685262.351, 6685237.925, 6685238.229, 6685243.763,
>> 6685261.336, 6685265.653, 6685237.925, 6685238.265, 6685244.458,
>> 6685264.123, 6685268.954, 6685237.925, 6685238.301, 6685245.153,
>> 6685266.91, 6685272.256, 6685237.925, 6685238.337, 6685245.848,
>> 6685269.698, 6685275.557, 6685237.925, 6685238.374, 6685246.543,
>> 6685272.485, 6685278.858, 6685237.925, 6685238.41, 6685247.238,
>> 6685275.273, 6685282.16, 6685237.925, 6685238.446, 6685247.933,
>> 6685278.06, 6685285.461, 6685237.925, 6685238.482, 6685248.628,
>> 6685280.847, 6685288.763, 6685237.925, 6685238.519, 6685249.323,
>> 6685283.635, 6685292.064, 6685237.925, 6685238.555, 6685250.018,
>> 6685286.422, 6685295.366, 6685237.925, 6685238.591, 6685250.713,
>> 6685289.21, 6685298.667, 6685237.925, 6685238.627, 6685251.408,
>> 6685291.997, 6685301.969, 6685237.925, 6685238.663, 6685252.104,
>> 6685294.784, 6685305.27, 6685237.925, 6685238.7, 6685252.799,
>> 6685297.572, 6685308.571, 6685237.925, 6685238.736, 6685253.494,
>> 6685300.359, 6685311.873, 6685237.925, 6685238.772, 6685254.189,
>> 6685303.147, 6685315.174, 6685237.925, 6685238.808, 6685254.884,
>> 6685305.934, 6685318.476, 6685237.925, 6685238.844, 6685255.579,
>> 6685308.722, 6685321.777, 6685237.925, 6685238.881, 6685256.274,
>> 6685311.509, 6685325.079, 6685237.925, 6685238.917, 6685256.958,
>> 6685314.243, 6685328.329, 6685328.371, 6685328.381, 6685237.925,
>> 6685238.953, 6685257.103, 6685314.318, 6685329.039, 6685331.19,
>> 6685331.718, 6685237.925, 6685238.989, 6685257.248, 6685314.392,
>> 6685329.749, 6685334.01, 6685335.056, 6685237.925, 6685239.025,
>> 6685257.393, 6685314.467, 6685330.458, 6685336.83, 6685338.394,
>> 6685237.925, 6685239.062, 6685257.538, 6685314.541, 6685331.168,
>> 6685339.65, 6685341.731, 6685237.925, 6685239.098, 6685257.683,
>> 6685314.616, 6685331.878, 6685342.469, 6685345.069, 6685237.925,
>> 6685239.134, 6685257.827, 6685314.691, 6685332.587, 6685345.289,
>> 6685348.407, 6685237.925, 6685239.17, 6685257.972, 6685314.765,
>> 6685333.297, 6685348.109, 6685351.744, 6685237.925, 6685239.207,
>> 6685258.117, 6685314.84, 6685334.007, 6685350.928, 6685355.082,
>> 6685237.925, 6685239.243, 6685258.262, 6685314.914, 6685334.717,
>> 6685353.748, 6685358.419, 6685237.925, 6685239.279, 6685258.407,
>> 6685314.989, 6685335.426, 6685356.568, 6685361.757, 6685237.925,
>> 6685239.315, 6685258.552, 6685315.063, 6685336.136, 6685359.387,
>> 6685365.095, 6685237.925, 6685239.351, 6685258.696, 6685315.138,
>> 6685336.846, 6685362.207, 6685368.432, 6685237.925, 6685239.388,
>> 6685258.841, 6685315.213, 6685337.555, 6685365.027, 6685371.77,
>> 6685237.925, 6685239.424, 6685258.986, 6685315.287, 6685338.265,
>> 6685367.847, 6685375.108, 6685237.925, 6685239.46, 6685259.131,
>> 6685315.362, 6685338.975, 6685370.666, 6685378.445, 6685237.925,
>> 6685239.496, 6685259.276, 6685315.436, 6685339.684, 6685373.486,
>> 6685381.783, 6685237.925, 6685239.532, 6685259.421, 6685315.511,
>> 6685340.394, 6685376.306, 6685385.12, 6685237.925, 6685239.569,
>> 6685259.566, 6685315.585, 6685341.104, 6685379.125, 6685388.458,
>> 6685237.925, 6685239.605, 6685259.71, 6685315.66, 6685341.813,
>> 6685381.945, 6685391.796, 6685237.925, 6685239.641, 6685259.855,
>> 6685315.735, 6685342.523, 6685384.765, 6685395.133, 6685237.925,
>> 6685239.677, 6685260, 6685315.809, 6685343.233, 6685387.585,
>> 6685398.471, 6685237.925, 6685239.713, 6685260.145, 6685315.884,
>> 6685343.943, 6685390.404, 6685401.808, 6685237.925, 6685239.75,
>> 6685260.29, 6685315.958, 6685344.652, 6685393.224, 6685405.146,
>> 6685237.925, 6685239.786, 6685260.435, 6685316.033, 6685345.362,
>> 6685396.044, 6685408.484, 6685237.925, 6685239.822, 6685260.579,
>> 6685316.107, 6685346.072, 6685398.863, 6685411.821, 6685237.925,
>> 6685239.858, 6685260.724, 6685316.182, 6685346.781, 6685401.683,
>> 6685415.159, 6685237.925, 6685239.895, 6685260.869, 6685316.256,
>> 6685347.491, 6685404.503, 6685418.497, 6685237.925, 6685239.931,
>> 6685261.014, 6685316.331, 6685347.844, 6685405.568, 6685435.199,
>> 6685496.716, 6685515.938, 6685531.122, 6685538.725, 6685553.729,
>> 6685561.236, 6685576.008, 6685583.392, 6685597.877, 6685604.396,
>> 6685615.747, 6685621.422, 6685632.578, 6685638.148, 6685649.116,
>> 6685654.585, 6685665.375, 6685670.748, 6685681.367, 6685686.649,
>> 6685697.108, 6685700.858, 6685705.372, 6685706.5, 6685237.925,
>> 6685239.967, 6685261.159, 6685316.406, 6685347.998, 6685405.648,
>> 6685435.365, 6685496.768, 6685515.979, 6685531.142, 6685538.766,
>> 6685553.749, 6685561.276, 6685576.028, 6685583.432, 6685597.895,
>> 6685604.426, 6685615.762, 6685621.452, 6685632.593, 6685638.178,
>> 6685649.131, 6685654.615, 6685665.39, 6685670.777, 6685681.382,
>> 6685686.677, 6685697.129, 6685704.901, 6685725.326, 6685732.83,
>> 6685742.362, 6685744.744, 6685237.925, 6685240.003, 6685261.304,
>> 6685316.48, 6685348.151, 6685405.728, 6685435.532, 6685496.82,
>> 6685516.02, 6685531.163, 6685538.806, 6685553.769, 6685561.316,
>> 6685576.048, 6685583.471, 6685597.912, 6685604.457, 6685615.777,
>> 6685621.483, 6685632.608, 6685638.207, 6685649.146, 6685654.644,
>> 6685665.404, 6685670.806, 6685681.396, 6685686.706, 6685697.15,
>> 6685704.957, 6685725.354, 6685735.473, 6685755.328, 6685763.892,
>> 6685778.202, 6685781.778, 6685237.925, 6685240.039, 6685261.448,
>> 6685316.555, 6685348.305, 6685405.808, 6685435.699, 6685496.872,
>> 6685516.061, 6685531.183, 6685538.847, 6685553.789, 6685561.356,
>> 6685576.068, 6685583.51, 6685597.93, 6685604.488, 6685615.792,
>> 6685621.513, 6685632.623, 6685638.237, 6685649.16, 6685654.673,
>> 6685665.419, 6685670.835, 6685681.41, 6685686.734, 6685697.171,
>> 6685705.012, 6685725.381, 6685735.527, 6685755.355, 6685765.213,
>> 6685784.568, 6685794.162, 6685813.063, 6685817.786, 6685237.925,
>> 6685240.076, 6685261.593, 6685316.629, 6685348.458, 6685405.888,
>> 6685435.866, 6685496.924, 6685516.102, 6685531.204, 6685538.887,
>> 6685553.81, 6685561.396, 6685576.088, 6685583.549, 6685597.947,
>> 6685604.519, 6685615.808, 6685621.543, 6685632.638, 6685638.267,
>> 6685649.175, 6685654.702, 6685665.433, 6685670.863, 6685681.425,
>> 6685686.762, 6685697.192, 6685705.068, 6685725.409, 6685735.581,
>> 6685755.381, 6685765.266, 6685784.594, 6685794.226, 6685813.15,
>> 6685822.566, 6685841.162, 6685864.006, 6685936.995, 6685965.278,
>> 6686000.893, 6686015.045, 6686033.095, 6686043.087, 6686061.378,
>> 6686071.633, 6686090.149, 6686100.591, 6686119.437, 6686130.086,
>> 6686149.209, 6686160.037, 6686179.812, 6686197.19, 6686236.419,
>> 6686258.818, 6686297.386, 6686306.803, 6685237.925, 6685240.112,
>> 6685261.738, 6685316.704, 6685348.612, 6685405.968, 6685436.033,
>> 6685496.976, 6685516.143, 6685531.224, 6685538.928, 6685553.83,
>> 6685561.436, 6685576.107, 6685583.588, 6685597.965, 6685604.549,
>> 6685615.823, 6685621.573, 6685632.653, 6685638.296, 6685649.19,
>> 6685654.732, 6685665.448, 6685670.892, 6685681.439, 6685686.791,
>> 6685697.213, 6685705.123, 6685725.436, 6685735.635, 6685755.408,
>> 6685765.319, 6685784.62, 6685794.278, 6685813.176, 6685822.617,
>> 6685841.226, 6685864.208, 6685937.07, 6685965.376, 6686000.93,
>> 6686015.094, 6686033.119, 6686043.136, 6686061.403, 6686071.684,
>> 6686090.174, 6686100.643, 6686119.463, 6686130.138, 6686149.235,
>> 6686160.09, 6686179.852, 6686197.297, 6686236.472, 6686258.926,
>> 6686297.413, 6686306.803, 6685237.925, 6685240.148, 6685261.883,
>> 6685316.778, 6685348.765, 6685406.048, 6685436.199, 6685497.028,
>> 6685516.184, 6685531.244, 6685538.969, 6685553.85, 6685561.476,
>> 6685576.127, 6685583.628, 6685597.982, 6685604.58, 6685615.838,
>> 6685621.604, 6685632.668, 6685638.326, 6685649.205, 6685654.761,
>> 6685665.462, 6685670.921, 6685681.453, 6685686.819, 6685697.234,
>> 6685705.179, 6685725.464, 6685735.689, 6685755.435, 6685765.372,
>> 6685784.646, 6685794.33, 6685813.201, 6685822.668, 6685841.289,
>> 6685864.41, 6685937.145, 6685965.474, 6686000.967, 6686015.144,
>> 6686033.144, 6686043.186, 6686061.429, 6686071.734, 6686090.2,
>> 6686100.694, 6686119.488, 6686130.19, 6686149.262, 6686160.142,
>> 6686179.891, 6686197.403, 6686236.526, 6686259.033, 6686297.44,
>> 6686306.803, 6685237.925, 6685240.184, 6685262.028, 6685316.853,
>> 6685348.918, 6685406.128, 6685436.366, 6685497.08, 6685516.226,
>> 6685531.265, 6685539.009, 6685553.87, 6685561.516, 6685576.147,
>> 6685583.667, 6685598, 6685604.611, 6685615.853, 6685621.634,
>> 6685632.683, 6685638.356, 6685649.219, 6685654.79, 6685665.477,
>> 6685670.95, 6685681.468, 6685686.848, 6685697.255, 6685705.234,
>> 6685725.491, 6685735.743, 6685755.462, 6685765.425, 6685784.672,
>> 6685794.382, 6685813.227, 6685822.719, 6685841.352, 6685864.612,
>> 6685937.22, 6685965.572, 6686001.004, 6686015.193, 6686033.169,
>> 6686043.236, 6686061.454, 6686071.784, 6686090.225, 6686100.745,
>> 6686119.514, 6686130.242, 6686149.288, 6686160.195, 6686179.931,
>> 6686197.51, 6686236.58, 6686259.141, 6686297.467, 6686306.803,
>> 6685237.925, 6685240.22, 6685262.173, 6685316.928, 6685349.072,
>> 6685406.208, 6685436.533, 6685497.132, 6685516.267, 6685531.285,
>> 6685539.05, 6685553.89, 6685561.556, 6685576.167, 6685583.706,
>> 6685598.017, 6685604.641, 6685615.869, 6685621.664, 6685632.698,
>> 6685638.386, 6685649.234, 6685654.82, 6685665.491, 6685670.979,
>> 6685681.482, 6685686.876, 6685697.276, 6685705.29, 6685725.518,
>> 6685735.798, 6685755.488, 6685765.477, 6685784.698, 6685794.433,
>> 6685813.253, 6685822.77, 6685841.415, 6685864.814, 6685937.295,
>> 6685965.67, 6686001.041, 6686015.242, 6686033.194, 6686043.286,
>> 6686061.479, 6686071.835, 6686090.25, 6686100.796, 6686119.54,
>> 6686130.294, 6686149.314, 6686160.248, 6686179.971, 6686197.617,
>> 6686236.633, 6686259.249, 6686297.494, 6686306.803, 6685237.925,
>> 6685240.257, 6685262.318, 6685317.002, 6685349.225, 6685406.288,
>> 6685436.7), Z = c(3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L,
>> 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3600L, 3568L, 3568L,
>> 3568L, 3568L, 3568L, 3568L, 3578L, 3578L, 3578L, 3578L, 3578L,
>> 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L,
>> 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L,
>> 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3578L,
>> 3578L, 3578L, 3578L, 3578L, 3578L, 3578L, 3576L, 3576L, 3576L,
>> 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L,
>> 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L,
>> 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L,
>> 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3576L, 3574L,
>> 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L,
>> 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L,
>> 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L,
>> 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L, 3574L,
>> 3574L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L,
>> 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L,
>> 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L,
>> 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L, 3572L,
>> 3572L, 3572L, 3572L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L,
>> 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L,
>> 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L,
>> 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L, 3570L,
>> 3570L, 3570L, 3570L, 3570L, 3570L, 3584L, 3584L, 3584L, 3584L,
>> 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L,
>> 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L,
>> 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L,
>> 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L, 3584L,
>> 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L,
>> 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L,
>> 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L,
>> 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L, 3582L,
>> 3582L, 3582L, 3582L, 3582L, 3580L, 3580L, 3580L, 3580L, 3580L,
>> 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L,
>> 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L,
>> 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L,
>> 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3580L, 3586L,
>> 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L,
>> 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L,
>> 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L,
>> 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L,
>> 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3586L, 3588L, 3588L,
>> 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L,
>> 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L,
>> 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L,
>> 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L,
>> 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3588L, 3590L,
>> 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L,
>> 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L,
>> 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L,
>> 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L,
>> 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L, 3590L,
>> 3590L, 3590L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L,
>> 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L,
>> 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L,
>> 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L,
>> 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L,
>> 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3592L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L,
>> 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3594L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L, 3596L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L, 3598L,
>> 3598L, 3598L, 3598L, 3598L, 4416L, 4416L, 4416L, 4416L, 4416L,
>> 4414L, 4414L, 4414L, 4414L, 4414L, 4412L, 4412L, 4412L, 4412L,
>> 4412L, 4410L, 4410L, 4410L, 4410L, 4410L, 4408L, 4408L, 4408L,
>> 4408L, 4408L, 4406L, 4406L, 4406L, 4406L, 4406L, 4404L, 4404L,
>> 4404L, 4404L, 4404L, 4402L, 4402L, 4402L, 4402L, 4402L, 4400L,
>> 4400L, 4400L, 4400L, 4400L, 4398L, 4398L, 4398L, 4398L, 4398L,
>> 4396L, 4396L, 4396L, 4396L, 4396L, 4394L, 4394L, 4394L, 4394L,
>> 4394L, 4392L, 4392L, 4392L, 4392L, 4392L, 4390L, 4390L, 4390L,
>> 4390L, 4390L, 4388L, 4388L, 4388L, 4388L, 4388L, 4386L, 4386L,
>> 4386L, 4386L, 4386L, 4384L, 4384L, 4384L, 4384L, 4384L, 4382L,
>> 4382L, 4382L, 4382L, 4382L, 4380L, 4380L, 4380L, 4380L, 4380L,
>> 4378L, 4378L, 4378L, 4378L, 4378L, 4376L, 4376L, 4376L, 4376L,
>> 4376L, 4374L, 4374L, 4374L, 4374L, 4374L, 4372L, 4372L, 4372L,
>> 4372L, 4372L, 4370L, 4370L, 4370L, 4370L, 4370L, 4368L, 4368L,
>> 4368L, 4368L, 4368L, 4366L, 4366L, 4366L, 4366L, 4366L, 4364L,
>> 4364L, 4364L, 4364L, 4364L, 4362L, 4362L, 4362L, 4362L, 4362L,
>> 4362L, 4362L, 4360L, 4360L, 4360L, 4360L, 4360L, 4360L, 4360L,
>> 4358L, 4358L, 4358L, 4358L, 4358L, 4358L, 4358L, 4356L, 4356L,
>> 4356L, 4356L, 4356L, 4356L, 4356L, 4354L, 4354L, 4354L, 4354L,
>> 4354L, 4354L, 4354L, 4352L, 4352L, 4352L, 4352L, 4352L, 4352L,
>> 4352L, 4350L, 4350L, 4350L, 4350L, 4350L, 4350L, 4350L, 4348L,
>> 4348L, 4348L, 4348L, 4348L, 4348L, 4348L, 4346L, 4346L, 4346L,
>> 4346L, 4346L, 4346L, 4346L, 4344L, 4344L, 4344L, 4344L, 4344L,
>> 4344L, 4344L, 4342L, 4342L, 4342L, 4342L, 4342L, 4342L, 4342L,
>> 4340L, 4340L, 4340L, 4340L, 4340L, 4340L, 4340L, 4338L, 4338L,
>> 4338L, 4338L, 4338L, 4338L, 4338L, 4336L, 4336L, 4336L, 4336L,
>> 4336L, 4336L, 4336L, 4334L, 4334L, 4334L, 4334L, 4334L, 4334L,
>> 4334L, 4332L, 4332L, 4332L, 4332L, 4332L, 4332L, 4332L, 4330L,
>> 4330L, 4330L, 4330L, 4330L, 4330L, 4330L, 4328L, 4328L, 4328L,
>> 4328L, 4328L, 4328L, 4328L, 4326L, 4326L, 4326L, 4326L, 4326L,
>> 4326L, 4326L, 4324L, 4324L, 4324L, 4324L, 4324L, 4324L, 4324L,
>> 4322L, 4322L, 4322L, 4322L, 4322L, 4322L, 4322L, 4320L, 4320L,
>> 4320L, 4320L, 4320L, 4320L, 4320L, 4318L, 4318L, 4318L, 4318L,
>> 4318L, 4318L, 4318L, 4316L, 4316L, 4316L, 4316L, 4316L, 4316L,
>> 4316L, 4314L, 4314L, 4314L, 4314L, 4314L, 4314L, 4314L, 4312L,
>> 4312L, 4312L, 4312L, 4312L, 4312L, 4312L, 4310L, 4310L, 4310L,
>> 4310L, 4310L, 4310L, 4310L, 4308L, 4308L, 4308L, 4308L, 4308L,
>> 4308L, 4308L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L,
>> 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L,
>> 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4306L,
>> 4306L, 4306L, 4306L, 4306L, 4306L, 4306L, 4304L, 4304L, 4304L,
>> 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L,
>> 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L,
>> 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L, 4304L,
>> 4304L, 4304L, 4304L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L,
>> 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L,
>> 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L,
>> 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L, 4302L,
>> 4302L, 4302L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L,
>> 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L,
>> 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L,
>> 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L, 4300L,
>> 4300L, 4300L, 4300L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L,
>> 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4298L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L, 4296L,
>> 4296L, 4296L, 4296L, 4296L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L, 4294L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L, 4292L,
>> 4292L, 4292L, 4292L, 4292L, 4292L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L, 4290L,
>> 4290L, 4288L, 4288L, 4288L, 4288L, 4288L, 4288L, 4288L)), .Names = c("X",
>> "Y", "Z"), row.names = c(NA, 1500L), class = "data.frame")
>>
>>
>> On Wed, Dec 11, 2013 at 4:24 PM, Ben Bolker <bbolker at gmail.com> wrote:
>>
>>> On 13-12-11 11:13 AM, Shane Carey wrote:
>>>> Hey,
>>>>
>>>> Nope, it crashes. R stops abruptly. I sent just 500 points to show the
>>>> structure of the data. There is in fact over 22,000 points and it is 3-D
>>>> data composing of x,y,z data
>>>>
>>>> Cheers
>>>
>>>   We really, really need a reproducible example to be able to help
>>> further.  Since the 'akima.li' object is the one causing the trouble
>>> perhaps you could post the results of save() somewhere publicly
>>> accessible?
>>>
>>>   Do you have the same problem if you skip the rgl.spheres() and
>>> rgl.bbox() steps?
>>>
>>>   I'm able to plot a 100x100 random surface, even one with lots of NAs,
>>> without any trouble.  So there must be something special about akima.li
>>> , and guessing is going to be very hard ...
>>>
>>>   (At some point it might be useful to know your platform, version of R,
>>> etc., but that will only really become useful if it turns out that other
>>> people can plot your akima.li surface without any trouble ...)
>>>
>>>   Ben Bolker
>>>
>>>>
>>>>
>>>> On Wed, Dec 11, 2013 at 3:24 PM, Ben Bolker <bbolker at gmail.com> wrote:
>>>>
>>>>> Shane Carey <careyshan <at> gmail.com> writes:
>>>>>
>>>>>>
>>>>>> it just keeps crashing on me,
>>>>>>
>>>>>> It seems to crash on this line
>>>>>>
>>>>>> rgl.surface(akima.li$x,akima.li$y,akima.li
>>> $z,color="green",alpha=c(0.5))
>>>>>
>>>>>   When I run the code you gave us along with the data you gave us,
>>>>> I get
>>>>>
>>>>> Error in rgl.surface(akima.li$X, akima.li$Y, akima.li$Z, color =
>>> "green",
>>>>>  :
>>>>>   rows < 2
>>>>>
>>>>>   Maybe this is what you mean by "crashing"?  (On this list, "crash"
>>>>> is usually reserved to mean that your instance of R stops abruptly ...
>>>>> the behaviour above would be called an "error").
>>>>>
>>>>>   Looking at the plot generated from the data you sent, all of
>>>>> the (X,Y) points lie on a 1-dimensional curve in the (X,Y) plane --
>>>>> it's not surprising that it's going to be hard to generate an
>>>>> interpolated surface from these points.
>>>>>
>>>>>   Or maybe I'm misunderstanding something.
>>>>>
>>>>> 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))
>>>>> # interp surface:
>>>>> rgl.surface(akima.li$X,akima.li$Y,akima.li
>>> $Z,color="green",alpha=c(0.5))
>>>>>
>>>>> ______________________________________________
>>>>> R-help at 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.
>>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>> --
>> Shane
>>
> 
> 
>



More information about the R-help mailing list