婁寧，馬健，楊永崇，張全文

(西安科技大學(xué) 測(cè)繪科學(xué)與技術(shù)學(xué)院，西安 710054)

0 引言

無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)的迅速發(fā)展大大推動(dòng)了實(shí)景三維建模技術(shù)的進(jìn)步[1-3]，這些實(shí)景三維模型數(shù)據(jù)作為單體建筑物仿真、保護(hù)、研究的基礎(chǔ)數(shù)據(jù)與依據(jù)，具有巨大的市場(chǎng)價(jià)值與應(yīng)用前景[4-6]。由于它們本身復(fù)雜的外形結(jié)構(gòu)和豐富的紋理類型，導(dǎo)致單體建筑物精細(xì)化實(shí)景三維建模的難度較高[7]，快速、高效、完整地收集建筑物高分辨率影像數(shù)據(jù)是實(shí)景三維建模的關(guān)鍵問題。在實(shí)際生產(chǎn)中對(duì)單體建筑物主要有以下3種三維建模方法。第1種是傳統(tǒng)的三維模型生產(chǎn)制作方法，利用3DS MAX、AutoCAD等軟件手工繪制建筑物的三維模型，該方法效率較低、模擬紋理的真實(shí)度不高[8]；第2種是利用三維激光掃描儀(LiDAR)獲取地物三維點(diǎn)云數(shù)據(jù)，運(yùn)用紋理映射等方法快速建立三維模型[9]，但該方法對(duì)建筑物頂部或某些側(cè)面部位數(shù)據(jù)獲取難度較大；第3種是無(wú)人機(jī)低空攝影測(cè)量通過(guò)搭載成像與非成像傳感器，利用航帶網(wǎng)式航線或簡(jiǎn)單的環(huán)繞飛行航線采集地物信息，建立實(shí)景三維模型[10-12]，但該方法并沒有依據(jù)建筑物實(shí)際情況計(jì)算飛行航線關(guān)鍵參數(shù)，設(shè)置針對(duì)不同建筑物的三維環(huán)繞式飛行航線，受主觀因素影響較大。因而，目前仍缺乏一種高效的單體建筑物精細(xì)化實(shí)景三維建模方法。

本文根據(jù)單體建筑物的自身?xiàng)l件、影像采集需求和單鏡頭多旋翼無(wú)人機(jī)的特點(diǎn)，計(jì)算飛行航線的關(guān)鍵參數(shù)，對(duì)不同單體建筑物的傾斜影像和正射影像采集方法進(jìn)行航線自動(dòng)規(guī)劃設(shè)計(jì)，通過(guò)三維環(huán)繞拍攝[13-16]，在低空完整收集建筑物高分辨率傾斜影像數(shù)據(jù)和正射影像數(shù)據(jù)，并按照Smart3D Capture軟件實(shí)景三維建模的相關(guān)技術(shù)流程，將傾斜影像數(shù)據(jù)和正射影像數(shù)據(jù)相融合，進(jìn)行實(shí)景三維重建，并從完整性、紋理細(xì)節(jié)、量測(cè)精度3個(gè)方面對(duì)實(shí)景三維模型進(jìn)行分析。

1 單體建筑物三維環(huán)繞式航線設(shè)計(jì)

基于單體建筑物進(jìn)行精細(xì)化實(shí)景三維建模，需要盡可能獲取建筑物外側(cè)全部的紋理信息，根據(jù)本文提出的三維環(huán)繞式飛行方案，對(duì)于單鏡頭多旋翼無(wú)人機(jī)來(lái)說(shuō)，需要采用高效合理的傾斜式三維環(huán)繞航線和正射式三維環(huán)繞航線，對(duì)單體建筑物進(jìn)行均勻環(huán)拍，以完整采集建筑物高分辨影像數(shù)據(jù)。

1.1 基本理論

