大地测量实习报告(共23页)

《大地测量实习报告(共23页)》由会员分享，可在线阅读，更多相关《大地测量实习报告(共23页)（23页珍藏版）》请在装配图网上搜索。

1、课程编号： 课程性质：必修大地测量计算与实习报告学院： 测绘学院 专业： 测绘工程 地点： 三校区，珞珈山 班级： 组号： 2 姓名： 学号： 200 教师： 向 东 2010年09月8日 至 2010年09月23日 实习目的与要求大地测量计算与实习是07级测绘工和一专业学生继大地测量学基础等理论课程学习结束的集中教学实习。通过该集中教学实习，不但使学生巩固课堂上学到的理论知识，而且应使学生的实际动手能力，仪器操作能力等得到较大的锻炼和提高。实习内容1精密水准测量部分1.1踏勘测区：了解水准路线的布设和点的位置，介绍水准点位置的选择方法。1.2填写水准点之记：每人完成一个水准点点之记的绘制。1

2、.3仪器检验：每人完成水准仪i角测定成果和一对水准标尺基辅差常数的测量定各一份1.4水准观测：按二等水准测量的精度要求，每人完成一个闭合环观测的水准测量观测。每人轮流进行观测，记录，扶尺，量距，打伞等工作，并按各组记录登记。1.5外业成果的检查计算：每 对自己观测的成果进行百分之百的检查，并完成全组的外业高差与概略高程表的计算一份。1.6控制点情况2大地测量计算部分具体任务要求：每人完成实测斜距化算至高斯投影平面边长的边长改正程序编制，并计算出结果每人完成高斯投影正反算公式的计算程序编制，并计算出结果每人完成大地主题正反算公式的计算程序编制，并计算出正确结果每人完成不同平面坐标转换计算程序编制

3、，并计算出结果i角的检校仪器：N3方法：I1ABI2观测者：日期：2010.9.7标尺：A,B记录者：时间：9:20呈像：清晰稳定检查者：仪器距近标尺距离D1= 5.6 m 仪器距远标尺距离D2= 41.3 m仪器站I1I2观测次序A尺读数a1B尺读数b1A尺读数a2B尺读数b21295 272101 406496 468305 2192 220 416 709 1983 269 317 483 1604 249 249 573 246中数295 230101 347496 558305 206高差（a-b）,mm964.53956.76方法：I1ABI2,AI1I2B=(a2-b2)-(a1

4、-b1)/2= -3.885 mm方法：AI1BI2=(a2-b2)-(a1-b1)= mmi=/(D2-D1)-1.6110-5(D1+D2)= -16.54校正：a2=a2-D2/(D2-D1)=b2=b2-D1/(D2-D1)=一对标尺零点差及水准尺尺常数测定2010年9月13日 观测者 记录者 计算者测回桩号No 水准尺读数No 水准尺读数基本分划（mm）辅助分划（mm）基辅差（mm）基本分划（mm）辅助分划（mm）基辅差（mm）I13079.119144.326065.213079.439144.146064.713079.149144.326065.143079.439144.13

5、6064.883079.19144.286064.593079.529144.526065.0022659.328724.486065.162659.128724.426065.302659.208724.356065.152659.788724.756064.972659.008724.116065.112659.288724.496065.2132376.238441.426065.192376.648441.826065.182376.208441.396065.192376.728441.826065.102376.218441.606065.392376.698441.806065.

6、11平均6065.136065.05II13063.089128.116065.033063.509128.306064.803062.789127.686064.903062.999128.456065.463063.019128.386065.373062.989128.306065.3222642.678707.946065.272642.908708.116065.212642.698708.346065.652642.888708.106065.212642.698708.316065.622642.888708.796065.8932360.418425.506065.092360

7、.328425.506065.182360.398425.256064.862360.328425.806065.482360.418425.406064.992360.488425.686065.20平均6065.206065.31III13067.639132.506064.873067.709132.486064.783067.429132.186064.763067.629132.486064.863067.249132.186064.903067.309132.206064.9022647.108712.106065.002647.458712.516065.062647.51871

