基于無人機(jī)影像序列的城市精細(xì)化三維模型精度評(píng)估

陳仁朋　吳熠文　余加勇　何曠宇　薛現(xiàn)凱　李鋒

摘? ?要：為滿足數(shù)字城市對(duì)精細(xì)化實(shí)景三維模型日益顯著的需求，解決傳統(tǒng)方法建模效率低的問題，提出一種基于無人機(jī)影像序列的城市精細(xì)化三維模型構(gòu)建方法，通過現(xiàn)場(chǎng)25架次的航拍試驗(yàn)，構(gòu)建了20個(gè)300 m×300 m測(cè)區(qū)范圍的精細(xì)化模型和1個(gè)圖書館建筑精細(xì)化模型;利用檢查點(diǎn)評(píng)估了全部精細(xì)化模型的精度，重點(diǎn)討論了5種航拍高度和4種控制點(diǎn)選取方案對(duì)精細(xì)化模型精度的影響. 分析結(jié)果表明：隨著航拍高度降低、控制點(diǎn)數(shù)量增加，構(gòu)建的精細(xì)化模型誤差減小、精度得到提高;測(cè)區(qū)范圍精細(xì)化三維模型最高精度達(dá)到平面±3.4 cm、高程±3.1 cm，圖書館單體建筑精細(xì)化模型精度達(dá)到平面±3.4 cm、高程±1.5 cm;基于無人機(jī)序列影像的三維模型構(gòu)建方法能夠滿足城市實(shí)景模型構(gòu)建的精度要求，具有重要工程實(shí)用價(jià)值.

關(guān)鍵詞：無人機(jī);傾斜攝影測(cè)量;三維模型;精度評(píng)估

中圖分類號(hào)：P237? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Method Accuracy Evaluations of Building Urban Detailed 3D Model Based on

the Unmanned Aerial Vehicle Image Sequences and Its Accuracy Evaluatios

CHEN Renpeng1，2，WU Yiwen1，2，YU Jiayong1，2，HE Kuangyu1，2，XUE Xiankai1，2，LI Feng1，2

（1. Key Laboratory of Building Safety and Energy Efficiency of Ministry of Education，Hunan University，Changsha 410082，China;

2. College of Civil Engineering，Hunan University，Changsha 410082，China）

Abstract： In order to satisfy the increasingly significant demand of the digital city for the detailed 3D model and solve the low efficiency problem of traditional modeling method， this paper summarizes a method of building the urban detailed model based on unmanned aerial vehicle （UAV） image sequences. Through the in-situ aerial photography test with 25 sorties， detailed models for 20 experiment area （300 m×300 m） and a library building were built. The accuracy of all detailed models were evaluated by field check points. The effect of 5 different flight heights and 4 different control points on the accuracy of the detailed models were discussed. The results show that the errors of the detailed models can be effectively decreased and the accuracy is improved by declining the flight height or increasing the number of control point. The highest accuracy of detailed models is ±3.4 cm in the plane and ±3.1 cm in the height， respectively. The accuracy of the detailed model for the library building is ±3.4 cm and ±1.5 cm， respectively. The method of building the urban detailed model using UAV image sequences can satisfy the accuracy requirement of the urban real scene model， which has great practical values in practical engineering.

Key words：Unmanned Aerial Vehicle（UAV）;oblique photogrammetry;3D reconstruction;accuracy evaluation

隨著數(shù)字城市對(duì)城市精細(xì)化實(shí)景三維模型的需求日益顯著，建立真實(shí)的三維景觀凸顯出較高的應(yīng)用前景和經(jīng)濟(jì)價(jià)值[1].? ?數(shù)字城市三維模型包括地形和建筑三維模型. 傳統(tǒng)航空攝影測(cè)量方法不僅影像質(zhì)量和分辨率相對(duì)較差[2]，而且構(gòu)建的三維建模側(cè)面紋理信息相對(duì)缺乏，需要依靠人工從地面拍攝照片來補(bǔ)充紋理信息，對(duì)于城市內(nèi)數(shù)以萬計(jì)的建筑物來說，此方法耗時(shí)耗力[3].

無人機(jī)（Unmanned Aerial Vehicle，UAV）具有操作容易、數(shù)據(jù)采集靈活、時(shí)效性高、云下飛行等優(yōu)

