王云川，段平，李佳，王明果，姚永祥

(1.云南師范大學(xué) 旅游與地理科學(xué)學(xué)院，云南 昆明 650500;2.云南省高校資源與環(huán)境遙感重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650500;3.云南省地理空間信息工程技術(shù)研究中心，云南 昆明 650500;4.云南省地質(zhì)科學(xué)研究所，云南 昆明 650501)

1 引 言

近年來(lái)，隨著多旋翼無(wú)人機(jī)的興起，利用無(wú)人機(jī)影像進(jìn)行三維建模越來(lái)越廣泛地應(yīng)用于水利建設(shè)[1]、消防救援[2]及智慧城市[3-5]等領(lǐng)域.多旋翼無(wú)人機(jī)具有成本低、可操作性高以及獲取的影像空間分辨率高等優(yōu)點(diǎn)，通過(guò)云臺(tái)相機(jī)與航線的配合，實(shí)現(xiàn)傾斜攝影測(cè)量的任務(wù).同時(shí)，目前已有許多技術(shù)成熟的三維建模軟件，如Context capture、Pix4d以及Agisoft Metashape等，可將正射影像與傾斜影像聯(lián)合平差生成點(diǎn)云以及實(shí)景三維模型.然而，受鏡頭性能的限制以及拍攝環(huán)境的影響，無(wú)人機(jī)在同一區(qū)域內(nèi)采集的影像在色調(diào)方面存在一定差異，若不經(jīng)過(guò)處理而直接進(jìn)行三維重建，必然會(huì)對(duì)實(shí)景三維模型的質(zhì)量造成影響，甚至導(dǎo)致建模失敗,因此對(duì)色差過(guò)大的影像進(jìn)行色調(diào)勻化就非常必要.

常用的經(jīng)典勻色算法有Mask勻光法、Wallis勻光法以及直方圖均衡化(HE，Histogram Equalization)等.Mask勻光算法又稱模糊正像勻光法，它將光照不均勻的影像看作是由理想狀態(tài)下光照分布均勻的影像和背景影像構(gòu)成，因此可通過(guò)從原始影像中減去模擬出來(lái)的背景影像，從而消除影像的色調(diào)不均勻現(xiàn)象[6].Mask勻光算法適應(yīng)性強(qiáng)，應(yīng)用較廣泛，但存在偏色、光暈效應(yīng)以及運(yùn)算效率低等不足.Wallis算法是將影像的灰度均值及方差映射到相應(yīng)定值，并且使不同影像區(qū)域的灰度方差和灰度均值都近似相等，使影像反差小的區(qū)域反差增大，使反差大的區(qū)域反差減小，從而使得整幅影像的亮度方差達(dá)到一致[7]，然而Wallis算法在處理單幅影像色調(diào)不均勻現(xiàn)象時(shí)，效果往往并不明顯，并且容易產(chǎn)生“塊效應(yīng)”.直方圖均衡化則是一種適應(yīng)性及效果都較好的圖像增強(qiáng)方法，其本質(zhì)是通過(guò)一種變換函數(shù)，使原始圖像以目標(biāo)直方圖進(jìn)行變換，通過(guò)改變映射函數(shù)來(lái)控制變換結(jié)果的精度和效果，對(duì)對(duì)比度低的圖像有較好的處理效果[8].直方圖均衡化雖然使圖像色階分布更加均勻，但實(shí)際上很難達(dá)到格式所允許的最大灰度變化范圍，因此對(duì)圖像進(jìn)行直方圖均衡化后，使用灰度拉伸使圖像的色階進(jìn)一步擴(kuò)展至整個(gè)灰度級(jí).灰度拉伸又稱為對(duì)比度拉伸，是一種簡(jiǎn)單的線性點(diǎn)運(yùn)算，灰度拉伸可以更加精確地控制直方圖的變化范圍，可以有選擇地拉伸灰度區(qū)間以改善圖像.因此對(duì)整體色調(diào)較暗及色階局部分布不規(guī)則不均勻的無(wú)人機(jī)影像，使用直方圖均衡化及灰度拉伸方法能將原影像灰度拉伸到整個(gè)灰度級(jí)并且均勻分布，從而提高圖像整體對(duì)比度[9].

