簡(jiǎn)燦良，趙彬彬，鄧敏＊，余麗鈺

（1.武漢大學(xué)測(cè)繪遙感信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武昌 430079；2.中南大學(xué)地理信息系，湖南 長(zhǎng)沙 410083；3.福建省基礎(chǔ)地理信息中心，福建 福州 350003）

0 引言

隨著空間數(shù)據(jù)量的急速增長(zhǎng)，數(shù)據(jù)生產(chǎn)和應(yīng)用部門對(duì)已有數(shù)據(jù)重用和質(zhì)量改善的需求越來越迫切，亟須一種有效的技術(shù)手段用于空間數(shù)據(jù)集成和更新［1］。目前，一些發(fā)達(dá)國(guó)家的測(cè)繪部門已將工作重點(diǎn)從數(shù)據(jù)采集、建庫(kù)轉(zhuǎn)移到更新與應(yīng)用方面，并在建立更新機(jī)制、利用遙感影像獲取變化信息、歷史數(shù)據(jù)存取等方面做了大量工作［2］。例如，美國(guó)地質(zhì)調(diào)查局測(cè)繪部2001年在其國(guó)家地圖計(jì)劃（The National Map）中規(guī)劃到2010年建立近實(shí)時(shí)數(shù)據(jù)更新機(jī)制，將數(shù)據(jù)的現(xiàn)勢(shì)性保持在數(shù)月甚至幾天之內(nèi)［3］；英國(guó)軍械測(cè)量局（Ordnance Survey）根據(jù)實(shí)際情況建立了推帚式和散點(diǎn)式相結(jié)合的更新機(jī)制［4］；日本GSI基于柵格采用先更新1∶2.5萬地形圖、再用1∶2.5萬地形圖更新1∶5萬地形圖的更新方法。

在多尺度地圖之間傳播更新，首先需要通過對(duì)空間目標(biāo)幾何、拓?fù)浜驼Z(yǔ)義進(jìn)行相似性度量，識(shí)別出同一地區(qū)不同來源、不同尺度地圖數(shù)據(jù)庫(kù)中的同一地物，即目標(biāo)匹配，進(jìn)而通過比較、分析和識(shí)別判斷兩幅地圖中的同名目標(biāo)是否發(fā)生變化，即變化探測(cè)。變化探測(cè)是一個(gè)復(fù)雜而繁瑣的過程［5］，近年來遙感平臺(tái)和傳感器的快速發(fā)展使得國(guó)內(nèi)外學(xué)者更關(guān)注遙感影像數(shù)據(jù)變化探測(cè)方法及評(píng)價(jià)等相關(guān)問題［6－8］，主要集中在影像處理和影像數(shù)據(jù)時(shí)間維的操作方法上，對(duì)矢量地圖數(shù)據(jù)變化探測(cè)方面的關(guān)注不夠，研究成果較為零散［9－14］，特別是對(duì)不同比例尺矢量地圖更新過程中的變化探測(cè)研究更少。鑒于此，本文研究不同比例尺地圖空間面目標(biāo)的變化探測(cè)方法，主要是在目標(biāo)匹配基礎(chǔ)上，探討多尺度地圖空間匹配目標(biāo)對(duì)之間的差異問題及其產(chǎn)生的原因，進(jìn)而討論空間面目標(biāo)的變化類型、判別依據(jù)及變化探測(cè)等解決方案。

1 空間目標(biāo)差異驅(qū)動(dòng)因素和變化類型

1.1 空間目標(biāo)表達(dá)差異驅(qū)動(dòng)因素

導(dǎo)致不同比例尺地圖空間中的同名目標(biāo)間差異的原因（即差異驅(qū)動(dòng)因素）主要包括：制圖綜合、實(shí)際變化和數(shù)據(jù)獲取誤差（圖1）。數(shù)據(jù)獲取誤差導(dǎo)致的差異主要體現(xiàn)在不同來源的地圖數(shù)據(jù)中，如不同部門同一地區(qū)同一時(shí)期的相同比例尺地圖數(shù)據(jù)中，而同一地物的差異主要來源于數(shù)據(jù)獲取誤差。下面主要探討制圖綜合和實(shí)際變化驅(qū)動(dòng)的差異。