點(diǎn)[4-5]，已成為主要的低空遙感平臺(tái)[6]. 無人機(jī)傾斜攝影測(cè)量技術(shù)是一項(xiàng)新興的對(duì)地觀測(cè)技術(shù)[7]，顛覆了傳統(tǒng)航空攝影測(cè)量技術(shù)，它通過無人飛行平臺(tái)搭載的多臺(tái)相機(jī)和傳感器，從一個(gè)垂直和四個(gè)傾斜角度獲取目標(biāo)的垂直和傾斜高質(zhì)量影像，不僅能夠獲取豐富的地物信息，而且通過定位、融合、建模技術(shù)，生成真實(shí)的三維城市模型[8].

目前國(guó)內(nèi)外研究者已開展了基于無人機(jī)傾斜攝影測(cè)量技術(shù)的城市三維建模方法研究[9-11]，然而對(duì)采用無人機(jī)超低空、低空航拍方法來構(gòu)建城市精細(xì)化三維模型的精度影響因素的研究較少[12-14]. 為此，本文提出一種基于無人機(jī)影像序列的城市精細(xì)化三維模型構(gòu)建方法，通過大量現(xiàn)場(chǎng)試驗(yàn)評(píng)估三維模型的平面和高程精度，重點(diǎn)分析航拍高度和控制點(diǎn)數(shù)量對(duì)三維模型精度的影響.

1? ?城市精細(xì)化模型構(gòu)建方法

1.1? ?無人機(jī)影像采集

通過無人機(jī)傾斜攝影測(cè)量技術(shù)，從不同相機(jī)位置獲取目標(biāo)區(qū)域多個(gè)不同角度影像，包括垂直拍攝影像和傾斜拍攝影像.

結(jié)合實(shí)際情況確定合理的飛行計(jì)劃，包括飛行高度、線路和速度、影像拍攝角度、相機(jī)參數(shù)、拍照頻率、航向和旁向影像重疊度等參數(shù)的設(shè)計(jì). 合理的飛行計(jì)劃能夠有效減少外業(yè)工作量，并提高后期處理效率和產(chǎn)品精度. 航拍實(shí)時(shí)記錄飛行高度、經(jīng)緯度、大地高、姿態(tài)等信息，用于后續(xù)的分析和整理.

1.2? ?影像畸變差糾正

畸變差主要來源除了相機(jī)鏡頭的機(jī)械和光學(xué)誤差以外，有時(shí)也可能源于A/D轉(zhuǎn)換時(shí)的電學(xué)誤差[15]. 畸變差破壞了出入射光線的平行狀態(tài)，使得像點(diǎn)偏離了理論位置，故不滿足三點(diǎn)共線條件. 相機(jī)畸變的糾正是徑向畸變和切向畸變的總和，即：

式中：xr和yr為徑向畸變;xd和yd為切向畸變. 通過標(biāo)定板利用相機(jī)檢校軟件對(duì)相機(jī)進(jìn)行標(biāo)定，獲得畸變參數(shù)，對(duì)影像進(jìn)行畸變差糾正.

1.3? ?影像勻色勻光

在影像獲取過程中，受快門速度、曝光度、天氣和光照等因素的影響，使得獲取的影像出現(xiàn)光照和色彩分布差異較大的現(xiàn)象，繼而影響三維模型的制作和應(yīng)用[16]. 基于Wallis勻光算法原理，選擇色彩和光照合適的模板影像，統(tǒng)計(jì)其均值和方差參數(shù)，借助Wallis濾波器對(duì)目標(biāo)影像灰度的線性分布進(jìn)行調(diào)整. 線性變換式為：

式中：g（x，y）為原始影像在（x，y）處的灰度值;f（x，y）為線性變換后影像在（x，y）處的灰度值;r0、r1分別為線性變換式中的加系數(shù)和乘系數(shù)[17].

1.4? ?影像區(qū)域網(wǎng)聯(lián)合平差

影像區(qū)域網(wǎng)聯(lián)合平差技術(shù)目的是獲得無人機(jī)不同拍攝位置和角度采集的影像上每個(gè)點(diǎn)的外方位元素. 通過由粗到精的金字塔匹配方法，結(jié)合POS信息得到外方位元素，對(duì)每級(jí)影像同名點(diǎn)進(jìn)行自動(dòng)匹配和自由網(wǎng)光束法平差，建立控制點(diǎn)坐標(biāo)、POS輔助信息、連接點(diǎn)線的影像自檢校區(qū)域網(wǎng)平差的誤差方程，經(jīng)過聯(lián)合解算，達(dá)到平差結(jié)果精度，其線性化誤差方程為：

