韓麗蘋(píng)，黃 榮

安徽省地質(zhì)礦產(chǎn)勘查局325 地質(zhì)隊(duì)，安徽 淮北 235000

目前，礦山治理監(jiān)測(cè)最常使用的是衛(wèi)星遙感技術(shù)，它的優(yōu)勢(shì)在于覆蓋面積大，成本低，但時(shí)效性、靈活性、分辨率偏低，主要適用于調(diào)查全域礦產(chǎn)開(kāi)采秩序問(wèn)題[1]。無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)是一項(xiàng)近年來(lái)不斷發(fā)展的高新技術(shù)，它通過(guò)無(wú)人機(jī)空中作業(yè)來(lái)獲得高分辨率、高重疊的影像數(shù)據(jù)，快速地建立帶有高精度坐標(biāo)的三維模型，并以模型為基礎(chǔ)矢量化輸出DSM、DOM、TDOM、DLG 等多種測(cè)繪數(shù)字產(chǎn)品[2]。無(wú)人機(jī)傾斜攝影技術(shù)機(jī)動(dòng)性強(qiáng)，分辨率高，面積小，經(jīng)費(fèi)投入低，能夠彌補(bǔ)衛(wèi)星遙感技術(shù)的不足，直觀地展現(xiàn)整個(gè)礦山治理區(qū)的范圍和治理狀況。露天礦山地形復(fù)雜多變，無(wú)人機(jī)傾斜攝影技術(shù)可以實(shí)時(shí)動(dòng)態(tài)對(duì)礦山治理進(jìn)行監(jiān)測(cè)，提供真實(shí)的三維模型，直觀反映礦山開(kāi)發(fā)利用現(xiàn)狀和周邊地貌環(huán)境。將無(wú)人機(jī)傾斜攝影技術(shù)應(yīng)用于露天礦山資源開(kāi)發(fā)的全生命周期中，有利于進(jìn)一步推動(dòng)礦山生態(tài)修復(fù)工作的開(kāi)展。研究小組以淮北市三五山土地綜合整治項(xiàng)目為例，對(duì)其精度進(jìn)行評(píng)價(jià)，淺析無(wú)人機(jī)技術(shù)在礦山生態(tài)治理監(jiān)測(cè)中的應(yīng)用。

1 無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)

無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)以無(wú)人機(jī)為航攝平臺(tái)，搭載多種航攝傳感器，獲取地面垂直和傾斜角度多個(gè)方向的影像二維相片或掃描的點(diǎn)云數(shù)據(jù)。數(shù)據(jù)獲取后，工作人員通過(guò)多光譜影像分析、正射影像制作、DSM 和DLG 數(shù)據(jù)制作等方法，結(jié)合機(jī)載POS 數(shù)據(jù)、像控測(cè)量成果及其他矢量數(shù)據(jù)，構(gòu)建地物與地貌的真實(shí)三維模型[3]。

主要操作流程包括：航線飛行設(shè)計(jì)、像控點(diǎn)布設(shè)與測(cè)量、無(wú)人機(jī)航攝、空三解算（約束網(wǎng)平差、刺點(diǎn)、影像接邊處理平差）和三維建模。

2 研究區(qū)及實(shí)驗(yàn)數(shù)據(jù)獲取

2.1 研究區(qū)概況

淮北市烈山區(qū)三五山土地綜合整治項(xiàng)目位于淮北市東南約9.0 km，總面積約63.14 hm2，分為三五山東部和三五山西部，其中東部面積約26.47 hm2，西部面積約36.67 hm2。研究區(qū)為廢棄采石宕口，區(qū)內(nèi)植被已被破壞，基礎(chǔ)條件較差，地形錯(cuò)綜復(fù)雜，垂直高差不一，坡度較大，最大垂直高差達(dá)82 m，選取該區(qū)域進(jìn)行實(shí)驗(yàn)具有一定的代表性。

2.2 數(shù)據(jù)來(lái)源

