面狀居民地匹配骨架線傅里葉變化方法

黃智深，錢海忠，郭 敏，3，劉海龍，王 驍

1.信息工程大學(xué) 地理空間信息學(xué)院，河南 鄭州 450052；2.南疆測勤隊，新疆 喀什 844200；3.61175部隊，江蘇 南京 210049

1 引 言

各級應(yīng)用部門經(jīng)常需要將多源空間數(shù)據(jù)進(jìn)行集成與融合，增強(qiáng)其共享性，避免數(shù)據(jù)的重復(fù)采集，以節(jié)省人力物力。多源空間數(shù)據(jù)之間往往存在差異性［1］，空間數(shù)據(jù)匹配作為空間數(shù)據(jù)更新與融合的關(guān)鍵技術(shù)之一，越來越受到重視。

城市居民地是空間目標(biāo)變化中最為活躍的要素之一，也是空間同名實體匹配中復(fù)雜性最強(qiáng)、最具挑戰(zhàn)性的內(nèi)容之一。目前，針對大比例尺城市面狀居民地匹配的研究已有很多，例如基于面狀居民地重疊面積的匹配［2］、基于面實體幾何形狀相似性的匹配［3－4，19］、基于知識的非空間屬性數(shù)據(jù)通過計算屬性項的相似度值而進(jìn)行的匹配［5］、基于面質(zhì)心進(jìn)行粗匹配并結(jié)合多種匹配檢驗規(guī)則（如面積、面密度等）進(jìn)行最終匹配［6］、基于模糊拓?fù)潢P(guān)系分類的匹配［7］、利用面實體重心距離與重疊面積的基于概率的匹配［8］、利用多級弦長函數(shù)以及中心距離函數(shù)對要素幾何形狀進(jìn)行多級描述并建立相似性度量模型進(jìn)行匹配等［18］。居民地多為不規(guī)則圖形，形態(tài)復(fù)雜，主要依據(jù)其面積重疊率、大小、方向等進(jìn)行匹配，不但算法復(fù)雜，且描述結(jié)果存在許多不確定性；同時，居民地作為大比例尺城市地圖的主要組成部分，其數(shù)據(jù)量大，直接對其進(jìn)行匹配操作，將導(dǎo)致整個匹配過程計算量大，匹配效率低等問題。

為解決大比例尺面狀居民地匹配過程中數(shù)據(jù)量大、形態(tài)復(fù)雜、不確定性強(qiáng)等問題，首先對更新前面狀居民地和更新后面狀居民地分別提取能反映該居民地輪廓形態(tài)主要特征的骨架線；然后通過對雙方骨架線進(jìn)行傅里葉變換與分析，來實現(xiàn)多源居民地之間的匹配。該方法的特點(diǎn)在于：①把面狀居民地轉(zhuǎn)化為能夠反映其輪廓主要形態(tài)特征的骨架線，一方面數(shù)據(jù)復(fù)雜性降低了，另一方面可以引入線要素的許多匹配算法，拓展了面要素匹配的技術(shù)途徑；② 傅里葉變換的變換系數(shù)可較好地描述線狀要素的形態(tài)特征，并將形態(tài)信息從空間域轉(zhuǎn)化到頻率域，特別對幾何形態(tài)相似性判斷具備較好的識別和區(qū)分能力，可很好地提高匹配的準(zhǔn)確率。

同時，由于把面狀居民地轉(zhuǎn)化為骨架線進(jìn)行匹配，研究對象作如下限定：① 匹配雙方為相同或相近比例尺面狀居民地；② 針對許多數(shù)據(jù)語義信息不完整的現(xiàn)實，主要研究幾何匹配，暫不考慮語義信息；③ 由于面狀居民地與其骨架線之間一一對應(yīng)，因此本文只針對1∶1居民地匹配情形展開研究，而針對1∶n和m∶n（這實際上是跨比例尺之間匹配）等匹配情形暫不考慮，這是因為試驗表明目前對多個居民地骨架線進(jìn)行合并，其結(jié)果難以反映這些居民地的形態(tài)分布與輪廓特征，后續(xù)擬在本文研究基礎(chǔ)上，通過空間拓?fù)潢P(guān)系約束展開進(jìn)一步研究。

