巫兆聰，胡忠文，張 謙，崔衛(wèi)紅

武漢大學(xué) 遙感信息工程學(xué)院，湖北 武漢430079

1 引 言

影像分割是面向?qū)ο笥跋穹治龅幕A(chǔ)，其結(jié)果將對(duì)后續(xù)的分析處理產(chǎn)生重要的影響。遙感圖像中的光譜信息是影像解譯最基本的特征，也是諸多分割算法的基礎(chǔ)，如區(qū)域增長(zhǎng)、分水嶺算法［1］等。影像紋理是地物的空間分布模式，也是進(jìn)行影像解譯的重要特征，被廣泛應(yīng)用于影像分割中［2－4］。地物在影像中表現(xiàn)為一定光譜、紋理與形狀特征的綜合體，采用單一特征進(jìn)行分割具有較大的局限性，多特征結(jié)合的影像分割方法已經(jīng)成為當(dāng)前研究的熱點(diǎn)。許多學(xué)者提出了結(jié)合光譜與紋理的分割方法［5］，如JSEG［6－7］、高斯混合模型［8－9］、馬爾可夫隨機(jī)場(chǎng)［10－11］等。這些算法針對(duì)內(nèi)容相對(duì)簡(jiǎn)單的自然影像和合成影像取得了一定的成果，但是對(duì)于數(shù)據(jù)量巨大且地物類型多樣的遙感影像仍難以取得較為理想的結(jié)果。在遙感影像分割中，多特征綜合應(yīng)用也一直是研究的熱點(diǎn)和難點(diǎn)。最具代表的FNEA［12］算法結(jié)合了光譜與幾何信息實(shí)現(xiàn)了多尺度的影像分割，但其未能有效結(jié)合影像的紋理信息，對(duì)紋理豐富地區(qū)仍難以達(dá)到理想效果。

本文針對(duì)目前遙感影像分割中多特征綜合應(yīng)用問(wèn)題，提出一種結(jié)合光譜、紋理與形狀結(jié)構(gòu)信息的分割方法。新算法以初始過(guò)分割為起點(diǎn)，采用基于鄰接圖的方法對(duì)過(guò)分割區(qū)域進(jìn)行合并獲取分割結(jié)果。在度量相鄰區(qū)域相似性時(shí)，分別以光譜直方圖距離和LBP（local binary patterns）直方圖距離度量光譜和紋理差異，并實(shí)現(xiàn)了光譜與紋理特征的自動(dòng)加權(quán)。在此基礎(chǔ)上，本文算法還利用相鄰區(qū)域的公共邊界進(jìn)行形狀約束。試驗(yàn)結(jié)果表明，本文算法能充分利用影像中的多種信息，實(shí)現(xiàn)準(zhǔn)確而高效的分割。

2 方 法

本文算法分割主要流程：首先對(duì)輸入影像進(jìn)行初始的分割，獲取最低層次的同質(zhì)區(qū)域，同時(shí)，提取影像的LBP紋理信息；然后再統(tǒng)計(jì)各區(qū)域的光譜與紋理信息，計(jì)算相鄰區(qū)域的光譜、紋理和形狀異質(zhì)性，并將所有區(qū)域都映射為區(qū)域鄰接圖中的結(jié)點(diǎn)。當(dāng)兩區(qū)域相鄰時(shí)，圖中對(duì)應(yīng)的節(jié)點(diǎn)用邊緣連接，邊緣權(quán)重即對(duì)應(yīng)兩節(jié)點(diǎn)合并的代價(jià)。最后將區(qū)域合并視為圖像逼近問(wèn)題，在鄰接圖的基礎(chǔ)上采用逐步迭代優(yōu)化的方式逐步合并邊緣權(quán)重最小的節(jié)點(diǎn)，獲得最終結(jié)果。整體流程如圖1。

圖1 分割流程圖Fig.1 Flowchart of the proposed approach

2.1 LBP紋理特征