式中：t和X分別表示影像定向未知數(shù)和坐標(biāo)未知數(shù)向量;A和B分別表示相應(yīng)的系數(shù)矩陣;l和Pφ分別表示影像坐標(biāo)觀測(cè)值向量及其權(quán)矩陣[18].

1.5? ?三維模型表面重建與紋理映射

在確定像機(jī)拍攝的外方位元素后，通過空中三角測(cè)量和自動(dòng)圖像匹配技術(shù)對(duì)場(chǎng)景進(jìn)行數(shù)字化重建，并輸出密集三維點(diǎn)云. 表面網(wǎng)格重建是曲面重建的重要前處理，根據(jù)區(qū)域的有限個(gè)點(diǎn)集將區(qū)域劃分為相連的三角面網(wǎng)絡(luò)，三角面的形狀和大小取決于不規(guī)則分布的測(cè)點(diǎn)的密度和位置，這樣能夠避免地形平坦時(shí)的數(shù)據(jù)冗余，又能按地形特征點(diǎn)表示數(shù)字高程特征.

最后進(jìn)行三維模型紋理映射，是增強(qiáng)模型視覺效果的關(guān)鍵步驟. 基于無人機(jī)傾斜攝影測(cè)量技術(shù)的紋理映射是多視圖紋理映射，即通過多幅來自不同視點(diǎn)的圖像作為紋理圖，具體實(shí)現(xiàn)方法為：數(shù)據(jù)獲取、可見性判斷、重新配置紋理分布、融合紋理邊界、紋理空洞修補(bǔ)等[19].

1.6? ?三維模型精度評(píng)估

根據(jù)CityGML標(biāo)準(zhǔn)對(duì)多細(xì)節(jié)層次（Levels of Detail，LOD）的劃分，針對(duì)建筑模型點(diǎn)位精度劃分為4個(gè)級(jí)別，如表1所示[20].

對(duì)于每個(gè)重構(gòu)的三維精細(xì)化模型，均對(duì)其平面和高程精度進(jìn)行評(píng)估. 首先采用自動(dòng)型全站儀高精度地獲取所有檢查點(diǎn)的實(shí)測(cè)三維坐標(biāo)（Xi，Yi，Zi），然后從重構(gòu)的三維模型中識(shí)別相應(yīng)檢查點(diǎn)的影像三維坐標(biāo)（XOi，YOi，ZOi）， 進(jìn)行誤差統(tǒng)計(jì)分析，采用式（4）（5）計(jì)算出三維模型在水平和高程方向的均方根誤差RMSEXY、RMSEZ：

式中：n為檢查點(diǎn)個(gè)數(shù). 并利用箱形圖對(duì)模型構(gòu)建誤差值進(jìn)行中誤差、平均值和離散性進(jìn)行統(tǒng)計(jì)分析. 箱形圖包含的數(shù)據(jù)點(diǎn)如圖1所示，當(dāng)上、下四分位數(shù)之差和上、下限值之差越小時(shí)，數(shù)據(jù)的離散性越小.

2? ?現(xiàn)場(chǎng)試驗(yàn)

2.1? ?試驗(yàn)區(qū)域

航拍試驗(yàn)區(qū)域位于長(zhǎng)沙市岳麓區(qū)中南大學(xué)新校區(qū)（東經(jīng)112.939 1°，北緯28.152 8°），以圖書館為中心并向四周延伸而形成300 m×300 m的正方形區(qū)域，如圖2所示. 試驗(yàn)區(qū)域內(nèi)包含建筑物、道路、橋梁、綠地等城市要素，可較好地模擬城市精細(xì)化三維模型構(gòu)建. 圖書館為試驗(yàn)區(qū)域內(nèi)最高建筑，其長(zhǎng)寬高分別為150 m、72 m、32 m.

2.2? ?試驗(yàn)設(shè)備

