王炳新

(遼寧省有色地質(zhì)一〇一隊(duì)有限責(zé)任公司，遼寧 撫順 113006)

為了創(chuàng)造我國(guó)良好的用地環(huán)境，推動(dòng)綠色發(fā)展和生態(tài)文明建設(shè)與環(huán)境保護(hù)，遏制違法用地的蔓延趨勢(shì)，我國(guó)《土地管理法》和《城市規(guī)劃法》對(duì)于土地使用有著明確的要求[1]，對(duì)于擅自開(kāi)發(fā)、使用、倒賣(mài)等非法用地要根據(jù)違法用地土方量進(jìn)行處罰，故準(zhǔn)確的違法用地土方量成為了公正處罰的基礎(chǔ)。

就土方量計(jì)算而言，對(duì)于大面積違法用地，若采用傳統(tǒng)手段測(cè)量土方量會(huì)耗費(fèi)大量人力和時(shí)間，測(cè)量效率較低導(dǎo)致項(xiàng)目周期長(zhǎng)、難度大，現(xiàn)研究采用無(wú)人機(jī)攝影測(cè)量技術(shù)對(duì)違法用地進(jìn)行數(shù)據(jù)采集，計(jì)算違法用地土方量，提高土方測(cè)量作業(yè)效率和精度[2]。

國(guó)內(nèi)外已有學(xué)者研究采用無(wú)人機(jī)技術(shù)計(jì)算土方量。李衛(wèi)軍等[3]人采用無(wú)人機(jī)計(jì)算礦山開(kāi)采的土方量采剝量，避免了人員下井的安全風(fēng)險(xiǎn)問(wèn)題，同時(shí)對(duì)礦山開(kāi)采生產(chǎn)成本監(jiān)控有著重要的作用。劉兆慧等[4]人采用無(wú)人機(jī)傾斜攝影測(cè)量作業(yè)的方式，通過(guò)航測(cè)影像獲取測(cè)區(qū)三維點(diǎn)云，計(jì)算測(cè)區(qū)實(shí)際土方量。李東升等[5]人采用無(wú)人機(jī)獲取露天測(cè)區(qū)影像，根據(jù)計(jì)算軟件生成三維點(diǎn)云數(shù)據(jù)，制作測(cè)區(qū)的DEM和DOM等數(shù)字產(chǎn)品，基于前期的數(shù)字產(chǎn)品計(jì)算露天挖填土方量。

1 傳統(tǒng)與新方法對(duì)比

隨著無(wú)人機(jī)技術(shù)的飛速發(fā)展和應(yīng)用普及，它在日常測(cè)繪工作中的作用越來(lái)越重要，有效推動(dòng)了測(cè)繪技術(shù)的快速發(fā)展，與常規(guī)GNSS實(shí)地測(cè)繪相比，在作業(yè)效率、成本、安全性等方面有著顯著優(yōu)勢(shì)。下面從土方量數(shù)據(jù)采集、計(jì)算數(shù)據(jù)方法等方面出發(fā)分析測(cè)繪方式的變革情況[6]，比較傳統(tǒng)測(cè)量技術(shù)和無(wú)人機(jī)技術(shù)在各個(gè)方面的優(yōu)劣勢(shì)。

(1)數(shù)據(jù)采集方面：傳統(tǒng)測(cè)量方式采用全站儀或人工手持RTK進(jìn)行實(shí)地?cái)?shù)據(jù)采集，此方法精度高但工作量大、耗時(shí)長(zhǎng)，故投入成本較高。無(wú)人機(jī)技術(shù)在效率方面有了質(zhì)的飛躍，無(wú)接觸式的面采集方式提升了數(shù)據(jù)采集的效率，此外，大大減小了人員投入和人工勞動(dòng)強(qiáng)度，但無(wú)人機(jī)作業(yè)受天氣影響較大[7]。

(2)數(shù)據(jù)處理方面：傳統(tǒng)手段多采用插點(diǎn)式三角網(wǎng)計(jì)算土方量，方法精度可靠，但整個(gè)過(guò)程需要人工參與。無(wú)人機(jī)技術(shù)基于航測(cè)內(nèi)業(yè)處理軟件，通過(guò)專(zhuān)業(yè)數(shù)據(jù)處理軟件獲取照片的位置、姿態(tài)、相機(jī)內(nèi)方位元素等信息，采用影像像對(duì)匹配生成地物三維點(diǎn)云，基于點(diǎn)云數(shù)據(jù)生成數(shù)字高程模型DEM和正射影像圖DOM，計(jì)算出勾選區(qū)域的土方量，此方法精度可靠，且自動(dòng)化程度高。

