鄧凱 朱世泉 王浩 岳梓晨 鄧虎

摘要：紙質(zhì)地形圖在掃描過(guò)程中，由于掃描儀器操作或紙質(zhì)地形圖質(zhì)量問(wèn)題，會(huì)出現(xiàn)底紋、色彩失真及孔洞、斷點(diǎn)等噪聲。為避免這些噪聲干擾掃描地形圖等高線的自動(dòng)矢量化，利用Photoshop對(duì)掃描后的電子地形圖進(jìn)行降噪處理，得到光滑連續(xù)的等高線柵格圖，繼而使用ArcMap的ArcScan工具自動(dòng)提取等高線，最后，從等高線質(zhì)量和處理效率角度出發(fā)，對(duì)提出的方法與手動(dòng)提取以及未經(jīng)過(guò)降噪處理的自動(dòng)提取技術(shù)方案進(jìn)行對(duì)比分析。結(jié)果表明：該方法不僅能夠保證等高線的數(shù)據(jù)質(zhì)量，作業(yè)效率也有顯著提升。

關(guān)鍵詞：地形圖，等高線，矢量化，Photoshop軟件，ArcMap軟件

中圖分類號(hào)：P208?文獻(xiàn)標(biāo)志碼：A?文章編號(hào)：2095-5383（2020）01-0005-08

Abstract：In the process of scanning paper topographic maps，due to the operation of scanners or the quality problems of paper topographic maps，there will be background lines，color distortion，holes，breakpoints and other noises.?In order to avoid these noises from interfering with the automatic vectorization of contour lines during scanning the topographic maps，a smooth and continuous contour grid was obtained by using Photoshop to reduce the noise of scanned electronic topographic maps，and then contour lines were automatically extracted by ArcMaps ArcScan tool.?Finally，from the point of view of contour quality and processing efficiency，the contour lines were processed.?The method proposed in this paper was compared with manual extraction and automatic extraction without de-noising.?The results show that this method can not only ensure the data quality of contour lines，but also significantly improve the efficiency of operation.

Keywords：topographic map，contour line，vectorization，Photoshop，ArcMap

等高線是能夠正確表示地貌形態(tài)的重要數(shù)據(jù)源之一。實(shí)際工程中，等高線數(shù)據(jù)來(lái)源有紙質(zhì)地形圖掃描后矢量化、測(cè)繪離散高程點(diǎn)生成等，其中紙質(zhì)地形圖掃描方式較為普遍。然而，由于掃描儀操作不當(dāng)或是紙質(zhì)地形圖本身質(zhì)量問(wèn)題，致使掃描得到的電子地形圖出現(xiàn)孔洞、斷點(diǎn)等噪聲及底紋、色彩失真等問(wèn)題，極大地影響了地形圖等高線自動(dòng)矢量化作業(yè)的效率與質(zhì)量。

針對(duì)該問(wèn)題，目前主要是利用工具軟件或是通過(guò)改進(jìn)算法來(lái)提高等高線處理的質(zhì)量和效率。在利用軟件操作方面，較早的有董敏等[1]用ArcGIS自動(dòng)矢量化提取等高線，薛濱瑞[2]用CorelDraw和ArcGIS制作地形圖等高線，但這兩種方案易導(dǎo)致矢量化線條變形，王麗等[3]用CASSCAN5_0軟件掃描地形圖矢量化，方式是單一的手動(dòng)追蹤，工作量大，效率低，不適合大批量的矢量化作業(yè)，何立恒等[4]利用AutoCAD提取矢量化地形圖要素，但該方案僅適用于1∶2 000、1∶1 000比例尺的地形圖，羅虎[5]使用Photoshop與ArcGIS提取地形圖等高線，該方法在降噪處理過(guò)程中單一地使用人機(jī)交互方式[6]剔除干擾要素，并未發(fā)揮Photoshop對(duì)顏色和圖像的處理優(yōu)勢(shì)。在算法研究方面，主要有爨瑩等[7]通過(guò)圖像智能跟蹤算法開(kāi)發(fā)出地質(zhì)圖件智能數(shù)字化系統(tǒng)，主要針對(duì)噪聲較小的彩色地形圖來(lái)提取交叉點(diǎn)畸變較小的拓?fù)涔羌?，劉人午等[8]根據(jù)柵格線體的形態(tài)特征設(shè)計(jì)通過(guò)對(duì)柵格線體的序貫推掃及截?cái)鄴呙璜@取截?cái)嗵幍闹行南袼氐膾呙杓?xì)化算法。Naccache等[9]提出細(xì)化算法來(lái)提取等高線，雖然速度快、效果好，但算法對(duì)于掃描地形圖的質(zhì)量要求較高，且實(shí)現(xiàn)方式成本較大，不利于普及。

