DEM格网

合集下载

DEM网格单元大小的确定

DEM网格单元大小的确定简单方法1由地形图上的等高线生成DEM时，DEM网格大小的粗略估计: CELL Size = Scale分母/ 纸张分辨率纸张分辨率为300bpi(一般为200bpi)，即一英寸纸张上面可以印刷300条线，以1:5万地形图为例：cell size = 50000/300 (inch) = 4.24 (meter)方法2地图比例尺，航空摄影测量、影像分辨率的关系带来的启示航摄规范（GB/T 15661-1995）中规定航摄仪有效使用面积内镜头分辨率“每毫米内不少于25 线对”。

根据物镜分辨率和摄影比例尺可以估算出航摄影像上相应的地面分辨率D，即D=M/R。

(其中M 为摄影比例尺分母，R 为镜头分辨率。

)根据航摄规范中“航摄比例尺的选择”的规定和以上公式，可得下表。

成图比例尺航摄比例尺影像地面分辨率（m）1：5000 1：10,000～1：20,000 0.4～0.81：10,000 1：20,000～1：40,000 0.8～1.61：2,5000 1：25,000～1：60,000 1.0～2.41：50,000 1：35,000～1：80,000 1.4～3.2补充：卫星影像分辨率的选择考虑不同比例尺成图对影像分辨率要求和对应规格商用卫星影像产品的稳定货源。

卫星QuickBird-2 IKONOS-2 SPOT-5 SPOT-4 Landsat-7最高分辩率(m) 0.61 1 2.5 10 15成图比例尺卫星影像（分辨率）1：5000～1：10,000 QuickBird(0.61m)IKONOS-2 (1m)1：25,000 QuickBird-2(0.61m)IKONOS-2 (1m)SPOT-5(2.5m)1：50,000 SPOT-5(2.5m)DEM生成方法- ANUDEM 模型水是地貌形成的主要侵蚀因素。

ANUDEM (Australian National University Digital Elevation Model) 采用了这一思想，使用地貌与水文数据作为插值约束条件，插值等高线高程。

第4章-DEM表面建模PPT课件

• 为保证相邻分块之间的平滑连接，相邻分块之间要有一定宽度的重叠。
• 不同的分块单元可用不同的内插函数，常用的内插方法有线性内插、双线性内插、多项式内插、样条函数内插、多层曲面叠加法等。
.
19
4.3.2 局部分块内插
• 1）线性内插和双线性内插 • 将分块单元内部的地形曲面视为平面 • 线性内插： H= ax+by+c
• DEM应首先考虑保凸性和逼真性！
.
8
4.1.3 DEM建立的一般步骤和方法
• （1）采用合适的空间模型构造空间结构。 • 即DEM格网化过程（形成格网）或三角剖分
过程（形成TIN） • （2）采用合适的属性域函数，属性值为高程。 • （3）在空间结构中进行采样，构造空间域函
数（即内插函数的确定） • （4）利用空间域函数进行分析（即求取格网
• 2）二元样条函数内插 • 样条曲面就是将一张具有弹性的薄板压定在各
个采样点上，而其他的地方自由弯曲。 • 数学上讲，为一个分段的低次多项式，多项式
次数一般不超过三阶。 • 二元样条函数首先对采样区进行分块，对每一
块用一个多项式进行拟合，为保证各个分块之间的平滑过渡，按照弹性力学条件设立分块之间的连续性条件，即公共边界上的导数连续。 • 特点：不仅保留了局部地形的细部特征，还能获取光滑的DEM。但由于多项式阶数较低，对误差响应不敏感，具有较好的保凸性和逼真性和平滑性。
• 整体内插函数通常是多项式，要求地形采样点个数大于或等于多项式的系数数目。此时没有唯一解，一般采用最小二乘法求解，即要求多项式曲面与地形采样点之间差值的平方和为最小，属曲面拟合插值或趋势面插值。
• 虽然任何复杂曲面都可以由多项式在任意精度上逼近，但在DEM内插中整体内插并不常用。

