yodaUtility/sgp4/SxP4Test.cpp

//
// SxP4Test.cpp 
// This sample code demonstrates how to use the C++ classes in order
// to determine satellite position and look angles.
//
// mfh 12/07/2003
//
#include "stdafx.h"

#include <stdio.h>
#include "cTle.h"
#include "cEci.h"
#include "cOrbit.h"
#include "cSite.h"


/////////////////////////////////////////////////////////////////////////////
// Test routine to output position and velocity information
void PrintPosVel(const cTle &tle)
{
   cOrbit       orbit(tle);
   cEci         eci;
   vector<cEci> Pos;

   // Calculate position, velocity
   // mpe = "minutes past epoch"
   for (int mpe = 0; mpe <= (360 * 4); mpe += 360)
   {
      // Get the position of the satellite at time "mpe"
      // The coordinates are placed into the local variable "eci".
      orbit.getPosition(mpe, &eci);
    
      // Push the coordinates object onto the end of the vector.
      Pos.push_back(eci);
   }

   // Print TLE data
   printf("%s\n",   tle.getName().c_str());
   printf("%s\n",   tle.getLine1().c_str());
   printf("%s\n\n", tle.getLine2().c_str());

   // Header
   printf("  TSINCE            X                Y                Z\n\n");

   // Iterate over each of the ECI position objects pushed onto the
   // position vector, above, printing the ECI position information
   // as we go.
   for (unsigned int i = 0; i < Pos.size(); i++)
   {
      printf("%8d.00  %16.8f %16.8f %16.8f\n",
               i * 360,
               Pos[i].getPos().m_x,
               Pos[i].getPos().m_y,
               Pos[i].getPos().m_z);
   }

   printf("\n                    XDOT             YDOT             ZDOT\n\n");

   // Iterate over each of the ECI position objects in the position
   // vector again, but this time print the velocity information.
   for (unsigned int i = 0; i < Pos.size(); i++)
   {
      printf("             %16.8f %16.8f %16.8f\n",
             Pos[i].getVel().m_x,
             Pos[i].getVel().m_y,
             Pos[i].getVel().m_z);
   }
}

//////////////////////////////////////////////////////////////////////////////
int main(int argc, char* argv[])
{
   // Test SGP4
   string str1 = "SGP4 Test";
   string str2 = "1 88888U          80275.98708465  .00073094  13844-3  66816-4 0    8";
   string str3 = "2 88888  72.8435 115.9689 0086731  52.6988 110.5714 16.05824518  105";

   cTle tle1(str1, str2, str3);

   PrintPosVel(tle1);

   printf("\n");

   // Test SDP4
   str1 = "SDP4 Test";
   str2 = "1 11801U          80230.29629788  .01431103  00000-0  14311-1       8";
   str3 = "2 11801  46.7916 230.4354 7318036  47.4722  10.4117  2.28537848     6";

   cTle tleSDP4(str1, str2, str3);

   PrintPosVel(tleSDP4);

   printf("\nExample output:\n");

   // Example: Define a location on the earth, then determine the look-angle
   // to the SDP4 satellite defined above.

   // Create an orbit object using the SDP4 TLE object.
   cOrbit orbitSDP4(tleSDP4);

   // Create an ECI object to hold the location of the satellite
   cEci eciSDP4;

   // Get the location of the satellite from the Orbit object. The 
   // earth-centered inertial information is placed into eciSDP4.
   // Here we ask for the location of the satellite 90 minutes after
   // the TLE epoch.
   orbitSDP4.getPosition(90.0, &eciSDP4);

   // Now create a site object. Site objects represent a location on the 
   // surface of the earth. Here we arbitrarily select a point on the
   // equator.
   cSite siteEquator(0.0, -100.0, 0); // 0.00 N, 100.00 W, 0 km altitude

   // Now get the "look angle" from the site to the satellite. 
   // Note that the ECI object "eciSDP4" contains a time associated
   // with the coordinates it contains; this is the time at which
   // the look angle is valid.
   cCoordTopo topoLook = siteEquator.getLookAngle(eciSDP4);

   // Print out the results. Note that the Azimuth and Elevation are
   // stored in the cCoordTopo object as radians. Here we convert
   // to degrees using rad2deg()
   printf("AZ: %.1f  EL: %.1f\n", 
          rad2deg(topoLook.m_Az), 
          rad2deg(topoLook.m_El));

}

