测绘程序实验报告

实验1 Visual C++.Net环境和程序设计初步1.掌握 VC++ 语言的基本语法；2.理解顺序结构、选择结构和循环结构程序设计的特点及应用；3.掌握对基于对话框的 MFC 应用程序设计方法；4.掌握一些简单算法。

5.编写一个方位角计算程序。

提示：先使用反正切函数计算，然后利用坐标增量的符号来判断所在的象限。

设计思路：在按钮下面添加程序。

X Y同时大于0在第一象限，方位角等于arctan(y/x)；X>0,Y<0在第二象限，方位角等于arctan(y/x)+90；X<0,Y<0在第三象限，方位角等于arctan(y/x)+180；X<0,Y>0在第四象限，方位角等于arctan(y/x)+270；界面设计：主要代码：计同一参考椭球下的三维地心坐标（笛卡儿坐标系）与大地坐标系转换的程序。

（提示：用 do…while 迭代，B、H 初始为 0 进行迭代，直到 H 的精度达到米）注意：东经 0~180（Y>0）,西经：0~-180(Y<0)式中，B、L、H 为椭球面上的大地纬度、大地经度、大地高；X、Y、Z 为空间直角坐标；N 为卯酉圈曲率半径，e 为椭球的偏心率，a 为椭球的长半径，b 为椭球的短半径。

（WGS84 椭球参数：长半径 a=6378137m，扁率α=1/）设计思路：在按钮下面设置主程序，按照指导书给的思路编辑公式界面设计：主要代码：写一个后方交会计算程序。

基本原理及计算公式若将 Pa、Pb、Pc 看成权，则 P 点的坐标即为三个已知点的加权平均值计算程序设计步骤（1）设计界面，用于输入 3 个已知点的坐标和三个观测角、和，以及用于输出待定点坐标的文本框（12 个）、静态标签框和 Button 按钮；（2）定义文本框控件变量（Value）；（3）根据已知点计算三个内角 A、B、C；（4）计算 Tan(α)、Tan(β)、Tan(γ)、Tan(A)、Tan(B)、Tan(C)；（5）计算 Pa、Pb、Pc；（6）计算待定点坐标 Xp、Yp。

界面要求:三个坐标输入框，两个角度输入框一个坐标结果输出框一个计算按钮，一个清除按钮，一个退出按钮。

2.由三角形三个边长求内角函数计算公式：设计思路：通过示例编辑框添加变量，在按钮下面添加程序。

先将α、β、的度分秒之转换成度利用三角形内角和等于180算出γ并将γ转换成弧度。

接着计算三角形的内角，判断P点是否在危险圆上若不在则计算P点坐标若在则弹出“该点在危险圆上”。

界面设计：主要代码：protected:HICON m_hIcon;写一个求任意多边形面积的程序。

提示：通过界面输入数据，并把数据保存在一个二维数组或一个一维的自定义结构体类型的数组中，然后再进行计算。

要求计算部分写成函数的形式，使计算程序与界面无关。

动态数组创建动态数组结构体的定义多边形面积计算原理及算法计算原理：面积计算的算法：提示：显示框用 Cedit 控件变量对每个输入的坐标用 CString str 临时变量格式化，然后用 CEdit 的控件变量插入设计思路：将实现写在按钮在下面。

添加顶点下面实现创建动态数组，确认按钮实现返回数组大小值并将添加的数据显示到显示框中去，通过计算按钮实现多边形面积的计算，清除按钮实现清除输入的数据，退出按钮退出程序。

