4等高線的自動(dòng)追蹤

1、等高線自動(dòng)追蹤等高線自動(dòng)追蹤DEM的表示的表示格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤TIN等高線追蹤等高線追蹤TIN和格網(wǎng)模型的應(yīng)用比較和格網(wǎng)模型的應(yīng)用比較數(shù)字地形表達(dá)20世紀(jì)中葉后，隨著計(jì)算機(jī)科學(xué)、現(xiàn)代數(shù)學(xué)和計(jì)算機(jī)圖形學(xué)的發(fā)展，各種數(shù)字地形表達(dá)方式得到迅猛的發(fā)展。計(jì)算機(jī)和計(jì)算機(jī)技術(shù)在測(cè)繪方面的應(yīng)用使得測(cè)繪學(xué)科逐步向數(shù)字化與自動(dòng)化、實(shí)時(shí)處理與多用途的方向發(fā)展。計(jì)算機(jī)技術(shù)在很大程度上改變了地圖制圖的生產(chǎn)方式，同時(shí)也改變著地圖產(chǎn)品的樣式和用圖概念。數(shù)字地形表達(dá)的方式可以分為兩大類，即數(shù)學(xué)描述和圖像描述。使用傅立葉級(jí)數(shù)和多項(xiàng)式來(lái)描述地形是常用的數(shù)學(xué)描述方法。規(guī)則格網(wǎng)、不規(guī)則格網(wǎng)、等高線則是圖像描述的常用方

2、式。2.1 2.1 DEM的表示的表示2.1 2.1 DEM的表示的表示數(shù)字地形表達(dá)測(cè)繪學(xué)中數(shù)字地形數(shù)字地形(地面地面)模型模型（DTM：Digital Terrain Model）是地形數(shù)字化的表達(dá)方式。數(shù)字地形模型是新一代的地形圖，地貌和地物不再用直觀的等高線和圖例符號(hào)在紙上表達(dá)，而是通過(guò)計(jì)算機(jī)儲(chǔ)存介質(zhì)中的大量密集的地面點(diǎn)的空間坐標(biāo)和地形屬性編碼，以數(shù)字的形式描述。1.通常根據(jù)不同的具體需要，將某些地面的特性信息與地形信息結(jié)合在一起，構(gòu)成有特殊用途的數(shù)字地面模型。 數(shù)字地面模型中所包含的地面特性信息類型一般可分為下列四組： 1地貌信息：高程、坡度、坡向、坡面形態(tài)及描述地表起伏情況的更為復(fù)雜

3、的地貌因子； 2基本地物信息：水系、交通網(wǎng)、居民點(diǎn)和工礦企業(yè)及境界線； 3主要的自然資源和環(huán)境信息：土壤、植被、地質(zhì)、氣候； 4主要的社會(huì)經(jīng)濟(jì)信息：人口、工農(nóng)業(yè)產(chǎn)值、經(jīng)濟(jì)活動(dòng)等。2.1 2.1 DEM的表示的表示數(shù)字高程模型-DEM傳統(tǒng)地形測(cè)圖中，地物用地物符號(hào)表達(dá)，地貌用高程點(diǎn)或等高線表達(dá)。從地形測(cè)繪的角度出發(fā)，一般僅把基本地形圖中的地形要素、特別是高程信息高程信息，作為數(shù)字地形模型數(shù)字地形模型（ DTM ）的內(nèi)容。即通常DTM的主要內(nèi)容是數(shù)字高程模型（DEM: Digital Elevation Model ）。從測(cè)繪學(xué)的角度，地面高程模型是高程Z關(guān)于平面坐標(biāo)X，Y兩個(gè)自變量的連續(xù)函數(shù)，D

4、EM只是地面高程模型的一個(gè)有限的離散表示。2.1 2.1 DEM的表示的表示數(shù)字高程模型-DEM表示法1.數(shù)學(xué)方法：用數(shù)學(xué)方法來(lái)表達(dá)，可以采用整體擬合方法，即根據(jù)區(qū)域所有的高程點(diǎn)數(shù)據(jù)，用傅立葉級(jí)數(shù)和高次多項(xiàng)式擬合統(tǒng)一的地面高程曲面。也可用局部擬合方法，將地表復(fù)雜表面分成正方形規(guī)則區(qū)域或面積大致相等的不規(guī)則區(qū)域進(jìn)行分塊搜索，根據(jù)有限個(gè)點(diǎn)進(jìn)行擬合形成高程曲面。 2.1 2.1 DEM的表示的表示數(shù)字高程模型-DEM表示法2.圖形方法：線模式線模式等高線是表示地形最常見(jiàn)的形式。其它的地形特征線也是表達(dá)地面高程的重要信息源，如山 脊 線 、 谷 底 線 、 海 岸 線 及 坡 度 變 換 線 等 。

