賀曉暉，陳楠

(福州大學(xué) 福建省空間信息工程研究中心，福州 350002)

1 引 言

地圖制圖綜合歷來是地圖學(xué)中最富挑戰(zhàn)性和創(chuàng)造性的研究領(lǐng)域[1]，一直受到專業(yè)人士的重視，其中，等高線的綜合占據(jù)非常重要的地位。等高線是地面上高程相等的點(diǎn)所連成的閉合曲線在水平面上的投影，是地形圖上用來表達(dá)地形起伏形態(tài)的主要手段之一。真實(shí)地貌形態(tài)的多樣性決定了等高線的復(fù)雜性，隨著地形圖比例尺的縮小，等高線的小細(xì)節(jié)不斷增多，產(chǎn)生數(shù)據(jù)冗余，影響地圖的閱讀，而且給地圖數(shù)據(jù)的存儲(chǔ)和處理帶來了不便。為了節(jié)約數(shù)據(jù)的存儲(chǔ)空間，加快空間數(shù)據(jù)更新速度，較好的解決辦法是在GIS數(shù)據(jù)庫中存儲(chǔ)單一大比例尺等高線數(shù)據(jù)，然后根據(jù)用戶需求派生各種小比例尺等高線數(shù)據(jù)，這就涉及到等高線的綜合過程。等高線綜合就是為滿足目標(biāo)比例尺下的應(yīng)用要求，對(duì)等高線進(jìn)行一定程度的選取、化簡(jiǎn)和位移等操作的總稱。由于等高線具有復(fù)雜、密集、多樣等特點(diǎn)，成為地圖綜合的重點(diǎn)和難點(diǎn)，因此也在地圖學(xué)領(lǐng)域引起了眾多學(xué)者的廣泛關(guān)注。

等高線是一類攜帶地理信息的特殊的曲線要素，與一般的幾何曲線要素既有相似之處又有較大區(qū)別。目前，很多等高線綜合算法都是依托曲線壓縮算法實(shí)現(xiàn)的，考慮到等高線之間特殊的空間拓?fù)潢P(guān)系，眾多學(xué)者還從整體角度出發(fā)研究成組等高線的自動(dòng)綜合方法。

2 單根等高線綜合方法

單根等高線綜合方法把等高線看作是由有序的離散點(diǎn)組成的一般的曲線要素，根據(jù)目標(biāo)比例尺保留滿足條件的曲線特征點(diǎn)，刪除冗余的曲線點(diǎn)，使得綜合后的等高線既能最大程度的表達(dá)目標(biāo)比例尺下的曲線形態(tài)特征又能最大限度的刪除冗余節(jié)點(diǎn)。根據(jù)算法的思想差異，單根等高線綜合算法大致可分為抽樣取點(diǎn)法、局部處理算法以及整體處理算法等，抽樣取點(diǎn)法由于無法保證綜合前后等高線的相似性而用得比較少。

2.1 局部處理算法

局部處理算法的思想是：與目標(biāo)點(diǎn)相鄰的兩點(diǎn)或多點(diǎn)共同參與目標(biāo)點(diǎn)取舍的計(jì)算，通過計(jì)算目標(biāo)點(diǎn)與相鄰點(diǎn)之間的距離、角度等數(shù)值來衡量目標(biāo)點(diǎn)對(duì)曲線局部區(qū)域的重要程度，從而決定目標(biāo)點(diǎn)的取舍。屬于局部處理算法的種類繁多，如垂距法、角度限差法、斜率比較法、光柵法、最小面積重復(fù)刪除法、Tobler法、McMaster-Jenks法等。此外，還有基于約束擴(kuò)展的局部處理法(該算法應(yīng)用角度、距離、點(diǎn)數(shù)來搜索相鄰點(diǎn)，如Opheim法、Lang法、Johannsen法、Roberge法、Deveau法)以及基于非約束擴(kuò)展的局部處理法(該算法利用線的復(fù)雜度、坐標(biāo)點(diǎn)的密度、開始點(diǎn)的定位來搜索相鄰點(diǎn)，如Reumann-Witka[2-6]法等)。

這類局部處理算法運(yùn)算簡(jiǎn)單，處理速度快。但也存在自身難以克服的一些缺陷，忽略了目標(biāo)點(diǎn)對(duì)于保持整條曲線形態(tài)特征的重要性，有時(shí)可能會(huì)將曲線的彎曲極值點(diǎn)刪除，造成圖形失真；當(dāng)曲線比較平緩時(shí)，各相鄰節(jié)點(diǎn)之間的數(shù)學(xué)特征差異比較小，這類算法容易造成綜合后曲線特征點(diǎn)的聚集現(xiàn)象。

2.2 整體處理算法

