dem粗差剔除方法研究-畢業(yè)論文

1、<p>　　南陽(yáng)師范學(xué)院20XX屆畢業(yè)生</p><p><b>　　畢業(yè)論文（設(shè)計(jì)）</b></p><p>　　題 目： DEM粗差剔除方法研究 </p><p>　　完 成 人： </p><p>　　班 級(jí)：

2、 </p><p>　　學(xué) 制： </p><p>　　專 業(yè)： 測(cè)繪工程 </p><p>　　指導(dǎo)教師： </p><p&

3、gt;　　完成日期： </p><p><b>　　目錄</b></p><p><b>　　摘要(1)</b></p><p><b>　　1引言(1)</b></p><p>　　1.1 選題意義(1)&l

4、t;/p><p>　　1.2 研究現(xiàn)狀(2)</p><p>　　2 DEM粗差探測(cè)與剔除方法研究(4)</p><p>　　2.1 基于三維可視化的DEM粗差探測(cè)與剔除(4)</p><p>　　2.2 基于高程信息的不規(guī)則分布數(shù)據(jù)粗差檢測(cè)與剔除(4)</p><p>　　2.3 檢測(cè)規(guī)則分布數(shù)據(jù)中粗差的算法(

5、5)</p><p>　　2.3.1基于參數(shù)統(tǒng)計(jì)的粗差檢測(cè)(5)</p><p>　　2.3.2基于參數(shù)統(tǒng)計(jì)的粗差檢測(cè)(6)</p><p>　　2.3.3基于主成分分析的粗差檢測(cè)(8)</p><p>　　2.4 檢測(cè)不規(guī)則分布數(shù)據(jù)中粗差的算法(9)</p><p>　　2.4.1基于點(diǎn)方式的算法(9)&l

6、t;/p><p>　　2.4.2基于粗差簇群的算法(10)</p><p>　　2.5 基于等高線拓?fù)潢P(guān)系的粗差檢測(cè)與剔除(10)</p><p>　　3基于趨勢(shì)面的粗差探測(cè)與剔出方法研究(11)</p><p>　　3.1 獲取數(shù)據(jù)采集后形成的原始DEM數(shù)據(jù)(11)</p><p>　　3.2 數(shù)據(jù)預(yù)處理(12

7、)</p><p>　　3.3粗差檢測(cè)及修正(13)</p><p>　　3.4 實(shí)驗(yàn)分析(14)</p><p>　　4總結(jié)與展望(17)</p><p><b>　　參考文獻(xiàn)(18)</b></p><p>　　Abstract(18)</p><p>　　D

8、EM粗差剔除方法研究</p><p>　　摘要: 數(shù)字高程模型（DEM）作為4D產(chǎn)品（DEM-數(shù)字高程模型、DOM-數(shù)字正射影像、DLG-數(shù)字線劃圖、DRG-數(shù)字柵格圖）之一，是一種對(duì)地球表面的數(shù)字化描述和模擬，是地球空間數(shù)據(jù)基礎(chǔ)設(shè)施的重要組成部分，是建立地形高程數(shù)據(jù)庫(kù)和各類GIS（地理信息系統(tǒng)）庫(kù),進(jìn)行地形定量分析等方面所必需的基礎(chǔ)數(shù)據(jù)。數(shù)字高程模型的應(yīng)用越來(lái)越廣泛，可用于地學(xué)分析、二維地理空間上連續(xù)分布并逐漸

9、變化的各種非高程屬性數(shù)據(jù)的建模與分析。DEM的數(shù)據(jù)來(lái)源方式有地面測(cè)量（利用自動(dòng)記錄的測(cè)距經(jīng)緯儀在野外實(shí)測(cè)）、現(xiàn)有地圖數(shù)字化（利用數(shù)字化儀對(duì)已有地圖上的信息如等高線、地形線等進(jìn)行數(shù)字化，目前常用的數(shù)字化儀有手扶跟蹤數(shù)字化儀與掃描數(shù)字化儀）、空間傳感器（利用GPS等進(jìn)行數(shù)據(jù)采集）、數(shù)字?jǐn)z影測(cè)量方法（這是DEM數(shù)據(jù)采集最常用的一種方法）。事實(shí)上,不管采用何種測(cè)量方法,測(cè)址數(shù)據(jù)總會(huì)包含各種各樣的誤差, 可將誤差分為三種,即系統(tǒng)誤差、偶然誤差和粗

10、差。同系統(tǒng)誤差、偶然誤差相比,粗差實(shí)際上是一種錯(cuò)誤,他們?cè)跍y(cè)量中出現(xiàn)的可能性一般比較小。但是與前兩種誤差相比,粗差對(duì)數(shù)字高程模型所反映的空間變化的扭曲更為嚴(yán)重,因而粗差是影響DEM質(zhì)量的</p><p>　　關(guān)鍵詞：數(shù)字高程模型 ； 粗差探測(cè) ； 粗差剔除</p><p><b>　　1引言</b></p><p><b>　　1.1

11、 選題意義</b></p><p>　　數(shù)字高程模型(Digital Elevation Model,簡(jiǎn)稱DEM)是以數(shù)字的形式按一定結(jié)構(gòu)組織在一起，表示實(shí)際地形特征空間分布的模型，也是地形形狀大小和起伏的數(shù)字描述，由一系列地面點(diǎn)X，Y位置及其相聯(lián)系的高程Z組成，是各種信息的載體，是地理信息系統(tǒng)(Geographic Information System，GIS)技術(shù)中最重要的內(nèi)容，是空間數(shù)據(jù)基礎(chǔ)設(shè)施

12、的重要組成部分，在生產(chǎn)中具有很高的利用價(jià)值。首先，它能夠反映區(qū)域內(nèi)的地形條件，為各用圖部門提供地形基礎(chǔ)，其次，用它制作沙盤，具有快速、簡(jiǎn)便、精確的優(yōu)點(diǎn)，可用于軍事指揮和模型演示；還可以用于農(nóng)業(yè)部門的農(nóng)田水利規(guī)劃，水利部門的洪水淹沒(méi)損失估算及水利建設(shè)的土方量計(jì)算，交通、建筑等各建設(shè)部門的選址、規(guī)劃，通訊部門信號(hào)覆蓋范圍的規(guī)劃與計(jì)算，地址、勘探部門地形的分析，各旅游景點(diǎn)的規(guī)劃及土地資源調(diào)查；也能用于二維地理空間上連續(xù)分布并逐漸變化的各種非高