2 基于骨架線傅里葉變換的面狀居民地匹配算法

本算法主要由以下4個步驟組成：

（1）面狀居民地主骨架線提取。對匹配雙方的居民地集分別提取骨架線，將2維居民地轉(zhuǎn)化為1維骨架線，降低數(shù)據(jù)的復(fù)雜性，提高計算效率。

（2）骨架線插值計算。包括插值數(shù)據(jù)準(zhǔn)備和插值計算實施。前者量化描述骨架線的幾何形態(tài)特征；后者對前者的量化描述進(jìn)行插值計算，為傅里葉變換打下基礎(chǔ)。

（3）離散傅里葉變換。對插值計算后的骨架線幾何形態(tài)數(shù)據(jù)進(jìn)行離散傅里葉變換，將骨架線形態(tài)信息從空間域轉(zhuǎn)到頻率域，更加充分地反映骨架線的幾何形態(tài)信息。

（4）計算匹配雙方的相關(guān)系數(shù)，并據(jù)此進(jìn)行匹配。根據(jù)傅里葉變換結(jié)果，計算更新前數(shù)據(jù)骨架線和更新后數(shù)據(jù)骨架線之間的相關(guān)系數(shù)，據(jù)此來確定面狀居民地之間的匹配關(guān)系。下文將詳細(xì)闡述。

2.1 居民地骨架線提取

目前提取骨架線的常見算法有：平行線切割中點(diǎn)連線法、內(nèi)側(cè)緩沖區(qū)法以及Delaunay三角網(wǎng)法［9－10］。平行線切割中點(diǎn)連線法簡單、直觀，但對復(fù)雜多邊形（如含島），往往得不到理想結(jié)果。內(nèi)測緩沖區(qū)法原理為：對多邊形內(nèi)側(cè)迭代做緩沖區(qū)，直至獲取小于面積要求的內(nèi)側(cè)緩沖區(qū)，將該緩沖區(qū)形心連線作為骨架線［12］，該方法與第一種方法存在相似的問題。Delaunay三角網(wǎng)法具有很好的幾何特性，能夠方便建立起不同空間目標(biāo)的空間鄰近關(guān)系［11］，并能較為詳細(xì)地反映面要素的幾何形態(tài)特征。

綜合考慮，本文選擇Delaunay方法來提取面要素骨架線。圖1是采用Delaunay方法提取面要素骨架線的原理圖，圖2是一個提取居民地主骨架線的例子。可以看出該方法能很好地反映居民地輪廓的主要形態(tài)特征，同時以主骨架線來代替居民地進(jìn)行匹配，可有效降低居民地匹配的復(fù)雜度。

圖1 采用Delaunay方法提取面要素骨架線原理圖Fig.1 Pick up skeleton－line with Delaunay technique

圖2 面狀居民地輪廓及其主骨架線提取例子Fig.2 Example of a polygon habitation and its main skeleton line

2.2 骨架線插值計算

為描述骨架線的幾何形態(tài)，把骨架線的每個直線段轉(zhuǎn)換為長度和角度信息，即骨架線每個直線段都轉(zhuǎn)化為一對坐標(biāo)值（長度，角度），這種坐標(biāo)點(diǎn)稱為幾何形態(tài)數(shù)據(jù)點(diǎn)。同時，原有骨架線幾何形態(tài)數(shù)據(jù)點(diǎn)的數(shù)量往往不能滿足匹配分析需求，需要對其進(jìn)行插值計算。詳細(xì)方法如下。

2.2.1 骨架線幾何形態(tài)描述

骨架線幾何形態(tài)描述主要包括骨架線起始點(diǎn)選擇、骨架線長度特征描述、骨架線角度特征描述以及幾何形態(tài)數(shù)據(jù)點(diǎn)生成4個部分。

2.2.1.1 骨架線起始點(diǎn)選取

進(jìn)行骨架線長度指標(biāo)和角度指標(biāo)描述之前，需首先確定骨架線上的某一頂點(diǎn)為起始點(diǎn)P0，選擇不同起始點(diǎn)P0，將會得到不同的幾何形態(tài)描述數(shù)據(jù)。本文對起始點(diǎn)P0選取的依據(jù)為

