/Kalman

Some Python Implementations of the Kalman Filter, from balzer82/Kalman

Primary LanguageJupyter Notebook

Multidimensional Kalman-Filter

See Vimeo for some Explanations. Or if you want to start with the basics, you might want to take a look at these Blogposts:

Some Python Implementations of the Kalman Filter

Kalman Filter Step

Kalman Filter with Constant Velocity Model

Situation covered: You drive with your car in a tunnel and the GPS signal is lost. Now the car has to determine, where it is in the tunnel. The only information it has, is the velocity in driving direction. The x and y component of the velocity (x˙ and y˙) can be calculated from the absolute velocity (revolutions of the wheels) and the heading of the vehicle (yaw rate sensor).

View IPython Notebook ~ See Vimeo

Kalman Filter with Constant Acceleration Model

in 2D

Situation covered: You have an acceleration sensor (in 2D: $\ddot x¨ and y¨) and try to calculate velocity (x˙ and y˙) as well as position (x and y) of a person holding a smartphone in his/her hand.

View IPython Notebook ~ See Vimeo

Second example is the same dynamic model but this time you measure the position as well as the acceleration. Both values have to be fused together with the Kalman Filter. Situation covered: You have an acceleration sensor (in 2D: x¨ and y¨) and a Position Sensor (e.g. GPS) and try to calculate velocity (x˙ and y˙) as well as position (x and y) of a person holding a smartphone in his/her hand.

View IPython Notebook ~ See Vimeo

in 3D

Third example is in 3D space, so the state vector is 9D. This model is for ball tracking or something else in 3D space.

Kalman 3D

View IPython Notebook

Adaptive Kalman Filter with Constant Velocity Model

Here the Measurement Covariance Matrix R is calculated dynamically via the maximum likelihood of the acutal standard deviation of the last measurements.

View IPython Notebook

Kalman Filter for Motorbike Lean Angle Estimation

Also know as the Gimbal Stabilization problem: You can measure the rotationrate, but need some validation for the correct lean angle from time to time, because simply an integration of the rotationrate adds up a lot of noise. There comes the vertical acceleration, which is a pretty good estimator for the angle in static situations. This Kalman Filter implementation fuses both together with some adaptive components.

View IPython Notebook

Extended Kalman Filter

EKF Filter Step

Extended Kalman Filter with Constant Turn Rate and Velocity (CTRV) Model

Situation covered: You have an velocity sensor which measures the vehicle speed (v) in heading direction (ψ) and a yaw rate sensor (ψ˙) which both have to fused with the position (x & y) from a GPS sensor.

State Vector

View IPython Notebook ~ See Vimeo

Extended Kalman Filter with Constant Heading and Constant Velocity (CHCV) Model

Situation covered: You have the position (x & y) from a GPS sensor and extimating the heading direction (ψ) and the velocity (v).

View IPython Notebook

Extended Kalman Filter with Constant Turn Rate and Acceleration (CTRA) Model

Situation covered: You have an acceleration and velocity sensor which measures the vehicle longitudinal acceleration and speed (v) in heading direction (ψ) and a yaw rate sensor (ψ˙) which all have to fused with the position (x & y) from a GPS sensor.

View IPython Notebook

License

CC-BY-SA2.0 Lizenz

You are free to:

  • Share — copy and redistribute the material in any medium or format for any purpose, even commercially.
  • Adapt — remix, transform, and build upon the material for any purpose, even commercially. The licensor cannot revoke these freedoms as long as you follow the license terms.

Under the following terms:

  • Attribution - You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
  • ShareAlike - If you remix, transform, or build upon the material, you must distribute your contributions under the same license as the original.

No additional restrictions - You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits.

Repo Stars History

Star History Chart