整體處理算法的思想是：從全局角度出發(fā)考慮目標(biāo)點(diǎn)對(duì)保持曲線形態(tài)的重要程度，以此來判斷目標(biāo)點(diǎn)是否保留。最經(jīng)典的算法就是Douglas-Peucker算法(簡(jiǎn)稱DP算法)[7-10]。1973年前后多人同時(shí)提出DP算法，該算法是對(duì)垂距法的改進(jìn)研究，運(yùn)用遞歸的思想實(shí)現(xiàn)曲線的綜合，是一種常用的矢量數(shù)據(jù)壓縮算法[11]，根據(jù)算法的實(shí)現(xiàn)過程可分為錨點(diǎn)前進(jìn)法和分而治之法。分而治之法是為了提高計(jì)算效率而對(duì)錨點(diǎn)前進(jìn)法的改進(jìn)。

Douglas-Peucker算法的“分而治之法”首先將曲線的首末點(diǎn)連成一條直線，計(jì)算中間所有點(diǎn)到該直線的距離，找出最大的距離dmax，與事先設(shè)定的限差閾值D進(jìn)行比較。若d

圖1 Douglas-Peucker算法的“分而治之法”示意圖

用DP算法化簡(jiǎn)后的曲線具有嚴(yán)格的保凹凸性，能夠保持曲線的整體形態(tài)特征，克服了局部處理算法存在的一些缺陷。但該算法也存在缺陷，在用于地圖比較復(fù)雜的等高線或者比例尺變化范圍較大的情況下，運(yùn)算效率不高，沒有考慮等高線之間的拓?fù)潢P(guān)系，容易造成化簡(jiǎn)后曲線之間的相交以及自身的相交；閾值的設(shè)定始終靠人為經(jīng)驗(yàn)或反復(fù)實(shí)驗(yàn)獲得，計(jì)算量大又費(fèi)時(shí)，還沒有實(shí)現(xiàn)自動(dòng)化；沒有考慮化簡(jiǎn)起始位置的影響，不同的起始點(diǎn)會(huì)得到不同的化簡(jiǎn)結(jié)果；曲線壓縮精度一般用位移矢量和面積偏差來衡量，該算法用垂向距離作為選取點(diǎn)的約束條件可以控制位移矢量，但無法控制面積偏差的大小。暫且不提DP算法存在的這些缺陷，其對(duì)曲線的綜合效果還是令人滿意的，保留了曲線的重要特征點(diǎn)，比較適合矢量曲線的壓縮，所以眾多學(xué)者還是熱衷于對(duì)該算法的改進(jìn)研究，在如何提高該算法的綜合效率[13-17]、如何解決該算法在曲線綜合中出現(xiàn)漏選特征點(diǎn)的問題[18-20]、如何改進(jìn)該算法只采用垂向距離作為選點(diǎn)的限制條件的缺陷[10，21]、如何改進(jìn)該算法以避免綜合后引起等高線的相交[22-24]等方面都有研究，針對(duì)其缺陷或者進(jìn)行算法改進(jìn)或者進(jìn)行一些后期處理，使得對(duì)曲線的綜合過程更加完善。費(fèi)立凡等[12，25-27]在分析2維Douglas-Peucker算法原理的基礎(chǔ)上，提出了3維Douglas-Peucker算法，之后又對(duì)其進(jìn)行改進(jìn)研究，引入“孤獨(dú)指數(shù)”的概念，并將該算法用于規(guī)則格網(wǎng)DEM和等高線的自動(dòng)綜合，初步實(shí)驗(yàn)表明，該算法保持了主要的地貌形態(tài)特征，綜合效果比較好。

2.3 基于自然規(guī)律的Li-Openshaw算法

Zhilin Li和Stan Openshaw將人類由近及遠(yuǎn)觀察物體的視覺原理用于研究曲線綜合，于1992年提出了Li-Openshaw算法[28]。該算法模擬了人類觀察物體的自然綜合過程，隨著人類遠(yuǎn)離空間目標(biāo)，某些細(xì)節(jié)特征將在視覺上無法表達(dá)出來甚至消失，此時(shí)就認(rèn)為被綜合掉了。實(shí)際上，人類的肉眼視覺分辨率有一定的限制，存在一種叫做模糊圓的最小可分辨單元，即人類無法感知小于該最小可分辨單元的空間細(xì)節(jié)，經(jīng)研究確定目標(biāo)比例尺最小可視目標(biāo)的限定范圍為0.6mm～0.7mm[29]。該算法的簡(jiǎn)化過程是在目標(biāo)比例尺下放置一個(gè)大小為0.6mm～0.7mm的柵格，然后每個(gè)柵格單元內(nèi)以某種規(guī)則選取一個(gè)新的節(jié)點(diǎn)代替柵格單元內(nèi)所有節(jié)點(diǎn)，從而實(shí)現(xiàn)曲線的綜合，通常情況下該節(jié)點(diǎn)不是原始曲線上的點(diǎn)，如圖2所示：

Li和Openshaw對(duì)復(fù)雜海岸線及河流在不同比例尺上的人工綜合結(jié)果與該算法自動(dòng)生成的結(jié)果進(jìn)行了比較，證明該算法對(duì)這類線狀地物的綜合非常有效，且運(yùn)行效率高，能有效地減少綜合后曲線的自相交。但用于等高線綜合時(shí)存在難以克服的缺陷：綜合過程采用統(tǒng)一大小的柵格單元，柵格單元大小一經(jīng)確定，每條曲線經(jīng)過的柵格內(nèi)都要記錄一個(gè)點(diǎn)，對(duì)于簡(jiǎn)單的曲線，會(huì)產(chǎn)生大量冗余點(diǎn)；對(duì)于較為復(fù)雜的曲線，有可能出現(xiàn)要素自身粘連的現(xiàn)象；最重要的是，用該算法綜合后的曲線點(diǎn)集大部分不是原始坐標(biāo)點(diǎn)，綜合質(zhì)量得不到保證。

