徐文學(xué)，楊必勝，魏 征，方莉娜

1.武漢大學(xué) 測繪遙感信息工程國家重點(diǎn)實(shí)驗(yàn)室，湖北武漢430079；2.武漢大學(xué)時(shí)空數(shù)據(jù)智能獲取技術(shù)與應(yīng)用教育部工程研究中心，湖北 武漢430079

1 引 言

激光掃描測量是近十幾年快速發(fā)展的一種新型測量技術(shù)，在三維城市建模、植被檢測、通信線路設(shè)計(jì)、城市規(guī)劃和災(zāi)害處理等領(lǐng)域都有廣泛的應(yīng)用前景［1］，現(xiàn)已逐漸成為三維城市數(shù)據(jù)模型獲取的一種重要方法。建筑物、道路和樹冠是城市空間中最主要、最重要的空間實(shí)體，其三維信息是重要的基礎(chǔ)地理信息。因此建筑物和樹冠的目標(biāo)提取對于城市制圖和三維建模也就具有極大的研究和實(shí)用價(jià)值［2］。

近年來，國內(nèi)外對基于LiDAR數(shù)據(jù)的建筑物目標(biāo)提取進(jìn)行了廣泛的研究。主要可以分為兩類：一類方法是只使用LiDAR數(shù)據(jù)提取建筑物目標(biāo)［3］；另一類方法是融合LiDAR數(shù)據(jù)與IKONOS影像［4］、航空影像［5］以及其他遙感影像或地面規(guī)劃設(shè)計(jì)圖［6］來提取建筑物。第1類方法，首先根據(jù)LiDAR數(shù)據(jù)得到DSM和DEM，然后利用兩者之差分離出地面點(diǎn)和非地面點(diǎn)，最后從非地面點(diǎn)中進(jìn)一步提取感興趣的建筑物目標(biāo)點(diǎn)云。這類方法的缺點(diǎn)是當(dāng)?shù)匦屋^復(fù)雜時(shí)生成的DEM精度較差，從而影響最終的提取結(jié)果。第2類方法是融合LiDAR數(shù)據(jù)與IKONOS影像［4］、航空影像［5］以及其他遙感影像或地面規(guī)劃設(shè)計(jì)圖［6］來提取建筑物。這類方法可以降低處理難度，但是對數(shù)據(jù)的同時(shí)獲取提出了更高的要求，并且地面規(guī)劃設(shè)計(jì)圖對新建房屋不能及時(shí)處理［7］。

在目標(biāo)提取問題中，引入空間信息和空間關(guān)系可以顯著提高算法的準(zhǔn)確性和可靠性。一種引入空間信息和空間關(guān)系的有效方法是使用統(tǒng)計(jì)幾何方法。該方法通過先驗(yàn)?zāi)Ｐ鸵肟臻g特性。這類方法中的一個(gè)主要模型是馬爾可夫隨機(jī)場（MRF）。MRF在很多問題中都得到了較好的效果［8］。但MRF是一種基于像素級(jí)隨機(jī)場的分析方法，對噪聲較敏感，且容易陷入局部最小值。而標(biāo)記點(diǎn)過程是一種基于高層隨機(jī)場的分析方法，可以很好的解決這些問題，該方法是一種面向?qū)ο蟮慕y(tǒng)計(jì)幾何方法。標(biāo)記點(diǎn)過程方法根據(jù)目標(biāo)的幾何特征建立吉布斯自由能變模型，通過目標(biāo)的一致性建立該模型的數(shù)據(jù)項(xiàng)，通過目標(biāo)的拓?fù)湫再|(zhì)等空間特性建立該模型的先驗(yàn)項(xiàng)，然后利用可逆跳轉(zhuǎn)馬爾可夫蒙特卡洛算法（RJMCMC）進(jìn)行采樣，并采用模擬退火算法進(jìn)行優(yōu)化求解［9］。近年來，基于標(biāo)記點(diǎn)過程的方法已被廣泛用于高分辨率影像處理中，并取得了較好的效果。文獻(xiàn)［10—11］都使用標(biāo)記點(diǎn)過程方法從DEM中提取建筑物目標(biāo)。但是，文獻(xiàn)［10］的方法參數(shù)設(shè)置復(fù)雜，針對不同的數(shù)據(jù)，需要對每一個(gè)參數(shù)進(jìn)行調(diào)試，處理起來非常麻煩。文獻(xiàn)［11］則使用目標(biāo)邊緣的梯度信息作為吉布斯自由能變模型的數(shù)據(jù)項(xiàng)，但當(dāng)目標(biāo)的邊緣不明顯或與其他目標(biāo)緊密相連時(shí)提取結(jié)果較差。文獻(xiàn)［12］在文獻(xiàn)［10］研究的基礎(chǔ)上，利用基于標(biāo)記點(diǎn)過程的方法從DEM數(shù)據(jù)中進(jìn)行建筑物三維重建。文獻(xiàn)［9］利用標(biāo)記點(diǎn)過程進(jìn)行車道線的識(shí)別和提取。文獻(xiàn)［13］采用標(biāo)記點(diǎn)過程方法用橢圓模型進(jìn)行樹冠提取，從而統(tǒng)計(jì)森林中的樹木數(shù)量。但該方法需要對遙感影像進(jìn)行預(yù)處理，把遙感影像轉(zhuǎn)化為二值影像。文獻(xiàn)［14］對直線、矩形和圓等模型進(jìn)行綜合，提出基于多標(biāo)記點(diǎn)過程的目標(biāo)提取方法，它包括7種隨機(jī)模型，能適應(yīng)多種應(yīng)用，如計(jì)數(shù)統(tǒng)計(jì)，紋理表達(dá)，路網(wǎng)提取等。但該方法中吉布斯自由能變模型的參數(shù)設(shè)置范圍較大，調(diào)試較困難。