本文參考了鄭華利[10]提出的像素空間關(guān)系信息與顏色信息相結(jié)合的地形圖分色算法，通過(guò)噪聲分析，利用Photoshop軟件在圖像處理中的優(yōu)勢(shì)，對(duì)掃描地形圖進(jìn)行降噪處理，并對(duì)地形圖中的等高線進(jìn)行了提取和優(yōu)化，最后在ArcMap軟件中進(jìn)行等高線自動(dòng)化處理。

1?數(shù)據(jù)來(lái)源與噪聲分析

1.1?數(shù)據(jù)來(lái)源

選取具有代表性的、不同復(fù)雜度的3幅1∶10 000柵格電子地形圖，均采用高斯克呂格西安1980坐標(biāo)系，采用分辨率為600×1 200，灰度系數(shù)在1.4～1.8之間，24位色彩模式掃描儀掃描處理，獲得的地形圖圖像清晰易辨。3幅地形圖等高線分布由稀疏到密集，矢量化工作由易到難，各地形圖的基本情況如表1所示，地形圖如圖1～3所示。

1.2?噪聲分析

1）底紋：底紋與掃描儀感光元件、透鏡組等本身特性[11]有關(guān)，客觀上很難消除。而且掃描儀的分辨率越低，底紋現(xiàn)象越嚴(yán)重。所以在掃描紙質(zhì)地形圖時(shí)，對(duì)設(shè)備與操作技術(shù)的要求不容忽視。底紋較容易判斷，參雜在白色背景區(qū)域內(nèi)的粉紅色、排布有規(guī)律且斷斷續(xù)續(xù)的線，就是底紋，如圖4所示。

2）與等高線相交的水系：水系是地形圖組成的基本要素，如圖5所示。根據(jù)自動(dòng)分色原理[10]，水系線體現(xiàn)為藍(lán)色，在二值化時(shí)，水系線若不去除，會(huì)被識(shí)別為等高線的一部分，給自動(dòng)矢量化等高線造成極為嚴(yán)重的影響?？筛鶕?jù)圖斑顏色特征與等高線明顯不同的特點(diǎn)，利用Photoshop的顏色識(shí)別功能去除。

3）與等高線相交的圖斑：經(jīng)緯線、道路、注記，這些圖斑之所以難以去除，是因?yàn)檫@些圖斑與等高線顏色相近，如圖6所示。該誤差可以采用人機(jī)交互剔除[12]或在ArcMap中用像素清除。

4）等高線斷裂：如圖7所示，紙質(zhì)地圖掃描成柵格圖像時(shí)，由于儀器操作和圖像顏色分布不均勻，都會(huì)使線條會(huì)出現(xiàn)孔洞、斷點(diǎn)等噪聲，但正確的等高線應(yīng)該是連續(xù)、獨(dú)立的?？赏ㄟ^(guò)在Photoshop中擴(kuò)展與收縮填補(bǔ)孔洞。

2?數(shù)據(jù)處理與結(jié)果分析

2.1?處理流程

為了方便后續(xù)對(duì)比分析工作，本次實(shí)踐選用3幅1∶10 000地形圖，并制定了A、B、C 3種方案，其中方案C為本文采用的技術(shù)方案，方案A、B用于對(duì)比，具體說(shuō)明如表2所示。