2.4 基于分形分析的綜合法

1973年，B.B.Mandelbrot首次提出了分維和分形幾何學(xué)的設(shè)想[30]。分形的原意是不規(guī)則的、支離破碎的，指部分與整體以某種方式相似的形體。分維是描述分形體的數(shù)學(xué)手段，可以反映出分形的不規(guī)則性和復(fù)雜性，越復(fù)雜的形體，分維數(shù)越大[31]。地圖上表示的是自然界中復(fù)雜的地理現(xiàn)象，往往呈現(xiàn)出不規(guī)則的變化規(guī)律，因此，分形幾何應(yīng)用于制圖綜合具有得天獨(dú)厚的優(yōu)勢(shì)，為該領(lǐng)域的發(fā)展提供了一種新思想。分形分析就是一種能夠描述圖形形狀結(jié)構(gòu)及其復(fù)雜度的方法，并且建立圖形形狀結(jié)構(gòu)與尺度之間的數(shù)量關(guān)系，而尺度與比例尺的對(duì)應(yīng)關(guān)系使我們可以建立圖形形狀結(jié)構(gòu)及其復(fù)雜度與尺度之間的變化關(guān)系，分形理論的尺度依賴性為研究地圖目標(biāo)隨比例尺的變化規(guī)律提供了一種新的思路。

分形幾何用于地圖綜合的研究可以追溯到Dutton在1980年提出的用分形方法控制地圖細(xì)節(jié)的思想，從此，越來越多的學(xué)者在該領(lǐng)域展開了研究[32-34]。理想情況下，分維是尺度變化下的不變量，因此可以根據(jù)綜合前后圖形分維值保持不變，以目標(biāo)比例尺、現(xiàn)在的比例尺和分維值動(dòng)態(tài)地求出與目標(biāo)比例尺對(duì)應(yīng)的圖形簡(jiǎn)化閾值，然后應(yīng)用DP算法化簡(jiǎn)圖形[35]。但地理現(xiàn)象比較復(fù)雜，變化多樣，地圖目標(biāo)的相似性都帶有統(tǒng)計(jì)意義，在尺度變化下其分維值不再是個(gè)定值。在地圖綜合中，隨著比例尺的縮小難免會(huì)丟失一部分信息，復(fù)雜的目標(biāo)會(huì)被化簡(jiǎn)，不重要的細(xì)節(jié)會(huì)被刪除，信息量會(huì)有所損失，即地圖結(jié)構(gòu)形態(tài)隨著目標(biāo)比例尺的變化發(fā)生改變，造成用于描述其信息量的衰減，反應(yīng)到分維數(shù)的變化上應(yīng)該也是衰減的，因此，在應(yīng)用單一分維進(jìn)行地圖目標(biāo)的化簡(jiǎn)時(shí)存在一定的局限性。隨后，許多學(xué)者從變維分析的角度出發(fā)，研究地圖目標(biāo)的自相似性隨比例尺的變化而衰減的規(guī)律[36]。

分形方法的引入使得人們可以充分利用分維值這個(gè)量化指標(biāo)研究曲線的綜合，在一定程度上克服了現(xiàn)有綜合算法的缺陷，但也存在一定的局限性。分形理論應(yīng)用于曲線綜合，仍然需要與DP算法相結(jié)合，綜合結(jié)果只是比單純使用DP算法多選了一些點(diǎn)。分形理論利用分維值去衡量曲線的復(fù)雜程度，但客觀世界的地理現(xiàn)象只是統(tǒng)計(jì)意義下的分形，不同的分形體可能會(huì)有相同的分?jǐn)?shù)維，而且分維值也是統(tǒng)計(jì)得到的；另外，由于使用的概念和所采用的標(biāo)度與尺碼不同，同一圖形的分?jǐn)?shù)維也會(huì)有不同的數(shù)值。因此，曲線綜合得到的結(jié)果不能保證保留了原始曲線的主要特征。

2.5 其他算法

小波分析是20世紀(jì)80年代后期發(fā)展起來的一門新興的應(yīng)用數(shù)學(xué)理論，由于其在時(shí)域、頻域都具有表征信號(hào)局部特征能力和多分辨率分析的“數(shù)學(xué)顯微鏡”特點(diǎn)，成為眾多學(xué)科領(lǐng)域解決自身難題的有力工具[37]。小波分析理論可提供在不同分辨率下分析表達(dá)信息的有效途徑，因此受到了地學(xué)領(lǐng)域?qū)W者的高度重視，其多分辨率分析原理與制圖綜合具有天然的聯(lián)系，進(jìn)入90年代，小波理論逐漸成熟并被廣泛應(yīng)用于數(shù)據(jù)壓縮領(lǐng)域，許多學(xué)者都將其應(yīng)用到數(shù)據(jù)的多比例尺分析和自動(dòng)制圖綜合研究中，于浩等[38-39]將小波多尺度分析應(yīng)用到DEM數(shù)據(jù)的壓縮簡(jiǎn)化，曲線要素是地圖綜合的重點(diǎn)，小波分析在線狀要素的簡(jiǎn)化方面也取得了很多重要的成果[40-44]。