由于激光掃描獲取的點(diǎn)云數(shù)據(jù)量大，難以直接對點(diǎn)云的幾何數(shù)據(jù)進(jìn)行分類和特征提取。故本文將機(jī)載激光掃描點(diǎn)云數(shù)據(jù)轉(zhuǎn)換成點(diǎn)云特征影像，針對激光掃描數(shù)據(jù)生成點(diǎn)云特征影像的算法已經(jīng)比較成熟［15－16］，采用文獻(xiàn)［16］中的方法生成點(diǎn)云特征影像，提出基于多標(biāo)記點(diǎn)過程的建筑物和樹冠目標(biāo)對象的自動(dòng)提取方法。該方法首先根據(jù)目標(biāo)的幾何特征建立吉布斯自由能變模型，通過目標(biāo)的一致性建立該模型的數(shù)據(jù)項(xiàng)，通過目標(biāo)的拓?fù)湫再|(zhì)等空間特性建立該模型的先驗(yàn)項(xiàng)，然后利用RJMCMC算法采樣，并采用模擬退火算法進(jìn)行優(yōu)化求解，實(shí)現(xiàn)建筑物目標(biāo)和樹冠目標(biāo)幾何對象的多目標(biāo)自動(dòng)提取，并同時(shí)得到建筑物和樹木的數(shù)量、中心點(diǎn)位置、目標(biāo)的大小以及目標(biāo)的方向角等信息。本文算法流程圖見圖1。

圖1 目標(biāo)提取流程圖Fig.1 The flow chart of object extraction based on marked point process

2 標(biāo)記點(diǎn)過程

在廣泛應(yīng)用于圖像分析的隨機(jī)方法中，標(biāo)記點(diǎn)過程具有組合全局信息來識(shí)別幾何形狀的優(yōu)點(diǎn)，該方法是一種面向?qū)ο蟮慕y(tǒng)計(jì)幾何方法。一個(gè)按一定統(tǒng)計(jì)規(guī)律在某空間中隨機(jī)分布的點(diǎn)集就形成一個(gè)隨機(jī)點(diǎn)過程（簡稱點(diǎn)過程）。隨機(jī)點(diǎn)過程是描述隨機(jī)點(diǎn)分布的隨機(jī)過程［17］。標(biāo)記點(diǎn)過程就是帶有輔助變量的點(diǎn)過程，這個(gè)輔助變量即為聯(lián)系于每一點(diǎn)的標(biāo)記。

對點(diǎn)云特征影像上的建筑物、樹木等不同目標(biāo)的提取過程可以看做隨機(jī)過程（空間點(diǎn)過程）的一種 —— 泊松過程（Poisson process）。需要處理的數(shù)據(jù)構(gòu)成的有界區(qū)域組成一個(gè)框架空間Λ，一個(gè)定義在框架空間Λ上具有密度估計(jì)μ和密度函數(shù)ρ的泊松過程X滿足［14］：對任意的有界區(qū)域B?R2，μ（B）＞0，滿足N（B）為均值為μ（B）的泊松分布；已知N（B），XB中的點(diǎn)相互獨(dú)立，且密度與ρ（u），u∈B成比例。因此可利用具有強(qiáng)度ν（·）的齊次泊松分布和框架空間Λ上的非負(fù)函數(shù)h（·）定義概率密度測度μ（·）：μ（B）＝