處理流程如圖8所示。第1步，將紙質(zhì)地形圖掃描得到電子?xùn)鸥竦匦螆D。第2步，利用Photoshop分離原始圖的色調(diào)，賦予地形圖中各類要素不同的特征顏色，重點(diǎn)要突出等高線的紅棕色，再通過(guò)對(duì)顏色的選取與編輯剔除非等高線圖斑噪聲，對(duì)于等高線線條質(zhì)量問(wèn)題，通過(guò)平滑曲線和調(diào)整閾值等操作優(yōu)化線條質(zhì)量。第3步，將符合要求的影像加載到ArcMap中，經(jīng)過(guò)配準(zhǔn)、創(chuàng)建要素等工作，激活Editor和Arcscan工具欄，即可進(jìn)行自動(dòng)矢量化。

2.2?Photoshop處理

2.2.1?等高線分色處理

根據(jù)自動(dòng)分色原理，由于地形圖上等高線與注記、水系等要素相比具有紅棕色特征，地形圖中紅棕色等高線的飽和度分量值相對(duì)黑色、藍(lán)色較高，據(jù)此可以有效提取等高線。紅棕色的RGB值為144、164、96，通過(guò)RGB的調(diào)整，增強(qiáng)等高線的棕色特征。但要徹底分離出等高線，還需淡化其他雜色。HSV是一種直觀而且比較符合人體視覺(jué)對(duì)色彩感知特征的色域空間，H是色調(diào)、S是飽合度、V是明度，通過(guò)判別各個(gè)分量的閾值，修改飽合度，描繪出等高線的輪廓。分色處理的效果如圖9所示。

2.2.2?底紋噪聲處理

底紋的產(chǎn)生主要是因?yàn)閽呙鑳x的元器件材料本身、掃描頭的機(jī)械抖動(dòng)及掃描儀中光源的光照均勻程度等因素。同時(shí)，掃描儀的光電穩(wěn)定性能直接影響所得到的RGB三分色單版圖像的灰度，掃描儀濾光片的非理想性也可引起顏色誤差。色階是表示圖像亮度強(qiáng)弱的指數(shù)標(biāo)準(zhǔn)，由于底紋通常顯示為粉紅色，色階的值越小，深色越深、淺色越淺，粉紅色逐漸向白色（背景色）變化，有利于去除底紋。故利用色調(diào)分離工具，調(diào)整底圖的色階來(lái)改變圖像的亮暗程度，從而達(dá)到去除底紋的目的。但色階值過(guò)小，等高線容易出現(xiàn)斷裂，所以，通常將色階設(shè)為2～4。

底紋去除效果如圖10所示。

2.2.3?非等高線噪聲處理

由于色彩的散射作用[11]，掃描地形圖圖像含有成千上萬(wàn)種不同顏色，但都是以黑、棕、綠、藍(lán)4色為主衍生而來(lái)，與其他圖像的分色相比，分色范圍更加明確。下面以水系線為例，說(shuō)明如何去除水系、注記等顏色與等高線明顯不同的圖斑。因水系用藍(lán)色表示，若采集到圖中藍(lán)色的樣本顏色，將其填充為白色，即讓水系要素與背景融為一體，最終只區(qū)別于背景色的只有等高線。在使用魔棒工具選取水系的顏色時(shí)，注意調(diào)整容差，容差越大，樣本選取的范圍也越大。在處理水系線時(shí)，將容差選為32～60。處理完水系后，有很多小且多的短微線以及線條斷裂或色調(diào)分離未完全產(chǎn)生的點(diǎn)，在后續(xù)步驟中通過(guò)Photoshop的選擇填充功能去除。

去除效果如圖11所示。

2.2.4?等高線線條質(zhì)量?jī)?yōu)化

在掃描過(guò)程中，由于圖像邊緣效應(yīng)[11]存在，地形圖上要素從一種顏色轉(zhuǎn)變到另一種顏色時(shí)會(huì)產(chǎn)生連續(xù)的過(guò)渡顏色，比如受掃描儀濾光片的實(shí)際透射性質(zhì)影響，等高線邊緣經(jīng)常表現(xiàn)為紅黑混合，通過(guò)色調(diào)分離使等高線呈現(xiàn)出與其他顏色特征相近的圖斑不一樣的像素，方便后續(xù)操作。