5、點(diǎn)模式點(diǎn)模式用離散采樣數(shù)據(jù)點(diǎn)建立DEM是DEM建立常用的方法之一。數(shù)據(jù)采樣可以按規(guī)則格網(wǎng)采樣，可以是密度一致的或不一致的；可以是不規(guī)則采樣，如不規(guī)則三角網(wǎng)、鄰近網(wǎng)模型等；也可以有選擇性地采樣，采集山峰、洼坑、隘口、邊界等重要特征點(diǎn)。在地理信息系統(tǒng)中，DEM最主要的三種表示模型是：規(guī)則格網(wǎng)模型，等高線模型和不規(guī)則三角網(wǎng)模型。2.1 2.1 DEM的表示的表示DEM表示法2.1 2.1 DEM的表示的表示點(diǎn)模式DEM的幾種表示模型規(guī)則格網(wǎng)模型不規(guī)則三角網(wǎng)模型(TIN: Triangulated Irregular Network)等高線模型層次模型 (Layer of Details,LOD):表

6、達(dá)多種不同精度水平的數(shù)字高程模型。*規(guī)則格網(wǎng)模型和TIN模型是最常用的模型2.1 2.1 DEM的表示的表示點(diǎn)模式DEM的幾種表示模型規(guī)則格網(wǎng)模型規(guī)則網(wǎng)格將區(qū)域空間切分為規(guī)則的格網(wǎng)單元，每個(gè)格網(wǎng)單元對(duì)應(yīng)一個(gè)高程值。數(shù)學(xué)上表示為一個(gè)矩陣，計(jì)算機(jī)實(shí)現(xiàn)中則是一個(gè)二維數(shù)組。 每個(gè)格網(wǎng)的高程值有兩種不同的解釋。格網(wǎng)柵格觀點(diǎn)認(rèn)為格網(wǎng)單元對(duì)應(yīng)的地面面積內(nèi)高程是均一的高度，這種DEM是一個(gè)不連續(xù)的函數(shù)。點(diǎn)柵格觀點(diǎn)認(rèn)為該網(wǎng)格單元的高程值是網(wǎng)格中心點(diǎn)的高程，這樣就需要用距離加權(quán)平均方法來(lái)計(jì)算每個(gè)點(diǎn)的高程。點(diǎn)柵格觀點(diǎn)比較常用。矩陣結(jié)構(gòu)的2.1 2.1 DEM的表示的表示點(diǎn)模式DEM的幾種表示模型TIN模型TIN模型

7、根據(jù)區(qū)域有限個(gè)點(diǎn)集將區(qū)域劃分為相連的三角面網(wǎng)絡(luò)，區(qū)域中任意點(diǎn)落在三角面的頂點(diǎn)、邊上或三角形內(nèi)。如果點(diǎn)不在頂點(diǎn)上，該點(diǎn)的高程值通常通過(guò)線性插值的方法得到。TIN是一個(gè)三維空間的分段線性模型，在整個(gè)區(qū)域內(nèi)連續(xù)但不可微。TIN的數(shù)據(jù)存儲(chǔ)方式不僅要存儲(chǔ)每個(gè)點(diǎn)的高程，還要存儲(chǔ)其平面坐標(biāo)、節(jié)點(diǎn)連接的拓?fù)潢P(guān)系，三角形及鄰接三角形等關(guān)系。頻繁使用指針。2.1 2.1 DEM的表示的表示點(diǎn)模式DEM的幾種表示模型等高線模型等高線模型表示高程，每一條等高線對(duì)應(yīng)一個(gè)已知的高程值，這樣一系列等高線集合和它們的高程值一起就構(gòu)成了一種地面高程模型 。等高線通常用二維的鏈表來(lái)存儲(chǔ)，其被存成一個(gè)有序的坐標(biāo)點(diǎn)對(duì)序列，可以認(rèn)為是