傳統(tǒng)手段采用全站儀或者GNSS-RTK是“點(diǎn)”數(shù)據(jù)采集方式，對(duì)堆放物實(shí)地進(jìn)行三維數(shù)據(jù)采集的密度遠(yuǎn)遠(yuǎn)沒(méi)有無(wú)人機(jī)數(shù)據(jù)處理后生成的三維點(diǎn)云的密度大[8]，所以在計(jì)算土石方量時(shí)，航測(cè)計(jì)算的土石方量要比傳統(tǒng)測(cè)量手段計(jì)算結(jié)果精度更高，如表1所示。

表1 傳統(tǒng)方式與航測(cè)方式效果對(duì)比表

由上述分析可知，與傳統(tǒng)測(cè)量手段相比，無(wú)人機(jī)測(cè)量技術(shù)在土方量計(jì)算上具有一定的優(yōu)越性。

2 無(wú)人機(jī)航測(cè)

無(wú)人機(jī)航測(cè)流程主要涵蓋幾個(gè)方面，首先，要制定本次航測(cè)項(xiàng)目方案、布設(shè)像控點(diǎn)、航測(cè)正式作業(yè)、成果質(zhì)量檢查、內(nèi)業(yè)影像處理、空三計(jì)算、生成DEM和DOM等成果，再進(jìn)行土方量計(jì)算[9]。無(wú)人機(jī)航測(cè)流程圖如圖1所示。

圖1 無(wú)人機(jī)航測(cè)流程圖

2.1 測(cè)區(qū)情況

本次測(cè)區(qū)違法用地占地面積約6 000 m2，測(cè)區(qū)為空曠的鹽堿地，周?chē)胁糠謽?shù)林等地物遮擋，測(cè)區(qū)主體坡勢(shì)平緩。

2.2 無(wú)人機(jī)設(shè)備

本項(xiàng)目采用大疆精靈4RTK無(wú)人機(jī)，精靈4RTK差分模塊優(yōu)異的定位精度保證了無(wú)人機(jī)在航測(cè)時(shí)數(shù)據(jù)的準(zhǔn)確性,其差分模塊的參數(shù)如表2所示。

表2 差分模塊參數(shù)表

2.3 像控點(diǎn)數(shù)學(xué)基礎(chǔ)

根據(jù)違法用地情況，在測(cè)區(qū)范圍內(nèi)均勻布置了6個(gè)像控點(diǎn)，點(diǎn)號(hào)為X1～X6，標(biāo)志均清晰可見(jiàn)，布設(shè)后對(duì)其采用RKT測(cè)量像控點(diǎn)三維坐標(biāo)。

坐標(biāo)系統(tǒng)方面，像控點(diǎn)平面采用2000國(guó)家大地坐標(biāo)系，高程采用1985國(guó)家高程基準(zhǔn)，投影為高斯克呂格3°帶投影，中央子午線為120°。

2.4 外業(yè)航測(cè)作業(yè)

首先，根據(jù)違法用地區(qū)范圍進(jìn)行航線規(guī)劃，本次采用“井”字航線飛行，可以保證更多的同名點(diǎn)，進(jìn)而提高三維點(diǎn)云的密度和精確度。同時(shí)對(duì)航高、航速、拍攝模式、相片重疊度、測(cè)光模式、云臺(tái)角度等參數(shù)進(jìn)行設(shè)置，具體參數(shù)設(shè)置情況如圖2所示，完成航線規(guī)劃與參數(shù)設(shè)置后，一鍵起飛無(wú)人機(jī)。飛行過(guò)程中關(guān)注外界環(huán)境風(fēng)力、信號(hào)強(qiáng)度、電池電量等因素，如遇緊急情況應(yīng)立即返航。

圖2 參數(shù)設(shè)置

本次航飛累計(jì)飛行2個(gè)架次，飛行時(shí)間35 min，外業(yè)踏勘像控點(diǎn)布置與采集30 min，外業(yè)無(wú)人機(jī)航飛共采集影像數(shù)據(jù)359張。

2.5 數(shù)據(jù)處理

數(shù)據(jù)處理環(huán)節(jié)主要包括影像檢查、影像預(yù)處理、空三計(jì)算及成果生成與輸出。

