毛崔磊　田毅　陳建平

摘要：針對廢棄礦山原地形地貌、堆棄物的規(guī)模、堆棄物的穩(wěn)定性未知等問題，以北京市密云某礦區(qū)為例，采用無人機技術(shù)獲取研究區(qū)航測數(shù)據(jù)，提出基于Agisoft PhotoScan的影像三維建模以及精度檢驗方法，利用Agisoft PhotoScan軟件進行定向、點云提取和正射影像圖的制作，生成了研究區(qū)的三維模型視圖。結(jié)果表明，該方法可以用于建立廢棄礦山的三維模型，縮短了建模的時間，有效降低了成本，為今后廢棄礦山地質(zhì)災(zāi)害信息獲取提供新的思路。

關(guān)鍵詞：無人機技術(shù)；Agisoft PhotoScan軟件；三維建模

中圖分類號： S127文獻標志碼：

文章編號：1002-1302（2017）16-0198-04

收稿日期：2017-02-08

基金項目：北京市教育委員會課題（編號：京教函〔2015〕282號）；國土資源部公益性行業(yè)科研專項（編號：201411007-3）。

作者簡介：毛崔磊（1992—），男，河南新鄉(xiāng)人，碩士，主要從事礦山資源環(huán)境研究。E-mail：1042274447@qqcom。

通信作者：田毅，博士，副教授，主要從事土地綜合整治、城市地下空間資源定量評價的教學(xué)與科研。E-mail：tianyi@cugbeducn。

社會經(jīng)濟的快速發(fā)展加快了對礦產(chǎn)資源的開發(fā)與利用，在帶給社會效益的同時，由于長期不當?shù)牡V產(chǎn)開發(fā)，對于礦山原地形地貌破壞較大，同時潛在的致災(zāi)隱患不斷增多，隨時可能發(fā)展成地質(zhì)災(zāi)害，常見的礦區(qū)地質(zhì)災(zāi)害為地面塌陷、邊坡失穩(wěn)、滑坡等，且礦區(qū)周邊常有居民點，一旦災(zāi)害發(fā)生將會造成不可估量的經(jīng)濟損失，甚至人員傷亡。

由于長期的開采，礦山廢棄物隨意堆棄，粗放的堆棄導(dǎo)致原地形原地貌未知，堆棄物的規(guī)模未知，堆棄物的穩(wěn)定性未知，以及可能導(dǎo)致潛在的礦山環(huán)境災(zāi)害未知，因此對于廢棄礦山的地形地貌的精準測量變得尤為重要。

近年來，隨著無人機航拍技術(shù)的不斷成熟，加之無人機低成本，機動性強，大比例尺高精度、維護簡單；在小區(qū)域和飛行困難地區(qū)高分辨率影像快速獲取方面具有明顯優(yōu)勢，如在災(zāi)害發(fā)生的區(qū)域進行快速的影像獲取，以及礦山地形圖的快速更新，使得無人機逐漸成為獲取更高精度的影像的主要技術(shù)手段之一。

本研究以北京密云某廢棄礦山為例，提出了利用無人機遙感技術(shù)獲取的影像來進行三維地形建模，并對廢棄礦山地貌精準測量方法進行探討[2-3]。

1三維重建方法及軟件

11研究區(qū)簡介

密云區(qū)位于北京市東北部，屬燕山山地與華北平原交接地，總地形為三面環(huán)山，中部低緩，西南開口的簸箕形。密云區(qū)為暖溫帶季風型大陸性半濕潤半干旱氣候。四季分明，干濕冷暖變化明顯。地理坐標西起116°39′33″E，東至117°30′25″E，東西長69 km；南起40°13′7″N，北至 40°47′57″N，是華北通往東北、內(nèi)蒙古的重要門戶。

密云區(qū)的礦產(chǎn)資源豐富，已探明儲量的礦產(chǎn)有25種，主要礦產(chǎn)有金、銀、鐵、鎢、鉛等。鐵礦儲量達95億t，在全國 2 000 多個縣中居19位，占北京地區(qū)總儲量的964%，石灰?guī)r的儲量也有近億t。

12原始數(shù)據(jù)的特點及來源

以北京密云某廢棄礦山為研究區(qū)域，通過無人機利用非量測相機對該地區(qū)進行航空攝影測量（所用sony a5100相機參數(shù)見表1，相機見圖1）。本研究區(qū)設(shè)計飛行相對高度為 250 m，絕對高度為450 m，無人機上搭載的相機焦距為 285 mm，測區(qū)范圍總面積達88 km2，獲得原始照片663張，同時保證航空攝影測量時60%側(cè)向重疊和80%航向重疊的規(guī)范要求。采用GeoFly-X8固定翼無人機，相對于多旋翼無人機其速度更快、效率更高、更節(jié)能，適合大范圍的高清影像的拍攝，可作業(yè)范圍甚至可達幾十平方公里，是農(nóng)業(yè)監(jiān)測、礦山監(jiān)測，災(zāi)情監(jiān)測的首選機型。該無人機飛行高度可達 1 200 m，影像精度為005～010 m，飛行速度35～70 km/h。

