高維寅，李 晨

（甘肅省有色金屬地質(zhì)勘查局天水礦產(chǎn)勘查院，甘肅 天水 741024）

伴隨著國(guó)際經(jīng)濟(jì)的發(fā)展，不同領(lǐng)域的信息化建設(shè)也呈現(xiàn)著飛速發(fā)展現(xiàn)狀。必須建立一套可以定期更新的礦山大比例尺地形圖數(shù)據(jù)庫(kù)，實(shí)時(shí)監(jiān)測(cè)到礦山周圍的土地利用情況。因此，一些相關(guān)的測(cè)繪行業(yè)迫切需求具有高分辨率的遙感影像技術(shù)，對(duì)礦山大比例尺地形圖的精度分析要求也越來(lái)越高。

礦山大比例尺地形圖是礦山設(shè)計(jì)、礦石生產(chǎn)以及礦產(chǎn)資源管理的基礎(chǔ)資料之一，同時(shí)也在礦山的可持續(xù)發(fā)展中起著非常重要的作用[1]。隨著多金屬礦山的深入開(kāi)采以及礦區(qū)地形的拓展，礦區(qū)地表正在發(fā)生著變化，如何將礦區(qū)地表變化反映在礦山大比例尺地形圖中，是學(xué)術(shù)界一個(gè)非常重視的問(wèn)題。礦山大比例尺地形圖的更新通常是通過(guò)全站儀或航空像片的方式來(lái)采集數(shù)據(jù)，這兩種方法也存在各自的缺點(diǎn)，具有較長(zhǎng)的更新周期并且還會(huì)花費(fèi)大量的成本，在應(yīng)用過(guò)程中，很難滿足礦山生產(chǎn)的實(shí)際需求[2]。隨著RTK測(cè)量技術(shù)的發(fā)展，尤其是多平臺(tái)RTK測(cè)量技術(shù)的發(fā)展和應(yīng)用，為礦山大比例尺地形圖精度的分析提供了更有效的手段。

1 礦區(qū)概況與選擇多平臺(tái)RTK測(cè)量數(shù)據(jù)

1.1 分析礦區(qū)概況

選擇某一東西長(zhǎng)度為42km、南北長(zhǎng)度為86km的礦區(qū)，每一個(gè)礦的占地總面積大約為379.5km2，更新比例為1:2000的地形圖395幅。整個(gè)礦區(qū)的地勢(shì)比較平坦，但是由于礦區(qū)地貌變化緩慢，因此，本文以地物為主，礦山大比例尺地形圖中的地物主要包括植被、土壤質(zhì)量、交通附屬、居民地垣柵、采礦管線以及附屬設(shè)施設(shè)備等。

1.2 選擇多平臺(tái)RTK測(cè)量數(shù)據(jù)

多平臺(tái)RTK測(cè)量數(shù)據(jù)的選擇對(duì)礦山大比例尺地形圖精度分析有一定影響，因此，需要選擇合理的空間分辨率作為數(shù)據(jù)源[3]。

2 礦山大比例尺地形圖精度分析的實(shí)現(xiàn)

2.1 礦山大比例尺地形圖精度分析流程

利用多平臺(tái)RTK測(cè)量技術(shù)分析礦山大比例尺地形圖精度的技術(shù)流程如圖1所示。

圖1 礦山大比例尺地形圖精度分析流程

礦山大比例尺地形圖精度分析方法是將高分辨率的礦山大比例尺地形圖影像數(shù)據(jù)與數(shù)字化線圖之間進(jìn)行配準(zhǔn)，通過(guò)將兩幅地形圖相互疊加，識(shí)別礦山地形圖中地物的變化情況，并對(duì)變化過(guò)程中的地物進(jìn)行一些有效的數(shù)字化處理，更新礦山中原有的地物，實(shí)現(xiàn)礦山大比例尺地形圖的精度分析。

2.2 礦山大比例尺地形圖精度分析方法

礦山區(qū)域內(nèi)采集一套比例為1:2000的地形圖，經(jīng)過(guò)一定的掃描處理、數(shù)字化處理之后，形成一套數(shù)字線劃圖。將數(shù)字線劃圖與大比例尺地形圖進(jìn)行配準(zhǔn)時(shí)，將糾正之后的礦山大比例尺地形圖插入到繪圖軟件中，與此同時(shí)，還要在另外一個(gè)窗口中打開(kāi)數(shù)字線劃圖，利用繪圖軟件中按照原點(diǎn)粘貼的功能，將數(shù)字線劃圖粘貼到礦山大比例尺地形圖的影像圖窗口中，使兩者之間相互疊加對(duì)齊，實(shí)現(xiàn)兩個(gè)圖像的配準(zhǔn)。

作為數(shù)據(jù)源的大比例尺地形圖數(shù)據(jù)，數(shù)據(jù)質(zhì)量直接影響礦山大比例尺地形圖的精度，因此，需要選擇質(zhì)量高的大比例尺地形圖數(shù)據(jù)。

在預(yù)處理好大比例尺地形圖數(shù)據(jù)某一區(qū)域上均勻分布地選擇20個(gè)地面控制點(diǎn)，利用多平臺(tái)RTK測(cè)量技術(shù)進(jìn)行糾正。

