实验三 高斯投正算.docVIP

实验三 高斯投影正算 地信050212 赵静 实验目的 1、高斯-克吕格投影正算的编程计算。 2、熟悉C++编程练习。 实验内容 设计一个高斯投影正算的VC程序，要求有单点和多点文件的输入，输出结果包含列表、文件及绘图输出，并且要求计算结果正确。 实验步骤 利用VC建立工程（MFC appwizard（exe））：Test_Gause 设计对话框窗体布局。如下图所示： 在classwizard中增加变量及消息映射函数。 设计选择椭球函数 /******************************************** 功能：选择计算所要用的椭球和投影带 共三个椭球，两种投影带 ********************************************/ void CTest_GauseDlg::choose() { this-UpdateData(); int m_index; int m_daiindex; m_index=((CComboBox*)GetDlgItem(IDC_COMBO_TUOQIU))-GetCurSel(); m_daiindex=((CComboBox*)GetDlgItem(IDC_COMBO_TUOYINGDAI))-GetCurSel(); flag_daicha=3*m_daiindex+3; switch (m_index) { case 0: //克氏 a0 = 111134.8610828; a2 = -16036.48022; a4 = 16.82805; a6 = -2.197E-02; a8 = 3E-05; c = 6399698.901782711; a = 6378245; square_e1 = 6.738525414683E-03; square_e = 6.693421622966E-03; q0 = 157046064.12328E-15; q2 = 25258869461.858E-12; q4 = 6572E-12; q6 = 120717.4265E-12; q8 = -1075.1509E-12; break; case 1: //1975 a0 = 111133.0046793; a2 = -16038.52818; a4 = 16.83263; a6 = -2.198E-02; a8 = 3E-05; c = 6399596.6519880105; a = 6378140; square_e1 = 6.739501819473E-03; square_e = 6.694384999588E-03; q0 = 157048687.47416E-15; q2 = 2526252791.9786E-12; q4 = 4356E-12; q6 = 120769.9608E-12; q8 = -1075.7700E-12; break; case 2: //wgs a0 = 111132.9525494; a2 = -16038.50840; a4 = 16.83260; a6 = -2.198E-02; a8 = 3E-05; c = 6399593.6258; a = 6378137; square_e1 = 6.73949674227E-03; square_e = 6.6943799013E-03; q0 = 157048761.142065E-15; q2 = 2526250855.8867E-12; q4 = 5362E-12; q6 = 120769.6828E-12; q8 = -1075.7667E-12; break; default: break; } } 正算公式设计 /******************************************** 功能：进行高斯投影正算并输出结果 ********************************************/ void CTest_GauseDlg::OnButtonZcount() { // TODO: Add your control notification handler code here this-choose(); double angle_to_ra