局部二進(jìn)制模式（local binary patterns，LBP）是由 Ojala［3］于1996年提出并逐步完善的一種紋理描述算子，該算子計(jì)算局部窗口內(nèi)像素的局部空間結(jié)構(gòu)與灰度反差，并采用統(tǒng)計(jì)分析的方法來(lái)描述紋理。該算子計(jì)算簡(jiǎn)單，并且紋理描述能力強(qiáng)，被廣泛應(yīng)用于紋理圖像分析中，如紋理圖像的分割［3，15］、紋理圖像分類［13－14］等。

LBP算子描述像素的局部結(jié)構(gòu)信息，以中心像素灰度為閾值，對(duì)鄰域像素二值分割為亮像元與暗像元，并與等大小的權(quán)重模板卷積得到結(jié)構(gòu)信息。其定義如下

式中

P表示半徑R時(shí)鄰域點(diǎn)的個(gè)數(shù)；gc為中心點(diǎn)的灰度值；gp為鄰域點(diǎn)灰度值。LBP描述的為紋理的結(jié)構(gòu)信息，對(duì)于紋理的強(qiáng)度以局部反差（local contrast，LC）進(jìn)行描述，其定義為

式中，sum1表是窗口內(nèi)亮像元灰度和，即sp為1的所有像元灰度值的和；n1為對(duì)應(yīng)像元的個(gè)數(shù)；sum2和n2分別表示暗像元的灰度和與個(gè)數(shù)。LBP和LC相結(jié)合將能同時(shí)描述紋理的結(jié)構(gòu)和強(qiáng)度信息，對(duì)于不同紋理具有良好的區(qū)分性。圖2所示為3×3窗口內(nèi)計(jì)算LBP與LC的示意圖。

圖2 LBP和LC的計(jì)算Fig.2 Calculation of LBP and LC values

原始LBP算法鄰域像素權(quán)重值和位置固定，這使LBP值對(duì)紋理的方向較為敏感。為此，文獻(xiàn)［13］中進(jìn)一步提出了旋轉(zhuǎn)不變的LBP，即不斷旋轉(zhuǎn)圓形鄰域得到一系列初始定義的LBP值，取其最小值作為該鄰域的LBP值。

式中，ROR（x，i）表示對(duì)P位二進(jìn)制數(shù)x進(jìn)行向右循環(huán)旋轉(zhuǎn)位移i次。當(dāng)鄰域?yàn)?×3時(shí)，P＝8，旋轉(zhuǎn)不變LBP模式僅有36種可能值。本文中根據(jù)LBP值出現(xiàn)的序號(hào)，將相應(yīng)的LBP值換算成序號(hào)值，方便后續(xù)計(jì)算［17］。

LBP紋理反應(yīng)局部的結(jié)構(gòu)信息，對(duì)于多光譜影像，紋理在不同波段上具有高度的一致性，為了減少算法計(jì)算量，將各波段影像求均值得到灰度影像，并利用此灰度影像計(jì)算影像的LBP／C紋理值。

2.2 紋理與光譜特征異質(zhì)性

區(qū)域LBP／C紋理特征表現(xiàn)為統(tǒng)計(jì)信息，即區(qū)域紋理直方圖之間的相似性。計(jì)算直方圖距離的方法有直方圖交叉法、交叉熵方法、G統(tǒng)計(jì)法等。其中G統(tǒng)計(jì)法是一種非參數(shù)的統(tǒng)計(jì)方法，不需要預(yù)期知道隨機(jī)變量的分布情況，因此被廣泛應(yīng)用于基于LBP紋理的影像分析中［3，14，17］。對(duì)于兩個(gè)隨機(jī)變量集合m和n，令f表示概率密度函數(shù)，G統(tǒng)計(jì)公式為

對(duì)于任意同質(zhì)區(qū)域，其直方圖累積概率值均為1，因此式（5）中第2項(xiàng)和第4項(xiàng)的值始終為常數(shù)，即

將式（6）和式（7）代入式（5）可得