孔洞、斷點(diǎn)和線線相交，使用人機(jī)交互方式剔除[12]。用魔棒選中線條，結(jié)合閾值設(shè)置與填充工具，通過(guò)“擴(kuò)展”和“收縮”改善線條的斷裂與相交：“擴(kuò)展”可以有規(guī)則地將等高線的孔洞、斷點(diǎn)擴(kuò)張，使孔洞左右的2個(gè)斷點(diǎn)聚攏在一起，但也會(huì)導(dǎo)致距離過(guò)近的等高線相交在一起，“收縮”可以讓原本粗的線條在一定程度上變細(xì)進(jìn)而線線分離。

但在轉(zhuǎn)換色階時(shí)，部分毛刺及斷線顏色變淺，難以觀察?？紤]到這些噪聲的像素大小與等高線相差甚遠(yuǎn)，需要使用色階工具進(jìn)行分離，設(shè)定一個(gè)色階閾值，以此值為標(biāo)準(zhǔn)，比閾值亮的像素轉(zhuǎn)換為白色，比閾值暗的像素轉(zhuǎn)換為黑色，將原來(lái)的彩色圖像轉(zhuǎn)換為高對(duì)比度的黑白圖像。

本文根據(jù)用到的地形圖中等高線的亮度判斷，其色階位于中間值附近。通過(guò)觀察毛刺、碎線以及等高線的邊緣變化，當(dāng)毛刺、碎線成為極易觀察的黑色，同時(shí)等高線斷裂不明顯的時(shí)候，此時(shí)閾值最為合適。經(jīng)過(guò)反復(fù)試驗(yàn)，將色階直方圖的閾值定義在中值128處比較合適，最后通過(guò)調(diào)整邊緣的命令，修改形狀的半徑大小，剔除細(xì)碎的像素。

由于文字、注記等因素干擾等高線線條邊緣走向，產(chǎn)生毛刺，在自動(dòng)矢量化工作中會(huì)產(chǎn)生短折線。本文通過(guò)平滑處理毛刺[13]，光滑等高線的邊緣，結(jié)合閾值設(shè)置與填充工具，使等高線達(dá)到如圖12所示的效果。

2.3?ArcMap矢量化處理

矢量化的產(chǎn)品都有屬性信息，因此沒(méi)有地理信息的柵格圖像是不能被矢量化的。首先通過(guò)地形圖左上角的標(biāo)注，找到該幅地形圖的分帶，在經(jīng)緯線交點(diǎn)處判斷坐標(biāo)并配準(zhǔn)。方案A與方案B中配準(zhǔn)時(shí)的點(diǎn)文件注意保存。由于方案C中的經(jīng)緯線和坐標(biāo)模糊不清，實(shí)行方案C時(shí)需要在配準(zhǔn)工作中要將方案A、B的點(diǎn)文件加載進(jìn)去。創(chuàng)建可供編輯的等高線要素集，給ArcScan提供可供編輯與保存數(shù)據(jù)的圖層。

自動(dòng)矢量化的ArcScan只能識(shí)別二值化的圖像，二值化就是在Spatial Analyst工具中選擇重分類，把圖中的背景色與等高線（包括未完全清理的干擾要素）劃分成如圖13所示的2個(gè)波段，將圖像的灰度值設(shè)為2個(gè)值，使整幅圖像顯出明顯的2種顏色的特征，突顯出目標(biāo)的輪廓。

由于在Photoshop圖中的文字、注記等地物符號(hào)清理得不徹底，且這些干擾因素的像素與等高線不同，干擾像元主要聚集在圖框和要素密集處，利用ArcScan選擇像素清理像元。

在ArcScan中生成要素后，3種方案的矢量化完成后最終效果如圖14所示。

