孟麗艷，張 旭，尹鵬飛，宋碧波

（1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院，河南 焦作 454000；2.河南恒旭力創(chuàng)測(cè)繪工程有限公司，河南 鄭州 450007；3.河南省測(cè)繪地理信息局，河南 鄭州 450007）

基于DOM與DLG模型的建筑物頂部紋理重建

孟麗艷1，張 旭2，尹鵬飛1，宋碧波3

（1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院，河南 焦作 454000；2.河南恒旭力創(chuàng)測(cè)繪工程有限公司，河南 鄭州 450007；3.河南省測(cè)繪地理信息局，河南 鄭州 450007）

提出一種基于DOM快速分割及與DLG模型自動(dòng)匹配的方式，通過(guò)利用區(qū)域生長(zhǎng)算法及形態(tài)學(xué)中的開(kāi)閉運(yùn)算，以Python腳本語(yǔ)言為開(kāi)發(fā)工具，實(shí)現(xiàn)了建筑物頂部紋理的快速重建。實(shí)驗(yàn)證明，這種方法簡(jiǎn)單、快速，適用于大范圍場(chǎng)景中建筑物三維快速重建。

DOM區(qū)域生長(zhǎng)法；開(kāi)閉運(yùn)算；Python語(yǔ)言；紋理自動(dòng)匹配

1 紋理快速構(gòu)建設(shè)計(jì)思想

針對(duì)建筑物立面紋理設(shè)計(jì)，本文對(duì)DOM進(jìn)行圖像分割并與DLG建立起來(lái)的模型進(jìn)行匹配。其設(shè)計(jì)思想是：首先選取簡(jiǎn)單直觀的區(qū)域分割算法，然后對(duì)分割后的影像進(jìn)行開(kāi)運(yùn)算和閉運(yùn)算，使分割后的影像輪廓變得光滑，再將處理好的正射影像利用坐標(biāo)匹配的方式，通過(guò)程序?qū)崿F(xiàn)紋理映射自動(dòng)化。與傳統(tǒng)的建筑物頂面紋理映射相比，該方法具有鮮明的特點(diǎn)：①表達(dá)了建筑物的真實(shí)信息，使建立的三維城市及虛擬城市更逼近真實(shí)地物；②采用批量的紋理匹配設(shè)計(jì)，適用于大范圍建筑三維模型的自動(dòng)重建；③與人工相比，它具有程序簡(jiǎn)單、速度快、現(xiàn)勢(shì)性強(qiáng)的特點(diǎn)。

2 建筑物紋理的映射

基于DOM圖像分割的建筑物頂面紋理的映射，實(shí)質(zhì)是在3DGIS平臺(tái)中對(duì)整片區(qū)域進(jìn)行影像分割和形態(tài)學(xué)處理，其核心是將處理好的正射影像自動(dòng)化或半自動(dòng)化映射到基于DLG建立的三維模型上。

2.1 正射影像的自動(dòng)分割

數(shù)字影像的表示是一個(gè)整數(shù)矩陣，其中的每個(gè)元素都表示其下標(biāo)所指位置處的電磁波強(qiáng)度，式（1）表示一個(gè)數(shù)字影像矩陣，每個(gè)值稱(chēng)為每個(gè)元素的像素，下標(biāo)m、n分別表示該像素的行列數(shù)。