此外，隨著等高線綜合算法的不斷改進(jìn)研究，最大通視原理也被用于研究等高線的化簡(jiǎn)，該算法利用最大通視性原理對(duì)等高線進(jìn)行區(qū)域劃分，在每個(gè)區(qū)域范圍內(nèi)對(duì)等高線進(jìn)行彎曲嵌套層次分析和等級(jí)劃分，然后根據(jù)曲線的等級(jí)和預(yù)先設(shè)定的閾值進(jìn)行等高線彎曲取舍[45]。

綜上所述，新的曲線綜合理論思想和算法不斷出現(xiàn)，算法在執(zhí)行效率和綜合質(zhì)量方面也在不斷改進(jìn)，各種算法的綜合原則不盡相同，處理方式也存在差異，但無論如何這些算法都是研究單根曲線要素的簡(jiǎn)化綜合，且每種地理要素都有其自身的特點(diǎn)，單純依靠統(tǒng)一的方法對(duì)所有線要素綜合是不現(xiàn)實(shí)的。作為地形圖上最復(fù)雜的曲線要素，應(yīng)該深入分析等高線之間的空間拓?fù)潢P(guān)系，根據(jù)等高線的特點(diǎn)制定特殊的綜合算法，才能使等高線的綜合過程更加合理。等高線的綜合不僅是單根曲線的綜合，還要考慮與相鄰等高線以及其他空間目標(biāo)之間的關(guān)系，因此單根曲線綜合法不適于直接應(yīng)用于等高線的簡(jiǎn)化。

3 成組等高線綜合方法

實(shí)際上，等高線的綜合比單根曲線的綜合要復(fù)雜的多。地形圖上，等高線是以成組的方式來表達(dá)地貌形態(tài)的，其簡(jiǎn)化不僅涉及到彎曲的判斷，還涉及到所綜合區(qū)域地形特征的判斷，以及與其他地理要素的空間關(guān)系。因此，將等高線以成組的方式進(jìn)行綜合更符合人們的實(shí)際思維方式，目前成組等高線的綜合主要有兩種方法。

3.1 等高線結(jié)構(gòu)化綜合法

地性線是描述地貌形態(tài)時(shí)的控制線，反映了等高線之間的內(nèi)在聯(lián)系，因此，可以用地性線來控制等高線的綜合，這與傳統(tǒng)制圖中手工綜合等高線的原理是一致的。等高線結(jié)構(gòu)化綜合法的主要思路是：首先根據(jù)原始數(shù)據(jù)提取地性線，建立地性線的層次樹；其次生成地性線的Voronoi圖作為同一層次中地性線差異的判斷指標(biāo)，通過樹結(jié)構(gòu)和Voronoi圖對(duì)地性線進(jìn)行綜合；最后以綜合后的地性線指導(dǎo)對(duì)應(yīng)等高線彎曲的取舍。

等高線結(jié)構(gòu)化綜合方法的重點(diǎn)和難點(diǎn)在于地性線的提取和評(píng)價(jià)。地性線是地貌形態(tài)的骨架線，控制著地貌起伏變化的形態(tài)，主要包括山谷線和山脊線，山谷線和山脊線是正交于等高線的曲線族，在地形圖表示為一組相關(guān)的等高線彎曲，由于等高線彎曲的復(fù)雜嵌套，基于DEM提出來的地形特征線也處于一種無組織、無層次、雜亂無章的狀態(tài)。等高線結(jié)構(gòu)化綜合的實(shí)質(zhì)是保留重要地形特征對(duì)應(yīng)的等高線彎曲組，舍去次要地形特征對(duì)應(yīng)的等高線彎曲組，故在等高線結(jié)構(gòu)化綜合之前，評(píng)價(jià)每條地形特征線所代表地形的重要程度顯得至關(guān)重要。目前對(duì)山谷進(jìn)行評(píng)價(jià)的指標(biāo)主要有長(zhǎng)度、谷間距、匯流累積量、等級(jí)等，事實(shí)上，這些指標(biāo)之間是相關(guān)的，不同的指標(biāo)適應(yīng)不同的數(shù)據(jù)類型。由于地形結(jié)構(gòu)的復(fù)雜性，有時(shí)候用單一指標(biāo)往往難以準(zhǔn)確的評(píng)價(jià)山谷線，因此可將多個(gè)指標(biāo)綜合運(yùn)用評(píng)價(jià)山谷的質(zhì)量，并建立山谷線的結(jié)構(gòu)樹。