由圖14可知，方案A手動(dòng)矢量化的線條與原始等高線吻合度較高，極少出現(xiàn)錯(cuò)誤，但效率低，方案B是未經(jīng)過(guò)降噪處理自動(dòng)提取的等高線，發(fā)生了距離較近的線條易相交、顏色不均勻的線條斷裂等情況、且生成了大量雜點(diǎn)，給后續(xù)的拓?fù)錂z查與糾錯(cuò)帶來(lái)大量難題，方案C效率高，質(zhì)量較方案B有了明顯提高，雖然在少數(shù)噪聲清理不完全的地方發(fā)生線線相交與線條斷裂，但整體上還原了原始圖等高線的走向。

3?效率與質(zhì)量分析

3.1?效率對(duì)比

為了證實(shí)效率有所提高，選用了3幅地形圖，適合自動(dòng)矢量化與手動(dòng)矢量化，設(shè)定每次矢量化的范圍為一整張地圖：在地形圖內(nèi)由同一起點(diǎn)開(kāi)始打點(diǎn)或框選，矢量化相同范圍的地形圖等高線，通過(guò)耗時(shí)長(zhǎng)短來(lái)比較3種方案的效率，數(shù)據(jù)如表3所示。

方案A需要操作人員的判斷和描點(diǎn)，且通常需要將數(shù)據(jù)框比例尺放大至1∶4 000以上，在等高線分布密集的區(qū)域甚至要將數(shù)據(jù)框比例尺增大到1∶2 000。數(shù)據(jù)框比例尺越大，每次操作之后矢量化范圍就越小。但方案A在后續(xù)的拓?fù)錂z查及處理中耗時(shí)較少，考慮到方案A省略了Photoshop預(yù)處理步驟，遇到等高線排布簡(jiǎn)單的地形圖時(shí)，方案A的手動(dòng)矢量化與拓?fù)錂z查及糾錯(cuò)步驟節(jié)省了更多的時(shí)間，因此遇到線條簡(jiǎn)單的地形圖時(shí)，方案A的效率比方案B高。

自動(dòng)矢量化一般不用考慮比例尺，只要保證選中的區(qū)域質(zhì)量符合要求，對(duì)于同樣大小的區(qū)域，自動(dòng)矢量化可以在選中的區(qū)域同時(shí)矢量化多條等高線，而手動(dòng)追蹤矢量化每次最多只能矢量化一條等高線其中的某一段，對(duì)于相同長(zhǎng)度一條等高線時(shí)，基于ArcMap所提供的算法支持，自動(dòng)矢量化可以通過(guò)計(jì)算機(jī)在該條等高線范圍內(nèi)立刻生成要素，手動(dòng)矢量化則需要操作人員根據(jù)等高線的走向逐個(gè)打點(diǎn)。

根據(jù)以上信息，在地形復(fù)雜的區(qū)域，自動(dòng)矢量化的速度是手動(dòng)追蹤矢量化的3倍以上。地形圖等高線分布越密集，自動(dòng)矢量化能夠節(jié)省的時(shí)間就越多，尤其在地形復(fù)雜且的地形圖Ⅲ中，方案C較方案A節(jié)省了近40 h，而對(duì)于地形簡(jiǎn)單的地形圖Ⅰ，方案B拓?fù)錂z查和糾錯(cuò)的時(shí)間與方案A的手動(dòng)追蹤打點(diǎn)相差不大，但隨著難度指標(biāo)增高，3種方案的效率優(yōu)勢(shì)差異愈加明顯。但方案B容易出錯(cuò)，比如線線相交未處理完畢，矢量化的等高線可能會(huì)出現(xiàn)交點(diǎn)，給拓?fù)錂z查及糾錯(cuò)工作帶來(lái)了大量難題。因此，方案C是效率最高的方案。

3.2?質(zhì)量對(duì)比

等高線的質(zhì)量指標(biāo)可以參考矢量化等高線與原始地形圖等高線的吻合程度[14]，依據(jù)湯國(guó)安等[15]提出的等高線套合積差A(yù)P可以評(píng)定該吻合程度。

