楊建強(qiáng)
(上海市政工程設(shè)計(jì)研究總院(集團(tuán))有限公司,上海市200092)
三維空間信息技術(shù)的發(fā)展非常迅速,傳統(tǒng)的二維平面圖對(duì)于三維空間的表達(dá)能力十分有限,難以滿足人們對(duì)三維空間信息的需求。鑒于此,國(guó)內(nèi)外眾多學(xué)者對(duì)于如何快速且有效地獲取空間三維信息并對(duì)其進(jìn)行三維模型的重建展開了深入的科學(xué)研究。三維模型重建常用的方法有基于攝影測(cè)量和基于三維激光掃描技術(shù)等。其中,使用三維激光掃描儀對(duì)建筑物進(jìn)行掃描獲得三維點(diǎn)云數(shù)據(jù),可以快速重建出建筑物的三維模型,但是對(duì)較高建筑物的屋頂卻無能為力。由于低空無人機(jī)具有成本低、操作靈活、受天氣影響小等特點(diǎn),尤其適合于小區(qū)域的航空影像獲取,所以在建筑物三維重建中正發(fā)揮著越來越重要的作用。使用無人機(jī)可以充分地獲取建筑物表面豐富的顏色、紋理等信息,特別是能夠拍攝到地面三維激光掃描儀掃描不到的屋頂?shù)葏^(qū)域。同時(shí)由于影像匹配算法的發(fā)展,無人機(jī)獲取到的影像可以生成建筑物帶色彩的三維點(diǎn)云,使建筑物的三維重建越來越方便與逼真。因此,研究基于無人機(jī)影像的建筑物三維模型重建方法,具有現(xiàn)實(shí)意義。
目前,國(guó)內(nèi)外對(duì)建筑物的三維重建都已經(jīng)取得了一些研究成果。邵振峰提出了一種基于航空立體影像對(duì)的人工地物三維提取和重建方法,并研究開發(fā)了一體化立體影像平臺(tái),利用該平臺(tái),可以實(shí)現(xiàn)部分較復(fù)雜房屋的三維重建[1]。翁姝提出了一種基于雷達(dá)數(shù)據(jù)和航攝圖像的區(qū)域建筑物三維重建技術(shù),但建筑物模型并不精細(xì)[2]。張龍?zhí)岢鲆环N結(jié)合自動(dòng)重建和交互建模的兩階段建模系統(tǒng),從多個(gè)視角下拍攝的建筑物圖像出發(fā),快捷有效地構(gòu)建出建筑物的三維真實(shí)感模型[3]。Pueschel等使用小型無人機(jī)拍攝的序列影像,以及地面控制點(diǎn)處拍攝的影像,經(jīng)過了影像處理、模型重建,以及紋理貼圖等步驟生成了帶紋理的且具有高分辨率的蘭登伯格城堡的三維模型[4]。
現(xiàn)以某建筑物為研究對(duì)象,提出了一種基于無人機(jī)影像的建筑物三維模型重建方法。首先使用無人機(jī)對(duì)建筑物進(jìn)行拍攝,影像匹配后獲得密集點(diǎn)云;然后使用掃描儀對(duì)無人機(jī)無法拍攝的一些區(qū)域進(jìn)行掃描,并將掃描點(diǎn)云與影像生成點(diǎn)云進(jìn)行配準(zhǔn);對(duì)點(diǎn)云預(yù)處理,再進(jìn)行Delaunay三角網(wǎng)構(gòu)建,最后使用3D Max構(gòu)建精細(xì)化模型,并對(duì)模型進(jìn)行質(zhì)量評(píng)估。
所研究的對(duì)象為某中學(xué)的一幢教學(xué)樓,表面色彩鮮艷,紋理豐富。教學(xué)樓墻面遮擋一般,可使用無人機(jī)拍攝屋頂和大部分的墻壁,對(duì)無人機(jī)無法飛行拍攝區(qū)域,采用地面激光掃描儀進(jìn)行掃描補(bǔ)充。
使用大疆無人機(jī)自帶相機(jī)對(duì)教學(xué)樓的屋頂采用航帶法連續(xù)拍攝,相機(jī)鏡頭垂直向下,與傳統(tǒng)航攝在航向和旁向重疊度的嚴(yán)格要求相比,輕型無人機(jī)的要求可適當(dāng)放寬,航向重疊度一般設(shè)置為60%~80%,最小不得小于53%;旁向重疊度一般設(shè)置為15%~60%,最小不小于8%[5]。由此拍攝到49張像素大小為4 000×3 000的屋頂序列影像,如圖1所示。然后再采用PhotoScan軟件對(duì)屋頂影像進(jìn)行密集匹配生成密集點(diǎn)云,如圖2所示。
圖1 教學(xué)樓屋頂影像圖示
圖2 教學(xué)樓屋頂密集點(diǎn)云圖示
以西立面為例,同樣采用航帶法對(duì)墻面進(jìn)行垂直拍攝,得到16張連續(xù)運(yùn)行,如圖3所示。
圖3 教學(xué)樓西立面影像圖示
使用PhotoScan軟件對(duì)墻面影像進(jìn)行密集匹配生成密集點(diǎn)云,但是由于無人機(jī)GPS定位并不精確,且無人機(jī)的位置難以使用全站儀進(jìn)行坐標(biāo)測(cè)量,因此采用外加控制點(diǎn)坐標(biāo)的方法求得墻面密集點(diǎn)云,控制點(diǎn)坐標(biāo)見表1所列。由此獲得的教學(xué)樓西墻面密集點(diǎn)云,以及控制點(diǎn)位置如圖4所示,較好地還原了墻體的特征形狀以及墻體的色彩信息。
表1 建筑物表面控制點(diǎn)坐標(biāo)表
1.3.1 墻面密集點(diǎn)云空洞修補(bǔ)
圖4 教學(xué)樓西立面密集點(diǎn)云圖示
墻面影像經(jīng)過軟件處理后獲得密集點(diǎn)云,能夠大體上反映出教學(xué)樓的形狀,但是由于遮擋,密集點(diǎn)云中還是存在一些空洞,給建筑物重建帶來一定困難,因此使用地面激光掃描點(diǎn)云來彌補(bǔ)空洞,并使用ICP算法進(jìn)行配準(zhǔn)。圖5的(紅)色框選范圍為點(diǎn)云空洞,圖6為激光掃描點(diǎn)云,圖7為配準(zhǔn)后的點(diǎn)云。
圖5 密集三維點(diǎn)云空洞圖示
1.3.2 建筑物屋頂點(diǎn)云與墻面點(diǎn)云拼接
由于無人機(jī)GPS定位不精確,因此生成的屋頂點(diǎn)云比例尺不對(duì),采取的解決方法是量測(cè)屋頂點(diǎn)云兩個(gè)角點(diǎn)A、B之間的水平距離DAB,并將其與前面配準(zhǔn)后的墻面頂部角點(diǎn)點(diǎn)云A'和B'之間的水平距離DA'B'進(jìn)行比較,確定放大的比例尺M(jìn)為1.0165,這樣就可以得到正確的屋頂點(diǎn)云。其中圖8為屋頂點(diǎn)云圖,圖9為墻面點(diǎn)云俯視圖。將建筑物屋頂點(diǎn)云和墻面點(diǎn)云同時(shí)導(dǎo)入到Geomagic軟件中,固定墻面點(diǎn)云,人工移動(dòng)屋頂點(diǎn)云,并將其與墻面點(diǎn)云精確拼接在一起,圖10為拼接后的某中學(xué)教學(xué)樓三維點(diǎn)云圖。
圖6 地面激光點(diǎn)云圖示
圖7 配準(zhǔn)后的墻面點(diǎn)云圖示
圖8 屋頂點(diǎn)云圖示
圖9 墻面點(diǎn)云俯視圖
圖10 教學(xué)樓三維點(diǎn)云圖示
經(jīng)過影像匹配生成密集點(diǎn)云,且與地面激光點(diǎn)云配準(zhǔn)后獲得的建筑物表面點(diǎn)云數(shù)據(jù)不能直接用于建筑物三維建模,還需要進(jìn)行點(diǎn)云數(shù)據(jù)預(yù)處理。包括點(diǎn)云壓縮、網(wǎng)格數(shù)據(jù)生成等[6]。點(diǎn)云壓縮方法有均勻采樣法、曲率估算法、平均點(diǎn)距值法等。本文使用Geomagic軟件提供的減少噪音和曲率法進(jìn)行點(diǎn)云數(shù)據(jù)的去噪與壓縮,曲率采用率為60%,壓縮效果如圖11所示。再采用Delaunay三角剖分法建立的墻面網(wǎng)格,如圖12所示。
圖11 點(diǎn)云數(shù)據(jù)去噪與壓縮圖示
對(duì)經(jīng)過點(diǎn)云預(yù)處理,以及三角網(wǎng)構(gòu)建后,得到了建筑物墻面的三角網(wǎng)格,同時(shí)建筑物的外形特征也能夠表現(xiàn)出來,但是這種方法得到的模型并不適用于有棱有角且外形規(guī)整的建筑物,因此使用3D Max對(duì)教學(xué)樓的模型進(jìn)行精細(xì)化構(gòu)建,實(shí)現(xiàn)了建筑物的三維模型重建。
2.2.1 墻體建模
在對(duì)建筑物的墻體進(jìn)行建模時(shí),需要先將這一部分凍結(jié),其他不需要的部位隱藏,再從其它三個(gè)視圖對(duì)其進(jìn)行建模。建模時(shí)需要開啟2.5維捕捉功能,使用“樣條線”對(duì)單墻面進(jìn)行重繪,然后再使用矩形對(duì)窗戶進(jìn)行繪制。在線的“修改器”面板下找到“附加”命令,將所有的窗戶和墻面附加在一起,運(yùn)用“擠出”命令,擠出厚度為200 mm,再在頂部視圖中將建好的模型移動(dòng)到真實(shí)的位置,對(duì)其余三面墻進(jìn)行同樣的操作,并移動(dòng)到合適的位置。
2.2.2 部件建模
建筑物部件主要包含了建筑物主體墻壁以外的門窗、雕花、斗拱、屋頂裝飾等部分。部件建模前需要開啟2.5維捕捉功能,并凍結(jié)需要構(gòu)建的部件。對(duì)門窗等規(guī)則的部件,使用“樣條線”里的矩形命令對(duì)其進(jìn)行描繪,再使用“擠出”命令,擠出門窗的厚度即可;對(duì)雕花、斗拱、屋頂裝飾而言,則需要先使用“樣條線”里的直線命令描繪,然后使用“Bezier角點(diǎn)”命令來調(diào)整線條不規(guī)則的區(qū)域,使得描繪出來的線條與實(shí)際需求盡可能一致,最后再擠出實(shí)際的厚度即可。圖13為建筑物的各個(gè)部件。
圖12 建筑物墻面點(diǎn)云封裝圖
圖13 建筑物部件圖示
建筑物的所有部件建模完成后,將所有的部件都移動(dòng)到正確的位置,生成建筑物的素模,再使用紋理貼圖功能將紋理貼到建筑物模型上,最后生成的建筑物模型如圖14所示。
圖14 教學(xué)樓樓體三維重建模型
建筑物的三維模型是利用影像或者點(diǎn)云數(shù)據(jù)進(jìn)行重建得到的模型,對(duì)于客觀存在的實(shí)體,使用儀器采集的空間數(shù)據(jù)往往是不確定的,在原始數(shù)據(jù)采集、數(shù)據(jù)預(yù)處理和模型重建階段都存在著影響模型精度的因素,這往往導(dǎo)致了重建后模型精度的損失。因此,任何重建的建筑物三維模型實(shí)際上只能是對(duì)建筑物實(shí)際形態(tài)的近似模擬,不可能完全真實(shí)反映實(shí)際情況,重建的建筑物模型需要進(jìn)行精度評(píng)價(jià)。為此,采用距離法進(jìn)行精度評(píng)價(jià),在教學(xué)樓的西立面,通過在重建后的墻面上選擇了8段特征線段,量測(cè)它們的距離。然后,使用全站儀在墻面上測(cè)量相同的點(diǎn)位坐標(biāo),計(jì)算其長(zhǎng)度,將其與特征線段的長(zhǎng)度比較,借以評(píng)定重建的精度。墻面特征點(diǎn)如圖15所示,特征線段的距離比較見表2所列。
圖15 墻面測(cè)距選點(diǎn)點(diǎn)位圖
由表2可知特征線的實(shí)測(cè)長(zhǎng)度與模型量測(cè)長(zhǎng)度存在一定的偏差,但偏差較小,最大偏差值為1.5 cm,平均誤差為10.6 mm,從距離比較結(jié)果來看,整體長(zhǎng)度值偏差都在1.5 cm以內(nèi),能滿足一般的建筑物三維建模的要求。
表2 特征線段距離表
本文以無人機(jī)影像為主要數(shù)據(jù)來源,以激光掃描點(diǎn)云數(shù)據(jù)作為補(bǔ)充,經(jīng)過影像匹配、密集點(diǎn)云生成、點(diǎn)云配準(zhǔn),以及使用3D Max建模等步驟實(shí)現(xiàn)了教學(xué)樓的精細(xì)化建模,并對(duì)其進(jìn)行了質(zhì)量評(píng)價(jià),結(jié)果表明三維重建質(zhì)量可靠,能夠滿足一般建筑物的建模要求。