13、程屬性數(shù)據(jù)的建模與分析上?？梢哉f(shuō)DEM數(shù)據(jù)具有廣泛的應(yīng)用潛力。目前，由于地理信息系統(tǒng)(GIS)的普及以及空間數(shù)據(jù)基礎(chǔ)設(shè)施的發(fā)展和建設(shè)，DEM作為數(shù)字線化圖(DLG)、數(shù)字高程模型(DEM)、數(shù)字正</p><p><b>　　1.2 研究現(xiàn)狀 </b></p><p>　　DEM是地理空間定位的數(shù)字?jǐn)?shù)據(jù)集合，它最初是美國(guó)麻省理工學(xué)院Miller教授為高速公路的自動(dòng)設(shè)計(jì)

14、于1956年提出來(lái)的，隨著各種相關(guān)技術(shù)的發(fā)展,特別是計(jì)算機(jī)技術(shù)在測(cè)繪方面的應(yīng)用使得測(cè)繪學(xué)科逐漸向數(shù)字化、實(shí)時(shí)處理與多用途的方向發(fā)展[1]。質(zhì)量控制是數(shù)字高程模型(Digital Elevation Model- DEM) 生產(chǎn)的重要環(huán)節(jié)之一。然而粗差對(duì)DEM 數(shù)據(jù)所造成的空間扭曲往往最為嚴(yán)重,有時(shí)能導(dǎo)致DEM及其產(chǎn)品嚴(yán)重失真,甚至完全不能使用,因此很有必要研究一些方法對(duì)DEM的粗差進(jìn)行探測(cè)和修正。要有效地探測(cè)DEM數(shù)據(jù)中存在的粗差,必須

15、在數(shù)據(jù)采集后形成的原始DEM 數(shù)據(jù)中進(jìn)行。原始的DEM數(shù)據(jù)結(jié)構(gòu)有規(guī)則和不規(guī)則兩類。從實(shí)用的角度來(lái)看,在規(guī)則格網(wǎng)的DEM數(shù)據(jù)中探測(cè)粗差相對(duì)簡(jiǎn)單一些,因而研究成果也相對(duì)豐富,如Hannah[2](1981)的基于坡度信息算法,FolicíSimon (1994)的統(tǒng)計(jì)參數(shù)法, 以及Lóper[3](1997) 的主成分分析法等等。如果原始DEM 為不規(guī)則數(shù)據(jù),要轉(zhuǎn)成規(guī)則格網(wǎng)點(diǎn),則需經(jīng)過(guò)數(shù)學(xué)內(nèi)插的方法處理。在此過(guò)程中,原始

16、數(shù)據(jù)中粗差點(diǎn)會(huì)影響到轉(zhuǎn)換后的多個(gè)格網(wǎng)點(diǎn),從而</p><p>　　2 DEM粗差探測(cè)與剔除方法研究</p><p>　　2.1 基于三維可視化的DEM粗差探測(cè)與剔除</p><p>　　通過(guò)地形表面的三維可視化建模來(lái)審查DEM中可疑數(shù)據(jù)點(diǎn)，從而剔除嚴(yán)重影響數(shù)據(jù)質(zhì)量的粗差或者說(shuō)錯(cuò)誤。DEM有著非常適宜于建立3維可視化的特點(diǎn)，采用DEM3維可視化技術(shù)，該方法可以交互式的

17、來(lái)檢查DEM中出現(xiàn)的可疑數(shù)據(jù)，剔除嚴(yán)重影響DEM數(shù)據(jù)質(zhì)量的粗差[4]。一般對(duì)于一個(gè)特定的研究區(qū)域，在三維透視圖上可疑點(diǎn)是否表現(xiàn)為粗差非常直觀，很容易據(jù)此作出正確地判斷。實(shí)際上，由于DEM有著非常適宜于建立三維可視化的特點(diǎn)，所以可以首先通過(guò)目視效果對(duì)粗差進(jìn)行檢測(cè)。通常粗差的地形很不自然，因此在實(shí)際應(yīng)用中，可以首先通過(guò)目視進(jìn)行粗差的檢測(cè)這種方法需要高效可靠的構(gòu)網(wǎng)技術(shù)、快速的交互相應(yīng)效率以及對(duì)異常值敏感的視化圖形，如線框透視圖、暈渲圖等常用的

18、可視化圖形，在技術(shù)層面上，操作的經(jīng)驗(yàn)和工作態(tài)度對(duì)結(jié)果也會(huì)有相當(dāng)大的影響。三維可視化的前提是要建立數(shù)字地面模型，為了保證所分析都基于原始數(shù)據(jù),可選的方法是直接利用原始數(shù)據(jù)建立不規(guī)則三角網(wǎng)絡(luò)模型(Tin)。該方法不利因素在于，一方面需要高效可靠的建模技術(shù)以及可視化處理的策略;另一面，它僅適用于較大粗差的判釋，對(duì)于中小粗差并不敏感。</p><p>　　2.2 基于高程信息的不規(guī)則分布數(shù)據(jù)粗差檢測(cè)與剔除</p&g

19、t;<p>　　呈散亂分布的數(shù)據(jù)點(diǎn)粗差探測(cè)技術(shù)在原理上與規(guī)則格網(wǎng)比較類似，但由于散亂分布的數(shù)據(jù)點(diǎn)的分布特征，坡度信息獲取比較困難，具體實(shí)現(xiàn)上有兩點(diǎn)不同:第一，窗口確定，在規(guī)則格網(wǎng)上采用3*3局部窗口是適宜的，但不規(guī)則分布點(diǎn)的鄰域范圍要進(jìn)行指定，一般可采用窗口尺寸或窗口區(qū)域的采樣點(diǎn)數(shù)量?jī)煞N方式確定。第二，一致性標(biāo)準(zhǔn)確定，規(guī)則格網(wǎng)上比較容易獲取坡度信息，而不規(guī)則分布上獲取坡度信息比較困難，因此，由于高程和坡度同是刻畫地形曲面連