設(shè)定手動(dòng)矢量化與原始地形圖等高線之間的AP=0，以手動(dòng)矢量化結(jié)果為準(zhǔn)，計(jì)算自動(dòng)矢量化等高線與手動(dòng)矢量化等高線之間的套合積差，評(píng)定自動(dòng)矢量化等高線與原始地形圖等高線的吻合程度。如果測(cè)得自動(dòng)矢量化生成等高線與手動(dòng)追蹤矢量化的套合積差在0.05以內(nèi)，則說(shuō)明其吻合度較高，矢量化等高線質(zhì)量合格。

單條等高線套合差：

（1）整幅地形圖的等高線套合積差：

（2）其中：Si為該條自動(dòng)矢量化等高線與手動(dòng)追蹤矢量化等高線復(fù)合后細(xì)碎多邊形的面積，M為相鄰2條間曲線的面積，m為等高線的條數(shù)。首先計(jì)算手動(dòng)追蹤矢量化與自動(dòng)矢量化生成的兩種等高線相交而產(chǎn)生的面積與相鄰的每3條等高線之間的面積之和的比值，再收集所有的比值計(jì)算出AP值。經(jīng)過(guò)提取半距等高線和分帶統(tǒng)計(jì)與計(jì)算，得到各方案AP值，如圖15所示。

由式（1）、（2）可知，P值表示的是每一條自動(dòng)矢量化等高線的吻合度，P值越小，吻合度越高，越能還原原始柵格地形圖的走向，AP值的大小可以反映出該矢量化等高線方案在全圖范圍內(nèi)的各種地形是否會(huì)產(chǎn)生大幅度波動(dòng)。

方案B在全圖范圍內(nèi)等高線受到干擾變形嚴(yán)重，即便在拓?fù)錂z查及糾錯(cuò)之后，偏離了原圖等高線的矢量化線條也不便糾正，而方案C的AP值均為0.03，矢量化線條與原始地形圖吻合度高，質(zhì)量好。

4?結(jié)論與討論

隨著地圖學(xué)的快速發(fā)展，如何高效提取地圖上有價(jià)值的信息，使掃描紙質(zhì)地形圖矢量化成為一項(xiàng)具有挑戰(zhàn)性的課題。本文分析了紙質(zhì)地形圖等高線矢量化的研究現(xiàn)狀，以快速高質(zhì)量提取等高線為目的，制訂了3種方案，分別在3種具有代表性的地形圖上進(jìn)行等高線的提取，最后通過(guò)測(cè)算矢量化作業(yè)耗時(shí)與套合積差來(lái)反映等高線提取的效率與質(zhì)量，得出了以下成果與結(jié)論：

1）地形圖等高線噪聲的根本來(lái)源在于儀器掃描地圖的操作以及地形圖原圖質(zhì)量問(wèn)題，可以通過(guò)優(yōu)化儀器性能、提高紙質(zhì)地形圖的印刷工藝使掃描得到的電子?xùn)鸥竦匦螆D的分辨率有所提高，為降噪工作減輕負(fù)擔(dān);

2）利用Photoshop與ArcMap自動(dòng)矢量化提取等高線的降噪技術(shù)較以往有了改進(jìn)，主要依據(jù)等高線飽和度分量值與其他要素的差別，通過(guò)對(duì)彩色電子?xùn)鸥竦匦螆D的分色來(lái)實(shí)現(xiàn)，可以在復(fù)雜地形圖環(huán)境下，有效地凸出等高線的輪廓，比目前較通用的人機(jī)交互剔除方式節(jié)省了大量時(shí)間，提高了矢量化等高線工作降噪步驟的效率;