圖1 不同比例尺地圖空間同名目標(biāo)之間的差異驅(qū)動(dòng)因素Fig.1 The driven factors of discrepancy between corresponding objects from maps of different scales

1.1.1 制圖綜合驅(qū)動(dòng) 地圖符號(hào)是對(duì)地理要素的抽象表達(dá)，不同要素使用不同地圖符號(hào)表達(dá)，同時(shí)亦與表達(dá)尺度直接相關(guān)。在不同比例尺地圖中，同一地物表達(dá)的詳細(xì)程度不同［15，16］，并最終通過制圖綜合的方法實(shí)現(xiàn)。尤其對(duì)大比例尺地圖數(shù)據(jù)進(jìn)行制圖綜合派生出中、小比例尺地圖數(shù)據(jù)時(shí)，綜合操作算子［17］（如選取、合并、夸大等）的運(yùn)用，使得不同比例尺地圖對(duì)同一要素的表達(dá)存在區(qū)別，顯然，這種差異主要由制圖綜合驅(qū)動(dòng)。如圖2所示，圖2a和圖2c分別為較大比例尺面目標(biāo)，圖2b和圖2d分別為與之相應(yīng)的較小比例尺面目標(biāo)。圖2a和圖2b匹配目標(biāo)之間形狀差異明顯，即較大比例尺面目標(biāo)西南方向的缺角經(jīng)綜合后被“簡(jiǎn)化”；圖2c和圖2d匹配目標(biāo)間數(shù)量差異明顯，即3個(gè)較大比例尺面目標(biāo)經(jīng)綜合后被“聚合”為1個(gè)較小比例尺面目標(biāo)。因此，制圖綜合驅(qū)動(dòng)的差異源于綜合操作算子的作用，具體體現(xiàn)在匹配目標(biāo)對(duì)包含的目標(biāo)數(shù)量、幾何形狀和空間位置等方面，某些綜合算子（如Collapse）亦可導(dǎo)致匹配目標(biāo)對(duì)在目標(biāo)類型、維度方面的差異。

圖2 制圖綜合驅(qū)動(dòng)的差異Fig.2 Discrepancy driven by cartographic generalization

1.1.2 實(shí)際變化驅(qū)動(dòng) 如圖3所示，圖3a和圖3c分別為較大比例尺面目標(biāo)，圖3b和圖3d分別為與之匹配的較小比例尺面目標(biāo)。圖3b中3號(hào)面目標(biāo)僅與圖3a中1號(hào)面目標(biāo)匹配，兩者之間形狀、大小等幾何特征均一致，故該地物未發(fā)生變化。而2號(hào)面目標(biāo)則無與之匹配的面目標(biāo)，其大小近似于1號(hào)面目標(biāo)，可推斷不屬于制圖綜合差異，而應(yīng)為新出現(xiàn)的地物，從而在現(xiàn)勢(shì)性較弱的較小比例尺地圖（圖3b）中無對(duì)應(yīng)面目標(biāo)。4、5號(hào)面目標(biāo)亦是一對(duì)匹配目標(biāo)，兩者的形狀、大小差異明顯，且差異大小超出表達(dá)尺度及制圖綜合所能引起的范圍，分析可知該處地物已發(fā)生變化。

圖3 實(shí)際變化驅(qū)動(dòng)的差異Fig.3 Discrepancy driven by real change

由此可見，實(shí)際變化驅(qū)動(dòng)的差異由現(xiàn)實(shí)世界地物發(fā)生實(shí)際變化所致，主要體現(xiàn)在匹配目標(biāo)對(duì)中大小、形狀和數(shù)量等幾何特征上，某些情況亦可導(dǎo)致目標(biāo)對(duì)在類型、維度和屬性方面的差異。不難看出，制圖綜合和實(shí)際變化都可導(dǎo)致匹配目標(biāo)出現(xiàn)各種差異，這些差異體現(xiàn)在目標(biāo)幾何特征（如形狀、大小、維度等）、語(yǔ)義特性（如屬性信息）上。變化探測(cè)的主要任務(wù)就是區(qū)分實(shí)際變化驅(qū)動(dòng)的差異和制圖綜合驅(qū)動(dòng)的差異，進(jìn)而為判斷變化類型提供依據(jù)。