8、2.216064.702647.428712.416064.992647.108712.306065.202647.758712.806065.0532364.738429.506064.772365.058429.756064.702364.628429.856065.232364.928429.506064.582364.628429.886065.262364.728429.896065.17平均6064.976064.90总中数6065.106065.09一对标尺零点不等差=0.01mmi角的检校仪器：N3方法：I1ABI2观测者：日期：2010.9.13标尺：A,B记录者：时间：12

9、:50呈像：清晰稳定检查者：仪器距近标尺距离D1= 5.7 m 仪器距远标尺距离D2= 42.4 m仪器站I1I2观测次序A尺读数a1B尺读数b1A尺读数a2B尺读数b21148 778156 801142 000 150 3592 776 739 030 3403 762 761 010 3284 772 749141 943 330中数148 772156 762141 996150 339高差（a-b）,mm-79.90-83.43方法：I1ABI2,AI1I2B=(a2-b2)-(a1-b1)/2= -1.765 mm方法：AI1BI2=(a2-b2)-(a1-b1)= mmi=/(D

10、2-D1)-1.6110-5(D1+D2)= -10.69校正：a2=a2-D2/(D2-D1)=b2=b2-D1/(D2-D1)=一对标尺零点差及水准尺尺常数测定2010年9月13日 观测者 记录者 计算者测回桩号No 水准尺读数No 水准尺读数基本分划(mm)辅助分划(mm)基辅差(mm)基本分划(mm)辅助分划(mm)基辅差(mm)111502.794518.623015.831502.414517.713015.31503.124518.713015.591502.414517.783015.371502.944518.533015.591502.614517.913015.32129

11、2.614307.993015.381291.984307.223015.241292.664307.843015.181292.024307.263015.241292.344307.883015.541291.924307.423015.531152.054168.143016.091151.674167.233015.561152.234167.883015.651151.724167.213015.491152.34167.793015.491151.74166.983015.28平均1315.8934331.4873015.5931315.3824330.7473015.364211

12、449.314465.063015.751449.494464.913015.421449.254465.053015.81449.224464.93015.681449.214465.013015.81449.294464.83015.5121239.324254.823015.51239.114254.493015.381239.44254.83015.41238.854254.53015.651239.254254.733015.481238.84254.453015.6531097.894114.033016.141099.354114.873015.521097.874113.983

13、016.111099.354114.953015.61097.94114.013016.111099.234114.893015.66平均1262.1564277.9433015.7881262.5214278.0843015.563311470.334485.913015.581470.284485.553015.271470.314485.913015.61470.234485.63015.371470.294485.833015.541470.254485.633015.3821258.654274.33015.651258.164273.613015.451258.554274.213

14、015.661258.294273.623015.331258.64274.223015.621258.124273.663015.5431118.174133.813015.641118.384133.73015.321118.264133.863015.61118.174133.73015.531118.24133.893015.691118.144133.83015.66平均1282.3734297.9933015.621282.2244297.6523015.428总中数1286.8074302.4743015.6671286.7094302.1613015.452一对标尺零点不等差=

15、0.二等水准点之记 线 点名： 标石断面图图幅编号标石类型经纬度标石质料所在地湖北武汉市武汉大学内地别土地使用单位武汉大学交通路线距广八路830m点位详细说明接管单位武汉大学保管人武汉大学选点单位埋石单位武汉大学维修单位武汉大学选点者埋石者维修者选点日期埋石日期维修日期备注基上高程（暗）高程系统基下高程（明）县（市）检查者_ 市（地）检查者_测段高差配赋表对集中实习的认识本次实习主要分为二等精密水准测量和大地测量解算编程。由于以前做过四等水准测量的实习，加上本次实习前老师发给我们的测量规范，我们对这次的二等精密水准测量上手很快。二等水准的精度要求比以前四等的要高很多，特别是测量方法步骤有一定的