本次航拍試驗(yàn)采用大疆經(jīng)緯M210-RTK型號(hào)無人機(jī)，搭載大疆禪思X5S相機(jī)，最大像素2 080萬，鏡頭焦距15 mm，試驗(yàn)設(shè)備如圖3所示. 無人機(jī)前方和下方配備立體視覺系統(tǒng)，上方配備紅外感知系統(tǒng)，可實(shí)現(xiàn)前方、下方和上方避障. M210-RTK無人機(jī)配備有實(shí)時(shí)動(dòng)態(tài)差分（Real-time Kinematic，RTK）模塊，可實(shí)時(shí)接收全球定位系統(tǒng)（Global Positioning System，GPS）和北斗導(dǎo)航衛(wèi)星系統(tǒng)（Beidou Satellite System，BDS）信號(hào)，利用地面基準(zhǔn)站進(jìn)行載波相位差分改正，可獲取平面1 cm、高程2 cm精度的無人機(jī)飛行定位信息.

測(cè)點(diǎn)（控制點(diǎn)和檢查點(diǎn)）采用瑞士Leica TS30型號(hào)自動(dòng)型全站儀測(cè)量，其精密測(cè)距標(biāo)稱精度為0.6 mm +1 ppm，免棱鏡測(cè)距的標(biāo)稱精度為2 mm +2 ppm，測(cè)角標(biāo)稱精度為±0.5″. 并采用天寶DiNi12電子水準(zhǔn)儀復(fù)核測(cè)點(diǎn)高程，其標(biāo)稱精度為每公里往返觀測(cè)標(biāo)準(zhǔn)誤差為±0.3 mm.

2.3? ?控制點(diǎn)和檢查點(diǎn)選取

控制點(diǎn)和檢查點(diǎn)共兩組，第1組分布在300 m×300 m試驗(yàn)區(qū)域內(nèi)地面，共25個(gè)測(cè)點(diǎn)，用于區(qū)域內(nèi)三維模型構(gòu)建及其精度評(píng)估;第2組分布在圖書館建筑外表面，共98個(gè)測(cè)點(diǎn)，用于單體建筑精細(xì)化模型構(gòu)建及其精度評(píng)估.

第1組25個(gè)測(cè)點(diǎn)均勻分布在300 m×300 m試驗(yàn)區(qū)域，并埋設(shè)直徑1 cm、長(zhǎng)度7 cm的不銹鋼測(cè)量標(biāo)志. 設(shè)計(jì)4種不同的三維模型構(gòu)建方案，根據(jù)均勻分布控制的原則，分別選取其中3、6、9、12個(gè)測(cè)點(diǎn)作為控制點(diǎn)，其余點(diǎn)均為精度檢查點(diǎn)，控制點(diǎn)均勻分布在試驗(yàn)區(qū)內(nèi)，如圖4所示. 按照一級(jí)導(dǎo)線和二等水準(zhǔn)測(cè)量的技術(shù)要求，采用Leica TS30型號(hào)自動(dòng)型全站儀和天寶DiNi11型號(hào)電子水準(zhǔn)儀精確測(cè)量地面測(cè)點(diǎn)的平面和高程坐標(biāo)，其測(cè)量精度均優(yōu)于±3 mm. 控制點(diǎn)標(biāo)志采用長(zhǎng)寬均為1 m的旗幟布制作，旗幟布上印制黑白相間的輻射狀圖案及編號(hào)，圖案中心預(yù)留直徑3 mm的孔洞，便于對(duì)準(zhǔn)不銹鋼標(biāo)志的中心，如圖5所示.

第2組98個(gè)測(cè)點(diǎn)分布在圖書館建筑外表面上，按照10 m、20 m、30 m 3個(gè)不同高度分層布置，如圖6所示. 在建筑物4個(gè)側(cè)面的10 m和30 m高度處各取2個(gè)測(cè)點(diǎn)作為控制點(diǎn)，共16個(gè)控制點(diǎn)，其余測(cè)點(diǎn)作為檢查點(diǎn). 選取圖書館建筑外表面具有較強(qiáng)反射性且具有明顯特征的點(diǎn)位作為測(cè)點(diǎn)，不需額外粘貼反射標(biāo)志，直接采用Leica TS30自動(dòng)型全站儀的免棱鏡功能進(jìn)行三維坐標(biāo)測(cè)量，平面和高程測(cè)量精度均優(yōu)于±5 mm.

2.4? ?影像采集

