1、一、DEM數據獲取方法：定義：地形圖指的是地表起伏形態和地物位置、形狀在水平面上的地物和地貌按水平投影的方法，并按照一定的比例縮繪到圖紙上，這種圖稱為地形圖。特點：（1） 具有統一的大地坐標系統的高程系統（2） 具有完整的比例尺系列和分幅編號系統：國家基本地形圖含1:5千、1:1萬、1：2：2.5/1:5萬、1:10萬、1:25萬、1:50萬、1:100萬8種比例地形圖。缺點：（1） 地形圖現勢性較差：紙質地形圖制作工藝復雜，更新周期比較長，一般不及時反映局部地形地貌的變化情況（2） 地形圖存儲介質單一，容易變形：傳統地形圖多為紙質存儲介質，存放環境（溫濕度）導致地形圖圖幅產生不同程度的變形，

2、這種變形表現在不同方向上的長度變形和圖幅面積上的變形（3） 地圖精度有限：地圖精度決定這地形圖對實際地形表達的可信度，與地形圖比例尺、等高線密度（由等高距表示），成圖方法有關。不同比例尺的地形圖，其所表示的幾何精度和內容詳細程度有很大的差別。在應用DEM的時候要考慮DEM分辨率、存儲格式、數據精度和可信度等因素。二、DEM數據采樣策略與采樣方法：采樣：確定在何處需要測量點的過程，這個過程有三個參數。決定：點的分布、點的密度和點的精度。1. 采樣數據的分布：由數據位置和結構（分布）來確定，指數據點的分布形態位置有地理坐標系統中經緯度或者網格坐標系統中坐標決定。結構的形式很多，因地形特征、設備、應

3、用的不同而不同。2. 數據的密度：是指采樣數據密集程度，與研究區域的地貌類型和地形復雜程度有關。用于刻畫地形形態所必須的最少的數據點。表示方式：相鄰的兩點之間的距離、單元面積內的點數、截止頻率（采樣數據所能表示的最高頻率）、單位線段上的點數等。采樣距離：相鄰兩點之間的距離，也稱采樣間隔。·通常數字加單位來表示，如采樣距離為20米，表示規格網分布的采樣數據·另一種表示法是單位面積內的點數，如每平方米500點，描述隨機分布的采樣數據·描述數據分布是沿等高線或特征等線狀分布采樣點，常用單位線段上的點數，如每米2點3. 數據的精度：是指數據點本身所具有的精確度，是數據獲取

4、過程中各種不同類型誤差的綜合反映采樣數據精度與數據源、數據的采集方法和數據采集的一起密切相關。野外測量、影像、地形圖掃描的精度從高到低激光掃描、干涉雷達的精度是非常高的攝影測量比GPS的精度要高，達到厘米級。地形圖的手扶跟蹤和掃描矢量化的精度都比較低三、DEM數據采集方法：1） 手扶跟蹤數字化步驟：定參考點（固定地圖）定控制點跟蹤采集（點方式，流方式）2） 掃描數字化/矢量化（詳見矢量化文檔）DEM（全球數字高程模型）生成方法 - ANUDEM 模型 水是地貌形成的主要侵蝕因素。ANUDEM (Australian National University Digital Elevation M

5、odel) 采用了這一思想，使用地貌與水文數據作為插值約束條件，插值等高線高程。大大減少了DEM中的凹陷數據錯誤，顯著提高了DEM在水文特征方面的質量。插值方法為遞歸有限元差分插值技術，擬合樣條曲面。ANUDEM插值處理方法： 等高線是最常見的高程信息表達方法，也最難適用各種通用插值方法進行處理，處理要點： (1) 地表曲面形態：等高線->等高線局部最大曲率->坡度最陡區域->山脊線河流徑網->確認輸出DEM的水文地貌特征/驗證DEM的精確度。 (2) 地表曲面插值：每個網格單元的插值計算，使用臨近等值線上的點。(3) 多分辨率逐級插值：低分辨率DEM->高分辨率

6、DEM->用戶定義分辨率DEM。每級分辨率插值操作，水系限制條件都用于DEM凹陷生成的控制，保留下的凹陷會紀錄在log文件中等高線插值的舉例 (arcinfo: topogrid) Arc: topogrid laodem20 20 TopoGrid: datatype contour TopoGrid: contour laocontour elev TopoGrid: stream laostream TopoGrid: enforce on TopoGrid: endArc:其中：laodem20 - 輸出DEM的文件名20 - 用戶指定分辨率Laocontour - 用于插值的等