本文結(jié)合建筑物和樹冠的形狀，分別使用矩形和橢圓作為建筑物和樹冠的標(biāo)記（圖2）。

圖2 矩形和橢圓的位置和標(biāo)記信息Fig.2 Position and marks of rectangle and ellipse

矩形和橢圓的位置和標(biāo)記信息如下

式中，x、y表示矩形和橢圓的中心點(diǎn)坐標(biāo)；l、w分別表示矩形長和寬的二分之一；a、b分別表示橢圓的長半軸和短半軸；θ表示矩形或橢圓的長軸與X軸之間的夾角。根據(jù)上述信息可以得到建筑物和樹冠的中心點(diǎn)位置、目標(biāo)的大小以及目標(biāo)的方向角等信息。圖3展示了一個(gè)點(diǎn)過程及其對應(yīng)的標(biāo)記點(diǎn)過程。圖3（a）中的“＋”表示框架空間中隨機(jī)分布的一系列點(diǎn)形成的點(diǎn)過程，圖3（b）則表示對點(diǎn)過程中的每個(gè)點(diǎn)分別添加矩形或橢圓標(biāo)記，從而形成標(biāo)記點(diǎn)過程。

圖3 點(diǎn)過程和標(biāo)記點(diǎn)過程Fig.3 A point process and a marked point process

3 吉布斯自由能變模型

本文考慮提取目標(biāo)的標(biāo)記點(diǎn)過程符合Poisson分布，其概率密度為可以通過兩種方式定義：一是在貝葉斯框架下，但是這種情況需要知道數(shù)據(jù)的DEM；二是通過吉布斯自由能變模型［11］。由于在復(fù)雜地形下獲取數(shù)據(jù)的DEM比較困難，且精度較低，因此本文采用第2種方式，通過吉布斯自由能變模型U（X）說明點(diǎn)過程，吉布斯函數(shù)的測度定義為

對于吉布斯自由能變模型U（X），本文主要考慮數(shù)據(jù)一致性因子Ud（X）和先驗(yàn)約束因子Up（X）。

本文主要目的是提取建筑物輪廓，同時(shí)可以提取樹冠輪廓。下面的說明皆以矩形目標(biāo)為例，同理可得到對橢圓目標(biāo)的操作。為均衡數(shù)據(jù)一致性因子和先驗(yàn)約束因子對吉布斯自由能變值的貢獻(xiàn)，本文把各類因子的取值范圍歸一化到［－1，1］之內(nèi)。

3.1 數(shù)據(jù)一致性因子

Ud（X）累積目標(biāo)框架X中每個(gè)目標(biāo)xi的局部吉布斯自由能變值

式中，Ud（xi）是目標(biāo)xi的吉布斯自由能變模型的數(shù)據(jù)一致性因子。該因子滿足目標(biāo)之間相互獨(dú)立且目標(biāo)的吉布斯自由能變值具有負(fù)值兩個(gè)重要條件［14］。本文認(rèn)為目標(biāo)內(nèi)部區(qū)域具有一致性，與外部區(qū)域具有差異性。Ud（xi）定義為

式中，dm表示目標(biāo)內(nèi)部區(qū)域和外部區(qū)域之間的馬氏距離，（min，σin）、（mout，σout）分別表示內(nèi)部區(qū)域和外部區(qū)域的平均像素值和標(biāo)準(zhǔn)偏差；n1和n2分別表示內(nèi)部區(qū)域和外部區(qū)域的面積（見圖4）。外部區(qū)域的寬度設(shè)置為固定值，一般根據(jù)影像的分辨率設(shè)置為2～4。馬氏距離是一種有效的計(jì)算兩個(gè)未知樣本集的相似度的方法，它考慮到各種特性之間的聯(lián)系，并且是尺度無關(guān)的。文獻(xiàn)［18］使用巴氏距離建立吉布斯自由能變模型的數(shù)據(jù)一致性因子，但是當(dāng)巴氏距離的第1項(xiàng)和第2項(xiàng)為同一數(shù)量級(jí)時(shí)，此計(jì)算方法會(huì)產(chǎn)生較大誤差。閾值d0是模型唯一需要調(diào)試的參數(shù)，調(diào)節(jié)數(shù)據(jù)擬合的敏感性，控制吉布斯自由能變模型的數(shù)據(jù)一致性因子的取值。當(dāng)矩形目標(biāo)和特征影像中的目標(biāo)匹配的較好時(shí)，吉布斯自由能變模型的數(shù)據(jù)一致性因子取值應(yīng)為負(fù)值。