同時(shí)，本文所闌述的方案也有不足之處，比如在等高線線條質(zhì)量?jī)?yōu)化過(guò)程中，閾值的調(diào)整可能額外生成雜點(diǎn)，但在本次實(shí)驗(yàn)中并未對(duì)矢量化操作結(jié)果造成嚴(yán)重影響，另外，擴(kuò)展、收縮操作是在全圖范圍內(nèi)同時(shí)實(shí)現(xiàn)的等量變化，不能針對(duì)某一處斷點(diǎn)或孔洞給出最優(yōu)的變化數(shù)值，導(dǎo)致對(duì)排列過(guò)于緊密、或者粗細(xì)分布不均勻的線條處理不當(dāng)。所以某些微小的噪聲可以依據(jù)實(shí)際情況添加手動(dòng)操作，或者在拓?fù)錂z查后逐個(gè)糾錯(cuò)。

如果有更好的方案還需依賴通過(guò)算法或其他更先進(jìn)的圖像處理軟件實(shí)現(xiàn)高效率高質(zhì)量的降噪工作。后續(xù)研究將立足于簡(jiǎn)化降噪工序的流程和提高降噪工序的效率與質(zhì)量，找到更適用于所有地形圖的矢量化方案。

參考文獻(xiàn)：

[1]董敏，孫寶生，陳川，等.基于ArcGIS的柵格圖像自動(dòng)矢量化的實(shí)現(xiàn)新方法[J].西部探礦工程，2011，23（2）：139-142.

[2]薛濱瑞，余正軍.基于CorelDraw和GIS軟件的等高線地形圖制作：以翠華山地質(zhì)公園地形圖為例[C]//中國(guó)地質(zhì)學(xué)會(huì)旅游地學(xué)與地質(zhì)公園研究分會(huì)第25屆年會(huì)暨張家界世界地質(zhì)公園建設(shè)與旅游發(fā)展戰(zhàn)略研討會(huì)，2010：5.

[3]王麗，薛志勤.淺談CASSCAN5.0軟件在地形圖掃描矢量化中的應(yīng)用[J].計(jì)算機(jī)光盤軟件與應(yīng)用，2014，17（6）：51-52.

[4]何立恒，鮑其勝.AutoCAD R14環(huán)境下的地形圖矢量化[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2002（2）：63-65.

[5]羅虎.Adobe Photoshop CC 2015結(jié)合ArcMap10.2.2在地形圖矢量化中的應(yīng)用[J].綠色科技，2018（14）：229-231.

[6]朱春節(jié)，韓雪培，封殿波.土地利用圖自動(dòng)矢量化方法探討[J].測(cè)繪通報(bào)，2008（7）：50-53.

[7]爨瑩，袁通路，徐永高，等.地質(zhì)等值線圖矢量化算法的設(shè)計(jì)與實(shí)現(xiàn)[J].系統(tǒng)仿真學(xué)報(bào)，2003（11）：1565-1566，1617.

[8]劉人午，李燕.掃描細(xì)化算法的地圖自動(dòng)矢量化研究[J].測(cè)繪學(xué)報(bào)，2012，41（2）：309-314.

[9]NACCACHE N J，SHINGHAL R.?SPTA：a proposed algorithm for thinning binary patterns[J].?IEEE Transactions on Systems，Man，and Cybernetics，1984，SMC-14（3）：409-418.

[10]鄭華利.彩色地形圖的自動(dòng)識(shí)別與矢量化研究[D].南京：南京理工大學(xué)，2004.

[11]周蓓蓓.彩色地形圖點(diǎn)狀地物符號(hào)的提取與識(shí)別[D].蘇州：蘇州大學(xué)，2008.

[12]紀(jì)洋，武文波，楊曉偉.等高線自動(dòng)矢量化的后處理[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào)（自然科學(xué)版），2008（S1）：37-39.

[13]韋建波，塔西甫拉提特依拜，田源，等.一種改善地形圖等高線自動(dòng)矢量化精度的方法[J].測(cè)繪通報(bào)，2008（6）：61-64.

[14]張曉偉，沈軍.利用等高線套合差評(píng)定DEM質(zhì)量的應(yīng)用研究[J].青海師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，29（2）：26-29.

[15]湯國(guó)安，陶旸，王春.等高線套合差及在DEM質(zhì)量評(píng)價(jià)中的應(yīng)用研究[J].測(cè)繪通報(bào)，2007（7）：65-67.