不规则三角网(TIN)

不规则三角网（TIN）Ⅰ 数字高程模型（DEM）地球表面高低起伏，呈现一种连续变化的曲面，这种曲面无法用平面地图来确切表示。

于是我们就利用一种全新的数字地球表面的方法——数字高程模型的方法，这种方法已被普遍广泛采用。

数字高程模型即DEM（Digital Elevation Model），是以数字形式按一定结构组织在一起，表示实际地形特征空间分布的模型，也是地形形状大小和起伏的数字描述。

DEM有三种主要的表示模型：规则格网模型，等高线模型和不规则三角网。

格网（即GRID）DEM在地形平坦的地方，存在大量的数据冗余，在不改变格网大小情况下，难以表达复杂地形的突变现象，在某些计算，如通视问题，过分强调网格的轴方向。

不规则三角网（简称TIN，即Triangulated Irregular Network）是另外一种表示数字高程模型的的方法（Peuker等，1978），它既减少了规则格网带来的数据冗余，同时在计算（如坡度）效率方面又优于纯粹基于等高线的方法。

不规则三角网能随地形起伏变化的复杂性而改变采样点的密度和决定采样点的位置，因而它能够避免地形起伏平坦时的数据冗余，又能按地形特征点如山脊，山谷线，地形变化线等表示数字高程特征。

Ⅱ TIN的基本知识在TIN中，满足最佳三角形的条件为：尽可能的保证三角形的三个角都是锐角，三角形的三条边近似相等，最小角最大化。

TIN 是基于矢量的数字地理数据的一种形式，通过将一系列折点（点）组成三角形来构建。

形成这些三角形的插值方法有很多种，例如Delaunay 三角测量法或距离排序法。

ArcGIS 支持Delaunay 三角测量方法。

TIN 的单位是英尺或米等长度单位，而不是度分秒。

当使用地理坐标系的角度坐标进行构建时，Delaunay 三角测量无效。

创建TIN 时，应使用投影坐标系（PCS）。

TIN 模型的适用范围不及栅格表面模型那么广泛，且构建和处理所需的开销更大。

获得优良源数据的成本可能会很高，并且，由于数据结构非常复杂，处理TIN 的效率要比处理栅格数据低。

第四章格网DEM建立

光滑性：光滑性是指曲线上切线方向变化的连续性，或者说曲线上曲率的连续性。
7
Beijing University of Civil Engineering and Architecture
4.1 从散点到地形统计表面
DEM 质量评价标准
光滑性：光滑性是指曲线上切线方向变化的连续性，或者说曲线上曲率的连续性。三方面相互独立又相互影响：曲面的逼真性与保凸性有关，保凸性显然会影响曲面的整体逼近性；而保凸性和光滑性常常矛盾，一个光滑性很好的逼近面可能保凸性较差。不同的应用领域对这些要求的重视程度也不一致，例如实际地形曲面一般是比较粗糙的，DEM 应首先满足保凸性和逼真性，而对于飞机、汽车等制造业而言，首先考虑的却是光滑性。
4.1 从散点到地形统计表面
DEM 质量评价标准
逼真性： MAX ∣f (x, y) − F(x, y) ∣ ≤σ ，则认为逼近面达到逼真性要求。
6
Beijing University of Civil Engineering and Architecture
4.1 从散点到地形统计表面
DEM 质量评价标准
4.2 规则格网DEM建立的基本思路
DEM内插分类方法
11ห้องสมุดไป่ตู้
Beijing University of Civil Engineering and Architecture
4.3 DEM内插数学模型
整体内插定义：在整个区域用一个数学函数来表达地形曲面。
12
Beijing University of Civil Engineering and Architecture
4.4 DEM建立过程
基于不规则分布采样点的DEM建立 1）直接法---邻域和领域内点的确定常用的邻域搜索区域包括搜索圆与搜索正方形两种初始半径经验公式

基于梯度方向的格网DEM的断裂特征线生成

