测绘程序实验报告.docx

专业 专业.专注 专业 专业.专注 H H 0145110615 ymh 2,1 | 昜 | E I 坐标 I i x- 示例编辑柱 i I F E i I Y= 示例编辑框 I E i S E E I | i a 象限 示例漏鉅框 方位角:示例前框: I ■ ■ ■ E I 3 S E E E I ~~| SS~~] i : * ■ 实验1 Visual C++.Net 环境和程序设计初步 1.掌握VC++.net语言的基本语法; 2?理解顺序结构、选择结构和循环结构程序设计的特点及应用； 掌握对基于对话框的 MFC应用程序设计方法； 掌握一些简单算法。 编写一个方位角计算程序。提示：先使用反正切函数计算，然后利用坐标增量 的符号来判断所在的象限。 设计思路：在按钮下面添加程序。XY同时大于0在第一象限，方位角等于 arctan(y/x) ； X0,Y0 在第二象限，方位角等于 arctan(y/x)+90 ； X0,Y0 在第三象 限，方位角等于arctan(y/x)+180 ； X0,Y0在第四象限，方位角等于 arctan(y/x)+270 ； 界面设 计： 主要代码：/ 0145110615 ymh 2.1Dlg.h : 头文件 protected : HICON m_hlco n; //生成的消息映射函数 virtual BOOL OnlnitDialog(); afx_msg void On SysComma nd(UINT nID, LPARAM lParam); afx_msg void On Pai nt(); afx_msg HCURSOR On QueryDragIco n(); DECLARE_MESSAGE_MAP() public : double x; double y; double Q; afx_msg void OnBn ClickedOk(); afx_msg void OnBn ClickedCa ncel(); afx_msg void OnBn ClickedButt on 1(); double A; }; // 0145110615 ymh 2.1Dlg.cpp : 实现文件 #include stdafx.h #include 0145110615 ymh 2.1.h #include 0145110615 ymh 2.1Dlg.h #inelude math.h #ifdef _DEBUG #define new DEBUG_NEW #en dif void CMy0145110615ymh21Dlg::0 nBn ClickedOk() // 计算 { // TODO:在此添加控件通知处理程序代码 UpdateData(TRUE); if(x0) {if(y0) {Q=1; A=atan(y/x);} 〃x大于y大于在第一象限 else{Q=4; A=atan(y/x)+270;} //x大于y小于在第四象限 } else if(y0) {Q=2; A=atan(y/x)+90;} //x大于y小于在第二象限 else{Q=3; A=atan(y/x)+180;} //x小于y小于在第三象限 UpdateData(FALSE); //On OK(); } void CMy0145110615ymh21Dlg::O nBn ClickedCa ncel() { // TODO:在此添加控件通知处理程序代码 UpdateData( true); x=0; y=0; Q=0; A=0; UpdateData( false); } void CMy0145110615ymh21Dlg::O nBn ClickedButto n1() { // TODO:在此添加控件通知处理程序代码 //清除//退出 //清除 //退出 } 运行结 果： 第一次做这个实验的时候真的觉得挺难得可是当我经过更难的实验的洗礼之后在返 回来看它真是顿时亲切了许多，其实当时是对这个软件不了解简单的东西把它想得很 复杂绕来绕去都不知道该怎么实现一些很简单的东西了 。 实验2.2 .设计同一参考椭球下的三维地心坐标 （笛卡儿坐标系）与大地坐标系转换的 程 序。 （提示：用do…while迭代，B、H初始为0进行迭代，直到H的精度达到 0.00001 米） 注意：东经 0~180 （ Y0），西经：0~-180（丫0） 式中，B、L、H为椭球面上的大地纬度、大地经度、大地高；X、丫、Z为空间直 角坐标；N为卯酉圈曲率半径，e为椭球的偏心率，a为椭球的长半径，b为椭球 的短半径。（WGS84 椭球参数：长半径 a=6378137m ，扁率a 的短半径。 =1/298.257223563) 设计思