大地测量计算与实习报告参考模板

大地测量计算与实习报告大地测量计算与实习报告大地测量计算与实习报告大地测量计算与实习报告1/201/20课程编号：20132024093 课程性质：必修 大地测量计算与实习实习报告学院： 测绘学院专业方向： 测绘工程（C组）实习地点： 武汉大学班级组号：六班第2组学生姓名： 学号： 指导教师：丁士俊、刘宗泉、苏新洲、史俊波2014年6月17日至2014年7月9日TOC\o"1-5"\h\z\o"CurrentDocument"一二等精密水准测量 31.1实习目的及意义 31.2实习任务 31.3已知高程 31.4作业依据 31.5测区概况及踏勘、选点 41.6使用的仪器和检验 51.7精密二等水准数据采集与数据概算 5\o"CurrentDocument"二大地测量计算 62.1实习目的 62.2.具体要求 62.3程序编程思想与基本框图 6\o"CurrentDocument"2.4大地测量计算成果汇总 14三实习总结 14四附录 16附录1数字水准仪的i角检验结果表 16附录2数字水准尺的零点差检验结果表 16附录3水准点点之记 16附录4自测二等水准测量外业观测高差与高程概算成果表 16附录5精密水准测量外业观测高差与高程概算成果表 16五实习成绩评定表 17一、二等精密水准测量1.1实习目的及意义了解水准路线的布设和点的位置，介绍水准点位置的选择方法。每人完成水准仪i角测定；一对水准标尺零点差及基辅差常数的测定；一个水准点点之记的绘制。除此之外，按二等水准测量的精度要求，每人完成一个闭合环观测的水准测量观测。每人轮流进行测量中的各种工作，并且对自己观测的成果进行百分之百的检查，同时完成全组的外业高差与概略高程表的计算一份。在进行测量的同时，每组应同时记录当日日程安排及工作量完成情况登记表。本次《大地测量计算与实习》是在我们完成《大地测量学基础》课程的基础上进行的，通过完成水准仪有关仪器检验和二等水准测量，我们锻炼并提高了仪器操作能力、实际动手能力和团队协作能力，同时巩固了在课堂上学到理论知识。在后期的编程计算过程中，提升了我们实际解决问题的能力，将学到的知识进一步发散，锻炼了编程能力，为后期的工作打下基础。1.2实习任务二等精密水准测量外业观测与概算的实习任务具体来说，有以下几点：踏勘测区：了解水准路线的布设和点的位置，介绍水准点位置的选择方法；填写水准点之记：每人完成一个水准点点之记的绘制；仪器检验：每人完成水准仪i角测定成果和一对水准尺零点差及基辅差常数的测定各一份；水准观测：按二等水准测量的精度要求，每人使用数字水准仪至少完成80站的观测和记簿，每个班级按照指定闭合环进行观测。每人轮流进行观测、记录、扶尺、量距、打伞等工作，并按各组记录登；每人使用光学精密水准仪至少完成10站的观测和记簿，每组构成闭合环；每人完成全组的外业高差与概略高程表的计算一份；在进行测量的同时，每组应同时记录当日日程安排及工作量完成情况登记表。1.3已知高程各个闭合环的高程已知点如下:点名已知高程（m）备注郑汉85133.040郑汉85（BM1）未名湖（BM8）未名湖119.888但是环二内没有已知咼程点，所以测的时候，设起始点B01高程为100.0m。