常見(jiàn)的圖像分割方法有邊界檢測(cè)法、閾值分割法、區(qū)域增長(zhǎng)分割法、基于神經(jīng)網(wǎng)的分割法、聚類(lèi)分割法等[6,7]，其中，區(qū)域分割算法原理簡(jiǎn)單直觀，可以進(jìn)行串行計(jì)算，對(duì)于較均勻的聯(lián)通目標(biāo)分割效果較好，適合于輪廓比較明顯的地物提取。區(qū)域生長(zhǎng)的基本思想是將具有相似性質(zhì)的像素集合起來(lái)構(gòu)成區(qū)域。先對(duì)每個(gè)需要分割的區(qū)域找一個(gè)種子像素作為生長(zhǎng)的起點(diǎn)，然后將重點(diǎn)像素周?chē)徲蛑信c種子像素有相同或相似性質(zhì)的像素合并到種子像素所在的區(qū)域中。將這些新像素當(dāng)作新的種子像素繼續(xù)進(jìn)行上面的過(guò)程，直到再?zèng)]有滿(mǎn)足條件的像素被包括進(jìn)來(lái)。這樣一個(gè)區(qū)域就長(zhǎng)成了。選取建筑物內(nèi)一個(gè)種子點(diǎn)的灰度值，并以此為中心或中心集合與各領(lǐng)域像素灰度值作比較，將小于某一閾值T的像素合并起來(lái)并組成新的種子集合，當(dāng)灰度差值發(fā)生劇烈變化時(shí)說(shuō)明區(qū)域不能進(jìn)一步擴(kuò)張，檢測(cè)到了建筑物的邊界。一般情況下，建筑物形狀規(guī)則、輪廓明顯、相似度高，利用區(qū)域增長(zhǎng)法具有較好的分割效果。

本文選取許昌市某區(qū)域的建筑物為實(shí)驗(yàn)對(duì)象，如圖1所示。利用區(qū)域生長(zhǎng)法分割建筑，初始種子點(diǎn)的選取非常重要。這里采取人為選取種子點(diǎn)的方法。區(qū)域生長(zhǎng)的相似性準(zhǔn)則選取為圖像的灰度值，其中相似的程度T=10，建筑物的提取結(jié)果如圖2。

圖1 建筑物圖

圖2 初步提取結(jié)果

2.2 數(shù)學(xué)形態(tài)學(xué)的開(kāi)閉運(yùn)算

圖2所示的建筑物里面有小的空洞，且邊緣是不光滑的偽邊緣，這里采取數(shù)學(xué)形態(tài)的處理方法得到光滑的建筑邊界。

設(shè)G為圖像集合，A、B、C、D為結(jié)構(gòu)元素，A^為A關(guān)于原點(diǎn)的影像，數(shù)學(xué)形態(tài)學(xué)運(yùn)算是用A對(duì)G進(jìn)行操作。結(jié)構(gòu)元素的原點(diǎn)可以包含在結(jié)構(gòu)元素中，也可以不包含。定義的基本運(yùn)算有腐蝕、膨脹、開(kāi)和閉。基本的運(yùn)算為（以結(jié)構(gòu)元素A為例）：

式（2）表示用A對(duì)G膨脹得到的集合；式（3）表示用A對(duì)G腐蝕得到的集合；式（4）表示開(kāi)運(yùn)算，G先被A腐蝕，再被A膨脹；式（5）表示閉運(yùn)算，G先被A膨脹，再被A腐蝕。

從幾何的角度看，膨脹是對(duì)圖像進(jìn)行濾波處理，填充圖像空洞，對(duì)圖像進(jìn)行“擴(kuò)張”的行為。腐蝕則是對(duì)圖像進(jìn)行濾波處理、分割圖像，消除不感興趣的區(qū)域。膨脹和腐蝕是一對(duì)偶操作，而開(kāi)運(yùn)算和閉運(yùn)算也是一對(duì)關(guān)于集合及影響的偶操作。為得到合理的運(yùn)算，本文經(jīng)多次運(yùn)算和對(duì)比，選取了合適的結(jié)構(gòu)元素：

其中，A和C是方向上的一對(duì)結(jié)構(gòu)元素，而B(niǎo)和D也是X、Y方向上的一對(duì)結(jié)構(gòu)元素。對(duì)圖2所示的建筑物用A和C進(jìn)行開(kāi)運(yùn)算，用B和D進(jìn)行閉運(yùn)算，得到如圖3所示的提取結(jié)果。

圖3 最終提取結(jié)果

2.3 批量貼圖

