DEM数据获取方法

1、一、DE嗷据获取方法：定义：地形图指的是地表起伏形态和地物位置、形状在水平面上的地物和地貌按水平投影的方法，并按照一定的比例缩绘到图纸上，这种图称为地形图。特点：(1) 具有统一的大地坐标系统的高程系统具有完整的比例尺系列和分幅编号系统：国家基本地形图含1:5千、1:1万、1:2:2.5/1:5万、1:10万、1:25万、1:50万、1:100万8种比例地形图。缺点：地形图现势性较差：纸质地形图制作工艺复杂，更新周期比较长，一般不及时反映局部地形地貌的变化情况地形图存储介质单一，容易变形：传统地形图多为纸质存储介质，存放环境(温湿度)导致地形图图幅产生不同程度的变形，这种变形表现在不同方向上的

2、长度变形和图幅面积上的变形地图精度有限：地图精度决定这地形图对实际地形表达的可信度，与地形图比例尺、等高线密度(由等高距表示)，成图方法有关。不同比例尺的地形图，其所表示的几何精度和内容详细程度有很大的差别。在应用DEM勺时候要考虑DE附辨率、存储格式、数据精度和可信度等因素。二、DE嗷据采样策略与采样方法：采样：确定在何处需要测量点的过程，这个过程有三个参数。决定：点的分布、点的密度和点的精度。采样数据的分布：由数据位置和结构（分布）来确定，指数据点的分布形态位置有地理坐标系统中经纬度或者网格坐标系统中坐标决定。结构的形式很多，因地形特征、设备、应用的不同而不同。采样数据的分布现岫分布二惟晚

3、峭幡河橙矩形格河分布柔惮肇蚩点n校正有形格向令布陶捧M点，特殊规则分布按六边指。布采样数据女一推分布剖面卞规则分布沿等高线栗祥数般点链表分布沿制型at等特征魄分布采样数制随机分布随机分布栗#点1. 数据的密度：是指采样数据密集程度，与研究区域的地貌类型和地形复杂程度有关。用于刻画地形形态所必须的最少的数据点。表示方式：相邻的两点之间的距离、单元面积内的点数、截止频率（采样数据所能表示的最高频率）、单位线段上的点数等。米样距离：相邻两点之间的距离，也称米样间隔。通常数字加单位来表示，如采样距离为20米，表示规格网分布的采样数据另一种表示法是单位面积内的点数，如每平方米500点，描述随机分布的采样

4、数据描述数据分布是沿等高线或特征等线状分布采样点，常用单位线段上的点数，如每米2点数据的精度：是指数据点本身所具有的精确度，是数据获取过程中各种不同类型误差的综合反映采样数据精度与数据源、数据的采集方法和数据采集的一起密切相关。野外测量、影像、地形图扫描的精度从高到低激光扫描、干涉雷达的精度是非常高的摄影测量比GP郭精度要高，达到厘米级。地形图的手扶跟踪和扫描矢量化的精度都比较低三、DE嗷据采集方法：1）手扶跟踪数字化步骤：定参考点（固定地图）定控制点跟踪采集（点方式，流方式）2）扫描数字化/矢量化（详见矢量化文档）DEM（全球数字高程模型）生成方法-ANUDEM模型水是地貌形成的主要侵蚀因素

5、。ANUDEM（AustralianNationalUniversityDigitalElevationModel）采用了这一思想，使用地貌与水文数据作为插值约束条件，插值等高线高程。大大减少了DE时的凹陷数据错误，显著提高了DE雁水文特征方面的质量。插值方法为递归有限元差分插值技术，拟合样条曲面。ANUDEM值处理方法：等高线是最常见的高程信息表达方法，也最难适用各种通用插值方法进行处理，处理要点：(1)地表曲面形态：等高线->等高线局部最大曲率->坡度最陡区域->山脊线/河流径网->确认输出DEM的水文地貌特征/验证DEM4)精确度。(2)地表曲面插值：每个网格单元