當(dāng)?shù)V山大比例尺地形圖經(jīng)過(guò)糾正處理之后，開(kāi)始對(duì)礦山大比例尺地形圖進(jìn)行精度分析，精度分析的過(guò)程就是利用多平臺(tái)RTK測(cè)量技術(shù)得到的數(shù)據(jù)對(duì)原有的礦山大比例尺地形圖進(jìn)行修測(cè)和補(bǔ)測(cè)。當(dāng)?shù)V山大比例尺地形圖與數(shù)字線劃圖進(jìn)行準(zhǔn)確的配準(zhǔn)之后，由于礦山大比例尺地形圖的影像分辨率比較高，可以清晰地判讀地形圖中不同地物在變化過(guò)程中的位置，根據(jù)礦山區(qū)域地物的變化，采集地物特征點(diǎn)，實(shí)現(xiàn)礦區(qū)內(nèi)地物的修測(cè)和補(bǔ)測(cè)，并結(jié)合人機(jī)交互的方式，分析礦山大比例尺地形圖的精度。

2.3 礦山大比例尺地形圖精度分析結(jié)果

為了使礦山大比例地形圖精度分析具有科學(xué)性、合理性以及可靠性，在礦區(qū)內(nèi)選取了分布均勻的30個(gè)檢查點(diǎn)，用于1:2000的礦山大比例地形圖精度的外業(yè)檢核。通過(guò)對(duì)礦山區(qū)域?qū)嵉販y(cè)量得到的坐標(biāo)數(shù)據(jù)和大比例地形圖上同名點(diǎn)的坐標(biāo)數(shù)據(jù)進(jìn)行比較，經(jīng)過(guò)精確的分析得到了礦山大比例尺地形圖精度分析結(jié)果，如表1所示。

表1 礦山大比例尺地形圖精度分析結(jié)果

26 號(hào) -0.879 1.543 1.776 27 號(hào) 0.698 0.223 0.733 28 號(hào) -0.681 1.394 2.183 29 號(hào) -0.211 -0.107 0.237 30 號(hào) 1.685 -1.260 2.104

通過(guò)計(jì)算得到：

由于礦山大比例尺地形圖檢查點(diǎn)的坐標(biāo)誤差通常是由檢查點(diǎn)在礦山大比例尺地形圖上誤差m1和實(shí)地測(cè)量得到的同名點(diǎn)坐標(biāo)誤差m2兩個(gè)部分組成，因此，可以得到：

在計(jì)算精度上分析，可以認(rèn)為檢查點(diǎn)在礦山大比例尺地形圖上誤差m1和實(shí)地測(cè)量得到的同名點(diǎn)坐標(biāo)誤差m2是等精度的，因此，通過(guò)計(jì)算得到了：

由于礦山大比例尺地形圖上單點(diǎn)定位的絕對(duì)精度要求控制在1.0m以內(nèi)，因此，計(jì)算得到的礦山大比例尺地形圖精度已經(jīng)完全達(dá)到了精度要求。

綜合以上分析，為了提高礦山大比例尺地形圖的成圖質(zhì)量，礦山開(kāi)采和勘查的各部門(mén)必須采取比較嚴(yán)格的互相檢測(cè)以及地形圖精度分析措施，在分析過(guò)程中每一個(gè)精度分析環(huán)節(jié)必須緊密相連在一起，非常容易造成地形圖精度分析誤差的累積現(xiàn)象，因此，分析過(guò)程中的每一個(gè)關(guān)鍵環(huán)節(jié)都應(yīng)該嚴(yán)格執(zhí)行地形圖精度的分析，避免由于精度不足而出現(xiàn)沒(méi)有必要的經(jīng)濟(jì)損失。

采用基于多平臺(tái)RTK測(cè)量的礦山大比例尺地形圖精度分析方法分析了礦區(qū)內(nèi)大約400km2內(nèi)部1:2000的大比例尺地形圖，具有作業(yè)時(shí)間短且精度高的特點(diǎn)，可以完全滿足用戶對(duì)礦山大比例尺地形圖精度的需求，為礦山大比例尺地形圖精度分析提供了全新的經(jīng)驗(yàn)。

3 結(jié)語(yǔ)

本文提出了基于多平臺(tái)RTK測(cè)量的礦山大比例尺地形圖精度分析，在礦山大比例尺地形圖精度分析中發(fā)現(xiàn)，多平臺(tái)RTK測(cè)量中存在很多判讀準(zhǔn)確率比較低的地物，對(duì)于無(wú)法通過(guò)多平臺(tái)RTK測(cè)量判讀的地物需要再一次進(jìn)行實(shí)地調(diào)查，從而判斷地形圖的類型，因此，在分析礦山大比例尺地形圖精度分析中，要盡量加大調(diào)整多平臺(tái)RTK測(cè)量在其中的占比。在今后的研究中，還要將更多的礦山測(cè)繪技術(shù)應(yīng)用到礦山大比例尺地形圖精度分析中，進(jìn)一步提升礦山大比例尺地形圖精度分析性能。