我想将其他空间坐标合并到现有的空间线中。
我的数据:
观测值不一定与直线上的坐标匹配,因此我将观测值对齐到直线上。这为观察提供了新的坐标:它在先前存在的线上的位置。
我想将此新坐标合并到原始空间线数据中的正确位置,以免破坏沿线的原始坐标的正确顺序。
我研究了sp
对象的操作,但是无法弄清楚如何生成我想要的数据结构。
# Pre-existing line
> head(grpCoords)
LONG LAT
[1,] 140.6615 -8.846777
[2,] 140.5811 -8.728320
[3,] 140.4897 -8.620410
[4,] 140.1017 -8.300586
[5,] 140.0029 -8.195508
[6,] 139.9833 -8.166504
# Observation
> xy
[,1] [,2]
[1,] 167 -45
# New coordinate for observation on the line
> dist2Line(xy, grpCoords)
distance lon lat
[1,] 4123324 150.4824 -10.63691
例如:对于代表直线的XY坐标:
XY[1,1] _ _ _ _ _ _ _ XY[1,9] _ _ XY[1,12]
位置:
LONG LAT
1 1
1 9
1 12
如果观察到的坐标是XY [1,8],我希望:
XY[1,1] _ _ _ _ _ _ XY[1,8] XY[1,9] _ _ XY[1,12]
合并后的数据为:
LONG LAT
1 1
1 8
1 9
1 12