首先，通過(guò)地形圖和衛(wèi)星影像等參考資料對(duì)測(cè)區(qū)地物進(jìn)行詳細(xì)分析，并通過(guò)實(shí)地勘察全方位了解建筑物的長(zhǎng)a1、寬b1、高c1，以此為基礎(chǔ)建立建筑物的外接長(zhǎng)方體，并求得外接長(zhǎng)方體的長(zhǎng)a2、寬b2、高c2，然后建立外接長(zhǎng)方體的外接半球面，如圖1(a)所示。最小圓柱體包圍盒，如圖1(b)所示。其半徑分別為R0與L0。計(jì)算方法如式(1)、式(2)所示。

(1)

(2)

圖1 外接半球面、傾斜航路半球面、最小圓柱體包圍盒、正射航路圓柱體

以外接半球面和最小圓柱體包圍盒為攝影基準(zhǔn)面，根據(jù)衡量影像分辨率的重要指標(biāo)地面分辨率(ground sample distance，GSD)，以正射影像的地面采樣距離GSD為準(zhǔn)計(jì)算攝影物距D0，如式(3)所示。

(3)

式中：pl為長(zhǎng)邊像元數(shù)量，即相機(jī)感光元件參數(shù)；FOV為鏡頭視角，鏡頭固定參數(shù)。

傾斜式航路半球面的半徑R可以由外接半球面的半徑、攝影物距經(jīng)過(guò)簡(jiǎn)單的幾何推算得到，如圖1(a)所示。正射式航路圓柱體的半徑L可由最小圓柱體包圍盒的半徑、攝影物距得到，如圖1(b)所示。其計(jì)算如式(4)、式(5)所示。

R=R0+D0

(4)

L=L0+D0

(5)

1.2 傾斜航線參數(shù)計(jì)算

依據(jù)最小外接半球面的半徑、航路半球面的半徑、航向重疊度、旁向重疊度、航向像元數(shù)量、旁向像元數(shù)量、長(zhǎng)邊像元數(shù)量、最大飛行速度，計(jì)算每條環(huán)繞航路的相機(jī)俯仰角、航高、環(huán)繞半徑、攝影基線長(zhǎng)度、攝影間隔角、攝影間隔時(shí)間、任務(wù)執(zhí)行時(shí)間。

(6)

初始航高H1、環(huán)繞半徑R1、俯仰角α1之間的關(guān)系如式(7)、式(8)所示。

H1=R×sinα1

(7)

R1=R×cosα1

(8)

無(wú)人機(jī)在飛行方向上相鄰影像重疊的比例航向重疊率Ox，用于計(jì)算航向攝影基線長(zhǎng)度B1，計(jì)算方法如式(9)所示。

(9)

式中：Px為航向像元數(shù)量，相機(jī)像素與航線平行方向的像元數(shù)量。

根據(jù)當(dāng)前航高和相機(jī)的水平視場(chǎng)角推算攝影間隔角β1，其計(jì)算如式(10)所示。

(10)

當(dāng)無(wú)人機(jī)采用等時(shí)間間隔的方法采集影像時(shí)，可由攝影基線長(zhǎng)度和無(wú)人機(jī)允許范圍內(nèi)的航速比值得到攝影間隔時(shí)間t，其計(jì)算如式(11)所示。

(11)

式中：Vmax為無(wú)人機(jī)最大飛行速度。

在進(jìn)行航線規(guī)劃時(shí)將航線長(zhǎng)度除以容許無(wú)人機(jī)最大飛行速度就可以得到任務(wù)執(zhí)行時(shí)間參數(shù)T，其計(jì)算如式(12)所示。

(12)

2)當(dāng)前航線參數(shù)計(jì)算。根據(jù)上一條環(huán)繞航線的航高和旁向重疊度計(jì)算當(dāng)前環(huán)繞航路相機(jī)中軸線的俯仰角α2，計(jì)算方法如式(13)所示。

αn=αn-1+α1(1-Oy)

(13)

式中：αn-1為上一條環(huán)繞航路的相機(jī)俯仰角；Oy為旁向重疊度。