16、不同。规范上说，对于二等水准，中午12点到2点是不宜进行测量的，因此为了保证一天的作业时间，我们早上七点钟就出发。这样我们大概二天就对珞珈山进行了往返测。本次水准测量每班三组人员，每组分别领到N3，007，电子测量仪器，这样我们一次性就学会了如何使用该三种水准仪器，按规范上说，对于光学仪器单站视距差不宜超过1米，累计视距差不宜超过3米;而对于电子仪器单站视距差不宜超过1.5米，累计视距差不宜超过6米。作业过程中我们都严格按照规范上的要求进行。实习过程中我们深深体会到了小组合作的力量，我们组一共六个成员，刚好二人扶尺，一人读数一人记录，二人量距，每当读数人员读好数并确定可以搬站时，量距的同学已经

17、提前量好了二到三站了，保持这样只提前量好二到三站其实有很大的意义：既可以使后尺和仪器无延迟的移动到相应的点上并调整，又可以避免出现量距差错需要重新量距却已经找不到量距人员影子的现象。一整天下来的外业测量使我们从这次实习中体会到测绘的辛苦与坚持不懈的精神！大地测量解算编程主要是在编程方面。但是在编程之前一定要把书本上关于解算的思路和原理弄清，否则在编程的时候在临时去看会降低速度和效率。在编程的时候会耗费大量的精力和时间，并且编绎下来也许会有很多的错误与警告，这时就需要耐心并静下心来回头检查程序，对于复杂的程序，可以把它分成几部分来编，进行计算的时候，通过正反算的数据可以检验编写的程序是否有误。由

18、于我没选修C#，因此没能把编好的程序封装模块化。二、在实习中的工作及在实习中的收获：这次实习让我们小组所有成员都得到了充分的锻炼，每个人都进行过各项任务的所有外业工作，并且都独立的完成了一套数据完整的内业检验和计算工作，每个组员的测绘技能都有较大的提高。通过这次实习测量我学会了在实际的野外测量工作中会在课本上所学到的知识不会完全一样，课本上所交给我们的都使在理想环境下的理论知识，真正的野外作业和书本上说的会有很大的不同。关于仪器的检验以及复杂地形的处理都会比课本上的有所提高。比如在我们测量过程中所遇到的地形起伏太大的问题，如果量距的成员在坡度较陡的路段能有效的控制视距，使得上，下丝读数没有超过

19、规范的要求，那么进度将会很快，虽然有时候量距人员没有控制得太好，但是我们没有人去责怪任何人，我们都是提醒一声然后量距的重新量距，量好后整平观测记录因此我们快速又成功的完成了水准测量的实习任务。实习中，每个组员的表现都十分出色，能力强的求快求准，做到精益求精，能力稍差的组员也能按照实习要求完成各项任务，没有人一味的为了实习进度而不顾实习的完成质量，只做个人比较擅长的事情，因为大家都知道这是我们整个小组所有组员的实习！三、工作量完成情况这次实习我们按照老师和任务书上的要求按质按量的完成了任务，主要成果包括：（1）、水准测量观测手薄；（2）、侧段高差计算；（3）、i角测定成果和水准尺零点差和基辅差常

20、数的测定成果；（4）、水准点点之记；（5）、水准测量每日安排日程及工作量完成情况登记表；（6）、实测斜距化算至高斯平面距离的程序并附有结果；（7）、高斯正反算程序并附有结果；（8）、大地主题正反算（包括白赛尔和高斯平均引数两种方法）程序及其结果；（9）、精密水准的计算表格。 四.编程思想按照书上的解算过程写出代码，这几个解算编程都比较复杂，因此基本上都是将一个解算拆成好几个部分各自编写成一段段的函数，通过主函数调用实现。而且编程过程中发现中间转换变量非常多，并且有些中间变量的名字代号都很相近，因此甚至将部分中间转换变量的转换独自写成一段函数，通过其它函数来调用实现，这样既方便又可以避免混淆。所