20、續(xù)性的指標(biāo)，在散亂數(shù)據(jù)分布的區(qū)域上，高程信息取代坡度成為一致性標(biāo)準(zhǔn)。在每一個(gè)窗口中，用高程信息計(jì)算統(tǒng)計(jì)指標(biāo)以及確定閡值，方法與規(guī)則格網(wǎng)類似。</p><p>　　2.3 檢測(cè)規(guī)則分布數(shù)據(jù)中粗差的算法</p><p>　　2.3.1 基于坡度信息的粗差檢測(cè)</p><p>　　坡度是地表的固有屬性，在局部連續(xù)空間的漸變模型上，坡度變化也是連的，因此可采用采樣點(diǎn)與周圍點(diǎn)

21、的坡度變化是否一致來(lái)檢測(cè)是否含有粗差，通以局部3*3窗口對(duì)每一采樣點(diǎn)進(jìn)行判斷。其基本思想是對(duì)每個(gè)表面上的點(diǎn)，在坡度上，高程或突變量引起的形狀不連續(xù)，可能被懷疑有誤差，通過(guò)坡度上每個(gè)點(diǎn)，應(yīng)用坡度逼近或改變量來(lái)計(jì)算，考慮坡度變化的相對(duì)值，并以這些相對(duì)值計(jì)算一個(gè)統(tǒng)計(jì)值為判斷該點(diǎn)合法性的閩值，使計(jì)算結(jié)果更為可靠。如表1所示，P點(diǎn)在高程矩陣中的行列號(hào)為(I，J)，它的相臨8個(gè)點(diǎn)1，2，3，4，6，7，8，9的行列號(hào)分別為:1點(diǎn)(I+1，J-l)，

22、2點(diǎn)(I+1，J)，3點(diǎn)(I+，J+1)，4點(diǎn)(I，J-l)，6點(diǎn)(I，J+1)，7點(diǎn)(I-l，J-l)，8點(diǎn)(I-1，J)，9點(diǎn)(I-l，J+1)。 </p><p>　　表1 P點(diǎn)坡度計(jì)算</p><p>　　該算法分三步進(jìn)行，首先以檢測(cè)點(diǎn)P的8個(gè)鄰域點(diǎn)分別計(jì)算I，J方向的坡度值，然后計(jì)算各個(gè)方向的坡度變化值DSC，根據(jù)每個(gè)數(shù)據(jù)點(diǎn)在同一方向上的兩個(gè)DSC值相加，其值用于計(jì)算均方根差(

23、RMSE)，如果坡度變化一致的話，則同一點(diǎn)在同一方向上的兩個(gè)DSC值和的絕對(duì)值將是很小的值(接近0)，反之如果坡度變化不一致的話，這個(gè)值將比較大。若某一數(shù)據(jù)點(diǎn)在行列方向上的DSC值都大于閉值(閉值為RMSE的K倍)，則可確信它含有粗差。對(duì)于K值，不同情況下，可以使用不同的值，如果DSC值分布比較均勻(此時(shí)RMSE值比較小)，K可以取大一些的值;反之(RMSE值較大)，則K值以取小一些的值。對(duì)判斷含粗差的點(diǎn)進(jìn)行改正，以保證數(shù)據(jù)質(zhì)量的提高[

24、5]。</p><p>　　2.3.2 基于參數(shù)統(tǒng)計(jì)的粗差檢測(cè)</p><p>　　假設(shè)超限誤差只是局部相關(guān)的，F(xiàn)elicisimo算法研究的對(duì)象是某點(diǎn)高程值，和其鄰域點(diǎn)內(nèi)插出的高程值之，之間的差值民，，如果使用雙線性內(nèi)插而言，只要使用鄰域的四個(gè)點(diǎn)就夠了，在這種情況下，對(duì)矩陣第i行，第j列交叉處點(diǎn)高程值可按下式計(jì)算[6]:</p><p>　　高程估值與DEM中高程

25、值之差為：</p><p>　　如果此過(guò)程應(yīng)用于DEM中所有點(diǎn)，可以得到高程差值的代數(shù)均值和標(biāo)準(zhǔn)偏差，假設(shè)服從均值為，標(biāo)準(zhǔn)偏差為（均從這個(gè)采樣中獲得）的高斯分布 ，可以用雙尾檢驗(yàn)來(lái)驗(yàn)證是否屬于該分布的集合。通過(guò)引進(jìn)t統(tǒng)計(jì)量，，可以看做一個(gè)標(biāo)準(zhǔn)化殘差，因從模型中獲得的數(shù)據(jù)量很大，故可假設(shè)服從分布，對(duì)于置信水平ɑ=0.001，統(tǒng)計(jì)量的臨界值為3.219，由此可進(jìn)行兩個(gè)假設(shè)的檢驗(yàn)，其中零假設(shè)為，而被選假設(shè)為，任何殘差，

26、使的數(shù)據(jù)點(diǎn)都被懷疑為含有粗差，但事實(shí)上，大的值并不能指出粗差，而僅僅是一個(gè)警告符號(hào)。</p><p>　　表2 3*3窗口計(jì)算單元</p><p>　　粗差的修正應(yīng)該緊緊伴隨著粗差的探測(cè)過(guò)程，在DEM的柵格矩陣中，由于不允許空格點(diǎn)的存在，一旦某個(gè)粗差被檢測(cè)就應(yīng)該加以改正，這里最簡(jiǎn)單的方法就是用臨近點(diǎn)的高程均值來(lái)代替可疑點(diǎn)的高程值，需要指出的是任何可疑數(shù)據(jù)點(diǎn)都不應(yīng)參與高程估值的計(jì)算。上述過(guò)

27、程可以迭代進(jìn)行，其計(jì)算方法如表2 以3*3窗口為計(jì)算單元，每一次迭代中統(tǒng)計(jì)量和都會(huì)變化，直到?jīng)]有超過(guò)臨界值的殘差出現(xiàn)。</p><p>　　2.3.3 基于主成分分析的粗差檢測(cè)</p><p>　　主成分分析是把多個(gè)指標(biāo)化為少數(shù)幾個(gè)綜合指標(biāo)的一種統(tǒng)計(jì)方法[7]。在實(shí)際的研究中，為了全面分析問(wèn)題，往往使用眾多有關(guān)的變量。但是，變量太多不但會(huì)增加計(jì)算的復(fù)雜性，而且也給合理地分析問(wèn)題和解釋問(wèn)題帶

