基于等高線的DEM生成方法研究與實現(xiàn)

一、等高線原始數(shù)據(jù)存儲分析

1、等高線數(shù)據(jù)以shape格式存儲，包括等高線和高程點兩個文件。

2、在等高線文件中，等高線以polyline格式存儲，高程點以point格式存儲。

3、其中存儲的polyline包括首尾重合和首尾不重合兩種類型。

4、首尾重合的polyline記錄了封閉的等高線。

5、首尾不重合的polyline記錄了被邊界裁剪后的等高線，同一等高線可能被分為多段。且存儲不具有連續(xù)性。

小結：

由于無邊界數(shù)據(jù)，故難以通過現(xiàn)有數(shù)據(jù)復原不封閉的等高線。從而采用根據(jù)等高線上的已知數(shù)據(jù)計算網(wǎng)格點上未知數(shù)據(jù)的方法方法生成DEM數(shù)據(jù)。

二、基于等高線的三維地形造型方法

目前基于等高線的三維地形造型方法主要包括如下三種模型：

1、均勻網(wǎng)格模型

實現(xiàn)方法：數(shù)據(jù)由一系列等間隔的等高線高程值表示，代表一塊方形網(wǎng)格地形；網(wǎng)格交叉點的高程值由相鄰等高線的高程值插值得到。

優(yōu)點：易于生成具有多分辨率的層次細節(jié)模型，易于實現(xiàn)地表特征物與地形間的自適應匹配。

缺點：不能直接利用等高線的高度信息，對于地形平坦地區(qū)網(wǎng)格浪費，對于起伏較大的地區(qū)由于精度不夠而丟失細節(jié)。

2、非均勻網(wǎng)格模型

實現(xiàn)方法：根據(jù)等高線上非均勻分布的采樣點在地平面上的投影點所作的水平和垂直線來構成地平面的非均勻分割。

優(yōu)點：充分利用了等高線上的高度信息，具有較高的精確性。

缺點：劃分的網(wǎng)格數(shù)較多，實時性較差。

3、三角模型

實現(xiàn)方法：直接利用原始離散采樣點表示地形，并將所有的采樣點連成連續(xù)三角形，三角形形狀和大小取決于采樣點的位置和密度。

優(yōu)點：數(shù)據(jù)冗余少，能夠突出地形變化，繪制的三D立體地形圖效果好。

缺點：數(shù)據(jù)結構復雜

小結：

綜合對比上述三種方法，結合項目不同比例尺下顯示不同精度的高程數(shù)據(jù)的需求，擬采用均勻網(wǎng)格模型。

三、插值算法

目前高程內(nèi)存算法主要有：窗口擬合法、移動擬合法和線性插值法。

1、窗口擬合法

計算以待插入點為圓心，以 R 為半徑的圓內(nèi)諸點的加權平均值來確定高程值, 精度由已知點到內(nèi)插點的距離來確定。

2、移動擬合法

移動擬合法是典型的逐點內(nèi)插法, 該方法對每個待插入點用一個多項式曲面來擬合該點附近的地表面。

3、線性插值

線性插值可分為剖面內(nèi)插和比例內(nèi)插，剖面內(nèi)插中每一個點的值由剖面上各等高線值內(nèi)插得到；比例內(nèi)插中每一個網(wǎng)格點的值由離該點最近的兩條等高線的值內(nèi)插得到。

小結：

綜合對比上述三種算法，線形插值算法實現(xiàn)簡單，且能夠滿足工程需要，故擬采用線形插值法。具體算法如下：

(1) 求取橫向和縱向兩個方向的插值，將二者的平均值作為該點高程值。

(2) 若待插入點恰好在等高線上，則直接賦值。否則通過插值得到。插值公式如下：

z = (z1 * d1 + z2 * d2) / (d1 + d2)

(3) 若待插入點兩側存在網(wǎng)格線與等高線的交點，且高度值不同，則進行內(nèi)插(P1點)。

(4) 若待插入點兩側與等高線相交，且高度值相等，則求取外側第二個交點，若存在第二個交點，高度值由距離較近的同側的兩個高度值外插得到(P2)；若兩側均不存在第二個交點，則運用第一組交點進行插值(P4)。

(5) 若待插入點只有一側與等高線相交(邊界)，則在存在交點一側求取兩個交點，高度值由同側的兩個高度值外插得到(P3)。若整條切割線上只存在一個交點，則整條切割線上的所有點均賦值為該等高線值(P5)。

(6) 若待插入點兩側均不與等高線相交，則本次不進行計算，在第二輪時運用周邊計算完畢的點進行插值(P6)。

四、數(shù)據(jù)制作流程

在運用插值方法計算網(wǎng)格點的高程值時，受周邊等高線的值和距離等高線的距離兩個因素的影響。當對數(shù)據(jù)分割后進行計算時，可能導致不同分塊邊界上的同一個點計算出不同的高度值，故采用先制作網(wǎng)格數(shù)據(jù)，后進行BlockSet、Block、Parcel分割的策略。數(shù)據(jù)制作流程如下：

1、解析 shape 文件，以數(shù)組結構保存所有等高線。

2、采用線性插值算法制作高程網(wǎng)格數(shù)據(jù)。

3、分割制作完畢的高程網(wǎng)格數(shù)據(jù)。

對分割完畢的高程網(wǎng)格數(shù)據(jù)進行三角化。