数字高程模型(DEM)内插程序设计总结

数字高程模型（DEM）内插程序设计总结测绘08-2廖小军3080208220【摘要】：数字高程模型(DEM)是构建虚拟地形环境的重要步骤之一,DEM的精度不仅影响地形可视化的效果,而且更重要的是制约着地形仿真和GIS分析与决策的可信度。随着现代测绘、GIS和VR等技术的发展,地形建模的应用范围越来越广阔。本次实习主要内容为在分析移动曲面拟合与加权平均算法特点的基础上,提出以移动曲面拟合法为主、加权平均法为辅两者相结合的方法,使它们能够相互取长补短,从而解决大区域DEM建模中存在的问题。【关键词】：数字高程模型移动曲面拟合加权平均DEM内插1.实验目的。掌握移动曲面法数字高程模型内插原理及其内插子程序的设计方法，了解其它逐点高程内插方法的基本原理。2.实验内容。根据提供的10个数据点的坐标（Xn,Yn,Zn）和待求点的平面坐标（Xp,Yp），要求利用移动二次曲面拟合法，由格网点P（Xp,Yp）周围的10个已知点内插出待求格网点P的高程，编制相应的程序并进行调试，最后解算出格网点P的高程并提交源程序代码。3.资料准备。已知数据点坐标编程计算点（110,110）上的高程。4.基本思路。5.源代码。//tggfhgfh.cpp:Definestheentrypointfortheconsoleapplication.//#includestdafx.h#includeSMatrix.hintmain(intargc,char*argv[]){printf(HelloWorld!\n);SMatrixa(6,1);//未知数（A,B,C,D,E,F）SMatrixx(10,1);//光标X坐标，已知SMatrixy(10,1);//光标y坐标，已知SMatrixz(10,1);//光标Z坐标，已知SMatrixM(10,6);//系数矩阵SMatrixP(10,10);//权阵SMatrixMTM(6,6);x[0][0]=102;x[1][0]=109;x[2][0]=105;x[3][0]=103;x[4][0]=108;x[5][0]=105;x[6][0]=115;x[7][0]=118;x[8][0]=116;x[9][0]=113;y[0][0]=110;y[1][0]=113;y[2][0]=115;y[3][0]=103;y[4][0]=105;y[5][0]=108;y[6][0]=104;y[7][0]=108;y[8][0]=113;y[9][0]=118;z[0][0]=15;z[1][0]=18;z[2][0]=19;z[3][0]=17;z[4][0]=21;z[5][0]=15;z[6][0]=20;z[7][0]=15;z[8][0]=17;z[9][0]=22;for(inti=0;i10;i++){M[i][0]=(x[i][0]-110)*(x[i][0]-110);M[i][1]=(x[i][0]-110)*(y[i][0]-110);M[i][2]=(y[i][0]-110)*(y[i][0]-110);M[i][3]=x[i][0]-110;M[i][4]=y[i][0]-110;M[i][5]=1;P[i][i]=1/((x[i][0]-110)*(x[i][0]-110)+(y[i][0]-110)*(y[i][0]-110));}MTM=M.T()*P*M;a=MTM.Invert()*M.T()*P*z;printf(a[0]:%f,,a[0][0]);printf(a[1]:%f,,a[1][0]);printf(a[2]:%f,,a[2][0]);printf(a[3]:%f,,a[3][0]);printf(a[4]:%f,,a[4][0]);printf(a[5]:%f\n,a[5][0]);return0;}6.计算结果。运行程序得待定点高程为：17.9701807.心得体会。通过这次上机实习学到了许多知识，特别是有关矩阵的运算方面的知识。虽然以前也学过了矩阵的运算，但由于运用较少，已经几乎忘了。孔子说过“温故而知新，可以为师矣。”正是这样一遍一遍的练习才会让我们学到更多的东西，每一次操作都有一定的收获。在这一次实习不仅学到的知识增加了，而且更懂得了怎样去学习才可以获得更多知识。在上机课上有什么不懂得地方可以问同学和老师，这是一个学习好的资源，只有不断的学习才可以使自己了解得更多对自己也是非常有益，因此，只有把握好机会就一定可以获得更多。机会也是自己把握得，错过了也就很难再找到。只有做好目前的任务也就可以更好地为以后服务，做好这一步对以后学习更加省力轻松。因此，在这次实习中不仅学的更多，而且也为以后地学习打下了坚实的基础。再次谢谢老师的辛勤指导。