圖4 目標(biāo)和數(shù)據(jù)的一致性測度Fig.4 Coherence measure between object and data

3.2 先驗(yàn)約束因子

Up（X）引入目標(biāo)之間的交互作用和目標(biāo)梯度等先驗(yàn)知識(shí)。本文限制目標(biāo)間的重疊，建立適合非重疊目標(biāo)的一般模型，并加入強(qiáng)結(jié)構(gòu)信息。

每一區(qū)域只能有一種標(biāo)記，對具有重疊區(qū)域的目標(biāo)進(jìn)行懲罰。定義兩個(gè)重疊目標(biāo)的懲罰因子的計(jì)算公式為（見圖5）

圖5 重疊目標(biāo)的量化Fig.5 Two overlapping objects and the quality

該先驗(yàn)知識(shí)能夠體現(xiàn)目標(biāo)之間的重疊關(guān)系，懲罰有重疊的目標(biāo)，通過該因子可以對簡單模型進(jìn)行組合，從而提取復(fù)雜目標(biāo)。

本文定義一個(gè)強(qiáng)制約束，對中心點(diǎn)距離非常近的目標(biāo)賦予非常大的吉布斯自由能變值，以杜絕同一區(qū)域具有不同的標(biāo)記。距離閾值h根據(jù)具體的數(shù)據(jù)確定（式6）。

考慮到目標(biāo)的邊緣梯度對目標(biāo)的確定具有較大作用，本文加入梯度約束。對矩形每條邊上的點(diǎn)搜索其八鄰域，計(jì)算每個(gè)點(diǎn)與其八鄰域內(nèi)的點(diǎn)的灰度差的絕對值的最大值，把該值作為該點(diǎn)的灰度梯度值，當(dāng)一條邊的平均灰度梯度值小于一定閾值時(shí)對該目標(biāo)進(jìn)行懲罰。計(jì)算公式為

式中，nxi表示每個(gè)目標(biāo)中小于給定灰度梯度閾值的邊的數(shù)目。該梯度閾值可以通過生成目標(biāo)的梯度圖像提前確定，不需要進(jìn)行反復(fù)調(diào)試。

綜上所述，計(jì)算吉布斯自由能變模型的先驗(yàn)因子的公式為

4 優(yōu)化算法

本算法的實(shí)質(zhì)是求解吉布斯自由能變模型U（X）的全局最小能量，實(shí)現(xiàn)實(shí)體對象與點(diǎn)云數(shù)據(jù)的最佳匹配，從而識(shí)別出點(diǎn)云數(shù)據(jù)中蘊(yùn)含的幾何對象。該方法克服了傳統(tǒng)方法只能得到局部最優(yōu)的不足。本文采用RJMCMC算法進(jìn)行采樣，使用模擬退火算法進(jìn)行優(yōu)化［9］。下面對RJMCMC算法和模擬退火算法兩個(gè)關(guān)鍵技術(shù)進(jìn)行說明。

4.1 RJMCMC算法

馬爾可夫蒙特卡洛（MCMC）算法是一種簡單且行之有效的貝葉斯計(jì)算方法，該迭代算法并不依賴于初始狀態(tài)空間［9］。文獻(xiàn)［19］提出了RJMCMC方法，該算法是一種特殊的Metropolis－Hastings算法，可以在不同維度之間可逆跳轉(zhuǎn)，被廣泛地應(yīng)用于最大化問題。決定跳轉(zhuǎn)的概率為計(jì)算得到的Green比。

RJMCMC算法對目標(biāo)的參數(shù)進(jìn)行修改，并且改變狀態(tài)空間的大小。該算法令q（Y｜X）為任意一個(gè)轉(zhuǎn)移核，假設(shè)已經(jīng)生成狀態(tài)空間X0、X1、…、Xi，為了生成狀態(tài)空間Xi＋1，進(jìn)行如下步驟：

（1）產(chǎn)生一個(gè)轉(zhuǎn)移核，并根據(jù)轉(zhuǎn)移核得到新的狀態(tài)空間Y～q（Y｜Xi）。

（2）計(jì)算 Green比r＝r（Xi，Y），其中

式中，f（Xi）和f（Y）分別表示狀態(tài)空間Xi和狀態(tài)空間Y中所有目標(biāo)的吉布斯自由能變值；q（Xi｜Y）和q（Y｜Xi）分別表示從狀態(tài)空間Xi轉(zhuǎn)移到狀態(tài)空間Y的概率和從狀態(tài)空間Y轉(zhuǎn)移到狀態(tài)空間Xi的概率。

