用Matlab 实现移动曲面拟合法生成DEM
标签： Matlab 移动曲面拟合法 DEM 分类： [E]【 MATLAB 】2006-12-22 10:02
用Matlab 实现移动曲面拟合法生成DEM
杜玉军
（武汉大学测绘工程0408班 200431610007 武汉 430079）
摘要：移动曲面拟合法是DEM格网点内插常用的一种方法，利用Matlab可以轻松实现该方法生成DEM。
关键字：移动曲面拟合法 DEM Matlab
1． 概述
为了获取规则格网DEM，内插是必不可少的过程。内插的方法很多，其中移动曲面拟合法由于其方法灵活、计算简便、精度较高、占用内存较少等诸多优点而经常被使用。
2． 实现原理
移动曲面拟合法是一种以待定点为中心的逐点内插法，它以每个待定点为中心，定义一个局部函数去拟合周围的数据点。其过程为：
（1） 对每个格网点，从数据点中检索出邻近的个（至少6个）数据点。
以待定点（）为圆心，以选定长为半径作圆，凡落入圆内的数据点都被采用。
Xpi=Xi-XYpi=Yi-Y
di2= Xpi2+Ypi2
di
标准
excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载
对话框
fid1=fopen(strcat(pathname,filename),'rt'); %以只读形式打开
if fid1==-1 %若没有选择文件则警告
msgbox('Input File or Path is not correct','Warning','warn');
break;
end
Dt=load(filename); %获取数据
fclose(fid1);
·函数文件2
function zp=GridZ(pt,x,y,s0)
%移动曲面拟合法内插(x,y)处的高程
%pt为数据点矩阵，s0为单点平均占用面积
N=length(pt); %点数
n0=8; %搜索点数
%初始化各值
n=1; m=1; %计数器
xp=0; yp=0; d2=0; %原点在(x,y)时数据点坐标和与(x,y)距离的平方
ptin.x=0; ptin.y=0; ptin.z=0; din2=0; %落入搜索范围内的数据点坐标值和与(x,y)距离的平方
ptout.x=0; ptout.y=0; ptout.z=0; dout2=0; %搜索区外的数据点坐标和距离
P=0; %权阵
for k=1:N
xp(k)=pt(k,1)-x; yp(k)=pt(k,2)-y;
d2(k)=xp(k)^2+yp(k)^2;
if d2(k)
本文档为【matlab离散点生成dem】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483