yodaUtility/sgp4/cTle.h

//
// cTle.h 
//
// This class will accept a single set of two-line elements and then allow
// a client to request specific fields, such as epoch, mean motion, 
// etc., from the set.
//
// Copyright 1996-2003 Michael F. Henry
//
#pragma once

#include "globals.h"
#include "stdafx.h"
/////////////////////////////////////////////////////////////////////////////
class cTle
{
public:
   cTle(string&, string&, string&);
   cTle(const cTle &tle);
   ~cTle();
   
   enum eTleLine
   {
      LINE_ZERO,
      LINE_ONE,
      LINE_TWO
   };

   enum eField
   {
      FLD_FIRST,
      FLD_NORADNUM = FLD_FIRST,
      FLD_INTLDESC,
      FLD_SET,       // TLE set number
      FLD_EPOCHYEAR, // Epoch: Last two digits of year
      FLD_EPOCHDAY,  // Epoch: Fractional Julian Day of year
      FLD_ORBITNUM,  // Orbit at epoch
      FLD_I,         // Inclination
      FLD_RAAN,      // R.A. ascending node
      FLD_E,         // Eccentricity
      FLD_ARGPER,    // Argument of perigee
      FLD_M,         // Mean anomaly
      FLD_MMOTION,   // Mean motion
      FLD_MMOTIONDT, // First time derivative of mean motion
      FLD_MMOTIONDT2,// Second time derivative of mean motion
      FLD_BSTAR,     // BSTAR Drag
      FLD_LAST       // MUST be last
   };

   enum eUnits
   {
      U_FIRST,
      U_RAD = U_FIRST,  // radians
      U_DEG,            // degrees
      U_NATIVE,         // TLE format native units (no conversion)
      U_LAST            // MUST be last
   };
   
   void Initialize();
   
   static int    CheckSum(const string&);
   static bool   IsValidLine(string&, eTleLine);
   static string ExpToDecimal(const string&);

   static void TrimLeft(string&);
   static void TrimRight(string&);
   
   double getField(eField fld,               // which field to retrieve
                   eUnits unit  = U_NATIVE,  // return units in rad, deg etc.
                   string *pstr = NULL,      // return ptr for str value
                   bool bStrUnits = false)   // 'true': append units to str val
                   const;
   string getName()  const { return m_strName; }
   string getLine1() const { return m_strLine1;}
   string getLine2() const { return m_strLine2;}

protected:
   static double ConvertUnits(double val, eField fld, eUnits units); 

private:
   string getUnits(eField) const;
   double getFieldNumeric(eField) const;

   // Satellite name and two data lines
   string m_strName;
   string m_strLine1;
   string m_strLine2;

   // Converted fields, in atof()-readable form
   string m_Field[FLD_LAST];

   // Cache of field values in "double" format
   typedef int FldKey;
   FldKey Key(eUnits u, eField f) const { return (u * 100) + f; }
   mutable map<FldKey, double>  m_mapCache;
};

///////////////////////////////////////////////////////////////////////////
//
// TLE data format
//
// [Reference: T.S. Kelso]
//
// Two line element data consists of three lines in the following format:
//
//  AAAAAAAAAAAAAAAAAAAAAA
//  1 NNNNNU NNNNNAAA NNNNN.NNNNNNNN +.NNNNNNNN +NNNNN-N +NNNNN-N N NNNNN
//  2 NNNNN NNN.NNNN NNN.NNNN NNNNNNN NNN.NNNN NNN.NNNN NN.NNNNNNNNNNNNNN
//  
//  Line 0 is a twenty-two-character name.
// 
//   Lines 1 and 2 are the standard Two-Line Orbital Element Set Format identical
//   to that used by NORAD and NASA.  The format description is:
//      
//     Line 1
//     Column    Description
//     01-01     Line Number of Element Data
//     03-07     Satellite Number
//     10-11     International Designator (Last two digits of launch year)
//     12-14     International Designator (Launch number of the year)
//     15-17     International Designator (Piece of launch)
//     19-20     Epoch Year (Last two digits of year)
//     21-32     Epoch (Julian Day and fractional portion of the day)
//     34-43     First Time Derivative of the Mean Motion
//               or Ballistic Coefficient (Depending on ephemeris type)
//     45-52     Second Time Derivative of Mean Motion (decimal point assumed;
//               blank if N/A)
//     54-61     BSTAR drag term if GP4 general perturbation theory was used.
//               Otherwise, radiation pressure coefficient.  (Decimal point assumed)
//     63-63     Ephemeris type
//     65-68     Element number
//     69-69     Check Sum (Modulo 10)
//               (Letters, blanks, periods, plus signs = 0; minus signs = 1)
//
//     Line 2
//     Column    Description
//     01-01     Line Number of Element Data
//     03-07     Satellite Number
//     09-16     Inclination [Degrees]
//     18-25     Right Ascension of the Ascending Node [Degrees]
//     27-33     Eccentricity (decimal point assumed)
//     35-42     Argument of Perigee [Degrees]
//     44-51     Mean Anomaly [Degrees]
//     53-63     Mean Motion [Revs per day]
//     64-68     Revolution number at epoch [Revs]
//     69-69     Check Sum (Modulo 10)
//        
//     All other columns are blank or fixed.
//          
// Example:
//      
// NOAA 6
// 1 11416U          86 50.28438588 0.00000140           67960-4 0  5293
// 2 11416  98.5105  69.3305 0012788  63.2828 296.9658 14.24899292346978

