Paparazzi UAS  v6.0_unstable-41-g5cfda00-dirty
Paparazzi is a free software Unmanned Aircraft System.
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
ins_int.c
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2008-2010 The Paparazzi Team
3  *
4  * This file is part of paparazzi.
5  *
6  * paparazzi is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2, or (at your option)
9  * any later version.
10  *
11  * paparazzi is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with paparazzi; see the file COPYING. If not, write to
18  * the Free Software Foundation, 59 Temple Place - Suite 330,
19  * Boston, MA 02111-1307, USA.
20  */
21 
29 #include "subsystems/ins/ins_int.h"
30 
31 #include "subsystems/abi.h"
32 
33 #include "subsystems/imu.h"
34 #include "subsystems/gps.h"
35 
36 #include "generated/airframe.h"
37 
38 #if USE_VFF_EXTENDED
40 #else
42 #endif
43 
44 #if USE_HFF
46 #endif
47 
48 #if defined SITL && USE_NPS
49 //#include "nps_fdm.h"
50 #include "nps_autopilot.h"
51 #include <stdio.h>
52 #endif
53 
54 #include "math/pprz_geodetic_int.h"
55 #include "math/pprz_isa.h"
56 #include "math/pprz_stat.h"
57 
58 #ifndef VFF_R_AGL
59 #define VFF_R_AGL 0.2
60 #endif
61 
62 #if USE_SONAR
63 #if !USE_VFF_EXTENDED
64 #error USE_SONAR needs USE_VFF_EXTENDED
65 #endif
66 
67 #ifdef INS_SONAR_THROTTLE_THRESHOLD
69 #endif
70 
71 #ifndef INS_SONAR_MIN_RANGE
72 #define INS_SONAR_MIN_RANGE 0.001
73 #endif
74 #ifndef INS_SONAR_MAX_RANGE
75 #define INS_SONAR_MAX_RANGE 4.0
76 #endif
77 #define VFF_R_SONAR_0 0.2
78 #ifndef VFF_R_SONAR_OF_M
79 #define VFF_R_SONAR_OF_M 0.2
80 #endif
81 
82 #endif // USE_SONAR
83 
84 #if USE_GPS
85 #ifndef INS_VFF_R_GPS
86 #define INS_VFF_R_GPS 2.0
87 #endif
88 
89 #ifndef INS_VFF_VZ_R_GPS
90 #define INS_VFF_VZ_R_GPS 2.0
91 #endif
92 #endif // USE_GPS
93 
95 #ifndef INS_MAX_PROPAGATION_STEPS
96 #define INS_MAX_PROPAGATION_STEPS 200
97 #endif
98 
99 #ifndef USE_INS_NAV_INIT
100 #define USE_INS_NAV_INIT TRUE
101 PRINT_CONFIG_MSG("USE_INS_NAV_INIT defaulting to TRUE")
102 #endif
103 
105 #define INS_BARO_MAX_INIT_VAR 1.f // variance threshold to set initial baro measurement
106 #ifndef INS_INT_BARO_ID
107 #if USE_BARO_BOARD
108 #define INS_INT_BARO_ID BARO_BOARD_SENDER_ID
109 #else
110 #define INS_INT_BARO_ID ABI_BROADCAST
111 #endif
112 #endif
113 PRINT_CONFIG_VAR(INS_INT_BARO_ID)
115 static void baro_cb(uint8_t sender_id, uint32_t stamp, float pressure);
116 
120 #ifndef INS_INT_IMU_ID
121 #define INS_INT_IMU_ID ABI_BROADCAST
122 #endif
124 static void accel_cb(uint8_t sender_id, uint32_t stamp, struct Int32Vect3 *accel);
125 
126 #ifndef INS_INT_GPS_ID
127 #define INS_INT_GPS_ID GPS_MULTI_ID
128 #endif
130 static void gps_cb(uint8_t sender_id, uint32_t stamp, struct GpsState *gps_s);
131 
135 #ifndef INS_INT_VEL_ID
136 #define INS_INT_VEL_ID ABI_BROADCAST
137 #endif
139 static void vel_est_cb(uint8_t sender_id,
140  uint32_t stamp,
141  float x, float y, float z,
142  float noise_x, float noise_y, float noise_z);
143 #ifndef INS_INT_POS_ID
144 #define INS_INT_POS_ID ABI_BROADCAST
145 #endif
147 static void pos_est_cb(uint8_t sender_id,
148  uint32_t stamp,
149  float x, float y, float z,
150  float noise_x, float noise_y, float noise_z);
151 
155 #ifndef INS_INT_AGL_ID
156 #define INS_INT_AGL_ID ABI_BROADCAST
157 #endif
159 static void agl_cb(uint8_t sender_id, uint32_t stamp, float distance);
160 
162 
163 #if PERIODIC_TELEMETRY
165 
166 static void send_ins(struct transport_tx *trans, struct link_device *dev)
167 {
168  pprz_msg_send_INS(trans, dev, AC_ID,
172 }
173 
174 static void send_ins_z(struct transport_tx *trans, struct link_device *dev)
175 {
176  pprz_msg_send_INS_Z(trans, dev, AC_ID,
178 }
179 
180 static void send_ins_ref(struct transport_tx *trans, struct link_device *dev)
181 {
182  if (ins_int.ltp_initialized) {
183  pprz_msg_send_INS_REF(trans, dev, AC_ID,
187  }
188 }
189 #endif
190 
191 static void ins_ned_to_state(void);
192 static void ins_update_from_vff(void);
193 #if USE_HFF
194 static void ins_update_from_hff(void);
195 #endif
196 
197 
198 void ins_int_init(void)
199 {
200 
201 #if USE_INS_NAV_INIT
203  ins_int.ltp_initialized = true;
204 #else
205  ins_int.ltp_initialized = false;
206 #endif
207 
208  /* we haven't had any measurement updates yet, so set the counter to max */
210 
211  // Bind to BARO_ABS message
212  AbiBindMsgBARO_ABS(INS_INT_BARO_ID, &baro_ev, baro_cb);
213  ins_int.baro_initialized = false;
214 
215  ins_int.vf_reset = false;
216  ins_int.hf_realign = false;
217 
218  /* init vertical and horizontal filters */
219  vff_init_zero();
220 #if USE_HFF
221  hff_init(0., 0., 0., 0.);
222 #endif
223 
227 
228 #if PERIODIC_TELEMETRY
232 #endif
233 
234  /*
235  * Subscribe to scaled IMU measurements and attach callbacks
236  */
237  AbiBindMsgIMU_ACCEL_INT32(INS_INT_IMU_ID, &accel_ev, accel_cb);
238  AbiBindMsgGPS(INS_INT_GPS_ID, &gps_ev, gps_cb);
239  AbiBindMsgVELOCITY_ESTIMATE(INS_INT_VEL_ID, &vel_est_ev, vel_est_cb);
240  AbiBindMsgPOSITION_ESTIMATE(INS_INT_POS_ID, &pos_est_ev, pos_est_cb);
241  AbiBindMsgAGL(INS_INT_AGL_ID, &agl_ev, agl_cb); // ABI to the altitude above ground level
242 }
243 
245 {
246 #if USE_GPS
247  if (GpsFixValid()) {
251  ins_int.ltp_initialized = true;
253  } else {
254  ins_int.ltp_initialized = false;
255  }
256 #else
257  ins_int.ltp_initialized = false;
258 #endif
259 
260 #if USE_HFF
261  ins_int.hf_realign = true;
262 #endif
263  ins_int.vf_reset = true;
264 }
265 
267 {
268 #if USE_GPS
269  if (GpsFixValid()) {
270  struct LlaCoor_i lla = {
272  .lon = state.ned_origin_i.lla.lon,
273  .alt = gps.lla_pos.alt
274  };
278  }
279 #endif
280  ins_int.vf_reset = true;
281 }
282 
283 void ins_int_propagate(struct Int32Vect3 *accel, float dt)
284 {
285  /* untilt accels */
286  struct Int32Vect3 accel_meas_body;
287  struct Int32RMat *body_to_imu_rmat = orientationGetRMat_i(&imu.body_to_imu);
288  int32_rmat_transp_vmult(&accel_meas_body, body_to_imu_rmat, accel);
289  stateSetAccelBody_i(&accel_meas_body);
290  struct Int32Vect3 accel_meas_ltp;
291  int32_rmat_transp_vmult(&accel_meas_ltp, stateGetNedToBodyRMat_i(), &accel_meas_body);
292 
293  float z_accel_meas_float = ACCEL_FLOAT_OF_BFP(accel_meas_ltp.z);
294 
295  /* Propagate only if we got any measurement during the last INS_MAX_PROPAGATION_STEPS.
296  * Otherwise halt the propagation to not diverge and only set the acceleration.
297  * This should only be relevant in the startup phase when the baro is not yet initialized
298  * and there is no gps fix yet...
299  */
301  vff_propagate(z_accel_meas_float, dt);
303  } else {
304  // feed accel from the sensors
305  // subtract -9.81m/s2 (acceleration measured due to gravity,
306  // but vehicle not accelerating in ltp)
307  ins_int.ltp_accel.z = accel_meas_ltp.z + ACCEL_BFP_OF_REAL(9.81);
308  }
309 
310 #if USE_HFF
311  /* propagate horizontal filter */
312  hff_propagate();
313  /* convert and copy result to ins_int */
314  ins_update_from_hff();
315 #else
316  ins_int.ltp_accel.x = accel_meas_ltp.x;
317  ins_int.ltp_accel.y = accel_meas_ltp.y;
318 #endif /* USE_HFF */
319 
321 
322  /* increment the propagation counter, while making sure it doesn't overflow */
325  }
326 }
327 
328 static void baro_cb(uint8_t __attribute__((unused)) sender_id, __attribute__((unused)) uint32_t stamp, float pressure)
329 {
330  if (pressure < 1.f)
331  {
332  // bad baro pressure, don't use
333  return;
334  }
335 
336  if (!ins_int.baro_initialized) {
337 #define press_hist_len 10
338  static float press_hist[press_hist_len];
339  static uint8_t idx = 0;
340 
341  press_hist[idx] = pressure;
342  idx = (idx + 1) % press_hist_len;
343  float var = variance_f(press_hist, press_hist_len);
344  if (var < INS_BARO_MAX_INIT_VAR){
345  // wait for a first positive value
346  ins_int.vf_reset = true;
347  ins_int.baro_initialized = true;
348  }
349  }
350 
352  float height_correction = 0.f;
354  // Calculate the distance to the origin
355  struct EnuCoor_f *enu = stateGetPositionEnu_f();
356  double dist2_to_origin = enu->x * enu->x + enu->y * enu->y;
357 
358  // correction for the earth's curvature
359  const double earth_radius = 6378137.0;
360  height_correction = (float)(sqrt(earth_radius * earth_radius + dist2_to_origin) - earth_radius);
361  }
362 
363  if (ins_int.vf_reset) {
364  ins_int.vf_reset = false;
365  ins_int.qfe = pressure;
366  vff_realign(height_correction);
368  }
369 
370  float baro_up = pprz_isa_height_of_pressure(pressure, ins_int.qfe);
371 
372  // The VFF will update in the NED frame
373  ins_int.baro_z = -(baro_up - height_correction);
374 
375 #if USE_VFF_EXTENDED
377 #else
379 #endif
380 
381  /* reset the counter to indicate we just had a measurement update */
383  }
384 }
385 
386 #if USE_GPS
387 void ins_int_update_gps(struct GpsState *gps_s)
388 {
389  if (gps_s->fix < GPS_FIX_3D) {
390  return;
391  }
392 
393  if (!ins_int.ltp_initialized) {
395  }
396 
397  struct NedCoor_i gps_pos_cm_ned;
398  ned_of_ecef_point_i(&gps_pos_cm_ned, &ins_int.ltp_def, &gps_s->ecef_pos);
399 
400  /* calculate body frame position taking BODY_TO_GPS translation (in cm) into account */
401 #ifdef INS_BODY_TO_GPS_X
402  /* body2gps translation in body frame */
403  struct Int32Vect3 b2g_b = {
404  .x = INS_BODY_TO_GPS_X,
405  .y = INS_BODY_TO_GPS_Y,
406  .z = INS_BODY_TO_GPS_Z
407  };
408  /* rotate offset given in body frame to navigation/ltp frame using current attitude */
409  struct Int32Quat q_b2n = *stateGetNedToBodyQuat_i();
410  QUAT_INVERT(q_b2n, q_b2n);
411  struct Int32Vect3 b2g_n;
412  int32_quat_vmult(&b2g_n, &q_b2n, &b2g_b);
413  /* subtract body2gps translation in ltp from gps position */
414  VECT3_SUB(gps_pos_cm_ned, b2g_n);
415 #endif
416 
418  struct NedCoor_i gps_speed_cm_s_ned;
419  ned_of_ecef_vect_i(&gps_speed_cm_s_ned, &ins_int.ltp_def, &gps_s->ecef_vel);
420 
421 #if INS_USE_GPS_ALT
422  vff_update_z_conf(((float)gps_pos_cm_ned.z) / 100.0, INS_VFF_R_GPS);
423 #endif
424 #if INS_USE_GPS_ALT_SPEED
425  vff_update_vz_conf(((float)gps_speed_cm_s_ned.z) / 100.0, INS_VFF_VZ_R_GPS);
427 #endif
428 
429 #if USE_HFF
430  /* horizontal gps transformed to NED in meters as float */
431  struct FloatVect2 gps_pos_m_ned;
432  VECT2_ASSIGN(gps_pos_m_ned, gps_pos_cm_ned.x, gps_pos_cm_ned.y);
433  VECT2_SDIV(gps_pos_m_ned, gps_pos_m_ned, 100.0f);
434 
435  struct FloatVect2 gps_speed_m_s_ned;
436  VECT2_ASSIGN(gps_speed_m_s_ned, gps_speed_cm_s_ned.x, gps_speed_cm_s_ned.y);
437  VECT2_SDIV(gps_speed_m_s_ned, gps_speed_m_s_ned, 100.f);
438 
439  if (ins_int.hf_realign) {
440  ins_int.hf_realign = false;
441  hff_realign(gps_pos_m_ned, gps_speed_m_s_ned);
442  }
443  // run horizontal filter
444  hff_update_gps(&gps_pos_m_ned, &gps_speed_m_s_ned);
445  // convert and copy result to ins_int
446  ins_update_from_hff();
447 
448 #else /* hff not used */
449  /* simply copy horizontal pos/speed from gps */
450  INT32_VECT2_SCALE_2(ins_int.ltp_pos, gps_pos_cm_ned,
452  INT32_VECT2_SCALE_2(ins_int.ltp_speed, gps_speed_cm_s_ned,
454 #endif /* USE_HFF */
455 
457 
458  /* reset the counter to indicate we just had a measurement update */
460 }
461 #else
462 void ins_int_update_gps(struct GpsState *gps_s __attribute__((unused))) {}
463 #endif /* USE_GPS */
464 
469 #if USE_VFF_EXTENDED
470 static void agl_cb(uint8_t __attribute__((unused)) sender_id, __attribute__((unused)) uint32_t stamp, float distance) {
471  if (distance <= 0 || !(ins_int.baro_initialized)) {
472  return;
473  }
474 
475 #if USE_SONAR
476  if (distance > INS_SONAR_MAX_RANGE || distance < INS_SONAR_MIN_RANGE){
477  return;
478  }
479 #endif
480 #ifdef INS_AGL_THROTTLE_THRESHOLD
481  if(stabilization_cmd[COMMAND_THRUST] < INS_AGL_THROTTLE_THRESHOLD){
482  return;
483  }
484 #endif
485 #ifdef INS_AGL_BARO_THRESHOLD
486  if(ins_int.baro_z < -INS_SONAR_BARO_THRESHOLD){ /* z down */
487  return;
488  }
489 #endif
490 
491 #if USE_SONAR
492  vff_update_agl(-distance, VFF_R_SONAR_0 + VFF_R_SONAR_OF_M * fabsf(distance));
493 #else
494  // TODO: this assumes that you will either have sonar or other agl sensor never both
495  vff_update_agl(-distance, VFF_R_AGL);
496 #endif
497  /* reset the counter to indicate we just had a measurement update */
499 }
500 #else
501 static void agl_cb(uint8_t __attribute__((unused)) sender_id, __attribute__((unused)) uint32_t stamp, __attribute__((unused)) float distance) {}
502 #endif
503 
505 static void ins_ned_to_state(void)
506 {
510 
511 #if defined SITL && USE_NPS
512  if (nps_bypass_ins) {
514  }
515 #endif
516 }
517 
519 static void ins_update_from_vff(void)
520 {
524 }
525 
526 #if USE_HFF
527 
528 static void ins_update_from_hff(void)
529 {
536 }
537 #endif
538 
539 
540 static void accel_cb(uint8_t sender_id __attribute__((unused)),
541  uint32_t stamp, struct Int32Vect3 *accel)
542 {
543  PRINT_CONFIG_MSG("Calculating dt for INS int propagation.")
544  /* timestamp in usec when last callback was received */
545  static uint32_t last_stamp = 0;
546 
547  if (last_stamp > 0) {
548  float dt = (float)(stamp - last_stamp) * 1e-6;
549  ins_int_propagate(accel, dt);
550  }
551  last_stamp = stamp;
552 }
553 
554 static void gps_cb(uint8_t sender_id __attribute__((unused)),
555  uint32_t stamp __attribute__((unused)),
556  struct GpsState *gps_s)
557 {
558  ins_int_update_gps(gps_s);
559 }
560 
561 /* body relative velocity estimate
562  *
563  */
564 static void vel_est_cb(uint8_t sender_id __attribute__((unused)),
565  uint32_t stamp __attribute__((unused)),
566  float x, float y, float z,
567  float noise_x, float noise_y, float noise_z)
568 {
569  struct FloatVect3 vel_body = {x, y, z};
570 
571  /* rotate velocity estimate to nav/ltp frame */
572  struct FloatQuat q_b2n = *stateGetNedToBodyQuat_f();
573  QUAT_INVERT(q_b2n, q_b2n);
574  struct FloatVect3 vel_ned;
575  float_quat_vmult(&vel_ned, &q_b2n, &vel_body);
576 
577  // abi message contains an update to the horizontal velocity estimate
578 #if USE_HFF
579  struct FloatVect2 vel = {vel_ned.x, vel_ned.y};
580  struct FloatVect2 Rvel = {noise_x, noise_y};
581 
582  hff_update_vel(vel, Rvel);
583  ins_update_from_hff();
584 #else
585  if (noise_x >= 0.f)
586  {
588  }
589  if (noise_y >= 0.f)
590  {
592  }
593 
594  static uint32_t last_stamp_x = 0, last_stamp_y = 0;
595  if (noise_x >= 0.f) {
596  if (last_stamp_x > 0)
597  {
598  float dt = (float)(stamp - last_stamp_x) * 1e-6;
599  ins_int.ltp_pos.x += lround(POS_BFP_OF_REAL(dt * vel_ned.x));
600  }
601  last_stamp_x = stamp;
602  }
603 
604  if (noise_y >= 0.f)
605  {
606  if (last_stamp_y > 0)
607  {
608  float dt = (float)(stamp - last_stamp_y) * 1e-6;
609  ins_int.ltp_pos.y += lround(POS_BFP_OF_REAL(dt * vel_ned.y));
610  }
611  last_stamp_y = stamp;
612  }
613 #endif
614 
615  // abi message contains an update to the vertical velocity estimate
616  vff_update_vz_conf(vel_ned.z, noise_z);
617 
619 
620  /* reset the counter to indicate we just had a measurement update */
622 }
623 
624 /* NED position estimate relative to ltp origin
625  */
626 static void pos_est_cb(uint8_t sender_id __attribute__((unused)),
627  uint32_t stamp __attribute__((unused)),
628  float x, float y, float z,
629  float noise_x, float noise_y, float noise_z)
630 {
631 
632 #if USE_HFF
633  struct FloatVect2 pos = {x, y};
634  struct FloatVect2 Rpos = {noise_x, noise_y};
635 
636  hff_update_pos(pos, Rpos);
637  ins_update_from_hff();
638 #else
639  if (noise_x >= 0.f)
640  {
642  }
643  if (noise_y >= 0.f)
644  {
646  }
647 #endif
648 
649  vff_update_z_conf(z, noise_z);
650 
652 
653  /* reset the counter to indicate we just had a measurement update */
655 }
bool nps_bypass_ins
static const struct usb_device_descriptor dev
Definition: usb_ser_hw.c:74
static void stateSetPositionNed_i(struct NedCoor_i *ned_pos)
Set position from local NED coordinates (int).
Definition: state.h:531
Event structure to store callbacks in a linked list.
Definition: abi_common.h:65
#define VECT3_SUB(_a, _b)
Definition: pprz_algebra.h:154
int32_t z
in centimeters
static abi_event gps_ev
Definition: ins_int.c:129
static uint32_t idx
static struct Int32RMat * orientationGetRMat_i(struct OrientationReps *orientation)
Get vehicle body attitude rotation matrix (int).
#define INS_INT_AGL_ID
ABI binding for AGL.
Definition: ins_int.c:156
static void baro_cb(uint8_t sender_id, uint32_t stamp, float pressure)
Definition: ins_int.c:328
bool vf_reset
request to reset vertical filter.
Definition: ins_int.h:53
static abi_event pos_est_ev
Definition: ins_int.c:146
void vff_update_z_conf(float z_meas, float conf)
Interface for extended vertical filter (in float).
Horizontal filter (x,y) to estimate position and velocity.
#define INS_SONAR_MIN_RANGE
Default AGL sensor minimum range.
Definition: ins_ekf2.cpp:58
Periodic telemetry system header (includes downlink utility and generated code).
static struct Int32RMat * stateGetNedToBodyRMat_i(void)
Get vehicle body attitude rotation matrix (int).
Definition: state.h:1119
#define INS_VFF_VZ_R_GPS
Definition: ins_int.c:90
static abi_event vel_est_ev
Definition: ins_int.c:138
struct InsInt ins_int
global INS state
Definition: ins_int.c:161
int32_t y
in centimeters
vector in East North Up coordinates Units: meters
float qfe
Definition: ins_int.h:62
#define POS_BFP_OF_REAL(_af)
#define QUAT_INVERT(_qo, _qi)
Definition: pprz_algebra.h:627
static void ins_update_from_vff(void)
update ins state from vertical filter
Definition: ins_int.c:519
#define INT32_VECT3_ZERO(_v)
void hff_realign(struct FloatVect2 pos, struct FloatVect2 vel)
Definition: hf_float.c:602
static void stateSetAccelBody_i(struct Int32Vect3 *body_accel)
Set acceleration in Body coordinates (int).
Definition: state.h:855
Main include for ABI (AirBorneInterface).
void vff_update_vz_conf(float vz_meas, float conf)
static void gps_cb(uint8_t sender_id, uint32_t stamp, struct GpsState *gps_s)
Definition: ins_int.c:554
void ins_reset_local_origin(void)
INS local origin reset.
Definition: ins_int.c:244
#define VECT2_ASSIGN(_a, _x, _y)
Definition: pprz_algebra.h:62
float z
z-position estimate in m (NED, z-down)
struct Imu imu
global IMU state
Definition: imu.c:108
#define INT32_SPEED_OF_CM_S_DEN
vector in Latitude, Longitude and Altitude
static void stateSetSpeedNed_i(struct NedCoor_i *ned_speed)
Set ground speed in local NED coordinates (int).
Definition: state.h:763
#define GPS_FIX_3D
3D GPS fix
Definition: gps.h:39
void ltp_def_from_lla_i(struct LtpDef_i *def, struct LlaCoor_i *lla)
static struct EnuCoor_f * stateGetPositionEnu_f(void)
Get position in local ENU coordinates (float).
Definition: state.h:719
int32_t z
Down.
static struct Int32Quat * stateGetNedToBodyQuat_i(void)
Get vehicle body attitude quaternion (int).
Definition: state.h:1113
#define press_hist_len
bool baro_initialized
Definition: ins_int.h:63
void ned_of_ecef_vect_i(struct NedCoor_i *ned, struct LtpDef_i *def, struct EcefCoor_i *ecef)
Rotate a vector from ECEF to NED.
#define INS_INT_POS_ID
Definition: ins_int.c:144
struct EcefCoor_i ecef
Reference point in ecef.
int32_t hmsl
Height above mean sea level in mm.
void ins_int_update_gps(struct GpsState *gps_s)
Definition: ins_int.c:387
static struct FloatQuat * stateGetNedToBodyQuat_f(void)
Get vehicle body attitude quaternion (float).
Definition: state.h:1131
int32_t alt
in millimeters above WGS84 reference ellipsoid
Vertical filter (in float) estimating altitude, velocity and accel bias.
Roation quaternion.
static abi_event agl_ev
The agl ABI event.
Definition: ins_int.c:158
abi_event baro_ev
Definition: ins_int.c:114
int32_t y
East.
static void stateSetAccelNed_i(struct NedCoor_i *ned_accel)
Set acceleration in NED coordinates (int).
Definition: state.h:986
#define INT32_VECT2_SCALE_2(_a, _b, _num, _den)
static void send_ins_z(struct transport_tx *trans, struct link_device *dev)
Definition: ins_int.c:174
int32_t hmsl
height above mean sea level (MSL) in mm
Definition: gps.h:94
float x
in meters
#define INS_INT_BARO_ID
Definition: ins_int.c:110
static void send_ins_ref(struct transport_tx *trans, struct link_device *dev)
Definition: ins_int.c:180
data structure for GPS information
Definition: gps.h:87
void ins_int_init(void)
Definition: ins_int.c:198
#define INT32_POS_OF_CM_NUM
void hff_update_vel(struct FloatVect2 vel, struct FloatVect2 Rvel)
Definition: hf_float.c:806
static void stateSetLocalOrigin_i(struct LtpDef_i *ltp_def)
Set the local (flat earth) coordinate frame origin (int).
Definition: state.h:457
Device independent GPS code (interface)
Paparazzi atmospheric pressure conversion utilities.
#define ACCEL_FLOAT_OF_BFP(_ai)
bool ltp_initialized
Definition: ins_int.h:41
int32_t x
North.
struct EcefCoor_i ecef_pos
position in ECEF in cm
Definition: gps.h:91
#define INS_INT_IMU_ID
ABI binding for IMU data.
Definition: ins_int.c:121
PRINT_CONFIG_MSG("USE_INS_NAV_INIT defaulting to TRUE")
static void ins_ned_to_state(void)
copy position and speed to state interface
Definition: ins_int.c:505
struct LlaCoor_i lla
Reference point in lla.
#define DefaultPeriodic
Set default periodic telemetry.
Definition: telemetry.h:66
int32_t lon
in degrees*1e7
#define INS_MAX_PROPAGATION_STEPS
maximum number of propagation steps without any updates in between
Definition: ins_int.c:96
static void pos_est_cb(uint8_t sender_id, uint32_t stamp, float x, float y, float z, float noise_x, float noise_y, float noise_z)
Definition: ins_int.c:626
void ins_init_origin_i_from_flightplan(struct LtpDef_i *ltp_def)
initialize the local origin (ltp_def in fixed point) from flight plan position
Definition: ins.c:39
void float_quat_vmult(struct FloatVect3 *v_out, struct FloatQuat *q, const struct FloatVect3 *v_in)
rotate 3D vector by quaternion.
static void accel_cb(uint8_t sender_id, uint32_t stamp, struct Int32Vect3 *accel)
Definition: ins_int.c:540
uint16_t f
Camera baseline, in meters (i.e. horizontal distance between the two cameras of the stereo setup) ...
Definition: wedgebug.c:204
Paparazzi fixed point math for geodetic calculations.
void int32_quat_vmult(struct Int32Vect3 *v_out, struct Int32Quat *q, struct Int32Vect3 *v_in)
rotate 3D vector by quaternion.
void vff_init_zero(void)
bool hf_realign
request to realign horizontal filter.
Definition: ins_int.h:48
void hff_propagate(void)
Definition: hf_float.c:479
float zdotdot
z-acceleration in m/s^2 (NED, z-down)
#define INS_SONAR_MAX_RANGE
Default AGL sensor maximum range.
Definition: ins_ekf2.cpp:64
Inertial Measurement Unit interface.
void vff_update_agl(float z_meas, float conf)
struct NedCoor_i ltp_pos
Definition: ins_int.h:56
void vff_propagate(float accel, float dt)
Propagate the filter in time.
struct HfilterFloat hff
Definition: hf_float.c:124
struct OrientationReps body_to_imu
rotation from body to imu frame
Definition: imu.h:49
void hff_update_pos(struct FloatVect2 pos, struct FloatVect2 Rpos)
Update position.
Definition: hf_float.c:711
struct NedCoor_i ltp_accel
Definition: ins_int.h:58
#define INS_BARO_MAX_INIT_VAR
default barometer to use in INS
Definition: ins_int.c:105
float baro_z
z-position calculated from baro in meters (z-down)
Definition: ins_int.h:61
struct LtpDef_i ned_origin_i
Definition of the local (flat earth) coordinate system.
Definition: state.h:166
#define VECT2_SDIV(_vo, _vi, _s)
Definition: pprz_algebra.h:104
uint32_t propagation_cnt
number of propagation steps since the last measurement update
Definition: ins_int.h:43
float variance_f(float *array, uint32_t n_elements)
Compute the variance of an array of values (float).
Definition: pprz_stat.c:139
void vff_realign(float z_meas)
float x
Definition: hf_float.h:41
float y
Definition: hf_float.h:45
Ins implementation state (fixed point)
Definition: ins_int.h:39
#define INT32_POS_OF_CM_DEN
void hff_init(float init_x, float init_xdot, float init_y, float init_ydot)
Definition: hf_float.c:274
float xdot
Definition: hf_float.h:42
struct VffExtended vff
float zdot
z-velocity estimate in m/s (NED, z-down)
void ned_of_ecef_point_i(struct NedCoor_i *ned, struct LtpDef_i *def, struct EcefCoor_i *ecef)
Convert a point from ECEF to local NED.
vector in North East Down coordinates
#define INS_INT_VEL_ID
ABI binding for VELOCITY_ESTIMATE.
Definition: ins_int.c:136
float xdotdot
Definition: hf_float.h:43
static abi_event accel_ev
Definition: ins_int.c:123
#define INT32_SPEED_OF_CM_S_NUM
float ydot
Definition: hf_float.h:46
Statistics functions.
static void send_ins(struct transport_tx *trans, struct link_device *dev)
Definition: ins_int.c:166
General stabilization interface for rotorcrafts.
int32_t stabilization_cmd[COMMANDS_NB]
Stabilization commands.
Definition: stabilization.c:32
void int32_rmat_transp_vmult(struct Int32Vect3 *vb, struct Int32RMat *m_b2a, struct Int32Vect3 *va)
rotate 3D vector by transposed rotation matrix.
void vff_update(float z_meas)
Definition: vf_float.c:189
int32_t x
in centimeters
struct LtpDef_i ltp_def
Definition: ins_int.h:40
void vff_update_baro(float z_meas)
rotation matrix
struct EcefCoor_i ecef_vel
speed ECEF in cm/s
Definition: gps.h:95
#define SPEED_BFP_OF_REAL(_af)
struct LlaCoor_i lla_pos
position in LLA (lat,lon: deg*1e7; alt: mm over ellipsoid)
Definition: gps.h:92
static float pprz_isa_height_of_pressure(float pressure, float ref_p)
Get relative altitude from pressure (using simplified equation).
Definition: pprz_isa.h:102
#define ACCEL_BFP_OF_REAL(_af)
void hff_update_gps(struct FloatVect2 *pos_ned, struct FloatVect2 *speed_ned)
Definition: hf_float.c:541
int32_t lat
in degrees*1e7
#define INS_VFF_R_GPS
Definition: ins_int.c:86
uint8_t fix
status of fix
Definition: gps.h:107
static void agl_cb(uint8_t sender_id, uint32_t stamp, float distance)
agl_cb This callback handles all estimates of the height of the vehicle above the ground under it Thi...
Definition: ins_int.c:501
#define INS_INT_GPS_ID
Definition: ins_int.c:127
float y
in meters
Rotation quaternion.
struct GpsState gps
global GPS state
Definition: gps.c:69
void ins_int_propagate(struct Int32Vect3 *accel, float dt)
Definition: ins_int.c:283
struct NedCoor_i ltp_speed
Definition: ins_int.h:57
void sim_overwrite_ins(void)
void ltp_def_from_ecef_i(struct LtpDef_i *def, struct EcefCoor_i *ecef)
struct State state
Definition: state.c:36
int8_t register_periodic_telemetry(struct periodic_telemetry *_pt, uint8_t _id, telemetry_cb _cb)
Register a telemetry callback function.
Definition: telemetry.c:46
void ins_reset_altitude_ref(void)
INS altitude reference reset.
Definition: ins_int.c:266
static void vel_est_cb(uint8_t sender_id, uint32_t stamp, float x, float y, float z, float noise_x, float noise_y, float noise_z)
Definition: ins_int.c:564
float ydotdot
Definition: hf_float.h:47
unsigned char uint8_t
Typedef defining 8 bit unsigned char type.
Definition: vl53l1_types.h:98
INS for rotorcrafts combining vertical and horizontal filters.
#define VFF_R_AGL
Definition: ins_int.c:59
unsigned int uint32_t
Typedef defining 32 bit unsigned int type.
Definition: vl53l1_types.h:78
#define GpsFixValid()
Definition: gps.h:43