即選擇橫坐標(biāo)最小的節(jié)點(diǎn)為起始點(diǎn)，且設(shè)定角度特征描述中以順時針方向為正。其優(yōu)點(diǎn)在于：可以確保骨架線的長度特征和角度特征均為正，且分析結(jié)果可視化圖形位于坐標(biāo)系的第一象限內(nèi)。

2.2.1.2 骨架線長度特征描述

對于骨架線長度特征，采用骨架線上各相鄰節(jié)點(diǎn)之間的長度在骨架線全長中所占比率的形式進(jìn)行描述。因為所有節(jié)點(diǎn)距之和與骨架線全長比值為1，所以該比率值又稱為歸一化距離［3］。其描述步驟為：① 從骨架線起始點(diǎn)開始，計算骨架線的全長及其任意兩個相鄰節(jié)點(diǎn)（設(shè)起點(diǎn)、終點(diǎn)也為節(jié)點(diǎn)）之間的長度；② 計算相鄰節(jié)點(diǎn)距離在骨架線全長中所占比率，用長度比率值代替相鄰節(jié)點(diǎn)間的實際長度值。具體實現(xiàn)公式骨架線全長

相鄰節(jié)點(diǎn)距離所占比率

這里需注意兩個方面：① 為方便后續(xù)數(shù)據(jù)的運(yùn)算及可視化，需將長度比率數(shù)組的第一個數(shù)據(jù)賦0；②前面各段的長度比率值需在后面的比率數(shù)據(jù)中累加，最后一個線段的長度比率值應(yīng)為1。

2.2.1.3 骨架線角度特征描述

對于骨架線角度特征，采用骨架線上各線段的實際坐標(biāo)方位角進(jìn)行描述。具體步驟為：① 從骨架線起始點(diǎn)出發(fā)，求得各線段的轉(zhuǎn)角φk（以順時針方向為正）；② 根據(jù)所求的轉(zhuǎn)角φk，計算骨架線上各線段的坐標(biāo)方位角θk（以順時針方向為正）。轉(zhuǎn)角φk計算公式

方位角計算公式：Δy＝y(tǒng)i＋1－yi，Δx＝xi＋1－xi

上述計算保證了θk均為正數(shù)。由于角度特征數(shù)據(jù)與骨架線上線段一一對應(yīng)，即使在一條骨架線的角度特征描述過程中出現(xiàn)多次θk值相同的情況，但該值的計算存在先后關(guān)系，可以避免一對多情形的出現(xiàn)。另外，在角度特征數(shù)據(jù)存儲時，從骨架線起始點(diǎn)出發(fā)，骨架線第一個線段的角度特征描述數(shù)據(jù)需進(jìn)行兩次記錄，即將其同時存儲在角度特征數(shù)組的第1位和第2位。這是因為，對于骨架線為直線這種特殊情況，其只有首端點(diǎn)和末端點(diǎn)，首端點(diǎn)和末端點(diǎn)之間只能計算一個方位角，該方位角與長度之間只能構(gòu)成一個幾何形態(tài)數(shù)據(jù)點(diǎn)，不能滿足對直線骨架線的表達(dá)。因此需要將骨架線角度特征數(shù)組的第2個數(shù)據(jù)賦予與第1個數(shù)據(jù)相同的值。

2.2.1.4 骨架線幾何形態(tài)數(shù)據(jù)點(diǎn)的生成

骨架線幾何形態(tài)數(shù)據(jù)點(diǎn)由長度特征數(shù)據(jù)和角度特征數(shù)據(jù)共同組成。由于這兩組數(shù)據(jù)數(shù)量相同，且生成時的起始點(diǎn)相同，因此可按照先后順序，以長度特征為橫軸，以角度特征為縱軸，一一對應(yīng)，得到一系列有序的坐標(biāo)對，這些坐標(biāo)對即為骨架線的幾何形態(tài)數(shù)據(jù)點(diǎn)。例如，圖3中骨架線A的幾何形態(tài)描述數(shù)據(jù)如表1所示。

表1 骨架線A幾何形態(tài)描述數(shù)據(jù)結(jié)果表Tab.1 The data describing the geometry of skeleton line Ashowing

圖3 骨架線A示意圖Fig.3 The diagram of skeleton line A

2.2.2 幾何形態(tài)數(shù)據(jù)點(diǎn)插值計算

