测绘工程实习任务1

《测绘工程实习任务1》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、摄影测量实习报告DEM内插法编程实验报告一、实习目的通过本次实习，使学生进一步理解移动曲面法数字高程模型内插原理，掌握移动曲面法数字高程模型内插原理及其内插子程序的设计方法，了解其它逐点高程内插方法的基木原理。掌握移动曲面法数字高程模型内插原理及其内插子程序的设计方法，了解其它逐点高程内插方法的基本原理。二、实习内容根据提供的io个数据点的坐标（x,y,z）和待求点的平面坐标（x,y）,要求利用移nnnpp动二次曲而拟合法，由格网点P（X,Y）周圉的10个已知点内插出待求格网点P的高程，pp编制相应的程序并进行调试，最后解算

2、出格网点P的高程并提交源程序代码。三、实习准备资料C知数据点坐标点号XYZ1102110152109113183105115194103103175108105216105108157115104208118108159116113171011311822编程计算点（110,110）上的高程。四、程序的代码#includeHstdafx.hH#include"SMatrix.h1'#include,,SingleImageResection.h,'intmain(intargc,char*argv[])SMatrixX(10

3、,l);SMatrixY(10,l);SMatrixZ(10,l);X[0][0]=102;X[l][0]=109;X⑵⑹=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[1JIO]=113;Y[2][0]=115;Y[3][0]=103;Y[4][0]=105;Y[5][0]=10&Y16J10]=104;Y[7][0]=10&Y[8][0]=113;Y[9][0]=1

4、1&系数矩阵系数矩阵系数矩阵输入已知数据点坐标输入己知数据点坐标输入已知数据点坐标输入已知数据点坐标输入已知数据点坐标输入已知数据点坐标输入己知数据点坐标输入已知数据点坐标输入已知数据点坐标输入已知数据点坐标输入己知数据点坐标输入己知数据点坐标输入已知数据点坐标输入已知数据点坐标输入已知数据点坐标输入己知数据点坐标输入已知数据点坐标输入已知数据点坐标输入已知数据点坐标输入己知数据点坐标ZLO][OJ=15;输入已知数据点坐标Z[l][0]=18;输入已知数据点坐标Z[2][0]=19;输入已知数据点坐标Z⑶[0]=17;输入

5、C知数据点坐标Z[4][0]=21;输入己知数据点坐标Z⑸[0J=15;输入已知数据点坐标Z[6][0]=20;输入已知数据点坐标Z[7][0]=15;输入已知数据点坐标Z[8][0]=17;输入C知数据点坐标Z[9][0]=22;输入己知数据点坐标for(inti=O;i<10;i++){xri][o]=xriiro]-no；Y[i][0]=Y[i][0]-110;}SMatrixM(10,6);for(i=0;i<10;i++){M[i][O]=X[i][O]*X[i][O];M[i][l]=X[i][O]*Y[i][O

6、];M[i][2]=Y[i][O]*Y[i][O]；M[i]⑶=X[i][O];M[i][4]=Y[i][O];⑸二1;SMatrixP(10,10);for(i=0;i<10;i++)P[i][i]=l/(X[i][0]*X[i][0]+Y[i][0]*Y[i][0]);x=(M.T()*P*M).Invert()*M.T()*P*Z;printfC*待定点的高程是：％f“,x⑸[0]);return0;HSMatrix.hHuSingleImageResection.hH2、输入己知点数据的坐标3、选择二次曲面作为拟

7、合曲面：Z=Ax2^Bxy^Cy-^-Dx+Ey^F列出误差方程:叫=戸:A+XXB+X'C+yD+EE+F-Z4、组方程，求出6个系数:5、计算出待求点高程。六、实验心得体会由于这次是实验是要用C++进行DEM内插法的编程，了解DEM内插法的特点，由于以前没写过类似的程序，感觉有些困难，但其编写有一定的规律特点，加上以前学习过c语言编程，所以仔细研究后还是能了解其代码的含义，程序编写过程中还要不断地调试。积极发现错误，细细改正，最后得出想要的结果。通过这次编程学习，主要培养的发现问题、思考问题和解决问题的能力，还要不断地学

8、习来完善自己。只有认真努力去做，才能把事情做好，编程更是如此，还要细心。