matlab曲面拟合的算法,用Matlab 实现移动曲面拟合法生成DEM

用Matlab实现移动曲面拟合法生成DEM

杜玉军

(武汉大学测绘工程0408班 200431610007  武汉  430079)

摘要：移动曲面拟合法是DEM格网点内插常用的一种方法，利用Matlab可以轻松实现该方法生成DEM。

关键字：移动曲面拟合法  DEM  Matlab

1．概述

为了获取规则格网DEM，内插是必不可少的过程。内插的方法很多，其中移动曲面拟合法由于其方法灵活、计算简便、精度较高、占用内存较少等诸多优点而经常被使用。

2．实现原理

移动曲面拟合法是一种以待定点为中心的逐点内插法，它以每个待定点为中心，定义一个局部函数去拟合周围的数据点。其过程为：

(1)           对每个格网点，从数据点中检索出邻近的n个(至少6个)数据点。

以待定点(X,Y)为圆心，以选定长R为半径作圆，凡落入圆内的数据点都被采用。

Xpi=Xi-X；Ypi=Yi-Y

di2= Xpi2+Ypi2

di

(2)           列立误差方程式。

选择二次曲面Z=Ax2+Bxy+Cy2+Dx+Ey+F为拟合面，则数据点pi对应的误差方程式为

vi=Xpi2A+XpiYpiB+Ypi2C+XpiD+YpiE+F-Zi

n个数据点列出的误差方程可写为：

v=MX-Z  X=[A B C…F]T

(3)           计算每一数据点的权。

本文选取Pi=1/di2定权