眾多學(xué)者從不同角度對(duì)基于地性線的結(jié)構(gòu)化綜合方法進(jìn)行了研究，著重于研究如何基于等高線提取地性線，如何以地性線控制等高線的綜合[46-50]，或者研究等高線之間的拓?fù)潢P(guān)系，如何根據(jù)等高線圖形之間的空間關(guān)系來建立某種樹結(jié)構(gòu)，以表達(dá)等高線之間的拓?fù)淇臻g關(guān)系[51-53]。另外，郭慶勝[54-55]提出了漸進(jìn)式曲線綜合的思想，采用從局部到整體的方式，以圖形的漸變實(shí)現(xiàn)圖形的突變，隨著比例尺的縮小，圖形的細(xì)節(jié)不斷刪除，從而實(shí)現(xiàn)目標(biāo)比例尺線圖形的綜合。

等高線結(jié)構(gòu)化綜合方法顧及到了等高線之間的空間關(guān)系，用地性線控制等高線的綜合過程，在一定程度上避免了綜合后等高線的相交現(xiàn)象。但如何準(zhǔn)確地提取及評(píng)價(jià)地性線是該算法的難點(diǎn)，而且綜合結(jié)果還不能令人滿意，沒有達(dá)到傳統(tǒng)手工綜合的質(zhì)量要求，并且存在過程復(fù)雜，效率不高等諸多問題。

此外，有學(xué)者研究了基于地貌高程帶的等高線綜合方法，該方法在相鄰兩條異高等高線所定義的高程帶內(nèi)插生成中間等高線以實(shí)現(xiàn)等高線的綜合[56-57]。

3.2 等高線的間接綜合法

等高線間接綜合法就是利用等高線與格網(wǎng)DEM可以相互派生的關(guān)系將輸入的等高線數(shù)據(jù)首先轉(zhuǎn)換為其他格式的一定分辨率的格網(wǎng)DEM數(shù)據(jù)，根據(jù)目標(biāo)比例尺的精度要求首先對(duì)格網(wǎng)DEM數(shù)據(jù)進(jìn)行綜合，然后利用簡(jiǎn)化后的格網(wǎng)DEM數(shù)據(jù)派生對(duì)應(yīng)比例尺的等高線數(shù)據(jù)，從而實(shí)現(xiàn)等高線的綜合。

格網(wǎng)DEM、不規(guī)則三角網(wǎng)DEM、等高線是數(shù)字高程模型最常見的3種表現(xiàn)形式。GRID模型具有數(shù)據(jù)規(guī)則排列的特點(diǎn)，對(duì)其簡(jiǎn)化的核心思想是將GRID模型看作柵格影像圖，采用濾波算子對(duì)每個(gè)柵格的灰度值進(jìn)行平滑，刪除地表貢獻(xiàn)較小的地形點(diǎn)，從而實(shí)現(xiàn)GRID模型的簡(jiǎn)化[58-59]。TIN簡(jiǎn)化算法的基本思想是從初始最細(xì)化水平層次開始，逐層對(duì)三角形頂點(diǎn)進(jìn)行消去或合并，重新構(gòu)造TIN，從而達(dá)到TIN的層次簡(jiǎn)化。TIN模型簡(jiǎn)化主要有3種算法：①Hoppe提出的漸進(jìn)網(wǎng)格的思想對(duì)TIN模型進(jìn)行簡(jiǎn)化[2，60]；②最大獨(dú)立點(diǎn)法對(duì)TIN模型進(jìn)行簡(jiǎn)化[61]；③基于頂點(diǎn)樹的TIN模型層次動(dòng)態(tài)簡(jiǎn)化方法，該算法是Luebke提出的基于緊致四叉樹空間劃分的三角形頂點(diǎn)合并與分裂方法，從而建立動(dòng)態(tài)的層次地形表現(xiàn)模型。此外，于江佩[62]運(yùn)用道格拉斯算法得到更抽象的TIN實(shí)現(xiàn)等高線的化簡(jiǎn)。

等高線間接綜合得到的中間數(shù)據(jù)可以根據(jù)需要用于地理信息的空間分析，綜合后的等高線不會(huì)出現(xiàn)相交問題，綜合過程簡(jiǎn)單。但是從綜合目的來看，該方法要經(jīng)過二維數(shù)據(jù)到三維數(shù)據(jù)，三維數(shù)據(jù)到二維數(shù)據(jù)的轉(zhuǎn)換，轉(zhuǎn)換過程中必然會(huì)造成信息量的額外損失。

綜上所述，成組等高線的綜合算法一直處于研究探討階段，在一定程度上克服了單根等高線綜合算法的某些缺陷，但等高線之間的拓?fù)淇臻g關(guān)系問題仍然未得到有效解決。在實(shí)際工作中，如果對(duì)綜合結(jié)果的精度要求比較高，適合選擇等高線結(jié)構(gòu)化綜合方法，等高線間接綜合法則比較適合精度要求不高的研究。用地性線控制等高線的綜合更符合人們的思維方式，但該方法還有很多問題需要解決，以后的研究重點(diǎn)應(yīng)該放在如何準(zhǔn)確提取地性線和評(píng)價(jià)地性線上。為了更好地解決等高線的綜合問題，在等高線綜合之前分析研究區(qū)域的地形特點(diǎn)以及合理建立等高線之間的空間關(guān)系顯得至關(guān)重要，還應(yīng)該進(jìn)一步考慮地貌等其他要素對(duì)等高線綜合的影響以及等高線與其他要素之間的空間關(guān)系。

