2 * Copyright © 2010 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; version 2 of the License.
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.
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.
19 * Track flight state from telemetry data stream
24 import altosui.AltosTelemetry;
25 import altosui.AltosGPS;
27 public class AltosState {
38 boolean ascent; /* going up? */
40 double ground_altitude;
51 double max_acceleration;
62 AltosGreatCircle from_pad;
67 double speak_altitude;
72 return System.currentTimeMillis() / 1000.0;
75 void init (AltosTelemetry cur, AltosState prev_state) {
78 double accel_counts_per_mss;
82 ground_altitude = AltosConvert.cc_barometer_to_altitude(data.ground_pres);
83 height = AltosConvert.cc_barometer_to_altitude(data.flight_pres) - ground_altitude;
85 report_time = aoview_time();
87 accel_counts_per_mss = ((data.accel_minus_g - data.accel_plus_g) / 2.0) / 9.80665;
88 acceleration = (data.ground_accel - data.flight_accel) / accel_counts_per_mss;
89 speed = data.flight_vel / (accel_counts_per_mss * 100.0);
90 temperature = AltosConvert.cc_thermometer_to_temperature(data.temp);
91 drogue_sense = AltosConvert.cc_ignitor_to_voltage(data.drogue);
92 main_sense = AltosConvert.cc_ignitor_to_voltage(data.main);
93 battery = AltosConvert.cc_battery_to_voltage(data.batt);
97 if (prev_state != null) {
99 /* Preserve any existing gps data */
100 npad = prev_state.npad;
101 gps = prev_state.gps;
102 pad_lat = prev_state.pad_lat;
103 pad_lon = prev_state.pad_lon;
104 pad_alt = prev_state.pad_alt;
105 max_height = prev_state.max_height;
106 max_acceleration = prev_state.max_acceleration;
107 max_speed = prev_state.max_speed;
109 /* make sure the clock is monotonic */
110 while (tick < prev_state.tick)
113 time_change = (tick - prev_state.tick) / 100.0;
115 /* compute barometric speed */
117 double height_change = height - prev_state.height;
119 baro_speed = (prev_state.baro_speed * 3 + (height_change / time_change)) / 4.0;
121 baro_speed = prev_state.baro_speed;
129 if (state == AltosTelemetry.ao_flight_pad) {
130 if (data.gps != null && data.gps.gps_locked && data.gps.nsat >= 4) {
133 /* filter pad position */
134 pad_lat = (pad_lat * 31.0 + data.gps.lat) / 32.0;
135 pad_lon = (pad_lon * 31.0 + data.gps.lon) / 32.0;
136 pad_alt = (pad_alt * 31.0 + data.gps.alt) / 32.0;
138 pad_lat = data.gps.lat;
139 pad_lon = data.gps.lon;
140 pad_alt = data.gps.alt;
144 ascent = (AltosTelemetry.ao_flight_boost <= state &&
145 state <= AltosTelemetry.ao_flight_coast);
147 /* Only look at accelerometer data on the way up */
148 if (ascent && acceleration > max_acceleration)
149 max_acceleration = acceleration;
150 if (ascent && speed > max_speed)
153 if (height > max_height)
155 if (data.gps != null) {
156 if (gps == null || !gps.gps_locked || data.gps.gps_locked)
158 if (npad > 0 && gps.gps_locked)
159 from_pad = new AltosGreatCircle(pad_lat, pad_lon, gps.lat, gps.lon);
162 gps_height = gps.alt - pad_alt;
168 public AltosState(AltosTelemetry cur) {
172 public AltosState (AltosTelemetry cur, AltosState prev) {