121影像特征

在實際飛行中，該區(qū)域1個架次完成，共獲取663張照片，其中有效照片658張。對應(yīng)的照片信息見圖2。該任務(wù)區(qū)內(nèi)包含建筑、道路、植被，各個類型的航拍原始圖片細節(jié)見圖3、圖4、圖5。

13三維建模軟件（Agisoft PhotoScan）

Agisoft PhotoScan是一款基于靜態(tài)圖像構(gòu)建高精度三維模型的軟件。依據(jù)最新的計算機多目視覺影像三維重建技術(shù)，就可以對具有影像重疊的照片進行處理；也可以通過給予的控制點生成真實坐標的三維模型。其在三維場景中的應(yīng)用相當廣泛，如考古（Agisoft Photoscan 和 Adobe Illustrator軟件在田野考古繪圖中的綜合應(yīng)用）、工業(yè)、軍事、醫(yī)學(xué)等諸多領(lǐng)域都可以使用該軟件進行處理。整個工作流程都是完全自動化的。

Agisoft PhotoScan區(qū)別于其他傳統(tǒng)航攝軟件的地方在于它使用的是運動信息中恢復(fù)三維場景結(jié)構(gòu)（structure from motion，SfM）的方法。SfM三維重建的主要過程如下：（1）利用SIFT算子從每張影像上提取特征點并獲取其對應(yīng)的Descriptor；（2）根據(jù)POS或其他約束關(guān)系選出可能具有重疊關(guān)系的像對；（3）對每一像對Descriptor進行匹配，并使用RANSAC算法進行粗差剔除，消除誤匹配；（4）根據(jù)計算機多目視覺原理，將各像對匹配的同名像點連接統(tǒng)一起來，在逐次平差迭代過程中，剔除粗差，估算出視覺場景中相機和由匹配點形成的稀疏點云的相對位置；（5）結(jié)合相機的GPS和像控點的位置坐標，通過Damped Least-Squares（DLS）算法來獲取相機和稀疏點云的真實空間位置。

2三維模型構(gòu)建

后期影像處理基于POS數(shù)據(jù)進行，經(jīng)過一系列的處理后，得到研究區(qū)的三維模型[4]。處理流程見圖6。

21原始數(shù)據(jù)的預(yù)處理

由于無人機類別及用途存在差異， 無人機所獲取的影像

POS信息的文件格式也有所區(qū)別。首先要將POS數(shù)據(jù)格式按照處理軟件的要求進行修改，使其能夠順利地導(dǎo)入到工程當中。將數(shù)據(jù)格式依次修改為緯度、經(jīng)度、海拔高度、航向、俯仰、橫滾。

22導(dǎo)入影像

在導(dǎo)入之前，檢查影像的質(zhì)量，刪除模糊、變形、色彩異常的影像，提高數(shù)據(jù)的精度。本次拍攝的區(qū)域較大，獲取的影像數(shù)量較多，分區(qū)塊依次單獨進行處理，對每一個塊逐一進行照片對齊、點云提取和立體模型的建立，最后再進行合并、提取紋理。分區(qū)塊進行處理的好處在于，既可以充分利用計算機的性能特點，還可以提高影像數(shù)據(jù)處理的效率。需要注意的是，導(dǎo)入的影像必須是依照航線且連續(xù)的影像，與POS數(shù)據(jù)一一對應(yīng)的數(shù)據(jù)，中間不能有遺漏。

23數(shù)據(jù)定向、點云提取

Agisoft photoscan根據(jù)空中三角測量的基本原理，根據(jù)多目視覺影像三維模型重建技術(shù)，自動解算出影像的坐標、姿態(tài)等信息，同時自動完成影像匹配、內(nèi)定向、相對定向及絕對定向等過程，處理過程中不再需要進行人工干預(yù)。后期處理所需要的基本數(shù)據(jù)是影像數(shù)據(jù)、POS數(shù)據(jù)和地面控制點數(shù)據(jù)，影像的拍攝姿態(tài)與位置是由無人機在飛行過程中根據(jù)航線的設(shè)計和飛行的姿態(tài)來決定。

24立體建模

利用Agisoft photoscan從無人機航拍的原始影像數(shù)據(jù)和位置信息提取帶有坐標信息的點云數(shù)據(jù)，可以快速重新構(gòu)建任務(wù)目標的三維模型，即線、面、體、空間等各種數(shù)據(jù)，還原任務(wù)目標真實的形態(tài)特性。

25模型賦予紋理

模型構(gòu)建完成后，需要將模型賦予材質(zhì)、著色、紋理，Agisoft photoscan軟件對于紋理的生成具有5種不同的模式，包括通用、正射影像、自適應(yīng)正射影像、球形、單一相機。根據(jù)研究區(qū)特點采用適合研究區(qū)域的自適應(yīng)正射影像模式，以達到符合實際的效果，保證質(zhì)量。

26成果輸出