首先，將外業(yè)采集的攜帶POS數(shù)據(jù)的照片導(dǎo)入到ContextCapture軟件中，檢查影像是否包含全部測(cè)區(qū)范圍，若出現(xiàn)遺漏則需要補(bǔ)測(cè)或重測(cè)相應(yīng)區(qū)域數(shù)據(jù)。完成影像數(shù)據(jù)檢查后，設(shè)置軟件自動(dòng)化處理參數(shù)，影像進(jìn)入預(yù)處理環(huán)節(jié)，預(yù)處理環(huán)節(jié)會(huì)完成影像的內(nèi)定向和外定向，再進(jìn)行空中三角測(cè)量[10]，ContextCapture軟件可自動(dòng)讀取影像攜帶的POS數(shù)據(jù)，然后，直接進(jìn)行空中三角測(cè)量計(jì)算，如圖3所示，大大提升了航測(cè)數(shù)據(jù)處理的自動(dòng)化程度，空中三角測(cè)量完成后，輸出測(cè)區(qū)的DEM和DOM，根據(jù)拼接好的正射影像圖檢查測(cè)區(qū)變形情況，若出現(xiàn)變形較大的情況應(yīng)及時(shí)檢查原始數(shù)據(jù)、處理過(guò)程等是否出問(wèn)題。最后，根據(jù)相鄰的重疊影像匹配測(cè)區(qū)的同名點(diǎn)，生成測(cè)區(qū)密集的三維點(diǎn)云坐標(biāo)，實(shí)現(xiàn)實(shí)景三維模型制作[11]。

圖3 空中三角計(jì)算

2.6 土方量計(jì)算

為準(zhǔn)確計(jì)算違法用地土方量，分別生產(chǎn)像控刺點(diǎn)和不刺點(diǎn)的數(shù)據(jù)格式為.osgb和.3mx模型產(chǎn)品，具體計(jì)算分析過(guò)程如下：

(1)對(duì)有植被覆蓋或樹(shù)木遮擋的沙堆等不規(guī)則堆積物進(jìn)行土方量計(jì)算時(shí)：將帶有像控的.osgb模型加載到清華山維的EPS三維測(cè)圖平臺(tái)，同時(shí)，將外業(yè)植被遮擋區(qū)域用GNSS-RTK采集到的數(shù)據(jù)導(dǎo)入到模型中，然后，對(duì)不規(guī)則沙堆在模型中裸露明顯地形特征點(diǎn)進(jìn)行三維坐標(biāo)及高程數(shù)據(jù)采集，利用采集的高程信息首先生成起算基準(zhǔn)面的三角網(wǎng)數(shù)據(jù)文件與地表三角網(wǎng)TIN數(shù)據(jù)文件[12]，通過(guò)土方量計(jì)算命令對(duì)不規(guī)則的沙堆進(jìn)行土石方量計(jì)算，如圖4所示。

圖4 三角網(wǎng)模型TIN

(2)對(duì)完全裸露的沙堆等不規(guī)則堆積物進(jìn)行土石方量計(jì)算時(shí)：將已加像控的.3mx數(shù)據(jù)模型加載到ContextCapture viewer軟件中，利用體積量算工具，對(duì)違法用地測(cè)區(qū)范圍進(jìn)行手動(dòng)勾繪，計(jì)算出占地面積與體積，同時(shí)對(duì)幾個(gè)明顯地物點(diǎn)進(jìn)行三維坐標(biāo)記錄。

通過(guò)ContextCapture viewer軟件，手動(dòng)勾劃其中一塊違法占用鹽堿地區(qū)域，其面積為516.530 m2，體積為880.854 m3，如圖5所示。照此步驟采集整個(gè)測(cè)區(qū)的總面積和土方量，采樣結(jié)果如表3所示。

圖5 土方量計(jì)算

表3 測(cè)區(qū)總采樣面積和土方量

3 精度驗(yàn)證

任何一個(gè)測(cè)繪項(xiàng)目成果需要經(jīng)過(guò)多級(jí)檢查，各項(xiàng)指標(biāo)、精度檢查滿足規(guī)范要求后，才能保證是準(zhǔn)確、合格的測(cè)繪成果，以提供給其他部門(mén)應(yīng)用。故需要對(duì)本次測(cè)繪的成果進(jìn)行檢查，檢查內(nèi)容主要包括：檢校點(diǎn)中誤差統(tǒng)計(jì)分析、成果質(zhì)量分析、坐標(biāo)較差分析等幾個(gè)環(huán)節(jié)。

