2 * Copyright © 2017 Keith Packard <keithp@keithp.com>
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, either version 2 of the License, or
7 * (at your option) any later version.
9 * This program is distributed in the hope that it will be useful, but
10 * WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 * General Public License for more details.
15 package org.altusmetrum.altoslib_13;
19 public class AltosFlightSeries extends AltosDataListener {
21 public ArrayList<AltosTimeSeries> series = new ArrayList<AltosTimeSeries>();
23 public double speed_filter_width = 4.0;
24 public double accel_filter_width = 1.0;
26 public int[] indices() {
27 int[] indices = new int[series.size()];
28 for (int i = 0; i < indices.length; i++)
30 step_indices(indices);
34 private double time(int id, int index) {
35 AltosTimeSeries s = series.get(id);
38 return Double.NEGATIVE_INFINITY;
40 if (index < s.values.size())
41 return s.values.get(index).time;
42 return Double.POSITIVE_INFINITY;
45 public boolean step_indices(int[] indices) {
46 double min_next = time(0, indices[0]+1);
48 for (int i = 1; i < indices.length; i++) {
49 double next = time(i, indices[i]+1);
54 if (min_next == Double.POSITIVE_INFINITY)
57 for (int i = 0; i < indices.length; i++) {
58 double t = time(i, indices[i] + 1);
66 public double time(int[] indices) {
67 double max = time(0, indices[0]);
69 for (int i = 1; i < indices.length; i++) {
70 double t = time(i, indices[i]);
77 public double value(String name, int[] indices) {
78 for (int i = 0; i < indices.length; i++) {
79 AltosTimeSeries s = series.get(i);
80 if (s.label.equals(name)) {
81 int index = indices[i];
84 if (index >= s.values.size())
85 index = s.values.size() - 1;
86 return s.values.get(index).value;
89 return AltosLib.MISSING;
92 public double value(String name, double time) {
93 for (AltosTimeSeries s : series) {
94 if (s.label.equals(name))
97 return AltosLib.MISSING;
100 public double value_before(String name, double time) {
101 for (AltosTimeSeries s : series) {
102 if (s.label.equals(name))
103 return s.value_before(time);
105 return AltosLib.MISSING;
108 public double value_after(String name, double time) {
109 for (AltosTimeSeries s : series) {
110 if (s.label.equals(name))
111 return s.value_after(time);
113 return AltosLib.MISSING;
116 public AltosTimeSeries make_series(String label, AltosUnits units) {
117 return new AltosTimeSeries(label, units);
120 public void add_series(AltosTimeSeries s) {
121 for (int e = 0; e < series.size(); e++) {
122 if (s.compareTo(series.get(e)) < 0){
130 public AltosTimeSeries add_series(String label, AltosUnits units) {
131 AltosTimeSeries s = make_series(label, units);
136 public void remove_series(AltosTimeSeries s) {
140 public boolean has_series(String label) {
141 for (AltosTimeSeries s : series)
142 if (s.label.equals(label))
147 public AltosTimeSeries state_series;
149 public static final String state_name = "State";
151 public void set_state(int state) {
153 if (state != AltosLib.ao_flight_pad && state != AltosLib.MISSING && state != AltosLib.ao_flight_stateless) {
154 if (state_series == null)
155 state_series = add_series(state_name, AltosConvert.state_name);
156 if (this.state() != state)
157 state_series.add(time(), state);
159 super.set_state(state);
162 public AltosTimeSeries accel_series;
163 public boolean accel_computed;
165 public static final String accel_name = "Accel";
167 public AltosTimeSeries vert_accel_series;
169 public static final String vert_accel_name = "Vertical Accel";
171 public void set_acceleration(double acceleration) {
172 if (acceleration == AltosLib.MISSING)
174 if (accel_series == null)
175 accel_series = add_series(accel_name, AltosConvert.accel);
177 accel_series.add(time(), acceleration);
178 accel_computed = false;
181 private AltosTimeSeries compute_accel() {
182 AltosTimeSeries new_accel_series = null;
184 if (speed_series != null) {
185 AltosTimeSeries temp_series;
186 if (accel_filter_width > 0) {
187 temp_series = make_series(speed_name, AltosConvert.speed);
188 speed_series.filter(temp_series, accel_filter_width);
190 temp_series = speed_series;
192 new_accel_series = make_series(accel_name, AltosConvert.accel);
193 temp_series.differentiate(new_accel_series);
195 return new_accel_series;
198 public void set_filter(double speed_filter, double accel_filter) {
199 this.speed_filter_width = speed_filter;
200 this.accel_filter_width = accel_filter;
202 AltosTimeSeries new_speed_series = compute_speed();
204 if (new_speed_series != null) {
205 speed_series.erase_values();
206 for (AltosTimeValue tv : new_speed_series)
207 speed_series.add(tv);
209 if (accel_computed) {
210 AltosTimeSeries new_accel_series = compute_accel();
211 if (new_accel_series != null) {
212 accel_series.erase_values();
213 for (AltosTimeValue tv : new_accel_series)
214 accel_series.add(tv);
219 public void set_received_time(long received_time) {
222 public AltosTimeSeries tick_series;
224 public static final String tick_name = "Tick";
226 public void set_tick(int tick) {
227 super.set_tick(tick);
228 if (tick_series == null)
229 tick_series = add_series(tick_name, null);
230 tick_series.add(time(), tick);
233 public AltosTimeSeries rssi_series;
235 public static final String rssi_name = "RSSI";
237 public AltosTimeSeries status_series;
239 public static final String status_name = "Radio Status";
241 public void set_rssi(int rssi, int status) {
242 if (rssi_series == null) {
243 rssi_series = add_series(rssi_name, null);
244 status_series = add_series(status_name, null);
246 rssi_series.add(time(), rssi);
247 status_series.add(time(), status);
250 public AltosTimeSeries pressure_series;
252 public static final String pressure_name = "Pressure";
254 public AltosTimeSeries altitude_series;
256 public static final String altitude_name = "Altitude";
258 public AltosTimeSeries height_series;
260 public double max_height = AltosLib.MISSING;
262 public void set_min_pressure(double pa) {
263 double ground_altitude = cal_data().ground_altitude;
264 if (ground_altitude != AltosLib.MISSING)
265 max_height = AltosConvert.pressure_to_altitude(pa) -
269 public static final String height_name = "Height";
271 public void set_pressure(double pa) {
272 if (pa == AltosLib.MISSING)
275 if (pressure_series == null)
276 pressure_series = add_series(pressure_name, AltosConvert.pressure);
277 pressure_series.add(time(), pa);
278 if (altitude_series == null)
279 altitude_series = add_series(altitude_name, AltosConvert.height);
281 if (cal_data().ground_pressure == AltosLib.MISSING)
282 cal_data().set_ground_pressure(pa);
284 double altitude = AltosConvert.pressure_to_altitude(pa);
285 altitude_series.add(time(), altitude);
288 private void compute_height() {
289 double ground_altitude = cal_data().ground_altitude;
290 if (height_series == null && ground_altitude != AltosLib.MISSING && altitude_series != null) {
291 height_series = add_series(height_name, AltosConvert.height);
292 for (AltosTimeValue alt : altitude_series)
293 height_series.add(alt.time, alt.value - ground_altitude);
296 if (gps_height == null && cal_data().gps_pad != null && cal_data().gps_pad.alt != AltosLib.MISSING && gps_altitude != null) {
297 double gps_ground_altitude = cal_data().gps_pad.alt;
298 gps_height = add_series(gps_height_name, AltosConvert.height);
299 for (AltosTimeValue gps_alt : gps_altitude)
300 gps_height.add(gps_alt.time, gps_alt.value - gps_ground_altitude);
304 public AltosTimeSeries speed_series;
306 public static final String speed_name = "Speed";
308 private AltosTimeSeries compute_speed() {
309 AltosTimeSeries new_speed_series = null;
310 AltosTimeSeries alt_speed_series = null;
311 AltosTimeSeries accel_speed_series = null;
313 if (altitude_series != null) {
314 AltosTimeSeries temp_series;
316 if (speed_filter_width > 0) {
317 temp_series = make_series(speed_name, AltosConvert.height);
318 altitude_series.filter(temp_series, speed_filter_width);
320 temp_series = altitude_series;
322 alt_speed_series = make_series(speed_name, AltosConvert.speed);
323 temp_series.differentiate(alt_speed_series);
325 if (accel_series != null && !accel_computed) {
327 if (orient_series != null) {
328 vert_accel_series = add_series(vert_accel_name, AltosConvert.accel);
330 for (AltosTimeValue a : accel_series) {
331 double orient = orient_series.value(a.time);
332 double a_abs = a.value + AltosConvert.gravity;
333 double v_a = a_abs * Math.cos(AltosConvert.degrees_to_radians(orient)) - AltosConvert.gravity;
335 vert_accel_series.add(a.time, v_a);
339 AltosTimeSeries temp_series = make_series(speed_name, AltosConvert.speed);
341 if (vert_accel_series != null)
342 vert_accel_series.integrate(temp_series);
344 accel_series.integrate(temp_series);
346 accel_speed_series = make_series(speed_name, AltosConvert.speed);
347 temp_series.filter(accel_speed_series, 0.1);
350 if (alt_speed_series != null && accel_speed_series != null) {
351 double apogee_time = AltosLib.MISSING;
352 if (state_series != null) {
353 for (AltosTimeValue d : state_series) {
354 if (d.value >= AltosLib.ao_flight_drogue){
355 apogee_time = d.time;
360 if (apogee_time == AltosLib.MISSING) {
361 new_speed_series = alt_speed_series;
363 new_speed_series = make_series(speed_name, AltosConvert.speed);
364 for (AltosTimeValue d : accel_speed_series) {
365 if (d.time <= apogee_time)
366 new_speed_series.add(d);
368 for (AltosTimeValue d : alt_speed_series) {
369 if (d.time > apogee_time)
370 new_speed_series.add(d);
374 } else if (alt_speed_series != null) {
375 new_speed_series = alt_speed_series;
376 } else if (accel_speed_series != null) {
377 new_speed_series = accel_speed_series;
379 return new_speed_series;
382 public AltosTimeSeries orient_series;
384 public static final String orient_name = "Tilt Angle";
386 private void compute_orient() {
388 if (orient_series != null)
391 if (accel_ground_across == AltosLib.MISSING)
394 AltosCalData cal_data = cal_data();
396 if (cal_data.pad_orientation == AltosLib.MISSING)
399 if (cal_data.accel_zero_across == AltosLib.MISSING)
402 AltosRotation rotation = new AltosRotation(accel_ground_across,
403 accel_ground_through,
405 cal_data.pad_orientation);
406 double prev_time = ground_time;
408 orient_series = add_series(orient_name, AltosConvert.orient);
409 orient_series.add(ground_time, rotation.tilt());
411 for (AltosTimeValue roll_v : gyro_roll) {
412 double time = roll_v.time;
413 double dt = time - prev_time;
416 double roll = AltosConvert.degrees_to_radians(roll_v.value) * dt;
417 double pitch = AltosConvert.degrees_to_radians(gyro_pitch.value(time)) * dt;
418 double yaw = AltosConvert.degrees_to_radians(gyro_yaw.value(time)) * dt;
420 rotation.rotate(pitch, yaw, roll);
421 orient_series.add(time, rotation.tilt());
427 public AltosTimeSeries kalman_height_series, kalman_speed_series, kalman_accel_series;
429 public static final String kalman_height_name = "Kalman Height";
430 public static final String kalman_speed_name = "Kalman Speed";
431 public static final String kalman_accel_name = "Kalman Accel";
433 public void set_kalman(double height, double speed, double acceleration) {
434 if (kalman_height_series == null) {
435 kalman_height_series = add_series(kalman_height_name, AltosConvert.height);
436 kalman_speed_series = add_series(kalman_speed_name, AltosConvert.speed);
437 kalman_accel_series = add_series(kalman_accel_name, AltosConvert.accel);
439 kalman_height_series.add(time(), height);
440 kalman_speed_series.add(time(), speed);
441 kalman_accel_series.add(time(), acceleration);
444 public AltosTimeSeries thrust_series;
446 public static final String thrust_name = "Thrust";
448 public void set_thrust(double N) {
449 if (thrust_series == null)
450 thrust_series = add_series(thrust_name, AltosConvert.force);
451 thrust_series.add(time(), N);
454 public AltosTimeSeries temperature_series;
456 public static final String temperature_name = "Temperature";
458 public void set_temperature(double deg_c) {
459 if (temperature_series == null)
460 temperature_series = add_series(temperature_name, AltosConvert.temperature);
461 temperature_series.add(time(), deg_c);
464 public AltosTimeSeries battery_voltage_series;
466 public static final String battery_voltage_name = "Battery Voltage";
468 public void set_battery_voltage(double volts) {
469 if (volts == AltosLib.MISSING)
471 if (battery_voltage_series == null)
472 battery_voltage_series = add_series(battery_voltage_name, AltosConvert.voltage);
473 battery_voltage_series.add(time(), volts);
476 public AltosTimeSeries apogee_voltage_series;
478 public static final String apogee_voltage_name = "Apogee Voltage";
480 public void set_apogee_voltage(double volts) {
481 if (volts == AltosLib.MISSING)
483 if (apogee_voltage_series == null)
484 apogee_voltage_series = add_series(apogee_voltage_name, AltosConvert.voltage);
485 apogee_voltage_series.add(time(), volts);
488 public AltosTimeSeries main_voltage_series;
490 public static final String main_voltage_name = "Main Voltage";
492 public void set_main_voltage(double volts) {
493 if (volts == AltosLib.MISSING)
495 if (main_voltage_series == null)
496 main_voltage_series = add_series(main_voltage_name, AltosConvert.voltage);
497 main_voltage_series.add(time(), volts);
500 public ArrayList<AltosGPSTimeValue> gps_series;
502 public AltosGPS gps_before(double time) {
503 AltosGPSTimeValue nearest = null;
504 for (AltosGPSTimeValue gtv : gps_series) {
508 if (gtv.time <= time) {
509 if (nearest.time <= time && gtv.time > nearest.time)
512 if (nearest.time > time && gtv.time < nearest.time)
523 public AltosTimeSeries sats_in_view;
524 public AltosTimeSeries sats_in_soln;
525 public AltosTimeSeries gps_altitude;
526 public AltosTimeSeries gps_height;
527 public AltosTimeSeries gps_ground_speed;
528 public AltosTimeSeries gps_ascent_rate;
529 public AltosTimeSeries gps_course;
530 public AltosTimeSeries gps_speed;
531 public AltosTimeSeries gps_pdop, gps_vdop, gps_hdop;
533 public static final String sats_in_view_name = "Satellites in view";
534 public static final String sats_in_soln_name = "Satellites in solution";
535 public static final String gps_altitude_name = "GPS Altitude";
536 public static final String gps_height_name = "GPS Height";
537 public static final String gps_ground_speed_name = "GPS Ground Speed";
538 public static final String gps_ascent_rate_name = "GPS Ascent Rate";
539 public static final String gps_course_name = "GPS Course";
540 public static final String gps_speed_name = "GPS Speed";
541 public static final String gps_pdop_name = "GPS Dilution of Precision";
542 public static final String gps_vdop_name = "GPS Vertical Dilution of Precision";
543 public static final String gps_hdop_name = "GPS Horizontal Dilution of Precision";
545 public void set_gps(AltosGPS gps, boolean new_location, boolean new_sats) {
546 super.set_gps(gps, new_location, new_sats);
547 AltosCalData cal_data = cal_data();
548 if (gps_series == null)
549 gps_series = new ArrayList<AltosGPSTimeValue>();
550 gps_series.add(new AltosGPSTimeValue(time(), gps));
553 if (sats_in_soln == null) {
554 sats_in_soln = add_series(sats_in_soln_name, null);
556 sats_in_soln.add(time(), gps.nsat);
557 if (gps.pdop != AltosLib.MISSING) {
558 if (gps_pdop == null)
559 gps_pdop = add_series(gps_pdop_name, null);
560 gps_pdop.add(time(), gps.pdop);
562 if (gps.hdop != AltosLib.MISSING) {
563 if (gps_hdop == null)
564 gps_hdop = add_series(gps_hdop_name, null);
565 gps_hdop.add(time(), gps.hdop);
567 if (gps.vdop != AltosLib.MISSING) {
568 if (gps_vdop == null)
569 gps_vdop = add_series(gps_vdop_name, null);
570 gps_vdop.add(time(), gps.vdop);
573 if (gps.alt != AltosLib.MISSING) {
574 if (gps_altitude == null)
575 gps_altitude = add_series(gps_altitude_name, AltosConvert.height);
576 gps_altitude.add(time(), gps.alt);
578 if (gps.ground_speed != AltosLib.MISSING) {
579 if (gps_ground_speed == null)
580 gps_ground_speed = add_series(gps_ground_speed_name, AltosConvert.speed);
581 gps_ground_speed.add(time(), gps.ground_speed);
583 if (gps.climb_rate != AltosLib.MISSING) {
584 if (gps_ascent_rate == null)
585 gps_ascent_rate = add_series(gps_ascent_rate_name, AltosConvert.speed);
586 gps_ascent_rate.add(time(), gps.climb_rate);
588 if (gps.course != AltosLib.MISSING) {
589 if (gps_course == null)
590 gps_course = add_series(gps_course_name, null);
591 gps_course.add(time(), gps.course);
593 if (gps.ground_speed != AltosLib.MISSING && gps.climb_rate != AltosLib.MISSING) {
594 if (gps_speed == null)
595 gps_speed = add_series(gps_speed_name, null);
596 gps_speed.add(time(), Math.sqrt(gps.ground_speed * gps.ground_speed +
597 gps.climb_rate * gps.climb_rate));
602 if (gps.cc_gps_sat != null) {
603 if (sats_in_view == null)
604 sats_in_view = add_series(sats_in_view_name, null);
605 sats_in_view.add(time(), gps.cc_gps_sat.length);
610 public static final String accel_along_name = "Accel Along";
611 public static final String accel_across_name = "Accel Across";
612 public static final String accel_through_name = "Accel Through";
614 public AltosTimeSeries accel_along, accel_across, accel_through;
616 public static final String gyro_roll_name = "Roll Rate";
617 public static final String gyro_pitch_name = "Pitch Rate";
618 public static final String gyro_yaw_name = "Yaw Rate";
620 public AltosTimeSeries gyro_roll, gyro_pitch, gyro_yaw;
622 public static final String mag_along_name = "Magnetic Field Along";
623 public static final String mag_across_name = "Magnetic Field Across";
624 public static final String mag_through_name = "Magnetic Field Through";
626 public AltosTimeSeries mag_along, mag_across, mag_through;
628 public void set_accel(double along, double across, double through) {
629 if (accel_along == null) {
630 accel_along = add_series(accel_along_name, AltosConvert.accel);
631 accel_across = add_series(accel_across_name, AltosConvert.accel);
632 accel_through = add_series(accel_through_name, AltosConvert.accel);
634 accel_along.add(time(), along);
635 accel_across.add(time(), across);
636 accel_through.add(time(), through);
639 private double accel_ground_along = AltosLib.MISSING;
640 private double accel_ground_across = AltosLib.MISSING;
641 private double accel_ground_through = AltosLib.MISSING;
643 private double ground_time;
645 public void set_accel_ground(double along, double across, double through) {
646 accel_ground_along = along;
647 accel_ground_across = across;
648 accel_ground_through = through;
649 ground_time = time();
652 public void set_gyro(double roll, double pitch, double yaw) {
653 if (gyro_roll == null) {
654 gyro_roll = add_series(gyro_roll_name, AltosConvert.rotation_rate);
655 gyro_pitch = add_series(gyro_pitch_name, AltosConvert.rotation_rate);
656 gyro_yaw = add_series(gyro_yaw_name, AltosConvert.rotation_rate);
658 gyro_roll.add(time(), roll);
659 gyro_pitch.add(time(), pitch);
660 gyro_yaw.add(time(), yaw);
663 public void set_mag(double along, double across, double through) {
664 if (mag_along == null) {
665 mag_along = add_series(mag_along_name, AltosConvert.magnetic_field);
666 mag_across = add_series(mag_across_name, AltosConvert.magnetic_field);
667 mag_through = add_series(mag_through_name, AltosConvert.magnetic_field);
669 mag_along.add(time(), along);
670 mag_across.add(time(), across);
671 mag_through.add(time(), through);
674 public void set_orient(double orient) {
675 if (orient_series == null)
676 orient_series = add_series(orient_name, AltosConvert.orient);
677 orient_series.add(time(), orient);
680 public static final String pyro_voltage_name = "Pyro Voltage";
682 public AltosTimeSeries pyro_voltage;
684 public void set_pyro_voltage(double volts) {
685 if (pyro_voltage == null)
686 pyro_voltage = add_series(pyro_voltage_name, AltosConvert.voltage);
687 pyro_voltage.add(time(), volts);
690 private static String[] igniter_voltage_names;
692 public String igniter_voltage_name(int channel) {
693 if (igniter_voltage_names == null || igniter_voltage_names.length <= channel) {
694 String[] new_igniter_voltage_names = new String[channel + 1];
697 if (igniter_voltage_names != null) {
698 for (; i < igniter_voltage_names.length; i++)
699 new_igniter_voltage_names[i] = igniter_voltage_names[i];
701 for (; i < channel+1; i++)
702 new_igniter_voltage_names[i] = AltosLib.igniter_name(i);
703 igniter_voltage_names = new_igniter_voltage_names;
705 return igniter_voltage_names[channel];
708 public AltosTimeSeries[] igniter_voltage;
710 public void set_igniter_voltage(double[] voltage) {
711 int channels = voltage.length;
712 if (igniter_voltage == null || igniter_voltage.length <= channels) {
713 AltosTimeSeries[] new_igniter_voltage = new AltosTimeSeries[channels];
716 if (igniter_voltage != null) {
717 for (; i < igniter_voltage.length; i++)
718 new_igniter_voltage[i] = igniter_voltage[i];
720 for (; i < channels; i++)
721 new_igniter_voltage[i] = add_series(igniter_voltage_name(i), AltosConvert.voltage);
722 igniter_voltage = new_igniter_voltage;
724 for (int channel = 0; channel < voltage.length; channel++)
725 igniter_voltage[channel].add(time(), voltage[channel]);
728 public static final String pyro_fired_name = "Pyro Channel State";
730 public AltosTimeSeries pyro_fired_series;
734 public void set_pyro_fired(int pyro_mask) {
735 if (pyro_fired_series == null)
736 pyro_fired_series = add_series(pyro_fired_name, AltosConvert.pyro_name);
737 for (int channel = 0; channel < 32; channel++) {
738 if ((last_pyro_mask & (1 << channel)) == 0 &&
739 (pyro_mask & (1 << channel)) != 0) {
740 pyro_fired_series.add(time(), channel);
743 last_pyro_mask = pyro_mask;
746 public void set_companion(AltosCompanion companion) {
749 public void finish() {
751 if (speed_series == null) {
752 speed_series = compute_speed();
753 if (speed_series != null)
754 add_series(speed_series);
756 if (accel_series == null) {
757 accel_series = compute_accel();
758 if (accel_series != null) {
759 add_series(accel_series);
760 accel_computed = true;
766 public AltosTimeSeries[] series() {
768 return series.toArray(new AltosTimeSeries[0]);
771 public AltosFlightSeries(AltosCalData cal_data) {