28、來(lái)困難。一般來(lái)說(shuō)，雖然每個(gè)變量都提供了一定的信息，但其重要性有所不同。實(shí)際上，在很多情況下，眾多變量間有一定的相關(guān)關(guān)系，人們希望利用這種相關(guān)性對(duì)這些變量加以“改造”，用為數(shù)較少的新變量來(lái)反映原變量所提供的大部分信息，通過(guò)對(duì)新變量的分析達(dá)到解決問(wèn)題的目的。主成分分析便是在這種降維的思想下產(chǎn)生的處理高維數(shù)據(jù)的統(tǒng)計(jì)方法。主成分分析的基本方法是通過(guò)構(gòu)造原變量的適當(dāng)?shù)木€性組合，以產(chǎn)生一系列互不相關(guān)的新變量，從中選出少數(shù)幾個(gè)新變量并使它們含有盡可能

29、多的原變量帶有的信息，從而使得用這幾個(gè)新變量代替原變量分析問(wèn)題和解決問(wèn)題成為可能。當(dāng)研究的問(wèn)題確定之后，變量中所含的“信息”的大小通常用該變量的方差來(lái)度量。主成分分析法進(jìn)行DEM粗差的探測(cè)包括以下幾個(gè)步驟:</p><p>　　(l)對(duì)于給定的大小為n*m的DEM，將其分別劃分為列主方向和行主方向度為w的帶狀區(qū)域，實(shí)施下列步驟直到滿足一定的準(zhǔn)則。</p><p>　　(2)處理列主方向的帶

30、狀區(qū)域。a.確定可能含有粗差的列;b.在每一列中確定差的位置，從而獲取備選粗差數(shù)組1。</p><p>　　(3)處理行主方向的帶狀區(qū)域。a.確定可能含有粗差的行;b.在每一行中確定差的位置，從而獲取備選粗差數(shù)組2。</p><p>　　(4)將兩套備選粗差進(jìn)行比較。</p><p>　　(5)提供計(jì)算準(zhǔn)則，改正所有的粗差。</p><p>

31、<b>　　(6)結(jié)束。</b></p><p>　　由于粗差數(shù)據(jù)影響統(tǒng)計(jì)量的計(jì)算，因而上述過(guò)程需要迭代進(jìn)行，在每一次代中，都可以得到一組“粗差候選序列”。如果確定為真正的粗差，就對(duì)相應(yīng)數(shù)據(jù)點(diǎn)加以修正，并從“候選序列”中去除，接著進(jìn)行下一步迭代過(guò)程，并設(shè)一定的條件來(lái)控制程序的進(jìn)程。</p><p>　　2.4 檢測(cè)不規(guī)則分布數(shù)據(jù)中粗差的算法</p>&l

32、t;p>　　規(guī)則格網(wǎng)DEM具有很多優(yōu)點(diǎn)，它的數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單，便于存儲(chǔ)和處理;但有不足之處，它對(duì)地表的描述沒(méi)有非規(guī)則DEM的精確，也無(wú)法顧及地形變化特征點(diǎn)和線。除此以外，在實(shí)際的生產(chǎn)中，規(guī)則格網(wǎng)DEM通常都是由非規(guī)則網(wǎng)DEM通過(guò)內(nèi)插而得到的。如果原始DEM為非規(guī)則數(shù)據(jù)，要轉(zhuǎn)成規(guī)則格網(wǎng)點(diǎn)在此過(guò)程中，原始數(shù)據(jù)中粗差點(diǎn)會(huì)影響到轉(zhuǎn)換后的多個(gè)格網(wǎng)點(diǎn)，從而增加格網(wǎng)粗差檢測(cè)的難度。因此，很有必要討論和驗(yàn)證基于非規(guī)則格DEM的粗差探算法。目前，最典型

33、的基于非規(guī)則DEM的粗差探測(cè)算法是李志林提出的點(diǎn)方算法[8]。</p><p>　　2.4.1 基于點(diǎn)方式的算法</p><p>　　首先，確定待定點(diǎn)P周圍的鄰域點(diǎn)范圍，然后計(jì)算窗口范圍內(nèi)所有點(diǎn)的平高程(或加權(quán)平均值)作為P點(diǎn)的估值，最后計(jì)算P點(diǎn)高程值與估值的高程差，如果高程差值大于閩值，則認(rèn)為P點(diǎn)含有誤差。</p><p><b>　　(l)鄰域點(diǎn)的范圍

34、</b></p><p>　　確定待定點(diǎn)P周圍的鄰域點(diǎn)范圍，可根據(jù)以P為中心的窗口指定，窗口的確定有三種方法，一種是定義窗口的尺寸，另一種是定義窗口覆蓋區(qū)域內(nèi)高程點(diǎn)的數(shù)量，還有一種方法是同時(shí)使用上述兩種方法來(lái)確定窗口的大小，通過(guò)計(jì)算區(qū)域內(nèi)點(diǎn)的數(shù)量和坐標(biāo)范圍確定一平均窗口。</p><p><b>　　(2)代表值的計(jì)算</b></p><

35、;p>　　在點(diǎn)方式算法中，把待測(cè)點(diǎn)鄰域點(diǎn)的平均高程作為該點(diǎn)的代表值。有兩種方法可計(jì)算鄰域點(diǎn)的平均高程，一種是簡(jiǎn)單的計(jì)算高程值的算術(shù)平均值，另一種是對(duì)每一個(gè)鄰域點(diǎn)賦以不同的權(quán)值。如果P的鄰域點(diǎn)包含粗差，用簡(jiǎn)單算術(shù)平均值法更加可信，計(jì)算速度也比較快。</p><p><b>　　(3)計(jì)算閉值</b></p><p>　　假設(shè)城是以第i個(gè)點(diǎn)為中心的鄰域點(diǎn)的算術(shù)平均值

