大地电磁测深实习报告.doc

精品学习资料范文 大地电磁测深实习报告 篇一：大地测量上机实习报告 题目： 大地测量学上机实习学 号： 姓 名： 班 级： 专 业： 测绘工程 课程名称： 大地测量学基础课程老师： 大地测量学上机实习报告 一、实习目的： 了解大地测量学的定义以及作用，让同学们充分了解大地测量学在现代社会的重要地位及其广泛的应用，让同学们体验大地测量的一些解题计算方法与运用。也让同学们对于复杂的测绘计算如何用计算机通过编程来解决,同时对于坐标转化和高斯投影的处理方法进行一定的了解。 二、实习任务： 1、 2、 空间坐标转换编程； 高斯投影正算编程； 3、 高斯投影反算编程； 三、实习时间与地点： 地点：核工楼410教室； 时间：每周四上午一二节课，共三周； 四、实习仪器： 电脑，C++编程软件； 5、实习步骤： 5.1实验方法步骤 1、先了解书本上关于大地主题解算的步骤； 2、根据公式用VC编写程序 3、对编写的程序代码进行调试运行 5.1.1大地坐标空间直角坐标系转换 X?(N?H)cosBcosL? 1）空间大地坐标系向空间直角坐标系的转换公式为：Y?(N?H)cosBsinL?? Z?[N(1?e2)?H]sinB?? 式中，N?a，a为椭球的长半轴，N为椭球的卯酉圈曲率半径a=6378.137km W W??e2sin2B，e2?a2?b2 a2，e为椭球的第一偏心率， b=6356.7523141km。 2）程序代码 #include stdio.h #include math.h #define PI 3.1415926535897932384626433832795 double rad(double a,double b,double c) {b=b+c/60; a=a+b/60; a=a/180*PI; return a;} void deg(double d,int *d1,int *d2,double *d3) {d=d/PI; *d1=(int)d; *d2=(int)((d-*d1)*60); *d3=((d-*d1)*60-*d2)*60;} void X1(double a1,double e) {double l,b,h,x,y,z,a,b1,c,n,w; printf( please input L:\n scanf( %lf%lf%lf , a, b1, l=rad(a,b1,c); b为椭球的短半轴 printf( please input B:\n scanf( %lf%lf%lf , a, b1, b=rad(a,b1,c); printf( please input H:\n scanf( %lf , w=sqrt((1-e*e*sin(b)*sin(b))); n=a1/w; x=(n+h)*cos(b)*cos(l); y=(n+h)*cos(b)*sin(l); z=(n*(1-e*e)+h)*sin(b); printf( X=%f\nY=%f\nZ=%f\n ,x,y,z);} void X2(double a1,double e) {double l,b,h,x,y,z,n,w,m,g,m1,g1,l3,b3;int l1,l2,b1,b2; printf( please input X:\n scanf( %lf , printf( please input Y:\n scanf( %lf , printf( please input Z:\n scanf( %lf , l=atan2(y,x); do{m=(z+a1*e*e*g/sqrt(1+g*g-e*e*g*g))/sqrt(x*x+y*y); m1=atan(m); g1=atan(g); m1=m1/PI*180*3600; g1=g1/PI*180*3600; g=m;} while(sqrt((m1-g1)*(m1-g1)) 0.0001); b=atan(m); w=sqrt(1-e*e