對于接受率r（Xi，Y）的計(jì)算需保證該跳轉(zhuǎn)是可逆跳轉(zhuǎn)。實(shí)際上，轉(zhuǎn)移核由不同建議的轉(zhuǎn)移核組成，不同的轉(zhuǎn)移核對應(yīng)不同的Green比r（Xi，Y）。本文中各種轉(zhuǎn)移核的Green比的計(jì)算見文獻(xiàn)［13，14，19］。

為保證馬爾可夫鏈可逆并收斂到最優(yōu)解，本文采用7種轉(zhuǎn)移核，分別為：生、滅、平移、旋轉(zhuǎn)、縮放、分裂和合并。均勻生滅轉(zhuǎn)移核保證馬爾可夫鏈?zhǔn)諗康阶顑?yōu)解，非跳躍變換轉(zhuǎn)移核和分裂合并轉(zhuǎn)移核加速馬爾可夫鏈的收斂［9］。為加速馬爾可夫鏈的收斂，本文對不同轉(zhuǎn)移核的選擇賦予不同的概率。實(shí)際中，一般對平移、旋轉(zhuǎn)和縮放3種轉(zhuǎn)移核賦予較大的概率［13］。

均勻生滅轉(zhuǎn)移核是指在當(dāng)前狀態(tài)空間中隨機(jī)添加或移除一個(gè)目標(biāo)，這一變換對應(yīng)于高維狀態(tài)空間和低維狀態(tài)空間之間的跳轉(zhuǎn)。向高維跳轉(zhuǎn)相當(dāng)于生的過程，而向低維跳轉(zhuǎn)則相當(dāng)于滅的過程。圖6（a）中黑色目標(biāo)表示在原來的狀態(tài)空間中隨機(jī)加入一個(gè)目標(biāo)，該目標(biāo)的位置和形狀參數(shù)都是隨機(jī)生成的，灰色目標(biāo)則表示在原來的狀態(tài)空間中隨機(jī)的移除一個(gè)目標(biāo)。

因?yàn)槠揭啤⑿D(zhuǎn)和縮放3種轉(zhuǎn)移核不改變目標(biāo)的數(shù)量，只是對當(dāng)前目標(biāo)的參數(shù)進(jìn)行修改，故稱之為非跳躍變換轉(zhuǎn)移核。這3種轉(zhuǎn)移核是為了使矩形更好的擬合目標(biāo)。圖6（b）中黑色目標(biāo)表示對選取目標(biāo)進(jìn)行平移、旋轉(zhuǎn)和縮放的結(jié)果，每一次縮放只對目標(biāo)的一條邊進(jìn)行修改。

分裂合并轉(zhuǎn)移核避免對提取的目標(biāo)探測不足或過探測。在有些情況下，狀態(tài)空間陷于局部吉布斯自由能變最小值，均勻生滅轉(zhuǎn)移核和非跳躍變換轉(zhuǎn)移核不能解決該問題，從而提出分裂合并轉(zhuǎn)移核（見圖6（c）），該轉(zhuǎn)移核可以有效地解決這一問題。

圖6 RJMCMC轉(zhuǎn)移核Fig.6 Transition kernels of RIMCMC

4.2 模擬退火算法

模擬退火算法是一種通用概率演算法，用來在一個(gè)大的搜尋空間內(nèi)找尋命題的最優(yōu)解［20］。該算法可以確保在任何初始狀態(tài)空間下收斂于全局最優(yōu)解。模擬退火算法具有漸近收斂性，已在理論上被證明是收斂于全局最優(yōu)解的全局優(yōu)化算法。模擬退火算法的應(yīng)用很廣泛，但其參數(shù)難以控制。實(shí)際應(yīng)用中，由于必須考慮計(jì)算復(fù)雜度的切實(shí)可行性等問題，常采用如下的降溫方式

式中，k為正的小于1的常數(shù)；t為降溫次數(shù)。

5 試驗(yàn)與分析

本文采用ISPRS的機(jī)載激光掃描測試數(shù)據(jù)進(jìn)行試驗(yàn)，其中包含樹木、建筑物、道路等目標(biāo)。分別對兩組不同的數(shù)據(jù)進(jìn)行試驗(yàn)。第1組數(shù)據(jù)的區(qū)域掃描范圍約為247m×179m，共146 235個(gè)激光點(diǎn)。點(diǎn)云及生成的0.25m分辨率的點(diǎn)云特征影像見圖7（a）、（b）。點(diǎn)云特征影像的大小為：988像素×715像素。第2組數(shù)據(jù)的區(qū)域掃描范圍約為275m×320m，共351 692個(gè)激光點(diǎn)。點(diǎn)云及生成的0.25m分辨率的點(diǎn)云特征影像見圖8（a）、（b）。點(diǎn)云特征影像的大小為：1098像素×1280像素。第1組數(shù)據(jù)的建筑物結(jié)構(gòu)相對比較簡單，第2組數(shù)據(jù)的建筑物結(jié)構(gòu)比較復(fù)雜。