如果幾何形態(tài)數(shù)據(jù)點(diǎn)數(shù)量過少，則不能滿足相似性分析需求，這就需要對其插值，進(jìn)一步突出骨架線主要幾何形態(tài)特征，特別是整體形態(tài)特征。為避免出現(xiàn)新的誤差，對待匹配雙方骨架線必須采用相同插值方法進(jìn)行插值。常用插值方法有拉格朗日插值、厄爾米特插值法以及樣條插值法等［13］。

拉格朗日插值法是高次多項式插值，插值結(jié)果光滑，但不能保證插值收斂性，容易引起新的誤差或不確定性；厄爾米特插值法屬于低次多項式插值，使用簡單，但是插值后曲線過于粗糙；三次樣條插值法也屬于低次多項式插值，但較厄爾米特插值法有較大改進(jìn)，可有效確保插值后骨架線的準(zhǔn)確性［14］。本文采用三次樣條插值法，方法如下［15］。

由于長度特征數(shù)據(jù)采用歸一化距離，其值域為x∈［0，1］，插值時首先將［0，1］區(qū)間進(jìn)行等分，然后以區(qū)間等分點(diǎn)處的坐標(biāo)值為插值點(diǎn)，以骨架線幾何形態(tài)數(shù)據(jù)點(diǎn)為樣條節(jié)點(diǎn)，進(jìn)行插值運(yùn)算。

在每個［xi－1，xi］區(qū)間上，三次樣條函數(shù)可以表示為以下形式

式中，hi－1＝xi－xi－1；Mi＝Sn（xi）。

Mi所滿足的方程

式中

邊界條件有如下兩種：

第一型插值條件

由此導(dǎo)出

自然樣條

由此得出M1＝0、Mn＝0，從而得出Mi滿足三對角方程。

以圖4中所示的待匹配骨架線為例，其幾何形態(tài)特征數(shù)據(jù)點(diǎn)的正切空間圖形如圖5所示，對圖5中的兩條待匹配骨架線的幾何形態(tài)特征圖形進(jìn)行三次樣條插值分析，將（0，1）區(qū)間等分成50份，如圖6所示。對比圖5和圖6可知，三次樣條插值后，骨架線幾何形態(tài)數(shù)據(jù)由直方圖轉(zhuǎn)化為曲線，不但豐富了幾何形態(tài)數(shù)據(jù)點(diǎn)，而且骨架線幾何特征的整體趨勢以及局部差別也更為直觀。

圖4 某更新前及更新后居民地骨架線ig.4 The skeleton lines’before and after updating

圖5 骨架線幾何形態(tài)數(shù)據(jù)的可視化結(jié)果Fig.5 Visualization of skeleton lines geometry data

圖6 對圖5進(jìn)行三次樣條插值后的結(jié)果Fig.6 The result of cubic spline interpolation to fig.5

2.3 離散傅里葉變換

對經(jīng)過三次樣條插值的骨架線幾何形態(tài)數(shù)據(jù)進(jìn)行離散傅里葉變換，其主要變換過程［15］如下。

令ωn＝e－i2π／n，對于N個離散數(shù)據(jù)點(diǎn)xk（k＝0，1，…，N－1），進(jìn)行如下變換

經(jīng)過離散傅里葉變換，將各個離散數(shù)據(jù)點(diǎn)xk分解為它的基本頻率的組合，也就是變換之后y的各個分量。為提高算法效率，采取離散傅里葉變換中的快速傅里葉變換，其變換過程如下

骨架線幾何形態(tài)數(shù)據(jù)經(jīng)過離散傅里葉變換，其結(jié)果由實數(shù)和虛數(shù)兩部分組成。圖6的插值結(jié)果經(jīng)離散傅里葉變換，其實數(shù)和虛數(shù)部分如圖7所示。

圖7 骨架線離散傅里葉變換后數(shù)據(jù)圖像Fig.7 The graphies of the skeleton line discreted with Fourier transform

2.4 相關(guān)系數(shù)計算

傅里葉變換后的骨架線幾何形態(tài)特征數(shù)據(jù)，只是將骨架線信息從空間域轉(zhuǎn)到頻率域，還不能反映骨架線之間的差異。為量化描述骨架線之間在頻率域的相似性，進(jìn)一步采用相關(guān)系數(shù)進(jìn)行描述。因為幾何形態(tài)特征數(shù)據(jù)為離散數(shù)據(jù)，因此采用離散相關(guān)系數(shù)進(jìn)行表達(dá)。離散相關(guān)系數(shù)的一般定義如下。

