Hi Sarah,
here is some of the data:
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), 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), 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)), .Names = c("X",
"Y", "Z"), row.names = c(NA, 500L), class = "data.frame")
and here is the code:
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:
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)
On Tue, Dec 10, 2013 at 3:48 PM, Sarah Goslee <sarah.gos...@gmail.com>wrote:
> Hi,
>
> On Tue, Dec 10, 2013 at 10:42 AM, Shane Carey <careys...@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 <careys...@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
>
--
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.