7、高線Elev-等高線cover的屬性表中，表示高程屬性的字段名Laostream-河流網絡coverEnforce on-清除凹陷2、 基于電源數據生成DEM的插值方法數據源：GPS測點，地形圖網格采點等。 Arc: items pntcov.pat COLUMN ITEMNAME WIDTH OUTPUT TYPE N.DEC ALTERNATENAME 9 PNTCOV# 4 5 B13 PNTCOV-ID 4 5 B21 YIELD 8 8 F 629 X_COORD 8 10 F 337 Y_COORD 8 11 F 3 Arc: q1) 反距離加權插值法 IDW 反距離加權法中距離的

8、指數，指數越小表面越光滑，常用范圍為0.53,缺省為2。 Arc: grid Grid: pntgrdidw=idw(pntcov,yield,#,2,sample,#,#,2.5) Running . 100% Grid: 2) 樣條插值法 Tension表示張力樣條插值法，Regularized表示規則樣條插值法 Grid: pntgrdspline = spline(pntcov,yield,tension,#,#,2.5) Running . 100% Grid: 3) 克呂格插值法 Kriging方法包括普通Kriging ( SPHERICAL, CIRCULAR, EXPONEN

9、TIAL, GAUSSIAN, LINEAR)和泛Kriging(UNIVERSAL1 UNIVERSAL2) Grid: pntcovkriging = kriging(pntcov,yield,#,#,#,SPHERICAL,SAMPLE,#,#,2.5) Estimating semi-variogram . Ordinary Kriging with the Following Model(s): SPHERICAL c0 = 0.003 c = 0.004 a = 208.237 sill = 0.007 Interpolating grid . 100% Grid: 3 常用DEM

10、數據資源 GTOPO30 DEM 1) GTOPO30數據特點 由USGS EROS完成于1996年, 歷時3年 GTOPO30數據源：DTM, DCW(VMAP0), USGS 1 degree DEM 覆蓋全球陸地 90N 90S, 小于1平方公里小島嶼忽略 數據庫大小：1.74 Giga (行21600 列43200) 分辨率: 30 arc seconds (約1公里) 高程范圍: -407 8,752 meter 數據格式：16bit二進制格式bil，MSB - big endian 2) Vertical accuracy (meters)Source L.E. at 90% RM

11、SE Estimation method- - - DTED 30 18 product specification DCW 160 97 calculated vs. DTED USGS DEM 30 18 product specification AMS maps 250 152 estimated from 500-meter interval IMW maps 50 30 estimated from 100-meter interval Peru map 500 304 estimated from 1,000-meter interval N.Z. DEM 15 9 estima

12、ted from 100-foot interval ADD highly variab wide le range of scales and intervals3) GTOPO30數據的獲取 整個數據庫分為33 tiles，中國占據4 tiles： E60N90 E100N90 E60N40 E100N40 Lat Lon Elev Tile Min Max Min Max Min Max Mean Std.Dev. - - - - E060N90 40 90 60 100 -152 7169 509 698E100N90 40 90 100 140 1 3877 597 455 E060

13、N40 -10 40 60 100 1 8752 1804 1892 E100N40 -10 40 100 140 -40 7213 692 910 4) GTOPO30數據文件解釋 (以E100N40為例) FileName Contents - - E100N40.DEM digital elevation model (BIP/BIL/BSQ) E100N40.HDR header file for DEM E100N40.DMW world file E100N40.STX statistics file E100N40.PRJ projection information file

14、E100N40.GIF shaded relief image E100N40.SRC source map E100N40.SCH header file for source map5) SRTM DEM數據的處理與拼接 #(ASCII DEM情況下, 投影參數手工定義) Arc: projectdefine grid n040e116 Project: projection geographic Project: datum wgs84 Project: zunits meters Project: units dd Project: spheroid wgs84 Project: parameters Arc: #DEM Mosaic Arc: grid Grid: mosaic Usage: (*) MOSAIC (<grid, ., grid>) Grid: bjdem = mosaic(0N39E116,N039E117,N040E116,N040E117) Mosaicing grids . 100% #Set Shadow/Water to NULL (fo