本文探討了基于直方圖均衡化的無(wú)人機(jī)影像色調(diào)勻化處理與三維建模方法.以云南師范大學(xué)某不規(guī)則建筑物為測(cè)區(qū)，利用無(wú)人機(jī)設(shè)備采集該建筑的正射以及傾斜影像，對(duì)該組影像進(jìn)行色調(diào)勻化處理，以勻色前及勻色后的影像為基礎(chǔ)，分別構(gòu)建實(shí)景三維模型，并對(duì)其三維模型質(zhì)量進(jìn)行對(duì)比分析.

2 無(wú)人機(jī)影像三維建模方法

無(wú)人機(jī)影像實(shí)景三維建模一般過(guò)程包括數(shù)據(jù)采集及預(yù)處理、空中三角測(cè)量、密集匹配、數(shù)字表面模型(digital surface model，DSM)構(gòu)建以及紋理映射等，其中空中三角測(cè)量作為三維模型構(gòu)建的關(guān)鍵步驟，往往會(huì)因?yàn)橛跋裆{(diào)過(guò)暗或色調(diào)差異過(guò)大而失敗，導(dǎo)致特征點(diǎn)的識(shí)別以及點(diǎn)云的構(gòu)建效果較差，因此對(duì)影像進(jìn)行色調(diào)勻化處理極為重要.結(jié)合三維建模流程，將其劃分為三個(gè)環(huán)節(jié)，即：數(shù)據(jù)采集及勻化處理、三維建模及質(zhì)量對(duì)比分析，總體流程如圖1所示.數(shù)據(jù)采集及勻化處理環(huán)節(jié)：通過(guò)無(wú)人機(jī)實(shí)拍獲取影像，以直方圖均衡化及灰度拉伸的方法對(duì)影像進(jìn)行色調(diào)勻化處理；三維建模環(huán)節(jié)：借助POS數(shù)據(jù)進(jìn)行空中三角測(cè)量，獲取影像的地面點(diǎn)坐標(biāo)，生成稀疏點(diǎn)云，隨后進(jìn)行密集匹配，構(gòu)建DSM模型，最后映射紋理生成實(shí)景三維模型；質(zhì)量對(duì)比分析環(huán)節(jié)：對(duì)生成的兩組實(shí)景三維模型從模型完整性以及幾何紋理細(xì)節(jié)兩個(gè)方面進(jìn)行對(duì)比分析.

圖1 勻化處理的無(wú)人機(jī)影像三維建模總體流程

2.1 數(shù)據(jù)獲取

選擇云南師范大學(xué)武之樓為采集實(shí)驗(yàn)區(qū)，該建筑結(jié)構(gòu)規(guī)整，四周植被茂盛.以大疆“DJI Phantom 4 Pro”無(wú)人機(jī)為數(shù)據(jù)采集設(shè)備，設(shè)置曝光間隔為2 s，相對(duì)航高為60 m，航向重疊度和旁向重疊度均設(shè)定為80%，影像分辨率為2 cm，獲取影像418張，共耗時(shí)15 min.拍攝期間天氣晴朗，光照較強(qiáng)，采集影像色調(diào)差異明顯.

2.2 影像色調(diào)勻化關(guān)鍵技術(shù)

影像拍攝時(shí)由于光線太強(qiáng)，使圖像對(duì)比度減弱，細(xì)節(jié)分辨不清，這樣的圖像直方圖灰度往往集中在某一色階范圍之內(nèi),此時(shí)需要將這些灰度拉伸到整個(gè)灰度級(jí)上，并使它們?cè)谥狈綀D中均勻分布，以達(dá)到色調(diào)勻化的目的.直方圖均衡化是灰度變換的一個(gè)重要應(yīng)用，是一種簡(jiǎn)單有效的圖像增強(qiáng)方法，因其運(yùn)算簡(jiǎn)單且能夠有效地提高圖像的整體對(duì)比度而被廣泛應(yīng)用.