1.2 面／面目標(biāo)變化類型

參照現(xiàn)實(shí)世界中地物實(shí)際變化情況可知，地物變化類型主要包括：1）出現(xiàn)：原地表不存在某類地物，現(xiàn)在新出現(xiàn)了該類地物，如原本空曠的地表新挖了一個(gè)水塘；2）消失：原地表存在某類地物，現(xiàn)在該類地物不再存在，如水塘被填平，成為空地；3）改變：原地表存在某類地物，現(xiàn)在該類地物被另外的地物所替代，如水塘被填平并新建了房屋。對(duì)于變化類型“改變”，可將其理解為“消失”和“出現(xiàn)”兩種變化的合成，即原有的地物消失后，又出現(xiàn)新地物，且新地物在大小、形狀、屬性等方面至少有一項(xiàng)不同于原有地物。3種變化類型分別反映了現(xiàn)實(shí)世界中地物“從無到有”、“從有到無”和“從某地物（群）到其他地物（群）”的變化過程。

若將上述面狀地物視為居民地，則上述3種變化類型分別對(duì)應(yīng)新建（即一片空地上出現(xiàn)一棟新建房屋）、拆除（舊建筑被拆除后不再建新建筑）和改建（將舊房屋拆除后再新建房屋，新建房屋在形狀、大小等方面至少有一項(xiàng)有別于舊建筑）。

2 面目標(biāo)匹配模式與變化類型判別

2.1 面／面目標(biāo)匹配模式

由于地理要素在不同比例尺地圖中的表達(dá)形式各異，以及點(diǎn)、線和面3類基本空間目標(biāo)之間的排列組合，以致匹配目標(biāo)對(duì)的構(gòu)成包括6種可能情況［18］（從較大比例尺到較小比例尺：點(diǎn)與點(diǎn)、線與點(diǎn)、線與線、面與點(diǎn)、面與線和面與面）。本文以面與面匹配目標(biāo)對(duì)為例，探討其變化探測(cè)問題?？紤]到匹配目標(biāo)對(duì)各自包含的目標(biāo)數(shù)量差異，面面匹配目標(biāo)對(duì)又分為1∶0、1∶1、0∶1、1∶M、N∶1和 N∶M 共6種匹配模式［19］（圖4），這些匹配模式不僅反映了不同時(shí)段地面目標(biāo)數(shù)目的變化，而且表達(dá)了制圖綜合操作中不同比例尺地圖目標(biāo)之間的對(duì)應(yīng)關(guān)系，亦傳遞了不同類型的變化信息。

2.2 面／面匹配目標(biāo)對(duì)差異度量指標(biāo)

常用的面目標(biāo)度量指標(biāo)包括大小、形狀和質(zhì)心等幾何特征，顧及幾何參數(shù)表達(dá)的繁簡(jiǎn)程度及計(jì)算量差異，這里按參數(shù)計(jì)算量大小依次選取邊界長(zhǎng)度差、面積差和重疊度3個(gè)指標(biāo)對(duì)匹配目標(biāo)的差異進(jìn)行度量（考慮到涉及兩個(gè)以上面目標(biāo)時(shí)，邊界長(zhǎng)度和面積差異均較大，實(shí)際意義不大，故邊界長(zhǎng)度差和面積差僅用于1∶0、1∶1和0∶1模式），進(jìn)而對(duì)差異驅(qū)動(dòng)因素進(jìn)行“由簡(jiǎn)到難”、“由快速到準(zhǔn)確”的層次判斷。

（1）設(shè)匹配面目標(biāo)對(duì)中較大、較小比例尺地圖面目標(biāo)邊界周長(zhǎng)分別為PN、PO，邊界周長(zhǎng)差異記為δP，其閾值記為εP，則δP可表達(dá)為：

（2）設(shè)匹配面目標(biāo)對(duì)中較大、較小比例尺地圖面目標(biāo)面積分別為AN、AO，面積差異記為δA，其閾值記為εA，則δA可表達(dá)為：