研究小組利用飛馬D2000 無(wú)人機(jī)搭載傾斜正射模塊，在研究區(qū)分別航飛3 個(gè)架次，每個(gè)架次間隔1 個(gè)月，3 次航飛分別在2022 年3月23 日、4 月24 日和5 月21 日。如圖1 所示，根據(jù)研究區(qū)地形，研究小組在縱向航線上按照研究區(qū)走向直線布設(shè)，平行于研究區(qū)邊界線的首末航線必須確保側(cè)視鏡頭能獲得研究區(qū)的有效影像；橫向上，選擇從東到西布設(shè)航線。行高109 m，航向重疊度67%，旁向重疊度47%，時(shí)速10.8 m/s，航向間距50 m。

圖1 測(cè)區(qū)航線布設(shè)圖

通過(guò)航飛，研究小組獲得影像照片1 671張（第1 期544 張，第2 期541 張，第3 期586 張）、3 期的POS 數(shù)據(jù)和RTK 采集的像控點(diǎn)坐標(biāo)。相機(jī)型號(hào)SONY ILCE-6000，影像大小6 000×4 000。

2.3 數(shù)據(jù)處理方法

數(shù)據(jù)處理包括空三解算和三維建模。無(wú)人機(jī)數(shù)據(jù)采集與地面數(shù)據(jù)采集同步進(jìn)行，研究小組采用Context Capture 軟件，結(jié)合影像的外方位元素，進(jìn)行全自動(dòng)空三解算，采用人工智能交互技術(shù)將數(shù)據(jù)進(jìn)行融合，生成一體化的三維模型[4]。三維模型生產(chǎn)完成后應(yīng)保證模型的完整性、位置精度、表現(xiàn)精度、屬性精度、現(xiàn)勢(shì)性等要求[5]，測(cè)繪航空攝影處理流程圖如圖2所示。

圖2 測(cè)繪航空攝影處理流程圖

3 實(shí)驗(yàn)結(jié)論和精度討論

3.1 三維模型成果

此次三維模型生成使用ArcGIS 軟件，將Context Capture 軟件生產(chǎn)的DSM 作為地形數(shù)據(jù)，DOM 作為影像數(shù)據(jù)，通過(guò)紋理映射技術(shù)將DOM 數(shù)據(jù)貼合在DSM 上，形成三維影像[6]。3 期模型對(duì)比如圖3 所示，模型與現(xiàn)場(chǎng)對(duì)比如圖4 所示。

圖3 3 期模型對(duì)比圖

圖4 模型與現(xiàn)場(chǎng)對(duì)比圖

整體上看，三維模型整體上較好地還原了礦山地物形態(tài)，礦山開(kāi)挖狀況與現(xiàn)場(chǎng)表現(xiàn)一致。治理區(qū)域原有的房屋建筑、植被樹(shù)木都能較完整地體現(xiàn)出來(lái)；局部上，模型色彩相對(duì)厚重，尤其是高低變化較大的區(qū)域，色調(diào)過(guò)度起伏較大，突出了礦山開(kāi)采的變化情況；空間形變上，區(qū)域內(nèi)的植被與礦山高低起伏一致，礦區(qū)廠房、礦坑、覆蓋痕跡等均在模型上有所體現(xiàn)。研究小組認(rèn)為，基于無(wú)人機(jī)傾斜攝影測(cè)量構(gòu)建的研究區(qū)三維模型可以直觀形象地表達(dá)礦區(qū)的地物特征，構(gòu)建的三維模型完整全面，細(xì)節(jié)層次清晰。

3.2 數(shù)據(jù)評(píng)價(jià)

研究小組以徠卡GS18 采集的檢查點(diǎn)數(shù)據(jù)為真實(shí)值，與ArcGIS 軟件提取的對(duì)應(yīng)影像高程值進(jìn)行對(duì)比，對(duì)生成的影像數(shù)據(jù)進(jìn)行檢驗(yàn)，精度統(tǒng)計(jì)原理如下。

式（1）至式（4）中，Δx、Δy、Δh分別為每個(gè)檢查點(diǎn)在X、Y、H 方向上的殘差值；n為檢查數(shù)量；mx為X 方向的中誤差；my為Y 方向的中誤差；ms為平面中誤差；mh為高程中誤差。