本文中紋理的描述綜合了紋理結(jié)構(gòu)與紋理強(qiáng)度，即采用LBP與LC的聯(lián)合概率分布直方圖進(jìn)行度量。為了減少計(jì)算復(fù)雜度，同時(shí)考慮紋理的可分性，將LC量化為8個(gè)等級(jí)，并與LBP構(gòu)成36×8的聯(lián)合概率直方圖［3］，其結(jié)合具有極強(qiáng)的紋理描述和區(qū)分能力。則相鄰區(qū)域的紋理異質(zhì)性定義為L(zhǎng)BP／C聯(lián)合直方圖之間的G統(tǒng)計(jì)值

對(duì)象的光譜特征量化等級(jí)越高，對(duì)地物的區(qū)分能力越強(qiáng)，但是較高的量化等級(jí)使面積有限的區(qū)域內(nèi)直方圖變得稀疏不穩(wěn)定。并且量化等級(jí)越高，存儲(chǔ)直方圖信息所需的內(nèi)存空間就越大，計(jì)算效率越低。適當(dāng)減少灰度值的量化等級(jí)，可以使直方圖變得緊密而提高穩(wěn)定性，但這樣會(huì)損失一定的光譜區(qū)分能力。因此實(shí)際計(jì)算中需要在區(qū)分能力和穩(wěn)定性之間均衡。在本文中，光譜直方圖被壓縮到32級(jí)?？紤]到影像中灰度值的分布范圍可能較小，因此對(duì)直方圖直接進(jìn)行線性的壓縮可能大大降低其區(qū)分能力。本文中，首先對(duì)原始影像先采用2倍標(biāo)準(zhǔn)差拉伸，然后再進(jìn)行線性的灰度壓縮，將較窄的灰度范圍較均勻的壓縮到32個(gè)灰度級(jí)中。在此基礎(chǔ)上統(tǒng)計(jì)區(qū)域的光譜直方圖，此時(shí)區(qū)域m和n之間的光譜距離定義為兩區(qū)域光譜直方圖的距離

式中，N表示影像的波段數(shù)；GiC（m，n）表示區(qū)域m和n在第i波段的直方圖距離。

2.3 光譜與紋理特征自動(dòng)加權(quán)

光譜特征與紋理特征的組合是實(shí)現(xiàn)光譜紋理分割的關(guān)鍵之一，最為常用的便是加權(quán)組合方式。設(shè)hC和hT分別為光譜與紋理異質(zhì)性，wC、wT為對(duì)應(yīng)的權(quán)重，則光譜與紋理加權(quán)組合形式如下

wC和wT的確定便成為光譜與紋理組合的關(guān)鍵問(wèn)題。文獻(xiàn)［16—17］等針對(duì)遙感影像分割中光譜與紋理特征加權(quán)組合提出了兩種半自動(dòng)方法。這兩種方法都采用相同的原則，即對(duì)于紋理占主導(dǎo)區(qū)域設(shè)置較高紋理權(quán)重，反之設(shè)置較大光譜權(quán)重。這兩種方法以區(qū)域的方差為依據(jù)，設(shè)定經(jīng)驗(yàn)閾值以判定區(qū)域是否屬于紋理主導(dǎo)區(qū)域，然后再根據(jù)判定結(jié)果選擇預(yù)先設(shè)置的權(quán)重。這種經(jīng)驗(yàn)閾值方法對(duì)于影像中紋理特征突出的區(qū)域具有良好的效果，而一些紋理細(xì)密的區(qū)域可能被歸入光譜主導(dǎo)的區(qū)域。同時(shí)，所有判定為紋理主導(dǎo)的區(qū)域所設(shè)置的權(quán)重均是固定的，難以反映不同區(qū)域光譜和紋理特征構(gòu)成成分的差異。