（3）設(shè)匹配面目標(biāo)對(duì)中面目標(biāo)并集的面積為AU，面目標(biāo)交集的面積為AI，兩面目標(biāo)重疊度記為δO，重疊度閾值記為εO，則δO可表達(dá)為：

鑒于兩個(gè)面目標(biāo)匹配的多種模式［19］及不同匹配模式所用判斷指標(biāo)的差異，下面按匹配模式分別探討面／面匹配目標(biāo)變化類型判別方法。

2.3 面／面匹配目標(biāo)變化類型判別

2.3.1 1∶0模式 如圖4，現(xiàn)勢(shì)性強(qiáng)的較大比例尺地圖面目標(biāo)無與之對(duì)應(yīng)的現(xiàn)勢(shì)性弱的較小比例尺地圖面目標(biāo)，根據(jù)式（1）、式（2）和式（3）計(jì)算得到邊界周長(zhǎng)和面積差異均為1，重疊度為0，此模式的匹配面目標(biāo)之間的差異非常大，但該差異的可能驅(qū)動(dòng)因素包括制圖綜合和實(shí)際變化。假設(shè)該面目標(biāo)的面積為AN，較小比例尺地圖所能表達(dá)的最小面積為AMin，若AN≤AMin，則該面狀地物未發(fā)生變化，差異由制圖綜合驅(qū)動(dòng)（如Selective Elimination）；若AN＞AMin，則該面狀地物已發(fā)生變化，差異由實(shí)際變化驅(qū)動(dòng)，其變化類型為“出現(xiàn)”。

2.3.2 1∶1模式 1∶1模式的匹配面目標(biāo)對(duì)之間的差異主要體現(xiàn)在形狀和大小方面，該差異可能由制圖綜合（如Simplification、Magnification等算子會(huì)使得匹配目標(biāo)對(duì)的形狀、大小、邊界周長(zhǎng)等特征發(fā)生差異，而Displacement算子則會(huì)使位置等發(fā)生差異）或?qū)嶋H變化驅(qū)動(dòng)。因此，此模式的匹配面目標(biāo)對(duì)需選用3種指標(biāo)進(jìn)行差異度量，根據(jù)式（1）、式（2）和式（3），分別計(jì)算邊界周長(zhǎng)差異、面積差異和重疊度，依次考察匹配對(duì)的邊界周長(zhǎng)差異、面積差異及重疊度是否在閾值范圍內(nèi)：若δP＞εP，則認(rèn)為差異由實(shí)際變化驅(qū)動(dòng)，地物已發(fā)生變化，變化類型為“改變”；若δP≤εP，且δA＞εA，則認(rèn)為差異由實(shí)際變化驅(qū)動(dòng)，地物已發(fā)生變化，變化類型為“改變”；若δP≤εP，且δA≤εA，又δO＜εO，則認(rèn)為差異由實(shí)際變化驅(qū)動(dòng)，該地物已發(fā)生變化，變化類型為“改變”；若δP≤εP，且δA≤εA，又δO≥εO，則認(rèn)為差異由制圖綜合驅(qū)動(dòng)，該地物未發(fā)生變化，如圖5所示。

2.3.3 0∶1模式 0∶1模式匹配目標(biāo)的差異程度類似于1∶0模式，不同的是現(xiàn)勢(shì)性較弱的較小比例尺地圖面目標(biāo)無與之對(duì)應(yīng)的現(xiàn)勢(shì)性較強(qiáng)的較大比例尺地圖面目標(biāo)，根據(jù)式（1）、式（2）和式（3）計(jì)算得到邊界周長(zhǎng)和面積差異均為1，重疊度為0，分析可知，該差異由實(shí)際變化驅(qū)動(dòng)，即地物已發(fā)生變化，變化類型為“消失”。

圖5 1∶1模式面／面匹配目標(biāo)變化層次判斷Fig.5 Hierarchical determination of corresponding areas for 1∶1 mode