8、一條帶有高程值屬性的簡(jiǎn)單多邊形或多邊形弧段。由于等高線模型只表達(dá)了區(qū)域的部分（位于等高線上的點(diǎn)）高程值，需要用插值方法來(lái)計(jì)算落在等高線外的其它點(diǎn)的高程。通常只使用外包最鄰近的兩條等高線的高程進(jìn)行插值。2.1 2.1 DEM的表示的表示點(diǎn)模式DEM的幾種表示模型層次模型層次模型是一種表達(dá)多種不同精度水平的數(shù)字高程模型。大多數(shù)層次模型是基于TIN模型的，通常TIN的數(shù)據(jù)點(diǎn)越多精度越高，數(shù)據(jù)點(diǎn)越少精度越低，但數(shù)據(jù)點(diǎn)多則要求更多的計(jì)算資源。所以如果在精度滿足要求的情況下，最好使用盡可能少的數(shù)據(jù)點(diǎn)。層次模型允許根據(jù)不同的任務(wù)要求選擇不同精度的TIN地形模型。2.1 2.1 DEM的表示的表示點(diǎn)模式DEM

9、的幾種表示模型間的相互轉(zhuǎn)換 格網(wǎng)TIN 格網(wǎng)等高線 TIN等高線 TIN的形成：Delaunay三角剖分方法2.1 2.1 DEM的表示的表示 DEM的建立與質(zhì)量2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤DEM規(guī)則格網(wǎng)模型的生成等高線點(diǎn)的追蹤2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤規(guī)則格網(wǎng)模型的生成 對(duì)于不規(guī)則分布的高程點(diǎn)，可以通過(guò)距離加權(quán)平均法或其它方法求格網(wǎng)點(diǎn)高程，從而生成格網(wǎng)模型。a.距離加權(quán)平均法 距離加權(quán)平均法基于一種假設(shè)：在一定區(qū)域內(nèi)的某點(diǎn)高程受周圍點(diǎn)高程的影響，其影響的大小與它們之間的距離成反比。 通常需要選擇410個(gè)點(diǎn)來(lái)確定待定點(diǎn)的高程。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法

10、等高線追蹤規(guī)則格網(wǎng)模型的生成a.距離加權(quán)平均法以待定高程的格網(wǎng)點(diǎn)為中心（圓心），在預(yù)先設(shè)定的搜索圓半徑范圍內(nèi)，從小到大進(jìn)行高程點(diǎn)的搜索。1.當(dāng)搜索圓內(nèi)的高程點(diǎn)數(shù)目介于預(yù)先設(shè)定的數(shù)目（通常為410）時(shí)，搜索成功，將搜索到的高程點(diǎn)用以格網(wǎng)點(diǎn)高程的距離加權(quán)平均計(jì)算。否則，搜索失敗，初始高程點(diǎn)數(shù)據(jù)有缺陷。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤規(guī)則格網(wǎng)模型的生成a.距離加權(quán)平均法22)0()0(yyixxiDi設(shè)格網(wǎng)點(diǎn)坐標(biāo)為(x0,y0)，搜索到的高程點(diǎn)坐標(biāo)為(xi, yi, zi) (4i10, z為高程)，高程點(diǎn)到格網(wǎng)點(diǎn)的距離為Di則格網(wǎng)點(diǎn)高程為：)/1()/(DiDiziz2.22.2格網(wǎng)法

11、等高線追蹤格網(wǎng)法等高線追蹤規(guī)則格網(wǎng)模型的生成b.其它求格網(wǎng)點(diǎn)高程的方法 改進(jìn)的距離加權(quán)方法 普通的距離加權(quán)平均法在原始高程數(shù)據(jù)點(diǎn)均勻分布的情況下是合理的。但是，有時(shí)原始數(shù)據(jù)的分布并不理想，搜索到的原始高程點(diǎn)可能集中在某一個(gè)方向，甚至是在一個(gè)十分狹小的方向內(nèi)，而在其它方向則沒(méi)有點(diǎn)。在這種情況下，需要擴(kuò)大搜索半徑，并分象限（扇面）去確定和選擇搜索點(diǎn)。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤規(guī)則格網(wǎng)模型的生成b.其它求格網(wǎng)點(diǎn)高程的方法多項(xiàng)式內(nèi)插線性：搜索最近3個(gè)高程點(diǎn)，插值函數(shù)雙線性：搜索最近4個(gè)高程點(diǎn)，插值函數(shù)其它多項(xiàng)式（已知高程點(diǎn)數(shù)目必須等于待定系數(shù)個(gè)數(shù)） *對(duì)格網(wǎng)點(diǎn)鄰近點(diǎn)的搜索，有時(shí)并不唯

12、一地以距離的遠(yuǎn)近作為取舍的依據(jù)。適當(dāng)考慮高程點(diǎn)的方位分布，使搜索點(diǎn)環(huán)布格網(wǎng)點(diǎn)四周為好。012zaaxay0123zaaxayaxy2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤A.等高線點(diǎn)位置的確定確定高程為z0的等高線點(diǎn)是否在格網(wǎng)邊上僅當(dāng)z0值介于邊的2個(gè)端點(diǎn)高程值之間時(shí)，等高線z0才通過(guò)該邊。等高線z0通過(guò)某邊的判別式L=(z0-z1)(z0-z2)L0時(shí)，等高線不過(guò)該邊。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤A.等高線點(diǎn)位置的確定確定高程為z0的等高線點(diǎn)在格網(wǎng)邊上的平面位置如果確認(rèn)等高線z0通過(guò)格網(wǎng)邊，則按等高線線性內(nèi)插線性內(nèi)插原則確定等高線點(diǎn)的平面位置。

13、關(guān)于格網(wǎng)邊界直線的參數(shù)方程是：x=x1+(x2-x1)t; y=y1+(y2-y1)t; z=z1+(z2-z1)t; 令z=z0，解得t0；將t0回代入?yún)?shù)方程解得(x0, y0)，其就是等高線點(diǎn)的平面位置。21210110112121();()xxyyxxzzyyzzzzzz2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤二義性問(wèn)題當(dāng)一個(gè)格網(wǎng)的4條邊上都有等高線點(diǎn)時(shí)，這時(shí)會(huì)產(chǎn)生等高線追蹤的二義性問(wèn)題（如僅12條邊有等高線點(diǎn)，則無(wú)此問(wèn)題。不可不可能能3 3條邊有等高線點(diǎn)條邊有等高線點(diǎn)）。如圖(a) ,(b)所示(z=50)，2種等高線追蹤都是可能合理，但真正正