然后根据以上数学分析判断该点是否属于断裂带的
在ｆ
＼
１０况情下，
，。．
（Ｌｔ且Ｌ北。（。
或者
边缘点。通过对所有点的遍历，获得断裂特征线的候选点，之所以这里得到的只称之为候选特征点，
是地学分析的基础，何获取ＤＭ的结构特征始终的结构特征。本文探讨的就是如何分析断裂带边缘如Ｅ
是其中基础而重要的问题。一般情况下，可以通过点，并生成其上下边缘线，为后续断裂特征面对象断面极值分析［１取局部地形特征点，以及最大梯的形成作准备。１获，２
地形断裂线
特征提取
断裂面对象
１引言
在ＤＥ结构特征线分析与连接过程中碰到断裂边Ｍ
地面高程模型（ＥＤＭ，ＤｇｔｌｌｖｔｎＭｏｅ）缘就采用新的处理方式，这样才能形成准确而完整ｉｉｅａｉｄ１ａＥｏ
摘
２０９）００２
要：本文按照断裂地形局部区域的变化特性建立数学模型，先获得表示断裂带边缘候选特征点，然后按
梯度最小方向搜索和局部区域特征线唯一性的地形特征超关系（ＣＲ原则连接分别得到地形断裂区域的上ＴＨ）下沿特征线并形成对偶。具体试验中，我们针对海底地形中台地区域的ＤＭ进行分析和提取，得了实际地Ｅ取形吻合得较满意的断裂特征线。进一步研究，我们将根据面向对象的思想，生成相应独立不规则但具有自适

规则格网DEM的解算及精度检验

双线性插值直接通过网格单元的四个顶点拟合曲面袁所采用的曲面方程为院
11
移移 H=f渊x袁y冤=
aij xiyi
j = 0i = 0
=a00+a10x+a01y+a11xy
=a0+a1x+a2y+a3xy
渊2冤
2.3 双三次多项式内插数学模型
双三次多形式内插通过网格单元周围的顶点依据某种数学方式
进行曲面拟合袁所采用的方程为院
11
移移 H=f渊x袁y冤=
aij xiyi
渊3冤
j = 0i = 0
式中的坐标是各个各网点的直角坐标遥
将当前格网单元与周围格网单元的高程值组成方程组联立求
解袁将点 P 坐标渊xp袁yp冤代入袁就可求出其高程 Hp[3]遥 2.4 反距离加权法内插数学模型
渊1冤根据实际插值要求袁计算 P 点渊Xp袁Yp冤的平面坐标遥渊2冤为了选择相邻的数据点袁使用点 P 为圆心袁R 为半径遥选择落在圆内的所有数据点遥当插值二次曲面时袁考虑到计算的方便性袁将坐标原点移动到 DEM 网格点 P渊Xp袁Yp冤袁由于反距离加权方法的系数数量需要 4~10 个袁因此所需的数据点数量为 4<n<10遥当数据点 Pi
确定权利的原则应该与数据点与等待点之间的距离 D 以及权的形式
渊1冤总结和分析现有的内插模型袁讨论阅耘酝内插模型的数学原理遥渊2冤探讨同一精度指标对于不同内插方法的差异性遥渊3冤编写 DEM 数据内插以及 DEM 精度评定程序遥
2 DEM 空间数据内插原理与方法
2.1 线性内插数学模型
线线性插值将网格单元分割成两个三角形袁每个三角形由三个顶点定义一个唯一平面袁插值过程为遥

格网 DEM 分析的主要应用。

答案：
（1）地形曲面拟合：DEM 最基础的应用是求 DEM 范围内任意点的高程，在此基础上进行地形属性分析。

由于已知有限个格网点的高程，可以利用这些格网点高程拟合一个地形曲面，推求区域内任意点的高程。

（2）立体透视图：绘制透视立体图是 DEM 的一个极其重要的应用。

透视立体图能更好地反映地形的立体形态，非常直观。

人们可以根据不同的需要，对于同一个地形形态作各种不同的立体显示，更好地研究地形的空间形态。