圖7 點(diǎn)云特征影像及多目標(biāo)提取結(jié)果Fig.7 Point cloud feature image and multi－objects extraction result

圖8 點(diǎn)云特征影像及多目標(biāo)提取結(jié)果Fig.8 Point cloud feature image and multi－objects extraction result

由于兩組數(shù)據(jù)屬于同一區(qū)域，本文對兩組數(shù)據(jù)的試驗(yàn)參數(shù)進(jìn)行相同的設(shè)置。根據(jù)建筑物及樹冠的大小，設(shè)定矩形的長和寬的取值范圍分別為［30，200］和［30，100］，橢圓的長短半軸的取值范圍為［16，32］，避免目標(biāo)中心點(diǎn)重疊的距離閾值設(shè)置為橢圓短半軸的最小值；對于模擬退火的參數(shù)，一般初始溫度為1000，最低溫度為0.000 1，降溫因子0.999 5，每一溫度下運(yùn)行次數(shù)100；泊松分布強(qiáng)度為提取目標(biāo)的大體數(shù)量，與最終提取目標(biāo)的數(shù)量無關(guān)，本文設(shè)置為70；通過對點(diǎn)云特征影像的梯度影像測試，梯度閾值為15時(shí)，目標(biāo)的邊緣較明確，且較完整。圖7（c）和圖8（c）分別展示了點(diǎn)云特征影像中建筑物和樹冠的提取結(jié)果，其中橢圓形目標(biāo)為樹冠，矩形目標(biāo)為建筑物（不考慮具有橢圓形等特殊形狀的建筑物和特殊形狀的樹木）。圖7（d）、圖8（d）和圖7（e）、圖8（e）則分別為提取建筑物和樹冠的點(diǎn)云以及提取建筑物和樹冠后的點(diǎn)云。試驗(yàn)結(jié)果見表1。

表1 試驗(yàn)結(jié)果Tab.1 Experimental result

由于第2組數(shù)據(jù)中建筑物的結(jié)構(gòu)比較復(fù)雜，因此第2組數(shù)據(jù)中目標(biāo)的馬氏距離要比第1組數(shù)據(jù)中目標(biāo)的馬氏距離小一些，因此數(shù)據(jù)一致性因子閾值d0的取值也就小一些。通過圖7（c）、（d）、（e）和圖8（c）、（d）、（e）可以看出，該方法對于比較規(guī)則的、可分離的建筑物和樹冠提取效果很好，且邊界精確；對于比較復(fù)雜的建筑物和樹冠也能粗略的提取出其輪廓，可根據(jù)該結(jié)果進(jìn)行進(jìn)一步的精細(xì)處理，確定目標(biāo)的具體邊緣。在第1組數(shù)據(jù)中可以提取出所有建筑物，并且提取的目標(biāo)的邊緣比較精確。對于第2組數(shù)據(jù)，結(jié)構(gòu)比較簡單的建筑物能得到較精確的提取結(jié)果，對比較復(fù)雜的建筑物也能粗略的提取出其輪廓。該方法也可以提取出大部分樹冠，但是少部分樹冠無法提出。總體來說，該試驗(yàn)結(jié)果達(dá)到了預(yù)期的目的。部分樹冠無法提取出的原因如下：目標(biāo)處于點(diǎn)云特征影像的邊緣（無法計(jì)算數(shù)據(jù)一致性自由能變值）、目標(biāo)較小或較大（超出提取目標(biāo)給定的大小范圍）、目標(biāo)形狀不規(guī)則、樹冠與建筑物或樹冠之間距離太近（不具有可分性）等因素。當(dāng)樹冠的形狀接近于矩形時(shí)，可能被錯(cuò)誤地提取為矩形目標(biāo)，這主要是因?yàn)楸疚闹械募妓棺杂赡茏兡Ｐ椭皇强紤]目標(biāo)內(nèi)部具有一致性，而沒有加入目標(biāo)的特征分析，僅依靠吉布斯自由能變模型不能對這類目標(biāo)進(jìn)行區(qū)分。對于建筑物和樹冠的詳細(xì)區(qū)分可以根據(jù)目標(biāo)點(diǎn)云的空間密度分析等方法進(jìn)行進(jìn)一步區(qū)分。本文采用正確率和完整率對提取目標(biāo)進(jìn)行評價(jià)，對兩組數(shù)據(jù)進(jìn)行目標(biāo)提取的正確率與完整率見表2。