為重構(gòu)試驗(yàn)區(qū)域內(nèi)地物和地貌的精細(xì)化三維模型，并評(píng)估三維模型的精度，共設(shè)計(jì)了5種影像采集方案，每個(gè)方案進(jìn)行1組垂直拍攝和4組傾斜拍攝. 5種方案中無人機(jī)飛行高度分別為80 m、100 m、120 m、140 m和160 m，對(duì)應(yīng)的影像分辨率和數(shù)量如表2所示. 采用DJI GS Pro軟件規(guī)劃無人機(jī)航拍路徑及控制無人機(jī)飛行，相機(jī)傾斜角設(shè)置為45°，航向重疊率和旁向重疊率均設(shè)置為70%. 共進(jìn)行了25架次航拍，共采集908張相片.

3? ?模型構(gòu)建與精度分析

3.1? ?模型構(gòu)建

對(duì)于300 m×300 m測(cè)區(qū)三維模型，本次試驗(yàn)設(shè)置了5個(gè)飛行高度和4種控制點(diǎn)方案，共重構(gòu)了20個(gè)三維模型;對(duì)于圖書館單體建筑精細(xì)化模型，本次試驗(yàn)利用飛行高度為80 m的影像序列進(jìn)行精細(xì)化三維模型構(gòu)建，重構(gòu)了1個(gè)三維模型. 首先對(duì)影像預(yù)處理、影像區(qū)域網(wǎng)聯(lián)合平差、導(dǎo)入控制點(diǎn)及坐標(biāo);然后按照空中三角測(cè)量生成密集三維點(diǎn)云、生成三維白模、表面網(wǎng)格重建生成三角網(wǎng)、紋理映射生成實(shí)景三維模型等步驟構(gòu)建三維模型，如圖7所示.

3.2? ?測(cè)區(qū)模型精度分析

對(duì)300 m×300 m測(cè)區(qū)范圍的整體三維模型進(jìn)行精度分析. 不同飛行高度、不同控制點(diǎn)數(shù)量下共構(gòu)建了20個(gè)三維模型. 計(jì)算出每個(gè)模型中所有檢查的誤差值，即實(shí)測(cè)三維坐標(biāo)與模型中識(shí)別三維坐標(biāo)值的差值.

不同飛行高度下，構(gòu)建的三維模型平面和高程精度如圖8所示. 飛行高度為80 m時(shí)，不同控制點(diǎn)數(shù)量（3、6、9和12個(gè)）對(duì)應(yīng)的三維模型平面方向均方根誤差分別為5.8 cm、4.9 cm、4.0 cm、3.4 cm，高程方向均方根誤差分別為12.1 cm、5.9 cm、5.1 cm、3.1 cm. 飛行高度為160 m時(shí)，不同控制點(diǎn)數(shù)量對(duì)應(yīng)的三維模型平面方向均方根誤差分別為14.9 cm、14.2 cm、12.5 cm、10.6 cm，高程方向均方根誤差分別為16.5 cm、10.0 cm、8.1 cm、7.0 cm.

當(dāng)無人機(jī)飛行高度增加時(shí)，各控制點(diǎn)方案的三維模型平面和高程誤差均逐漸增大，精度逐步降低. 當(dāng)飛行高度由160 m減小至140 m時(shí)，平面均方根誤差減小較多，飛行高度對(duì)模型精度的水平精度影響較為明顯. 當(dāng)飛行高度為80 m時(shí)，平面和高程精度最高. 原則上要求無人機(jī)飛行高度不低于建筑物1.5倍高度，在此前提下，當(dāng)飛行高度越低時(shí)，模型的平面和高程精度越高.

在不同控制點(diǎn)數(shù)量下構(gòu)建的三維模型精度如圖9所示. 當(dāng)控制點(diǎn)數(shù)量為3時(shí)，不同飛行高度（80 m、100 m、120 m、140 m、160 m）對(duì)應(yīng)的平面均方根誤差分別為5.8 cm、9.0 cm、10.1 cm、10.9 cm、14.9 cm，高程均方根誤差分別12.1 cm、14.4 cm、15.8 cm、16.4 cm、16.5 cm. 當(dāng)控制點(diǎn)個(gè)數(shù)為12時(shí)，不同飛行高度對(duì)應(yīng)的平面均方根誤差分別為3.4 cm、4.8 cm、6.2 cm、6.8 cm、10.6 cm，高程均方根誤差分別3.1 cm、5.4 cm、6.7 cm、6.9 cm、7.0 cm. 隨著控制點(diǎn)數(shù)量的增加，三維模型的平面和高程誤差均逐漸減小，精度越來越高. 因此，控制點(diǎn)數(shù)量越多時(shí)，模型的平面和高程精度越高. 當(dāng)控制點(diǎn)數(shù)量由3個(gè)增加至6個(gè)時(shí)，高程均方根誤差明顯減小，控制點(diǎn)數(shù)量對(duì)模型高程精度影響顯著.

