[测绘程序设计_01]调MFC实现地心坐标与大地坐标相互转换界面

一、界面设计

二、参数定义与界面控件

（1）
const double PI = 4 * atan(1.0);//声明π；
const double a = 6378137;//椭球长半径
const double α = 1 / 298.257223563;//椭球扁率
const double e_2 = 2 * α - α * α;//第一偏心率

三、代码

（2）
//大地坐标转笛卡尔坐标主程序
void CMFCExp4Dlg::OnBnClickedGeo2car()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData(TRUE);
	Geo2Car(L, B, H, X, Y, Z);
	UpdateData(FALSE);
}

（3）
//笛卡尔坐标转大地坐标主程序
void CMFCExp4Dlg::OnBnClickedCar2geo()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData(TRUE);
	Car2Geo(X, Y, Z, L, B, H);
	UpdateData(FALSE);
}

（4）
//大地坐标转笛卡尔坐标函数
void CMFCExp4Dlg::Geo2Car(double L,double B,double H,double &X,double &Y,double &Z)
{
	L  = Dms2Deg(L);
	B = Dms2Deg(B);
	double N = a / (sqrt(1 - e_2 *sin(B)*sin(B)));
	X = (N + H)*cos(B)*cos(L);
	Y = (N + H)*cos(B)*sin(L);
	Z = (N*(1 - e_2) + H)*sin(B);
}

（5）
//笛卡尔坐标转大地坐标函数
void CMFCExp4Dlg::Car2Geo(double X,double Y,double Z,double &L,double &B,double &H )
{
	H = 0; // 设定高程初始值进行迭代
	double dH = 0;
	double N, pcs;
	B = atan(Z / (sqrt(X*X + Y * Y)*(1 - e_2)));
	do // do...while语句迭代计算
	{
		N = a / sqrt(1 - e_2*sin(B)*sin(B));
		B = atan(Z / (sqrt(X*X + Y * Y)*(1 - e_2*N / (N + H)))); 
		H = sqrt(X*X + Y * Y) / cos(B) - N;
		pcs = H - dH;
		dH = H;
	} while (fabs(pcs) >= 1e-5); // 判断精度是否符合要求
	L = atan2(Y , X);//大地经度

	//弧度转度分秒输出
	L  = Deg2Dms(L);
	B = Deg2Dms(B);
}

（6）
//清空按钮，初始化操作
void CMFCExp4Dlg::OnBnClickedClear()
{
	// TODO: 在此添加控件通知处理程序代码
	L = 0;
	B = 0;
	H = 0;
	X = 0;
	Y = 0;
	Z = 0;
	UpdateData(FALSE);
}

（7）
//度分秒转度
double CMFCExp4Dlg::Dms2Deg(double dDms)
{
	//用于存放度、分、秒三个值的变量
	int iDegree, iMin;
	double dSec;
	iDegree = int(dDms); //截取整数部分，即度
	iMin = int((dDms - iDegree) * 100); //截取分
	dSec = ((dDms - iDegree) * 100 - iMin) * 100;//获得秒
	dDms = iDegree + double(iMin) / 60 + dSec / 3600; //把分和秒转换成度，再相加，得到角度

	dDms = dDms * PI / 180;//化为弧度
	return dDms;//返回弧度
	
}

（8）
//度转度分秒
double CMFCExp4Dlg::Deg2Dms(double dDeg)
{
	//用于存放度、分、秒三个值的变量
	int iDegree, iMin;
	double dSec;

	dDeg = dDeg * 180 / PI;//弧度转化为角度
	
	double dTmp;//临时变量
	iDegree = int(dDeg);//截取度
	dTmp = (dDeg - iDegree) * 60;//把度的小数部分转换为分
	iMin = int(dTmp);//截取分
	dSec = (dTmp - iMin) * 60;//把分的小数部分转换为秒
	dDeg = iDegree + double(iMin) / 100 + dSec / 10000;//合并为度分秒形式
	return dDeg;//返回度分秒
}