通過試驗(yàn)分析可知，數(shù)據(jù)一致性因子閾值d0用于歸一化目標(biāo)的數(shù)據(jù)一致性因子的取值。馬氏距離的計(jì)算為目標(biāo)內(nèi)部區(qū)域和外部區(qū)域的均值差的平方與協(xié)方差的比，是一種相對關(guān)系，與點(diǎn)云特征影像的分辨率和提取目標(biāo)的大小無關(guān)；與目標(biāo)邊緣的規(guī)則度以及與目標(biāo)相鄰的區(qū)域是否有其他目標(biāo)存在相關(guān)。針對本文處理的數(shù)據(jù)，d0的取值范圍為［1，1.5］。

表2 提取目標(biāo)的正確率與完整率Tab.2 The correctness and completeness of the extraction object

本文算法只有一個(gè)參數(shù)需要調(diào)試，并把目標(biāo)的梯度信息作為先驗(yàn)信息，可以克服文獻(xiàn)［10—11］方法的不足，能夠有效地提取出建筑物和樹冠。同時(shí)把吉布斯自由能變模型的數(shù)據(jù)一致性因子歸一化到［－1，1］內(nèi)，解決了文獻(xiàn)［12］方法中吉布斯自由能變模型參數(shù)設(shè)置閾值范圍過大、調(diào)試比較困難的問題。另外，本文算法不需要對點(diǎn)云特征影像進(jìn)行預(yù)處理，消除了對影像進(jìn)行預(yù)處理時(shí)造成目標(biāo)及其邊緣模糊的影響，提取結(jié)果更加精確。

與所有基于采樣的方法一樣，該方法的缺點(diǎn)是運(yùn)行時(shí)間相對較長。該方法的效率和很多因素相關(guān)：數(shù)據(jù)大小、提取目標(biāo)大小、目標(biāo)間差異、目標(biāo)數(shù)量、提取結(jié)果精度、背景復(fù)雜度、模擬退火中參數(shù)設(shè)置等。提高算法的效率，引入布朗運(yùn)動(dòng)等一些高效的優(yōu)化算法是下一階段的重點(diǎn)研究內(nèi)容。

利用本文中的算法，可以同時(shí)得到建筑物和樹木的數(shù)量、中心點(diǎn)位置、目標(biāo)的大小以及目標(biāo)的方向角等信息，可以用于估計(jì)城區(qū)的建筑物和植被綠化面積、確定和標(biāo)注建筑物和樹冠的位置和大小，對于城建部門掌握管理城市的建筑物及綠化狀況、地圖測繪人員豐富城市地圖信息和三維建模等方面有著非常重要的使用價(jià)值。

6 結(jié) 論

本文利用多標(biāo)記點(diǎn)過程方法，從根據(jù)機(jī)載激光掃描數(shù)據(jù)生成的點(diǎn)云特征影像中自動(dòng)提取建筑物和樹冠。本文算法的優(yōu)點(diǎn)為：①本方法是一種基于對象的統(tǒng)計(jì)幾何方法，可以較好地消除噪聲對目標(biāo)提取結(jié)果的影響；②在點(diǎn)云特征影像中無需進(jìn)行預(yù)處理，直接提取目標(biāo)；③構(gòu)建新的、能準(zhǔn)確提取目標(biāo)的吉布斯自由能變模型。試驗(yàn)證明該方法可以效地提取建筑物和樹冠及其相應(yīng)信息，并具有較強(qiáng)的穩(wěn)健性。缺點(diǎn)為比較費(fèi)時(shí)。下一步可以利用蛇形算法等方法跟蹤出不規(guī)則建筑物和樹冠的具體形狀，使提取結(jié)果更加精確。同時(shí)，提高算法的效率也是本方法需要解決的問題之一。

［1］ XIA Bing，SU Guozhong，ZHOU Mei，et al.Algorithms

［2］ LI Yunfan，MA Hongchao.New Method for Building Planar Objects Extraction from LiDAR Data［J］.Computer Engineering and Applications，2011，47（10）：5－7.（李云帆，馬洪超.從LiDAR數(shù)據(jù)中提取建筑物平面目標(biāo)的新方法［J］.計(jì)算機(jī)工程與應(yīng)用，2011，47（10）：5－7.）