當(dāng)飛行高度為80 m、控制點(diǎn)數(shù)量為12時(shí)，平面均方根誤差為3.4 cm，高程均方根誤差為3.1 cm，三維模型的平面精度和高程精度達(dá)到最高. 所有三維模型的水平及高程精度均滿足精細(xì)化模型精度要求，故在滿足模型精度的情況下可盡量減少控制點(diǎn)數(shù)量（但至少為3個(gè)）或提高飛行高度.

3.3? ?測(cè)區(qū)模型誤差離散性分析

對(duì)300 m×300 m測(cè)區(qū)范圍的整體三維模型進(jìn)行兩組誤差離散性分析，第1組選取控制點(diǎn)數(shù)量為12時(shí)構(gòu)建的5種不同飛行高度的三維模型;第2組選取飛行高度為80 m時(shí)構(gòu)建的4種不同控制點(diǎn)方案的三維模型.

在第1組5個(gè)三維模型中，控制點(diǎn)數(shù)量均為12個(gè)，飛行高度由80 m 增加到160 m，誤差箱形圖如圖10所示. 當(dāng)飛行高度為160 m時(shí)，平面誤差上限值為16.0 cm，下限值為1.7 cm，上四分位數(shù)為13.8 cm，下四分位數(shù)為6.9 cm，中位數(shù)為10.0 cm，平均數(shù)為9.7 cm;高程誤差上限值為10.6 cm，下限值為2.2 cm，上四分位數(shù)為9.5 cm，下四分位數(shù)為2.7 cm，中位數(shù)為4.7 cm，平均數(shù)為5.9 cm. 隨著飛行高度的增加，模型平面和高程誤差的上、下四分位數(shù)之差和上、下限值逐步增大，離散性逐漸增大，精度降低. 當(dāng)飛行高度為80 m時(shí)，誤差離散性最小，精度最高;當(dāng)飛行高度由120 m降低至100 m時(shí)，平面誤差和高程誤差離散性明顯減小;當(dāng)飛行高度小于120 m時(shí)，誤差離散性較小;當(dāng)飛行高度大于120 m時(shí)，誤差離散性較大.

在第2組4個(gè)三維模型中，飛行高度均為80 m，控制點(diǎn)數(shù)量由3個(gè)增加至12個(gè)，誤差箱形圖如圖11所示. 隨著選取控制點(diǎn)數(shù)量的增加，三維模型誤差的上、下四分位數(shù)之差和上、下限值逐步減小，離散性逐步減小，精度提高. 當(dāng)控制點(diǎn)個(gè)數(shù)為12個(gè)時(shí)，離散性達(dá)到最小，精度最高;當(dāng)控制點(diǎn)數(shù)量由3個(gè)增加至6個(gè)時(shí)，平面誤差和高程誤差離散性明顯減小. 故當(dāng)控制點(diǎn)數(shù)量小于6個(gè)時(shí)，誤差離散性較大;當(dāng)控制點(diǎn)數(shù)量大于6個(gè)時(shí)，誤差離散性較小.

3.4? ?單體建筑物精細(xì)化模型精度分析

選取無人機(jī)飛行高度為80 m、建筑物外表面控制點(diǎn)數(shù)量為12個(gè)的實(shí)施方案，重構(gòu)圖書館建筑精細(xì)化三維模型，對(duì)其精度進(jìn)行分析. 建筑物外表面檢查點(diǎn)數(shù)量共86個(gè)，分布在建筑物10 m、20 m、30 m 3個(gè)高度層次，如圖12所示.如表3所示，三維模型在X軸、Y軸、高程方向的三維模型平面精度基本相同，高程精度隨著位置增高而提高，各方向各高度的均方根誤差都低于3.4 cm. 最后，計(jì)算全部檢查點(diǎn)的均方根誤差，其平面和高程標(biāo)準(zhǔn)差分別為±3.4 cm、

±1.5 cm，該精度遠(yuǎn)高于精細(xì)化建筑模型精度要求. 上述精度分析表明，采用無人機(jī)超低空、低空飛行航拍方法可以獲取高精度的地物和地貌精細(xì)化三維模型.

4? ?結(jié)? ?論