經(jīng)過(guò)區(qū)域生長(zhǎng)法的圖像分割和形態(tài)學(xué)的開(kāi)閉運(yùn)算后，得到可以直接用于和DLG模型進(jìn)行匹配的建筑物紋理，經(jīng)過(guò)運(yùn)算得到大量的建筑物頂面紋理，而DLG模型想要搜索紋理進(jìn)行，就必須賦予每個(gè)建筑物紋理一個(gè)記錄索引的屬性，用以表明建筑物紋理與DLG模型的對(duì)應(yīng)關(guān)系（如圖4）。索引信息使用坐標(biāo)的對(duì)應(yīng)關(guān)系，生成的紋理文件使用對(duì)應(yīng)要素外框最大坐標(biāo)“x-y” 進(jìn)行命名（如表1），在3D Max 中讀取每個(gè)頂面的坐標(biāo)信息，與紋理文件命名一致的，判斷為相對(duì)應(yīng)的頂部貼圖，自動(dòng)貼在模型頂面上，從而完成一定區(qū)域內(nèi)的紋理圖形的映射。

圖4 索引信息的建立

表1 影像坐標(biāo)命名

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

本文以Windows XP為平臺(tái)，以Python腳本語(yǔ)言為開(kāi)發(fā)工具，對(duì)許昌市某區(qū)域的數(shù)據(jù)進(jìn)行處理，對(duì)DOM影像的裁切與DLG數(shù)字模型的匹配進(jìn)行實(shí)驗(yàn)。其中，實(shí)現(xiàn)了快速的影像分割，在3D Max平臺(tái)上建立了材質(zhì)庫(kù)，并使用Python腳本語(yǔ)言實(shí)現(xiàn)影像和模型的信息索引。以下是信息索引的部分代碼：

s=#()

for i=1 to selection.count do

(s[i]=selection[i])

for x=1 to s.count do (select s[x]

n1=0

n2=0

for m=1 to s[x].numverts do

(n1=n1+s[x].vertices[m].pos.x

n2=n2+s[x].vertices[m].pos.y)

通過(guò)影像信息的索引實(shí)現(xiàn)了建筑物頂部紋理的自動(dòng)映射。在圖5所示的模型中，能看到利用傳統(tǒng)建模方法得到的模型，不具備真實(shí)特性。而利用DOM影像裁切的自動(dòng)化建模方法，可以實(shí)現(xiàn)批量及自動(dòng)化，得到的模型紋理較為真實(shí)。

P237.3

B

1672-4623（2015）03-0168-02

10.3969/j.issn.1672-4623.2015.03.059城市建筑物建模主要分為幾何建模和紋理建模[1,2]，在建筑物紋理映射中主要包括建筑物的頂部紋理映射和立面紋理映射。傳統(tǒng)的頂面紋理映射大多采取手工貼圖的方式，但隨著建模范圍的擴(kuò)大，用Photoshop等繪圖軟件對(duì)DOM進(jìn)行手工裁切，然后利用3D Max等三維平臺(tái)進(jìn)行貼圖，工作量巨大，費(fèi)時(shí)費(fèi)力。很多研究人員力圖采用其他方法來(lái)快速獲取建筑物紋理，其中包括近些年發(fā)展起來(lái)的機(jī)載LiDAR技術(shù)。Alharthy等利用多次回波信息對(duì)LiDAR數(shù)據(jù)進(jìn)行提??；黃先鋒提出逐步剔除非建筑物點(diǎn)的方法[3-5]。但是這些方法對(duì)點(diǎn)云的密度要求高，數(shù)據(jù)處理量龐大。還有一些學(xué)者采用低空傾斜攝影測(cè)量的方法構(gòu)建模型，但是生成的模型不夠精細(xì)，不能體現(xiàn)細(xì)節(jié)。本文提出一種基于DOM的頂部紋理自動(dòng)匹配機(jī)制，將真實(shí)的正射影像映射到建筑物模型上，可有效解決快速建模中紋理映射問(wèn)題。

2014-04-18。

項(xiàng)目來(lái)源：國(guó)家自然科學(xué)基金資助項(xiàng)目（41001226）。