為此，本文采用了一種全自動(dòng)的權(quán)重設(shè)置方法，依據(jù)不同兩區(qū)域的影像內(nèi)容動(dòng)態(tài)調(diào)整光譜和紋理特征的權(quán)重。首先統(tǒng)計(jì)兩近鄰區(qū)域的光譜特征分布情況，然后以光譜特征所表述的信息自動(dòng)調(diào)整紋理與光譜的權(quán)重。對(duì)于紋理粗糙的區(qū)域，由于紋理結(jié)構(gòu)的影響，影像灰度值比較分散或集中于某幾個(gè)值附近，使得直方圖中最大概率值較小。在紋理細(xì)密的區(qū)域，灰度分布較為集中，灰度頻率最大值較高。圖3所示為某真彩色航空影像上林地與農(nóng)田區(qū)域及其光譜直方圖。兩區(qū)域光譜均值相近，但紋理特征存在較大的差異，在光譜直方圖上反映為離散程度和峰值的不同。

由以上分析可知，當(dāng)兩區(qū)域紋理結(jié)構(gòu)均比較細(xì)密時(shí)，兩光譜直方圖最大值均比較高，此時(shí)設(shè)置較大的光譜權(quán)重。反之，只要任一區(qū)域紋理粗糙則需加大紋理權(quán)重。據(jù)此，以兩區(qū)域光譜直方圖情況設(shè)置光譜與紋理權(quán)重。定義

式中，N為波段總數(shù)；maxf（）i表示第i波段的灰度概率最大值；k表示的對(duì)當(dāng)前區(qū)域各波段概率密度最大值的加權(quán)。分別對(duì)相鄰的兩區(qū)域統(tǒng)計(jì)k1和k2。由于光譜直方圖已經(jīng)在一定程度上反映部分紋理信息，因此適當(dāng)增大光譜的權(quán)重，取k1和k2最小值的平方根為光譜權(quán)重

式中，wC為光譜權(quán)重；wT為紋理權(quán)重。當(dāng)相鄰兩區(qū)域中任意一個(gè)區(qū)域由紋理占主導(dǎo)時(shí)，wC將較小，這時(shí)兩個(gè)區(qū)域的區(qū)分主要依靠其紋理的差別。當(dāng)k1和k2均較大時(shí)，兩區(qū)域的紋理特征均不明顯，以光譜信息作為區(qū)分兩區(qū)域的主要特征。

圖3 不同紋理區(qū)域及其光譜直方圖（已壓縮為32級(jí)）Fig.3 Examples of regions with different textures and their histograms（32bins）

2.4 結(jié)合形狀特征的異質(zhì)性度量

遙感影像上地物是具有一定形狀信息的地理實(shí)體，在對(duì)象合并中，僅利用地物的光譜與紋理信息往往使得圖斑輪廓蜿蜒曲折，較難形成具有實(shí)際意義的對(duì)象，加以幾何約束將有助于獲取具有實(shí)際意義的對(duì)象。初始分割往往受噪聲、紋理、細(xì)微邊緣的影響，使得對(duì)象破碎，輪廓復(fù)雜。在合并中，優(yōu)先合并公共邊較長(zhǎng)的相鄰區(qū)域?qū)⒂兄讷@取結(jié)構(gòu)更為緊湊的對(duì)象［20］。以式（11）為基礎(chǔ)，引入公共邊長(zhǎng)度的影響，定義新的異質(zhì)性準(zhǔn)則為

式中，lmn為相鄰區(qū)域m和n的公共邊長(zhǎng)度；λ為公共邊的影響系數(shù)，當(dāng)λ＝0時(shí)，lλmn＝1，表示公共邊不影響區(qū)域異質(zhì)性度量；當(dāng)λ不為0時(shí)，公共邊越長(zhǎng)，異質(zhì)性越小。

2.5 基于鄰接圖的區(qū)域合并