2.3.4 1∶M、N∶1和N∶M 模式 1∶M 和 N∶M 模式的匹配面目標(biāo)對(duì)在實(shí)際生產(chǎn)中很少出現(xiàn)，其差異可能由制圖綜合或?qū)嶋H變化驅(qū)動(dòng)。N∶1模式的匹配面目標(biāo)對(duì)出現(xiàn)較多，在6種匹配模式中比重最大，其差異的可能驅(qū)動(dòng)因素包括制圖綜合和實(shí)際變化?？紤]到涉及多個(gè)面目標(biāo)時(shí)，邊界周長(zhǎng)差異和面積差異的實(shí)際意義不大，故1∶M、N∶1和N∶M3種模式的匹配對(duì)用“重疊度”指標(biāo)考察：根據(jù)式（3）計(jì)算重疊度大小，若δO≥εO，則認(rèn)為差異由制圖綜合（如Simplification、Aggregation或 Selective Elimination等算子）驅(qū)動(dòng)，該地物未發(fā)生變化；否則認(rèn)為差異由實(shí)際變化驅(qū)動(dòng)，該地物已發(fā)生變化，變化類型為“改變”。

3 匹配面目標(biāo)對(duì)變化探測(cè)實(shí)驗(yàn)及結(jié)果分析

3.1 面／面匹配目標(biāo)對(duì)變化探測(cè)實(shí)驗(yàn)

以兩種不同比例尺的居民地?cái)?shù)據(jù)作為實(shí)驗(yàn)對(duì)象，如圖6a和圖6b分別為實(shí)驗(yàn)區(qū)1∶2 000（現(xiàn)勢(shì)性較好）和1∶10 000（現(xiàn)勢(shì)性較差）居民地?cái)?shù)據(jù)，面目標(biāo)數(shù)分別為732和181。兩種比例尺地圖中匹配目標(biāo)對(duì)共計(jì)550對(duì)，包括1∶0模式369對(duì)，1∶1模式98對(duì)，0∶1模式5對(duì)，N∶1模式78對(duì)，另兩種模式（即1∶M 和N∶M）均為0對(duì)。參照相關(guān)規(guī)范［20，21］規(guī)定的圖上最小尺寸（即?。瓷鲜龇椒ǚ謩e對(duì)各種模式（即分別?。┑钠ヅ淠繕?biāo)對(duì)進(jìn)行考察得到變化探測(cè)結(jié)果如表1。從表1中可以看出：1）1∶0模式匹配目標(biāo)對(duì)中共有93個(gè)1∶2 000比例尺面目標(biāo)（圖上面積大于等于1mm2）為新出現(xiàn)面目標(biāo)，即其變化類型為“出現(xiàn)”，另外276個(gè)面目標(biāo)（圖上面積小于1mm2）未發(fā)生變化，匹配目標(biāo)間的差異由制圖綜合驅(qū)動(dòng)；2）1∶1模式匹配目標(biāo)對(duì)中13個(gè)面目標(biāo)已發(fā)生變化，其變化類型為“改變”，另外85個(gè)面目標(biāo)未發(fā)生變化，匹配目標(biāo)間的差異由制圖綜合驅(qū)動(dòng)；3）0∶1模式匹配目標(biāo)對(duì)中的5個(gè)面目標(biāo)均已發(fā)生變化，其變化類型為“消失”；4）N∶1模式匹配目標(biāo)對(duì)中的8對(duì)面目標(biāo)已發(fā)生變化，變化類型為“改變”，其余70對(duì)面目標(biāo)未發(fā)生變化，其差異由制圖綜合驅(qū)動(dòng)。

3.2 實(shí)驗(yàn)結(jié)果分析

圖6 1∶2 000與1∶10 000居民地匹配目標(biāo)變化探測(cè)實(shí)驗(yàn)Fig.6 Change detection between corresponding areas from residential maps of 1∶2 000 and 1∶10 000 scales

表1 1∶2 000、1∶10 000比例尺地圖面目標(biāo)匹配對(duì)各模式變化類型與頻次的統(tǒng)計(jì)Table 1 Summary of detecting results between corresponding areas from residential maps of 1∶2 000 and 1∶10 000 scales

表2 1∶2 000、1∶10 000匹配面目標(biāo)差異度量樣本Table 2 Samples of discrepancy measuring between corresponding areas in maps of 1∶2 000 and 1∶10 000 scales

