Merge branch 'master' of ssh://git.gag.com/scm/git/fw/altos
[fw/altos] / altoslib / AltosFlightStats.java
1 /*
2  * Copyright © 2011 Keith Packard <keithp@keithp.com>
3  *
4  * This program is free software; you can redistribute it and/or modify
5  * it under the terms of the GNU General Public License as published by
6  * the Free Software Foundation; version 2 of the License.
7  *
8  * This program is distributed in the hope that it will be useful, but
9  * WITHOUT ANY WARRANTY; without even the implied warranty of
10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
11  * General Public License for more details.
12  *
13  * You should have received a copy of the GNU General Public License along
14  * with this program; if not, write to the Free Software Foundation, Inc.,
15  * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
16  */
17
18 package org.altusmetrum.altoslib_7;
19
20 import java.io.*;
21
22 public class AltosFlightStats {
23         public double           max_height;
24         public double           max_gps_height;
25         public double           max_speed;
26         public double           max_acceleration;
27         public double[] state_speed = new double[AltosLib.ao_flight_invalid + 1];
28         public double[] state_accel = new double[AltosLib.ao_flight_invalid + 1];
29         public int[]            state_count = new int[AltosLib.ao_flight_invalid + 1];
30         public double[] state_start = new double[AltosLib.ao_flight_invalid + 1];
31         public double[] state_end = new double[AltosLib.ao_flight_invalid + 1];
32         public int              serial;
33         public int              flight;
34         public int              year, month, day;
35         public int              hour, minute, second;
36         public double           lat, lon;
37         public double           pad_lat, pad_lon;
38         public boolean          has_flight_data;
39         public boolean          has_gps;
40         public boolean          has_gps_sats;
41         public boolean          has_gps_detail;
42         public boolean          has_flight_adc;
43         public boolean          has_battery;
44         public boolean          has_rssi;
45         public boolean          has_imu;
46         public boolean          has_mag;
47         public boolean          has_orient;
48         public int              num_ignitor;
49
50         double landed_time(AltosStateIterable states) {
51                 AltosState state = null;
52
53                 for (AltosState s : states) {
54                         state = s;
55                         if (state.state == AltosLib.ao_flight_landed)
56                                 break;
57                 }
58
59                 if (state == null)
60                         return 0;
61
62                 double  landed_height = state.height();
63
64                 state = null;
65
66                 boolean above = true;
67
68                 double  landed_time = -1000;
69
70                 for (AltosState s : states) {
71                         state = s;
72
73                         if (state.height() > landed_height + 10) {
74                                 above = true;
75                         } else {
76                                 if (above && state.height() < landed_height + 2) {
77                                         above = false;
78                                         landed_time = state.time;
79                                 }
80                         }
81                 }
82                 if (landed_time == -1000)
83                         landed_time = state.time;
84                 return landed_time;
85         }
86
87         double boost_time(AltosStateIterable states) {
88                 double boost_time = AltosLib.MISSING;
89                 AltosState      state = null;
90
91                 for (AltosState s : states) {
92                         state = s;
93                         if (state.acceleration() < 1)
94                                 boost_time = state.time;
95                         if (state.state >= AltosLib.ao_flight_boost && state.state <= AltosLib.ao_flight_landed)
96                                 break;
97                 }
98                 if (state == null)
99                         return 0;
100
101                 if (boost_time == AltosLib.MISSING)
102                         boost_time = state.time;
103                 return boost_time;
104         }
105
106
107         public AltosFlightStats(AltosStateIterable states) throws InterruptedException, IOException {
108                 double          boost_time = boost_time(states);
109                 double          end_time = 0;
110                 double          landed_time = landed_time(states);
111
112                 year = month = day = AltosLib.MISSING;
113                 hour = minute = second = AltosLib.MISSING;
114                 serial = flight = AltosLib.MISSING;
115                 lat = lon = AltosLib.MISSING;
116                 has_flight_data = false;
117                 has_gps = false;
118                 has_gps_sats = false;
119                 has_flight_adc = false;
120                 has_battery = false;
121                 has_rssi = false;
122                 has_imu = false;
123                 has_mag = false;
124                 has_orient = false;
125                 for (AltosState state : states) {
126                         if (serial == AltosLib.MISSING && state.serial != AltosLib.MISSING)
127                                 serial = state.serial;
128                         if (flight == AltosLib.MISSING && state.flight != AltosLib.MISSING)
129                                 flight = state.flight;
130                         if (state.battery_voltage != AltosLib.MISSING)
131                                 has_battery = true;
132                         if (state.main_voltage != AltosLib.MISSING)
133                                 has_flight_adc = true;
134                         if (state.rssi != AltosLib.MISSING)
135                                 has_rssi = true;
136                         end_time = state.time;
137
138                         if (state.pressure() != AltosLib.MISSING)
139                                 has_flight_data = true;
140
141                         int state_id = state.state;
142                         if (state.time >= boost_time && state_id < AltosLib.ao_flight_boost)
143                                 state_id = AltosLib.ao_flight_boost;
144                         if (state.time >= landed_time && state_id < AltosLib.ao_flight_landed)
145                                 state_id = AltosLib.ao_flight_landed;
146                         if (state.gps != null && state.gps.locked) {
147                                 year = state.gps.year;
148                                 month = state.gps.month;
149                                 day = state.gps.day;
150                                 hour = state.gps.hour;
151                                 minute = state.gps.minute;
152                                 second = state.gps.second;
153                         }
154                         max_height = state.max_height();
155                         max_speed = state.max_speed();
156                         max_acceleration = state.max_acceleration();
157                         max_gps_height = state.max_gps_height();
158
159                         if (0 <= state_id && state_id < AltosLib.ao_flight_invalid) {
160                                 double acceleration = state.acceleration();
161                                 double speed = state.speed();
162                                 if (acceleration != AltosLib.MISSING && speed != AltosLib.MISSING) {
163                                         state_accel[state_id] += acceleration;
164                                         state_speed[state_id] += speed;
165                                         state_count[state_id]++;
166                                 }
167                                 if (state_start[state_id] == 0.0)
168                                         state_start[state_id] = state.time;
169                                 if (state_end[state_id] < state.time)
170                                         state_end[state_id] = state.time;
171                         }
172                         if (state.pad_lat != AltosLib.MISSING) {
173                                 pad_lat = state.pad_lat;
174                                 pad_lon = state.pad_lon;
175                         }
176                         if (state.gps != null && state.gps.locked && state.gps.nsat >= 4) {
177                                 lat = state.gps.lat;
178                                 lon = state.gps.lon;
179                                 has_gps = true;
180                                 if (state.gps.cc_gps_sat != null)
181                                         has_gps_sats = true;
182                                 if (state.gps.course != AltosLib.MISSING)
183                                         has_gps_detail = true;
184                         }
185                         if (state.imu != null)
186                                 has_imu = true;
187                         if (state.mag != null)
188                                 has_mag = true;
189                         if (state.orient() != AltosLib.MISSING)
190                                 has_orient = true;
191                         if (state.ignitor_voltage != null && state.ignitor_voltage.length > num_ignitor)
192                                 num_ignitor = state.ignitor_voltage.length;
193                 }
194                 for (int s = AltosLib.ao_flight_startup; s <= AltosLib.ao_flight_landed; s++) {
195                         if (state_count[s] > 0) {
196                                 state_speed[s] /= state_count[s];
197                                 state_accel[s] /= state_count[s];
198                         } else {
199                                 state_speed[s] = AltosLib.MISSING;
200                                 state_accel[s] = AltosLib.MISSING;
201                         }
202                         if (state_start[s] == 0)
203                                 state_start[s] = end_time;
204                         if (state_end[s] == 0)
205                                 state_end[s] = end_time;
206                 }
207         }
208 }