采用基于RAG的合并方法是將初始分割的每一個(gè)區(qū)域均視為圖中的一個(gè)節(jié)點(diǎn)，當(dāng)兩區(qū)域相鄰時(shí)，對(duì)應(yīng)的節(jié)點(diǎn)便以一具有權(quán)重的邊緣進(jìn)行連接［21］，通過(guò)對(duì)圖中節(jié)點(diǎn)的合并來(lái)實(shí)現(xiàn)區(qū)域的合并。文獻(xiàn)［19］提出了一種逐步迭代優(yōu)化的分割模型，該模型將影像分割視為影像逼近問(wèn)題，在每一步的優(yōu)化過(guò)程中力求逼近誤差最小。在基于鄰接圖的區(qū)域合并中，逐步優(yōu)化問(wèn)題轉(zhuǎn)化為尋找合并代價(jià)最小的區(qū)域進(jìn)行合并，合并代價(jià)即對(duì)應(yīng)于圖中邊緣的權(quán)重，其定義為

式中，C（m，n）表示區(qū)域m和n的合并代價(jià)值，即對(duì)整個(gè)圖像逼近問(wèn)題所帶來(lái)誤差；Sm和Sn表示兩區(qū)域的面積；h（m，n）表示兩區(qū)域的異質(zhì)性。將h1（m，n）和h（m，n）逐次代入可得

上式所示合并代價(jià)不僅融合了區(qū)域的光譜、紋理與形狀特征，還充分考慮了區(qū)域的面積特征。分析式（17）可得知，在區(qū)域異質(zhì)性相同的情況下，區(qū)域的面積越小，合并對(duì)于整個(gè)圖像逼近帶來(lái)的誤差也越小，其合并也代價(jià)越小。在區(qū)域合并過(guò)程中優(yōu)先合并面積較小的區(qū)域，形成較大的區(qū)域，還有助于增強(qiáng)區(qū)域光譜和紋理統(tǒng)計(jì)特征的穩(wěn)定性。在鄰接圖中逐次選擇合并代價(jià)最小的節(jié)點(diǎn)進(jìn)行合并，并更新相關(guān)的邊緣權(quán)重，重復(fù)這一過(guò)程，直至合并代價(jià)最小值超過(guò)某一閾值或滿足其他設(shè)定條件時(shí)合并停止，獲得最終分割結(jié)果。

3 試驗(yàn)結(jié)果與分析

為驗(yàn)證本文算法效果，采用1幅TerraSar－X單極化影像和兩幅QuickBird真彩色影像進(jìn)行試驗(yàn)分析。影像原始大小為400像素×400像素，TerraSar－X影像空間分辨率為3m，QuickBird影像為RGB波段組合影像，空間分辨率2.44m。本試驗(yàn)影像分割目的是滿足地表覆蓋分類需求，因此需要盡可能將同一類型的地物分割在一起，如城市、森林、水體等。因區(qū)域合并可以逐步改善過(guò)分割效果，但無(wú)法處理欠分割的情況，因此在初始分割階段需保證不出現(xiàn)欠分割為宜。本文采用分水嶺分割算法獲取初始分割結(jié)果，并將面積小于50的區(qū)域合并到其光譜均值距離最小的相鄰區(qū)域中。分水嶺算法能較好地區(qū)分影像中的同質(zhì)區(qū)域，產(chǎn)生密集的過(guò)分割圖斑，保證后續(xù)處理的精度。

3.1 試驗(yàn)1