6、的插值计算，使用临近等值线上的点。多分辨率逐级插值：低分辨率DEM-流分辨率DEM-翔户定义分辨率DEM每级分辨率插值操作，水系限制条件都用于DEMH陷生成的控制，保留下的凹陷会纪录在log文件中等高线插值的举例(arcinfo:topogrid)Arc:topogridlaodem2020TopoGrid:datatypecontourTopoGrid:contourlaocontourelevTopoGrid:streamlaostreamTopoGrid:enforceonTopoGrid:endArc:其中：laodem20-输出DEM勺文件名20-用户指定分辨率Laocontour-

7、用于插值的等高线Elev-等高线cover的属性表中，表示高程属性的字段名Laostream-河流网络coverEnforceon-清除凹陷二、基于电源数据生成DEM勺插值方法数据源：GPSM点，地形图网格采点等。Arc:itemspntcov.patCOLUMNITEMNAMEWIDTHOUTPUTTYPEN.DECALTERNATENAME9PNTCOV#45B13PNTCOV-ID45B21YIELD88F629X_COORD810F337Y_COORD811F3Arc:q反距离加权插值法IDW反距离加权法中距离的指数，指数越小表面越光滑，常用范围为0.53,缺省为2。Arc:gridG

8、rid:pntgrdidw=idw(pntcov,yield,#,2,sample,#,#,2.5)Running.100%Grid:样条插值法Tension表示张力样条插值法，Regularized表示规则样条插值法Grid:pntgrdspline=spline(pntcov,yield,tension,#,#,2.5)Running.100%Grid:克吕格插值法Kriging方法包括普通Kriging(SPHERICAL,CIRCULAR,EXPONENTIAL,GAUSSIAN,LINEAR)乏Kriging(UNIVERSAL1UNIVERSAL2)Grid:pntcovkrigi

9、ng=kriging(pntcov,yield,#,#,#,SPHERICAL,SAMPLE,#,#,2.5)Estimatingsemi-variogram.OrdinaryKrigingwiththeFollowingModel(s):SPHERICALc0=0.003c=0.004a=208.237sill=0.007Interpolatinggrid.100%Grid:三.常用DEM据资源GTOPO30DEM1)GTOPO30据特点由USGSEROS成于1996年，历时3年GTOPO澈据源：DTM,DCW(VMAP0),USGS1degreeDEM覆盖全球陆地90N-90S,小于1平方

10、公里小岛屿忽略数据库大小：1.74Giga(行21600列43200)分辨率：30arcseconds(约1公里)高程范围：-4078,752meter数据格式：16bit二进制格式bil,MSB-bigendian2)Verticalaccuracy(meters)SourceL.E.at90%RMSEEstimationmethodDTED3018productspecificationDCW16097calculatedvs.DTEDUSGSDEM3018productspecificationAMSmaps250152estimatedfrom500-meterintervalIMWm

11、aps5030estimatedfrom100-meterintervalPerumap500304estimatedfrom1,000-meterintervalN.Z.DEM159estimatedfrom100-footintervalADDhighlyvariabwidelerangeofscalesandintervalsGTOPO30据的获取整个数据库分为33tiles，中国占据4tiles:E60N90E100N90E60N40E100N40LatLonElevTileMinMaxMinMaxMinMaxMeanStd.Dev.E060N90409060100-152716950

12、9698E100N90409010014013877597455E060N40-1040601001875218041892E100N40-1040100140-407213692910GTOPO30据文件解释(以E100N4,例)FileNameContentsE100N40.DEMdigitalelevationmodel(BIP/BIL/BSQ)E100N40.HDRheaderfileforDEME100N40.DMWworldfileE100N40.STXstatisticsfileE100N40.PRJprojectioninformationfileE100N40.GIFshad

13、edreliefimageE100N40.SRCsourcemapE100N40.SCHheaderfileforsourcemapSRTMDE畋据的处理与拼接#(ASCIIDEM情况下，投影参数手工定义)Arc:projectdefinegridn040e116Project:projectiongeographicProject:datumwgs84Project:zunitsmetersProject:unitsddProject:spheroidwgs84Project:parametersArc:#DEMMosaicArc:gridGrid:mosaicUsage:(*)MOSAIC(<grid,.,grid>)Grid:bjdem=mosaic(0N39E116,N039E117,N040E116,N040E117)Mosaicinggrids.100%#SetShadow/WatertoNULL(foreasydisplayonly)Grid:bjdemp=setnull(b