Monitoring of an athlete's kinematics both in training and in competition is important in the development and implementation of new approaches towards performance improvement as well as injury analysis and prevention.
Motion sensing devices are frequently used in order to determine the motion of an athlete. For example, such devices may sense motion parameters such as acceleration, angular rates, velocity, stride distance, total distance, speed, stride rate, and the like, for use in the training and evaluation of athletes, and the rehabilitation of the injured.
There are a number of solutions that measure kinematic parameters in one plane (X/Y) and the orientation (pitch) of an athlete's foot. These systems provide valuable insight into the biomechanics of motion, but fail to resolve the full 6D movement of the athlete's foot. 6D—in the context of stride based kinematics—representing both the position (X, Y, Z) as well as the orientation (pitch, roll, yaw) of the athlete's foot.
These designs, having focused on XY-plane stride kinematics, are implemented as single foot solutions. This assumes left/right symmetry, which for some metrics is safe, but for many, is an invalid assumption. Metrics like stride rate, velocity, even contact time (to some degree) will tend to be highly symmetric. However, pronation velocity, pronation angle, even pitch at footstrike, among others, can be radically different between an athlete's right and left sides. The disclosed detachable measurement system may be optionally implemented as either single (right or left) or both feet—providing full 6D space/orientation kinematic parameters in each combination.
Embodiments of the present invention provide a system for determining athletic kinematic characteristics. The system includes an inertial sensor, a processing system, and a wireless transceiver. The inertial sensor may be coupled with a user's footwear in order to generate one or more signals corresponding to the motion of the user's foot/feet. The processing system is in communication with the inertial sensor and is programmed to use the one or more signals to determine one or more kinematic characteristics. The present invention measures various parameters about each individual stride rather than assuming a given fixed rate. The stride based kinematic characteristics may include, but are not limited to, pitch, roll, yaw, vertical position, horizontal position, horizontal velocity, vertical velocity, distance traveled, foot strike, foot strike classification, toe off, contact time, stride rate, stride length, rate of pronation, maximum pronation, rate of plantarflexion and dorsiflexion, swing velocity, and pitch-roll signature.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not necessarily restrictive of the invention claimed. The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the invention and together with the general description, serve to explain the principles of the invention.
The following detailed description of embodiments of the invention references the accompanying drawings. The embodiments are intended to describe aspects of the invention in sufficient detail to enable those skilled in the art to practice the invention. Other embodiments can be utilized and changes can be made without departing from the scope of the claims. The following detailed description is, therefore, not to be taken in a limiting sense. The scope of the present invention is defined only by the appended claims, along with the full scope of equivalents to which such claims are entitled.
In this description, references to “one embodiment”, “an embodiment”, or “embodiments” mean that the feature or features being referred to are included in at least one embodiment of the technology. Separate references to “one embodiment”, “an embodiment”, or “embodiments” in this description do not necessarily refer to the same embodiment and are also not mutually exclusive unless so stated and/or except as will be readily apparent to those skilled in the art from the description. For example, a feature, structure, method, etc. described in one embodiment may also be included in other embodiments, but is not necessarily included. Thus, the present technology can include a variety of combinations and/or integrations of the embodiments described herein.
The information to permit stride based kinematic analysis is obtained via a suitable Motion Processing Unit (MPU) 20—comprised of sensors, preferably a 3D Accelerometer 21, a 3D Gyroscope 22, and (optionally) a 3D Compass 23 as shown in
As shown in
The Corrected Quaternion 32 shown in
Step 1: Convert Gyroscope 22 angular rate to a quaternion representation 38, where w(t) is the angular rate and q(t) is the normalized quaternion.
dq(t)/dt=½w(t)*q(t)
Step 2: Convert Accelerometer data to world coordinates. This means using the Quaternion above to get the appropriate coordinate system in world-frame motion. Here Ab(t) is in the body coordinates of the device 1, while Aw(t) is in world-frame.
A
w(t)=q(t)*Ab(t)*q(t)′
Step 3: Create an acceleration measurement feedback quaternion 39 as below.
qf(t)=[0Awy(t)−Awx(t)0]*gain
Step 4: Once converted to world coordinates, accelometer feedback and gain is used to generate a feedback quaternion which is then added to the previous quaternion along with the gyro generated quaternion. The result is a Corrected Quaternion 32 that will track the gyroscope measured data, but will drift towards the accelerometer measurement, according to the value chosen for gain. Similarly, compass data can be added to the yaw component of the quaternion in order to correct for drift in yaw.
Shown in
We associate a quaternion with a rotation around an axis by the expressions:
q
0=cos(α/2)
q
1=sin(α/2)cos(βx)
q
2=sin(α/2)cos(βy)
q
3=sin(α/2)cos(βz)
where α is a simple rotation angle (the value in radians of the angle of rotation) and cos(βx), cos(βy) and cos(βz) are the “direction cosines” locating the axis of rotation (Euler's Theorem). From this we can derive the following rotation matrix:
q02+q12−q22−q32 2(q1q2−q0q3) 2(q0q2−q1q3) 2(q1q2+q0q3) q02−q12+q22−q322(q2q3−q0q1) 2(q1q3−q0q2) 2(q0q1+q2q3) q02−q12−q22+q32
Pitch 35, Roll 36, and Yaw 37 can thus be computed by the following equations.
Θ=a tan 2(2(q0q1+q2q3), 1−2(q12+q22))
Φ=arcsin(2(q0q2−q3q1))
Ψ=a tan 2(2(q0q3+q1q2), 1−2(q22+q32))
Having determined device orientation, it is now possible to determine device Position and Velocity (X,Y,Z) 33 by integrating the Gravity Corrected Accelerations (X,Y,Z) 30 as follows:
These equations are integrated once to determine horizontal and vertical velocity, and twice to determine the stride length, and the vertical displacement of the foot. While the above calculations show corrections for the pitch (X/Y) axis, it is also understood that similar corrections may be made for both roll and yaw axes as well.
Referring to
Step 1: Locate the pitch gyro peak=max(pitch gyro data) since last detected pitch gyro peak.
Step 2: Determine Foot Strike (
Step 3: Determine Toe Off (
Step 4: Determine Maximum Pronation Angle (
Steps 5-N: Continue locating all other Stride Based Metrics—including, but not limited to:
The above kinematic metrics being recorded in Data Storage Memory 28 and optionally transmitted in real time via Wireless Transceiver 29 using, for example, wireless protocols such as ANT, ANT+, or Bluetooth Low Energy (BT Smart), as shown in
It is understood that in ideal (laboratory) environments, the sensors used to collect the kinematic parameters described above can operate with few error sources. That the data is ‘accurate’ as a result of the constrained environmental and operational settings. However, when the device is preferably used in non-laboratory settings, such as training and competition, the system must be capable of maintaining accuracy in order to continue to correctly determine the same high quality kinematic metrics as disclosed above. In order to do so, it is required that the device limitations be well understood, and compensated for accordingly.
The output of rate gyroscopes is rotational rate, and to obtain a relative change in angle, a single integration on the gyro outputs must be performed. Error in gyro bias (the output of the gyro when rotation is zero) leads to an error that increases with integration time. Methods must be taken to compensate for these bias errors, which are caused by drift due to time and temperature, and by noise.
Common methods of compensation involve the use of other sensors, such as accelerometers for tilt angle, and compasses for heading. Alternately, changes in bias may be sensed when the device is not moving (i.e. pause during a run). No motion is detected by looking at peak deviation in gyro output during a relatively short timeframe, such as two seconds. If the peak-to-peak signal is below a predetermined threshold, it is determined that the device is stationary, and the average gyro output during that time becomes the new bias setting.
Note that accelerometers and compass sensors also have bias drift, but since accelerometers provide tilt angle directly (without integration) by measuring gravity, and since compass sensors provide heading information directly by measuring the earth's magnetic field, bias errors in these sensors are not integrated when providing tilt angle or heading. However, when double integrating the output of an accelerometer to provide distance or when single integrating its output to provide velocity, the bias errors of the accelerometer become important.
Magnetic sensors (also known as compass sensors) are used to determine heading (yaw orientation) using magnetic north as a reference. The value of compass sensors is that they provide absolute heading information using a known reference (magnetic north). This is in contrast with gyros, which provide relative outputs that can accurately detect how far a device has rotated. Additionally, the compass sensors are typically only used for rotational information around the yaw axis, while gyros provide information around the X, Y, and Z axes (pitch, roll, and yaw).
Magnetic sensors respond to more than just the earth's magnetic field (typically ranging from 30 microteslas to over 60 microteslas). They also respond to interference, such as RF signals (caused by cell phones, radio towers, etc.) and to magnetic fields caused by magnets, such as those in cell phones and headphones. Compasses are often used in combination with gyroscopes, where the gyroscopes provide a heading signal for faster motions, and the filtered compass output provides a heading signal with a longer time constant to be used for bias and heading compensation. Additionally, since the earth's magnetic field is not perfectly parallel to the surface of the earth, its angle varies with position on the Earth, accelerometers are used in conjunction with compass sensors to provide tilt compensation.
Another source of error may arise from the arbitrary mounting angle of the detachable motion sensor 10. While it is possible to vertically align the +Y axis as shown in
Zero offset correction of depth is one of the first considerations in analyses of diving behaviour data from time-depth recorders (TDRs). Pressure transducers in TDRs often “drift” over time due to temperature changes and other factors, so that recorded depth deviates from actual depth over time at unpredictable rates.
For diving animals, such as marine mammals and seabirds, the problem of zero offset correction is simplified by the cyclical return to or from the surface as study animals perform their dives throughout the deployment period, thereby providing a reference for calibration (The short period where the foot is flat on the ground during each stride is the equivalent in kinematic stride analysis).
The method consists of recursively smoothing and filtering the input time series using moving quantiles. It uses a sequence of window widths and quantiles, and starts by filtering the time series using the first window width and quantile in the specified sequences. The second filter is applied to the output of the first one, using the second specified window width and quantile, and so on. In most cases, two steps are sufficient to detect the surface signal in the time series: the first to remove noise as much as possible, and the second to detect the surface level. Depth is corrected by subtracting the output of the last filter from the original.
Using the above dual filter technique, the ‘corrupted’ roll and yaw data can be recursively filtered as depicted in
Inherent in the biomechanics of humans is intrinsic asymmetry which can manifest itself in different ways which may adversely affect performance and even lead to injury. The ability of the disclosed invention to measure and record the motion of an athlete can provide valuable insight into these asymmetries when the motion sensing system 10 is affixed to both the athlete's left and right feet. Information particularly between Foot Strike (
When both right and left data are to be simultaneously recorded, the motion sensing system 10 on the left foot may be preferably designated as a slave device, forwarding its stride based metrics to the master device on the right foot, which will aggregate the data from the two systems, then record and/or transmit the information via wireless interface.
Using the kinematic metrics collected by the system, it is possible to compute a metric that can be used to represent the intensity (runScore) of an activity. Specifically, using an equation of the form:
runScore=a*Pace+b*StrideRate+c*PronationExcursion+d*Maximum Pronation Velocity+e*ImpactGs+f*BrakingGs+ . . .
This intensity metric can then be used to quickly visualize the ‘stress’ of a given run (such as
The intensity formula may also be expanded to further include other non-kinematic metrics, such as physiological parameters like: heart rate, heart rate variability, oxygen consumption, and perceived exertion.
Using the data collected by the system, it is possible to interpret plots (such as
Again, using the kinematic data collected by the system, it is possible to visualize the change of kinematic parameters (such as ImpactGs, BrakingGs, Maximum Pronation Excursion) on a given pair of shoes as mileage increases (such as
Using the data collected by the system, it is possible to aggregate kinematic metrics from a large population of users. Enabling specific demographic comparisons to be made, such as: age group, weight, competitive level, type of terrain, length of run, and average pace. Such aggregate data can then be used to enable injury correlations to the collected kinematic data, looking for trends in individual and combinations of metrics, such as ImpactGs and Pronation Velocity. The aggregate data can also be gathered for specific events which have a large number of participants (such as Boston and NYC Marathons), where the mean and variance of key kinematic metrics can be compared over the course of that specific event (shown in
As described above, the motion system shown in
The device may be battery powered; this requires that the primary components and associated circuits possess low-power consumption characteristics.
The sensor is mounted on the foot or shoe and will thus be subjected to large impact forces and abuse. It is necessary that the sensors be rugged and durable to be able to survive in this environment.
The linearity, repeatability and noise levels must be such that the accuracy of measurement is acceptable for the application.
The motion processing units used in the development work of this invention are manufactured by InvenSense (part no.'s MPU-9150 and MPU-9250). These devices are constructed using MEMS techniques to build the transducers into a silicon chip. This accounts for the small size, low power consumption and accuracy of the devices.
The invention described herein is not limited to the above mentioned sensor family. Other MEMS accelerometers, gyroscopes, and compasses are currently produced or are under development by different manufacturers and could be considered for this purpose.
The integrated application processor and wireless transceiver used in the development work of this invention is manufactured by Nordic Semiconductor (part no.'s nRF51422, nRF51822, and nRF51922). These devices comprise an ARM Cortex-MO level microcontroller with 256 kB of embedded flash program memory and 16 kB of RAM.
The data storage memory used in the development of this invention is manufactured by Macronix (part no. MX25L25635EZNI). This device is a Serial Flash containing 256 Mbit (32 Mbyte) of non-volatile storage for data storage and retention.
Although the invention has been described with reference to various exemplary embodiments illustrated in the attached drawing figures, it is noted that equivalents may be employed and substitutions made herein without departing from the scope of the invention as recited in the claims. Having thus described embodiments of the invention, what is claimed as new and desired to be protected by patent includes the following:
This application claims benefit of Provisional Application 61/890,299 filed Oct. 13, 2013 entitled “Detachable Wireless Motion System for Human Kinematic Analysis”.
Number | Date | Country | |
---|---|---|---|
20160100801 A1 | Apr 2016 | US |
Number | Date | Country | |
---|---|---|---|
61890299 | Oct 2013 | US |