【正文】 37 ?椀最甀爀攀...吀梔攀挀漀渀琀攀渀琀猀漀昀渀攀眀昀攀愀琀甀爀攀攙愀琀愀猀攀琀.椀最甀爀攀....椀最椀琀愀氀挀漀渀琀漀甀爀氀椀渀攀猀.椀最甀爀攀...吀梔攀爀攀猀甀氀琀猀漀昀攙椀最椀琀愀氀挀漀渀琀漀甀爀氀椀渀攀猀.椀最甀爀攀...吀梔攀愀琀琀爀椀戀甀琀攀猀漀昀挀漀渀琀漀甀爀氀椀渀攀猀.椀最甀爀攀....。 感謝我的父母和家人對(duì)我多年來在物質(zhì)及精神上的幫助和支持，是他們的不斷鼓勵(lì)使得我能夠順利完成學(xué)業(yè)。從張老師身上，我看到了一絲不茍的治學(xué)態(tài)度和兢兢業(yè)業(yè)的科學(xué)研究精神，這些都將影響我今后的學(xué)習(xí)和工作，指導(dǎo)我在今后的道路上不斷前進(jìn)。岳老師淵博的知識(shí)、嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度、在學(xué)術(shù)上積極創(chuàng)新的開拓精神和在事業(yè)上銳意進(jìn)取的執(zhí)著精神都令我欽佩。張加龍，男，1981年生，湖北襄樊人，碩士，助教，西南林學(xué)院森林經(jīng)理教研室教師；主要擔(dān)任課程為地理信息系統(tǒng)原理、 GIS集成與開發(fā)、專業(yè)英語等課程，研究方向?yàn)?3S系統(tǒng)集成與開發(fā)、資源信息系統(tǒng)管理與應(yīng)用。[1]沈小樂,郭志佳, [J].測繪與空間地理信息,2007,02(30):140~141 參考文獻(xiàn) [2]張凱選,潘夢清 , DEM方法的研究 [J].測繪工程 ,2007,03(16):16 [3] DEM [C].武漢大學(xué) ,2005,4:5~7 [4]王曉東. DEM精度檢查與質(zhì)量評(píng)估系統(tǒng)研究 [C].吉林大學(xué) ,2006,6:2 [5]劉學(xué)軍, Delaunay三角剖分與修改算法 [J].測繪學(xué)報(bào) ,2001,02(10):82~87 [6]江帆. DEM表面建模與精度評(píng)估方法研究 [C].中國人民解放軍信息工程大學(xué) ,2006,7:7 [7]胡鵬,吳艷蘭 , [J].地理信息科學(xué) ,2003,03:64~69 [8]湯國安, [M].科學(xué)出版社 ,2001,7:121~122 [9]湯國安,劉學(xué)軍, [M].北京:科學(xué)出版社,2005,8:18~23 [10]吳立新, [M].北京:科學(xué)出版社 ,2003 [11]黃力,邱滌珊 ,王峰. GIS中 DEM的表示方法[J].計(jì)算機(jī)與信息技術(shù) ,1999 [12] [M].北京:科學(xué)出版社 ,2001 [13] GIS的若干理論與技術(shù)[M].武漢測繪科技大學(xué)出版社 ,1999 [14][M].北京:北京高等教育出版社 ,2001 [15]鄔倫,劉瑜,——原理、方法和應(yīng)用[ M].北京:科學(xué)出版社 ,2002 [16][ M].解放軍出版社,2001 [17]邵春麗,胡鵬,[J].測繪科學(xué). 2004,29（6）:68~71 [18][C].河海大學(xué) ,2007,6:3 [19]湯國安,[M].北京:科學(xué)出版社, 2006 [20]楊曉云 , (DEM)的質(zhì)量控制及精度分析 [J].廣西工學(xué)院學(xué)報(bào),2006,02(15):37~38 [21] DEM的精度評(píng)價(jià)方法探討 [J].遙感信息,2006,02(13):16~28 [22] [J].測繪通報(bào) ,2002,02(3):18~20 [23]李志林, [M].武漢大學(xué)出版社 ,2001 [24] 1:10000,1:50000數(shù)字高程模型（CH/T10082001） [25]Zhu C Q, SHI W Z etal. The impact of a random error at the TIN nodes on the accuracy of the DEM under linear interpolation[J].International Journal of Remote Sensing， 2005,26(24):5493~5508 [26]Tang Guoan. A Research on the accuracy of digital elevation models[M].Beijing: Science Press,2000 指導(dǎo)教師簡介岳彩榮，男，漢族， 1964年 12月生于云南建水，博士，教授，碩士研究生導(dǎo)師，現(xiàn)任西南林學(xué)院資源學(xué)院森林經(jīng)理學(xué)教研室主任。因分析建立 DEM整體精度理論模型。 （3）實(shí)際地形地貌復(fù)雜多變，如何建立適應(yīng)各種地貌形態(tài)特征的高精度的 DEM，更逼真的數(shù)字再現(xiàn)實(shí)際地貌形態(tài)。 （1）地形特征線的提取工作是手動(dòng)提取的，應(yīng)該研究算法實(shí)現(xiàn)自動(dòng)提取，盡量減少人工干預(yù)。 論文需要深入的地方 DEM的研究是當(dāng)前計(jì)算機(jī)圖形學(xué)、地理信息系統(tǒng)、虛擬現(xiàn)實(shí)等領(lǐng)域的熱點(diǎn)內(nèi)容之一，特別是在 GIS中，數(shù)字高程的精度研究更有著特殊的意義。 5 VB開發(fā)了一個(gè)自動(dòng)生成 TIN的程序，從而進(jìn)一步理解了 TIN的生成方式。解決三角形的邊穿越等高線的方法是把等高線當(dāng)作特征線處理，構(gòu)成約束 TIN，本論文也用到了該方法。三角形的邊穿越等高線是指滿足 DT三角剖分準(zhǔn)則的三角形與等高線相交。本論文特征點(diǎn)的提取都是手工的，這種方法雖然簡單，但工作量大，耗時(shí)多。 3 在查閱文獻(xiàn)時(shí)，發(fā)現(xiàn)是有方法可以自動(dòng)提取等高線上特征點(diǎn)。 本論文完成的工作有： 1 實(shí)驗(yàn)區(qū)大部分地勢比較陡峭，等高線密集，通過采用先數(shù)字化等高線的計(jì)曲線后數(shù)字化等高線間曲線的方法，減少粗差的產(chǎn)生。如果精度得不到保證， DEM信息將變得毫無意義。對(duì) DEM數(shù)據(jù)質(zhì)量的科學(xué)度量及精確表達(dá)，己成為 DEM數(shù)據(jù)的生產(chǎn)者和使用者非常關(guān)注的問題。逼近誤差，屬于函數(shù)誤差它更能反應(yīng)誤差的來源。如果精度評(píng)定的結(jié)果達(dá)不到上述標(biāo)準(zhǔn)，說明數(shù)字化等高線有誤，要重新數(shù)字化等高線。 假設(shè)檢查點(diǎn)的高程為 Z k (k=1，2，…，n)，在建立 DEM之后，由 DEM內(nèi)插出這些點(diǎn)的高程為 R k，則 DEM的精度為： σDEM = n 1 ∑ ? n k ZRk( k)2 （4）表 2 1:1萬 DEM精度標(biāo)準(zhǔn)（中國國家測繪局 1998） [24] Table2 1:1 million DEM accuracy standards (Chinese State Bureau of Surveying and Mapping 1998) 地形類別地形圖基本等高距（m）地面坡度（度） DEM格網(wǎng)間距（m）格網(wǎng)點(diǎn)高程中誤差（m）一級(jí)二級(jí)三級(jí)平地 1 ≤2 丘陵地 2~6 山地 5 6~25 高山地 10 ≥25 根據(jù) TIN內(nèi)插出 DEM，可得知這些點(diǎn)的內(nèi)插高程值。 （2）用檢查點(diǎn)法做定量評(píng)價(jià)檢查點(diǎn)法：即事先將檢查點(diǎn)按格網(wǎng)或任意形式進(jìn)行分布，對(duì)生成的 DEM在這些點(diǎn)處進(jìn)行檢查。結(jié)果如下所示： 圖 (紅色代表回放等高線，紫色代表原始等高線) Contour lines intervals由上圖可見，通過 DEM內(nèi)插生成的等高線和原始等高線在等高線密集區(qū)域，擬合程度比較好，在等高線稀疏區(qū)域，擬合程度比較差。因此，需要對(duì)研究區(qū)域進(jìn)行裁剪，運(yùn)用 ArcToolboxAnalysist ToolsExtractClip命令，新建點(diǎn)層（points），選擇進(jìn)行精度評(píng)價(jià)需要的特征點(diǎn)，結(jié)果如下所示： Figure Accuracy of regional plans 精度評(píng)定區(qū)域圖 精度評(píng)定方法 本論文選擇兩種方法進(jìn)行評(píng)價(jià)，一種用等高線回放法做定性評(píng)價(jià)，一種用檢查點(diǎn)法做定量評(píng)價(jià)。由于影響 DEM精度的多樣性，在考察 DEM的精度時(shí)，不僅要 DEM的單點(diǎn)誤差，還要考慮 DEM在山區(qū)、平緩地區(qū)和破碎地區(qū)的整體形狀，使 DEM不僅在單點(diǎn)的精度達(dá)到相當(dāng)?shù)乃?，而且整個(gè) DEM的形狀和實(shí)際地形保持一致。對(duì) DEM的內(nèi)插模型的檢查則要復(fù)雜一些，從數(shù)學(xué)的角度而言，可從逼近程度、外推能力、平滑效果、惟一性、計(jì)算時(shí)間等方面進(jìn)行比較檢查和評(píng)價(jià)，但在實(shí)際應(yīng)用中，無法對(duì)內(nèi)插模型的這些特性進(jìn)行檢查，更為主要的是，大量的實(shí)踐表明，影響 DEM精度的主要因素取決于原始數(shù)據(jù)的精度和內(nèi)插模型顧及地形特征與否，而與內(nèi)插的其它因素并無明顯的關(guān)系。對(duì) DEM原始數(shù)據(jù)的精度進(jìn)行檢查的實(shí)質(zhì)是檢查數(shù)據(jù)中是否含有誤差(包括系統(tǒng)誤差、偶然誤差和粗差)。 影響 DEM精度的主要因素 （1）地形表面的特征 （2）DEM原始數(shù)據(jù)的分布和密度 （3）DEM表面建模方法 （4）DEM表面的自身特性 （5）矢量數(shù)據(jù)狀況 DEM精度檢查的內(nèi)容 DEM的精度檢查主要包括以下內(nèi)容[21]： 1 DEM原始的數(shù)學(xué)基礎(chǔ) 2 DEM數(shù)據(jù)起止點(diǎn)坐標(biāo)的正確性 3 DEM原始數(shù)據(jù)的精度 4 DEM的高程值有效范圍區(qū)是否正確 5 DEM的內(nèi)插模型 6 DEM產(chǎn)品的精度 7 DEM的元數(shù)據(jù)文件是否正確 這些內(nèi)容中，對(duì)于 DEM原始的數(shù)學(xué)基礎(chǔ)、DEM數(shù)據(jù)起止點(diǎn)坐標(biāo)的正確性、 DEM的高程值有效范圍區(qū)的正確性、DEM的元數(shù)據(jù)文件的正確性等問題的檢查一般都比較容易，并非特別困難。 DEM誤差源主要有兩種：一是對(duì)實(shí)際地形表面采樣所引起的誤差，即原始數(shù)據(jù)的誤差；另一種是數(shù)據(jù)重采樣引起的誤差，即由 DEM內(nèi)插算法引起的誤差。 實(shí)驗(yàn)區(qū) DEM圖 DEM of experimental area 5質(zhì)量檢查與精度評(píng)定 DEM誤差及誤差源 DEM是通過對(duì)自然界中空間實(shí)體的各種直接或間接觀測得到的。（5）TIN轉(zhuǎn)換成規(guī)則格網(wǎng) DEM選擇 3DAnalysistConversion—TIN to Raster生成格網(wǎng)間距為 * DEM。若在建立 TIN時(shí)考慮山脊點(diǎn)，山谷點(diǎn)，鞍點(diǎn)這些特征點(diǎn)，出現(xiàn)平三角形現(xiàn)象會(huì)減少，從而建立高質(zhì)量的TIN（（2））。本論文實(shí)驗(yàn)樣區(qū)為宜良縣某村委會(huì)地形圖數(shù)據(jù)。轄五鎮(zhèn)七鄉(xiāng)，總?cè)丝?38萬。 Chr(10) Loop Close 3 ②判斷新加入點(diǎn)是否在已生成三角形的外接圓內(nèi) eps = InCircle = False If Abs(y1 y2) eps And Abs(y2 y3) eps Then MsgBox INCIRCUM F Points are coincident !! Exit Function End If If Abs(y2 y1) eps Then m2 = (x3 x2) / (y3 y2) mx2 = (x2 + x3) / 2 my2 = (y2 + y3) / 2 xc = (x2 + x1) / 2 yc = m2 * (xc mx2) + my2 ElseIf Abs(y3 y2) eps Then m1 = (x2 x1) / (y2 y1) mx1 = (x1 + x2) / 2 my1 = (y1 + y2) / 2 xc = (x3 + x2) / 2 yc = m1 * (xc mx1) + my1 Else m1 = (x2 x1) / (y2 y1) m2 = (x3 x2) / (y3 y2) mx1 = (x1 + x2) / 2 mx2 = (x2 + x3) / 2 my1 = (y1 + y2) / 2 my2 = (y2 + y3) / 2 xc = (m1 * mx1 m2 * mx2 + my2 my1) / (m1 m2) yc = m1 * (xc mx1) + my1 End If dx = x2 xc dy = y2 yc rsqr = dx * dx + dy * dy r = Sqr(rsqr) dx = xp xc dy = yp