14、確的方式應(yīng)該只有一種。一條等高線完全包含在一個(gè)矩形單元內(nèi)非常少見(jiàn)。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤二義性問(wèn)題 通過(guò)在格網(wǎng)的中心位置補(bǔ)測(cè)一個(gè)高程點(diǎn)，構(gòu)成4個(gè)三角形來(lái)解決二義性問(wèn)題。 可以認(rèn)為距離近的連線方式是正確的追蹤方向； 參考地性線的走向來(lái)決定。1.縮小格網(wǎng)尺寸，重建格網(wǎng)模型。（1）（2）（3）2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤微量調(diào)整格網(wǎng)點(diǎn)高程問(wèn)題如果一個(gè)網(wǎng)格點(diǎn)的高程數(shù)值恰好等于要提取的等高線的高程數(shù)值（整數(shù)），會(huì)使等高線追蹤的判斷過(guò)程變得復(fù)雜，并且會(huì)生成不閉合的等高線，一般的解決辦法是將部分（高

15、程為整數(shù)）或全部的網(wǎng)格點(diǎn)的數(shù)值增加一個(gè)小的偏移量偏移量 （ = =最小高程度量值的1/10或1/100）。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤追蹤算法1）設(shè)立格網(wǎng)標(biāo)志數(shù)組，其初始值為0。其元素與格網(wǎng)對(duì)應(yīng)，凡搜索處理過(guò)的格網(wǎng)將標(biāo)志置為1，以后不再被搜索。當(dāng)?shù)雀呔€的高程數(shù)值改變，格網(wǎng)標(biāo)志數(shù)組應(yīng)該清零。DEM的格網(wǎng)尺寸一般保證相應(yīng)等高距的等高線不會(huì)在一個(gè)格網(wǎng)內(nèi)閉合。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤追蹤算法2）按格網(wǎng)點(diǎn)在計(jì)算機(jī)中的順序，依次判斷格網(wǎng)邊中是否有等高線穿過(guò)。直至搜索到等高線與格網(wǎng)邊的第一個(gè)交點(diǎn)（該點(diǎn)

16、為搜索起點(diǎn)）。線性內(nèi)插該點(diǎn)的平面坐標(biāo)，同時(shí)標(biāo)記該格網(wǎng)邊所在的一個(gè)格網(wǎng)作為起始格網(wǎng)（起始格網(wǎng)的選擇也是對(duì)搜索方向的選擇），另一個(gè)格網(wǎng)為二次搜索格網(wǎng)。在搜索過(guò)程中，如果格網(wǎng)點(diǎn)高程為等高距的整數(shù)倍，則對(duì)其作高程值微量調(diào)整。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤B.等高線點(diǎn)位置的追蹤追蹤算法3）依次對(duì)該格網(wǎng)的另外三條邊進(jìn)行搜索，尋找等高線在該格網(wǎng)的離去邊，也就是相鄰格網(wǎng)的進(jìn)入邊，并內(nèi)插等高線點(diǎn)的平面坐標(biāo)。在搜索過(guò)程中，一定要對(duì)另外三條邊全部進(jìn)行搜索判斷。如果一個(gè)格網(wǎng)的四條邊都有等高線穿過(guò)，則要按照預(yù)設(shè)的方式來(lái)解決二義性問(wèn)題。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤

