gis-lecture6-spatialgrid.ppt

VI 空间数据分析 6.1 空间数据的量算 6.2 空间数据的叠置分析 6.3 空间数据的邻域分析 6.4 TIN与泰森多边形 6.5 数字高程模型(DEM) 6.6 栅格GIS几何代数语言 6.7 优化问题 – 最短路径 6.5数字高程模型(DEM) DTM模型(Digital Terrain Model)数字地形模型 50年代由MIT摄影测量实验室提出，是用数字形式描述地形表面的模型。实质上这是对地面形态和属性信息的数字表达。 DEM模型(Digital Elevation Model)数字高程模型 当DTM模型中数字属性为高程时称DEM模型, 即数字高程模型。DEM模型是DTM模型的一种特例。 格网数据结构 格网数据结构是典型的栅格数据结构，可采用栅格矩阵及其压缩编码的方法表示。通常，其数据包括三部分： 1）元数据 描述DEM数据的数据，如数据表示的时间，边界，测量单位，投影参数. 2）数据头 DEM数据的起点坐标，坐标类型，格网大小，行列数等. 3）数据体 行列数分布的数据阵列。 6.5.1 DEM数据的获取 1）航空航天遥感数据源 航空照片、卫星影像通过立体象对，用摄影测量的方法建立空间立体模型量取密集的高程数据。 激光断面扫描仪(Laser Scanner)直接用于测定地面高程，从而建立数字高程模型。它向地面发射高频激光波束并接受反射波，精确地记录波束传输时间。传感器的位置和姿态参数由GPS和IMU (Inertial Measuring Units). 合成孔径雷达（SAR）能够全天候、全天时获取大范围、高分辨率图像，利用干涉SAR（InSAR）的精确相位信息可以提取地表的三维信息即生成数字高程图（DEM）和检测高精度的地表高度变化。 2）以地形图为数据源 主要用比例尺不大于己于1:1万的国家近期地形图为数据源从中量取等密度地面点集的高程数据，建立DEM。数据采集通常用手工、数字化仪几及扫描仪。 3）以地面实測记录为数据源 对小范围的大比例尺（如1:5000）的DEM数据可用GPS，电子測速仪（全站仪），測距经纬仪等获取数据。 6.5.2 数据的插值方法 1）网格尺寸的确定 网格尺寸的确定是DEM数据网格化的重要问题，它关系到数据模型的精度。 网格尺寸的确定首先应分析地形的形态特征。 一般地说，采样点的密度基本上确定了网格点的密度。在采集优选点的情况下，网格点数应大于或接近于样点数，取 n N 2n， N 为网格点数，n为样点数。 此外，DEM的网格大小还可以根据地图比例尺和人眼分辨率/地图纸张印刷分辨率来初步计算。如何算？例如：1:5万地形图+人眼分辨率300dpi,DEM网格大小？ 2) 内插方法 数据插值的核心问题是确定内插方法，任何一种内插方法都是基于原始函数光滑性，也就是说认为邻近数据点之间存在很大的相关性。 插值的方法很多，如按距离加权法、多项式内插法、样条函数内插法、多面函数法等等。大量的实验证明，由于实际地形的非平稳性，不同的内插方法对DEM的精度并无显著影响，主要取决于原始采样点的密度和分布。简单而常用的为线性内插法和双线性多项式内插法。 DEM常用的插值算法(ArcView/ArcInfo)： 1）距离加权法 距离愈远的点对待插点的影响愈小； 2）趋势面插值法 用通过采样点的空间曲面来模拟地形表面，常用二次多项式来拟合 3）样条插值法 用距离函数样条法，或分片样条法进行插值 4）克吕格插值法 最优线性无偏内插估计量插值法考虑采样点分布，但计算量大。 IDW Z0 = ?[zi (1/dik)] /? (1/dik) Trend Surface (arcinfo 最高12次多项式) Zx,y = b0 +b1x + b2y Zx,y = b0 + b1x + b2y + b3x2 + b4xy + b5y2 Zx,y = b0 + b1x + b2y + b3x2 + b4xy + b5y2 + b6x3 + b7x2y + b8xy2 + b9y3 Splines Q(x,y) = ?(Aidi2logdi) + a + bx + cy di2 = (x-xi) 2 +(y-yi) 2 – the distance between the point to be interpolated and the control points a + bx + cy - local trend function di2logdi – the minimum curvature surface function ?(Aidi2logdi) + a + b x