-
Notifications
You must be signed in to change notification settings - Fork 1
/
pass.cpp
103 lines (87 loc) · 2.54 KB
/
pass.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
/********************************************************************************
* File Name : pass.cpp
* Title :
* Programmer: Ron Hashimoto
* Belonging :
* Date : 2014.11.1
* Language : C++
*********************************************************************************
* class to calculate next AOS, LOS and TCA
*
********************************************************************************/
#include "pass.h"
PassFinder::PassFinder (void)
{
departure_unix_ = 0.0;
}
void PassFinder::setDepartureTime (double departure_unix)
{
departure_unix_ = departure_unix;
}
void PassFinder::getDepartureTime (double* departure_unix) const
{
*departure_unix = departure_unix_;
}
Pass PassFinder::findFirst(SpacecraftTracker* tracker, double begin_unix, double end_unix, double interval) const
{
static const Pass error_pass = {-1., -1., -1.};
Pass pass = {begin_unix, begin_unix, begin_unix};
if (!tracker || begin_unix < 0.0 || end_unix < 0.0 || begin_unix > end_unix) {
return error_pass;
}
if (tracker->setTargetTime (pass.AOS) != 0) {
return error_pass;
}
double elevation,azimuth;
tracker->getSpacecraftDirection (&elevation, &azimuth);
if (elevation > 0.0) {
while (elevation > 0.0) {
pass.AOS -= interval;
if (pass.AOS < departure_unix_) {
pass.AOS = departure_unix_;
break;
}
if (tracker->setTargetTime (pass.AOS) != 0) {
return error_pass;
}
tracker->getSpacecraftDirection (&elevation, &azimuth);
}
}
while (elevation < 0.0) {
pass.AOS += interval;
if (pass.AOS > end_unix) {
return error_pass;
}
if (tracker->setTargetTime (pass.AOS) != 0) {
return error_pass;
}
tracker->getSpacecraftDirection (&elevation, &azimuth);
}
pass.LOS = pass.AOS;
double max_elevation = 0.0;
while (elevation > 0.0) {
pass.LOS += interval;
if (pass.LOS > end_unix) {
return error_pass;
}
if (tracker->setTargetTime (pass.LOS) != 0) {
return error_pass;
}
tracker->getSpacecraftDirection (&elevation, &azimuth);
if (max_elevation < elevation) {
max_elevation = elevation;
pass.TCA = pass.LOS;
}
}
return pass;
}
std::vector<Pass> PassFinder::findAll(SpacecraftTracker* tracker, double begin_unix, double end_unix, double interval) const
{
std::vector<Pass> passes;
Pass pass = findFirst(tracker, begin_unix, end_unix, interval);
while (pass.AOS > 0) {
passes.push_back(pass);
pass = findFirst(tracker, pass.LOS, end_unix, interval);
}
return passes;
}