表2為實(shí)驗(yàn)結(jié)果的部分采樣，其中無填充面目標(biāo)和斜線填充面目標(biāo)比例尺分別為1∶2 000和1∶10 000，各模式匹配目標(biāo)對(duì)變化探測(cè)示例如下：1）1∶0模式：168號(hào)目標(biāo)面積 A168為1 012.735m2，大于AMin（圖上1mm2），變化類型為“出現(xiàn)”，即“新建”房屋，166號(hào)目標(biāo)面積A166為16.533m2，小于AMin，未變化；2）1∶1模 式：263號(hào)目標(biāo)對(duì)周長(zhǎng)差異δP為83.4%，大于εP（15%），變化類型為“改變”，即“改建”房屋，157號(hào)目標(biāo)對(duì)周長(zhǎng)差異δP為12.4%，小于εP（15%），面積δA差異為3.9%，小于εA（20%），重疊度δO為92.6%，大于εO（75%），未變化；3）0∶1模式：82號(hào)面目標(biāo)變化類型為“消失”，即“拆除”房屋；4）N∶1模式：179號(hào)目標(biāo)對(duì)重疊度δO為31.5%，小于εO（75%），變化類型為“改變”，即“改建”房屋，14號(hào)目標(biāo)對(duì)重疊度δO為87.5%，大于εO，未變化。

從面／面匹配目標(biāo)變化探測(cè)實(shí)驗(yàn)結(jié)果看，效果良好，達(dá)到了預(yù)期目標(biāo)，能準(zhǔn)確、有效地區(qū)分不同模式的匹配面目標(biāo)之間的實(shí)際變化差異和制圖綜合差異。總體而言，上述方法取得了較理想的變化探測(cè)結(jié)果，但亦存在有待深入研究之處，如實(shí)驗(yàn)結(jié)果與比例尺跨度和閾值大小的關(guān)系，同時(shí)，受到某種匹配模式的現(xiàn)實(shí)出現(xiàn)概率很小和可獲得的有限實(shí)驗(yàn)數(shù)據(jù)的局限，未能獲得某種匹配模式的目標(biāo)對(duì)的變化探測(cè)結(jié)果。

4 結(jié)語(yǔ)

目標(biāo)匹配是實(shí)現(xiàn)多尺度地圖數(shù)據(jù)合并、持續(xù)更新的關(guān)鍵技術(shù)和重要環(huán)節(jié)之一，其目的是利用已更新的現(xiàn)勢(shì)性強(qiáng)的較大比例尺地圖數(shù)據(jù)更新現(xiàn)勢(shì)性弱的較小比例尺地圖數(shù)據(jù)，從而在不同比例尺地圖數(shù)據(jù)間傳播變化，實(shí)現(xiàn)多尺度地圖協(xié)同更新［22］。目標(biāo)匹配技術(shù)的應(yīng)用已深入到地理信息科學(xué)領(lǐng)域的空間數(shù)據(jù)質(zhì)量改善和評(píng)價(jià)、多源空間數(shù)據(jù)集成（或融合）、多尺度空間數(shù)據(jù)庫(kù)的維護(hù)和更新、基于位置服務(wù)的導(dǎo)航等諸多方面［23，24］。本文從空間目標(biāo)的面積等幾何特征出發(fā)，對(duì)多尺度矢量地圖空間目標(biāo)變化探測(cè)問題進(jìn)行探討，分析了多尺度地圖空間匹配面目標(biāo)之間的差異產(chǎn)生原因，并將差異驅(qū)動(dòng)因素歸納為制圖綜合和實(shí)際變化兩類，明確了變化探測(cè)的目的和任務(wù)，對(duì)多尺度空間面目標(biāo)的變化類型進(jìn)行了歸納、分類，提出了針對(duì)面面匹配目標(biāo)的差異驅(qū)動(dòng)因素判別方法和變化探測(cè)方法，實(shí)驗(yàn)亦證明了該方法的可行性和實(shí)用性，為多尺度地圖更新提供了技術(shù)方法。

［1］ 陳軍，李志林，蔣捷，等.基礎(chǔ)地理數(shù)據(jù)庫(kù)的持續(xù)更新問題［J］.地理信息世界，2004，2（5）：1－5.