界面设计：主要代码："#include"05(3)"#ifdef _DEBUG#define new DEBUG_NEW#endifd % lf %lf \r\n"),t,x,y);(str);px[t-1]=x;py[t-1]=y;UpdateData(FALSE);}void CMy073Dlg::OnBnClickedButton5(){掌握面向对象编程基本思想2. 掌握 VC++ 中创建类3. 掌握建立和使用对象4. 掌握运算符号重载5. 理解类的继承和多态性二、实验内容1.设计一个角度类。

要求该类具有度分秒至度的换算、度至度分秒的换算、度与弧度的换算等功能。

提示：设置一个角度大小属性，并设定该属性为缺省属性；另设一个状态属性，表示当前设置的角度大小的形式；度分秒、度、弧度间的相互转换的方法；定义运算符号（加、减）方法，使得角度类能够像一种普通的数据类型样的方便使用。

设计思路：设计一个角度类在类的头文件里申明度到度分秒的转换、度分秒到度的转换、度到弧度的转换、运算符的重载，在类的.cpp文件中写出具体的函数实现，在按钮下面调用类的各个函数来实现角度的转换。

界面设计：主要代码：DECLARE_MESSAGE_MAP()public:double dmstodgree(double dmg);编制简单的 Cass 数据文件进行数据整理的程序。

整理后的数据文件中要求无重复点数据，且数据按点号大小的升序进行排序。

要求整理后的数据按与原始数据文件同样的格式保存为另外一个文件示例数据“民用园燃气.dat”文件数据格式：总点数点号，编码，X,Y,H…例如：要求：a.用SaveFileDialog和OpenFileDialog控件获取文件打开或保存的文件名。

b.自定义一个测量点数据结构体，其元素包括：点号，编码，X,Y,Hc.用文本框显示原始数据和整理后的数据d..按编程规范进行编码a. 读数据、保存数据、判断一个点数据是否已经存在、排序、在文本框中显示文件内容等可以分别定义成一个子过程或函数。

c.编程技巧打开文件获得总点数根据点数调整数组的大逐行读取数据，判断该点是否已经存在数组中，若不存在，则把该数据存放在数组中根据删除重复点后的总点数，重新调整数组大小按点号大小的升序排序按原数据格式输出到另外一个文件设计思路：在按钮下面设置主程序。

通过打开文件按钮将原文件打开并显示到示例编辑框内，通过保存文件按钮实现对文件分行、删除重复点、排序、保存更改后的文件并把它显示到示例编辑框中等操作。

界面设计：主要代码：xt"));at"));umber=_ttoi(strTmp[0]);oder=strTmp[1];pPoint[i].X=_tstof(strTmp[2]);pPoint[i].Y=_tstof(strTmp[3]);pPoint[i].H=_tstof(strTmp[4]);}if(strTmp!=NULL)umber=pPoint[j].number)delete [] pPoint;umber>pPoint[j].number){t=pPoint[i].number;pPoint[j].number=t;pPoint[i].number=pPoint[j].number;}}}umber,pPoint[i].coder,pPoint[i].X,pPoint[i].Y,pPoint[i].H);}xiugai=xiugai+strOutPut;列数据为一变形监测点的 24 期位移监测结果（分别为X,Y,H），编制程序绘制出该点的变形曲线图，每个方向一个位移序列图。

设计思路：用类DeformationCurve实现原始文件的读入、文件的拆分、画图等功能，调用类的相关函数来实现相应的功能。