當(dāng)前環(huán)繞航線的航高、環(huán)繞半徑、攝影基線長(zhǎng)度、攝影間隔角、攝影間隔時(shí)間、任務(wù)執(zhí)行時(shí)間可通過(guò)式(7)～式(12)進(jìn)行循環(huán)計(jì)算，若當(dāng)前環(huán)繞航路滿足式(14)～式(16)則終止計(jì)算，說(shuō)明當(dāng)前航線下的相機(jī)視場(chǎng)已完全覆蓋建筑物的頂部，否則繼續(xù)執(zhí)行式(6)的后續(xù)步驟。

Hn>C2

(14)

(15)

(16)

1.3 正射航線參數(shù)計(jì)算

依據(jù)建筑物最小圓柱體包圍盒的半徑、航路圓柱體的半徑、重疊度、像元數(shù)量等已知參數(shù)，計(jì)算三維環(huán)繞航線的俯仰角、航高、環(huán)繞半徑、攝影基線長(zhǎng)度、攝影間隔角、攝影間隔時(shí)間、任務(wù)執(zhí)行時(shí)間等參數(shù)。

1)初始航線參數(shù)計(jì)算。正射式三維環(huán)繞航線的初始云臺(tái)俯仰角為α1，初始航高為H1，環(huán)繞半徑為R1，攝影基線長(zhǎng)度為B1，攝影間隔角為β1。

α1=0°

(17)

(18)

R1=L

(19)

B1=(1-Ox)×GSD×Px

(20)

(21)

攝影間隔時(shí)間、任務(wù)執(zhí)行時(shí)間通過(guò)式(11)～式(12)進(jìn)行計(jì)算。

2)當(dāng)前航線參數(shù)計(jì)算。根據(jù)上一條環(huán)繞航路的航高和旁向重疊度計(jì)算當(dāng)前環(huán)繞航路的航高，計(jì)算方法如式(22)所示。

Hn=Hn-1+2H1×(1-Oy)

(22)

當(dāng)前環(huán)繞航線的云臺(tái)俯仰角重復(fù)式(17)，環(huán)繞半徑、攝影基線長(zhǎng)度和間隔角循環(huán)式(19)～式(21)，攝影間隔時(shí)間、任務(wù)執(zhí)行時(shí)間通過(guò)式(11)～式(12)進(jìn)行計(jì)算，直到當(dāng)前航線的航高滿足式(23)時(shí)終止計(jì)算。說(shuō)明當(dāng)前航線下的相機(jī)視場(chǎng)角已完全拍攝到建筑物側(cè)面頂部的影像數(shù)據(jù)。

Hn>C2-H1

(23)

2 環(huán)繞航路數(shù)學(xué)模型的設(shè)計(jì)與實(shí)現(xiàn)

在以往的無(wú)人機(jī)傾斜攝影環(huán)繞飛行過(guò)程中，由于沒有專門針對(duì)單體建筑精細(xì)化建模影像數(shù)據(jù)采集的三維環(huán)繞航線規(guī)劃軟件，要求飛手根據(jù)航空攝影測(cè)量的基礎(chǔ)知識(shí)和以往的飛行經(jīng)驗(yàn)人為控制無(wú)人機(jī)的飛行航線參數(shù)。在進(jìn)行航攝的過(guò)程中，有很多因素導(dǎo)致影像數(shù)據(jù)質(zhì)量下降。若是能根據(jù)建筑物的自身?xiàng)l件，通過(guò)數(shù)學(xué)模型的方式自動(dòng)規(guī)劃出針對(duì)不同種類單體建筑物的三維環(huán)繞式航線，以確立環(huán)繞立體采集的最優(yōu)方式，使得無(wú)人機(jī)航攝的外業(yè)流程實(shí)現(xiàn)自動(dòng)化或半自動(dòng)化，這對(duì)于影像數(shù)據(jù)的航攝質(zhì)量和航攝效率會(huì)具有較大的意義。