設(shè)｛gk｝和｛hk｝為兩個未經(jīng)傅里葉變換的離散數(shù)組，k＝1，2，…，n，則｛gk｝與｛hk｝的離散相關(guān)系數(shù)Cgh的一般計算公式如下［17］

在離散傅里葉變換的基礎(chǔ)上求取相關(guān)系數(shù)，根據(jù)相關(guān)離散定理［15］，如果｛gk｝和｛hk｝經(jīng)過離散傅里葉變換后分別為｛Gk｝、｛Hk｝，則Cgh等于｛Gk｝·｛Hk｝＊的逆離散傅里葉變換，其中星號表示復(fù)共軛。其中，逆離散傅里葉變換的一般過程如下

根據(jù)上述分析結(jié)果，設(shè)未經(jīng)傅里葉變換之前的待匹配雙方骨架線的幾何形態(tài)描述數(shù)據(jù)分別為｛gk｝和｛hk｝，對其進(jìn)行離散傅里葉變換的結(jié)果分別為｛Gk｝、｛Hk｝。令DFT和DFT－1分別代表離散傅里葉變換和逆離散傅里葉變換，則

令Z＝｛G｝·｛H｝＊，其中｛H｝＊為表示｛H｝的復(fù)共軛。對Z進(jìn)行逆傅里葉變換，其結(jié)果即為離散相關(guān)系數(shù)，相關(guān)公式如下

相關(guān)系數(shù)的值域為［－1，1］區(qū)間，當(dāng)Cgh＞0時，表示兩要素同向相關(guān)，否則表示異向相關(guān)，Cgh的絕對值越接近1，表示兩數(shù)組的相似性越高；越接近0，表示兩數(shù)組的相似性越低。由于在進(jìn)行骨架線幾何形態(tài)描述時統(tǒng)一了骨架線描述的起始點(diǎn)，從而保證了待匹配雙方骨架線幾何形態(tài)描述數(shù)據(jù)的同向性，所以離散相關(guān)系數(shù)的取值區(qū)間為［0，1］。

3 試驗驗證及分析

3.1 試驗驗證

以某城市同一區(qū)域1∶1萬的更新前面狀居民地和更新后面狀居民地作為試驗數(shù)據(jù)，進(jìn)行算法驗證。詳細(xì)匹配流程如圖8所示。

圖8 基于骨架線傅里葉變換的居民地匹配流程圖Fig.8 The flowchart of habitation matching based on their skeleton lines by the Fourier transform

（1）面狀居民地骨架線提?。禾崛【用竦毓羌芫€［16］，圖9為更新前居民地及其骨架線，圖10為更新后居民地及其骨架線。

（2）插值計算：分別對圖9和圖10中的骨架線進(jìn)行幾何形態(tài)描述，把幾何坐標(biāo)轉(zhuǎn)化為“長度－角度”坐標(biāo)系下的幾何形態(tài)數(shù)據(jù)，然后對其采用三次樣條插值法進(jìn)行插值計算。

圖9 更新前居民地及對其提取骨架線的結(jié)果Fig.9 Block data before the update，and the result of getting their sketch lines

圖10 更新后居民地及對其提取骨架線的結(jié)果Fig.10 Updated block data，and the result of getting their sketch lines

（3）傅里葉變換及相關(guān)系數(shù)計算：根據(jù)插值結(jié)果，分別對更新前居民地骨架線和更新后居民地骨架線的幾何形態(tài)數(shù)據(jù)進(jìn)行離散傅里葉變換，并計算相關(guān)系數(shù)。將相關(guān)系數(shù)與指定閾值比較，若相關(guān)系數(shù)大于指定閾值，則認(rèn)為雙方骨架線之間幾何形態(tài)吻合，滿足匹配條件，進(jìn)而認(rèn)定骨架線所對應(yīng)的居民地之間匹配成功；反之，匹配失敗。