完成上述步驟后，就可以根據(jù)實際的需要將數(shù)字高程模型（DEM）、數(shù)字正射影像（DOM）、數(shù)字表面模型（DSM）導(dǎo)出，數(shù)字高程模型（DEM）即地形表面形態(tài)的數(shù)字化表達，它是用一組有序數(shù)值陣列形式表示地面高程的一種實體地面模型。數(shù)字正射影像（DOM）是對航空（或航天）相片進行數(shù)字微分糾正和鑲嵌，按一定圖幅范圍裁剪生成的數(shù)字正射影像集。它是同時具有地圖幾何精度和影像特征的圖像。數(shù)字表面模型（DSM）指物體表面形態(tài)以數(shù)字表達的集合。其中DEM必須是高程信息，是地表的模擬，DSM可以是地物表面的模擬，包括植被表面、房屋的表面，對DSM進行加工，去掉房屋、植被等信息，可以形成DEM。成果輸出時，根據(jù)需求調(diào)整像元大小、投影坐標的類型，選擇分幅圖或整幅圖輸出。輸出的格式常用為TIFF/GeoTIFF（tif）格式，也可導(dǎo)出KML文件，在google earth上打開，數(shù)字地表模型（DSM）也可導(dǎo)出為PDF或3DS等多種格式。三維模型見圖7。

三維模型的精度決定了其可以滿足的不同應(yīng)用需求，為檢驗軟件生成的三維模型，在實地采集了32個控制點，其中7個作為匹配點，對三維模型進行匹配，將控制點輸入軟件，坐標為wgs84，匹配模型見圖8，導(dǎo)入的GPS點與軟件匹配后的點進行比較，計算得到Z值誤差數(shù)據(jù)（表2）。

從圖9可以看出，誤差總體控制在05 m左右，由于受實際條件的限制，以及天氣原因的影響，基本上達到了預(yù)期的效果，與實際外業(yè)測量的真實坐標差距不是很大，因此構(gòu)建的三維模型具有較高的精度，可以滿足實際需求[5]。

3結(jié)論

本研究以廢棄礦山為例，通過無人機獲取的原始數(shù)據(jù)，運用Agisoft PhotoScan軟件建立了相應(yīng)的三維模型，實現(xiàn)了廢棄礦山的精準測量，試驗結(jié)果表明，基于Agisoft PhotoScan可以快速得到研究區(qū)的三維模型，并結(jié)合精度檢驗，可以滿足相應(yīng)的需求，為廢棄礦山的地貌精準測量提供了幫助，相比與傳統(tǒng)的三維建模方式縮短了時間，降低了成本，具有更好的應(yīng)用前景[6-7]。

同時需要注意GCP點應(yīng)布置在高對比度、可清楚地確定出中心位置的目標。這樣才能提高影像坐標和空間坐標的匹配度與可靠性。實際實施中受一些客觀條件的限制未能將GCP點做到最優(yōu)化。下一步需要優(yōu)化GCP點位置，從而提高三維模型的精度[8-10]。

今后將研究局部微地貌的表達，地形的緩陡程度，將三維數(shù)據(jù)與實地數(shù)據(jù)進行對比分析，進一步驗證該方法的可行性與快捷性。

參考文獻：

王俊強，韓宏強，孫建霖 基于無人機低空攝影高分辨率三維地形景觀的實現(xiàn)[J] 測繪技術(shù)裝備，2013（4）：10-13

王鈞超 礦山環(huán)境評價方法綜述與治理經(jīng)驗探討[J] 中國礦業(yè)，2014（增刊1）：103-105

[3]王來貴，潘一山，趙娜 廢棄礦山的安全與環(huán)境災(zāi)害問題及其系統(tǒng)科學(xué)研究方法[J] 渤海大學(xué)學(xué)報（自然科學(xué)版），2007，28（2）：97-101

[4]李秀全，陳竹安，張立亭 非量測相機影像三維模型構(gòu)建及精度檢驗[J] 測繪科學(xué)，2016，41（6）：144-147

[5]李秀全，陳竹安，張立亭 基于Agisoft PhotoScan的無人機影像快速拼接在新農(nóng)村規(guī)劃中的應(yīng)用[J] 湖北農(nóng)業(yè)科學(xué)，2016，55（3）：743-745

[6]黎富忠 無人機航拍技術(shù)在三維建模上的應(yīng)用[J] 廣西水利水電，2016（4）：21-23

[7]王鶴，劉軍，王秋玲 利用無人機影像進行滑坡地形三維重建[J] 測繪與空間地理信息，2015（12）：68-71

[8]李文慧，楊斌，黃永璘，等 無人機遙感在三維地形建模中的應(yīng)用初探[J] 氣象研究與應(yīng)用，2008，29（4）：38-41

[9]魏占玉，Ramon A，何宏林，等 基于SfM方法的高密度點云數(shù)據(jù)生成及精度分析[J] 地震地質(zhì)，2015，37（2）：636-648

[10]趙云景，龔緒才，杜文俊，等 PhotoScan Pro軟件在無人機應(yīng)急航攝中的應(yīng)用[J] 國土資源遙感，2015，27（4）：179-182