测绘程序设计(VS2008)实验报告--高斯投影正反算.docx

《测绘程序设计（VC.net）》上机实验报告（Visual C++.Net）班级：学号：姓名：序号：二零一一年五月实验7常用测量程序设计实验目的：巩固类的创建与使用;掌握数组参数的传递；掌握常用测绘程序设计的技巧。实验内容：编写高斯投影正反算程序。设计思路：这次的实验目的是实现高斯正反算。需要考虑以下几个问题投影方式即分带的方式，又要考虑椭球参数的类型，可通过单选组框和相应函数来完成此功能，分别是int SetProjectType(int m)和void SetParameter(int m,double a,double b)界面设计界面设计很简单，如下：其中每个复选框的第一个选项的属性中Group属性设置为TRUE主要代码：文件名：Ccehui1201_shiyan7_yang_GaussDlg.hvoid Ccehui1201_shiyan7_yang_GaussDlg::SetParameter(double a,double b);//设置椭球参数double Ccehui1201_shiyan7_yang_GaussDlg::SetProjectType(double LL);//设置分带方式double Deg(double dDms);文件名：Ccehui1201_shiyan7_yang_GaussDlg.cpp代码：constdouble PI=4*atan(1.0);在BOOL Ccehui1201_shiyan7_yang_GaussDlg::OnInitDialog()添加代码：《CheckRadioButton(IDC_RADIO1,IDC_RADIO2,IDC_RADIO3);CheckRadioButton(IDC_RADIO4,IDC_RADIO5,IDC_RADIO6);》//获得分带方式返回中央子午线经度double Ccehui1201_shiyan7_yang_GaussDlg::SetProjectType(double LL){UpdateData(TRUE);int n; //记录分带带号double L2;//经度L2=Deg(LL);if(IsDlgButtonChecked(IDC_RADIO1)) //6度带{n=int(L2/6)+1;L0=6*n-3;}elseif(IsDlgButtonChecked(IDC_RADIO2)) //3度带{n=int((L2+1.5)/3);L0=3*n;} elseif(IsDlgButtonChecked(IDC_RADIO3)) //自主分带L0=L0; return L0;}void Ccehui1201_shiyan7_yang_GaussDlg::SetParameter(double a,double b){if(IsDlgButtonChecked(IDC_RADIO4)) //克拉索夫斯基椭球{a=6378245.0;b=6356863.0187730473;//e=sqrt(0.006693421622966);}elseif(IsDlgButtonChecked(IDC_RADIO5)) //1975国际协议椭球{a=6378140.0;b=6356755.2881575287;//e=sqrt(0.006694384999588);} elseif(IsDlgButtonChecked(IDC_RADIO6)) //WGS-84椭球{a=6378137.0;b=6356752.3142;//e=sqrt(0.0066943799013);} }//正算按钮void Ccehui1201_shiyan7_yang_GaussDlg::OnBnClickedOk(){// TODO: 在此添加控件通知处理程序代码double N;//卯酉圈曲率半径double t;doubleη;double X; //当前点到赤道的子午线长度double A0,A2,A4,A6,A8;//计算X的参数double RadB;doubleρ;ρ=180*3600/PI;//定义弧秒值，近似于″double a,b,e1,e2; //椭球参数SetParameter(a,b);e1=sqrt(a*a-b*b)/a;e2=sqrt(a*a-b*b)/b;double l;//经度L同当前带中央子午线的差值，单位为秒L0=SetProjectType(L);double RadL;RadL= Deg(L);l=(L-L0)*3600;RadB= Deg(B)*PI/180;N=a/sqrt(1-e1*e1*sin(RadB)*sin(RadB));t=tan(RadB);η=e2*cos(RadB);A0=1+