項(xiàng)目檢查點(diǎn)共計(jì)27 個(gè)，研究小組將3 期檢查數(shù)據(jù)對(duì)比分析。1 期27 個(gè)點(diǎn)平面誤差主要分布在0.033～0.241 m 之間，其中平面誤差在0.033～0.137 m 之間的點(diǎn)共有24 個(gè)，占比89%；0.137～0.241 m 之間的點(diǎn)共有3 個(gè)，占比11%；高程誤差分布在0.054～0.123 m之間，其中高程誤差在0.054～0.09 m 之間的點(diǎn)共有13 個(gè)，占比48%；其中高程誤差在0.09～0.123 m 之間的點(diǎn)共有14 個(gè)，占比52%。平面誤差最大為0.241 m，中誤差為±0.118 m；高程誤差最大為0.123 m，中誤差為±0.089 m。

2 期27 個(gè)點(diǎn)平面誤差主要分布在0.07～0.134 m 之間，其中平面誤差在0.07～0.102 m之間的點(diǎn)共有16 個(gè)，占 比59%；0.102～0.134 m 之間的點(diǎn)共有11 個(gè)，占比41%；高程誤差分布在0.01～0.2 m 之間，其中高程誤差在0.01～0.11 m 之間的點(diǎn)共有20 個(gè)，占比74%；其中高程誤差在0.11～0.2 m之間的點(diǎn)共有7 個(gè)，占比26%。平面誤差最大為0.134 m，中誤差為±0.102 m，高程誤差最大為 0.2 m，中誤差為±0.093 m。

3 期27 個(gè)點(diǎn)平面誤差主要分布在0.059～0.135 m之間，其中平面誤差在0.059～0.097 m 之間的點(diǎn)共有11 個(gè)，占比41%；0.097～0.135 m 之間的點(diǎn)共有16 個(gè)，占比59%；高程誤差分布在0.05～0.128 m之間，其中高程誤差在0.05～0.089 m 之間的點(diǎn)共有21 個(gè)，占比78%；其中高程誤差在0.89～0.128 m 之間的點(diǎn)共有6 個(gè)，占比22%。平面誤差最大為0.135 m，中誤差為±0.101 m，高程誤差最大為 0.128 m，中誤差為±0.084 m。

根據(jù)《三維地理信息模型數(shù)據(jù)產(chǎn)品規(guī)范》（CH/T9015-2012），比例尺為1 ∶500 的平面精度應(yīng)控制在0.3 m 以內(nèi)，高程精度應(yīng)控制在0.5 m 內(nèi)。根據(jù)實(shí)驗(yàn)結(jié)果，研究區(qū)平面中誤差最大為0.118 m，最大誤差0.241 m，數(shù)據(jù)＜0.3 m 的標(biāo)準(zhǔn)；高程中誤差最大為0.093 m，最大誤差0.2 m，數(shù)據(jù)＜0.5 m 的標(biāo)準(zhǔn)。由此可見(jiàn)，此次實(shí)驗(yàn)結(jié)果符合《三維地理信息模型數(shù)據(jù)產(chǎn)品規(guī)范》精度Ⅰ級(jí)要求。

4 礦山地質(zhì)環(huán)境治理監(jiān)測(cè)實(shí)例應(yīng)用分析

研究小組通過(guò)無(wú)人機(jī)傾斜攝影技術(shù)，按流程操作獲得3 期實(shí)際數(shù)據(jù)，利用ArcGIS 軟件和南方CASS 軟件，對(duì)符合要求的DSM 和DOM 影像進(jìn)行處理，將標(biāo)準(zhǔn)范圍線套合在新生成的DOM 影像中。通過(guò)3 期數(shù)據(jù)對(duì)比分析、數(shù)據(jù)平面量測(cè)，研究小組獲得需要重點(diǎn)監(jiān)測(cè)、治理的露天礦山的開(kāi)采范圍、開(kāi)采面積、開(kāi)采標(biāo)高、回填信息、礦山地質(zhì)環(huán)境情況等基礎(chǔ)信息，并判斷研究區(qū)開(kāi)采東南方向局部擴(kuò)張，西部擴(kuò)大，在研究區(qū)范圍界線內(nèi)，不存在越界開(kāi)采情況。

