王海鹏,孙海龙,鲁红峰,高永钊
辽宁省地质勘查院,大连 金州 116100
区域化探工作过程中涉及到大量的化探采样点的编号工作,如果是规则网化探工作,可根据规律性的采样点分布进行批量采样点编号。但在实际工作过程中无论是土壤测量还是水系沉积物测量常常布设不规则网采样点,这样数千甚至上万个采样点的编号工作就耗费了大量的人力和时间,而且容易出错。笔者在实际工作过程中在MAPGIS条件下总结一套较好的方法,将采样点的样号属性挂接到采样点的MAPGIS点文件中,从而完美实现不规则网化探采样点自动编号工作,并为后续的采样点的处理工作带来很大的方便。
水系或土壤测量一般可采用地形图定点。在工作的范围内划出长宽各为0.5 km的方格网,以四个方格(1 km2)作为采样大格。大格的编号顺序自左而右再自上而下。每个大格中有四个面积为0.25km2的小格,编号顺序自左而右自上而下标号a,b,c,d。在每一小格中采集的第一号样品为1,第二号样品标号为2,以此类推。每个采样点根据其所处的位置按上述顺序进行编号(图1)。按此原则区域性采样点可预先设计并标绘在地形图上。但在实际工作中我们常常按地形、水系及地质等综合条件进行采样点布设,采样点可能数千甚至上万个,布设完成之后的编号工作量非常大,容易出错。
化探采样点编号涉及到大格区编号、小格区编号及采样点本身编号(图1),如:002(大格区编号)a(小格区编号)1(采样点编号)。本文论及的不规则网化探采样点自动编号总体思路、流程见图2。
按照化探采样点布设规则,巧妙地利用点位坐标(X、Y)的升降排序来实现大格区、小格区的自动编号。
图 1 化探采样点布设Fig.1 Geochemical sampling points
图 2 采样点编号思路流程Fig.2 Numbering of sampling points
首先在MAPGIS条件下生成“大格区.wp”文件的lable点文件,其次将lable点坐标位置赋属性到lable点文件内[1],输出lable点属性到EXCELE中,然后按照X升序、Y降序原则进行重新排序,同时增加大格编号字段,按照重新排序后的次序批量对大格区进行编号(表1)。接下来以ID号为索引字段,利用MAPGIS属性挂接[1]功能将EXCELE中大格编号字段挂接到“大格区.wp”文件中的“大格编号”属性字段内,利用MAPGIS中根据属性标注示的功能将“大格区.wp”文件中大格编号属性字段自动标注到采样点布设图上,这样就实现了大格编号的自动化标注。
表 1 大格区编号Table 1 Numbering of the big network
首先在MAPGIS采样点布设文件中将第1个大格内的4个小格标注a,b,c,d,按照小格固定的行列间距进行阵列复制点最后使得所有的小格内具有相应的标注(图3),存储点文件为“小格标注.wt”。利用MAPGIS点文件注释赋属性功能将小格区内的点注释赋为相应点的属性(图4),再利用MAPGIS中lable点合并功能将“小格标注.wt”中“abcd”的标注属性合并到“小格区.wp”文件“小格编号”属性字段内(图4),这样就实现了小格编号的自动化。
在MAPGIS“空间分析”模块下,将“化探采样.wt”文件分别与“大格区.wp”、“小格区.wp”进行交集的空间分析[2],将“大格编号”、“小格编号”属性字段合并到“化探采样.wt”文件内(图5)。然后将采样点坐标位置赋到“化探采样.wt”文件属性内。
输出“化探采样.wt”属性到EXCELE中,以大格编号为第1升序,小格编号为第2升序、X坐标为第3升序、Y坐标为第4降序原则进行重新排序(表2),同时增加“过渡编号”、“样品编号”、“图面样品编号”“最终样品编号”4个属性字段,然后按照重新排序后的次序依次对化探采样点进行编号。下面以第1行为例进行具体的剖析:“过渡编号”内公式[3]为“=D2&E2”,“样品编号”内公式为“=IF(F1=F2,1+G1,1)”,“图面样品编号”内公式为“=E2&G2”“最终样品编号”内公式为“=F2&G2”,利用EXCELE填充功能将第1行的公式内容向下充填到最后一行完成相应的计算。然后以ID号为索引字段,利用MAPGIS属性挂接功能将EXCELE中“图面样品编号”及“最终样品编号”挂接到“化探采样.wt”文件内,这样就完成了所有化探采样点的自动编号工作,我们可以根据不同的需求把不同的编号标注到采样点布设图上。由于“化探采样.wt”文件内属性字段内容齐全,可以根据不同的需要进行相应的处理工作,方便快捷。
表 2 化探采样点样品自动编号Table 2 Automatic numbering of geochemical exploration sample points
图 3 小格编号Fig.3 Numbering of the small network
图 4 小格编号属性挂接Fig.4 Attributes articulation of the small network
图 5 化探采样点空间分析Fig.5 Spatial analysis of geochemical exploration sample points
本文讨论MAPGIS条件下利用lable点的生成及合并、点位置转属性、空间分析、属性挂接等功能,在充分理解化探采样点布设规则的前提下,对“化探采样.wt”文件中X、Y坐标进行不同级别的排序,结合EXCELE相关公式的编制来完美地实现化探采样点的自动编号工作,方法实用高效。对于其它类似的工作同样具有很强的适用性。
[1] 吴信才,MAPGIS地理信息系统[M].北京:电子工业出版社,2004.
[2] 刘俊长,龚红雷,张玉领,等.用MAPGIS实现区域化探数据的空间分析[J]物探与化探.2008.32(6):690.
[3] 张 强,刘 飚,Excele2007与VBA编程从入门到精通[M].北京:清华大学出版社.2010.