為驗(yàn)證本文算法，本試驗(yàn)將試驗(yàn)影像分別以JSEG算法［6］、FNEA 算法［12］進(jìn)行對(duì) 比 試 驗(yàn)。試驗(yàn)平臺(tái)CPU為Intel E7200雙核，2GB內(nèi)存，Windows XP SP3操作系統(tǒng)。JSEG作為光譜－紋理分割的經(jīng)典算法，常用于評(píng)判其他顏色－紋理分割算法的優(yōu)劣，該算法也成功應(yīng)用于遙感影像分割中［6］；而FNEA算法是世界上首個(gè)面向?qū)ο笥跋穹治鲕浖Cognition的核心算法，被廣泛應(yīng)用于各種遙感影像分析中。因而本文選取這兩種標(biāo)志性的算法作為參考衡量本文算法效果。本文算法合并時(shí)設(shè)置形狀參數(shù)為0.5，試驗(yàn)中當(dāng)剩余區(qū)域的個(gè)數(shù)達(dá)到設(shè)定閾值時(shí)停止合并，通過(guò)多組試驗(yàn)，最終設(shè)定SAR影像和QuickBird影像區(qū)域合并的停止閾值均為24。采用JSEG算法和FNEA算法時(shí)，同樣調(diào)整參數(shù)進(jìn)行多次試驗(yàn)，并采用目視解譯的方式選擇整體效果最佳的分割結(jié)果作為對(duì)比。圖4（a）和（d）JSEG參數(shù)分別為：q＝255，l＝2，m＝0.25和q＝379，l＝2，m＝0.4。圖4（b）和（e）FNEA算法參數(shù)分別為shape＝0.4，compactness＝0.55，scale＝140和shape＝0.4，compactness＝0.55，scale＝120。

圖4 分割效果對(duì)比圖Fig.4 Segmentation results of different method

3種分割算法效果如圖4所示，圖中比較明顯的過(guò)分割或欠分割以箭頭指示。圖4（a）、（d）所示為JSEG算法結(jié)果，兩圖中區(qū)域個(gè)數(shù)分別為31和39。JSEG算法針對(duì)遙感影像分割存在一定的局限，對(duì)于紋理粗糙地區(qū)過(guò)分割，而對(duì)于紋理相對(duì)細(xì)密地區(qū)則易出現(xiàn)欠分割。圖4（b）和（e）所示為FNEA算法分割結(jié)果，兩圖中區(qū)域個(gè)數(shù)分別為24和25。FNEA算法采用了光譜特征和形狀特征進(jìn)行分割，但由于其沒(méi)有考慮影像中紋理特征，在紋理密集分布的地方往往出現(xiàn)過(guò)分割現(xiàn)象，并且對(duì)象邊界彎曲復(fù)雜的情況，如圖4（b）、（e）中箭頭處居民地；而對(duì)于一些光譜相近紋理存在一定差別的區(qū)域容易出現(xiàn)欠分割，如圖4（e）中箭頭所指處耕地和水體。圖4（c）、（f）所示為本文算法結(jié)果，兩圖中區(qū)域個(gè)數(shù)均為24。從圖中可以看出本文算法能夠較好度量不同地物的光譜與紋理信息，反映地物之間的差別，較好地實(shí)現(xiàn)不同類型地物的分割。從圖中還可以看出本文算法在能較為精確地獲取不同類型地物之間的邊界，保證后續(xù)分析的精度。

表1給出了3種不同算法針對(duì)兩幅試驗(yàn)圖片的運(yùn)行時(shí)間。由于FNEA算法采用eCognition軟件實(shí)現(xiàn)，該軟件在運(yùn)行時(shí)間大于1s時(shí)統(tǒng)計(jì)精確到秒，因而表1中FNEA算法時(shí)間僅供參考。從表1中可以看出，本文算法的運(yùn)行效率最高。這充分體現(xiàn)了所采用的LBP紋理的高效，同時(shí)也體現(xiàn)了基于區(qū)域鄰接圖的區(qū)域合并算法的高效性。

3.2 試驗(yàn)2

為了分析形狀參數(shù)對(duì)分割結(jié)果的影響，采用一幅QuickBird影像，以分水嶺獲取的初始結(jié)果，再分別采用形狀參數(shù)0、0.5和1進(jìn)行區(qū)域合并，通過(guò)調(diào)整合并停止閾值，使得3組試驗(yàn)結(jié)果中區(qū)域數(shù)目相同。試驗(yàn)結(jié)果如圖5所示。

表1 分割時(shí)間比較Tab.1 Running times of different algorithms for the two images s

圖5 不同形狀參數(shù)下分割結(jié)果及局部放大圖像Fig.5 Segmentation results using different shape parameters

