(2024-02-13, 07:30 PM)Sparhawk76 Wrote: (2023-04-18, 02:30 AM)mgrouch Wrote: Take a look at this link https://bareboat-necessities.github.io/m...-math.html

It should give you idea how to calculate it. There are code parts in python and links to

rust code as well.

Thanks for that info, I've tried reading through that a few times but unfortunately the math is beyond my level of education. Gives me enough info for some insight, but no idea how to practically apply it.

From a bit more research I have done, I have confirmed my guess that it is supposed to be output in meters and the sentence format (enough info that if I had heave data, I could use node-red to output a working sentence):

https://www.manualslib.com/manual/163824...l?page=520

@Sparhawk76

@seandepagnier

I’ve built a heave sensor now. It uses accelerometer and gyroscope and

runs on esp32. It integrates with SignalK and ever more with BBN OS.

https://github.com/bareboat-necessities/...riod-esp32
BBN Boat Heave Sensor using IMU

Boat Heave Sensor on m5stack atomS3.

m5stack atomS3 is ESP32 microcontroller with built-in IMU MPU6886 (accelerometer and gyroscope)

Use arduino IDE to compile and upload sketch to esp32

Check required libraries and version in .github/workflows/build.yaml

Produces NMEA-0183 XDR sentences over USB (See NmeaXDR.h). Baud rate: 115200.

Mount atomS3 with LCD facing up.

Estimating Boat Heave using IMU

Estimate vessel heave (vertical displacement) in ocean waves using IMU on esp32

The method for estimating wave height and heave from a moving boat implemented here is using the following on-line algorithm:

Sample MPU6886 3D acceleration and 3D gyroscope (angular velocities) measurements at about 250 Hz.

Estimate attitude and get attitude quaternion using Mahony agorithm. Using acceleration and gyroscope is enough. No magnetometer required because we only interested vertical acceleration for the next steps.

Double integrate vertical acceleration into vertical displacement using specially designed Kalman filter which corrects for integral drift in wave and corrects for the constant accelerometer bias.

Estimate observed heave frequency with Aranovskiy on-line filter. The correction for accelerometer bias is important for this step.

Smooth frequency produced by Aranovskiy filter with Kalman smoother.

Use another specially designed Kalman filter knowing the frequency and fusing model with trochoidal wave model to double integrate vertical acceleration. Assuming convergence of frequency this method would give real-time phase correction of heave compared to the first Kalman method. Doppler effect due to a boat moving relating to waves has no impact on displacement amplitude.

Implementation Notes

Rolling min/max algorithm with window of about couple wave periods samples to produce wave height measurement. Algorithm:

https://github.com/lemire/runningmaxmin from Daniel Lemire paper, and improvements from:

https://github.com/EvanBalster/STL_mono_wedge
References

Alexey A. Bobtsov, Nikolay A. Nikolaev, Olga V. Slita, Alexander S. Borgul, Stanislav V. Aranovskiy: The New Algorithm of Sinusoidal Signal Frequency Estimation. 11th IFAC International Workshop on Adaptation and Learning in Control and Signal Processing, 2013

Sharkh S. M., Hendijanizadeh2 M., Moshrefi-Torbati3 M., Abusara M. A.: A Novel Kalman Filter Based Technique for Calculating the Time History of Vertical Displacement of a Boat from Measured Acceleration, Marine Engineering Frontiers Volume 2, 2014

Daniel Lemire, Streaming Maximum-Minimum Filter Using No More than Three Comparisons per Element. Nordic Journal of Computing, 13 (4), pages 328-339, 2006.

Trochoidal Wave Wikipedia