常用测量程序设计.docVIP

《测绘程序设计（VC.net）》 上机实验报告 （Visual C++.Net） 班 级： 学 号： 姓 名： 序 号： 二零一零年四月 实验7 常用测量程序设计 一、实验目的 巩固类的创建与使用 掌握数组参数的传递 掌握常用测绘程序设计的技巧 二、实验内容 编写高斯投影正、反算程序 设计思路 基于对话框的设计，添加两个头文件分别是Angle和Gauss ,在角度类中添加角度类的函数，角度类的函数我用了老师的例子，因为我以前写的那个角度类，定义的存贮角度类型的值是一个变量，引用起来不好用。然后就是在Gauss的头文件和cpp中添加代码，我一共定义了四个函数，分别是高斯正算，求X的一个函数，高斯反算，求Bf的一个函数。定义完类之后，就在界面上引用函数即可。 主要代码 Gauss.h #pragma once class Gauss { /*double a; double b;*/ //double e=sqrt(a*a-b*b)/a; //double _e=sqrt(a*a-b*b)/b; public: /*Gauss(double a1,double b1);*/ Gauss(void); ~Gauss(void); public: void CGaussForward(double dB,double dL,double dCML,int type,double x,double y); double ComptX(double dB); void CGaussInverse(double x, double y, double dCML,int type, double dB, double dL); double ComptBf(double x); void ellipd(int type); }; Gauss.cpp #include StdAfx.h #include Gauss.h #include math.h #include Angle.h const double PI=4.0*atan(1.0); double a=6378245.0; double b=6356863.0187730473;//初始化 double e=sqrt(a*a-b*b)/a; double _e=sqrt(a*a-b*b)/b; /*Gauss::Gauss(double a1=6378245.0,double b1=6356863.0187730473) { a=a1; b=b1; }*/ //Gauss::Gauss(double a1, double b1) //{ a1=6378245.0; // b1=6356863.0187730473; // a=a1; // b=b1; // /*double e=sqrt(a*a-b*b)/a; //double _e=sqrt(a*a-b*b)/b;*/ //} Gauss::Gauss(void) { } Gauss::~Gauss(void) { } void Gauss::ellipd(int type) { if(type==1) //克拉索夫斯基椭球 { a=6378245.0; b=6356863.0187730473; //e=sqrt(0.006693421622966); } else if(type==2) //1975国际协议椭球 { a=6378140.0; b=6356755.2881575287; //e=sqrt(0.006694384999588); } else if(type==3) //WGS-84椭球 { a=6378137.0; b=6356752.3142; //e=sqrt(0.0066943799013); } } void Gauss::CGaussForward(double dB,double dL,double dCML,int type,double x,double y) { CAngle B(dB); CAngle L(dL); CAngle CML(dCML); CAngle _L(L-CML); //计算经差 dB=B(RAD); ellipd(type);//设置椭球参数 //double fengdai(type);//投影带 double l=_L(RAD); double N=a/sqrt