本文提出一種基于無人機(jī)影像序列的城市精細(xì)化模型構(gòu)建方法，通過無人機(jī)影像采集、影像畸變糾正、影像勻色勻光、影像區(qū)域聯(lián)合平差、三維表面重建和紋理映射等步驟實(shí)現(xiàn)模型重構(gòu)，能完全真實(shí)反映城市地形地貌，并通過大量現(xiàn)場(chǎng)無人機(jī)航拍試驗(yàn)和內(nèi)業(yè)建模工作，對(duì)構(gòu)建的精細(xì)化模型進(jìn)行了精度分析，得出以下結(jié)論：

1）當(dāng)飛行高度為80 m且控制點(diǎn)數(shù)量為12個(gè)時(shí)，

平面均方根誤差為±3.4 cm，高程均方根誤差為±3.1 cm，三維模型的平面和高程精度達(dá)到最高;當(dāng)飛行高度小于120 m時(shí)，且控制點(diǎn)數(shù)量大于6個(gè)時(shí)，誤差離散性較小.

2）通過對(duì)單體建筑物精度進(jìn)行評(píng)估，檢查點(diǎn)所在的高度越高，高程精度有所提高，平面精度基本一致. 飛行高度為80 m且控制點(diǎn)個(gè)數(shù)為12個(gè)的三維模型的平面均方根誤差為±3.4 cm，高程均方根誤差為±1.5 cm.

綜上所述，基于無人機(jī)序列影像的三維模型構(gòu)建方法可達(dá)到厘米級(jí)精度，這不僅能夠滿足城市實(shí)景模型構(gòu)建的精度要求，而且操作容易、數(shù)據(jù)采集靈活、時(shí)效性高，具有重要工程實(shí)用價(jià)值. 在今后的工作中，將從影像重疊率、拍攝傾斜角度、控制點(diǎn)空間分布、天氣狀況等精度影響因素進(jìn)行綜合分析.

參考文獻(xiàn)

[1]? ? 康志忠. LiDAR數(shù)據(jù)與數(shù)字影像的配準(zhǔn)和地物提取研究[D]. 武漢： 武漢大學(xué)遙感信息工程學(xué)院，2004：5—60.

KANG Z Z. Registration and feature extraction of LiDAR data and digital image[D]. Wuhan：School of Remote Sensing and Information Engineering，Wuhan University，2004：5—60.（In Chinese）

[2]? ? IMMERZEEL W W，KRAAIJENBRINK P D A，SHEA J M，et al.? High-resolution monitoring of Himalayan glacier dynamics using unmanned aerial vehicles[J]. Remote Sens，2014，150：93—103.

[3]? ? 王建強(qiáng)，鐘春惺，江麗鈞，等. 基于多視航空影像的城市三維建模方法[J]. 測(cè)繪科學(xué)，2014，39（3）：70—74.

WANG J Q，ZHONG C X，JIANG L J，et al. 3D city modeling construction based on multi-view aerial images[J]. Science of Surveying & Mapping，2014，39（3）：70—74.（In Chinese）

[4]? ? NEITZE F，KLONOWSK J. Mobile 3D mapping with a low-cost UAV system[J]. Remote Sensing and Spatial Information Sciences，2011，38：1—6.