36、，耳為從與第i個(gè)點(diǎn)的高程值的差值，即：</p><p>　　對(duì)DEM中所有的點(diǎn)，可以得到一系列的值，計(jì)算均值和標(biāo)準(zhǔn)偏差[9]：</p><p>　　其中，為由中心點(diǎn)的窗口范圍內(nèi)數(shù)據(jù)點(diǎn)計(jì)算的高程估值，為點(diǎn)的高程測(cè)量值。假設(shè)服從均值為，標(biāo)準(zhǔn)偏差為的正態(tài)分布，則檢測(cè)粗差的閾值為的倍（為常數(shù)）。閾值確定后，對(duì)DEM任一數(shù)據(jù)點(diǎn)，如果,則認(rèn)為含有粗差。一旦數(shù)據(jù)點(diǎn)被檢測(cè)超限，用其估值來(lái)代替可疑的高程值，

37、迭代進(jìn)行，直到?jīng)]有超過(guò)閾值的高程較差出現(xiàn)。</p><p>　　2.4.2 基于粗差簇群的算法</p><p>　　基于點(diǎn)方式的粗差檢測(cè)算法是針對(duì)檢測(cè)數(shù)據(jù)中僅存在單個(gè)粗差的情況，而事實(shí)上DEM中的粗差經(jīng)常以一種排列緊湊、數(shù)據(jù)巨大的簇群方式存在—這在自動(dòng)相關(guān)技術(shù)獲取的數(shù)據(jù)中經(jīng)常存在。與點(diǎn)方式算法相同，首先需要定義一以為中心的窗口，將窗口中的第一點(diǎn)從窗口中移去，從窗口中剩余的點(diǎn)計(jì)算新的“代表值

38、”即平均值，然后計(jì)算并記錄這個(gè)平均值與移去的數(shù)據(jù)點(diǎn)的值之差，此過(guò)程迭代進(jìn)行，直到窗口中所有的點(diǎn)都通過(guò)檢驗(yàn)[10]。</p><p>　　假設(shè)在窗口中有M個(gè)點(diǎn)，那么通過(guò)下式可計(jì)算M個(gè)差值：</p><p>　　式中，是窗口所有剩余數(shù)據(jù)點(diǎn)的平均值，P是窗口中所有數(shù)據(jù)點(diǎn)的平均值。余下的過(guò)程與點(diǎn)方式檢測(cè)粗差的算法相同，也是M個(gè)值用來(lái)計(jì)算一個(gè)統(tǒng)計(jì)值，并使用該統(tǒng)計(jì)值生成閩值，如果某一差值超過(guò)了這個(gè)閩值

39、，則認(rèn)為這個(gè)數(shù)據(jù)點(diǎn)含有粗差而需要將其進(jìn)行剔除。</p><p>　　2.5 基于等高線拓?fù)潢P(guān)系的粗差檢測(cè)與剔除</p><p>　　眾所周知，相鄰等高線的高程值之間的關(guān)系有且僅有三種:遞增、遞減或相等，根據(jù)這些關(guān)系，可對(duì)等高線的高程值是否有錯(cuò)作一判斷。在等高線地形圖上由于存在等高線密集、注記的壓蓋、斷崖地形等情況，常常造成等高線的不連續(xù)有時(shí)丟失的情形，因此不能僅僅依靠等高距來(lái)確定可疑處是否

40、錯(cuò)誤。在對(duì)所有的可疑處自動(dòng)檢測(cè)后，應(yīng)當(dāng)對(duì)每個(gè)可疑處根據(jù)等高線的關(guān)系由人工進(jìn)行校驗(yàn)并進(jìn)行修改，剔除粗差。</p><p>　　以上所提出的幾種探測(cè)方法各有千秋,而我們應(yīng)在實(shí)踐中根據(jù)自己的DEM的數(shù)據(jù)源、數(shù)據(jù)特點(diǎn)、現(xiàn)有的條件選擇合適的粗差檢測(cè)與剔除方法。</p><p>　　3 基于趨勢(shì)面的粗差剔除方法研究 </p><p>　　3.1 獲取原始DEM數(shù)據(jù)</p&

41、gt;<p>　　DEM的數(shù)據(jù)來(lái)源方式有地面測(cè)量（利用自動(dòng)記錄的測(cè)距經(jīng)緯儀在野外實(shí)測(cè)）、現(xiàn)有地圖數(shù)字化（利用數(shù)字化儀對(duì)已有地圖上的信息如等高線、地形線等進(jìn)行數(shù)字化，目前常用的數(shù)字化儀有手扶跟蹤數(shù)字化儀與掃描數(shù)字化儀）、空間傳感器（利用GPS等進(jìn)行數(shù)據(jù)采集）、數(shù)字?jǐn)z影測(cè)量方法（這是DEM數(shù)據(jù)采集最常用的一種方法）。DEM數(shù)據(jù)采集的方式有以上幾種方法，在實(shí)際生產(chǎn)中，究竟要采用哪種方法和生產(chǎn)工藝，主要是取決于DEM的分辨率、精度

42、要求、應(yīng)用范圍的大小以及成本和速度快慢等綜合考慮。現(xiàn)把常用的DEM數(shù)據(jù)采集方法以及各自特性的比較列于表3。</p><p>　　表3 DEM采集方法及各自特性的比較</p><p>　　3.2 DEM 數(shù)據(jù)預(yù)處理</p><p>　　DEM數(shù)據(jù)預(yù)處理是DEM 內(nèi)插之前的準(zhǔn)備工作，它是整個(gè)數(shù)據(jù)處理的一部分，一般包括數(shù)據(jù)格式的轉(zhuǎn)換、坐標(biāo)系統(tǒng)的變換、數(shù)據(jù)的編輯、柵格數(shù)據(jù)的

43、矢量化及數(shù)據(jù)分塊等內(nèi)容。</p><p>　　3.2.1 格式轉(zhuǎn)換</p><p>　　由于數(shù)據(jù)采集的硬、軟件系統(tǒng)各不相同，因而數(shù)據(jù)的格式可能也不相同。常用的代碼有ASCII碼、BCD碼及二進(jìn)制碼，每一次記錄的各項(xiàng)內(nèi)容及每項(xiàng)內(nèi)容的類型位數(shù)也可能各不相同，要根據(jù)相應(yīng)的DEM內(nèi)插軟件的要求，將各種數(shù)據(jù)轉(zhuǎn)換成該軟件所要求的數(shù)據(jù)格式[11]。</p><p>　　3.2.2