1	//
2	// SxP4Test.cpp
3	// This sample code demonstrates how to use the C++ classes in order
4	// to determine satellite position and look angles.
5	//
6	// mfh 12/07/2003
7	//
8	#include "stdafx.h"
9
10	#include <stdio.h>
11	#include "cTle.h"
12	#include "cEci.h"
13	#include "cOrbit.h"
14	#include "cSite.h"
15
16
17	/////////////////////////////////////////////////////////////////////////////
18	// Test routine to output position and velocity information
19	void PrintPosVel(const cTle &tle)
20	{
21	cOrbit orbit(tle);
22	cEci eci;
23	vector<cEci> Pos;
24
25	// Calculate position, velocity
26	// mpe = "minutes past epoch"
27	for (int mpe = 0; mpe <= (360 * 4); mpe += 360)
28	{
29	// Get the position of the satellite at time "mpe"
30	// The coordinates are placed into the local variable "eci".
31	orbit.getPosition(mpe, &eci);
32
33	// Push the coordinates object onto the end of the vector.
34	Pos.push_back(eci);
35	}
36
37	// Print TLE data
38	printf("%s\n", tle.getName().c_str());
39	printf("%s\n", tle.getLine1().c_str());
40	printf("%s\n\n", tle.getLine2().c_str());
41
42	// Header
43	printf(" TSINCE X Y Z\n\n");
44
45	// Iterate over each of the ECI position objects pushed onto the
46	// position vector, above, printing the ECI position information
47	// as we go.
48	for (unsigned int i = 0; i < Pos.size(); i++)
49	{
50	printf("%8d.00 %16.8f %16.8f %16.8f\n",
51	i * 360,
52	Pos[i].getPos().m_x,
53	Pos[i].getPos().m_y,
54	Pos[i].getPos().m_z);
55	}
56
57	printf("\n XDOT YDOT ZDOT\n\n");
58
59	// Iterate over each of the ECI position objects in the position
60	// vector again, but this time print the velocity information.
61	for (unsigned int i = 0; i < Pos.size(); i++)
62	{
63	printf(" %16.8f %16.8f %16.8f\n",
64	Pos[i].getVel().m_x,
65	Pos[i].getVel().m_y,
66	Pos[i].getVel().m_z);
67	}
68	}
69
70	//////////////////////////////////////////////////////////////////////////////
71	int main(int argc, char* argv[])
72	{
73	// Test SGP4
74	string str1 = "SGP4 Test";
75	string str2 = "1 88888U 80275.98708465 .00073094 13844-3 66816-4 0 8";
76	string str3 = "2 88888 72.8435 115.9689 0086731 52.6988 110.5714 16.05824518 105";
77
78	cTle tle1(str1, str2, str3);
79
80	PrintPosVel(tle1);
81
82	printf("\n");
83
84	// Test SDP4
85	str1 = "SDP4 Test";
86	str2 = "1 11801U 80230.29629788 .01431103 00000-0 14311-1 8";
87	str3 = "2 11801 46.7916 230.4354 7318036 47.4722 10.4117 2.28537848 6";
88
89	cTle tleSDP4(str1, str2, str3);
90
91	PrintPosVel(tleSDP4);
92
93	printf("\nExample output:\n");
94
95	// Example: Define a location on the earth, then determine the look-angle
96	// to the SDP4 satellite defined above.
97
98	// Create an orbit object using the SDP4 TLE object.
99	cOrbit orbitSDP4(tleSDP4);
100
101	// Create an ECI object to hold the location of the satellite
102	cEci eciSDP4;
103
104	// Get the location of the satellite from the Orbit object. The
105	// earth-centered inertial information is placed into eciSDP4.
106	// Here we ask for the location of the satellite 90 minutes after
107	// the TLE epoch.
108	orbitSDP4.getPosition(90.0, &eciSDP4);
109
110	// Now create a site object. Site objects represent a location on the
111	// surface of the earth. Here we arbitrarily select a point on the
112	// equator.
113	cSite siteEquator(0.0, -100.0, 0); // 0.00 N, 100.00 W, 0 km altitude
114
115	// Now get the "look angle" from the site to the satellite.
116	// Note that the ECI object "eciSDP4" contains a time associated
117	// with the coordinates it contains; this is the time at which
118	// the look angle is valid.
119	cCoordTopo topoLook = siteEquator.getLookAngle(eciSDP4);
120
121	// Print out the results. Note that the Azimuth and Elevation are
122	// stored in the cCoordTopo object as radians. Here we convert
123	// to degrees using rad2deg()
124	printf("AZ: %.1f EL: %.1f\n",
125	rad2deg(topoLook.m_Az),
126	rad2deg(topoLook.m_El));
127
128	}
129