2.2.1 直方圖均衡化

直方圖均衡化的中心思想就是把原始圖像的灰度直方圖從某個(gè)比較集中的區(qū)間變換成在更大區(qū)間內(nèi)均勻分布的形式，以達(dá)到增強(qiáng)圖像整體對(duì)比度的效果.變換公式如下：

s=T(r)

(1)

式中：s為經(jīng)直方圖均衡化后的圖像灰度級(jí)，r為原圖像的灰度級(jí)，T(r)為變換函數(shù).

直方圖均衡化就是根據(jù)直方圖對(duì)原像素點(diǎn)經(jīng)變換函數(shù)T(r)產(chǎn)生一個(gè)新的像素點(diǎn)，其中變換函數(shù)T(r)應(yīng)滿足以下兩個(gè)條件：

(1)在0≤r≤1時(shí)，T(r)為單調(diào)遞增函數(shù)；

(2)在0≤r≤1時(shí)，0≤T(r)≤1.

條件(1)保證均衡化后圖像的灰度級(jí)依然保持從黑到白的排列次序，條件(2)保證均衡化后圖像的灰度級(jí)在允許的范圍內(nèi)[10].綜合以上兩個(gè)條件，因?yàn)槔鄯e分布函數(shù)(CDF，Cumulative Distribution Function)是單調(diào)增函數(shù)，并且值域在0-1，因此直方圖均衡化中采用累積分布函數(shù)作為變換函數(shù)，即：

(2)

式中：sk為均衡化后的像素灰度級(jí)，rk為原元素灰度級(jí)，nj表示當(dāng)前灰度級(jí)的像素?cái)?shù)量，N表示圖像像素總數(shù).

直方圖均衡化一般首先統(tǒng)計(jì)原圖像的灰度直方圖，并計(jì)算出像素?cái)?shù)量、灰度分布頻率和灰度累積分布頻率等參數(shù)；然后根據(jù)變換函數(shù)(2)計(jì)算灰度映射表；最后根據(jù)映射表通過(guò)映射關(guān)系生成具有新灰度的圖像.

2.2.2 灰度拉伸

經(jīng)過(guò)直方圖均衡化之后，圖像的灰度級(jí)的分布范圍雖然更大更均勻，但可能依然沒(méi)有填充至整個(gè)區(qū)間，因此通過(guò)灰度拉伸使其充滿整個(gè)灰度級(jí)范圍，以達(dá)到對(duì)比度大幅增強(qiáng)的目的.灰度拉伸公式如下：

(3)

式中：ck為拉伸后的像素灰度級(jí)，sk為拉伸前的像素灰度級(jí)，smax和smin分別為拉伸前圖像的最大灰度值和最小灰度值，cmax和cmin分別為要拉伸到的灰度空間的最大灰度值和最小灰度值.

灰度拉伸一般首先統(tǒng)計(jì)圖像的灰度直方圖，并確定拉伸前的灰度值以及想要拉伸到的灰度空間的最大和最小灰度值等參數(shù)；然后根據(jù)拉伸公式(3)計(jì)算灰度映射表；最后根據(jù)映射表通過(guò)映射關(guān)系生成具有新灰度的圖像.

2.2.3 實(shí)驗(yàn)結(jié)果