44、 坐標(biāo)變換</p><p>　　若采集的數(shù)據(jù)一般要轉(zhuǎn)換到地面坐標(biāo)系。地面坐標(biāo)系一般采用國(guó)家坐標(biāo)系，也可采用局部坐標(biāo)系。</p><p>　　3.2.3 數(shù)據(jù)編輯</p><p>　　將采集的數(shù)據(jù)用圖形方式顯示在計(jì)算機(jī)屏幕上（或展繪在數(shù)控測(cè)圖儀上），作業(yè)人員采用圖形交互方式進(jìn)行數(shù)據(jù)編輯。包括剔除錯(cuò)誤的、過(guò)密的與重復(fù)的點(diǎn)，發(fā)現(xiàn)某些需要補(bǔ)測(cè)，對(duì)斷面掃面數(shù)據(jù)，還要進(jìn)行掃描的

45、系統(tǒng)誤差的糾正。</p><p>　　3.2.4 柵格數(shù)據(jù)轉(zhuǎn)換為矢量數(shù)據(jù)</p><p>　　由地圖掃描數(shù)字化獲取的掃描影像是一灰度陣列，首先經(jīng)過(guò)二值化處理，再經(jīng)過(guò)濾波或形態(tài)處理（利用數(shù)學(xué)形態(tài)進(jìn)行各種運(yùn)算），并進(jìn)行邊緣跟蹤，獲得等高線上按順序排列的點(diǎn)坐標(biāo)，即矢量數(shù)據(jù)。</p><p>　　3.2.5 數(shù)據(jù)分塊</p><p>　　由于數(shù)據(jù)采集

46、方式不同，數(shù)據(jù)的排列順序也不同，例如等高線是按各條等高線采集的先后順序排列的。但在內(nèi)插DEM時(shí)，待定點(diǎn)常常只與其周圍的數(shù)據(jù)點(diǎn)有關(guān)，為了能在大量的數(shù)據(jù)點(diǎn)中迅速地找的所需要的數(shù)據(jù)點(diǎn)，必須將其進(jìn)行分塊。在某些軟件中，需要將數(shù)據(jù)點(diǎn)劃分成計(jì)算單元，每個(gè)計(jì)算單元之間有一定的重疊度，以保證單元的連續(xù)性。分塊的方法是先將整個(gè)區(qū)域分成等間隔得格網(wǎng)（通常比DEM格網(wǎng)大），然后將數(shù)據(jù)點(diǎn)按格網(wǎng)分成不同的類，通常有交換法和鏈指針?lè)ā?lt;/p><

47、;p>　　3.2.6 子邊界的提取</p><p>　　根據(jù)離散的數(shù)據(jù)點(diǎn)內(nèi)插規(guī)則格網(wǎng)DEM，通常是將地面看做是一個(gè)光滑的連續(xù)曲面，但是地面上存在著各種各樣的斷裂帶，如陡崖、絕壁以及各種人工地物，如路堤等，使地面并不光滑，這就需要將地面分成若干區(qū)域，即子區(qū)，是每一個(gè)子區(qū)的表面呈一連續(xù)光由特征線滑曲面。這些子區(qū)的邊界由特征線（如斷裂線）與區(qū)域的邊界線組成。確定每一子區(qū)的邊界可以采用專門的數(shù)據(jù)結(jié)構(gòu)或利用圖論等多種

48、方法來(lái)解決。</p><p>　　3.3 基于趨勢(shì)面的粗差檢測(cè)與剔除方法</p><p>　　基于趨勢(shì)面的粗差探測(cè)基礎(chǔ)是地形所具有的自相關(guān)性，即地形變化符合一定的自然趨勢(shì)，表現(xiàn)為連續(xù)空間的漸變模型，并且這種連續(xù)變化可以用一種平滑的數(shù)學(xué)曲面—趨勢(shì)面來(lái)加以描述。對(duì)粗差的檢測(cè)，可以通過(guò)模型誤差即實(shí)際觀測(cè)值與趨勢(shì)面計(jì)算值(模型值)之差來(lái)判斷其是否屬于異常數(shù)據(jù)，因此趨勢(shì)面分析的一個(gè)典型應(yīng)用就是揭示研

49、究區(qū)域中不同于總趨勢(shì)的最大偏離部分。因此，可以采用趨勢(shì)面分析找出偏離趨勢(shì)超過(guò)一定閩值的異常數(shù)據(jù)可疑點(diǎn)。應(yīng)用趨勢(shì)面分析進(jìn)行粗差探測(cè)需要注意兩點(diǎn)，一是趨勢(shì)面函數(shù)的確定，另一個(gè)是閾值問(wèn)題，即當(dāng)觀測(cè)值和計(jì)算值相差多大時(shí)，該點(diǎn)可被懷疑為粗差點(diǎn)，理論上任何復(fù)雜的曲面都可用高階多項(xiàng)式去逼近，但用高階多項(xiàng)式的解本身不穩(wěn)定，同時(shí)多項(xiàng)式系數(shù)的物理意義也不清楚，可能會(huì)導(dǎo)致不符合實(shí)際地形起伏的趨勢(shì)。關(guān)于閾值問(wèn)題，通常按統(tǒng)計(jì)方法，假設(shè)誤差頻率呈正態(tài)分布，取三倍的

50、中誤差為極限值，但畢竟設(shè)定閩值的人為因素很大，且現(xiàn)實(shí)世界中的變化也可導(dǎo)致不正確的判斷，造成粗差點(diǎn)的遺漏或錯(cuò)選。趨勢(shì)面可有各種不同的形式，在現(xiàn)實(shí)生活中，我們通常選擇由下式構(gòu)成的最小二乘趨勢(shì)面:</p><p>　　據(jù)處理區(qū)域的形味大小，可以靈活地選擇不同階次的多項(xiàng)式，對(duì)大而復(fù)雜的區(qū)域采用較高階次。根據(jù)統(tǒng)計(jì)規(guī)律，常用2倍或3倍中誤差作為極限誤差，即模型誤差大于極限誤差的觀測(cè)數(shù)據(jù)被認(rèn)為是粗差。趨勢(shì)面分析技術(shù)的特點(diǎn)是能將

