面向对象程序设计——陈士城-副本.docVIP

面向对象程序设计 开发环境：VC++6.0。 主要工作：在测量中会涉及到几种坐标系统的转换，本程序实现了一下几类坐标系的转换（包括：北京-54坐标，国家-80坐标，WGS-84坐标）。 大地坐标与高斯平面坐标的转换。 将大地坐标转为高斯平面坐标，按照高斯投影正算公式进行。将高斯平面坐标转为大地坐标，按照高斯投影反算公式进行。 （1）、由大地坐标求高斯平面坐标： 高斯投影坐标正算公式（由求）： （2）、由高斯平面坐标求大地坐标： 高斯投影坐标反算公式（由求）： 大地坐标与空间直角坐标坐标的转换。 （1）、由大地坐标求空间直角坐标： （2）、由空间直角坐标求大地坐标： 式中，,为该点的卯酉圈曲率半径；，、分别为该大地坐标系对应椭球的长半轴和第一偏心率。 空间直角坐标与空间直角坐标的转换。 应用的是布尔莎公式： 公式存在7个参数：3个平移参数，3个旋转参数，1个尺度变化参数。附1 程序设计界面 附2 主要程序代码 // DlgCoord.cpp : implementation file ///////////////////////////////////////////////////////////////////////////// BOOL CDlgCoord::OnInitDialog() { CDialog::OnInitDialog(); m_comOld.SetCurSel(0); //默认椭球类型 m_comNew.SetCurSel(0); CWnd* pWnd; pWnd=GetDlgItem(IDC_COMBO_GEO_T0); pWnd-ShowWindow(SW_HIDE); pWnd=GetDlgItem(IDC_COMBO_GEO_T1); pWnd-ShowWindow(SW_HIDE); return TRUE; // return TRUE unless you set the focus to a control // EXCEPTION: OCX Property Pages should return FALSE } void CDlgCoord::OnRadioGeoOld() //源 大地坐标 { CWnd* pWnd; pWnd=GetDlgItem(IDC_COMBO_GEO_T0); pWnd-ShowWindow(SW_SHOW); m_comGeoOld.SetCurSel(0); m_strX0=; //置0 m_strY0=; m_strH0=; pWnd=GetDlgItem(IDC_STATIC_X0);//标签标题改变 pWnd-SetWindowText(B=); pWnd=GetDlgItem(IDC_STATIC_Y0); pWnd-SetWindowText(L=); pWnd=GetDlgItem(IDC_STATIC_H0); pWnd-SetWindowText(H=); } …… 其余几个单选按钮类似 void CDlgCoord::OnButtonCoordOk() //坐标转换 按钮 { UpdateData(); //高斯正算 大地-平面 if(m_CoordTypeOld==1m_CoordTypeNew==2) { if(m_strX0.GetLength()==0||m_strY0.GetLength()==0) {AfxMessageBox(请输入数据！);m_ctrX0.SetFocus();} else GaussForward(); } //高斯反算 平面-大地 …… //大地坐标-空间直角坐标 …… //空间直角坐标-大地坐标 …… //空间直角坐标-空间直角坐标 …… } void CDlgCoord::GaussForward()//高斯正算 { UpdateData(); double B,Bbak,L,l,ll,L0,x,y,cos2B,p2; double N,a0,a4,a6,a3,a5; p2=206264.806247; B=stod(m_strX0); // 化为 秒 L=stod(m_strY0); //化为 秒 Bbak=B; //备份B B=B/p2; //化为弧度值 L0=stod(g_center);//获得中央子午线 // AfxMessageBox(g_center); l=L-L0; //化为 秒 l=l/p2;//化为弧度 ll=l*l; //l相乘 cos2B=cos(B)*cos(B); //国家-80 大地坐标