KJ135
/
sgp4
mirror of http://172.16.0.12:5000/LAMPSARToolSoftware/sgp4.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Formatting update

feature/19
Daniel Warner 2011-12-14 11:48:22 +00:00
parent 5f1583a6fc
commit 0fc2d8f57f
2 changed files with 16 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Eci.cpp
View File

@ -17,8 +17,8 @@ void Eci::ToEci(const Julian& date, const CoordGeodetic &g)
/* /*
* take into account earth flattening * take into account earth flattening
*/ */
const double c = 1.0 / const double c = 1.0
sqrt(1.0 + kF * (kF - 2.0) * pow(sin(g.latitude), 2.0)); / sqrt(1.0 + kF * (kF - 2.0) * pow(sin(g.latitude), 2.0));
const double s = pow(1.0 - kF, 2.0) * c; const double s = pow(1.0 - kF, 2.0) * c;
const double achcp = (kXKMPER * c + g.altitude) * cos(g.latitude); const double achcp = (kXKMPER * c + g.altitude) * cos(g.latitude);
@ -54,8 +54,8 @@ CoordGeodetic Eci::ToGeodetic() const
// 180 >= lon < 180 // 180 >= lon < 180
const double lon = fmod(theta - date_.ToGreenwichSiderealTime(), kPI); const double lon = fmod(theta - date_.ToGreenwichSiderealTime(), kPI);
const double r = sqrt((position_.x * position_.x) + const double r = sqrt((position_.x * position_.x)
(position_.y * position_.y)); + (position_.y * position_.y));
static const double e2 = kF * (2.0 - kF); static const double e2 = kF * (2.0 - kF);

23
Observer.cpp
View File

@ -1,11 +1,10 @@
#include "Observer.h" #include "Observer.h"
#include "Globals.h"
/* /*
* calculate lookangle between the observer and the passed in Eci object * calculate lookangle between the observer and the passed in Eci object
*/ */
CoordTopographic Observer::GetLookAngle(const Eci &eci) { CoordTopographic Observer::GetLookAngle(const Eci &eci)
{
/* /*
* update the observers Eci to match the time of the Eci passed in * update the observers Eci to match the time of the Eci passed in
@ -31,21 +30,23 @@ CoordTopographic Observer::GetLookAngle(const Eci &eci) {
double sin_theta = sin(theta); double sin_theta = sin(theta);
double cos_theta = cos(theta); double cos_theta = cos(theta);
double top_s = sin_lat * cos_theta * range.x + double top_s = sin_lat * cos_theta * range.x
sin_lat * sin_theta * range.y - + sin_lat * sin_theta * range.y - cos_lat * range.z;
cos_lat * range.z; double top_e = -sin_theta * range.x
double top_e = -sin_theta * range.x + + cos_theta * range.y;
cos_theta * range.y; double top_z = cos_lat * cos_theta * range.x
double top_z = cos_lat * cos_theta * range.x + + cos_lat * sin_theta * range.y + sin_lat * range.z;
cos_lat * sin_theta * range.y +
sin_lat * range.z;
double az = atan(-top_e / top_s); double az = atan(-top_e / top_s);
if (top_s > 0.0) if (top_s > 0.0)
{
az += kPI; az += kPI;
}
if (az < 0.0) if (az < 0.0)
{
az += 2.0 * kPI; az += 2.0 * kPI;
}
double el = asin(top_z / range.w); double el = asin(top_z / range.w);
double rate = range.Dot(range_rate) / range.w; double rate = range.Dot(range_rate) / range.w;