大地测量实习报告.doc

课程编号： 课程性质：必修大地测量计算与实习实 习 报 告学 院： 测绘学院 专业方向： 测绘工程 A方向 实习地点： 武汉大学 班级组号： 10班第小组 学生姓名： 学 号： 指导教师：苏新洲 刘宗泉 丁士俊向 东 黄海兰 等2012年 6 月18 日 至 2012年 7 月4 日目录封面1目录.2第一章：技术报告31.实习目的和意义32.实习任务33.测区概况44.已知高程45.作业依据（包括各种限差和要求）46.踏勘、选点67.使用的仪器和仪器检验68. 精密二等水准数据采集和外业概算过程中的有关情况79.精密二等水准测量外业高程概算成果910. 大地测量计算课程设计911.大地测量计算成果汇总16第二章：实习总结.21第三章：附录26封底.62第一章技术报告一：实习目的和意义通过完成水准仪有关仪器检验和二等水准测量，达到以下几点目的：1、巩固加深所学理论知识，培养理论联系实际、实际动手能力；2、初步掌握精密水准仪仪器主要检验方法和精密水准测量的观测程序，记录计算和高差计算等工作；3、熟练掌握常用水准测量仪器；4、熟练掌握高程概算的相关计算及其改正。本次大地测量计算与实习课程是在我们完成大地测量学课程的基础上进行的，这次实习，主要是为了提高我们同学对仪器操作能力、实际动手能力以及小组成员之间的团队协作能力，巩固在课堂上学到理论知识。在后期的编程计算过程中，提升了我们发现问题、提出问题、解决问题的能力，将学到的知识进一步扩展与发散，加强自己的计算机编程能力，提升创新能力，为以后同学们以后参加工作以及生产打下坚实的基础。本次实习要求学生提高认识，加强纪律，保证安全，培养团队精神，听从指挥，严守规范，认真操作，增强实践能力。二 实习任务本次实习共有两项任务一 二等精密水准测量外业观测与概算（1）踏勘测区：了解水准路线的布设和点的位置和水准点位置的选择方法。（2）填写水准点之记：每人完成一个水准点点之记的绘制。（3）仪器检验：每人完成水准仪i角测定成果和一对水准标尺零点差及基辅差常数的测定各一份。（4）水准观测：每人使用数字水准仪至少完成80站的观测和记簿，每小组完成珞珈山环，樱园环和武测环3个闭合水准环中的2个环线；（5）每人使用光学水准仪至少完成8站的观测和记簿，每组构成闭合环。（6）外业观测数据检查与概算。二 大地测量计算（1）高斯投影正反算的计算程序编制及计算结果；（2）实测斜距化算至高斯投影平面边长的边长改正的程序编制及计算结果；（3）大地主题正反算的计算程序编制及计算结果；（4）实测数据进行外业高差与概略高程表的编制。三 测区概况本次实习我们小组的测区为武测环和珞珈山环。武测环，起点是一号教学楼前郑汉85点，要求围绕武测校区观测一圈构成一个闭合环线。全部测段都位于水泥路面或柏油马路上，地质坚硬，测区中部分都有树荫。测段坡度很缓，高差较小，观测进度较快。珞珈山环路以教务部前水准点为已知点，以围绕珞珈山一周为水准环路，全路线位于硬化路面，测段坡度起伏较大，高差大，路线相对较长，但测段中水准点分布较为密集，对与分段测量较为便利,且较为利于检核测量中是否出现超限情况。我们还在国际软件学院的校园内完成了小组一共70站的光学水准仪测量，小组内部构成闭合环，路线平整，易于测量。四 已知高程武测环：郑汉85，高程为133.040m樱园环：未名湖，119.888m珞珈山环：教务部，126.157m五 作业依据水准测量方法：使用水准仪和水准尺，根据水平视线确定两点之间的高差，从而由已知点的高程推求未知点的高程。 各类测站的选择：在所选择的闭合回路上选取4个普通水准标石作为固定点，将整个闭合回路分为4个测段，每个测段中设置偶数个测站，站与站之间距离应适中。 施测方法及使用的仪器：二等水准测量中使用了两种仪器，一是007自动安平水准仪，二是科力达DL-201电子水准仪。按照国家二等水准测量技术要求并参照城市测量规范CJJ8-9进行施测。 A：限差：1水准仪的i角限差为152标尺的零点不等差为0.10mm。仪器类型光学水准仪数字水准仪备注前后视距差1.0m1.5m任一测站前后视距累积差3.0m6.0m测站观测顺序和方法奇数站后前前后奇数站后前前后左边表格均为往测，光学仪器返测时与往测相反，数字水准仪返测与往测相同偶数站前后后前偶数站前后后前闭合差往返高差不符值 4mm ( k为测段长度，单位为km,0.1km按0,1km算)环闭合差 4mm ( F为环线长度，单位为km )B：要求：1.仪器保护问题：在实习过程中中，测量人员应做到，人不离仪器太远，一般人应在仪器旁一米左右的范围；2. 保持后尺尺垫不动，扶后尺的同学应注意保证其尺垫不动。由于尺垫起着传递高程的作用因此在测量的过程中如果后尺的尺垫移动将会影响整个测量的精度，应从固定点重新测起；3. 记录人应时刻注意前后视距差，以便在测量的过程中及时的进行调整；4. 注意各个值得限差，如果超限就及时进行重测；5.注意奇数站和偶数站时观测尺的顺序，返测时进行刚好相反或交换前后尺；6.抚尺同学如果在观测过程中移动了标尺应及时报告读数的人员，进行重读；7.在老师给我们分配的路线中进行选点，尽量应选固定水准点。在进行测量前利用测绳量取大约相等的两段距离，然后利用油漆或其他标记物记下点的位置，以便以后测量时进行利用。在画仪器点和架尺点是建议用不同的符号，最好写上点号，以便掌握奇偶数站；8.在提前画好的点上架尺，固定点不要尺垫，而在自己选取的点上放尺垫，然后架尺。架尺时首先将两根竹竿交叉顶在尺上，然后调整竹竿让标尺上的圆水准器的水泡居中；9.在提前画好的点上架仪器，首先大体上整平仪器底座，然后调整脚架，使圆水准器上的水泡大体居中，然后利用微调螺旋使水泡居中；10.架好仪器后让望远镜对准尺子，调好视线和焦距进行测量。在往测中按照奇数站后前前后，偶数站前后后前的顺序进行测量；而返测时刚好相反或交换前后尺按照往测得顺序进行测量；11.同一测站观测不能两次调焦。12.测站高差超限，在本测站发现立即重测，迁站后应从水准点后者间歇点开始重测。测段往返高差超限，应进重测。当测段需重测时，应认真分析超限原因，选取可靠性较小的进行重测。六 踏勘、选点因为实习的路线已经确定，我们经过比较和筛选，发现已知刚点距离较远且为了防止一旦发生重测事故时距离过长，我们决定尽可能的多用已知刚点，多侧段的策略。七 使用的仪器和仪器检验7.1 使用的仪器水准仪（DL-201和007）及脚架一套、水准标尺2根、尺垫2个、测绳一根、背包一个、记录板一个、竹竿四根，水准测量观测手簿、i角检测表格、水准点之记、一对水准尺零点差及基辅差常数的测定表格，铅笔、小刀、油漆等。7.2 仪器检验：1. i角检验在较平坦的地方选定适当距离（取20.6m）的两点A、B，并用尺垫确定不动。置水准仪于A点附近，如图：此时测量正确的高差hAB，然后将水准仪置于B点的附件，在线段AB的延长线上， 在测得的高差hAB中将有i角的影响，按上面的公式有：i=( hAB- hAB) /(2SAB)2. 水准尺零点差检验：在距离水准仪约20m30m的等距离处选择3个点，使得各个点之间存在着约20cm的高差。此项检验进行3个测回，每一测回中，分别在3个点上一次安置一对标尺，分别用仪器进行读数，测回间应变换仪器高。由于使用的是数字水准仪，因而设置重复测量次数为5次，每测回每点连续观测4次。水准尺零点差即两次读数中数之差。八 精密二等水准数据采集和外业概算过程中的有关情况8.1 水准路线图武测环珞珈山环图中用箭头标明我所完成的测段。观测日期与观测时段在观测记录薄中记载详细，数据记录规范、清晰。8.2人员安排本小组共有成员7人，在施测水准测量的时候，人员分配情况是：两人为跑尺员（前、后尺各一人），一人观测，一人记录，两人负责量距，另外一人负责后勤和适当轮换跑尺。8.3具体作业方法在选取了合适的水准路线和固定点之后，开始进行第一测站的观测，将水准尺立于固定点上作为后视，水准仪放置杂水准路线附近合适位置，然后在施测路径前进方向上取仪器与后尺大致相等距离放置尺垫，在尺垫上树立前尺。随后观测员对水准仪进行整平。观测过程为往测时，奇数站按“后前前后”的顺序对后尺前尺进行读数；偶数站则按照“前后后前”读数；返测时，则与按照奇数站“前后后前”，偶数站则按照“后前前后”。在该测站完毕后，通知后尺移站，此时前一站的前视点变为后一站的后视点，按照奇偶工作程序完成该站的测量，直到完成该测段为止。九 精密二等水准测量外业高程概算结果详见附录十 大地测量计算课程设计10.1课程设计目的及主要内容在学生自主编程计算的过程中，提高发现问题，分析问题，解决问题的对策与方法的能力；开拓视野，把基本理论与实践紧密结合在一起，对理论知识进一步消化与吸收，激发和培养学生创新思维的能力。（1）高斯投影正反算的计算程序编制及计算结果；（2）实测斜距化算至高斯投影平面边长的边长改正的程序编制及计算结果；（3）大地主题正反算的计算程序编制及计算结果。10.2 编程语言简介输出x，y根据正算公式计算高斯坐标x，y计算正算公式所需的各个系数A2，A3，A4，A5，A6计算卯酉圈半径N输入B,L,L0 我采用的是C+语言，C+语言功能齐全。具有各种各样的数据类型，并引入了指针概念，可使程序效率更高。而且计算功能、逻辑判断功能也比较强大。 10.3 任务框图高斯投影正算输入x，y，L0高斯投影反算计算底点纬度Bf和对应的卯酉圈半径Nf计算反算公式所需的各个系数B2，B3，B4，B5，B6根据反算公式计算大地经纬度B，L输出B，L输出高斯投影平面边长D12利用参考椭球面的距离D和正常高及高程异常计算出高斯投影平面边长D12利用大地方位角和实测斜距将实测斜距化算至参考椭球面的距离D利用两点经差l和两点纬度计算出大地方位角A12利用高斯坐标反算出两点经差l输入两点高斯坐标x1,y1,x2,y2，两点大地纬度B1,B2，两点正常高h1,h2，实测斜距S12，高程异常值实测斜距改化大地主题正算计算框图输入B1,L1,A12,S输出A21,B2,L2将球面元素换算到椭球面上，计算出B2，L2将椭球面元素投影到球面上，计算球面长度解算球面三角形，计算出A21输入B1,B2,L1,L2大地主题反算计算框图输出A12，A21，S将球面元素换算到椭球面上，计算出S和A21采用迭代方法计算球面对应经差和球面长度及A12输入B1,B2,L1,L210.4 编程思想及基本数学模型1.高斯投影正反算正算是指：由大地坐标（L,B）求得高斯平面坐标（x,y）的过程。反算是指：由高斯平面坐标（x,y）求得大地坐标（L,B）的过程。正算：高斯投影必须满足的三个条件：（1）中央子午线投影后为直线；（2）中央子午线投影后长度不变；（3）投影具有正性性质，即正性投影条件。由第一个条件可知，中央子午线东西两侧的投影必然对称于中央子午线。设在托球面上有P1 ，P2，且对称于中央子午线。其大地坐标为（l,B）,(-l,B)则投影后的平面坐标一定为P1（x,y）,P2（x,-y）由第二个条件可知，位于中央子午线上的点，投影后的纵坐标x应该等于投影前从赤道量至该点的子午弧长。相应计算公式：X=6367558.496*B1/-a0-0.5+(a4+a6*l*l)*l*l*N*sinB*cosBY=1+(a3+a5*l*l)*l*l*l*N*cosB式中l =(L-L0)”/”N=6399698.902-21562.267-(108.973-0.612*cosB*cosB)*cosB*cosB*cosB*cosB a0=32.140.404-135.3302-(0.7092-0.0040*cosB*cosB)*cosB*cosB*cosB*cosB a4=(0.25+0.00252*cosB*cosB)*cosB*cosB-0.04166 a6=(0.166*cosB*cosB-0.084)*cosB*cosBa3=(0.3333333+0.001123*cosB*cosB)*cosB*cosB-0.1666667 a5=0.0083-0.1667-(0.1968+0.0040*cosB*cosB)*cosB*cosB-0.1666667 反算：在高斯投影坐标反算时，原面是高斯平面，投影面是椭球面，则有如下三个条件：（1）x坐标轴投影成中央子午线，是投影的对称轴；（2）x轴上的长度投影保持不变; （3）正性投影条件。相应计算公式：B=Bf-1-(b4-0.12*Z*Z)*Z*Z*Z*Z*b2*”l =1-(b3-b5*Z*Z)*Z*Z*Z*”L=L0+l 式中Bf=+50221746+293622+(2350+22*cos*cos)*cos*cos*cos*cos*(10+e-10)sin*co s*”=x/6367558.4969*”Z=y/(Nf*cosBf) Nf=6399698.902-20562.267-(108.973-0.612*cosBf*cosBf)*cosBf*cosBf *cosBf*cosBfb2=(0.5+0.003369*cosBf*cosBf)*sinBf*cosBfb3=0.333333-(0.166667-0.001123*cosBf*cosBf) *cosBf*cosBfb4=0.25+(0.16161+0.00562*cosBf*cosBf) *cosBf*cosBfb5=0.2-(0.1667-0.0088*cosBf*cosBf) *cosBf*cosBf2.实测斜距化算到高斯平面距离假设两个大地点在椭球面上沿法线的投影点间的大地线的长度为S 由于在椭球面上两点间大地线长度与相应法截线长度之差是极微小的，可以忽略不计，则可以将两点间的发截线长度认为是该两点间的大地线长度。并且，两点间的发截线的长度与半径等于其起始点曲率半径的圆弧长相差也很小，则所求的大地线长度可以认为是半径。本题的思路是，先将斜距归算成为大地线长，再将大地线长归算为高斯平面投影，但是若要计算出大地线的前提是要知道两点间的大地方位角，大地方位角与坐标方位角，子午线收敛角以及方向改化角之间有一定的数量关系，后三个角均可以由已知条件求出。所用到的计算公式如下：A12=apha12+r1-delta12;S=S0*sqrt(1-(H2-H1)/S0*(H2-H1)/S0)/(1+H1/RA)*(1+H2/RA) D=(1+ym*ym/2/Rm/Rm)*S3.大地主题正反算大地主题解算：知道某些大地元素推求另一些大地元素的过程。正解是指：已知某点P1的大地坐标（L2,B2），且知该点到另一点P2(L2,B2)的大地线长及其大地方位角A12，计算P2点的大地坐标（L2,B2）和大地线在P2点的反方位角A21的过程。反解是指：已知P1和P2的大地坐标（L1,B1）和P2（L2,B2）计算P1至P2的大地线长，正反方位角A12、A21的过程。大地主题解算的基本思想：将椭球面上的大地元素按照高斯平均引数条件投影到辅助球面上，继而在球面上进行大地主题解算，最后在将球面上的计算结果换算到椭球面上。其关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式，同时解决在球面上进行大地主题解算的方法。正算思路，通过已知条件首先求出各辅助量，然后计算deltaB,deltaL,deltaA的初值，再由初值算出Bm,Lm,Am,再通过以上三值重新计算deltaB,deltaL,deltaA,将这三个数值与之间序列的三个数值求差，再一定的限差要求范围内进行迭代，直到满足限差条件，最后求出B2,L2,A2所用到的计算公式：deltaB = S * Cos(A) / M;deltaL0 = S *Sin(A) / (N * Cos(B);deltaA0 = S * Sin(A) * Tan(B) / N;Bm = B + deltaB0 / 2;Lm = L + deltaL0 / 2;Am = A + deltaA0 / 2;deltaB = S * Cos(Am) * (1 + S * S * (Sin(Am) * Sin(Am) * (2 + 3 * tm * tm + 2 * gm2) + 3 * Cos(Am) * Cos(Am) * gm2 * (tm * tm - 1 - gm2 - 4 * gm2 * tm * tm) / 24 / Nm / Nm) / Mm;deltaL = S * Sin(Am) * (1 + S * S * (tm * tm * Sin(Am) * Sin(Am) - Cos(Am) * Cos(Am) * (1 + gm2 - 9 * gm2 * tm * tm) / 24 / Nm / Nm) / Nm / Cos(Bm);deltaA = S * Sin(Am) * tm * (1 + S * S*(Cos(Am) * Cos(Am) * (2 + 7 * gm2 + 9 * gm2 * tm * tm + 5 * gm2 * gm2) + Sin(Am) * Sin(Am) * (2 + tm * tm + 2 * gm2) / 24 / Nm / Nm) / Nm;B2 = B + deltaB;L2 = L + deltaL;A2=A+deltaA0+PI(或者A+deltaA0-PI）反算思路：与正算类似，通过已知条件先算出中间要用的一些辅助量，此处，用U和V来表示S*SIN(Am)和S*COS(Am)的数值，最终可解算出Am。具体步骤是：先求出r01,r21,r03,s10,s12,s30,t01,t21,t03这些系数，再代入到UV的表达式中，求出UV，再通过判断条件判断T与Am取值，最后即可求出S与A1，A2的值。所用到的公式：r01 = Nm * Cos(Bm);r21 = Nm * Cos(Bm) * (1 + gm2 - 9 * gm2 * tm * tm + gm2 * gm2) / 24 / Pow(Vm, 4);r03 = -Nm * Pow(Cos(Bm), 3) * tm * tm;s01 = Nm / Vm / Vm;s12 = -Nm * Pow(Cos(Bm), 2) * (2 + 3 * tm * tm + 2 * gm2) / 24 / Vm / Vm;s30 = Nm * (gm2 - tm * tm * gm2) / 8 / Pow(Vm, 6);t01 = tm * Cos(Bm);t21 = Cos(Bm) * tm * (2 + 7 * gm2 + 9 * tm * tm * gm2) / 24 / Pow(Vm, 4);t03 = Pow(Cos(Bm), 3) * tm * (2 + tm * tm + 2 * gm2) / 24;V = s01 * deltaB + s12 * deltaB * Pow(deltaL, 2) + s30 * Pow(deltaB, 3);U = r01 * deltaL + r21 * deltaB * deltaB * deltaL + r03 * Pow(deltaL, 3);deltaA = t01 * deltaB + t21 * Pow(deltaB, 2) * deltaL + t03 * Pow(deltaL, 3);C=Abs(V/U);S = U / Sin(Am);A1 = Am - deltaA / 2;A2 = Am + deltaA / 2 十一 大地测量计算成果汇总我选择十班第一组4号数据高斯投影正反算：原始数据 3度带投影，第42号投影带，X=4575017.3855 Y=437854.163计算结果 正算反算实测斜距化算至高斯平面边长大地主题解算（1） 白塞尔大地主题正算（2）高斯引数法正算大地主题反算（1）白塞尔大地主题反算（2） 高斯引数发反算附录五 大地测量计算源程序高斯投影正算/ ZhengsuanDlg.cpp : implementation file/#include stdafx.h#include gaosi.h#include ZhengsuanDlg.h#include math.h#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CZhengsuanDlg dialogCZhengsuanDlg:CZhengsuanDlg(CWnd* pParent /*=NULL*/): CDialog(CZhengsuanDlg:IDD, pParent)/AFX_DATA_INIT(CZhengsuanDlg)m_tuoqiu = -1;/AFX_DATA_INITvoid CZhengsuanDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CZhengsuanDlg)DDX_Radio(pDX, IDC_RADIO1, m_tuoqiu);/AFX_DATA_MAPBEGIN_MESSAGE_MAP(CZhengsuanDlg, CDialog)/AFX_MSG_MAP(CZhengsuanDlg)ON_BN_CLICKED(IDC_BUTTON1, OnButton1)/AFX_MSG_MAPEND_MESSAGE_MAP()/ CZhengsuanDlg message handlersvoid CZhengsuanDlg:OnButton1() / TODO: Add your control notification handler code hereUpdateData();if(m_tuoqiu=-1)MessageBox(请选择你需要的椭球！);/return -1;double e2; double c;switch(m_tuoqiu)case 0:e2=0.006738525414683; /第二偏心率 c=6399698.901782711; /极点处的子午线曲率半径break;case 1:e2=0.006739501819473; /第二偏心率 c=6399596.6519880105; /极点处的子午线曲率半径break;case 2:e2=0.00673949674227; /第二偏心率 c=6399593.6258; /极点处的子午线曲率半径break;case 3: e2=0.00673949677548; /第二偏心率 c=6399593.6259; /极点处的子午线曲率半径 break;default:break;double X=0; /自赤道起的子午线弧长double dX; /微分 double M; /子午圈曲率半径double N; /卯酉圈曲率半径double pi; /double l; /距中央子午线的经度差（以秒为单位）double t; /B的tandouble yita;double rou; /以秒为单位double x,y; /坐标 int i,j=200000; /微分数double B1=0; /赤道的纬度 double B,m_B1,m_B2,m_B3;double dB;double m_L1,m_L2,m_L3;double L0; /中央子午线经度char ch130,ch230,ch330,ch430,ch530,ch630,ch730; CString str8,str9; GetDlgItem(IDC_EDIT1)-GetWindowText(ch1,30);GetDlgItem(IDC_EDIT2)-GetWindowText(ch2,30);GetDlgItem(IDC_EDIT3)-GetWindowText(ch3,30);GetDlgItem(IDC_EDIT4)-GetWindowText(ch4,30);GetDlgItem(IDC_EDIT5)-GetWindowText(ch5,30);GetDlgItem(IDC_EDIT6)-GetWindowText(ch6,30);GetDlgItem(IDC_EDIT7)-GetWindowText(ch7,30);m_B1=atof(ch1);m_B2=atof(ch2);m_B3=atof(ch3);m_L1=atof(ch4);m_L2=atof(ch5);m_L3=atof(ch6);L0=atof(ch7);l=(m_L1-L0)*3600+m_L2*60+m_L3;/角度与弧度的转化pi=4*atan(1);B=pi*(m_B1+m_B2/60.0+m_B3/3600.0)/180.0;rou=180*3600/pi;/求rou/积分得XdB=B/j;for (i=0;iSetWindowText(str8); GetDlgItem(IDC_EDIT9)-SetWindowText(str9);高斯投影反算/ FansuanDlg.cpp : implementation file/#include stdafx.h#include gaosi.h#include FansuanDlg.h#include math.h#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CFansuanDlg dialogCFansuanDlg:CFansuanDlg(CWnd* pParent /*=NULL*/): CDialog(CFansuanDlg:IDD, pParent)/AFX_DATA_INIT(CFansuanDlg)x = 0.0; y = 0.0;L0 = 0;m_B1 = 0;m_B2 = 0;m_B3 = 0.0;m_l1 = 0;m_l2 = 0;m_l3 = 0.0;m_tuoqiu = -1;/AFX_DATA_INITvoid CFansuanDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);AFX_DATA_MAP(CFansuanDlg)DDX_Text(pDX, IDC_EDIT1, x);DDX_Text(pDX, IDC_EDIT2, y);DDX_Text(pDX, IDC_EDIT3, L0);DDX_Text(pDX, IDC_EDIT4, m_B1);DDX_Text(pDX, IDC_EDIT5, m_B2);DDX_Text(pDX, IDC_EDIT6, m_B3); DDX_Text(pDX, IDC_EDIT7, m_l1);DDX_Text(pDX, IDC_EDIT8, m_l2);DDX_Text(pDX, IDC_EDIT9, m_l3);DDX_Radio(pDX, IDC_RADIO1, m_tuoqiu);/AFX_DATA_MAPBEGIN_MESSAGE_MAP(CFansuanDlg, CDialog)/AFX_MSG_MAP(CFansuanDlg)ON_BN_CLICKED(IDC_BUTTON1, OnButton1)/AFX_MSG_MAPEND_MESSAGE_MAP()/ CFansuanDlg message handlersvoid CFansuanDlg:OnButton1() / TODO: Add your control notification handler code hereUpdateData();if(m_tuoqiu=-1)MessageBox(请选择你需要的椭球！);/return -1;double e2; double c;switch(m_tuoqiu)case 0:e2=0.006738525414683; /第二偏心率 c=6399698.901782711; /极点处的子午线曲率半径break;case 1:e2=0.006739501819473; /第二偏心率 c=6399596.6519880105; /极点处的子午线曲率半径break;case 2:e2=0.00673949674227; /第二偏心率 c=6399593.6258; /极点处的子午线曲率半径break;case 3: e2=0.00673949677548; /第二偏心率 c=6399593.6259; /极点处的子午线曲率半径 break;default:break;double x,y; /坐标double Bf=0; /地点纬度double dBf; double pi=4*atan(1);double dx;double Mf; /底点纬度对应的子午圈曲率半径double Nf; /底点纬度对应的平行圈曲率半径double tf; /tan(Bf)double yitaf;double B,l; /纬度和经差int m_B1,m_B2; /B的度和分double m_B3; /B的秒int m_l1,m_l2; /l的度和分double m_l3; /l的秒int L0; /中央子午线经度 int i,j=10000000; /L0=111;/x=3404139.8389;/y=528680.5710; y=y-500000; dx=x/j;for (i=0;ij;i+)Mf=c/pow(sqrt(1+e2*pow(cos(Bf),2),3);dBf=dx/Mf;Bf=Bf+dBf;Mf=c/pow(sqrt(1+e2*pow(cos(Bf),2),3);Nf=c/sqrt(1+e2*pow(cos(B),2); tf=tan(Bf);yitaf=sqrt(e2)*cos(Bf);B=Bf-(tf*y*y)/(2*Mf*Nf)+(tf*(5+3*tf*tf+pow(yitaf,2)-9*yitaf*yitaf*tf*tf)*pow(y,4)/(24*Mf*pow(Nf,3)-(tf*(61+90*tf*tf+45*pow(tf,4)*pow(y,6)/(720*Mf*pow(Nf,5);l=y/(Nf*cos(Bf)-(1+2*tf*tf+yitaf*yitaf)*pow(y,3)/(6*pow(Nf,3)*cos(Bf)+(5+28*tf*tf+24*pow(tf,4)+6*pow(yitaf,2)+8*pow(yitaf,2)*pow(tf,2)*pow(y,5)/(120*pow(Nf,5)*cos(Bf);/弧度化成角度m_B1=int(B*180)/pi);m_B2=int(B*180*60)/pi-m_B1*60);m_B3=(B*180*3600)/pi-m_B1*3600-m_B2*60;m_l1=int(l*180)/pi);m_l1=m_l1+L0;m_l2=int(l*180*60)/pi-m_l1*60);m_l3=(l*180*3600)/pi-m_l1*3600-m_l2*60实测斜距化算至高斯平面/ XiejugaihuaDlg.cpp : implementation file/#include stdafx.h#include 大地课程设计.h#include XiejugaihuaDlg.h#include string.h#include math.h#define pi 3.1415926535897932384626433#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CXiejugaihuaDlg dialogCXiejugaihuaDlg:CXiejugaihuaDlg(CWnd* pParent /*=NULL*/): CDialog(CXiejugaihuaDlg:IDD, pParent)/AFX_DATA_INIT(CXiejugaihuaDlg)x1 = 0.0;y1 = 0.0;m_b1du = 0.0;m_b1fen = 0.0;m_b1miao = 0.0;H1 = 0;x2 = 0.0;y2 = 0.0;m_b2du = 0.0;m_b2fen = 0.0;m_b2miao = 0.0;H2 = 0;D = 0.0;yichang = 0;Dp = 0.0;m_tuoqiu = -1;/AFX_DATA_INITvoid CXiejugaihuaDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CXiejugaihuaDlg)DDX_Text(pDX, IDC_EDIT1, x1);DDX_Text(pDX, IDC_EDIT2, y1);DDX_Text(pDX, IDC_EDIT3, m_b1du);DDX_Text(pDX, IDC_EDIT4, m_b1fen);DDX_Text(pDX, IDC_EDIT5, m_b1miao);DDX_Text(pDX, IDC_EDIT6, H1);DDX_Text(pDX, IDC_EDIT7, x2);DDX_Text(pDX, IDC_EDIT8, y2);DDX_Text(pDX, IDC_EDIT9, m_b2du);DDX_Text(pDX, IDC_EDIT10, m_b2fen);DDX_Text(pDX, IDC_EDIT11, m_b2miao);DDX_Text(pDX, IDC_EDIT12, H2);DDX_Text(pDX, IDC_EDIT13, D);DDX_Text(pDX, IDC_EDIT14, yichang);DDX_Text(pDX, IDC_EDIT15, Dp);DDX_Radio(pDX, ID