4 等高線綜合方法存在問題分析

等高線本身的復(fù)雜性給實(shí)現(xiàn)自動(dòng)綜合帶來了很多困難，多年來人們提出了很多自動(dòng)綜合算法，并且對(duì)算法本身存在的缺陷也在不斷改進(jìn)。單根等高線綜合算法沒有顧及到等高線之間的空間關(guān)系。成組等高線綜合方法把等高線看成一群線狀目標(biāo)，已經(jīng)初步考慮到了等高線之間的空間關(guān)系，但忽略了地貌的其他要素對(duì)等高線綜合產(chǎn)生的影響。為此，龍毅等[63]根據(jù)地形與河流之間的空間依賴性提出一種等高線簇與河網(wǎng)協(xié)同的制圖綜合方法。在以后的工作中，等高線綜合算法還存在很多問題需要進(jìn)一步研究解決。

(1)等高線綜合的原則。對(duì)于以正地形為主的地貌形態(tài)，擴(kuò)脊消谷；對(duì)于以負(fù)地形為主的地貌形態(tài)，擴(kuò)谷消脊。但目前大多數(shù)研究都是通過刪除山谷線分支達(dá)到等高線綜合的目的，并沒有考慮山脊線，且山脊線分級(jí)的研究較少，無法判斷山脊線對(duì)應(yīng)等高線彎曲的重要程度，嚴(yán)格意義上來說，這樣做不太科學(xué)。對(duì)于以負(fù)地形為主的地貌形態(tài)，能否通過刪除山脊線分支達(dá)到等高線綜合的目的還尚未見研究；山脊線和山谷線共同控制等高線簡(jiǎn)化過程能否提高綜合結(jié)果的精度也缺乏研究。

(2)目前很多學(xué)者都是以某一特定區(qū)域?yàn)檠芯繉?duì)象探討等高線綜合方法，但是地形類型不同，對(duì)應(yīng)等高線疏密、等高距、起伏度等都存在差異，這些方法是否能適用于所有地形的等高線綜合涉及到綜合算法的普遍適用性問題，是以后要研究的一個(gè)重要問題。

(3)等高線綜合就是對(duì)等高線的化簡(jiǎn)、取舍、移位、夸大等操作的總稱，但大多數(shù)研究都集中于等高線圖形的化簡(jiǎn)，而如何對(duì)等高線進(jìn)行移位以及如何對(duì)某些部分夸大卻鮮有研究。等高線綜合不僅僅是圖形的化簡(jiǎn)，在綜合時(shí)如何處理與其他地形要素之間的空間關(guān)系也非常重要，目前這類問題的研究也比較少。

(4)大多數(shù)研究者都以1∶1萬等高線綜合為1∶5萬等高線為例研究自動(dòng)綜合的算法，但是這些算法是否適用于1∶1萬綜合至1∶10萬、1∶20萬、1∶50萬、1∶100萬，是否適用于1∶5萬綜合至1∶10萬、1∶20萬、1∶50萬等都缺少研究，同時(shí)也涉及到這些方法能否廣泛使用以及能否圓滿地解決等高線的綜合問題。

(5)學(xué)者們已經(jīng)提出了很多等高線綜合方法，但是綜合精度是否滿足要求以及如何對(duì)綜合結(jié)果進(jìn)行評(píng)價(jià)尚缺乏研究。

上述這些問題都是以后等高線綜合工作的中要重點(diǎn)探討的問題。

參考文獻(xiàn)：

[1] 王家耀.空間數(shù)據(jù)自動(dòng)綜合研究進(jìn)展及趨勢(shì)分析[J].測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2008，25(1)：1-7.

[2] 王文迭，賈建紅.等高線的自動(dòng)綜合[J].北京測(cè)繪，2010，(3)：22-30.

[3] 劉可晶.一種改進(jìn)的矢量曲線數(shù)據(jù)壓縮算法[J].甘肅科學(xué)學(xué)報(bào)，2005，17(3)：112-115.

[4] 胡鵬，黃杏元，華一新.地理信息系統(tǒng)教程[M].武漢：武漢大學(xué)出版社，2002.

[5] MAHESWARI V，WHYATT J D.Line generalisation by repeated elimination of points[J].Cartographic Journal，1993，30(1)：46-51.

[6] REUMANN K，WITKAM A P M.Optimizing curve segmentation in computer graphics[C].Proceedings of the International Computing Symposium，1974，467-472.

[7] URS R.An iterative procedure for the polygonal approximation of plane curves[J].Computer Graphics and Image Processing，1972，1(3)：244-256.

[8] DAVID H D，THOMAS K P.Algorithms for the reduction of the number of points required to represent a digitized line or its caricature[J].Cartographica：The International Journal for Geographic Information and Geovisualization，1973，10(2)：112-122.