（3）通视分析：通视分析有着广泛的应用背景。

典型的例子是观察哨所的设定、森林中火灾监测点的设定、无线发射塔的设定等。

通视问题可以分为五类：a）已知一个或一组观察点，找出某一地形的可见区域；b）欲观察到某一区域的全部地形表面，计算最少观察点数量；c）在观察点数量一定的前提下，计算能获得的最大观察区域；d）以最小代价建造观察塔，要求全部区域可见；e）在给定建造代价的前提下，求最大可见区。

根据问题输出维数的不同，通视可分为点的通视，线的通视和面的通视。

（4）流域特征地貌提取与地形自动分割：是进行流域空间模拟的基础技术。

主要包括两个方面：a）流域地貌形态结构定义，定义能反映流域结构的特征地貌，建立格网 DEM 对应的微地貌特征；b）特征地貌自动提取和地形自动分割算法。

1
（5）计算地形属性：DEM 派生的地形属性数据可以分为单要素属性和复合属性二种。

前者可由高程数据直接计算得到，如坡度因子，坡向。

后者是由几个单要素属性按一定关系组合成的复合指标，用于描述某种过程的空间变化，这种组合关系通常是经验关系，也可以使用简化的自然过程机理模型。

2。

格网DEM

格网DEM的优点：数据结构简单，便于管理；有利于地形分析，以及制作立体图。格网DEM的缺点：格网点高程的内插会损失精度，格网过大会损失地形的关键特征，如山峰、洼地、山脊等；如不改变格网的大小，不能适用于地形起伏程度不同的地区；地形简单的地区存在大量数据冗余。
格网DEM是DEM最常用的形式，其数据的组织形式类似于图像栅格数据，只是每个像元的值是高程值。即格网DEM是一种高程矩阵，其高程数据可直接由解析立体测图仪获取，也可由规则或不规则的离散点数据内插产生。
100 110 120 140 110 105 90 120 115 130 135 120 110 100 135 120 120 130 130 120 110 145 130 115 20 120 115 118

什么是4D（DRG、DLG、DOM、DEM）数据

什么是4D（DRG、DLG、DOM、DEM）数据1，DOM ，利⽤数字⾼程模型对扫描处理的数字化的航空相⽚、遥感影像抄，经逐个像元纠正，按图幅范围裁切⽣成的影像数据。

百DOM 是需要DEM进⾏⼆次加⼯的，也是4D产品中最为⾼级的产品。

2，DEM ，通过等⾼线、或航空航天影像建⽴以表达地⾯⾼程起伏形态的数字集合。

DEM数据为基础数据。

3，DRG，是纸制地形图的栅格形式的数字化产品，可与DOM、DEM集成派⽣出新的可视信息。

4，DLG，利⽤航空航天影像通过对影像进⾏识别和⽮度量化，建⽴基础地理要素分层存储的⽮量数据集，既包括空间信息也包括属性信息，可⽤于各专业信息系统的空间定位基础。

⼀、 DOM (图)：利⽤数字⾼程模型对扫描处理的数字化的航空相⽚、，经逐个像元纠正，按图幅范围裁切⽣成的影像数据，它的信息⽐较直观，具有良好的可判读性和可量测性，从中可直接提取⾃然地理和社会经济信息。

在SAR图像处理中，往往需借助DEM数据来解决RD定位导致的斜距成像⼏何失真。

因此，求解X,Y,Z考虑了三个⽅程。

即距离公式、多普勒频率公式和地球坐标公式。

也就是说DOM是需要DEM 进⾏⼆次加⼯的，也是4D产品中最为⾼级的产品。

DEM (数字⾼程模型) ：通过等⾼线、或影像建⽴以表达地⾯⾼程起伏形态的数字集合。

⽬前可得到的有90m的SRTM，和30m的Aster GDTM数据。

前者采⽤InSAR技术获取，后者则是⾼分辨率⽴体摄影测量技术。