這里，相關(guān)系數(shù)閾值的設(shè)定對匹配結(jié)果有直接影響。相關(guān)系數(shù)取值區(qū)間為［0，1］，相關(guān)系數(shù)越接近1，表明待匹配雙方骨架線的幾何形態(tài)差異越?。环粗浇咏?，說明待匹配雙方骨架線的幾何形態(tài)差異越大。相關(guān)系數(shù)閾值是待匹配雙方骨架線基于傅里葉變換的幾何形態(tài)相似程度標(biāo)準(zhǔn)。例如設(shè)置閾值為0.9時，本算法的各類匹配信息如表2所示（該閾值大小的設(shè)置可根據(jù)實際情況有所不同）。

相關(guān)系數(shù)閾值為0.9時，匹配結(jié)果如下：共有81塊更新前居民地參與匹配，成功匹配71塊，無法匹配10塊（其中已拆除居民地8塊）。共有77塊更新后居民地參與匹配，成功匹配71塊，無法匹配6塊（其中新增居民地4塊）。更新前居民地匹配總成功率為：87.7%（71／81），實際參與匹配成功率為：97.3%（71／（81－8））；更新后居民地匹配總成功率為：92.2%（71／77），實際參與匹配成功率為：97.26%（71／（77－4））。試驗表明本方法具有較高的匹配正確率。

表2 更新前居民地和更新后居民地匹配的各類信息計算結(jié)果（局部）Tab.2 Statistic matching information of block data before update and updated habitation（part）

3.2 試驗對比分析

面要素匹配方法大致可分為兩類［21］，一類是簡單形狀描述法，如包圍盒法、弦長面積描述法、面積周長比法等；另一類是詳細(xì)輪廓描述法，如矩描述子、小波變換法等。前一類方法借助面要素輪廓的某些特征進(jìn)行概略描述，沒有區(qū)分出主要形態(tài)信息和次要形態(tài)信息；后一類方法可非常詳細(xì)地反映居民地輪廓信息，但對細(xì)節(jié)信息過于敏感，不但匹配時計算量較大、匹配過程復(fù)雜，而且存在不確定性。相對已有算法而言，本方法的優(yōu)勢在于：重點(diǎn)對居民地的主要形態(tài)信息進(jìn)行分析，有效降低了匹配過程中的不確定性，提高了匹配準(zhǔn)確率。

為進(jìn)一步驗證本方法科學(xué)性，采用同一數(shù)據(jù)、不同方法的匹配結(jié)果進(jìn)行比較（仍以圖9、10為例）。圖11是本方法匹配結(jié)果，圖12是骨架線緩沖區(qū)法匹配結(jié)果；圖13是矩描述子法匹配結(jié)果，圖14是弦長面積法匹配結(jié)果。圖11至圖14中，高亮色顯示的居民地為未匹配成功的居民地。圖11和圖12中居民地內(nèi)的折線為提取的骨架線。

圖11 采用本方法進(jìn)行居民地匹配結(jié)果Fig.11 Matching result of habitations using the method termed in this paper

圖12 采用骨架線緩沖區(qū)法進(jìn)行居民地匹配結(jié)果Fig.12 Matching result of habitations by the method of skeleton lines′buffers

圖13 采用矩描述子法進(jìn)行居民地匹配結(jié)果Fig.13 Matching result of habitations with moment describe method

圖14 采用弦長面積法進(jìn)行居民地匹配結(jié)果Fig.14 Matching result of habitations by chord length and area method

不同匹配結(jié)果對比如表3所示。從表3可知，本方法匹配成功率優(yōu)于另外3種匹配算法，進(jìn)一步證明了本方法的科學(xué)性和優(yōu)越性。同時，圖12采用的方法為骨架線緩沖區(qū)匹配方法，該方法實際上是本文方法的前期研究成果，即對面狀居民主骨架線進(jìn)行緩沖區(qū)分析，并依據(jù)骨架線緩沖區(qū)之間的面積疊置率來判斷是否匹配成功；其缺陷在于把骨架線還原為緩沖區(qū)后，骨架線的形態(tài)特征被弱化，不確定性增強(qiáng)了，影響了匹配正確率。而本方法對居民地主骨架線采用傅里葉變換的變換系數(shù)來表達(dá)，并將形態(tài)信息從空間域轉(zhuǎn)化到頻率域，特別對幾何形態(tài)相似性判斷具備較好的識別和區(qū)分能力，從而有效提高了匹配準(zhǔn)確率。