21、写的函数格式基本都是如下格式：精度类型 函数名（函数变量）一系列运算;返回函数值;控制输出精度的命令有：cout.precision();和coutsetiosflags(ios:fixed);coutsetprecision(4);五计算公式1.高斯投影 和 投影换带高斯正算:x=X+NsinBcosBl2/22+NsinBcos3B(5-t2+92+44)l4/244+NsinBcos5B(61-58t2+t4)l6/7206y=NcosBl/+Ncos3B(1-t2+2)l3/63+Ncos5B(5-18t2+t4+142-582t2)l5/1205高斯反算:B=Bf-tfy2/2MfN

22、f+tf3(5+3tf2+f2_9f2tf2)y4L=y/NfcosBf (1+2tf2+f2)y3/6Nf3cosBf + (5+28tf2+24tf4)y5/120Nf5cosBf2.大地主题解算高斯平均引数法 正算大地主题解算高斯平均引数法 反算大地主题解算白赛尔 正算大地主题解算白赛尔 反算第一次趋近时，取=0=L+实习编程一：实测斜距化算至高斯平面距离利用高斯正反算程序算出4，5号点的X , Y坐标：4号点坐标5号点坐标再计算4号点至5号点的大地方位角：A12=202度56分27.972秒最后得实测斜距化算至高斯平面距离：实习编程二：高斯正反算所用数据为第六组：B=35.L=112.

23、正算：反算：实习编程三：大地主题解算高斯平均引数所用数据为第六组：B1=42. ， L1=136. ， A12=1.4943 ， S12=90000正算：反算：白赛尔方法所用数据为第六组：B1=42. ，L1=136. ，A12=1.4943 ，S12=正算：反算：平面坐标转换4参数计算部分主要代码一实测距离化算至高斯平面距离void main() double B1,B2,A,D,H1,H2,y1,y2; double b13,b23,a13; coutb10b11b12b20b21b22a10a11a12DH1H2y1y2; B1=b10+b11/60+b12/3600; B2=b20+b