两者相似之处都需要两幅图像，⽽且精确配准。

需要有⼀定的基线长度，需在⼀定范围内取值。

不同之处，前者是利⽤波的相⼲性原理求得，后者则是光直线传播所产⽣的共线⽅程。

DEM数据为基础数据。

DRG (数字栅格地图) ：数字栅格地图是纸制地形图的栅格形式的数字化产品，可与DOM、DEM集成派⽣出新的可视信息。

该类型数据主要是将已有的纸质地图进⾏栅格化，然后配准，⽬前这类图很少⽤到，多⽤⾼分辨率的影像来取代，或者就是将主要地物进⾏⽮量化表征和存储，⽬前⼤多数的都⽀持这⼀功能。

第4章格网DEM建立

等高线构建TIN法
• 等高线内插法、等高线离散化法相比，等高线构建TIN无论在效率还是内插精度上都是最优的 .
不规则三角网TIN的基本概念
TIN的基本元素
数据和类型
• 无约束数据域 • 约束数据域 • 无约束TIN • 约束TIN
TIN的三角剖分准则
局部优化过程，local optimal procedure ，LOP
• 保证相邻分块之间的平滑连接，相邻分块之间要有一定宽度的重叠
分块单元常用的内插数函数
• 线性内插、双线性内插
• 多项式内插 • 样条函数 • 多层曲面叠加法
线性内插、双线性内插
双线性插值
线性插值(假设：高程在直线上呈比例变化)
B C
| AC | f c f a ( f b f a) | AB | 基本公式：
不规则点双线性插值 C E D P B A 规则点双线性插值
A
A
E
B
P
D F C
样条插值
思路：利用样本点，采用一定的数学函数拟合出一个光滑曲面，使该曲面曲率最小，且通过所有的样本点。相当于扭曲一条橡皮让它通过所有的点，同时曲率最小。需具备满足弹性力学光滑性条件，但地形不具备（不理想）
保凸性
逼真性
光滑性
DEM质量评价标准
• 高质量的逼近曲面应该具有较好的保凸性和逼真性，如果实际曲面是光滑的，则逼近面也应该是光滑的 • 不同的应用领域对这些要求的重视程度也不一致（DEM和飞机、汽车等制造业）
DEM建立的一般步骤与方法
• 建模的目的是对复杂的客体进行简化和抽象，并把对客体（源域，DEM中为地形起伏）的研究转移到对模型的研究来

DEMDOMDLG生产流程

摄影负片
象片扫描 GPS外业控制象片参数
GPS自动空中三角测量
单个象对定向或安置（内定向，相对定向，绝对定向）
自动相关获取DEM
检查编辑修改
多模型DEM镶嵌
质量检测与元数据文件记录 DEM建库或刻盘
b.从现有地形图获取数据的流程

如果采用原有的地图，或使用未经改造的模拟精密立体测图仪获得的等高线，可以扫描，先进行矢量化，得到矢量化的等高线，再内插成DEM。该种方法由于使用的设备相对简单，所获得的DEM精度，与前面所述的解析测图仪方法相当，可以得到比较好的精度，它不失为目前DEM生产的重要方法之一。该种方法的生产工艺如下所示：
格网DEM的缺点

1）在地形平坦的地方，存在大量的数据冗余； 2）在不改变格网大小的情况下，难以表达复杂地形的突变现象；
2.等值线模型
3.不规则三角网DEM
不规则三角网（TIN）
1 2 3 4 5 6 7 8 X X X X X X X X Y Y Y Y Y Y Y Y Z Z
3 1 2 顶点邻接三角形
1
4
1 1 1 2 5 4 4 3
5 4 2 3 6 5 7 4
6 5 4 4 8 8 8 7
2 1
5 3
X 6 2
Z Z Z Z Z Z
6 5 1
2 5 4 8 7
2
3
3 4
X 4 3 1 2 6 4
X 8 X 6 5 8 7 7 X X
6
5 6 7
8
7
8
点文件
三角网的一种存储方式
三角形文件