51、問(wèn)題簡(jiǎn)單化、局部化，能找出大部分的可疑數(shù)。</p><p><b>　　3.4 實(shí)驗(yàn)分析 </b></p><p>　　根據(jù)提供的 10 個(gè)數(shù)據(jù)點(diǎn)的坐標(biāo)和待求點(diǎn)的平面坐標(biāo)，利用移動(dòng)二次曲面擬合法，由格網(wǎng)點(diǎn)周圍的 10 個(gè)已知點(diǎn)內(nèi)插出待求格網(wǎng)點(diǎn)的高程，編制相應(yīng)的程序并進(jìn)行調(diào)試，最后解算出格網(wǎng)點(diǎn)P的高程并提交源程序代碼[12]。</p><p>&

52、lt;b>　　操作過(guò)程：</b></p><p>　　1、首先要導(dǎo)入三個(gè)程序矩陣模塊"stdafx.h"</p><p>　　"SMatrix.h" "SingleImageResection.h"</p><p>　　2、輸入已知點(diǎn)數(shù)據(jù)的坐標(biāo)</p><p>　

53、　3、選擇二次曲面作為擬合曲面：</p><p><b>　　列出誤差方程：</b></p><p>　　4、組方程，求出6個(gè)系數(shù)：</p><p>　　5、計(jì)算出待求點(diǎn)高程。</p><p>　　表4 已知數(shù)據(jù)點(diǎn)坐標(biāo)</p><p>　　編程計(jì)算點(diǎn)（110,110）上的高程。</p>

54、;<p>　　運(yùn)用C++進(jìn)行DEM內(nèi)插法的編程程序的代碼:</p><p>　　#include "stdafx.h"</p><p>　　#include "SMatrix.h"</p><p>　　#include "SingleImageResection.h"</p>&

55、lt;p>　　int main(int argc, char* argv[])</p><p><b>　　{</b></p><p>　　SMatrix X(10,1); 系數(shù)矩陣</p><p>　　SMatrix Y(10,1); 系數(shù)矩陣</p><p>　　SMatrix Z(10,1); 系

56、數(shù)矩陣</p><p>　　X[0][0] = 102; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[1][0] = 109; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[2][0] = 105; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[3][0] = 103; // 輸

57、入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[4][0] = 108; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[5][0] = 105; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[6][0] = 115; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[7][0] = 118;

58、 // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[8][0] = 116; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　X[9][0] = 113; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[0][0] = 110; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[1][0] =

59、113; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[2][0] = 115; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[3][0] = 103; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[4][0] = 105; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[5][0

60、] = 108; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[6][0] = 104; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[7][0] = 108; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[8][0] = 113; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Y[

61、9][0] = 118; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[0][0] = 15; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[1][0] = 18; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[2][0] = 19; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>

62、　　Z[3][0] = 17; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[4][0] = 21; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[5][0] = 15; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[6][0] = 20; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p

63、>　　Z[7][0] = 15; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[8][0] = 17; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　Z[9][0] = 22; // 輸入已知數(shù)據(jù)點(diǎn)坐標(biāo) //</p><p>　　for(int i = 0 ;i <10;i++)</p><p&

64、gt;<b>　　{</b></p><p>　　X[i][0] = X[i][0] -110;</p><p>　　Y[i][0] = Y[i][0] -110;</p><p><b>　　}</b></p><p>　　SMatrix M(10,6);</p><p&g

65、t;　　for(i = 0 ;i <10;i++)</p><p><b>　　{</b></p><p>　　M[i][0] =X[i][0] *X[i][0] ;</p><p>　　M[i][1] =X[i][0] *Y[i][0] ;</p><p>　　M[i][2] =Y[i][0] *Y[i][0]

66、;</p><p>　　M[i][3] =X[i][0];</p><p>　　M[i][4] =Y[i][0] ;</p><p>　　M[i][5] = 1 ;</p><p><b>　　}</b></p><p>　　SMatrix P(10,10);</p><p&g

67、t;　　for(i = 0 ;i <10;i++)</p><p><b>　　{</b></p><p>　　P[i][i] = 1/(X[i][0] *X[i][0]+Y[i][0] *Y[i][0]);</p><p><b>　　}</b></p><p>　　SMatrix x(6

68、,1);</p><p>　　x = (M.T() * P * M).Invert() * M.T() * P * Z;</p><p>　　printf("待定點(diǎn)的高程是：%f\n",x[5][0]);</p><p><b>　　return 0;</b></p><p><b>　　}

69、</b></p><p><b>　　運(yùn)行結(jié)果截圖</b></p><p>　　按上述方法計(jì)算模型誤差即實(shí)際觀測(cè)值與趨勢(shì)面計(jì)算值(模型值)之差來(lái)判斷其是否屬于異常數(shù)據(jù)。采用趨勢(shì)面分析找出偏離趨勢(shì)超過(guò)一定閾值的異常數(shù)據(jù)可疑點(diǎn)。粗差的修正應(yīng)該緊緊伴隨著粗差的探測(cè)過(guò)程，在DEM的柵格矩陣中，由于不允許空格點(diǎn)的存在，一旦某個(gè)粗差被檢測(cè)就應(yīng)該加以改正，這里最簡(jiǎn)單的方

70、法就是用臨近點(diǎn)的高程均值來(lái)代替可疑點(diǎn)的高程值，需要指出的是任何可疑數(shù)據(jù)點(diǎn)都不應(yīng)參與高程估值的計(jì)算。對(duì)于不同的數(shù)據(jù)窗口，分別利用最小二乘法進(jìn)行移動(dòng)曲面擬合，得到各個(gè)數(shù)據(jù)窗口內(nèi)高程殘差值列向量，以此為基礎(chǔ)計(jì)算驗(yàn)后殘差值方差式，并以k值作為檢測(cè)粗差的閉值，由于同一個(gè)DEM高程點(diǎn)在不同的擬合區(qū)域中計(jì)算出不同的殘差值，而這些殘差值對(duì)同一DEM高程點(diǎn)是否為粗差的判斷并不一致，即在一個(gè)擬合面內(nèi)認(rèn)為該點(diǎn)是粗差，而在另一個(gè)擬合面內(nèi)認(rèn)為是正常點(diǎn)，因此統(tǒng)計(jì)每