本文利用C#軟件產(chǎn)品設(shè)計(jì)流程，對(duì)三維環(huán)繞航線的計(jì)算系統(tǒng)進(jìn)行需求分析和功能分析，并根據(jù)這些分析建立了數(shù)學(xué)模型，數(shù)學(xué)模型實(shí)現(xiàn)的技術(shù)路線如圖2所示。模型建立方法如下：首先根據(jù)1.1節(jié)的步驟，構(gòu)造傾斜航路半球面和正射航路圓柱體；然后根據(jù)1.2節(jié)的步驟，依據(jù)已知參數(shù)和終止條件，計(jì)算傾斜式三維環(huán)繞航路飛行參數(shù)；其次根據(jù)1.3節(jié)的步驟，依據(jù)航路圓柱體的半徑、高度及其終止條件，計(jì)算正射航路的飛行參數(shù)；最后在C#編程語(yǔ)言中將三維環(huán)繞航路數(shù)學(xué)模型的代碼進(jìn)行編譯和運(yùn)行測(cè)試。如圖3所示，其顯示界面主要分為2個(gè)部分：輸入?yún)?shù)部分包括單體建筑物幾何參數(shù)和實(shí)景三維模型的質(zhì)量要求指標(biāo)參數(shù)；輸出參數(shù)部分用于顯示傾斜攝影三維環(huán)繞航路和正射三維環(huán)繞航路的飛行參數(shù)。

圖2 三維環(huán)繞航路數(shù)學(xué)模型實(shí)現(xiàn)技術(shù)路線

3 建模實(shí)例及實(shí)驗(yàn)分析

3.1 實(shí)驗(yàn)對(duì)象和設(shè)備

實(shí)驗(yàn)挑選大疆精靈4(Phantom4)無(wú)人機(jī)為傾斜攝影航攝平臺(tái)，相片分辨為3 648×4 864，鏡頭視角為84°。實(shí)景三維建模的實(shí)驗(yàn)對(duì)象是一棟單體建筑物。該建筑的幾何參數(shù)為長(zhǎng)6.6 m、寬3.6 m、高4 m。采用本文設(shè)計(jì)的環(huán)繞航路數(shù)學(xué)模型自動(dòng)計(jì)算出無(wú)人機(jī)傾斜攝影三維環(huán)繞航路和正射三維環(huán)繞航路的飛行參數(shù)。軟件采用Smart3D Capture實(shí)景建

模大師，它具有強(qiáng)大的數(shù)據(jù)運(yùn)算能力，能還原出接近真實(shí)的毫米級(jí)模型。

3.2 數(shù)據(jù)采集及處理

利用本文設(shè)計(jì)的三維環(huán)繞航路數(shù)學(xué)模型，輸入建筑物的幾何參數(shù)長(zhǎng)寬高、相機(jī)視場(chǎng)角、像元數(shù)量、飛行速度、影像空間分辨率0.1 cm、航向重疊度80%、旁向重疊度80%，輸出傾斜攝影和正射三維環(huán)繞飛行的參數(shù)(圖3)。無(wú)人機(jī)則按照上述參數(shù)生成的三維環(huán)繞航線采集建筑物的傾斜影像數(shù)據(jù)和正射影像數(shù)據(jù)。按照Smart3D Capture軟件實(shí)景建模的相關(guān)技術(shù)流程，將傾斜影像數(shù)據(jù)和正射影像數(shù)據(jù)融合處理，依據(jù)影像的內(nèi)、外方位元素獲取多視角影像的同名點(diǎn)坐標(biāo)，進(jìn)行影像密集匹配[17]，生成建筑物的密集點(diǎn)云并構(gòu)建不規(guī)則三角網(wǎng)TIN，然后進(jìn)行紋理映射[18]，最終生成實(shí)景三維模型，如圖4所示。

圖3 三維環(huán)繞航路飛行參數(shù)圖

圖4 相機(jī)位置連接點(diǎn)、密集點(diǎn)云和實(shí)景三維模型

3.3 實(shí)驗(yàn)分析