第二种是点栅格观点：认为该网格单元的数值是网格中心点的高程或该网格单元的平均高程值，这样就需要用一种插值方法来计算每个点的高程。

基于加权最小二乘曲面拟合的规则格网DEM建立

本文使用方程为 Z = x exp( - 10- 4 x2 - 10- 4 y2 ) 的曲面来模拟地表。因为该曲面既有地势变化较快区域 ( 山峰、山谷 ) , 也有地势变化平缓的区域 ( 缓坡 ) 。有利于检验算法的实用性。
图 4 加权拟合曲面
44
海洋测绘
第 28卷
由图 4可以看出, 曲面拟合度较好, 地形起伏较大和起伏平缓处的地形都可以得到较好的模拟, 地形得以逼真的再现。 4 4 有无权值拟合的精度比较
(di )-u 2
i= 1 j= 1
i= 1
i∃ j
( 11) 同理:
nn
! ! !f
( x, y ) !y
=
i= 1
(
d
i
)-
u-
2
d
u j
(
y
-
yi ) u ( z
-
j= 1
zi )
i+ 1
! / n ( di ) -u 2 i= 1
( 12)
由式 (
11)、( 12)可知当
P
在 D i邻域中,
( 8)
可得法方程组的矩阵形式:
A∀ u = B
( 9)
式中, 系数矩阵 A = [ aij ]m #m, 常数项 B = ( b1, , bm )T; 未知项 u= ( u1, , um )T。
观察可知系数矩阵为一对称矩阵。由采样数据计算各项后, 便可用高斯消去法等经典方法求解出拟合表达式的系数 u1, u2, , um。
假定要计算网格点 ( a, b )上的曲面高度 f ( a, b ), 为此需要求出一个多项式 p ( x, y ), 一般二次多项式: p (x, y ) = C00 + C10X + C 01 Y + C20X 2 + C11X Y + C02 Y2

格网DEM的几种插值方法实践

，一，、２，
｛＋（
）
）
一
１
＋１ｋｚ＋ｚ一１ｐ（）Ｉ’ （一）（４Ｚ）＋一｝，、ｌ，
：
Ｐ＝ ” Ｐ＝ｆｋ／，或、ｌ来求取距离权重。如果ａ
考虑采样点方向的影响，应以公式（）改化权则１来
作者简介：武广臣（９９一）男，１７，黑龙江人，宁科技学院辽资源与建筑工程学院讲师，硕士．
第１４卷
第ｌ期
辽宁科技学院学报
Ｖ１１Ｎ．ｏ．４ｏ１
２１年３０２月
ＪＵＮＬＯＩＯＩＧＩＳＩＴＦＳＩＮＥＡＤＴＣＮＬＧＭｒ２ｌＯＲＡＦＬＡＮＮＴＵＥＯＣＥＣＮＥＨＯＯＹＮＴａＱ２
验表明：三次插值法具有良好的连续性和较高的地形仿真度。双关键词：最邻近插值法；线性内插法；三次插值法双双
中图分类号：Ｂ２Ｔ２
文献标识码：Ａ
格网ＤＭ插值是ＧＳ空间数据插值的一个方ＥＩ
Ｐ＝（ｉ（＋ｔｉ）１）
二出！２
ｄ
．ቤተ መጻሕፍቲ ባይዱ
—＿ ■一
＝
１，有时采样点并不规则（２，时需根据矩形）但图）此或圆形对邻近区域的采样点进行插值计算，用的采
方法多为顾及地形特征的加权平均值内插，可利用

基于算术编码的格网DEM压缩技术综述

