王海鵬，孫海龍，魯紅峰，高永釗

遼寧省地質(zhì)勘查院，大連 金州 116100

0 引言

區(qū)域化探工作過程中涉及到大量的化探采樣點的編號工作，如果是規(guī)則網(wǎng)化探工作，可根據(jù)規(guī)律性的采樣點分布進行批量采樣點編號。但在實際工作過程中無論是土壤測量還是水系沉積物測量常常布設(shè)不規(guī)則網(wǎng)采樣點，這樣數(shù)千甚至上萬個采樣點的編號工作就耗費了大量的人力和時間，而且容易出錯。筆者在實際工作過程中在MAPGIS條件下總結(jié)一套較好的方法，將采樣點的樣號屬性掛接到采樣點的MAPGIS點文件中，從而完美實現(xiàn)不規(guī)則網(wǎng)化探采樣點自動編號工作，并為后續(xù)的采樣點的處理工作帶來很大的方便。

1 化探采樣點布設(shè)規(guī)則

水系或土壤測量一般可采用地形圖定點。在工作的范圍內(nèi)劃出長寬各為0.5 km的方格網(wǎng)，以四個方格(1 km2)作為采樣大格。大格的編號順序自左而右再自上而下。每個大格中有四個面積為0.25km2的小格，編號順序自左而右自上而下標號a，b，c，d。在每一小格中采集的第一號樣品為1，第二號樣品標號為2，以此類推。每個采樣點根據(jù)其所處的位置按上述順序進行編號（圖1）。按此原則區(qū)域性采樣點可預先設(shè)計并標繪在地形圖上。但在實際工作中我們常常按地形、水系及地質(zhì)等綜合條件進行采樣點布設(shè)，采樣點可能數(shù)千甚至上萬個，布設(shè)完成之后的編號工作量非常大，容易出錯。

2 總體思路、流程

化探采樣點編號涉及到大格區(qū)編號、小格區(qū)編號及采樣點本身編號（圖1），如：002（大格區(qū)編號）a（小格區(qū)編號）1（采樣點編號）。本文論及的不規(guī)則網(wǎng)化探采樣點自動編號總體思路、流程見圖2。

3 采樣大格、小格編號

按照化探采樣點布設(shè)規(guī)則，巧妙地利用點位坐標（X、Y）的升降排序來實現(xiàn)大格區(qū)、小格區(qū)的自動編號。

圖 1 化探采樣點布設(shè)Fig.1 Geochemical sampling points

圖 2 采樣點編號思路流程Fig.2 Numbering of sampling points

3.1 大格編號

首先在MAPGIS條件下生成“大格區(qū).wp”文件的lable點文件，其次將lable點坐標位置賦屬性到lable點文件內(nèi)[1]，輸出lable點屬性到EXCELE中，然后按照X升序、Y降序原則進行重新排序，同時增加大格編號字段，按照重新排序后的次序批量對大格區(qū)進行編號（表1）。接下來以ID號為索引字段，利用MAPGIS屬性掛接[1]功能將EXCELE中大格編號字段掛接到“大格區(qū).wp”文件中的“大格編號”屬性字段內(nèi)，利用MAPGIS中根據(jù)屬性標注示的功能將“大格區(qū).wp”文件中大格編號屬性字段自動標注到采樣點布設(shè)圖上，這樣就實現(xiàn)了大格編號的自動化標注。

表 1 大格區(qū)編號Table 1 Numbering of the big network

3.2 小格編號

首先在MAPGIS采樣點布設(shè)文件中將第1個大格內(nèi)的4個小格標注a，b，c，d，按照小格固定的行列間距進行陣列復制點最后使得所有的小格內(nèi)具有相應的標注（圖3），存儲點文件為“小格標注.wt”。利用MAPGIS點文件注釋賦屬性功能將小格區(qū)內(nèi)的點注釋賦為相應點的屬性（圖4），再利用MAPGIS中l(wèi)able點合并功能將“小格標注.wt”中“abcd”的標注屬性合并到“小格區(qū).wp”文件“小格編號”屬性字段內(nèi)（圖4），這樣就實現(xiàn)了小格編號的自動化。

4 采樣點編號

4.1 空間分析

在MAPGIS“空間分析”模塊下，將“化探采樣.wt”文件分別與“大格區(qū).wp”、“小格區(qū).wp”進行交集的空間分析[2]，將“大格編號”、“小格編號”屬性字段合并到“化探采樣.wt”文件內(nèi)（圖5）。然后將采樣點坐標位置賦到“化探采樣.wt”文件屬性內(nèi)。

4.2 采樣點編號

輸出“化探采樣.wt”屬性到EXCELE中，以大格編號為第1升序，小格編號為第2升序、X坐標為第3升序、Y坐標為第4降序原則進行重新排序（表2），同時增加“過渡編號”、“樣品編號”、“圖面樣品編號”“最終樣品編號”4個屬性字段，然后按照重新排序后的次序依次對化探采樣點進行編號。下面以第1行為例進行具體的剖析：“過渡編號”內(nèi)公式[3]為“=D2&E2”，“樣品編號”內(nèi)公式為“=IF(F1=F2,1+G1,1)”，“圖面樣品編號”內(nèi)公式為“=E2&G2”“最終樣品編號”內(nèi)公式為“=F2&G2”，利用EXCELE填充功能將第1行的公式內(nèi)容向下充填到最后一行完成相應的計算。然后以ID號為索引字段，利用MAPGIS屬性掛接功能將EXCELE中“圖面樣品編號”及“最終樣品編號”掛接到“化探采樣.wt”文件內(nèi)，這樣就完成了所有化探采樣點的自動編號工作，我們可以根據(jù)不同的需求把不同的編號標注到采樣點布設(shè)圖上。由于“化探采樣.wt”文件內(nèi)屬性字段內(nèi)容齊全，可以根據(jù)不同的需要進行相應的處理工作，方便快捷。

表 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

5 結(jié)論

本文討論MAPGIS條件下利用lable點的生成及合并、點位置轉(zhuǎn)屬性、空間分析、屬性掛接等功能，在充分理解化探采樣點布設(shè)規(guī)則的前提下，對“化探采樣.wt”文件中X、Y坐標進行不同級別的排序，結(jié)合EXCELE相關(guān)公式的編制來完美地實現(xiàn)化探采樣點的自動編號工作，方法實用高效。對于其它類似的工作同樣具有很強的適用性。

[1] 吳信才，MAPGIS地理信息系統(tǒng)[M].北京：電子工業(yè)出版社，2004.

[2] 劉俊長，龔紅雷，張玉領(lǐng)，等.用MAPGIS實現(xiàn)區(qū)域化探數(shù)據(jù)的空間分析[J]物探與化探.2008.32（6）：690.

[3] 張 強，劉 飚，Excele2007與VBA編程從入門到精通[M].北京：清華大學出版社.2010.