從對(duì)比3組不同形狀參數(shù)下分割結(jié)果可以看出，當(dāng)形狀參數(shù)為0時(shí)，區(qū)域合并以對(duì)象的光譜和紋理特征為主導(dǎo)。城區(qū)中紋理復(fù)雜，且存在部分植被干擾，被分割為多個(gè)小塊，并且各個(gè)對(duì)象區(qū)域之間輪廓線較為復(fù)雜。隨著形狀參數(shù)的增大，對(duì)象公共邊長(zhǎng)度對(duì)區(qū)域異質(zhì)性的影響逐步加大，公共邊較長(zhǎng)的對(duì)象被合并在一起的優(yōu)先級(jí)升高，從而得到的對(duì)象輪廓更為緊致。如圖5（b）、（c）所示，伴隨形狀參數(shù)逐步增大，城區(qū)里原本邊界復(fù)雜的對(duì)象由于公共邊較長(zhǎng)而被優(yōu)先合并，城區(qū)分割效果逐步改善。但是當(dāng)形狀參數(shù)過(guò)大時(shí)，可能會(huì)導(dǎo)致某些線狀地物與其相鄰地物合并，引起欠分割，如圖5（f）中箭頭所指處植被，被合并到水體中。

4 結(jié) 論

本文方法實(shí)現(xiàn)了光譜與紋理信息的自動(dòng)加權(quán)組合，并以區(qū)域公共邊進(jìn)行空間約束進(jìn)行影像分割，為多特征結(jié)合的遙感影像分割提供了一條有效途徑。試驗(yàn)結(jié)果表明本文方法能結(jié)合光譜、紋理與形狀結(jié)構(gòu)信息，改善分割質(zhì)量，特別是紋理信息較為豐富的區(qū)域。與其他算法相比，具有分割結(jié)果準(zhǔn)確，效率高等優(yōu)點(diǎn)，具有較高的實(shí)用價(jià)值。為進(jìn)一步提高遙感影像分割的自動(dòng)化、智能化水平，在本文方法的基礎(chǔ)上還可進(jìn)一步引入其他特征（如邊緣、上下文關(guān)系等），同時(shí)，在區(qū)域合并停止閾值方面也有待進(jìn)一步的研究。

［1］ XIAO Pengfeng，F(xiàn)ENG Xuezhi，ZHAO Shuhe，et al.Segmentation of High－resolution Remote Sensed Imagery Based on Phase Congruency［J］.Acta Geodaetica et Cartographica Sinica，2007，36（2）：146－151.（肖鵬峰，馮學(xué)智，趙書河，等.基于相位一致的高分辨率遙感圖像分割方法［J］.測(cè)繪學(xué)報(bào)，2007，36（2）：146－151.）

［2］ GAETANO R，SCARPA G，POGGI G.Hierarchical Texture－based Segmentation of Multiresolution Remote Sensing Images［J］.IEEE Transactions on Geoscience and Remote Sensing，2009，47（7）：2129－2141.

［3］ OJALA T，PIETIKAINEN M.Unsupervised Texture Segmentation Using Feature Distribution［J］.Pattern Recognition，1999，32：477－486.

［4］ YU Peng，ZHANG Zhenlong，HOU Zhiqun.Textured Image Segmentation Based on Gaussian Markov Random Field Mixture Model［J］.Acta Geodaetica et Cartographica Sinica，2006，35（8）：224－228.（余鵬，張震龍，侯至群.基于高斯馬爾科夫隨機(jī)場(chǎng)混合模型的紋理圖像分割［J］.測(cè)繪學(xué)報(bào)，2006，35（8）：224－228.）

［5］ TRIAS－SANZ R，STAMON G，LOUCHET J.Using Colour，Texture，and Hierarchical Segmentation for High－resolution Remote Sensing［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2008，63：156－168.

［6］ DENG Y，MANJUNATH B S.Unsupervised Segmentation of Color－texture Regions in Image and Video［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2001，23（8）：800－810.