１．装备学院重点实验室，北京１１１０４６２装备指挥技术学院研究生管理大队，．北京１１１０４６
１ｈｅａｆｈＡｃｄｍｙｏＥｕｐｎ，ｅｉ０４６Ｃｉａ．ｅｙＬｂｏｔｅａｅｆｑｉｍｅｔＢｉｎ１１１，ｈｎＴＫｊｇ２Ｃｍｐｎｆｏｔｒｄａｅｎｇｍｅｔｔｅｃｄｍｙｆｑｉｍｅｔｏｍａｄ＆Ｔｃｎｌｇ，ｅｉｇ１１１，ｈｎ．ｏａｙｏＰｓａｕｔａｅｎ，ｈａｅＥｕｐｎｍｎｇＭａＡｏＣｅｈｏｏｙＢｉｎ０４６Ｃｉａｊ
ＣｍｕｒｎｉｅｉｄｐｌａｏｓｏｐｔｇｅｒｇｎＡｐｉｔＭ压缩技术综述Ｅ
邹鹏卞燕山，觳，曹
ＺＯＵｎｇ．ＡＮｎｈｎ．ＰｅＢＩＹａｓａ２ＣＡＯ。Ｒｕｉ
１引言
数字高程模型（Ｍ）ＤＥ是三维地形可视化系统最重要的数据源之一，是进行三维空间数据处理、地形特征分析的核心数据，用格网（ｉ）者不规则三角网（Ｉ来表示ＤＭ数常Ｇｒ或ｄＴＮ）Ｅ
ＺＯＵｅｇＢＩＰｎ，ＡＮａｓａ，ＹｎｈｎＣＡＯｉＳｒｅｆＧｒｄＤＥＭｏＲｕ．ｕｖｙｏｉｃｍｐｅｓｎｂｓｄｏｒｔｒｓｉａｅｎａｉｏｈｍｅｉｏｉｇＣｏｕｅｇｎｅｉｇａｄｔｃｄｎ．ｍｐｔｒＥｎｉｅｒｎｎｃ
拉格朗日预测和最小二乘预￣ｓ类，三类算法进行了对比分析。指出了算术编码算法在实际运行中存在的问题，其未来，－对这ｌ对

3 格网DEM的建立

3 格网DEM的建立
3 格网DEM的建立
31
从散点到统计表面描述
2
规则格网DEM建立的思路
3
格网DEM内插数学模型
4
格网DEM 建立过程
3.1 从散点到统计表面
❖DEM数学机理
地形起伏 H 看作是平面位置(x, y)的函数即
H f (x, y)
是一个单值光滑连续的函数； ❖地形变化复杂多样，地形表面存在的各种
（3）逐点内插计算简单，比较灵活，是较为常用的一类DEM内插方法。逐点内插方法的主要问题是内插点邻域的确定，它不仅影响到DEM内插精度，也影响到内插速度。
（4）通过建立TIN进而生成格网DEM的方法是一种被人们广为接受的DEM 建立方案，主要原因是TIN能完全覆盖研究区域表面，可适应各种分布数据，并能方便的处理各种地形结构线、断裂线。
的内插计内算插范；围
局部内插方法逐点内插方法
内插曲面与
纯二维内插
❖DEM概念提参考出点关至系今，经曲面过拟合内多插年的发展和完善，
已经提出多种高程内插方法线性插值
❖DEDMEM内内插插的根本是对地多项形式内插曲面特双征线性的插值认识和理解，具内体插函到数性方质法上，则是内插点邻域范
围的确定、权值确定方法（自相高次关多项程式插度值）、
3.4 格网DEM建立过程
主要讨论以下三个方面的内容： ❖3.4.1 基于不规则采样点的格网DEM建立 ❖3.4.2 基于规则采样点格网DEM建立 ❖3.4.3 基于等高线分布采样点的格网DEM建
立
3.4.1 基于不规则采样点的格网DEM建立
❖质基上于是不数规据则规分则布不化采规采处样样则理点点分布过的程格，网DEM建立：实
3.4.1 基于不规则采样点的格网DEM建立

基于DEM格网的改进型A＊路径搜索算法

