[R-sig-Geo] Kriging of values distributed on line

[zbynek.janoska at centrum.cz]

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.