17、B.等高線點(diǎn)位置的追蹤追蹤算法4）進(jìn)入相鄰格網(wǎng)，重復(fù)第（3）步，直至該格網(wǎng)不存在等高線的離去邊（此時(shí)等高線為開(kāi)曲線）或該格網(wǎng)為起始格網(wǎng)（此時(shí)等高線為閉曲線）為止。5）對(duì)于開(kāi)曲線，將已搜索到的等高線點(diǎn)順序倒過(guò)來(lái)；并回到搜索起點(diǎn)，由二次搜索格網(wǎng)開(kāi)始進(jìn)行搜索，直至到達(dá)搜索邊界。6）當(dāng)一條等高線的等高線點(diǎn)全部跟蹤完，將其光滑輸出。然后繼續(xù)格網(wǎng)的搜索，直至全部格網(wǎng)處理完，再改變等高線高程，重復(fù)以上過(guò)程，直到完成全部等高線的繪制為止。2.22.2格網(wǎng)法等高線追蹤格網(wǎng)法等高線追蹤等高線點(diǎn)的追蹤等高線點(diǎn)的插值、光滑 由于DEM格網(wǎng)模型尺寸的原因，一些DEM的分辨率受到限制。追蹤到的等高線點(diǎn)在平滑輸出時(shí)可能視覺(jué)

18、效果較差，或者和實(shí)際地形相差較大。因此有必要在追蹤到的等高線點(diǎn)間進(jìn)行一些高程點(diǎn)的增補(bǔ)，然后再作曲線光滑。 高程點(diǎn)的增補(bǔ)和光滑有如下要求：l曲線要通過(guò)等高線追蹤點(diǎn)l曲線在等高線追蹤點(diǎn)處的一階（或二階）導(dǎo)數(shù)連續(xù)l相鄰的等高線追蹤點(diǎn)間的曲線沒(méi)有多余擺動(dòng)l同一等高線自身不能相交2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成等高線點(diǎn)的追蹤2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成對(duì)于TIN模型的建構(gòu)，有如下基本要求生成的TIN必須具有唯一性力求最佳的三角形形狀，三角形形狀以等邊為最優(yōu)，避免小的銳角和大的鈍角保證最鄰近的點(diǎn)構(gòu)成三角形，即三角形的3邊之和最小Delaunay三角網(wǎng)

19、生成法就是大致滿足上述要求的一種常用方法。其在地形擬合方面十分出色。2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成Voronoi(沃羅努瓦)圖由一組連續(xù)多邊形(又稱泰森多邊形)組成，多邊形的邊界是由連接兩鄰點(diǎn)線段的垂直平分線組成。圖中紅點(diǎn)為原始高程數(shù)據(jù)點(diǎn)，虛黑線為Voronoi多邊形的邊界 線 。 具 有 公 共 邊 線 的Voronoi多邊形為鄰接多邊形。點(diǎn)與其所在多邊形為包含關(guān)系。2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成Delaunay三角形是由每每3 3個(gè)兩兩相鄰的個(gè)兩兩相鄰的VoronoiVoronoi多邊形多邊形所包含的點(diǎn)連接而成的。三角形的外接圓圓心是

20、多邊形的公共頂點(diǎn)。 2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成對(duì)于給定的初始點(diǎn)集P，Delaunay三角網(wǎng)有以下特性：只要不超過(guò)4個(gè)鄰域點(diǎn)共圓，Delaunay三角網(wǎng)是唯一的三角網(wǎng)的外邊界構(gòu)成了點(diǎn)集P的凸多邊形“外殼”；沒(méi)有任何點(diǎn)在三角形的外接圓內(nèi)部。反之，如果一個(gè)三角網(wǎng)滿足此條件，它就是Delaunay三角網(wǎng)。Delaunay三角網(wǎng)并不能理想地滿足構(gòu)建TIN三角網(wǎng)的3點(diǎn)要求，但總體上它趨于最佳，為構(gòu)建TIN三角網(wǎng)的合適選擇。為了進(jìn)一步優(yōu)化三角網(wǎng)的生成，可以采用改進(jìn)的Delaunay三角網(wǎng)生成方法。2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成改進(jìn)的Delaunay