4 結(jié) 論

本文提出了一種大比例尺城市面狀居民地匹配方法，對同一區(qū)域大比例尺異源居民地進(jìn)行等級化處理，得到能夠反映面狀居民地主要形態(tài)特征的骨架線，然后對骨架線進(jìn)行幾何形態(tài)描述、插值計算、傅里葉變換以及相關(guān)系數(shù)計算等操作，并依據(jù)相關(guān)系數(shù)的大小對骨架線進(jìn)行匹配判斷，最終根據(jù)骨架線的匹配完成面狀居民地的匹配。該方法有效降低匹配數(shù)據(jù)的復(fù)雜性和可能產(chǎn)生的不確定性因素，提高幾何形態(tài)相似性的識別與區(qū)分能力，同時匹配速度、正確率也得到提高。

表3 采用不同方法居民地匹配情況對比表Tab.3 Comparing results of matching habitations with different matching methods

［1］ GUO Li.Theory and Method Research on Multi－sources Geospatial Vector Data Fusion［D］.Zhengzhou：Information Engineering University，2008：51－52.（郭黎.多源地理空間矢量數(shù)據(jù)融合理論與方法研究［D］.鄭州：信息工程大學(xué)，2008：51－52.）

［2］ GUO Li，ZHENG Haiying，WANG Hao.Study for Area Feature Matching Technique Based on Area Similarity［J］.Hydrographic Surveying and Charting，2009，29（3）：12－15.（郭黎，鄭海鷹，王豪.面狀矢量空間數(shù)據(jù)匹配技術(shù)研究［J］.海洋測繪，2009，29（3）：12－15.）

［3］ FU Zhongliang，SHAO Shiwei，TONG Chunya.Multi－scale Area Entity Shape Matching Based on Tangent Space［J］.Computer Engineering，2010，36（17）：216－220.（付仲良，邵世維，童春芽.基于正切空間的多尺度面實體形狀匹配［J］.計算機(jī)工程，2010，36（17）：216－220.）

［4］ HAO Yanling，TANG Wenjing，ZHAO Yuxin，et al.Areal Feature Matching Algorithm Based on Spatial Similarity［J］.Acta Geodaetica et Cartographica Sinica，2008，27（3）：501－506.（郝燕玲，唐文靜，趙玉新，等.基于空間相似性的面實 體 匹 配 算 法 研 究 ［J］.測 繪 學(xué) 報，2008，27（3）：501－506.）

［5］ COBB M，CHUNG M，F(xiàn)OLEY H.A Rule－based Approach for the Conflation of Attributed Vector Data ［J］.Geo Informatica，1998，2（1）：7－35.

［6］ YUAN S，TAO C.Development of Conflation Components［C］∥Proceedings of Geo－informatics‘99Conference.Ann Arbor：［s.n.］，1999：1－13.

［7］ ZHANG Qiaoping，LI Deren，GONG Jianya.A Real Feature Matching among Urban Geographic Databases［J］.Journal of Remote Sensing，2004，8（2）：107－112.（張橋平，李德仁，龔健雅.城市地圖數(shù)據(jù)庫面實體匹配技術(shù)［J］.遙感學(xué)報，2004，8（2）：107－112.）

［8］ TONG Xiaohua，DENG Susu，SHI Wenzhong.A Probabilistic Theory－based Matching Method ［J］.Acta Geodaetica et Cartographica Sinica.2007，36（2）：210－217.（童小華，鄧愫愫，史文中.基于概率的地圖實體匹配方法［J］.測繪學(xué)報，2007，36（2）：210－217.）

［9］ CAI Mengyi，TIAN Desen.Newly Organized Cartography Tutorial［M］.Beijing：Higher Education Press，2000：109－110.（蔡孟裔，田德森.新編地圖學(xué)教程［M］.北京：高等教育出版社，2000：109－110.）

［10］ HUANG Zhishen，QIAN Haizhong.Dimension Decreaseoriented Habitation Matching Method ［J］.Journal of Geomatics Science and Technology，2012，29（1）：75－78.（黃智深，錢海忠.基于降維技術(shù)的面狀居民地匹配方法［J］.測繪科學(xué)技術(shù)學(xué)報，2012，29（1）：75－78.）