［7］ LIU T T，ZHANG L P，LI P X，et al.Region Segmentation of Multi－spectral Remote Sensing Images Using Improved JSEG Algorithm［J］.Journal of Remote Sensing，2009，13（1）：30－34.

［8］ CARSON C，BELONGIE S，GREENSPAN H，et al.Blobworld：Image Segmentation Using Expectation－maximization and Its Application to Image Querying［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2002，24（8）：1026－1038.

［9］ PERMUTER H，F(xiàn)RANCOS J，JERMYN I.A Study of Gaussian Mixture Models of Color and Texture Features for Image Classification and Segmentation［J］.Pattern Recognition，2006，39：695－706.

［10］ KATO Z，PONG T C.A Markov Random Field Image Segmentation Model for Colour Texture Images［J］.Image and Vision Computing，2006，24：1103－1114.

［11］ XIA Y，F(xiàn)ENG D，ZHAO R.Adaptive Segmentation of Texture Images by Using the Coupled Markov Random Field Model［J］.IEEE Transactions on Image Processing，2006，15（11）：3559－3566.

［12］ BAATZ M，SCHAPE A.Multiresolution Segmentation：An Optimization Approach for High Quality Multiscale Image Segmentation［J］.Journal of Photogrammetry and Remote Sensing，2000，58（3－4）：12－23.

［13］ OJALA T，PIETIKAINEN M，HARWOOD D.A Comparative Study of Texture Measures with Classification Based on Feature Distributions［J］.Pattern Recognition，1996，29（1）：51－59.

［14］ OJALA T，PIETIKAINEN M，MAENPAA T.Multiresolution Gray－scale and Rotation Invariant Texture Classification with Local Binary Patterns［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2002，24（7）：971－987.

［15］ CHEN K M，CHEN S Y.Color Texture Segmentation Using Feature Distributions［J］.Pattern Recognition Letters，2002，23：755－771.

［16］ HU X，VINCENT T，BJORN P.Automatic Segmentation of High－resolution Satellite Imagery by Integrating Texture，Intensity，and Color Features［J］.Photogrammetric Engineering and Remote Sensing，2004，71：1399－1406.

［17］ WANG A，WANG S，LUCIEER A.Segmentation of Multispectral High－resolution Satellite Imagery Based on Integrated Feature Distributions［J］.International Journal of Remote Sensing，2010，31（6）：1471－1483.

［18］ ZHENG Yongbin，HUANG Xinsheng，F(xiàn)ENG Songjiang.An Image Match Algorithm Based on Combination of SIFT and the Rotation Invariant LBP［J］.Journal of Computer－Aided Design ＆ Computer Graphics，2010，22（2）：286－292.（鄭永斌，黃新生，豐松江.SIFT和旋轉(zhuǎn)不變LBP相結(jié)合的圖像匹配算法［J］.計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2010，22（2）：286－292.）

［19］ BEAULIEU J M，GOLDBERG M.Hierarchy in Picture Segmentation：A Stepwise Optimization Approach ［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence.1989，11（2）：150－163.

［20］ TAN Yumin，HUAI Jianzhu，TANG Zhongshi.Edge－guided Segmentation Method for Multiscale and High Resolution Remote Sensing Image［J］.Journal of Infrared and Millimeter Waves，2010，29（4）：312－315.（譚玉敏，槐建柱，唐中實(shí).一種邊界引導(dǎo)的多尺度高分辨遙感圖像分割方法［J］.紅外與毫米波學(xué)報(bào)，2010，29（4）：312－315.）

［21］ TAN Yumin，HUAI Jianzhu，TANG Zhongshi.Objectoriented Adjacency Graph Partition Algorithm for Remote Sensing Image Segmentation［J］.Journal of Dalian Maritime University，2009，35（2）：81－83.（譚玉敏，槐建柱，唐中實(shí).基于鄰接圖的面向?qū)ο筮b感圖像分割算法［J］.大連海事大學(xué)學(xué)報(bào)，2009，35（2）：81－83.）