[9] ALAN S.Topologically consistent line simplification with the Douglas-Peucker algorithm[J].Cartography and Geographic Information Science，1999，26(1)：7-18.

[10] 楊得志，王杰臣，閭國(guó)年.矢量數(shù)據(jù)壓縮的Douglas-Peucker算法的實(shí)現(xiàn)與改進(jìn)[J].測(cè)繪通報(bào)，2002，(7)：18-22.

[11] 楊建宇，楊崇俊，明冬萍，等.WebGIS系統(tǒng)中矢量數(shù)據(jù)的壓縮與化簡(jiǎn)方法綜述[J].計(jì)算機(jī)工程與應(yīng)用，2004，40(32)：36-38.

[12] 費(fèi)立凡，何津，馬晨燕，等.三維Douglas-Peucker算法及其在DEM自動(dòng)綜合中的應(yīng)用研究[J].測(cè)繪學(xué)報(bào)，2006，35(3)：278-284.

[13] NIRWA A，DELP E J.On detecting dominant points[J].Pattern Recognition，1991，24(5)：441-451.

[14] 黃培之.具有預(yù)測(cè)功能的曲線矢量數(shù)據(jù)壓縮方法[J].測(cè)繪學(xué)報(bào)，1995，24(4)：316-320.

[15] 王晏民.矢量曲線的特征點(diǎn)提取[J].測(cè)繪工程，2002，11(2)：8-10.

[16] 夏仁波，劉偉軍，王越超.一種改進(jìn)的基于DP原理的分段輪廓采樣算法[J].計(jì)算機(jī)工程與應(yīng)用，2004，40(21)：21-24.

[17] 劉穎，翟京生.海岸線形態(tài)的表達(dá)與自動(dòng)綜合[J].測(cè)繪與空間地理信息，2005，28(4)：78-81.

[18] TADEUSZ C.A numerical method for generalizing the linear elements of large-scale maps，based on the example of rivers[J].Cartographica：The International Journal for Geographic Information and Geovisualization，2000，37(1)：49-56.

[19] 劉曉紅，李樹軍.矢量數(shù)據(jù)壓縮的角度分段道格拉斯算法研究[J].四川測(cè)繪，2005，28(2)：51-52.

[20] 張洪升，王金娜.矢量曲線自動(dòng)綜合算法的研究[J].測(cè)繪與空間地理信息，2011，34(3)：133-134.

[21] 王杰臣，閭國(guó)年.曲線矢量數(shù)據(jù)綜合的微凹凸消除算法[J].測(cè)繪通報(bào)，2002，(1)：17-19.

[22] 張勝，朱才連，鐘世明.Douglas-Peucker算法的改進(jìn)及應(yīng)用[J].武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版)，2005，29(5)：671-674.

[23] 楊平，胡鵬，邵春麗.一種拓?fù)浔３值恼劬€簡(jiǎn)化算法研究[J].地理空間信息，2006，4(6)：41-43.

[24] 趙永清，謝傳節(jié)，喬玉良，等.基于最值點(diǎn)的道格拉斯—普克壓縮算法[J].教育技術(shù)導(dǎo)刊，2008，7(11)：60-62.

[25] 何津，費(fèi)立凡.再論三維Douglas-Peucker算法及其在DEM綜合中的應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2008，33(2)：160-163.

[26] 黃麗娜，費(fèi)立凡.采用三維Douglas-Peucker算法的等高線綜合[J].測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2009，26(6)：444-448.

