matlab求坡向,基于杨赤中插值法的坡度、坡向计算的实现与分析

第 28卷第 1期 2008年 1月 海 　洋 　测 　绘 HYDROGRA PH IC SURV EYING AND CHARTING V ol28, N o11 Jan. , 2008 收稿日期 : 20070424; 修回日期 : 20070514 基金项目 : 广西区教育厅科研面上项目 (2006293)。 作者简介 : 夏元平 (19822) ,男 ,江西临川人 ,硕士研究生 ,主要从事 GIS的数据处理和空间分析研究。 基于杨赤中插值法的坡度、坡向计算的实现与分析 夏元平 ,陈天伟 ,张佃友 (桂林工学院 ,广西 桂林 　541004) 　　摘要 : 应用杨赤中滤波推估法 ,将地形碎部点高程观测值看作一种包含规律性和随机性变化的复合变量 ,求其二项系数加权游动平均 ,建立估值数学模型 ,求得各格网点高程估值。将离散的数据经插值计算转换为 DEM数据 ,并根据其网格顶点来分析、推算和显示网格单元的坡度和坡向。 关键词 : 杨赤中滤波推估法 ;DEM;坡度 ;坡向 中图分类号 : P224　　　文献标识码 : A　　　文章编号 : 167123044 (2008) 0100092041 　引 　言 坡度表示地表面的倾斜程度 ,而坡向反映斜坡所面对的方向 ,它们通常与确定的点有关。在地理信息系统中 ,坡度和坡向一般在特定的数字高程模型上通过计算得到。迄今为止 ,有多种插值算法可建立 DEM 网格数据。杨赤中滤波推估法作为一种插值算法 ,应用了中国古代数学杨辉三角这一民族遗产 ,几十年来为中国矿床统计几何学科发展做出了重要的贡献。有学者经过研究证明杨赤中法构造DEM能取得稳定的精度和较好的效果 [1 ] ,本文研究将其应用于构造格网 DEM数据并对坡度坡向进行计算分析和显示的实现方法。 2　杨赤中滤波推估法 211　杨赤中滤波 将高程观测值看作是一种包含规律性变化和随机性变化的复合变量 ,通过逐次加权游动平均 ,逐步消除随机性变化成分 ,显现出基本变化成分 [2 ] 。滤波实际上是加权游动平均 ,即在一定范围内取若干观测值 (按一定循序 ) 进行加权平均 ,逐个观测值向前游动 ,然后再将新的平均值当作观测值 ,再次求取加权游动平均值 ,并继续逐次地进行。有关衡量随机性、基本变化的指标主要有如下几种。 2111　综合随机方差 它反映了复合变量的综合随机性变化 , 每次滤波之后可求得综合随机方差 S2为 : S 2 =ζ2 - ξ2 (1) 式中 ,ζ2为各原始观测值的平方的平均值;ξ2为各加 权游动平均值的平方的平均值。 2112　基本方差 Dζ =ξ2 - ξ2 (2) 式中 ,ξ2为原始观测值的平均值的平方。 基本方差 Dζ反映了基本变化的特性 ,不属于对 复合变量产生随机影响的变化分量。 213　总方差 Dζ =ζ2 - ξ2 (3) 　　可见 :Dζ =S2 +Dζ 2114　广义杨赤中滤波的权系数 [3 ] 由于高程碎部点的分布一般是不规则的 ,此时符合广义杨赤中滤波的情形 ,描述杨赤中滤波的权系数的数学模型为 : p(d) = Ae- B d2 (4) 式中 , p ( d)为权系数; A、B 为待定参数; d = lj / l 0 ( lj 为抽样点与区间中心点的距离; l 0 为平均抽样间 距 )。 广义杨赤中滤波是以半径略大于 2l 0 的圆为基 本游动区间进行逐遍滤波计算的 ,利用最小二乘回归法来拟合函数式 (4)得到 A、B 参数值为 : 一维 : A = 2, B = 01693145; 二维 : A = 4, B =