[5]? ? 羅隆福，李冬，鐘杭. 基于改進(jìn)RRT的無人機(jī)電力桿塔巡檢路徑規(guī)劃[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，45（10）：80—86.

LUO L F，LI D，ZHONG H. Path planning of unmanned aircraft inspection for electric towers based on advanced RRT algorithm[J]. Journal of Hunan University（Natural Sciences），2018，45（10）：80—86.（In Chinese）

[6]? ? 李德任. 論“互聯(lián)網(wǎng)+”天基信息服務(wù)[J]. 遙感學(xué)報(bào)，2016，20（5）：708—715.

LI D R. The “Internet Plus” space-based information services[J]. Journal of Remote Sensing，2016，20（5）：708—715.（In Chinese）

[7]? ? 吳熠文，陳仁朋，余加勇，等. 無人機(jī)傾斜攝影測(cè)量技術(shù)及其工程應(yīng)用研究進(jìn)展[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，45（S1）：167—172.

WU Y W，CHEN R P，YU J Y，et al. Research and development of unmanned aerial vehicle-based oblique photogrammetry and its engineering applications[J]. Journal of Hunan University（Natural Sciences），2018，45（S1）：167—172.（In Chinese）

[8]? ? 譚金石，黃正忠. 基于傾斜攝影測(cè)量技術(shù)的實(shí)景三維建模及精度評(píng)估[J]. 現(xiàn)代測(cè)繪，2015，38（5）：21—24.

TAN J S，HUANG Z Z. Three dimensional true modeling and precision evaluation base on oblique photogrammetry technology[J]. Modern Surveying & Mapping，2015，38（5）：21—24.（In Chinese）

[9]? ? AJAYI O G，SALUBI A A，ANGBAS A F，et al. Generation of accurate digital elevation models from uav acquired low percentage overlapping images[J]. International Journal of Remote Sensing，2017，38（8）：3113—3134.

[10]? COVENEY S，ROBERTS K. Lightweight UAV digital elevation models and orthoimagery for environmental applications：data accuracy evaluation and potential for river flood risk modelling[J]. International Journal of Remote Sensing， 2017， 38（8）： 3159—3180.

[11]? 張純斌，楊勝天，趙長(zhǎng)森，等. 小型消費(fèi)級(jí)無人機(jī)地形數(shù)據(jù)精度驗(yàn)證[J]. 遙感學(xué)報(bào)，2018，22（1）：185—195.

ZHAGN C B，YANG S T，ZHAO C S，et al. Topographic data accuracy verification of small consumer UAV[J]. Journal of Remote Sensing，2018，22（1）：185—195.（In Chinese）

[12]? HARWIN S，LUCIEER A，OSBORN J. The impact of the calibration method on the accuracy of point clouds derived using unmanned aerial vehicle multi-view stereopsis[J]. Remote Sensing，2015，7（9）：11933—11953.

[13]? LUCIEER A，TURNER D，KING D H，et al. Using an Unmanned Aerial Vehicle （UAV） to capture microtopography of antarctic moss beds[J]. International Journal of Applied Earth Observation & Geoinformation，2014，27（4）：53—62.

[14]? LUCIEER A，JONG S M D，TURNER D. Mapping landslide displacements using Structure from Motion （SfM） and image correlation of multi-temporal UAV photography[J]. Progress in Physical Geography，2014，38（1）：97—116.

[15]? 徐秋輝. 無控制點(diǎn)的無人機(jī)遙感影像幾何校正與拼接方法研究[D]. 南京：南京大學(xué)地理與海洋科學(xué)學(xué)院，2013：12—56.

XU Q H. A method of geometric correction and mosaic of unmanned aerial vehicle remote sensing image without ground control points[D]. Nanjing：School of Geography and Ocean Science，Nanjing University，2013：12—56.（In Chinese）

[16]? 夏孟涵. 航空影像無控拼接及勻光勻色方法研究[D]. 武漢：武漢大學(xué)遙感信息工程學(xué)院，2017：24—68.

XIA M H. Research on pose-free stitching and color consistency correction for aerial images[D]. Wuhan：School of Remote Sensing and Information Engineering，Wuhan University，2017：24—68.（In Chinese）

[17]? 王博，張永軍，陳奇. 一種改進(jìn)的自適應(yīng)Wallis濾波算法[J]. 測(cè)繪科學(xué)，2015，40（5）：95—98.

WANG B，ZHANG Y J，CHEN Q. An improved adaptive algorithm of Wallis filter[J]. Science of Surveying and Mapping，2015，40（5）：95—98.（In Chinese）

[18]? 李德仁，鄭肇葆. 解析攝影測(cè)量學(xué)[M]. 北京：測(cè)繪出版社，1992：109—120.

LI D R，ZHENG Z B. Analytical photogrammetry[M]. Beijing：Surveying and Mapping Publishing House，1992：109—120.（In Chinese）

[19]? 萬剛. 無人機(jī)測(cè)繪技術(shù)及應(yīng)用[M]. 北京：測(cè)繪出版社，2015：206—209.

WAN G. Unmanned aerial vechicle surveying and mapping technology and application[M]. Beijing：Surveying and Mapping Publishing House，2015：206—209.（In Chinese）

[20]? 高山，陳思. 城市三維建模技術(shù)與標(biāo)準(zhǔn)研究[J]. 測(cè)繪通報(bào)，2013（3）：95—97.

GAO S，CHEN S. Research on technique and standard of 3D city modeling[J] .Bulletin of Surveying and Mapping，2013（3）：95—97.（In Chinese）