#include?“StdAfx.h“

#include?
#include?“GaussProjection.h“


#define?Pi?3.141592653589793

GaussProjection::GaussProjection（）
{
	
}

GaussProjection::~GaussProjection（）
{
	
}

/*************************************************
*??Function:???????AngleToRadian（）
*??Description:???????????
*??parameter:??????Angle（DD.MMSSS）
*??Return:?????????Radian
**************************************************/
double?GaussProjection::AngleToRadian（double?dfm1）
{
	//角度化弧度參數
	double?p_d?=?Pi?/?180;
	double?p_f?=?p_d?/?60;
????double?p_m?=?p_f?/?60;
	
	//彌補存儲缺陷
	double?dfm?=?dfm1;
	
	double?du?=?floor（dfm）;
	double?temp?=?dfm?-?du;
	double?fe?=?floor（temp?*?100）;
	double?mi?=?（temp?*?100?-?fe）?*?100;
	
	double?radian?=?du?*?p_d?+?fe?*?p_f?+?mi?*?p_m;
	
????return?radian;
}


/*************************************************
*??Function:???????RadianToAngle（）
*??Description:???????????
*??parameter:??????Radian
*??Return:?????????Angle（DD.MMSSS）
**************************************************/
double?GaussProjection::RadianToAngle（double?radian）
{
	double?flag;
	if?（radian?>?0）????{?flag?=?1;}
	else
	{
		flag?=?-1;
		radian?=?fabs（radian）;
	}
	
	double?angle_tmp?=?radian?*?180?/?Pi;	
	double?Du?=?floor（angle_tmp）;
	double?Fe?=?floor（（angle_tmp?-?Du）?*?60）;
	double?Mi?=?（（angle_tmp?-?Du）?*?60?-?Fe）?*?60;
	
	if?（fabs（Mi-60）<0.0001）??????//秒60進1???分60進1
	{
		Mi?=?0;
		Fe?=?Fe?+?1;
		if?（fabs（Fe?-?60）?		{
			Fe?=?0;
			Du?=?Du?+?1;
		}
	}???????????	
	double?angle?=?（Du?+?Fe?/?100?+?Mi?/?10000）?*?flag;
????return?angle;
}



void?GaussProjection::Nekrasov_ForWard（double?Bdouble?Ldouble?Lodouble?&xdouble?&y）
{
	//輔助量
	double?cosB?=?cos（B）;
	double?sinB?=?sin（B）;
	double?cosB_2?=?cosB?*?cosB;
	double?l?=?L?-?Lo?*?Pi?/?180;
	double?ll?=?l?*?l;
	
	//計算系數
	double?N?=?6399698.902?-?（21562.267?-?（108.973?-?0.612?*?cosB_2）?*?cosB_2）?*?cosB_2;
	double?a0?=?32140.404?-?（135.3302?-?（0.7092?-?0.0040?*?cosB_2）?*?cosB_2）?*?cosB_2;
	double?a4?=?（0.25?+?0.00252?*?cosB_2）?*?cosB_2?-?0.04166;
	double?a6?=?（0.166?*?cosB_2?-?0.084）?*?cosB_2;
	double?a3?=?（0.3333333?+?0.001123?*?cosB_2）?*?cosB_2?-?0.1666667;
	double?a5?=?0.0083?-?（0.1667?-?（0.1968?+?0.0040?*?cosB_2）?*?cosB_2）?*?cosB_2;
	
	//計算高斯平面坐標值
	x?=?6367558.4969?*?B?-?（a0?-?（0.5?+?（a4?+?a6?*?ll）?*?ll）?*?ll?*?N）?*?cosB?*?sinB;
	y?=?（1?+?（a3?+?a5?*?ll）?*?ll）?*?l?*?N?*?cosB?+?500000;
}


void?GaussProjection::Inte1975_Forward（double?Bdouble?Ldouble?Lodouble?&xdouble?&y）
{
	//輔助量
	double?cosB?=?cos（B）;
	double?sinB?=?sin（B）;
	double?cosB_2?=?cosB?*?cosB;
	double?l?=?L?-?Lo?*?Pi?/?180;
	double?ll?=?l?*?l;
	
	//計算系數
	double?N?=?6399596.596?-?（21565.045?-?（108.996?-?0.603?*?cosB_2）?*?cosB_2）?*?cosB_2;
	double?a0?=?32144.5189?-?（135.3646?-?（0.7034?-?0.0041?*?cosB_2）?*?cosB_2）?*?cosB_2;
	double?a4?=?（0.25?+?0.00253?*?cosB_2）?*?cosB_2?-?0.04167;
	double?a6?=?（0.167?*?cosB_2?-?0.083）?*?cosB_2;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-05-06?22:07??Projection\
?????目錄???????????0??2012-04-19?22:10??Projection\Debug\
?????文件???????21689??2012-04-19?20:19??Projection\Debug\GaussProjection.obj
?????文件???????20031??2012-04-19?20:12??Projection\Debug\MainFrm.obj
?????文件??????122965??2012-04-20?20:55??Projection\Debug\Projection.exe
?????文件??????347340??2012-04-20?20:55??Projection\Debug\Projection.ilk
?????文件???????23273??2012-04-19?22:10??Projection\Debug\Projection.obj
?????文件?????5499652??2012-04-19?20:12??Projection\Debug\Projection.pch
?????文件??????484352??2012-04-19?22:10??Projection\Debug\Projection.pdb
?????文件????????8036??2012-04-19?20:12??Projection\Debug\Projection.res
?????文件???????15030??2012-04-19?20:12??Projection\Debug\ProjectionDoc.obj
?????文件???????28692??2012-04-20?20:55??Projection\Debug\ProjectionView.obj
?????文件??????105734??2012-04-19?20:12??Projection\Debug\StdAfx.obj
?????文件??????214016??2012-04-20?20:55??Projection\Debug\vc60.idb
?????文件??????364544??2012-04-19?22:10??Projection\Debug\vc60.pdb
?????文件????????5279??2012-05-06?22:07??Projection\GaussProjection.cpp
?????文件?????????534??2012-04-19?20:19??Projection\GaussProjection.h
?????文件????????2511??2012-04-19?19:18??Projection\MainFrm.cpp
?????文件????????1581??2012-04-19?19:18??Projection\MainFrm.h
?????文件???????45400??2012-04-19?19:49??Projection\Projection.aps
?????文件????????3127??2012-05-06?22:07??Projection\Projection.clw
?????文件????????4281??2012-04-19?19:18??Projection\Projection.cpp
?????文件????????4781??2012-04-19?22:13??Projection\Projection.dsp
?????文件?????????545??2012-04-19?19:18??Projection\Projection.dsw
?????文件????????1400??2012-04-19?19:18??Projection\Projection.h
?????文件???????66560??2012-05-06?22:07??Projection\Projection.ncb
?????文件???????48640??2012-05-06?22:07??Projection\Projection.opt
?????文件?????????794??2012-04-20?20:55??Projection\Projection.plg
?????文件???????13205??2012-04-19?19:49??Projection\Projection.rc
?????文件????????1822??2012-04-19?19:18??Projection\ProjectionDoc.cpp
?????文件????????1519??2012-04-19?19:18??Projection\ProjectionDoc.h
............此處省略11個文件信息