1.4作业依据国家测绘局，《国家一，二等水准测量规范》2006-05-24,测绘出版社，2010.各项限差：水准仪的i角限差为15〃；标尺的零点不等差为0.10mm；(3)仪器类型光学水准仪数字水准仪备注前后视距差Wl.OmV1.5m任——测祜前后视V60m|丄以1口IJiJ丿口17b距累积差W3.0m视线长度V50nm视线长度 I1—ii.50UI视线咼度N0.3mV2.80m光学为下丝读数且叩0.55m.测站观测顺序和奇数站：后刖刖后奇数站：j—— 、，• 、，•j——后刖刖后左边表格均为往测，光学仪器返测时与往测相反，数字水准仪返测与往测相冋方法闭合差偶数站：刖后后刖偶数站：刖后后刖!D段长度，单位为km,V0.1km按0,1km算)往返咼差不符彳I4振mm(k为测环闭合差4仰mm(F为环线长度，单位为km)1.5测区概况及踏勘、选点数字水准仪测区：大地测量计算与实习报告大地测量计算与实习报告用呻—W7成廿B03W3W心.」忏|*护用呻—W7成廿B03W3W心.」忏|*护"專”Z光学水准仪测区:食龍1诂剤中.側題咛~|脂忠学部敦四fjjJLl-1-J•讶园夢工13松也工I虫栋謂I食龍1诂剤中.側題咛~|脂忠学部敦四fjjJLl-1-J•讶园夢工13松也工I虫栋謂I：服性匸］株維姓我苴琶陽_芯顼华；敦-楼摘丁"I棺，.f=«H.呵：：十「；忤卜忠眼•八椅：Q :□后勤此円I•垂2力宝_1'永蛉塑暧I：①［坚臼製懒怎）「雪安塑像DJ酮广场是傅愤"杵 g匚 ］「飢但半醐善.「部活，加卩心.牧工“栋i'i'V•悻本次实习的路线已经确定为环三，踏勘选点的过程比较简单。此次二等水准精密测量在樱园环（环3）和信息学部进行，路上车辆较多，所以测量过程中要特别注意仪器安全。测量时经常经过马路，较危险，需要兼顾人身和仪器安全。樱园环距离长，上下坡较多，并且人流量与车流量大，容易对实习中的人员和仪器造成干扰。选点时，我们尽量选在道路两旁，通视较好，地面状况良好，距离适当的地方。1.6使用的仪器和检验本次我组使用的仪器是南方的电子水准仪DL-2007和007自动安平水准仪以及配套水准尺，精度较高。我们按照<<国家一.二等水准测量规范"附录的方法求对电子水准仪进行了两项检验（1）数字水准仪的i角检验；（2）数字水准尺的零点差检验，完全符合测量规范要求。（结果另附）下面将检测方法列出：（1）数字水准仪的i角检验方法在一平坦的场地上用钢卷尺依次量取一直线AI1I2B,其中I1，I2为安置仪器处，A，B为立标尺处，在线段AI1I2B上使得AI1=I2B，在I1，I2处先后安置仪器，仔细整平仪器后分别在A,B标尺上各照准读数基本分划四次。将仪器温度与环境温度充分平衡，并开机预热后方可进行检测。（2） 数字水准尺的零点差检验方法在距水准仪20至30米的地方依次打三个尺桩，使桩顶间高差约20cm。此项检验共进行3个测回，每一个测回分别在3个桩上安置一对标尺，每测回每桩连续观测4次。。1.7精密二等水准数据采集与数据概算小组观测日期：2014年6月28号——2012年7月9号仪器检验时间：数字水准仪的i角检验（2014年6月29号早上8点到12点）数字水准尺的零点差检验（2014年6月29号下午14点到17点）往返测时段：（1） 数字二等精密水准：往测：武大樱园环（2014年6月30号）返测：武大樱园环（2014年7月2号）往测：武大樱园环（2014年7月3号）返测：武大樱园环（2014年7月6号）（2） 光学二等水准（2014年7月1号）二、大地测量计算2.1实习目的与要求通过大地测量计算的课程设计，特别是根据数学模型公式，自己动手编程，计算出正确结果，从而对大地测量中主要计算公式有更深入的了解和掌握。2.2实习任务 （1） 每人完成实高斯投影正反算程序编制，并计算结果；（2） 每人完成大地主题正反算公式的计算程序编制，并计算出正确的结果；（3） 每人完成不同平面之间的坐标系之间转换程序编制；（4） 根据提供数据进行二等精密水准外业概算的计算工作，并与其他组的数据进行平差计算2.3大地测量计算编程思想与基本框图高斯投影正反算和大地主题正反算编程，我采用的是C#语言，C#是一种安全的、稳定的语言，设计窗体程序特别方便简单。它在继承C和C++强大功能的同时去掉了一些它们的复杂特性。与C++不一样的是C#不支持对指针的访问，保证了程序的健壮性，但是不适应于对于指针的操作。该程序共有2个窗体组成，第一个窗体为总窗体，其余的分别是高斯正反算和大地主题解算。不同平面坐标系之间的转换编程，我米用的是Matlab编程。Matlab对矩阵计算具有得天独厚的优势，仅需几十行代码即可完成矩阵运算。窗体截图（1）高斯投影正反算正算是指：由大地坐标（L,B）求得高斯平面坐标（x,y）的过程。反算是指：由高斯平面坐标（x,y）求得大地坐标（L,B）的过程。将克拉索夫斯基椭球带入计算式，可得到正算公式：x=6367558.4969^-（aD-[0.5+（a4+I3）]3]lzN}sinBcasB、 y=[l+R+W邛屮阪汨其中：I（L-LaFN=6399^98.902-[21562.267-(108.973-0.612caS2E)coK2B]coK2EaD=5214&.404-[135.3302-(0.7092-D.CI04Ocas2B)coK2B]cOK2B%=(0,25+0,OO252cos2B}coE=B-0,04166afi=(0,166cG5=B-0,084}ca52Bag=(fl.3333333+0.001123cas2B)ais2B-0,16666667a5=0.0083-[0.1667-(0.1968+O.O&40cas2B)cOs2B]cos£E大地测量计算与实习报告大地测量计算与实习报告反算相应计算公式：p=Bf-[1-翌旳巧理扩l=[l-(b5-b5Z2)Z2]ZVf[ L=La+I其中：Bf=P+(50221746+[293622+{2350+22cos2p)cos2p]cos2p)l&-lcsinpcospp,f§=63675504969PZ=y/(NfcosBf)Mf=的99S的&902-[21562.267-(108.973-0.612cos2Bf)cas2Bf]cas2BEb2=(0.5+O.OOSS^^cas^BfJsiuBfCosBfba=0,333333-〔0,166667-0,001123cas=Bf)coE=BEb4=0.25+(0.16161+。.胸562a*盼亦卫再bE=0.2-(0,16667-0.0088cos2Bf)cos=Bf程序结构图：开始计算中央子午线经度Lo；所在带号；经差l选择椭球参数：e,a等；分带方法：3开始计算中央子午线经度Lo；所在带号；经差l选择椭球参数：e,a等；分带方法：3°或6°带输入大地元素：B,L由r的前两位(带号)计算中央子午线经度L0计算x,y关于经差计算x,y关于经差l的级数展开系数：m0；m1 P168(4—363)(4—365)计算B,l关于x,y的级数展开系数：n,n,n1 2 3等。P170(4—381)(2)大地主题解算大地主题解算：知道某些大地元素推求另一些大地元素的过程。正算：已知某点P1的大地坐标(L2,B2),且知该点到另一点P2(L2,B2)的大地线长及其大地方位角A12，计算P2点的大地坐标(L2,B2)和大地线在P2点的反方位角A21.的过程。反算：已知P1和P2的大地坐标(L1,B1)和卩2(L2,B2)计算P1至P2的大地线长，正反方位角A12、A21的过程。高斯平均引数算法基本思想：首先把勒让德級数在P1点展开改在大地线长度中点M展开，以使級数公式项数减少，收敛块，精度高，其次将M点用大地线两端平均纬度及平均方位角相对应的m点代替，并迭代计算，直到限差小于规定值。 白塞尔方法基本思想：将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后在将球面上的计算结果换算到椭球面上。其关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式，同时解决在球面上进行大地主题解算的方法。白塞尔大地主题解算的步骤：按椭球面上的已知值计算球面相应值，即实现椭球面向球面的过渡；在球面上解算大地问题；按球面上得到的数值计算椭球面上的相应数值，即实现从圆球向椭球的过渡。大地主题解算的数学模型请参照《大地测量学基础》教材(1)高斯平均引数正算计算公式(s<200km)(1)计算辅助量公式e2e'2= 1一e2N= 1一e2sin2Ba(1一e2)M=\i(1-e2sin2B)3t=tgB n2=e'2C0S2B(2)计算AB">△〃、AA”的初值p"AB=S-cosA

0M1

1(3)计算B、mAL'= S-sinA0NcosB111p'AA"=——S-sinAtanB=AL'sinB0N 1 1 0 11L、Amm1 1B=B+—AB”m1 2 0L=L+—AL'>m 1 2 0A=A+1AAm 1 2 0(4)再次计算AB”、AL"、AA"S-cosA{1+-^^[sin2A(2+3t2+2q2)~ ’ m mmAB“工Mm 24N2mm+3cos2An2(t2—1-n2—4n212)]Jmmm mmmAL"=NcosBm mp" IS2 SsinA〈1+ [t2sin2A—cos2A(1+n2—9n212)]mI24N2mm mmmmmAA"=-P-SsinAN mmm+sin2A(2+12+2q2)]}mmmS2t{1+[cos2A(2+7n2+9n212+5n4)24N2 m mmmmm(5)重复计算(3)，直到计算满足AB"一…B„£i+1 i…L一…Lv£，i+1 iAA”一AA”v£i+1 i 如按孤度计算可取£=1x10一io,按角度计算可取£=0.0001'o(6)计算B、L、A的最后值2 2 2B=B+AB2 1iL=L+AL21 iA=A+AA+1800(A>1800取+,Av1800取一21 i 1 12)高斯平均引数反算公式(Sv200km)(1)B=i(B+B) AB=B-Bm2 1 2 2 1AL=L-L21(2)U=SsinA=rAL+rAB2AL+rAL3m01 21 03V=ScosA=sAB+sABAL2+sAB330m1012AA=t AL+1AB2AL+1 AL01 21 03式中各系数:r=NcosB01m mr03Nmcos3B12mm24s01N mV2mr21s12cosBm 24V4m(1+n2一9门212门4)m mmmN m—24V2mcos2B(2+3t2+2q2)m m ms3001=tcoss3001=tcosBm mt2124V4m1cosB(2+7门2+9t2门2)mmmN/ 、—亠(n2一12n2)8V6m1t=—cos3Bt(2+12+2^2)03 24mmmm注：这里对教材公式中相应系数进行了修正与改进。大地测量计算与实习报告大地测量计算与实习报告SsinA U13)A=arctg m=arctg-ScosA Vmarctan(4)T兀——farctan4SsinA mScosAm1-C1++C|涸>|AL||ab|<|al|ScosA mSsinAm兀/2当AB>0,

当AB<0,

当AB<0,

当AB>0,

当AB=0,AL>0AL>0AL<0AL<0AL>0⑸S=SsinAsinAmU

sinAmScosAS= mcosAmVcosAm1A=A+-AA土兀2m2程序结构图：采用高斯平均引数算法开始反算K=1；2正算r选择椭球参数：/e，a；大地元素:/开始反算K=1；2正算r选择椭球参数：/e，a；大地元素:/B1，L1，B2，L2计算：U；V；AAP137(4—226)(4—228)计算：AB0；AL'。；AA0P136 (4—221)(3)坐标系之间转换的转换参数和转换精度计算模型如下:Ax+(1+m)cos0—sin0sin0cos0新坐标旧坐标构建向量：X=[abcd]

a=Ax、b=△yc=(1+m)cos9

d=(1+m)sin9由已知数据得到：V=[v1………VJB=[B11B48]L=[L「 得到V1=Bx—L大地测量计算与实习报告大地测量计算与实习报告2.4大地测量计算成果（1）高斯投影正反算计算结果：参考基准：克拉索夫斯基椭球②采用3°带，第42带投影③采用的是第二组第4号数据。巳知数据x=4585822.5452y=425000L0=126°正算x=4585822.5464y=42500.3050反算B=41.24115693L=125.06112054

F-^rrnl专氏崩注 ¥投影分苛鎚度L41£4专氏崩注 ¥投影分苛鎚度L41£415.9393摘球辭度B高斯投窮正兵苴大地主题解算坐标案转換正算¥带V12561L.2054正牌结果：s -1505022.5464y ^SD00.305说明：1、 输入角度值时，语以•童分秒'格式输入铀入聖标时,请輸入闻察大地坐标(X,Y),印F值前加帯号且如£00困H.梧&不.巳两位蛇时前面g-投黒鶴投黒鶴松肯 *4ES5022.5452D0l：4ES5022.5452D0l：.3OB0结果：反算纬度为41?史4'U.S国3”

金度为.12尸小".£054”反算(2)大地主题解算计算结果①采用高斯平均引数算②参考椭球为克拉索夫斯基椭球③采用5—6班第2组第4号数据。正算反算巳知数据计算结果巳知数据计算结果B1=43.18356874B2=44.35140096B1=43.18356874S=142000.0004L1=132.10122676L2=132.13376639L1=132.10122676A12=L49430007A12=1.4943A21=181.52055565B2=44.35140096A21=181.52055543S=142000L2=132.13376639(3)不同平面坐标系之间的转换，计算四个转换参数，评定转换精度采用Matlab进行编程计算，结果见下图。NewQpe-nSav*■*■■**■F=iPrintFIl£[nden?.EDITNe-wNtw Open阿Scrpt ■ ■faio.mUfx.=[21S3L.S51O;117363.2200:126133.03□□:120264.期純;101659.490Q]332060B207.173027IBS.0930-30775.459(1.43210.3240NewQpe-nSav*■*■■**■F=iPrintFIl£[nden?.EDITNe-wNtw Open阿Scrpt ■ ■faio.mUfx.=[21S3L.S51O;117363.2200:126133.03□□:120264.期純;101659.490Q]332060B207.173027IBS.0930-30775.459(1.43210.3240ri£ht=[S4S5eeO.250Z077.GS3OSloSao.04ZQfile<!■*►S$3 ►DrfQdlfi：件*1257481748►FileRerv*CammandWindow25055^3.5170□10747.0770:25375S4.4TL0491542.3L7QDld=ri£ht:□ld_K-ald(,1)Dl(Ly=oldS:,2):new™left;y2=：491B7.7330:132799r2230]:yl-[2524.192.3650；522075.4290：；B=[-10□-1-100-1-10□-1-old_z(1)-gld_y(l)-old_z⑵-old_y(2)-old_z⑴-old_y(3)fLd_y⑴■old_xl：L):fLU⑵:old_xl：2):-oLrLj'O):01d_x(3);13SEJSJ□3W3EI旬三;国»faio四堂粉为Zkx=-2^365701.口50=199,^=-^33212.253161fl(IS)—0,97064099GC19)—0,01