為檢驗(yàn)本文設(shè)計(jì)的三維環(huán)繞航路采集的影像數(shù)據(jù)所建單體精細(xì)化實(shí)景三維模型的可行性，對(duì)某單體建筑物進(jìn)行了實(shí)例建模，并從實(shí)景三維模型數(shù)據(jù)的完整性、紋理細(xì)節(jié)及量測(cè)精度3個(gè)方面進(jìn)行了質(zhì)量分析。

三維實(shí)景模型的完整性可由模型紋理不存在缺失和變形的面積與三維模型表面積的比值作為評(píng)價(jià)模型完整率的標(biāo)準(zhǔn)。經(jīng)計(jì)算，無(wú)人機(jī)單相機(jī)建模的完整率為0.91，具有較優(yōu)的完整性。在紋理細(xì)節(jié)方面，從圖4(c)可知，實(shí)景三維模型的細(xì)節(jié)基本沒有紋理缺失，色調(diào)均勻且能夠真實(shí)地反映建筑物的表面信息，具有較完整且清晰真實(shí)的紋理細(xì)節(jié)。由于實(shí)景三維建模實(shí)現(xiàn)了傾斜加正射的三維環(huán)繞式拍攝，該結(jié)果有效減少了觀測(cè)盲區(qū)且完整采集了建筑物的高分辨率影像數(shù)據(jù)。

對(duì)單體建筑物的仿真、保護(hù)和研究工作，更注重的是建筑本身的尺寸、結(jié)構(gòu)和紋理信息，因而在三維重建的精度更應(yīng)該關(guān)注相對(duì)精度。本文對(duì)模型的相對(duì)精度采用檢查點(diǎn)外業(yè)實(shí)測(cè)值與內(nèi)業(yè)模型量測(cè)值對(duì)比的方法進(jìn)行評(píng)估，將外業(yè)實(shí)測(cè)的10個(gè)距離和高差與內(nèi)業(yè)直接在實(shí)景三維模型上的量測(cè)值進(jìn)行對(duì)比檢驗(yàn)。圖5(a)顯示了距離量測(cè)值與實(shí)測(cè)值的差值，圖5(b)顯示了高差量測(cè)值與實(shí)測(cè)值的差值。從表1的統(tǒng)計(jì)結(jié)果可知，點(diǎn)與點(diǎn)之間的距離精度、其模型量測(cè)值與外業(yè)實(shí)測(cè)值的差值范圍是(0.002 4±0.001 5) m；點(diǎn)與點(diǎn)之間的高差精度、其模型量測(cè)值與外業(yè)實(shí)測(cè)值

圖5 三維模型量測(cè)精度分析

的差值范圍是(0.002 7±0.002 8) m。本實(shí)驗(yàn)的分析結(jié)果表明，基于本文設(shè)計(jì)的三維環(huán)繞式航線采集的影像數(shù)據(jù)所制作的實(shí)景三維模型尺寸精度滿足規(guī)范要求，所建模型有完整性好、分辨率高、紋理真實(shí)、幾何精度可靠的特性。

表1 實(shí)景三維模型精度統(tǒng)計(jì) m

4 結(jié)束語(yǔ)

本文針對(duì)單體建筑物真三維精細(xì)化建模問題，提出了以單鏡頭多旋翼無(wú)人機(jī)為載體的三維環(huán)繞式航線自動(dòng)規(guī)劃方法。該方法針對(duì)單體建筑物的不同特征進(jìn)行分析，自動(dòng)設(shè)計(jì)不同的三維環(huán)繞式無(wú)人機(jī)飛行方案，提升了單體建筑高分辨率影像數(shù)據(jù)采集的完整性、質(zhì)量和效率，所建模型具有完整性好、紋理真實(shí)、幾何精度可靠的特性。同時(shí)，本文所述方法對(duì)于單體建筑物中“回”字型等特殊形態(tài)的單體建筑內(nèi)部精細(xì)化建模存在局限性。該情況將作為本文未來(lái)進(jìn)一步研究的方向，并尋求“回”字型等單體建筑物精細(xì)化建模的最優(yōu)方式。