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	kusanagi	1.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