［11］ QIAN Haizhong，WU Fang，ZHU Kunpeng，et al.A Generalization Method of Street Block Based on Dimensionreducing Technique［J］.Acta Geodaetica et Cartographica Sinica，2007，36（1）：102－108.（錢海忠，武芳，朱鯤鵬，等.一種基于降維技術(shù)的街區(qū)綜合方法［J］.測繪學(xué)報，2007，36（1）：102－108.）

［12］ LIU Xiufang，YANG Yongping，LUO Ji，et al.Research on Extracting Polygon Skeleton by Using Inner Buffering Algorithm ［J］.Hydrographic Surveying and Charting，2010，30（5）：46－48.（劉秀芳，楊永平，羅吉，等.基于內(nèi)側(cè)緩沖區(qū)算法的多邊形骨架線提取模型［J］.海洋測繪，2010，30（5）：46－48.）

［13］ CHEN Jie.A Valuable Book of MATLAB［M］.Beijing：Publishing House of Electronics Industry，2010：139－170.（陳杰.MATLAB寶典（第二版）［M］.北京：電子工業(yè)出版社，2010：139－170.）

［14］ LIU Zhenjun.A Valuable Book of Matlab in Scientific Calculation and Visualization［M］.Beijing：Publishing House of Electronics Industry，2009：191－210.（劉正君.Matlab科學(xué)計算與可視化仿真寶典［M］.北京：電子工業(yè)出版社，2009：191－210.）

［15］ HE Guangyu.Frequently－used Numerical Arithmetic Aggregation of Visual C＋＋［M］.Beijing：Science Press，2002：99－100.（何光愈.Visual C＋＋常用數(shù)值算法集［M］.北京：科學(xué)出版社，2002：99－100.）

［16］ WANG Xin.Research on Identical Entity Geometric Matching in Multi－Source Spatial Data［D］.Zhengzhou：Information Engineering University，2008：10－11.（王馨.多源空間數(shù)據(jù)同名實體幾何匹配方法研究［D］.鄭州：信息工程大學(xué)，2008：10－11.）

［17］ XU Jianhua.Mathematical Methods in Contemporary Geography［M］.Beijing：Higher Education Press，2002：37－38.（徐建華.現(xiàn)代地理學(xué)中的數(shù)學(xué)方法［M］.北京：高等教育出版社，2002：37－38.）

［18］ AN Xiaoya，SUN Qun，XIAO Qiang，et al.A Shape Multilevel Description Method and Application in Measuring Geometry Similarity of Multi－scale Spatial Data［J］.Acta Geodaetica et Cartographica Sinica，2011，40（4）：495－501.（安曉亞，孫群，肖強(qiáng)，等.一種形狀多級描述方法及在多尺度空間數(shù)據(jù)幾何相似性度量中的應(yīng)用［J］.測繪學(xué)報，2011，40（4）：495－501.）

［19］ TANG Luliang，LI Qingquan，YANG Bisheng.Shape Similarity Measuring for Multi－resolution Transmission of Spatial Datasets over the Internet［J］.Acta Geodaetica et Cartographica Sinica，2009，38（4）：336－340.（唐爐亮，李清泉，楊必勝.空間數(shù)據(jù)網(wǎng)絡(luò)多分辨率傳輸?shù)膸缀螆D形相似性度量［J］.測繪學(xué)報，2009，38（4）：336－340.）

［20］ XU Junkui，WU Fang，WEI Huifeng.Research on Spatial Cognition Character of Areal Settlement Matching Process［J］.Journal of Geomatics Science and Technology，2012，29（4）：303－307.（許俊奎，武芳，魏慧峰.面狀居民地匹配的空間認(rèn)知特點(diǎn)研究［J］，測繪科學(xué)技術(shù)學(xué)報，2012，29（4）：303－307.）

［21］ ZHAI Renjian.Research on Automated Matching Methods for Multi－Scale Vector Spatial Data Based on Global Consistency Evaluation［D］.Zhengzhou：Information Engineering University，2011：82－87.（翟仁健.基于全局一致性評價的多尺度矢量空間數(shù)據(jù)匹配方法研究［D］.鄭州：信息工程大學(xué)，2011：82－87.）