［2］ HEIPKE C.Requirements for modern GIS［A］.CHEN J，JIANG J.Proceedings of ISPRS Commission II Symposium［C］.2002.173－179.

［3］ USGS，2001，http：／／nationalmap.usgs.gov，2013－05－02.

［4］ MURRAY K.Anewgeo－infromation Framework for Great Britain［C］.FIG XXII International Congress，2002.

［5］ GONG J Y，SUI H G，MA G R，et al.A review of multi－temporal remote sensing data change detection algorithms［A］.The International Archives of the Photogrammetry，Remote Sensing and Sensing and Spatial Information Sciences［C］.Vol.XXXVII.Part B7.2008.757－762.

［6］ TIAN M，WAN S，YUE L.A novel approach for change detection in remote sensing image based on saliency map［A］.International Conference on Computer Graphics，Imaging and Visualization［C］.2007.397－402.

［7］ SONG Y，JIA X X.Remote sensing change detection based on growing hierarchical self－organization map［A］.Environmental Science and Information Application Technology（ESIAT）［C］.2010.

［8］ 杜培軍，柳思聰.融合多特征的遙感影像變化檢測(cè)［J］.遙感學(xué)報(bào)，2012，16（4）：663－677.

［9］ BADARD T.Towards a generic updating tool for geographic databases［A］.GIS／LIS′98［C］.1998.352－363.

［10］ ANDERS K H，BOBRICH J.MRDB Approach for Automatic Incremental Update［C］.ICA Workshop on Generalization and Multiple Representation，2004.

［11］ 張保鋼，袁燕巖.城市大比例尺地形圖數(shù)據(jù)庫(kù)中地物變化的自動(dòng)發(fā)現(xiàn)［J］.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2005，30（7）：640－642.

［12］ 陳軍，胡云崗，趙仁亮，等.道路數(shù)據(jù)縮編更新的自動(dòng)綜合方法研究［J］.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2007，32（11）：1022－1027.

［13］ 吳建華，傅仲良.數(shù)據(jù)更新中要素變化檢測(cè)與匹配方法［J］.計(jì)算機(jī)應(yīng)用，2008，28（6）：1612－1615.

［14］ ARMENAKIS C，CYR I.PAPANIKOLAOU E.Change detection methods for the revision of topographic databases［A］.Proceedings of the Joint International Symposium on Geospatial Theory，Processing and Applications［C］.2002.

［15］ 李志林.地理空間數(shù)據(jù)處理的尺度理論［J］.地理信息世界，2005，3（2）：1－5.

［16］ 艾廷華，成建國(guó).對(duì)空間數(shù)據(jù)多尺度表達(dá)有關(guān)問題的思考［J］.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2005，30（5）：377－382.

［17］ LI Z，YAN H，AI T，et al.Automated building generalization based on urban morphology and Gestalt theory［J］.International Journal of Geographical Information Science，2004，18（5）：513－534.

［18］ 徐楓，鄧敏，趙彬彬，等.空間目標(biāo)匹配方法的應(yīng)用分析［J］.地球信息科學(xué)學(xué)報(bào)，2009，11（5）：657－663.

［19］ 趙彬彬，鄧敏，徐震，等.多尺度地圖面目標(biāo)匹配的統(tǒng)一規(guī)則研究［J］.武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2011，36（8），991－994.

［20］ 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB／T 20257.2—2006.國(guó)家基本比例尺地形圖圖式第2部分：1∶5 000 1∶10 000地形圖圖式［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2006.

［21］ 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局、中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).GB／T 12343.1—2008.國(guó)家基本比例尺地圖編繪規(guī)范第1部分：1∶25 000 1∶50 000 1∶100 000地形圖編繪規(guī)范［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2008.

［22］ 趙彬彬.多尺度矢量地圖空間目標(biāo)匹配方法及其應(yīng)用研究［D］.中南大學(xué)，2011.

［23］ 田原，蔡華，鄔倫，等.基于局部相似的地籍宗地圖斑匹配方法［J］.地理與地理信息科學(xué)，2011，27（5）：25－28.

［24］ 吳建華.顧及環(huán)境相似的多特征組合實(shí)體匹配方法［J］.地理與地理信息科學(xué)，2010，26（4）：1－6.