從圖2勻色前后直方圖的對(duì)比中可以發(fā)現(xiàn)，勻色之前影像的直方圖中，灰度值在前半部分(具體灰度值范圍為6-30以及73-115處)有兩個(gè)小高峰，尤其是前一個(gè)的高峰峰值較高，在0-128范圍內(nèi)的灰度值的百分位達(dá)到了70%，而在129-255范圍內(nèi)的灰度值百分位僅有30%，這使得影像整體上偏暗；勻色處理后影像的直方圖變得較為平緩，影像的灰度值被均勻地分為3份，每一份所占百分比都近似33%，在0-33灰度范圍內(nèi)仍有部分峰值較高的灰度值，說(shuō)明影像有些部分仍然較暗，但整體上已有較好地改觀.

圖2 影像勻色前(左)及勻色后(右)直方圖對(duì)比

從圖3影像本身的對(duì)比中可以看到，勻色后影像整體色調(diào)偏亮.在細(xì)節(jié)處，仔細(xì)觀察影像的左上角，如圖4所示，路旁有一條白色的砌磚，在原圖上由于樹(shù)影的遮擋幾乎不可見(jiàn)，而勻色后該砌磚清晰可見(jiàn)；在影像的右側(cè)，如圖5所示，建筑物的右側(cè)植被被建筑物的影子所遮擋，在原圖上幾乎是整片漆黑，看不到任何可用信息，經(jīng)勻色處理后，雖然整體色調(diào)較暗，但已經(jīng)可以識(shí)別出一定的植被信息.

圖3 影像勻色前(左)及勻色后(右)對(duì)比

圖4 影像細(xì)節(jié)處勻色前(左)及勻色后(右)對(duì)比

圖5 影像細(xì)節(jié)處勻色前(左)及勻色后(右)對(duì)比

2.3 實(shí)景三維建模方法

空中三角測(cè)量是無(wú)人機(jī)實(shí)景三維建模的關(guān)鍵步驟，將攝影瞬間POS系統(tǒng)記錄的值作為影像的初始外方位元素，并結(jié)合少量的實(shí)測(cè)控制點(diǎn)，采用光束法進(jìn)行區(qū)域網(wǎng)整體平差解算，實(shí)現(xiàn)多視角聯(lián)合空中三角測(cè)量，并將通過(guò)影像匹配得到的連接點(diǎn)及部分地面控制點(diǎn)納入已知控制點(diǎn)的大地坐標(biāo)系中，來(lái)獲得每張影像的外方位元素和加密點(diǎn)的大地坐標(biāo)，以此生成稀疏點(diǎn)云.然后，對(duì)每個(gè)立體像對(duì)進(jìn)行密集匹配，從而形成密集點(diǎn)云.通過(guò)點(diǎn)云生成不規(guī)則三角網(wǎng)(Triangulated Irregular Network ，TIN)構(gòu)建數(shù)字表面模型(Digital Surface Model，DSM).最后，對(duì)三維模型和紋理影像配準(zhǔn)，建立空間地物點(diǎn)到各影像的投影關(guān)系，篩選出效果最優(yōu)的目標(biāo)影像，將其反投影到模型的三角面片上，實(shí)現(xiàn)三維模型的紋理貼合，最終形成實(shí)景三維模型[11].

使用Context Capture軟件重建該建筑的三維模型，依據(jù)上述三維模型重建方法，由航攝影像生成點(diǎn)云，之后生成點(diǎn)云的TIN來(lái)構(gòu)建DSM，最后進(jìn)行紋理映射，最終生成實(shí)景三維模型.

3 三維模型重建質(zhì)量對(duì)比分析

在模型中選取具有代表性的建筑紋理區(qū)域，主要通過(guò)模型的完整性和幾何紋理細(xì)節(jié)來(lái)對(duì)比勻色前后影像的建模效果.從圖6建筑物正立面的對(duì)比中發(fā)現(xiàn)，由于拍攝時(shí)陽(yáng)光從東北面照射且強(qiáng)度較大，使屋頂和正立面的色差較大，這一點(diǎn)在未經(jīng)勻色處理的三維模型中表現(xiàn)得比較明顯，整個(gè)正立面幾乎被黑影籠罩，紋理細(xì)節(jié)難以分辨.反觀經(jīng)勻色處理后的三維模型，整體色調(diào)顯得更為均勻，屋頂和正立面之間的過(guò)渡平緩舒適，整體性較強(qiáng)，紋理更加清晰.