71、個(gè)點(diǎn)在相關(guān)的擬合面內(nèi)被懷疑為粗差的百分比（percent），當(dāng)percent大于某一限值，則認(rèn)為點(diǎn)i含有粗差。中位數(shù)算法是以殘差的中值作為檢驗(yàn)統(tǒng)計(jì)量來(lái)檢測(cè)粗差，均值算法是以殘差的平均值作為檢驗(yàn)統(tǒng)計(jì)量來(lái)檢測(cè)粗差，雖然它們都是對(duì)同一問(wèn)題的不同解決方式，但中位數(shù)算法比均值算法更為簡(jiǎn)單。由</p><p><b>　　4 總結(jié)與展望</b></p><p>　　認(rèn)真總結(jié)所做的

72、工作以及在論文撰寫中的遇到的困難，認(rèn)為在DEM的粗差探測(cè)中仍然有一些問(wèn)題有待進(jìn)一步研究與解決:本文在數(shù)據(jù)處理方面較為不足，如果能夠?qū)B可視化編程語(yǔ)言和MATLAB工具相結(jié)合，充分利用MATLAB的運(yùn)算功能和VB開(kāi)發(fā)界面方便的特點(diǎn)進(jìn)行混合編程，即用VB設(shè)計(jì)界面作為主程序，調(diào)用MATLAB子程序，數(shù)據(jù)處理將會(huì)簡(jiǎn)單、方便。今天，DEM已作為一個(gè)獨(dú)立的產(chǎn)品而存在，并越來(lái)越廣泛地被用來(lái)代替?zhèn)鹘y(tǒng)地圖中等高線對(duì)地形的描述，成為地理信息系統(tǒng)的核心數(shù)據(jù)

73、庫(kù)以及地學(xué)分析的基礎(chǔ)數(shù)據(jù)。目前DEM粗差探測(cè)領(lǐng)域中，針對(duì)于規(guī)則格網(wǎng)DEM算法研究成果較為豐富，這與規(guī)則格網(wǎng)形式的DEM具有簡(jiǎn)單的數(shù)據(jù)結(jié)構(gòu)和便利的存儲(chǔ)方式有很大的關(guān)系。相反，由于非規(guī)則形式的DEM數(shù)據(jù)結(jié)構(gòu)相對(duì)復(fù)雜，這方面的算法設(shè)計(jì)并不很多，較為實(shí)用的是基于點(diǎn)方式的算法。點(diǎn)方式算法是通過(guò)中心點(diǎn)的內(nèi)插值計(jì)算檢驗(yàn)量，然后進(jìn)行粗差的判釋。無(wú)論是規(guī)則格網(wǎng)DEM算法還是非規(guī)則格網(wǎng)DEM算法來(lái)剔除粗差，這些方法從理論上講都是將粗差歸入函數(shù)模型來(lái)實(shí)現(xiàn)粗差

74、的探測(cè)。在現(xiàn)代測(cè)量數(shù)據(jù)處理中如何消除粗差的影響，特別是自動(dòng)化觀測(cè)水平高，數(shù)據(jù)量大的情況下，就顯得越來(lái)越重要</p><p><b>　　參 考 文 獻(xiàn)</b></p><p>　　[1] 柯正誼，何建邦，池天河.數(shù)字地面模型[M」.中國(guó)科學(xué)技術(shù)出版社，1993：103-215.</p><p>　　[2] 李志林,朱慶.數(shù)字高程模型[M].武漢

75、:武漢大學(xué)出版社, 2001：77-83.</p><p>　　[3] 陳再輝，路曉峰.基于自適應(yīng)抗差最小二乘的DEM數(shù)據(jù)粗差剔除[J].海洋測(cè)繪，2006，26(6):15-17.</p><p>　　[4] ?？?，易思蓉，韓春華.基于坡度RMSE與3維可視化的格網(wǎng)DEM粗差檢測(cè)與剔除[J].測(cè)繪通報(bào)，2007，6(9):23-25.</p><p>　　[5]

76、韓玲.格網(wǎng)DEM中基于坡度信息的粗差檢測(cè)與剔除方法實(shí)驗(yàn)[J].長(zhǎng)安大學(xué)學(xué)報(bào)，2003，25(l):74-78.</p><p>　　[6] 黃宏波，梁鑫，楊曉云，等.基于參數(shù)統(tǒng)計(jì)的DEM粗差探測(cè)算法[J].測(cè)繪工程2008，8（2）：12-15.</p><p>　　[7] 王貴滿,王東華.DEM的粗差探測(cè)與生產(chǎn)質(zhì)量控制方法[J].測(cè)繪與空間地理信息.2011，2（1）：1-7.</

77、p><p>　　[8] 楊曉云,顧利亞,岑敏儀.基于不同大小窗口的移動(dòng)曲面擬合法探測(cè)不規(guī)則DEM粗差的一種方法[J].測(cè)繪學(xué)報(bào),2005，34 (2):14-15.</p><p>　　[9] 楊元喜.自適應(yīng)抗差最小二乘估計(jì)[J].測(cè)繪學(xué)報(bào),1996,25(3): 206-211.</p><p>　　[10] 黃幼才.數(shù)據(jù)探測(cè)與抗差估計(jì)[M].北京:測(cè)繪出版社,199

78、0：180-185.</p><p>　　[11] 徐士良.數(shù)值方法與計(jì)算機(jī)實(shí)現(xiàn)[M].北京:清華大學(xué)出版社,2006：33-56.</p><p>　　[12] 傅祖蕓,趙梅娜,丁巖,等,譯. C語(yǔ)言數(shù)值算法程序大全[M].北京:電子工業(yè)出版社,1995：98-115.</p><p>　　DEM Gross Error Picking Method Resear

79、ch</p><p>　　Abstract: Digital elevation model (DEM) as one of 4 D products (DEM, DOM - numbers are projective like, DLG-digital line drawing, DRG - digital raster graph), is a kind of digital description and

80、 simulation of the surface of the earth. The earth is an important part of spatial data infrastructure , it can be used to establish a terrain database and all kinds of GIS (geographic information system), the basis of q