通過(guò)生成高程DEM 數(shù)據(jù)，研究小組建立了研究區(qū)前后開(kāi)挖對(duì)比影像圖（見(jiàn)圖5），形象直觀地展布了整個(gè)礦山環(huán)境治理區(qū)的范圍與地形概況，提升了監(jiān)測(cè)與評(píng)價(jià)礦山地質(zhì)環(huán)境治理的效果。如圖5 所示，研究區(qū)中心存在超挖現(xiàn)象（黃色圖斑），研究區(qū)南北方向存在漏挖、少挖的情形（藍(lán)色圖斑）。

圖5 研究區(qū)前后開(kāi)挖對(duì)比影像圖

5 技術(shù)難點(diǎn)

攝影測(cè)量過(guò)程中，由于地物遮擋、空間位置起伏變化、分辨率的差異，獲得的數(shù)據(jù)當(dāng)中含有非常多的粗差，從而生成大量冗余信息，單獨(dú)利用一種匹配方式無(wú)法獲得建模所需數(shù)據(jù)。

因此，研究小組通過(guò)實(shí)地勘驗(yàn)，按照平均150 m 的間隔布設(shè)像控點(diǎn)，并保證其均勻分布在測(cè)區(qū)范圍內(nèi)。對(duì)于地形起伏異常較大、大面積植被及面狀水域特征點(diǎn)非常少的區(qū)域，應(yīng)酌情增加像控點(diǎn)布設(shè)密度；單個(gè)攝影測(cè)量區(qū)域像控點(diǎn)的布設(shè)，其外圍的像控點(diǎn)連線要覆蓋整個(gè)區(qū)域。調(diào)整后的布設(shè)方式可以有效辨析地物特性，多層次獲取影像信息，分析多視影像特征，將影像上不同視角的二維特征轉(zhuǎn)化為三維特征，實(shí)現(xiàn)三維模型的影像匹配。

6 結(jié)語(yǔ)

研究小組利用無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)對(duì)淮北市烈山區(qū)三五山土地綜合整治項(xiàng)目進(jìn)行監(jiān)測(cè)數(shù)據(jù)分析，得到以下結(jié)論。

（1）無(wú)人機(jī)傾斜攝影測(cè)量平面中誤差最大為0.118 m，最大誤差0.241 m；高程中誤差最大為0.93 m，最大誤差0.2 m。平面中誤差均小于0.3 m，高程中誤差均小于0.5 m，滿足《三維地理信息模型數(shù)據(jù)產(chǎn)品規(guī)范》（CH/T9015-2012）中規(guī)定的Ⅰ級(jí)精度誤差標(biāo)準(zhǔn)。

（2）在礦山治理環(huán)境監(jiān)測(cè)應(yīng)用中，研究小組通過(guò)建立實(shí)景三維模型對(duì)比，認(rèn)為研究區(qū)開(kāi)采東南方向局部擴(kuò)張，西部擴(kuò)大；在研究區(qū)范圍界線內(nèi)，不存在越界開(kāi)采情況；研究區(qū)中心存在超挖現(xiàn)象，研究區(qū)南北方向存在漏挖少挖的情形。

（3）該項(xiàng)目的生產(chǎn)實(shí)踐證明，無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)能夠在保證精度滿足要求的前提下真實(shí)完整地再現(xiàn)研究區(qū)治理現(xiàn)狀。下一步應(yīng)綜合考慮設(shè)備、環(huán)境、算法等影響因素，選取更加精準(zhǔn)、高效、復(fù)雜的精度評(píng)價(jià)因子，建立更嚴(yán)謹(jǐn)?shù)木仍u(píng)價(jià)體系，讓自動(dòng)生成的三維模型更加直觀全面地反映實(shí)際生產(chǎn)情況。