一、目标
一个点集里的点两两连线,比如有4个点,则连3+2+1=6条线;有38个点,则连37+36+35+···+3+2+1=703条线,其原始坐标数据及最终样式类似于以下这种:
→
图一:原始数据及最终效果
二、总体思路
使用arcgis中的XY To Line工具,如下图所
其中,Start X Field、Start Y Field、End X Field、End Y Field分别为起点和终点的坐标值,Line Type为要构建的大地测量线的类型,包括GEODESIC(默认)、GREAT_CIRCLE、RHUMB_LINE、NORMAL_SECTION这四种构造方法,ID为输入表中的字段,此字段和坐标值包含在输出中,可用于将输出要素与输入表中的记录相连接。Spatial_Reference为空间参考。
可以看到,要准备的数据是包含起点和终点坐标值、ID的一张表。
三、表数据准备
(一)坐标值数据准备
首先,是要将Start X、Start Y合并,为了便于区分,将StationName也合并到一个单元格。即在XYOld列中,即D2单元格输入“=A2&","&B2&","&C2”,并下拉至列尾;
其次,将D列复制到E列(XYOldCopy);
(二)坐标值排列组合
对D列和E列进行排列组合。右键左下方sheet1名称(本例为PointstoTest)→查看代码→在弹出的对话框中粘贴如下代码:
1 Sub两列组合()2 'EXCEL880.COM产品,2017-7-12
3 Dimn1, n2, arr, brr, i, j, k4 n1 =[D65536].End(xlUp).Row5 n2 =[E65536].End(xlUp).Row6 arr = Range("D2:D" &n1)7 brr = Range("E2:E" &n2)8 ReDim crr(1 To n1 * n2, 0)9 Dimdic10 Set dic = CreateObject("scripting.dictionary")11 For i = 1 To UBound(arr)12 For j = i + 1 To UBound(brr)13 Key = arr(i, 1) & "," & brr(j, 1)14 'If Not dic.Exists(Key) And Not dic.Exists(brr(i, 1) & "," & arr(j, 1)) Then
15 '字典检测是否双向重复
16 dic.Add Key, ""
17 k = k + 1
18 crr(k, 0) =Key19 'End If
20 Next
21 Next
22 [F2:F65536].Clear23 [F2].Resize(k) =crr24 End Sub
点击弹出框上方的“运行”→“运行子过程/用户窗体”。截图如下
对F列以","逗号为分隔符进行分列并将最终结果粘贴至新的表格,添加OIDMINE列用于XY To Line工具的ID输入字段,结果如下图所示:
四、在arcgis中使用XY To Line工具运行即可。
在arcgis中,ArcToolbox→Data Management Tools→Features→XY To Line,并输入相应参数即可。