界面设计：因为是基于当文本的所以没有界面主要代码：class DeformationCurve : public CWnd{DECLARE_DYNAMIC(DeformationCurve)public:DeformationCurve();virtual ~DeformationCurve();protected:DECLARE_MESSAGE_MAP()public:CString *DeformationCurve::Splitstring(CString str,char split,int &isubstr);void DeformationCurve::ReadData(double *&X,double *&Y,double *&H);void DeformationCurve::Draw(CDC* pDC,CRect& rect);};xt"));f"), pDC->TextOutW(dOrgX-60,dEndY+i*dy,str);}for(int i=0;i<=3;i++){(_T("%.3f"), pDC->TextOutW(dOrgX-120,dEndY+i*dy,str);}for(int i=0;i<=3;i++){(_T("%.3f"), pDC->TextOutW(dOrgX=180,dEndY+i*dy,str);}();=dOrgX+i*dx;ptX[i].y=dOrgY-int((X[i]/*3*dy);}pDC->Polyline(ptX,24);=dOrgX+i*dx;ptY[i].y=dOrgY-int((Y[i]+/*3*dy);}pDC->Polyline(ptY,24);=dOrgX+i*dx;ptH[i].y=dOrgY-int((H[i]*3*dy));}pDC->Polyline(ptH,24);trName=pstrData[0];m_pKnownPoint[i].strID=pstrData[0];m_pKnownPoint[i].H=_tstof(pstrData[1]);m_pKnownPoint[i].flag=1;trName=pstrData[i];m_pUnknownPoint[i].strID=pstrData[i];m_pUnknownPoint[i].H=0;lag=0;pBackObj=SearchPointUsingID(pstrData[0]);pF rontObj=SearchPointUsingID(pstrData[1]);ObsValue=_tstof(pstrData[2]);ist=_t stof(pstrData[3]);trID==ID){return &m_pKnownPoint[i];}}return NULL;}trID==ID){return &m_pUnknownPoint[i];}}return NULL;}LevelControlPoint* AdjustLevel::SearchPointUsingID(CString ID){LevelControlPoint* cp;cp=SearchKnownPointUsingID(ID);if(cp==NULL)cp=SearchUnknownPointUsingID(ID);return cp;}void AdjustLevel::ApproHeignt(void)lag!=1){pFrontObj->strID==m_pUnknownPoint[i].strID)&& m_pDhObs[j].cpBackObj->flag==1 ){ =m_pDhObs[i].cpBackObj->H - m_pDhObs[i].ObsValue;*/m_pUnknownPoint[i].H=m_pDhObs[j].cpBackObj->H +m_pDhObs[j].HObsValue;m_pUnknownPoint[i].flag=1;break;}}if(m_pUnknownPoint[i].flag!=1)pBackObj->strID==m_pUnknownPoint[i].strID) && m_pDhObs[j].cpFrontObj->flag==1 ){ =m_pDhObs[j].cpFrontObj->H-m_pDhObs[j].HObsValue;/*m_pUnknownPoint[i].H=m_pDhObs[i].cpFrontObj->H+m_pDhObs[i].ObsValue;*/m_pUnknownPoint[i].flag=1;break;}}}}if(i==m_iUnknownPointCount-1)lag!=1)ist);p(i,i)=value;}return p;}void AdjustLevel::FormErrorEquation(CMatrix &B, CMatrix &L){(m_iDhObsCount,m_iUnknownPointCount);(m_iDhObsCount,1);for(int i=0;i<m_iDhObsCount;i++)pBackObj->strID);tmpFront=SearchPointUsingID(m_pDhObs[i].cpFrontObj->strID);trID==tmpBack->strID)trID==tmpFront->strID)bsValue-(m_pDhObs[i].cpBackOb j->H-m_pDhObs[i].cpFrontObj->H);*/L(i,0)=m_pDhObs[i].HObsValue-(m_pDhObs[i].cpFrontObj->H -m_pDhObs[i].cpBackObj->H);(_T("%.3f"),L(i,0));L(i,0)=_tstof(tmp);L(i,0)=L(i,0)*1000;+=x(i,0);xt"));xt"));if()==IDCANCEL) return;CString strFileName=();setlocale(LC_ALL,"");CStdioFile sf;if(!(strFileName, CFile::modeCreate|CFile::modeWrite)) return;(LevleContent);();UpdateData(FALSE);}void CIndircLelveDlg::OnBnClickedComputelevel(){f\r\n"), [i].strID,[i].H);LevleContent+=Temp;}(_T("单位权中误差：%.1f mm\r\n"),r0*1000);LevleContent+=Temp;LevleContent+=_T("未知点高程中误差(mm):\r\n");for(int i=0;i< ;i++){();(_T("%s,%.1f\r\n"),[i].strName,Qx[i]*1000);LevleContent+=Temp;}UpdateData(false);}void CIndircLelveDlg::OnBnClickedSavelevleresult(){xt"));if()==IDCANCEL) return;CString strFileName=();setlocale(LC_ALL,"");CStdioFile sf;if(!(strFileName, CFile::modeCreate|CFile::modeWrite)) return;(LevleContent);();UpdateData(FALSE);}二、实验结果打开文件数据：平差结果：四、实验心得这从实验是我们测绘程序设计的最后一次实验，虽然这个学期我们做了好几次相关的实验，但是我却发现自己学的东西也越来越模糊，感觉很多内容都不理解。