《大地测量实习报告(武汉大学版本).docx》由会员分享,可在线阅读,更多相关《大地测量实习报告(武汉大学版本).docx(25页珍藏版)》请在第一文库网上搜索。
1、课程编号:课程性质:专业必修大地测量课程设计总结报告学院:测绘学院专业:测绘工程地点:武汉大学班级:2008级XX班组号:17姓名:XXX学号:XXXXX教师:郭际明丁仕俊苏新洲刘宗泉2011年06月27日至2011年07月15日目录封面错误!未定义书签。技术报告错误!未定义书签。1 .实习目的和意义错误!未定义书签。2 .实习任务错误!未定义书签。3 .测区概况错误!未定义书签。4 .已知高程错误!未定义书签。5 .作业依据(包括各种限差和要求)错误!未定义书签。6 .踏勘、选点错误!未定义书签。7 .使用的仪器和仪器检验48 .精密二等水准数据采集和外业概算过程中的有关情况49 .精密二等
2、水准测量外业高程概算成果410 .大地测量计算所采用的编程语言、各个任务的框图和编程思想以及基本数学模型等411 .大地测量计算成果汇总8实习报告1.1 技术报告1.1.1 实习的目的和意义本次大地测量计算与实习是在我们完成大地测量学课程的基础上进行的,通过这次实习,我们锻炼并提高了仪器操作能力、实际动手能力和团队协作能力,同时巩固了在课堂上学到理论知识。同时在后期的编程计算过程中,提升了我们发现问题、提出问题、解决问题的能力,将学到的知识进一步发散,提升自己的计算机能力,加强创新能力,为后期的工作打下基础。1.1.2 实习任务本次实习共有两项任务(1)二等精密水准测量外业观测与概算(2)大地
3、测量1.1.3 测区概况本次我们17小组抽取的测区是珞珈山测区,这段水准线路环绕武汉大学文理学部珞珈山一周,较大二上学期进行的三、四等水准测量线路较短,整个测段除了从教务部至珞珈山庄和教务部至校医院地势起伏较大外,其余测段都较为平坦,地质为混凝土和柏油马路,由于部分测段在居民生活区,来往车辆较多,道路为典型的盘山公路。卜.图标注的是我的测段。1.1.4 已知高程本次我组使用的己知高程点为教务部高程点,高程为126.157m1.1.5 作业依据国家测绘局,国家一、二等水准测量规范2006-05-24测绘出版社,2010仪器的检验水准仪的i角限差为15(2)标尺的零点不等差为0.10mm观测方式。
4、仪器类型光学水准仪数字水准仪备注前后视距差1.0m1.5m任一测站前后视距累积差3.0m6.0m测站观测顺序和方法奇数站后前前后奇数站后前前后左边表格均为往测,光学仪器返测时与往测相反,数字水准仪返测与往测相同偶数站前后后前偶数站前后后前闭合差往返高差不符值4kmm(k为测段长度,单位为km,40.1km按OzIkm算)环闭合差4Fmm(F为环线长度,单位为km)1.1.6 踏勘、选点本次实习的路线已经确定,踏勘选点的过程比较简单,最后我们组选取了教务部、政管院、研究总院指示牌、校区警务室、老同志干休所以及校医院六个点。1.1.7 使用的仪器和检验本次我组使用的仪器是TrimbIe的DiNi电
5、子水准仪以及配套水准尺,精度很高。我们进行了水准仪i角检验和水准尺零点差检验(结果另附),完全符合测量规范要求。1.1.8 精密二等水准数据采集与数据概算水准线路图已在测绘概况中绘出,观测日期与观测时段在观测记录薄中记载详细,数据记录规范、清晰。1.1.9 详见附录61.1.10 程序部分本次大地测量采用的是C#语言,C#是一种安全的、稳定的、简单的、优雅的语言,由C和C+衍生出来的面向对象的编程语言。它在继承C和C+强大功能的同时去掉了一些它们的复杂特性。C#综合了VB简单的可视化操作和C+的高运行效率,以其强大的操作能力、优雅的语法风格、创新的语言特性和便捷的面向组件编程的支持成为.NET
6、开发的首选语言,并且C#成为ECMA与ISO标准规范。C#看似基于C+写成,但又融入其它语言如PascaKJavaVB等。与C+不一样的是C#不支持对指针的访问,保证了程序的健壮性,但是不适应于对于指针的操作。此次程序共有四个窗体,一个类组成,类是白塞尔大地主题解算的算法。第一个窗体为总窗体,其余的分别是高斯正反算和实测斜距归算至高斯平面和白塞尔大地主题结算和一个GaussCacuIate的类。程序主体一、高斯投影正反算正算是指:由大地坐标(1,B)求得高斯平面坐标(x,y)的过程。反算是指:由高斯平面坐标(x,y)求得大地坐标(1,B)的过程。正算:高斯投影必须满足的三个条件:(1),中央子
7、午线投影后为直线。(2),中央子午线投影后长度不变。(3),投影具有正性性质,即正性投影条件。由第一个条件可知,中央子午线东西两侧的投影必然对称于中央子午线。设在托球面上有P1,P2,且对称于中央子午线。其大地坐标为(I,B),(I,B)则投影后的平面坐标一定为P1(x,y),P2(xz-y).由第二个条件可知,位于中央子午线上的点,投影后的纵坐标X应该等于投影前从赤道量至该点的子午弧长。计算公式:1当将克拉索夫斯基椭球带入计算式,可得到正算公式:(B”X=6367558.4969-a0-0.5+(a4+a612)1212NsinBcosBy=1+(a3+a512)I21NcosB其中:1(1
8、-10)PN=6399698.902-21562.267-(108.973-0.612cos2B)cos2Bcos2Ba0=32140.404-135.3302-(0.7092-0.0040cos2B)cos2Bcos2Ba4=(0.25+0.00252cos2B)cos2B-0.04166a6=(0.166cos2B-0.084)cos2Ba3=(0.3333333+0.001123cos2B)cos2B-0.16666667a5=0.0083-0.1667-(0.1968+0.0040cos2B)cos2Bcos2B2.反算公式为:p=Bf-1-(b4-0.12Z2)Z2Z2b2p1=1-
9、(b3-b5Z2)Z2Z2p,(1=10+1其中:Bf=+50221746+293622+2350+22cos2)cos2cos210-10sincosp=p,p6367558.4969PZ=y(NfcosBf)Nf=63995698.902-21562.267-(108.973-0.612cos2Bf)cos2Bfcos2Bfb2=0.5+0.003369cos2Bf)sinBfcosBfb3=0.333333-(0.166667-0.001123cos2Bf)cos2Bfb4=0.25+(0.16161+0.00562cos2Bf)cos2Bfb5=0.2-(0.16667-0.0088c
10、os2Bf)cos2Bf这道题目由于己经给出了相当精确地电算公式,因此过程比较简单,二实测斜距化算到高斯平面距离假设1、2两个大地点在椭球面上沿法线的投影点r和,间的大地线的长度为s,由于在椭球面上两点间大地线长度与相应法截线长度之差是极微小的,可以忽略不计,则可以将两点间的发截线长度认为是该两点间的大地线长度。并且,两点间的发截线的长度与半径等于其起始点曲率半径的圆弧长相差也很小,则所求的大地线长度可以认为是半径。其计算如下:S=D*1-(h2-h1)D*(h2-h1)D(1+h1Ra)*(1+h2Ra)这个题目的思想是先利用题目所给的XY坐标求出其大地纬度和大地经度,具体投影可以用三度带和
11、六度带,即先利用高斯投影正反算求出经纬度,再利用第三题白塞尔大地主题结算求出在一号的大地方位角,具体流程图如下。三大地主题正反算大地主题解算:知道某些大地元素推求另一些大地元素的过程。正解是指:已知某点P1的大地坐标(12,B2),且知该点到另一点P2(12,B2)的大地线长及其大地方位角A12,计算P2点的大地坐标(12,B2)和大地线在P2点的反方位角A21.的过程。反解是指:已知P1和P2的大地坐标(11,B1)和P2(12zB2)计算P1至P2的大地线长,正反方位角AI2、A21的过程。大地主题解算的基本思想:将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上,继而在球面上进行大地
12、主题解算,最后在将球面上的计算结果换算到椭球面上。其关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式,同时解决在球面上进行大地主题解算的方法。白塞尔的三个投影条件:(1)椭球面大地线投影到球面上为大圆弧;(2)大地线和大圆弧上相应点的方位角相等;(3)球面上任意一点的纬度等于球面上相应点的归化纬度。大地主题解算的步骤:(1)按椭球面上的已知值计算球面相应值,即实现椭球面向球面的过渡;(2)在球面上解算大地问题;(3)按球面上得到的数值计算椭球面上的相应数值,即实现从圆球向椭球的过渡。此次选用的是白塞尔大地主题解算。白塞尔大地主题正算公式:中间量:1-sinB11-e2cosB1W1
13、=1-e2sinzB1,Sinu1=-,cosu1=W1W1辅助函数:COSU1COsA1sinAn=cosu1sinA,cot1=:snu12cot1cot21-1sin21=,cos2=cot2i+1cot2+1球面长度:0=S-(B+Ccos21)sin21-A,sin2(1+)=0+B+5Ccos2(1+0)JA经差改正数:=+sin2(1+0)sin21sinA0终点大地坐标及大地方位角:SinU2=sinu1cos+cosu1cosA1sinsinu2=arctan|.,11e251sin2U2sinA1sin入=arctan:-Cosu1cossnu1sncosA112=1+入一
14、rCosu1SinA1A2=arctan:Cosu1coscosA1-snu1sn白塞尔大地主题反算公式:辅助计算:W1=y/1e2sin2B1,W2=y/1-e2sin2B2sinB11e2sinB21e2SinU1瓯,SinU2=-cosB1CosB2Cosu1=COSU2=W1W21=12-1a1=SinUISinU2,a2=cosu1COSU2b1=COSUISinU2,b2=SinU1COSU2采用迭代法同时计算起点大地方位角、球面长度及经差入=1+5第一次趋近,取6=0:p=CoSU2sin入,q=b1-b2c0s.A1=arctan-qsin=psinA1+qcosA1,cos=a1+a2c0s.sin=arctancossinA0=COSUiSinA1,x=2a1-cos2A0cos=xsinsinA0将计算得到的再带回计算经差,直到最后两次6相同或小于给定的允许值。大地线长y=(cos2A0-2x2)cosS=A+(Bx+Cy)sin反方位角A2=arctan(COSu1Sin入bcos.-b2三、大地测量计算成果(