24、21/60+b22/3600; A=a10+a11/60+a12/3600; double W=sqrt(1-e2*sin(B1)*sin(B1); double N=a/W; double R= N/(1+e2*cos(B1)*cos(B1)*cos(A)*cos(A); double m=(H2-H1)*(H2-H1)/(D*D); double n=(1+H1/R)*(1+H2/R); double S=D*sqrt(1-m)/n)+D*D*D/(24*R*R); coutsetiosflags(ios:fixed); coutsetprecision(4); cout先将实测斜距换算至

25、椭球面上得距离S:nS=Sendl; double ym=(y1+y2)/2; double bm=(B1+B2)/2; double Rm=N/(1+e2*cos(bm)*cos(bm)*cos(A)*cos(A); double d=S*(1+ym*ym/(2*Rm*Rm)+(y2-y1)*(y2-y1)/(24*Rm*Rm)+ym*ym*ym*ym/(24*Rm*Rm*Rm*Rm); cout然后将椭球面上距离换算至高斯平面距离d:nd=dendl;二高斯投影正反算（其中正反算函数均按书上公式敲出来的，没意思，故没放上来）void Degree_change(double point)/

26、角度转换 int du,fen;double miao;du=int(point/3600.0);fen=int(point-du*3600.0)/60.0);miao=(point-du*3600.0-fen*60.0);coutdu度fen分miaoi) return (6*(i+1)-3)*3600/p;else return (6*i-3)*3600/p;/返回弧度void main()double x,y,B,L;int choose;cout请选择高斯投影类型：1、高斯投影正算. 2、高斯投影反算. choose;if(choose=1) double b3,l3;cout请按顺序

27、输入B、L的值(B、L度分秒之间用空格隔开)b0b1b2l0l1l2;B=(3600.0*b0+60.0*b1+b2)/p;/将输入的B变为弧度制L=(3600.0*l0+60.0*l1+l2)/p;/将输入的L变为弧度制LB_xy(B,L); if(choose=2)cout请按顺序输入x、y的值.xy; xy_BL(x,y);三 大地主题高斯平均引数int i=1;/控制循环迭代的代码doHSW(mB);HSV(mB);HSM(W);HSN(W);HSt(mB);HSn(mB);zBi=V*V/N*p*S*cos(mA/p)*(1+S*S/(24*N*N)*(sin(mA/p)*sin(m

28、A/p)*(2+3*t*t+2*n)+3*n*cos(mA/p)*cos(mA/p)*(-1+t*t-n-4*t*t*n);zLi=p/N*S/cos(mB/p)*sin(mA/p)*(1+S*S/(24*N*N)*(sin(mA/p)*sin(mA/p)*t*t-cos(mA/p)*cos(mA/p)*(1+n-9*t*t*n);zAi=p/N*S*sin(mA/p)*t*(1+S*S/(24*N*N)*(cos(mA/p)*cos(mA/p)*(2+7*n+9*t*t*n+5*n*n)+sin(mA/p)*sin(mA/p)*(2+t*t+2*n); mB=B+0.5*zBi;mA=A+0

29、.5*zAi; i+;while(abs(zBi-1-zBi-2)0.|abs(zLi-1-zLi-2)0.|abs(zAi-1-zAi-2)0.);void main()/主函数double change10;double anti_change12;int choose;coutchoose;if(choose=1)cout请输入B,L,A的度分秒和S（以空格键作为分隔）: n;for(int i=0;ichangei;double B=change0*3600+change1*60+change2;double L=change3*3600+change4*60+change5;doub

30、le A=change6*3600+change7*60+change8;double S=change9;Gauss(B,L,A,S);else if(choose=2)cout请输入B1,L1,B2,L2的度分秒（以空格键作为分隔）: n;for(int i=0;ianti_changei;double B1=anti_change0*3600+anti_change1*60+anti_change2;double L1=anti_change3*3600+anti_change4*60+anti_change5;double B2=anti_change6*3600+anti_chang

31、e7*60+anti_change8;double L2=anti_change9*3600+anti_change10*60+anti_change11;f_Gauss(B1,L1,B2,L2);elsecout您输入的有误!endl;大地主题白塞尔方法除了正反算函数是按书上公式一步一步敲出来的外，其余控制迭代循环代码，主函数代码和确定方位角大小的代码均与“大地主题高斯平均引数法”大同小异四平面坐标转换4参数计算void matrix(int bX,int cY)/计算二个矩阵相乘int i,j,k; for(i=0;iX;i+) for(j=0;jY;j+) for(k=0;kY;k+)

32、aij+=bik*ckj; void transpose(int bXY,int b_tYX)/计算转置矩阵int i,j;for(i=0;iX;i+)for(j=0;jY;j+)b_tji=bij;int inverse(int n,double a)/计算逆矩阵int *is,*js,i,j,k,l,u,v;double d,p;is=(int *)malloc(n*sizeof(int);js=(int *)malloc(n*sizeof(int);for(k=0;k=n-1;k+)d=0.0;for(i=k;i=n-1;i+)for(j=k;jd) d=p;isk=i;jsk=j;if

33、(d+1.0=1.0)free(is);free(js);return(0);if(isk!=k)for(j=0;j=n-1;j+)/行交换u=k*n+j;v=isk*n+j;p=au;au=av;av=p;if(jsk!=k) for(i=0;i=n-1;i+)/列交换u=i*n+k;v=i*n+jsk;p=au;au=av;av=p;l=k*n+k;al=1.0/al;for(j=0;j=n-1;j+)if(j!=k) u=k*n+j;au=au*al;for(i=0;i=n-1;i+)/消元计算if(i!=k)for(j=0;j=n-1;j+)if(j!=k)u=i*n+j;au=au-ai*n+k*ak*n+j;for(i=0;i=0;k-)/恢复if(jsk!=k)for(j=0;j=n-1;j+)/行交换u=k*n+j;v=jsk*n+j;p=au;au=av;av=p;if(isk!=k)for(i=0;i=n-1;i+)/列交换u=i*n+k;v=i*n+isk;p=au;au=av;av=p;free(is); free(js);return(l);/正常返回