#ifndef _GETPOSITION_H
#define _GETPOSITION_H

#include "Time.h"
#include "Coordinate.h"
#include "Rinex.h"
#include <list>

using namespace std;
const double _DAY_IN_SECOND=86400;
const double _HOUR_IN_SECOND=3600;
const double _MINUTE_IN_SECOND=60;
const double c=2.99792458e+08;//光速
const double GM=3.986005e+14;//重力常數(shù)
const double we=7.2921151467e-05;//地球自轉(zhuǎn)角速度
const double F=-4.442807633e-10;
const double f1=1575.42e+6;



typedef struct tagUtilParam 
{
	double n;//平均角速度
	double E;//偏近點(diǎn)角
	double u0;
	double vk;//真近點(diǎn)角
	double uk;//經(jīng)過(guò)改正的升交角距
	double r;//經(jīng)過(guò)改正的向徑
	double i;//經(jīng)過(guò)改正的軌道傾角
	double L;//改正后的升交點(diǎn)經(jīng)度	
} UtilParam;
typedef UtilParam *PUtilParam;



typedef struct tagSatVelocity //衛(wèi)星速度
{         
	double xv;
	double yv;
	double zv;
} SatVel;
typedef SatVel	*PSatVel;

//得到最靠近所給時(shí)刻的歷元單位
list<PGMNREC>::iterator GetBestGMNREC(list<GMNREC>& navRecord,
				   int& nPRN,PCOMMONTIME pctEpoch);

//得到有關(guān)計(jì)算衛(wèi)星空間位置，發(fā)射時(shí)刻改正，衛(wèi)星速度有關(guān)的參數(shù)
//避免代碼重復(fù)
void GetUtilParameter(list<PGMNREC>& navRecord,int& nPRN,PCOMMONTIME pctEpoch,
					  PUtilParam pParam);

//計(jì)算衛(wèi)星在笛卡爾坐標(biāo)系中的位置
void GetOrbNClk(list<PGMNREC>& navRecord,int& nPRN,PCOMMONTIME pctEpoch, 
				PCRDCARTESIAN pcrdOrb/*,double* pdSVClkBias*/);


//計(jì)算衛(wèi)星C/A碼信號(hào)發(fā)射時(shí)刻的改正
void GetSVClkBias(list<PGMNREC>& navRecord,int& nPRN,PCOMMONTIME pctEpoch,
				  double* pdSVClkBias,double *detj);

//計(jì)算衛(wèi)星的速度
void GetSatVelocity(list<PGMNREC>& navRecord,int& nPRN,PCOMMONTIME pctEpoch,
					PSatVel psatv);




typedef struct tagSITEV    //測(cè)站速度
{
	double xv;
	double yv;
	double zv;
}SITEV;
typedef SITEV *PSITEV;


typedef struct tagPPONERESULT{
	COMMONTIME        epoch;//歷元時(shí)刻
	int               sat_num;//參與解算的衛(wèi)星數(shù)
	CRDCARTESIAN      crd;  /////////////             /////////////                  ///////////     測(cè)站坐標(biāo)
	double            clk_bias;//接收機(jī)鐘差
	double            PDOP;

	SITEV             m_sitev;//測(cè)站速度
	int               flag;//是否能進(jìn)行速度解算，0代表不能，1代表能
}PPONERESULT;

typedef PPONERESULT *PPPONERESULT;

  

void PPOne(GMOHDR& gmoh,/*觀測(cè)值頭文件*/
		   list<PGMOREC>::iterator pgmor,/*某一歷元觀測(cè)值 數(shù)據(jù)記錄*/

		   list<PGMNREC>&  navRecord,/*導(dǎo)航電文文件*/
		   PPPONERESULT presult);

#endif