［3］ KIM K，SHAN J.Building Roof Modeling from Airborne Laser Scanning Data Based on Level Set Approach［J］.7ISPRS Journal of Photogrammetry and Remote Sensing，2011，66（4）：1－14.

［4］ SOHN G，DOWMAN I.Data Fusion of Highresolution Satellite Imagery and LiDAR Data for Automatic Building Extraction［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2007（62）：43－63.

［5］ CHENG Liang，GONG Jianya，LI Manchun，et al.3D Building Model Reconstruction from Multi－view Aerial Images and LiDAR Data［J］.Acta Geodaetica et Cartographica Sinica，2009，38（6）：494－501.（程亮，龔健雅，李滿春，等.集成多視航空影像與LiDAR數(shù)據(jù)重建3維建筑物模型［J］.測繪學(xué)報(bào)，2009，38（6）：494－501.）

［6］ VOSSELMAN G，DIJKMAN S.3DBuilding Model Reconstruction from Point Clouds and Ground Plan［J］.International Archives and Remote Sensin， 2001，4（3W4）：37－43.

［7］ ZENG Qihong，MAO Jianhua，LI Xianhua，et al.Building Reconstruction from Airborne LiDAR Points Cloud Data［J］.Geomatics and Information Science of Wuhan University，2011，36（3）：321－324.（曾齊紅，毛建華，李先華，等.機(jī)載LiDAR點(diǎn)云數(shù)據(jù)的建筑物重建研究［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2011，36（3）：321－324.）

［8］ LIU Linlin.The Application of Marked Point Process in Road Extraction［D］.Beijing：Institute of Electronics，Chinese Academy of Sciences，2007.（劉琳琳.標(biāo)記點(diǎn)過程在道路提取中的應(yīng)用［D］.北京：中國科學(xué)院電子學(xué)研究所，2007.）

［9］ TOURNAIRE O，PAPARODITIS N.A Geometric Stochastic Approach Based on Marked Point Processes for Road Mark Detection from High Resolution Aerial Images［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2009，64（6）：621－631.

［10］ ORTNER M，DESCOMBES X，ZERUBIA J.Building Outline Extraction from Digital Elevation Models Using Marked Point Processes［J］.International Journal of Computer Vision，2007，72（2）：107－132.

［11］ TOURNAIRE O，BREDIF M，BOLDO D，et al.An Efficient Stochastic Approach for Building Footprint Extraction from Digital Elevation Models［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2010，65（4）：317－327.

［12］ LAFARGE F，DESCOMBES X.ZERUBIA J，et al.Automatic Building Extraction from Dems Using an Object Approach and Application to the 3D－city Modeling［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2008，63（3）：365－381.

［13］ PERRIN G，DESCOMBES X，ZERUBIA J.Point Processes in Forestry：An Application to Tree Crown Detection［R］.Sophia Antipolis：INRIA，2005.

［14］ LAFARGE F，GIMELFARB G，DESCOMBES X.Geometric Feature Extraction by a Multi－marked Point Process［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2010，32（9）：1597－1609.

［15］ LI Bijun ，F(xiàn)ANG Zhixiang，REN Juan.Extraction of Building’s Feature from Laser Scanning Data［J］.Geomatics and Information Science of Wuhan University，2003，28（1）：65－70.（李必軍，方志祥，任娟.從激光掃描數(shù)據(jù)中進(jìn)行建筑物特征提取研究［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2003，28（1）：65－70.）

［16］ YANG Bisheng，WEI Zheng，LI Qingquan，et al.A Classification－oriented Method of Feature Image Generation for Vehicle－borne Laser Scanning Point Clouds［J］.Acta Geodaetica et Cartographica Sinica，2010，39（5）：540－545.（楊必勝，魏征，李清泉，等.面向車載激光掃描點(diǎn)云快速分類的點(diǎn)云特征圖像生成方法［J］.測繪學(xué)報(bào)，2010，39（5）：540－545.）

［17］ SNYDER D L.Random Point Processes［M］.New York：John Wiley ＆Sons，1975.

［18］ XUAN Guorong，CHAI Peiqi.Feature Selection Bhattacharyya Distance［J］.Pattern Recognition and Artificial I，1996，9（4）：324－329.（宣國榮，柴佩琪.基于巴氏距離的特征選擇［J］.模式識(shí)別與人工智能，1996，9（4）：324－329.）

［19］ GREEN P.Reversible Jump Markov Chains Monte Carlo Computation and Bayesian Model Determination［J］.Biometrika，1995，82（4）：711－732.

［20］ LAM J.An Efficient Simulated Annealing Schedule［D］.New Haven：Yale University，1988.