ｏｔｅａｇｒｔｍｒａｉａｅｙｐａｔｃｌｐｌａｉｎｆｈｌｏｉｈａｅｖｌｔｄｂｒｃｉａｐｉｔｏ．ｄａｃ
Ｋｅｏｄ：ＤＥＧｒ；ｒａｐ —ｅｏｋＧＩ；ｐｔ— ｎｉｇｙｗｒｓＭ— ｉｏｄｏｏｎｔｒ；Ｓａｈｆｄｎ；Ｄｉｓａａｏｉｍ；Ａｌｏｉｍｄｔｗｉｊｔｌｒｈｋｒｇｔａｇｒｈｔ
ＡｂｔａｔＴｅｌｔａｈＦｉｄｎｒｂｅｂｓｄｏＭ — ｉ，ｔｅｄｆｅｅｃｎａｎｓｅｗｅｎＤＥＭ — ｉｎａｐ — ｅ — ｓｒｃ：ｏｄａｈＰｔ — ｎｉｇｐｏｌｍａｅｎＤＥＧｒｗｉｄｈｉｒｎｅａｄｓｍｅｅｓｔｅｆｂＧｒｄａｄｒｄｔｏｎｔｏｏ
ａｇｒｔｍｏｓｄｒｇｔｅｇｏｒｐｙｓａｅｃａａｔｒｓｍｐｏｅ，ａｄａｇｅｄｌｍｐｓｄＦｎｌ，ｔｅｆａｉｉｔｎａｉｉｌｏｉｃｎｉｅｎｅｇａｈｐｃｈｒｃｅｓｉｅｌｙｄｈｉｈｎｒｅｙｒｅｉｉｏｅ．ｉａｌｕｓｙｈｓｌａｄｖｌｔｅｂｉｙｄｙ
摘要：解决ＤＭ格网上的路径搜索问题，析了ＤＭ格网和道路拓扑网络的共性与区别，结合ＤＭ格网的特点对比为Ｅ分Ｅ并Ｅ
了Ｄｉｓａｊｔ算法和Ａ算法的优缺点，出了基于ＤＭ格网的改进型Ａ路径搜索算法。该算法充分考虑ＤＭ格网上的路径ｋｒ提ＥＥ可达条件，用移动窗口法获取搜索网格候选集，并根据地理空间分布特征选取Ａ算法的估价函数，同时引入贪婪准则缩采

DEM格网——精选推荐

DEM格⽹绘制规则格⽹和左右格⽹对⾓线⼀．绘制规则格⽹：x=ones(3,3)for i=1:3x(:,i)=i;endxy=1:3;plot(x,y,'R*-')hold on;plot(y,x,'R*-')⼆．绘制左对⾓线：m=[1,2,3];plot(m,n,'R*-')n=[2,3];m=[1,2];plot(m,n,'R*-')m=[2,3];n=[1,2];plot(m,n,'R*-')三．绘制右对⾓线：m=[1,2,3];n=[3,2,1];plot(m,n,'R*-')n=[2,1];plot(m,n,'R*-')m=[3,2];n=[2,3];plot(m,n,'R*-')四．同时绘制左右对⾓线：x=ones(3,3)for i=1:3x(:,i)=i;endy=1:3;plot(x,y,'R*-')plot(y,x,'R*-') %左斜线m=[1,2,3];n=[1,2,3];plot(m,n,'R*-') n=[2,3];m=[1,2];plot(m,n,'R*-') m=[2,3];n=[1,2];plot(m,n,'R*-') %右斜线m=[1,2,3];n=[3,2,1];plot(m,n,'R*-') n=[2,1];m=[1,2];plot(m,n,'R*-') m=[3,2];n=[2,3];plot(m,n,'R*-')Delaunay三⾓形算法实现⼀．理论依据：1.找出离散点集中相距最近的两点，连接这两点形成TIN的初始基线；2.沿基线的固定⼀侧搜寻第三点，⽣成第⼀个Delauney三⾓形；3.以三⾓形的两条新边作为新的基线；4.依次重复2,3直⾄所有的基线处理完毕；所以在算法中分别⽤了⼏个函数分别是⽤来获取第三个点，根据余弦值来判断三⾓形的，构建三⾓⽹的。