圖6 三維模型勻色前(左)及勻色后(右)對(duì)比

圖7為建筑物的俯視圖，由于建筑物較為高大，導(dǎo)致建筑物中部區(qū)域都被墻體的影子遮擋，在未經(jīng)勻色處理的建筑物三維模型中(圖7左)，不僅存在色調(diào)太暗而難以分辨紋理細(xì)節(jié)的區(qū)域，還存在較多的空洞，這是因?yàn)樵加跋襁^(guò)暗，三維重建時(shí)難以識(shí)別同名點(diǎn)以及點(diǎn)云的構(gòu)建所致.

圖7 三維模型勻色前(左)及勻色后(右)對(duì)比

放大對(duì)比圖7中的區(qū)域1及區(qū)域2，得圖8及圖9.從區(qū)域1的對(duì)比中發(fā)現(xiàn)，在未經(jīng)勻色處理的建筑物三維模型中(圖8左)，中間部分基本都是黑色，建筑物紋理細(xì)節(jié)無(wú)法辨認(rèn)，并且有少量空洞存在.而經(jīng)過(guò)勻色處理的三維模型中(圖8右)，未存在建筑物缺失的情況，其中地上的植被、石板以及藍(lán)色的臺(tái)子都清晰可見(jiàn)，因此區(qū)域1中經(jīng)勻色處理過(guò)的建筑物三維模型從紋理細(xì)節(jié)以及模型完整性方面都優(yōu)于未經(jīng)勻色處理的三維模型；從區(qū)域2的對(duì)比中可以看出，由于色調(diào)太暗，若不經(jīng)勻色處理(圖9左)，建筑物立面的整個(gè)墻體以及地面完全呈大片空洞狀態(tài)，極大地破壞了建筑物的完整性.經(jīng)勻色處理后(圖9右)，建筑物墻體不僅沒(méi)有空洞，其紋理細(xì)節(jié)也不存在丟失問(wèn)題，因?yàn)橹脖坏钠ヅ潆y度本身就較大且色調(diào)過(guò)暗，在地面處依然存在一部分區(qū)域無(wú)法構(gòu)建出完整的模型，但是模型的整體性以及細(xì)節(jié)的表現(xiàn)能力已經(jīng)遠(yuǎn)遠(yuǎn)優(yōu)于原三維模型.

圖8 三維模型勻色前(左)及勻色后(右)對(duì)比

圖9 三維模型勻色前(左)及勻色后(右)對(duì)比

4 結(jié) 語(yǔ)

針對(duì)無(wú)人機(jī)影像在曝光、飽和度以及色調(diào)等方面存在的問(wèn)題，對(duì)無(wú)人機(jī)影像進(jìn)行色調(diào)勻化處理，并以云南師范大學(xué)(呈貢校區(qū))為研究區(qū)，分別對(duì)未經(jīng)色調(diào)勻化處理及經(jīng)色調(diào)勻化處理后的無(wú)人機(jī)影像進(jìn)行三維建模實(shí)驗(yàn)并作對(duì)比分析，結(jié)果表明：由于曝光過(guò)強(qiáng)而使原始影像色調(diào)差異較大時(shí)，對(duì)其直接進(jìn)行三維模型構(gòu)建，其實(shí)景三維模型色調(diào)差異大，亮暗處的過(guò)渡生硬，紋理細(xì)節(jié)難以分辨并且存在較大的空洞，完整性較差；對(duì)影像進(jìn)行色調(diào)勻化處理使其色階均勻分布，勻化處理后的建筑物三維模型色調(diào)均勻，紋理細(xì)節(jié)顯著且模型完整性較高.