21、三角網(wǎng)生成方法最大化最小角原則每?jī)蓚€(gè)相鄰的三角形所構(gòu)成的凸四邊形中，在對(duì)角線互換后，6個(gè)內(nèi)角中的最小角不再增大，則原三角網(wǎng)較優(yōu)。2.3 2.3 TIN等高線追蹤等高線追蹤TIN模型的生成將一系列高程點(diǎn)組成的點(diǎn)集P轉(zhuǎn)成TIN，最常用的方法是用Delaunay三角剖分方法。生成過(guò)程分兩步完成：利用P中點(diǎn)集的平面坐標(biāo)產(chǎn)生Delaunay三角網(wǎng)；具體算法參見(jiàn)文獻(xiàn)2D-Delaunay 三角網(wǎng)格的數(shù)據(jù)結(jié)構(gòu)與遍歷1.給Delaunay三角形中的節(jié)點(diǎn)賦予高程值。（需要內(nèi)插嗎？）TIN模型生成方法 間接構(gòu)網(wǎng)：先生成泰森多邊形，再構(gòu)建三角網(wǎng) 直接構(gòu)網(wǎng)：離散點(diǎn)直接參與構(gòu)建三角網(wǎng)，如三角形生長(zhǎng)法2.3 2.3 TI

22、N等高線追蹤等高線追蹤等高線點(diǎn)位置的確定設(shè)高程為z的等高線點(diǎn)，通過(guò)三角形邊的兩個(gè)端點(diǎn)的三維坐標(biāo)分別為(x1,y1,z1)和(x2,y2,z2)，則等高線點(diǎn)的平面坐標(biāo)為：211121211121()()zzxxxxzzzzyyyyzzzz2.3 2.3 TIN等高線追蹤等高線追蹤等高線點(diǎn)的追蹤按記錄的三角形順序搜索，其過(guò)程如下：對(duì)給定的等高線高程h，將其與所有三角網(wǎng)的網(wǎng)點(diǎn)高程zi (i=1, 2, , n)，進(jìn)行比較，若zi=h，則將zi加上（或減）一個(gè)微小正數(shù)。(微量調(diào)整)設(shè)立三角形標(biāo)志數(shù)組，其初始值為0，每一元素與一個(gè)三角形對(duì)應(yīng)，凡處理過(guò)的三角形將標(biāo)志置為1，以后不再處理。等高線高程改變后，

23、三角形標(biāo)志數(shù)組應(yīng)清零。1.按順序判斷每一個(gè)三角形的三邊中的兩條邊是否有等高線穿過(guò)。直至搜索到等高線與網(wǎng)邊的第一個(gè)交點(diǎn)，稱該點(diǎn)為搜索起點(diǎn)。2.3 2.3 TIN等高線追蹤等高線追蹤等高線點(diǎn)的追蹤在三角形的另兩邊上搜索等高線在該三角形的離去邊（相鄰三角形的進(jìn)入邊），并內(nèi)插其平面坐標(biāo)。進(jìn)入相鄰三角形，重復(fù)第4步，直至沒(méi)有離去邊（此時(shí)等高線為開(kāi)曲線）或相鄰三角形即搜索起點(diǎn)所在的三角形（此時(shí)等高線為閉曲線）時(shí)為止。對(duì)于開(kāi)曲線，將已搜索到的等高線點(diǎn)順序倒過(guò)來(lái)，并回到搜索起點(diǎn)向另一方向搜索，直至到達(dá)邊界。4.當(dāng)一條等高線全部跟蹤完后，將其光滑輸出，方法與矩形格網(wǎng)等高線的繪制相同。然后繼續(xù)三角形的搜索，直至全部三角形處理完，再改變等高線高程，重復(fù)以上過(guò)程，直到完成全部等高線的繪制為止。2.4 2.4 TINTIN和格網(wǎng)模型的應(yīng)用比較和格網(wǎng)模型的應(yīng)用比較格網(wǎng)模型和TIN是應(yīng)用最廣的兩種數(shù)字地形表示方法。從數(shù)據(jù)存儲(chǔ)結(jié)構(gòu)、地形分辨率表達(dá)以及計(jì)算機(jī)運(yùn)算處理上，它們相比較，