MATLAB插值在测绘中的应用.docVIP

《现代测绘数据处理方法》课程 课间实验报告 MATLAB插值在测绘中的应用 姓 名： 学 号 班 级： 专业 级 班 指导教师： 实验日期：20XX年09月20日 测量与空间信息处理实验 一、实验目的及所用软件版本 1、实验目的 （1）掌握插值的基本原理； （2）利用Matlab进行地形数据的插值运算； （3）将插值结果以可视化的方式进行显示； （4）不同插值方法的对比。 2、实验软件所用版本 实验软件Matlab7.0+Window7 二、实验内容及问题背景 1、实验内容 基于matlab的插值在测绘中的应用 2、实验内容所涉及的问题背景 在某山区测得一些地点的高程如下表所示，平面区域1200≤≤4000，1200≤≤3600。试画出该山区的地貌图和等高线图，并对几种插值方法进行比较。 坐标 高程 1200 1600 2021 2400 2800 3200 3600 4000 坐 标 3600 1480 1500 1550 1510 1430 1300 1200 980 3200 1500 1550 1600 1550 1600 1600 1600 1550 2800 1550 1200 1100 1550 1600 1550 1380 1070 2400 1550 1200 1100 1350 1450 1200 1150 1010 2021 1390 1500 1500 1400 900 1100 1060 950 1600 1320 1450 1420 1400 1300 700 900 850 1200 1130 1250 1280 1230 1040 900 500 700 三、实验原理与数学模型 四、实验过程记录（含基本步骤、程序清单及异常情况记录等） 第1步，先利用原始三维坐标，画出粗糙的地表高程分布曲面图和等高线图。 （1）根据原始数据画出地貌曲面图 %x代表横坐标 x=1200:400:4000; %y代表纵坐标 y=3600:-400:1200; %Z是已知坐标点的高程值 Z=[1480 1500 1550 1510 1430 1300 1200 980; 1500 1550 1600 1550 1600 1600 1600 1550; 1550 1200 1100 1550 1600 1550 1380 1070; 1550 1200 1100 1350 1450 1200 1150 1010; 1390 1500 1500 1400 900 1100 1060 950; 1320 1450 1420 1400 1300 700 900 850; 1130 1250 1280 1230 1040 900 500 700]; mesh(x,y,Z); xlabel(x); ylabel(y); zlabel(Z:山的高度); 结果如下： 图1 原始地貌曲面图 备注：注意平面位置（）和对应关系。 （2）根据原始地貌曲面图画出原始地貌等高线图。 在继续添加，利用少数已知高程点勾绘出等高线，结果如下图所示： figure; [c,h]=contour(x,y,Z,[100:100:1600]); xlabel(x); ylabel(y); clabel(c,h); 结果如下： 图2 原始地貌等高线图 （3）原始地貌曲面图和等高线图对比及叠加。 将原始地貌曲面图和原始地貌等高线图放在一起，可以发现非常粗糙，如图3所示。 原始地貌曲面图 原始地貌等高线图 图3 原始地貌曲面图和等高线图 利用meshc函数将原始地貌曲面图和等高线图叠加如图4所示。 figure meshc(x,y,Z); 图4 原始地貌曲面图和等高线图叠加 第2步，将原始地貌数据插值，并画出曲面图和等高线图 （1）将原始地貌数据插值，得插值后地貌曲面图，如图5所示。 [X,Y]=meshgrid(x,y); [Xi,Yi]=meshgrid(1200:50:4000,3600:-50:1200); Zi=interp2(X,Y,Z,Xi,Yi,cubic); figure; mesh(Xi,Yi,Zi); xlabel(X); ylabel(Y); zlabel(山的高度); 图5 原始地貌插值后曲面图 （2）根据插值后地貌曲面图绘制等高线图，如图6所示。 figure; [c,h]=contour(Xi,Yi,Zi,[100:100:1600]); clabel(c,h); xlabel(X);