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

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

0 引言

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

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

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

2 總體思路、流程

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

3 采樣大格、小格編號(hào)

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

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

圖 2 采樣點(diǎn)編號(hào)思路流程Fig.2 Numbering of sampling points

3.1 大格編號(hào)

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

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

3.2 小格編號(hào)

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

4 采樣點(diǎn)編號(hào)

4.1 空間分析

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

4.2 采樣點(diǎn)編號(hào)

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

表 2 化探采樣點(diǎn)樣品自動(dòng)編號(hào)Table 2 Automatic numbering of geochemical exploration sample points

圖 3 小格編號(hào)Fig.3 Numbering of the small network

圖 4 小格編號(hào)屬性掛接Fig.4 Attributes articulation of the small network

圖 5 化探采樣點(diǎn)空間分析Fig.5 Spatial analysis of geochemical exploration sample points

5 結(jié)論

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

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

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

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