1	//
2	// cTle.h
3	//
4	// This class will accept a single set of two-line elements and then allow
5	// a client to request specific fields, such as epoch, mean motion,
6	// etc., from the set.
7	//
8	// Copyright 1996-2003 Michael F. Henry
9	//
10	#pragma once
11
12	#include "globals.h"
13	#include "stdafx.h"
14	/////////////////////////////////////////////////////////////////////////////
15	class cTle
16	{
17	public:
18	cTle(string&, string&, string&);
19	cTle(const cTle &tle);
20	~cTle();
21
22	enum eTleLine
23	{
24	LINE_ZERO,
25	LINE_ONE,
26	LINE_TWO
27	};
28
29	enum eField
30	{
31	FLD_FIRST,
32	FLD_NORADNUM = FLD_FIRST,
33	FLD_INTLDESC,
34	FLD_SET, // TLE set number
35	FLD_EPOCHYEAR, // Epoch: Last two digits of year
36	FLD_EPOCHDAY, // Epoch: Fractional Julian Day of year
37	FLD_ORBITNUM, // Orbit at epoch
38	FLD_I, // Inclination
39	FLD_RAAN, // R.A. ascending node
40	FLD_E, // Eccentricity
41	FLD_ARGPER, // Argument of perigee
42	FLD_M, // Mean anomaly
43	FLD_MMOTION, // Mean motion
44	FLD_MMOTIONDT, // First time derivative of mean motion
45	FLD_MMOTIONDT2,// Second time derivative of mean motion
46	FLD_BSTAR, // BSTAR Drag
47	FLD_LAST // MUST be last
48	};
49
50	enum eUnits
51	{
52	U_FIRST,
53	U_RAD = U_FIRST, // radians
54	U_DEG, // degrees
55	U_NATIVE, // TLE format native units (no conversion)
56	U_LAST // MUST be last
57	};
58
59	void Initialize();
60
61	static int CheckSum(const string&);
62	static bool IsValidLine(string&, eTleLine);
63	static string ExpToDecimal(const string&);
64
65	static void TrimLeft(string&);
66	static void TrimRight(string&);
67
68	double getField(eField fld, // which field to retrieve
69	eUnits unit = U_NATIVE, // return units in rad, deg etc.
70	string *pstr = NULL, // return ptr for str value
71	bool bStrUnits = false) // 'true': append units to str val
72	const;
73	string getName() const { return m_strName; }
74	string getLine1() const { return m_strLine1;}
75	string getLine2() const { return m_strLine2;}
76
77	protected:
78	static double ConvertUnits(double val, eField fld, eUnits units);
79
80	private:
81	string getUnits(eField) const;
82	double getFieldNumeric(eField) const;
83
84	// Satellite name and two data lines
85	string m_strName;
86	string m_strLine1;
87	string m_strLine2;
88
89	// Converted fields, in atof()-readable form
90	string m_Field[FLD_LAST];
91
92	// Cache of field values in "double" format
93	typedef int FldKey;
94	FldKey Key(eUnits u, eField f) const { return (u * 100) + f; }
95	mutable map<FldKey, double> m_mapCache;
96	};
97
98	///////////////////////////////////////////////////////////////////////////
99	//
100	// TLE data format
101	//
102	// [Reference: T.S. Kelso]
103	//
104	// Two line element data consists of three lines in the following format:
105	//
106	// AAAAAAAAAAAAAAAAAAAAAA
107	// 1 NNNNNU NNNNNAAA NNNNN.NNNNNNNN +.NNNNNNNN +NNNNN-N +NNNNN-N N NNNNN
108	// 2 NNNNN NNN.NNNN NNN.NNNN NNNNNNN NNN.NNNN NNN.NNNN NN.NNNNNNNNNNNNNN
109	//
110	// Line 0 is a twenty-two-character name.
111	//
112	// Lines 1 and 2 are the standard Two-Line Orbital Element Set Format identical
113	// to that used by NORAD and NASA. The format description is:
114	//
115	// Line 1
116	// Column Description
117	// 01-01 Line Number of Element Data
118	// 03-07 Satellite Number
119	// 10-11 International Designator (Last two digits of launch year)
120	// 12-14 International Designator (Launch number of the year)
121	// 15-17 International Designator (Piece of launch)
122	// 19-20 Epoch Year (Last two digits of year)
123	// 21-32 Epoch (Julian Day and fractional portion of the day)
124	// 34-43 First Time Derivative of the Mean Motion
125	// or Ballistic Coefficient (Depending on ephemeris type)
126	// 45-52 Second Time Derivative of Mean Motion (decimal point assumed;
127	// blank if N/A)
128	// 54-61 BSTAR drag term if GP4 general perturbation theory was used.
129	// Otherwise, radiation pressure coefficient. (Decimal point assumed)
130	// 63-63 Ephemeris type
131	// 65-68 Element number
132	// 69-69 Check Sum (Modulo 10)
133	// (Letters, blanks, periods, plus signs = 0; minus signs = 1)
134	//
135	// Line 2
136	// Column Description
137	// 01-01 Line Number of Element Data
138	// 03-07 Satellite Number
139	// 09-16 Inclination [Degrees]
140	// 18-25 Right Ascension of the Ascending Node [Degrees]
141	// 27-33 Eccentricity (decimal point assumed)
142	// 35-42 Argument of Perigee [Degrees]
143	// 44-51 Mean Anomaly [Degrees]
144	// 53-63 Mean Motion [Revs per day]
145	// 64-68 Revolution number at epoch [Revs]
146	// 69-69 Check Sum (Modulo 10)
147	//
148	// All other columns are blank or fixed.
149	//
150	// Example:
151	//
152	// NOAA 6
153	// 1 11416U 86 50.28438588 0.00000140 67960-4 0 5293
154	// 2 11416 98.5105 69.3305 0012788 63.2828 296.9658 14.24899292346978
155