Hello, perhaps I am asking for something obvious, but I would like to perform interpolation of values, which are distributed on the line. my dataset consists of two vectors: stationing (distance from origin) and values. I would like to krige the values to obtain smooothed curve with variance estimates.
First thing I would thing of doing, is to transform stationing into two dimensional points (x=stationing, y=zeros) and to perform kriging as usual. However, I am not sure, if this procedure is correct, and I can not find a way to perform kriging of points on the line using packages gstat or rgeo, which I usually use. Any help or direction where to look would be appreciated. Below is example of my data Thanks, Zbynek > dput(linedata[1:50,]) structure(list(dist = c(0, 11.2800709216506, 22.0504005359196, 34.490657768126, 47.3296608447357, 62.2641846066267, 76.2888058138145, 92.6607330062524, 109.936149065093, 121.724278687636, 134.66298032256, 151.482314435605, 171.706062851761, 196.025602320903, 218.497807375147, 240.858487150145, 290.944014001843, 336.594862843896, 368.801228674272, 385.537118253217, 411.595514213105, 451.732280400594, 497.242828628585, 520.615676123105, 542.359870750547, 621.08089716662, 693.351499765505, 733.295961208736, 767.154638087923, 812.663318578618, 838.17880407268, 866.77880407254, 894.949888538801, 947.5878976196, 964.187897619693, 976.363693281579, 986.73480701161, 998.052492287232, 1027.60735652313, 1050.68570135112, 1063.82381384872, 1075.09501609169, 1086.18374911208, 1115.63353688717, 1139.86657219042, 1159.27301441961, 1179.15365827232, 1198.9213086133, 1225.76176315483, 1272.772400249), val = c(NA, NA, 0.908160651409972, 0.871501182963617, 0.918368640158739, 0.535216008816336, 0.949021011317396, 0.973433687159185, 0.976070428632721, 0.802086031449843, 0.961906084900447, 0.967220590046406, 0.97692016190632, 0.73837288344605, 0.803375801850693, 0.96357872833355, 0.944930467022756, 0.523655970604211, 0.00738153445264014, 0.0768875597281344, 0.445633626058785, 0.920272743496457, 0.486084005381625, 0.0267197415075827, 0.161488140793818, 0.533492819032941, 0.780736244753444, 0.117253233499661, 0.349786014519948, 0.853346494457695, 0.615928981211447, 0.38897304499115, 0.887739425294155, 0.748869234397016, 0.204990476432674, 0.963361133937722, 0.914181073306763, 0.878624741707031, 0.907694604916518, 0.95382032345456, 0.964164531499797, 0.962783600058143, 0.90797701844016, 0.890780300868731, 0.931906872733242, 0.856026387027286, 0.68907649322984, 0.638278881901186, 0.681647064000946, 0.649443287250305)), .Names = c("dist", "val"), row.names = c(NA, 50L), class = "data.frame") -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Kriging-of-values-distributed-on-line-tp7581559.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo