空间分析-仪器信息网.ppt

第十四讲 空间分析与查询（五） 中山大学 遥感与地理信息工程系 数字高程模型 数字地面模型(Digital Terrain Model, 简称DTM)是对地理空间起伏变化的连续表面的数字表示形式。 一般情况下，被描述的连续表面是地形面，即高程值在地理空间上的变化，称为数字高程模型(DEM, Digital Elevation Model)。被描述的连续表面也可以是地理空间上的地价、污染负荷量、绿化率、降雨量、气温、人口密度、建筑物密度等等。 介绍两种数字高程模型，即高程矩阵(GRID)与不规则三角网(TIN) 。 高程矩阵(GRID)的生成方法 空间内插是通过使用现有数据点的变量值估计一个非样点位置的变量值的方法。 插值的方法有很多种，其中最重要的有线性内插法、双线性插值、移动拟合法、趋势面插值、样条插值及克里金(Kriging)插值，这些方法有的是整体拟合，有的是局部拟合。 线性内插法 用线性曲面方程Zp=a0+a1x+a2y拟合待定点附近的地形面（例如TIN模型中的三角平面可视为线性曲面）。可用3个数据点求解待定系数，然后根据线性曲面方程计算待定点的高程值。 双线性插值 用曲面方程Zp=a0+a1x+a2y+a3xy拟合待定点附近的地形面，也可以表示为 Zp=(ax+b)(cy+d)。 即认为待定点附近的地形面的高程Zp在x轴方向呈线性变化，同时在y轴方向上也呈线性变化。可用待定点附近的4个数据点来计算双线性曲面函数的待定系数，再根据曲面函数计算待定点的高程值。 移动拟合法 （1） 移动拟合法是用二次多项式来拟合地面高程，即Zp=Ax2+Bxy+Cy2+Dx+Ey+F，通常把坐标原点平移到待定点之后利用上式，采用待定点为圆心的圆周内6个数据点求解函数的待定系数，然后根据曲面方程计算Zp 移动拟合法（2） 当数据点多于6个时，可按最小二乘法求解待定系数。根据数据点距待定点远近不同，利用“距离越远对插值点的影响越小”的思想，可赋以适当的权值，权重可采用W=1/di2或W=[(R-di)/di]2或W=exp(-di2/R2)等形式（其中di为距离，R为采样圆半径），该方法称为按距离加权最小二乘内插法。 移动拟合法（3） 首先建立误差方程式： 再根据最小二乘法原理，按vwvT=min的方法，建立方程式求解待定系数。 移动拟合法的结果显示误差分布和地形坡度大小没有特别显著的关系，在地性线地区无法显示出这些转折地地形。 趋势面(trend)方法 趋势面分析是一种整体内插法，这种方法假设地形表面的一般趋势与高程值的局部变化无关，并利用曲面方程来模拟待定点附近地形表面的一般趋势。 通常使用的是一次、二次、三次趋势面，过高次的趋势面不利于反映空间趋势，并可能存在趋势面的“畸变”。其中，二次趋势面可用待定点附近的6个数据点来计算方程式系数。若取用的点数多于6个，则要用测量平差方法来拟合。 趋势面(trend)方法产生的DEM难以反映地形起伏的特性。 双三次多项式（样条函数）插值（1） 即用双三次多项式拟合地形表面。 双三次多项式（样条函数）插值（2） 该曲面模型有16个待定系数，需要方程组成的方程组来求解。通常用4个数据点（规则格网的4个顶点）的4个函数值组成的4x4方程组求解，见图4-11。这4个函数值是高程Z、x方向斜率R、y方向斜率S、及xy方向扭矩T： 双三次多项式（样条函数）插值（3） 其中Z保证曲面通过格网的4个数据点，R、S、T保证曲面在这4个数据点光滑连续。 双三次多项式内插法是规则格网插密的常用方法之一。这种方法通过一系列曲面片段来拼接地形表面，最终得到一个一阶、二阶连续的表面。该方法属于局部插值，计算负担中等；对于平滑表面拟合效果最好，对于起伏的表面拟合效果最差。 克里金(Kriging)插值 Kriging插值方法在空间相关范围分析的基础上，用相关范围内的采样点来估计待插点属性值。 （1）数据检验与分析，删去明显偏离实际的采样数据点。 （2）数据预处理。 （3）绘制方差图，了解空间变量的集聚范围与方向。 （4）克里金插值估计。 相对来说，克里金插值则能较好地反映各种地形变化，但克里金方法的计算量很大，因此在对大面积区域大数据量内插时，这是一个不能不考虑的因素。? DEM应用举例 1、地形图－景观描述，地形如坡度、方向等的属性。 2、市政工程－土木工程、建筑和道路规化。 3、规划－视觉影响分析、建筑以及其它物体的定位。 4、测深学/水道学－近岸建筑物的设计，小比例尺图表 5、采矿工程－亚表面模拟和抽取规划。 6、地质学上的模型－亚表面模型，断层模型。