3.1 檢校點(diǎn)中誤差分析

本次共選取了6個(gè)檢校點(diǎn)作為項(xiàng)目精度檢查點(diǎn)，點(diǎn)名分別列為J1～J6，這6個(gè)檢校點(diǎn)均勻分布整個(gè)測(cè)區(qū)且布設(shè)于特征地物上，便于內(nèi)業(yè)判讀與測(cè)量。

從外業(yè)GNSS-RTK測(cè)量手簿中導(dǎo)出檢校點(diǎn)坐標(biāo)和誤差，分析本次項(xiàng)目檢校點(diǎn)的平面和高程中誤差，以檢查檢校點(diǎn)的測(cè)量精度情況，檢校點(diǎn)的中誤差統(tǒng)計(jì)結(jié)果如表4所示。

由表4可知，檢校點(diǎn)的平面精度、高程精度和總中誤差均符合國(guó)家規(guī)范《1:500 1:1 000 1:2 000地形圖航空攝影測(cè)量?jī)?nèi)業(yè)規(guī)范》(GB/T 7930-2008)中相關(guān)要求，故本次項(xiàng)目檢校點(diǎn)達(dá)到了項(xiàng)目精度要求，可以作為本次項(xiàng)目精度的判斷標(biāo)準(zhǔn)。

表4 檢校點(diǎn)中誤差統(tǒng)計(jì)表/m

3.2 成果質(zhì)量分析

測(cè)區(qū)成果輸出后，需要對(duì)本次成果進(jìn)行精度檢查，檢查標(biāo)準(zhǔn)應(yīng)符合本次項(xiàng)目精度要求，對(duì)測(cè)區(qū)的平面精度和高程精度應(yīng)該按照《1:500 1:1 000 1:2 000地形圖航空攝影測(cè)量?jī)?nèi)業(yè)規(guī)范》(GB/T 7930-2008)中關(guān)于1：500地形測(cè)繪相關(guān)規(guī)定進(jìn)行精度檢查[13]。

通過(guò)軟件生成的空中三角計(jì)算報(bào)告得到像控點(diǎn)的平差精度，精度滿足限差要求后，對(duì)明顯地物點(diǎn)進(jìn)行檢核，檢核平面點(diǎn)330多個(gè)，最大平面誤差為0.030 m，檢核高程點(diǎn)270多個(gè)，最大高程差為0.027 m，均滿足1:500地形圖精度要求[14]。

3.3 檢校點(diǎn)坐標(biāo)較差

航測(cè)精度驗(yàn)證方式是根據(jù)檢校點(diǎn)的坐標(biāo)位置，從航測(cè)中選取相同的點(diǎn)位坐標(biāo)，將兩種坐標(biāo)成果在同坐標(biāo)系下進(jìn)行坐標(biāo)較差比較，比較結(jié)果如表5所示。

由表5可知，本次航測(cè)坐標(biāo)點(diǎn)和實(shí)地檢查坐標(biāo)點(diǎn)較差差值均在3 cm以下，滿足《低空數(shù)字航空攝影規(guī)范》(CH/Z 3005-2010)中對(duì)于航測(cè)精度的要求[15]。

表5 檢校點(diǎn)與航測(cè)點(diǎn)精度對(duì)比表

4 結(jié) 語(yǔ)

本文基于無(wú)人機(jī)攝影測(cè)量技術(shù)計(jì)算違法用地土方量，分析了傳統(tǒng)手段和無(wú)人機(jī)新測(cè)繪方式在數(shù)據(jù)采集和處理等方面的優(yōu)缺點(diǎn)。研究了航測(cè)內(nèi)外業(yè)整個(gè)流程，從數(shù)據(jù)采集、數(shù)據(jù)處理分析再到土方量計(jì)算成套的技術(shù)體系。驗(yàn)證分析了無(wú)人機(jī)技術(shù)在土石方量計(jì)算中的可行性，減小了人工勞動(dòng)強(qiáng)度，大大提高了生產(chǎn)效率，最后，通過(guò)精度驗(yàn)證分析，證明了無(wú)人機(jī)技術(shù)成果的精度符合相關(guān)規(guī)范要求，滿足項(xiàng)目生產(chǎn)精度，可以作為一種可靠的測(cè)繪成果為執(zhí)法部門(mén)提供基礎(chǔ)數(shù)據(jù)，該方法經(jīng)過(guò)多此理論推敲和實(shí)踐驗(yàn)證，值得應(yīng)用推廣。