James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 1 | // This script provides a basic utility for de-batching the IMUValues |
| 2 | // message. See imu_plotter.ts for usage. |
| 3 | import * as configuration from 'org_frc971/aos/configuration_generated'; |
| 4 | import * as imu from 'org_frc971/frc971/wpilib/imu_batch_generated'; |
| 5 | import {MessageHandler, TimestampedMessage} from 'org_frc971/aos/network/www/aos_plotter'; |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 6 | import {Point} from 'org_frc971/aos/network/www/plotter'; |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 7 | import {Table} from 'org_frc971/aos/network/www/reflection'; |
| 8 | import {ByteBuffer} from 'org_frc971/external/com_github_google_flatbuffers/ts/byte-buffer'; |
James Kuszmaul | f8355ff | 2021-12-19 22:08:45 -0800 | [diff] [blame] | 9 | import {Long} from 'org_frc971/external/com_github_google_flatbuffers/ts/long'; |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 10 | |
| 11 | import Schema = configuration.reflection.Schema; |
| 12 | import IMUValuesBatch = imu.frc971.IMUValuesBatch; |
| 13 | import IMUValues = imu.frc971.IMUValues; |
| 14 | |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 15 | const FILTER_WINDOW_SIZE = 100; |
| 16 | |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 17 | export class ImuMessageHandler extends MessageHandler { |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 18 | // Calculated magnitude of the measured acceleration from the IMU. |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 19 | private acceleration_magnitudes: Point[] = []; |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 20 | constructor(private readonly schema: Schema) { |
| 21 | super(schema); |
| 22 | } |
James Kuszmaul | f8355ff | 2021-12-19 22:08:45 -0800 | [diff] [blame] | 23 | private readScalar(table: Table, fieldName: string): number|Long|null { |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 24 | return this.parser.readScalar(table, fieldName); |
| 25 | } |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 26 | addMessage(data: Uint8Array, time: number): void { |
| 27 | const batch = IMUValuesBatch.getRootAsIMUValuesBatch( |
| 28 | new ByteBuffer(data) as unknown as flatbuffers.ByteBuffer); |
| 29 | for (let ii = 0; ii < batch.readingsLength(); ++ii) { |
| 30 | const message = batch.readings(ii); |
| 31 | const table = Table.getNamedTable( |
| 32 | message.bb as unknown as ByteBuffer, this.schema, 'frc971.IMUValues', |
| 33 | message.bb_pos); |
| 34 | if (this.parser.readScalar(table, "monotonic_timestamp_ns") == null) { |
| 35 | console.log('Ignoring unpopulated IMU values: '); |
| 36 | console.log(this.parser.toObject(table)); |
| 37 | continue; |
| 38 | } |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 39 | const time = message.monotonicTimestampNs().toFloat64() * 1e-9; |
| 40 | this.messages.push(new TimestampedMessage(table, time)); |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 41 | this.acceleration_magnitudes.push(new Point( |
| 42 | time, |
| 43 | Math.hypot( |
| 44 | message.accelerometerX(), message.accelerometerY(), |
| 45 | message.accelerometerZ()))); |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 46 | } |
| 47 | } |
| 48 | |
| 49 | // Computes a moving average for a given input, using a basic window centered |
| 50 | // on each value. |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 51 | private movingAverageCentered(input: Point[]): Point[] { |
| 52 | const num_measurements = input.length; |
| 53 | const filtered_measurements = []; |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 54 | for (let ii = 0; ii < num_measurements; ++ii) { |
| 55 | let sum = 0; |
| 56 | let count = 0; |
| 57 | for (let jj = Math.max(0, Math.ceil(ii - FILTER_WINDOW_SIZE / 2)); |
| 58 | jj < Math.min(num_measurements, ii + FILTER_WINDOW_SIZE / 2); ++jj) { |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 59 | sum += input[jj].y; |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 60 | ++count; |
| 61 | } |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 62 | filtered_measurements.push(new Point(input[ii].x, sum / count)); |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 63 | } |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 64 | return filtered_measurements; |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 65 | } |
| 66 | |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 67 | getField(field: string[]): Point[] { |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 68 | // Any requested input that ends with "_filtered" will get a moving average |
| 69 | // applied to the original field. |
| 70 | const filtered_suffix = "_filtered"; |
| 71 | if (field[0] == "acceleration_magnitude") { |
James Kuszmaul | 0d7df89 | 2021-04-09 22:19:49 -0700 | [diff] [blame] | 72 | return this.acceleration_magnitudes; |
James Kuszmaul | ac2b6b4 | 2021-03-07 22:38:06 -0800 | [diff] [blame] | 73 | } else if (field[0].endsWith(filtered_suffix)) { |
| 74 | return this.movingAverageCentered(this.getField( |
| 75 | [field[0].slice(0, field[0].length - filtered_suffix.length)])); |
| 76 | } else { |
| 77 | return super.getField(field); |
James Kuszmaul | 5f5e123 | 2020-12-22 20:58:00 -0800 | [diff] [blame] | 78 | } |
| 79 | } |
| 80 | } |