[27] 黃麗娜，費(fèi)立凡.采用3D D-P算法的等高線三維綜合實(shí)驗(yàn)研究[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2010，35(1)：55-58.

[28] LI Z L，OPENSHAW S.Algorithms for automated line generalization based on a natural principle of objective generalization[J].International Journal of Geographical Information Systems，1992，6(5)：373-389.

[29] 朱文博.基于多分辨率DEM模型的等高線綜合方法研究[D].太原：太原理工大學(xué)，2002.

[30] BENOIT B M.How long is the coast of Britain[J].Science，1967，156(3775)：636-638.

[31] 朱華，姬翠翠.分形理論及其應(yīng)用[M].北京：科學(xué)出版社，2011.

[32] BARBARA P B.Scale-dependence and self-similarity in cartographic lines[J].Cartographica：The International Journal for Geographic Information and Geovisualization，1989，26(1)：79-100.

[33] PAUL A L，BATTY M.Fractal measurement and line generalization[J].Computers & Geosciences，1989，15(2)：167-183.

[34] KEITH C C.Scale-based simulation of topographic relief[J].Cartography and Geographic Information Science，1988，15(2)：173-181.

[35] 王橋，毋河海.地圖信息的分形描述與自動(dòng)綜合研究[M].武漢：武漢測(cè)繪科技大學(xué)出版社，1998.

[36] 李雯靜，毋河海.地圖綜合中的地圖目標(biāo)自相似性分形衰減研究[J].測(cè)繪科學(xué)，2005，30(3)：21-23.

[37] 李晨釗，舒紅.等高線小波綜合中的閾值選取方法[J].測(cè)繪信息與工程，2006，31(5)：45-46.

[38] 于浩，楊勤科，張曉萍，等.基于小波多尺度分析的DEM數(shù)據(jù)綜合及尺度轉(zhuǎn)換[J].地理與地理信息科學(xué)，2009，25(4)：12-16.

[39] 于浩，楊勤科，張曉萍，等.基于小波多尺度分析的DEM數(shù)據(jù)綜合研究[J].測(cè)繪科學(xué)，2008，33(3)：93-95.

[40] 吳紀(jì)桃，王橋.小波分析在GIS線狀數(shù)據(jù)圖形簡(jiǎn)化中的應(yīng)用研究[J].測(cè)繪學(xué)報(bào)，2000，29(1)：71-75.

[41] 王玉海，朱長(zhǎng)青，游雄，等.基于B樣條小波的等高線數(shù)據(jù)簡(jiǎn)化[J].測(cè)繪科學(xué)，2003，28(2)：23-25.

[42] 朱長(zhǎng)青，王玉海，李清泉，等.基于小波分析的等高線數(shù)據(jù)壓縮模型[J].中國(guó)圖象圖形學(xué)報(bào)，2004，9(7)：841-845.

[43] 王明常，王宇，邢立新，等.多分辨率小波分析在線狀要素簡(jiǎn)化中的應(yīng)用[J].地理與地理信息科學(xué)，2005，21(5)：38-40.

[44] 黃娟，程耀東.多分辨率小波分析在GIS線狀要素簡(jiǎn)化中的應(yīng)用[J].唐山學(xué)院學(xué)報(bào)，2010，23(6)：13-16.

[45] 朱強(qiáng)，武芳，翟仁健.基于通視性原理的等高線化簡(jiǎn)算法研究[J].中國(guó)圖象圖形學(xué)報(bào)，2009，14(2)：359-364.

[46] 費(fèi)立凡.地形圖等高線成組綜合的試驗(yàn)[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，1993，8(S1)：7-21.

[47] 毛可標(biāo)，陳向東.地形結(jié)構(gòu)線自動(dòng)生成方法研究[J].測(cè)繪科學(xué)，1995，(3)：12-18.

[48] 余生晨，劉洪.山脊線與山谷線的計(jì)算機(jī)自動(dòng)檢測(cè)[J].中國(guó)圖象圖形學(xué)報(bào)，1999，4(8)：693-696.

[49] 曲均浩，程久龍，崔先國(guó).垂直剖面法自動(dòng)提取山脊線和山谷線[J].測(cè)繪科學(xué)，2007，32(5)：30-31.

[50] 楊英偉，羅娟.聯(lián)合Delaunay三角網(wǎng)的等高線群綜合研究[J].地理空間信息，2010，(2)：34-36.

[51] 吳凡，粟衛(wèi)民.顧及地形特征的等高線拓?fù)淇臻g關(guān)系表達(dá)[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版)，2006，39(3)：140-144.

[52] 吳凡，粟衛(wèi)民.基于網(wǎng)絡(luò)模型的等高線群拓?fù)潢P(guān)系表示[J].遙感學(xué)報(bào)，2007，11(2)：145-151.

[53] 朱強(qiáng)，武芳，錢海忠，等.基于空間認(rèn)知的等高線樹的構(gòu)建[J].信息工程大學(xué)學(xué)報(bào)2011，12(4)：458-462.

[54] 郭慶勝.線狀要素圖形綜合的漸進(jìn)方法研究[J].武漢測(cè)繪科技大學(xué)學(xué)報(bào)，1998，23(1)：52-56.

[55] 郭慶勝，于江佩，徐崢.等高線圖形綜合策略的研究[J].測(cè)繪信息與工程，2006，30(6)：16-17.

[56] HIND T，PARROT J F，ALVAREZ R.DEM generation by contour line dilation[J].Computers & Geosciences，1999，25(7)：775-783.

[57] 李志林，朱慶.數(shù)字高程模型[M].武漢：武漢測(cè)繪科技大學(xué)出版社，2000.

[58] SANJAY C，SUN P.SLOM：A new measure for local spatial outliers[J].Knowledge and Information Systems，2006，9(4)：412-429.

[59] HANS-PETER K，ZIMEK A.Angle-based outlier detection in high-dimensional data[C].Proceedings of the 14th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining，ACM，2008：444-452.

[60] 張海堂，羅睿，郭建星，等.基于三角網(wǎng)漸進(jìn)式簡(jiǎn)化的等高線多尺度綜合[J].測(cè)繪信息與工程，2004，29(5)：11-13.

[61] 朱文博，畢如田，郭磐石，等.一種基于多分辨率模型簡(jiǎn)化算法的等高線自動(dòng)綜合方法研究[J].山西農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2008，28(3)：332-337.

[62] 于江佩.基于ARCGIS的等高線與河流的自動(dòng)綜合[D].武漢：武漢大學(xué)，2005.

[63] 龍毅，曹陽，沈婕，等.基于約束D-TIN的等高線簇與河網(wǎng)協(xié)同綜合方法[J].測(cè)繪學(xué)報(bào)，2011，40(3)：379-385.