System and method for measuring mobility metrics of a joint using magnetometer-free inertial measurement units

Abstract

A method and a system include a first inertial measurement unit (IMU) is positioned on a first member of a joint and a second IMU is positioned on a second member of the joint, where the first IMU and the second IMU each lack a magnetometer unit. A processing device receives first IMU sensor data and second IMU sensor data. The first IMU sensor data is inputted into a quaternion computation algorithm (QCA) to compute a first quaternion. The first quaternion and the second IMU sensor data are inputted into the QCA to compute a one second quaternion, where the inputting of the first quaternion imposes yaw constraints. A calibration parameter algorithm determines a mobility metric of the joint based on the first and the second quaternions, and an orientation and position of the first IMU on the first member and of the second IMU on the second member.

1 . A system, comprising:
 at least one first inertial measurement unit (IMU) is positioned on a first member of a joint;   at least one second inertial measurement unit (IMU) is positioned on a second member of the joint;
 wherein each of the at least one first IMU and the at least one second IMU lacks a magnetometer unit; 
 wherein the at least one first IMU and the at least one second IMU respectively output at least one first IMU sensor data and at least one second IMU sensor data; 
   at least one processing device configured to determine at least one mobility metric of the joint by executing programming instructions that cause the at least one processing device to:
 receive, from the at least one first IMU, the at least one first IMU sensor data; 
 receive, from the at least one second IMU, the at least one second IMU sensor data; 
 input the at least one first IMU sensor data into a quaternion computation algorithm to compute at least one first quaternion; 
 input the at least one first quaternion and the at least one second IMU sensor data into the quaternion computation algorithm to compute at least one second quaternion;
 wherein the inputting of the at least one first quaternion imposes a yaw constraint that reduces an uncertainty in:
 at least one first yaw axis value in a coordinate frame of the at least one first IMU, and 
 at least one second yaw axis value in a coordinate frame of the at least one second IMU; and 
 
 
 utilize a calibration parameter algorithm to determine the at least one mobility metric of the joint based on:
 the at least one first quaternion, 
 the at least one second quaternion, 
 at least one first orientation and at least one first position of the at least one first IMU on the first member of the joint, and 
 at least one second orientation and at least one second position of the at least one second IMU on the second member of the joint. 
 
   
     
     
         2 . The system according to  claim 1 , wherein the at least one mobility metric of the joint further comprises a joint position, a joint angle, a step length, a stride length, a cadence, or any combination thereof. 
     
     
         3 . The system according to  claim 1 , further comprising a memory;
 wherein the memory is configured to store a configuration database; and   wherein the at least one processing device is configured to store system variables, preconfigured nominal frequencies, sensor information, user information, or any combination thereof in the configuration database.   
     
     
         4 . The system according to  claim 1 , further comprising a mobile device;
 wherein the at least one first IMU is configured to wirelessly transmit the at least one first IMU sensor data to the mobile device;   wherein the at least one second IMU is configured to wirelessly transmit the at least one second IMU sensor data to the mobile device;   wherein the mobile device comprises the at least one processing device; and   wherein the at least one processing device is configured to determine the at least one mobility metric of the joint based at least in part on the at least one first IMU sensor data and the at least one second IMU sensor data.   
     
     
         5 . The system according to  claim 4 , wherein the at least one processing device is configured to receive the at least one first IMU sensor data and the at least one second IMU sensor data over a wireless communication network. 
     
     
         6 . The system according to  claim 4 , wherein the at least one processing device is configured to synchronize the at least one first IMU sensor data and the at least one second IMU sensor data by calculating a local master time difference to compensate for an initial time delay between times that the at least one first IMU and the at least one second IMU start to collect data. 
     
     
         7 . The system according to  claim 6 , wherein the local master time difference is based at least in part on a master clock of the mobile device. 
     
     
         8 . The system according to  claim 1 , wherein the at least one first IMU comprises at least one first accelerometer and at least one first gyroscope;
 wherein the at least one first IMU sensor data comprises at least one first accelerometer sensor data and at least one first gyroscope sensor data;   wherein the at least one second IMU comprises at least one second accelerometer and at least one second gyroscope;   wherein the at least one second IMU sensor data comprises at least one second accelerometer sensor data and at least one second gyroscope sensor data; and   wherein the at least one processing device is configured to use the quaternion computation algorithm to:
 calculate from the at least one first accelerometer sensor data and the at least one first gyroscope sensor data, an estimation of the at least one first quaternion;
 wherein the at least one first quaternion represents at least one first pitch axis value, at least one first roll axis value, and the at least one first yaw axis value in the coordinate frame of the at least one first IMU; 
 
 calculate from the at least one second accelerometer sensor data and the at least one second gyroscope sensor data, an estimation of the at least one second quaternion;
 wherein the at least one second quaternion represents at least one second pitch axis value, at least one second roll axis value, and the at least one second yaw axis value in the coordinate frame of the at least one second IMU; and 
 
 impose the yaw constraint between the at least one first yaw axis value and the at least one second yaw axis value to improve:
 the estimation of the at least one first quaternion and 
 the estimation of the at least one second quaternion. 
 
   
     
     
         9 . The system according to  claim 1 , wherein the at least one processing device is configured to use the calibration parameter algorithm to compute the at least one mobility metric of the joint by calibrating training data during movements of the joint to determine:
 the at least one first orientation and the at least one first position of the at least one first IMU on the first member of the joint, and   the at least one second orientation and the at least one second position of the at least one second IMU on the second member of the joint.   
     
     
         10 . The system according to  claim 1 , wherein the joint is a knee joint of a leg, the first member is a thigh of the leg, and the second member is a calf of the leg. 
     
     
         11 . A method, comprising:
 receiving, by at least one processing device, from at least one first inertial measurement unit (IMU), at least one first IMU sensor data;
 wherein the at least one first IMU is positioned on a first member of a joint; 
   receiving, by the at least one processing device, from at least one second inertial measurement unit (IMU), at least one second IMU sensor data;
 wherein the at least one second IMU is positioned on a second member of the joint; 
 wherein each of the at least one first IMU and the at least one second IMU lacks a magnetometer unit; 
 wherein the at least one first IMU and the at least one second IMU respectively output the at least one first IMU sensor data and the at least one second IMU sensor data; 
 inputting, by the at least one processing device, the at least one first IMU sensor data into a quaternion computation algorithm to compute at least one first quaternion; 
 inputting, by the at least one processing device, the at least one first quaternion and the at least one second IMU sensor data into the quaternion computation algorithm to compute at least one second quaternion;
 wherein the inputting of the at least one first quaternion imposes a yaw constraint that reduces an uncertainty in:
 at least one first yaw axis value in a coordinate frame of the at least one first IMU and 
 at least one second yaw axis value in a coordinate frame of the at least one second IMU; and 
 
 
 utilizing, by the at least one processing device, a calibration parameter algorithm to determine at least one mobility metric of the joint based on:
 the at least one first quaternion, 
 the at least one second quaternion, 
 at least one first orientation and at least one first position of the at least one first IMU on the first member of the joint, and 
 at least one second orientation and at least one second position of the at least one second IMU on the second member of the joint. 
 
   
     
     
         12 . The method according to  claim 11 , wherein the at least one mobility metric of the joint further comprises a joint position, a joint angle, a step length, a stride length, a cadence, or any combination thereof. 
     
     
         13 . The method according to  claim 11 , further comprising storing, by the at least one processing device, system variables, preconfigured nominal frequencies, sensor information, user information, or any combination thereof in a configuration database in a memory. 
     
     
         14 . The method according to  claim 11 , further comprising a mobile device;
 wherein the at least one first IMU is configured to wirelessly transmit the at least one first IMU sensor data to the mobile device;   wherein the at least one second IMU is configured to wirelessly transmit the at least one second IMU sensor data to the mobile device;   wherein the mobile device comprises the at least one processing device; and   further comprising determining, by the at least one processing device, the at least one mobility metric of the joint based at least in part on the at least one first IMU sensor data and the at least one second IMU sensor data.   
     
     
         15 . The method according to  claim 14 , further comprising receiving, by the at least one processing device, the at least one first IMU sensor data and the at least one second IMU sensor data over a wireless communication network. 
     
     
         16 . The method according to  claim 14 , further comprising synchronizing, by the at least one processing device, the at least one first IMU sensor data and the at least one second IMU sensor data by calculating a local master time difference to compensate for an initial time delay between times that the at least one first IMU and the at least one second IMU start to collect data. 
     
     
         17 . The method according to  claim 16 , wherein the local master time difference is based at least in part on a master clock of the mobile device. 
     
     
         18 . The method according to  claim 11 , wherein the at least one first IMU comprises at least one first accelerometer and at least one first gyroscope;
 wherein the at least one first IMU sensor data comprises at least one first accelerometer sensor data and at least one first gyroscope sensor data;   wherein the at least one second IMU comprises at least one second accelerometer and at least one second gyroscope;   wherein the at least one second IMU sensor data comprises at least one second accelerometer sensor data and at least one second gyroscope sensor data; and   
       further comprising:
 calculating, by the at least one processing device, from the at least one first accelerometer sensor data and the at least one first gyroscope sensor data using the quaternion computation algorithm, an estimation of the at least one first quaternion;
 wherein the at least one first quaternion represents at least one first pitch axis value, at least one first roll axis value, and the at least one first yaw axis value in the coordinate frame of the at least one first IMU; 
 
 calculating, by the at least one processing device, from the at least one second accelerometer sensor data and the at least one second gyroscope sensor data using the quaternion computation algorithm, an estimation of the at least one second quaternion;
 wherein the at least one second quaternion represents at least one second pitch axis value, at least one second roll axis value, and the at least one second yaw axis value in the coordinate frame of the at least one second IMU; and 
 
 imposing, by the at least one processing device, using the quaternion computation algorithm, the yaw constraint between the at least one first yaw axis value and the at least one second yaw axis value to improve:
 the estimation of the at least one first quaternion and 
 the estimation of the at least one second quaternion. 
 
 
     
     
         19 . The method according to  claim 11 , wherein the at least one processing device is configured to use the calibration parameter algorithm to compute the at least one mobility metric of the joint by calibrating training data during movements of the joint to determine:
 the at least one first orientation and the at least one first position of the at least one first IMU on the first member of the joint, and   the at least one second orientation and the at least one second position of the at least one second IMU on the second member of the joint.   
     
     
         20 . The method according to  claim 11 , wherein the joint is a knee joint of a leg, the first member is a thigh of the leg, and the second member is a calf of the leg.