Source: https://patents.google.com/patent/WO2008143738A1/en
Timestamp: 2019-10-16 16:28:16
Document Index: 723034734

Matched Legal Cases: ['Application No. 60', 'Application No.\n60', 'Application No.\n60', 'Application No. 60', 'Application No. 60', 'Application No. 60']

WO2008143738A1 - Newtonian physical activity monitor - Google Patents
WO2008143738A1
WO2008143738A1 PCT/US2008/003745 US2008003745W WO2008143738A1 WO 2008143738 A1 WO2008143738 A1 WO 2008143738A1 US 2008003745 W US2008003745 W US 2008003745W WO 2008143738 A1 WO2008143738 A1 WO 2008143738A1
PCT/US2008/003745
WO2008143738A8 (en
Ultimate Balance, Inc.
2007-05-18 Priority to US93074907P priority Critical
2007-05-18 Priority to US60/930,749 priority
2007-06-11 Priority to US93421607P priority
2007-06-11 Priority to US60/934,216 priority
2007-06-18 Priority to US60/936,117 priority
2007-06-18 Priority to US93611707P priority
2008-03-21 Application filed by Ultimate Balance, Inc. filed Critical Ultimate Balance, Inc.
2008-11-27 Publication of WO2008143738A1 publication Critical patent/WO2008143738A1/en
2009-04-09 Publication of WO2008143738A8 publication Critical patent/WO2008143738A8/en
230000037120 immobility Effects 0 description 2
230000037081 physical activity Effects 0 abstract description title 33
230000000704 physical effects Effects 0 abstract description title 33
230000036544 posture Effects 0 description 43
An improved apparatus and methods of posture and physical activity monitoring. The apparatus is physically mountable to or associated with an object or person, includes a multi-axis accelerometer, and derives measurements of posture and of acceleration. Methods are disclosed which provide improved estimations of posture, acceleration, energy expenditure, movement characteristics and physical activity, detect the influence of externally- caused motion, and permit automatic calibration of the apparatus in the field.
TITLE OF THE INVENTION NEWTONIAN PHYSICAL ACTIVITY MONITOR
This application claims benefit of U.S. Provisional Patent
Application No. 60/930,749 filed May 18, 2007 entitled PHYSICAL
ACTIVITY MONITOR, U.S. Provisional Patent Application No.
60/934,216 filed June 11, 2007 entitled PHYSICAL ACTIVITY AND POSTURE MONITOR, and U.S. Provisional Patent Application No.
60/936,117 filed June 18, 2007 entitled PHYSICAL ACTIVITY AND
POSTURE MONITOR.
The present invention relates generally to the field of physical activity monitoring, and more specifically to apparatus and methods of monitoring and measuring posture and physical activity using acceleration data for research, training, monitoring, sports, leisure, therapeutic and rehabilitation purposes . Uni-axial and tri-axial accelerometers have been used for many years in wearable devices that estimate physical activity.
The tri-axial devices provide acceleration data for three orthogonal axes, plus the vector magnitude. Many conventional systems integrate the accelerometer data over a time window (an "activity count"), and use linear or higher-order formulas or look-up tables to estimate Energy Expenditure (EE) , expressed in units such as METS (multiple of basal metabolic rate), or calories. The offset and scaling factors of the conversion equations are often determined by statistical curve-fitting, using population-sample accelerometer data collected concurrently with reference EE measurements, such as from oxygen consumption or doubly-labeled water. The applicant has recognized that existing systems have shortcomings, including one or more of the following:
3) Existing systems that use a low-pass filter to suppress the acceleration signal perform the filtering function on both the direction and magnitude of the raw signal . This leads to inaccuracies in the resulting calculation of acceleration, as the filter can allow feed-through of either component of the raw signal in either or both directions. 4) Existing systems which subtract a calibrated value of gravity from the total signal to determine the magnitude of acceleration do not correct the calculation for differences in direction between gravity and the user's acceleration, causing errors in the resulting value of acceleration.
6) Existing systems require precise attachment, aligned with the wearer's vertical, anterior-posterior, or medio-lateral axes, and do not correct for alignment errors between the sensor axes and those of the wearer, or for changes in the alignment during the measurement process, or for changes in the wearer's posture during the measurement process. 7) Existing systems do not have a means to automatically detect the influence of modes of conveyance such as elevators, automobiles, or airplanes on measured acceleration.
8) Existing systems correlate step counts or total dynamic acceleration (either uni-axial, tri-axial vector, or tri-axial sum-of-axes) to Energy Expenditure (EE) , and do not discriminate between physical activity related to changes in potential energy
(EEp, i.e. displacement against gravity), and physical activity related to changes in kinetic energy (EEK, i.e. acceleration against inertia) . Since the relationships between acceleration and these two types of energy are different, this leads to measurement errors. In addition, separate measures of these two types of physical activity can be useful for purposes of monitoring, analysis, and feedback to the user.
10) Existing systems analyze the frequency or pattern of acceleration to estimate activity levels. However, pattern analysis is prone to misinterpretation, especially in daily- living situations where movements can be complex and unpredictable .
12) Existing systems use other positional sensors in addition to a single multi-axis linear accelerometer, such as rotational (gyroscopes or additional accelerometers) , magnetic, or locational sensors (GPS or altimeter) , or physiological sensors, such as heart-rate, muscular activity, ventilation
(breathing), or skin-temperature sensors, to improve the measurement of posture or activity. However, the addition of these other sensors increases the cost and complexity of the system, and makes it difficult to incorporate the system as software into a device such as a cell phone or wearable music player, which already incorporates a single multi-axis accelerometer but not these other sensor types . It would therefore be desirable to have improved apparatus and methods of monitoring and measuring posture and physical activity using multi-axis accelerometer sensors that would avoid the drawbacks of the above-described conventional systems.
In accordance with the present invention, an improved system and method of acquiring and/or processing data to monitor and measure posture and intensity of physical activity is disclosed. In one embodiment, an apparatus for simultaneously monitoring the posture and acceleration of an object includes a 3 -axis sensor, a memory, and at least one processor. The apparatus is configured to be attached to, held against, inserted into, or otherwise physically associated with the object to be monitored. The 3 -axis sensor is configured to sense apparent acceleration along each of a first axis, a second axis, and a third axis. The memory is operative to store data representing the sensed magnitudes of apparent acceleration along each of the three axes, and the processor is operative to process the data stored in the memory. Specifically, the processor uses a method to estimate both the actual gravity vector and the acceleration vector based upon the sensor signals. In this way, the apparatus provides two separate time-based streams of data that can be used individually or together to accurately determine the type and intensity of physical activity, and for any other suitable purpose.
In another embodiment, the acceleration data and the posture data are used to calculate one or more of the total EE, EEP/ and EEK of the system. In another embodiment, EEK in a plurality of directions are distinguished from each other.
In another embodiment, posture data, acceleration data, EE, EEK, EEK in one or more directions, EEP, and other inputs are used alone or in combination to evaluate activity patterns, determine activity types, or to evaluate movement efficiency, physical impairment or for other diagnostic, rehabilitative or training purposes .
In another embodiment, the accuracy of the estimation of posture and acceleration is improved by use of a feedback loop whose control signals are measures of quality derived from the posture and acceleration outputs .
In another embodiment, sensor errors are measured by analyzing the system outputs when the user is determined to be at rest. In another embodiment, the acceleration sensors are offset to correct the sensor errors, or to compensate for the effect of gravity, e.g., in order to make optimal use of the dynamic range of the A/D converter, or for any other suitable system design purposes .
In other embodiments, one or more of the above-described features may be implemented in whole or in part in software and incorporated into devices such as cell phones, personal digital assistants (PDAs) , personal music players, special purpose monitoring devices, implanted medical devices such as pacemakers, or any other suitable portable, wearable, or insertable device, and may make use of existing hardware or software features of the device to implement portions of the activity and posture monitoring system.
The various elements may be included in a single housing, or distributed among a plurality of housings, with suitable functional wired or wireless inter-housing communication,- or incorporated into articles to be worn such as apparel or clothing, or eyeglasses, sunglasses, wristwatches, or any other suitable accessories .
Fig. 2 is an illustration of certain force vectors and their relationships to the posture of a typical user; Fig. 3 is geometric model of the posture and acceleration signals that are sensed and calculated by the system of Fig. 1;
Fig. 4 is a schematic diagram of a method to discriminate between posture and acceleration in the system,- Fig. 5 is a flow diagram of a method of calibrating the input sensors automatically based on activity and posture patterns;
Fig. 6a is an illustration of a runner wearing an exemplary embodiment of a device according to this invention, showing certain axes of the runner and a reference orientation of the device;
Fig. 6b is a geometric representation of the magnitude of the runner's Kinetic Energy Expenditure (EEK) that is sensed and calculated by the system of Fig. 1, for various horizontal directions relative to the axes shown in Fig. 6a; and
The disclosures of U.S. Provisional Patent Application No. 60/930,749 filed May 18, 2007 entitled PHYSICAL ACTIVITY MONITOR, U.S. Provisional Patent Application No. 60/934,216 filed June 11, 2007 entitled PHYSICAL ACTIVITY AND POSTURE MONITOR, and U.S. Provisional Patent Application No. 60/936,117 filed June 18, 2007 entitled PHYSICAL ACTIVITY AND POSTURE MONITOR, are incorporated herein by reference in their entirety.
Activity type: "Activity type" is a qualitative description of physical activity, e.g., walking, running, swimming, vacuuming floors, washing dishes, climbing stairs, etc. Energy Expenditure: "Energy Expenditure (EE)" is a measure of the amount of energy that is expended during the physical effort required to perform an activity. For example, EE can be expressed in footsteps (an inaccurate measure) , calories of metabolic activity (a more accurate measure) , or any other suitable units of expended energy. For purposes of this disclosure, EE includes activity-related energy expenditure, but does not include basal or dietary metabolic activity. Posture: The "posture" of an object or person is the direction and magnitude of the tilt of the object or person in 3- dimensional space, relative to a reference direction in the horizontal plane. For example, as shown in Fig. 3, angles ΘG- 320 and φG< 322 associated with an exemplary apparent gravity vector G' 334 define the posture of an object or person aligned with the orientation of G' 334, relative to direction Y 304 in the XY plane. The posture of an individual is often meant to refer to the posture of the individual's spine, but can refer to the posture of the head or of an extremity. Activity monitoring: "Activity monitoring" can include capturing data indicative of either the activity type or the EE or both. One example of data indicative of activity type is upper body posture over time. Posture data can provide useful information in activity monitoring: e.g., whether an individual is recumbent or upright, during what time intervals, and at what times of day or night. Patterns of posture data can also be combined with patterns of measured EE for more complete activity monitoring.
Acceleration: "Acceleration" is the direction and magnitude of the acceleration of a user. References to a user may constitute references to a person or object being monitored.
Apparent gravity: "Apparent gravity" refers to the total measurement made by a sensor that is sensitive to the sum of acceleration and gravity, e.g., an accelerometer . Gravity vector or actual gravity vector: "Gravity vector" or "actual gravity vector" is a vector representing the gravity that is felt by an object, relative to a reference frame of the object. Since the gravity vector on earth always points towards the center of the earth, and is thus normal to the horizontal plane, changes in the direction of the gravity vector relative to a reference frame of the object result from changes in posture of the object. It is assumed in the embodiments described in this application that the object is at sea level on earth, and that the length of the gravity vector is "standard gravity", i.e., Ig « 9.8m/sec2. Corrections for other locations can be made using techniques known in the art.
Component of a vector: A "component" of a first vector refers to a second vector which can be added to a third vector to equal the first vector. In Fig. 3, vectors G 332 and μ 342 are components of vector G' 334.
Characteristic of a vector: A "characteristic" of a vector refers to the direction or magnitude of the vector, or to a statistic (such as the average value over a time interval) of the characteristic. In Fig. 3, ΘG- 320, φG- 322, and pG- 324 are characteristics of vector G' 334.
Potential Energy Expenditure (EEP) : "Potential Energy Expenditure" EEP is that portion of EE that is related to changes in vertical position that are resisted by gravity.
Kinetic Energy Expenditure (EEK) : "Kinetic Energy Expenditure" EEK is that portion of EE that is related to acceleration in any direction, which is resisted by the inertia of the object. Affected Mass (M) : "Affected mass" M is the mass of the body or body part of a person and objects to which the device is attached and to which the motion-related energy expenditure is applied. One example is the case of a single device used to measure the caloric expenditure of a runner wearing weights, for which M is the total mass of the runner and the weights. Another example is the use of multiple sensors on the trunk and extremities, to more accurately measure the energy expenditure in those different parts of the body, for which the total mass of the wearer is distributed among energy-expenditure formulas applied to the outputs of each of the respective sensors. For purposes of this disclosure, no distinction is made between mass and weight . Fig. 1 depicts an illustrative embodiment of a physical activity monitoring device 100, in accordance with the present invention. The physical activity monitoring device 100 provides proper operation and accuracy when physically associated in substantially any direction and orientation relative to the user. In addition, the device 100 provides visible, audible (e.g., alarms, tones, human speech) , and/or tactile feedback to the user.
In one embodiment, the physical activity monitor 100 includes an X-axis sensor 102, a Y-axis sensor 104, a Z-axis sensor 106, a signal multiplexer and analog-to-digital converter
(A/D) 110, a data processor 112, a program and data memory 120, a user display and controls 122, a real-time clock 124, an audio output converter 130, a speaker 132, one or more tactile vibration outputs 134, a PC/network interface 150, a wireless networking antenna 152, a wired networking connector 154, a battery charger connector 160, and a battery and power supply 162. The data processor 112 may selectively power down functional components during idle periods, such as powering down the sensors 102, 104, 106 between samples, to conserve power. When charging, the battery charger connector 160 is connected to a battery charger (not shown) . When in a network configuration, the
PC/network interface 150 is connected to a personal computer (PC; not shown) , or a local area network (LAN) , a wide area network
(WAN) , or any other suitable network, through either the wireless networking antenna 152 or the wired networking connector 154.
It is noted that alternative embodiments of the physical activity monitoring device 100 may include a subset of the functional components illustrated in Fig. 1. It is also noted that in the presently disclosed embodiment, the X, Y, and Z-axis sensors 102, 104, 106 are accelerometers oriented within the device 100 so that their respective sensitive axes, namely, the X-axis 302, the Y-axis 304, and the Z-axis 306 (see Fig. 3), are positioned at 90° relative to one another. In other embodiments, the X, Y, and Z-axis sensors may employ any other suitable technique for sensing gravity and acceleration, and may be oriented at any other suitable angle relative to one another. For example, each of the X, Y, and Z-axis sensors 102, 104, 106 may be a suspended proof-mass MEMS accelerometer such as the ADXL103 device sold by Analog Devices Inc., Norwood, Massachusetts, U.S.A. Alternatively, the X, Y, and Z-axis sensors 102, 104, 106 may be implemented using a single three-axis accelerometer such as the ADXL330 device sold by Analog Devices Inc. Alternatively, any other suitable set of accelerometers may be used. The outputs of the accelerometers are coupled to the multiplexer-A/D converter 110 to produce digitized accelerometer data for processing by the processor 112.
The X7 Y, and Z-axis sensors 102, 104, 106 sense gravity and acceleration along the X, Y, and Z axes 302, 304, 306, respectively, by measuring the projection of the apparent gravity vector G' 334 on their respective axes that is the sum of gravity G 332 at the location of the device 100 and the acceleration μ 342 (see Fig. 3) of the device 100.
The outputs of the X, Y, and Z-axis sensors 102, 104, 106 can be filtered in the analog or digital domain using any suitable hardware or software components to reject noise or signals at frequencies beyond those of interest for an intended application. For example, autogenous human trunk movement occurs at frequencies less than approximately 10 Hz. Filtering and sampling of outputs of the sensors 102, 104, 106 can be designed to ensure that relevant motion data can be captured without attenuation or aliasing. Additionally, data samples can be aggregated to reduce data storage requirements .
Fig. 3 is a geometric model 300 that illustrates the data format and operation of one embodiment of the physical activity monitoring device 100. Specifically, the geometric model 300, which is constructed with reference to both the spherical and Cartesian coordinate systems, may be used to quantify changes in the posture of the person or object to which the device is attached, and acceleration in the various directions relative to the horizontal plane. The geometric model 300 includes an entire unit sphere, of which only one octal portion is shown for clarity. The radius of the sphere is equal to 1 g, i.e., the magnitude of gravity on Earth. X-axis 302 represents the sensitive axis of the X-axis sensor 102, Y-axis 304 represents the sensitive axis of the Y-axis sensor 104, and Z-axis 306 represents the sensitive axis of the Z-axis sensor 106. In addition, Fig. 3 depicts one possible direction of an exemplary apparent gravity vector G' 334, where G' is the sum of the actual gravity vector G 332 and of the acceleration vector μ 342. Because G always points in the same direction, i.e., toward the center of the earth, the device 100 can determine changes in the posture and acceleration of the person or object to which it is attached by decomposing changes in the apparent gravity vector G' into its component parts, G and μ.
It is noted that the length of the vector G is equal to Ig, i.e., the radius of the unit sphere, and that G' and μ can be of any length and direction relative to G. It is also noted that use of the spherical coordinate system is for illustrative purposes only, and that all of the inputs, outputs, and formulas expressed herein can be represented and expressed in any other suitable coordinate system. Within the geometric model 300, direction (ΘG- 320 and φG- 322) and magnitude pG- 324 of the apparent gravity vector G' 334 can be determined by applying known trigonometric formulas to the outputs ax 312, ay 310, and az 314 of the X-axis, Y-axis, and Z- axis sensors sampled quasi-simultaneously. Methods of distinguishing the acceleration (μ) of device 100 from its posture (G) are described below with reference to Fig. 4. In a first variation on this method, feedback is not used and blocks 414 and 450 are not required. Signals representing the angles ΘG< (402 & 320) and cpG. (404 & 322) are filtered by filter components 410 and 412 respectively, to suppress just the effects of the acceleration vector μ 342 from the apparent gravity vector G' 334. Filter components 410 and 412 may be low-pass, band-pass, high-pass or combinations thereof, may be sampled or continuous, may use any suitable architecture such as finite- impulse response
(FIR) , infinite- impulse response (HR) , and frequency domain
(such as Discrete Cosine Transform (DCT) or Fast Fourier Transform (FFT) ) , and may use any suitable synthesized technique such as look-ahead, predictive, wavelet, or Kalman filtering. Further, any suitable filter architectures and values for the filter components 410 and 412 can be chosen to perform this step if there are known differences between the characteristic variations of G and μ. For example, if the majority of the magnitude variation of the acceleration vector μ is in a sufficiently higher frequency band than the directional variation of the posture, i.e., of the actual gravity vector G, then a 15- tap, 0.5 Hz FIR low-pass filter with a 6 Hz sampling rate will attenuate a periodic acceleration vector μ with a period of 1 second by 33 dB, while attenuating a periodic G with a period of 9 seconds by less than 1 dB .
The previous example improves on those conventional systems that filter out the G signal from the μ signal, by instead filtering out the effect of μ on the direction of G, imposing a Ig length to G, and determining μ arithmetically (see Fig. 4) . In many applications, G is considerably stronger than μ and varies in direction over a lower frequency band and range, and thus the presently disclosed embodiment provides a more accurate output. Additionally, the constant and known length of G (pG = Ig) is used in the presently disclosed embodiment to further refine the measurement process, and for other purposes such as calibrating the sensors during use . As shown in Fig. 4, the filters 410 and 412 generate output signals corresponding to angles ΘG and φG, respectively, which define the actual gravity vector G. Next, assuming that PG = Ig, signals representing the angles ΘG and φG and the known value Ig of the length pG of G are converted to Cartesian coordinates by the converter 420. Similarly, the signals representing the angles ΘG< 402 and φG' 404, and the length pG- 406, are converted to Cartesian coordinates by the converter 422. Finally, a representation of the acceleration vector μ 440 is obtained at a summation node 430 by subtracting the actual gravity vector G from the apparent gravity vector G' .
In a variation on this method, a filter control signal 451 is used to optimize the filters 410 and 412, with the optional addition of a filter 414 and a feedback function 450. In this variation, the characteristics of the filters 410, 412 and 414 are adjusted by the filter control signal 451 to minimize deviations from target values of characteristics of G and μ. For example, target values in some applications could be pG = 1 and avg(pμ) = 0, since the magnitude of the gravity vector is a constant, Ig, and the average acceleration of a human body performing physical activity is zero if its velocity tends towards a constant value. In such applications, an example of an optimization rule could be to reduce the bandwidth of low-pass filters 410, 412, and 414 if pG ≠ 1 (feed-through of μ into G), and increase the bandwidth if avg(pμ) ≠ 0 (feed-through of G into μ) . In other variations on this method, components of G or μ, such as the vertical component of μ, can be optimized.
It is noted that in some applications, a delay can occur between data acquisition and required output to the user. This can allow a range of recursive, successive-approximation or look- ahead filters, or signal synthesis techniques, to be used in variations of this method. For example, in cases where the frequency band of variations in G and μ overlap, the filters 410, 412, and 414 can be implemented as signal synthesizers, and synthesize an approximation of G with a first constraint of PG = 1. The approximate G sequence can then be optimized recursively until average (pμ) = 0. An additional exemplary quality metric is that μ often oscillates along a stable direction and its inverse over several cycles when the user is active in physical exercise.
A method to calibrate the accuracy of the X, Y, and Z sensors and the signal path through the multiplexer-A/D converter 110 is described below, with reference to Figs. 1, 3, and 5. For an individual engaged in typical activities, immobility of the unit for more than several seconds is significant and indicative that the unit has been removed or that the user is in a resting
(not necessarily supine) position, and thus G is equal to Ig and μ is equal to 0. These periods can be determined by the lack of significant change in values of the direction (ΘG- and φG- ) and length (pc) of G' . During these periods, the gain and/or offset of the X, Y, and Z-axis sensors 102, 104, and 106 can be adjusted to zero out the error in the computed value of pG. Any suitable calibration rule(s) can be applied based on sensor characteristics and the initial factory calibration of the system. For example, the ADXL330 (Analog Devices, Inc.) exhibits almost no gain error or drift over time or temperature when referenced to the power supply. The offsets of the three channels, however, can vary and drift significantly, but the drift slopes can be characterized in an initial calibration process and a weighted correction can be applied to each axis during each field-calibration opportunity. During periods of immobility, the accelerometer sample rate can also be reduced to conserve battery power. Fig. 5 is a flow diagram depicting the operation of a system, in accordance with the present invention, which recalibrates the accelerometers when the system is at rest (e.g., when the user is immobile or has removed the system) . This ensures that the accuracy of the accelerometers and A/D is maximized, and also serves as a periodic diagnostic on the proper operation of the system. As depicted in step 510, the system is monitored for occurrence of a significant period of time (e.g., several seconds) during which sensor values are constant, i.e., substantially unchanging, signaling that the system is at rest or is subject to a constant linear velocity (μ « 0) . Specifically, the multi-axis sensor measures the apparent gravity vector G' , during a specified time interval, to obtain a plurality of substantially unchanging samples of G' . Next, as depicted in step 520, using the samples of G', the magnitude pG. of G' is measured. The sensor calibration is then adjusted to set pG- to Ig.
In a variation of this calibration method, the user can be instructed to place the device in one or more reference positions, or to manipulate the device through a sequence of movements and/or rotations, and the values of calibration adjustments can be determined by monitoring changes in values of the direction (ΘG> and φG. ) and/or length (pG<) of G' .
In another variation of this calibration method, the offsets of each of the three signals from the three accelerometers can be adjusted to better use the full range of A/D converter 110. In three parallel feedback loops, the X, Y, and Z characteristics of the calculated value of G are used by offset feedback 170 to generate respectively three offsets to each of said signals, so that the input signals to A/D converter 110 are at the optimal point within the scale of A/D converter 110 when μ 342 = 0, for any orientation of G. The optimal point may be the mid-scale point of the A/D converter 110 for all three axes. Alternatively, any other suitable point (s) may be chosen, for example, the mid- scale point for the most horizontal axis, a -0.5g offset for the middle axis, and a -Ig offset for the most vertical axis, since the maximum vertical acceleration with gravity is likely to be a free fall, but the maximum vertical acceleration against gravity- can be several g's.
In an alternative embodiment, the EE of the system and the various components of EE can be calculated. Fig. 3 shows vector μv 344, which is the vertical component of the acceleration vector μ
342. μv 344 is the projection of μ 342 on G 332, and can be determined from μ and G using known trigonometric identities.
The magnitude of vector μv 344 (pμV) can be used, together with the affected mass M, the gravity g, and the sampling interval t, to calculate the change in Potential Energy Expenditure (EEP) . In one variant of this method, a formula derived from the Newtonian work equation W = F • d can be used, assuming that pμV is constant over the sampling interval t, as follows :
EEp = X • M • g • pμV • t2 (1)
Similarly, the magnitude of vector μ 342 (pμ) can be used, together with the affected mass M and the sampling interval t, to calculate the change in Kinetic Energy Expenditure (EEK) . In one variant of this method, a formula derived from the Newtonian work equation W = F • d can be used, assuming that pμ is constant over the sampling interval t, as follows:
EEK = M • M • pμ 2 • t2 (2)
It can be observed that the EEP formula (1) above is linear with pμv/ whereas the EEK formula (2) is quadratic with pμ . It can also be observed that μv contributes to both EEP and EEK. In this and other variants of the method, separating out the estimation of EEp from EEK allows the different relationships between acceleration and EE to be derived and used for the various components of EE. The Kinetic Energy Expenditure (EEK) can be calculated independently for horizontal and vertical directions, using known trigonometric identities to calculate μH (where μH is the projection of μ on the horizontal plane normal to G), and by substituting μv and μH respectively for μ in Formula (2) above to calculate the horizontal and vertical components of EEK.
In certain applications, the horizontal component of EEK can be further decomposed into component values corresponding to the horizontal projections of the anterior-posterior and medio- lateral axes of the user. This can be done, for example, by constraining the attachment of the device in known alignment with the anterior-posterior or medio- lateral axes of the user.
Fig. 6a is an illustration of a runner wearing an exemplary embodiment of a device 605 according to the present invention. As shown in Fig. 6a, a reference orientation 610 is associated with the device 605. The reference orientation 610 can be aligned with one of the sensed acceleration axes 302, 304, 306 (see Fig. 3) of the device, or can be any other orientation that is defined with respect to the axes 302, 304, 306. The device 605 and its orientation 610 are independent of the axes of the runner, and thus the device 605 can be attached in any position and any orientation. In a preferred embodiment, the device is attached to the belt of the runner. Fig 6a also shows the projection of the runner's anterior-posterior axis 620, 630 and medio- lateral axis 640, 650 on the horizontal plane.
Fig. 6b is a geometric model 600 that illustrates a method, in accordance with the present invention, to distinguish the anterior-posterior axis 620, 630 from the medio-lateral axis 640, 650 in applications where the horizontal component of EEK peaks in one or more known directions, such as in the anterior-posterior direction for rowing or running, and if the attachment of the device 605 (see Fig. 6a) is not in known alignment with the anterior-posterior or medio-lateral axes of the user. Geometric model 600 is a plane, intersecting the origin of unit sphere 300 (see Fig. 3) , and normal to the vector G 332 (see Fig. 3) . Plane 600 is thus the horizontal plane.
In the embodiment shown in Fig. 6b, plane 600 has been divided into 16 equal sectors such as sector 670. In other embodiments, a larger or smaller number of sectors can be used. The sector boundaries are defined with respect to a reference direction 615, which is the projection of the reference orientation 610 of the device 605 on the horizontal plane 600. Reference orientation 610 can be arbitrarily chosen among available reference orientations of device 605 (such as any axis of device 605) to minimize the sensitivity of the reference direction 615 to changes in the gravity vector G 332, for example, by choosing a reference orientation 610 which is close to horizontal. During a period of physical activity, the accumulated horizontal component of EEK can be decomposed into sixteen horizontal directions, relative to reference direction 615, using known trigonometric identities. As shown in Fig. 6b, the magnitude of EEK in the direction of sector 670 is represented by the distance between the origin and line 660. The magnitude of EEK in the other fifteen directions is shown in a similar manner. In the example shown in Fig. 6b, EEK peaks approximately along an axis running inside sectors 680 and 690, with a secondary peak approximately at right angles to sectors 680 and 690, i.e., sectors 640 and 650. This pattern of a major peak in one direction and a secondary peak at right angles to the major peak is typical of the horizontal kinetic energy expenditure of a runner. If the physical activity is known a priori to be running, it is highly likely that components 680 and 690 are aligned with the anterior-posterior axis 620, 630 of the user, and thus EEK can be estimated for every direction relative to the user, without a- priori knowledge of the alignment between device 605 and the axes of the user. These formulas (1) , (2) for calculating EEP, EEK, the components of EEK/ and the total EE may be modified to compensate for the effects of human physiology, complex motion, affected mass, extremity motion, device attachment, friction, particular types of physical activity, or for other factors. For example, potential energy losses and kinetic energy gains can be ignored in environments where no autogenous energy is expended during periods of downward vertical motion in quasi-free fall. Additionally, a sample of individuals from a target population can be tested in a controlled environment, and the data can then be correlated to independent and simultaneous measures of actual energy expenditure, such as from oxygen consumption or doubly- labeled water. Using known statistical techniques, offset and scaling factors for the formulas (1) , (2) that yield a best-fit correlation with the set of measured energy data points can be determined. The offset and scaling factors can then be used to modify the formulas (1) , (2) and to provide improved measurement accuracy in subsequent use. Offset and scaling factors can be derived for the total EE or for any of its components or sub- components. Different offset and scaling factors for the various components can compensate for the difference in the effect of the various factors enumerated above on the different components of EE.
These formulas (1) , (2) may also be modified by including posture data, posture patterns, or activity-type information, as additional variables in the correlation.
Values or patterns of the relationships between EE or its components, posture, and time can be used to analyze activity patterns, movement efficiency and impairment. One example is running efficiency. If a runner's anterior-posterior horizontal EEK contributes to moving her forward, and that runner's EEP and vertical EEx contribute primarily to up-and-down motion, then ratios of these energy expenditure components or any other suitable energy expenditure components may be related to running efficiency and may be used as immediate or time-delayed training feedback.
Fig. 7 is a flow diagram depicting the operation of a system, in accordance with the present invention, which determines that the user is subjected to an externally-applied acceleration, such as is caused by traveling in an elevator, automobile, airplane, or other means of transportation. As depicted in step 710, the system is monitored for occurrence of a significant period of time (e.g., several seconds) during which μv remains greater than a low threshold (i.e., persistent acceleration against gravity) . Under such conditions, it is highly likely that the user is being lifted upwards by an elevator or airplane. Similarly, the system is monitored for occurrence of a significant period of time (e.g., several seconds) during which μH in any one direction remains greater than a low threshold (i.e., persistent acceleration in any one horizontal direction) . Under such conditions, it is highly likely that the user is being moved horizontally in an automobile or other form of transportation. Next, as depicted in step 720, the system flags the event or alerts the user that the measurement may be compromised by externally-applied acceleration. Finally, as depicted in step 730, the system accepts a user input to confirm the presence of the externally-caused acceleration, when the transportation is ended, or in the case of a false alarm. In certain applications, these systems can be applied using an accelerometer with two sensing axes, if the wearer's motion is similarly limited to two directions. A runner, for example, moves her trunk in the vertical and anterior-posterior directions, but her trunk motion in the medio- lateral direction is limited. For such an application, a two-axis sensor whose sensor plane is aligned with the runner's sagittal plane can be used. The sagittal plane is a plane that divides the body into left and right portions.
In alternative embodiments, other sensing devices and/or techniques may be used along with or in place of accelerometers .
Systems that include these methods can also include one-way or two-way communication features. One example of a use of such a system is to remotely monitor an elderly person living by herself, to use posture or EE data to determine changes in health, mobility or lifestyle, and to use the communication features to communicate remotely with the individual .
1. An apparatus for monitoring posture and acceleration of a user, said apparatus being mountable to said user, comprising: a multi-axis sensor operative to measure an apparent gravity vector, the measured apparent gravity vector having a first component corresponding to an actual gravity vector and a second component corresponding to an acceleration vector,- at least one memory operative to store data representative of the measured apparent gravity vector; and at least one processor operative to process the data stored in said at least one memory, said at least one processor being operative : to calculate a representation of the measured apparent gravity vector; to calculate a representation of said actual gravity vector, said actual gravity vector having a predetermined magnitude equal to 1 g; to subtract said representation of said actual gravity vector from said representation of the measured apparent gravity vector to obtain a representation of said acceleration vector; and to generate at least one indication of one or more of the posture of said user based at least upon a direction of said representation of said actual gravity vector, and the acceleration of said user based at least upon said representation of said acceleration vector.
7. A method of monitoring posture and acceleration of a user, comprising the steps of: measuring, by a multi-axis sensor mounted to said user, an apparent gravity vector, the measured apparent gravity vector having a first component corresponding to an actual gravity vector and a second component corresponding to an acceleration vector, said actual gravity vector having a predetermined magnitude equal to 1 g; subtracting said actual gravity vector from the measured apparent gravity vector to obtain said acceleration vector; and generating at least one indication of one or more of the posture of said user based at least upon a direction of said actual gravity vector, and the acceleration of said user based at least upon said acceleration vector.
12. A method of monitoring energy expended by a user, comprising the steps of: measuring an acceleration vector representing an acceleration of said user, said acceleration vector having a vertical component; performing one or more of a first calculating step and a second calculating step, wherein said first calculating step includes calculating a first amount of energy expended by said user based at least upon the measured acceleration vector, said first amount of energy corresponding to a change in kinetic energy of said user, and wherein said second calculating step includes calculating a second amount of energy expended by said user based at least upon said vertical component of said acceleration vector, said second amount of energy expended corresponding to a change in potential energy of said user; and generating an indication of an amount of total energy expended by said user based upon one or more of said calculated first and second amounts of energy.
14. The method of claim 12 wherein said acceleration vector has at least one horizontal component, and wherein said first calculating step includes: calculating said first amount of energy expended by said user based at least upon one or more of said vertical component and said at least one horizontal component of said acceleration vector, said first amount of energy corresponding to at least one change in kinetic energy of said user in each of one or more of a vertical direction and at least one horizontal direction.
15. The method of claim 14 further including determining at least one characteristic of motion of said user based at least in part upon one or more of said at least one change in kinetic energy of said user in each of one or more of said vertical direction and said at least one horizontal direction and said change in potential energy of said user over a specified time interval .
16. A method of calibrating an accelerometer , said accelerometer including a multi-axis sensor, said method comprising the steps of: measuring, by said multi-axis sensor during a specified time interval, an apparent gravity vector, said measured apparent gravity vector having a magnitude and a direction, thereby obtaining a plurality of samples of said apparent gravity vector during said specified time interval, said multi-axis sensor having at least one associated axis, wherein said at least one axis of said multi-axis sensor has at least one associated offset and at least one associated gain, one or more of the at least one offset and the at least one gain of the at least one axis of said multi-axis sensor being proportional to at least one adjustable parameter; determining, from said plurality of samples of said apparent gravity vector, whether at least one characteristic of the measured apparent gravity vector is substantially unchanging during said specified time interval; and in the event said at least one characteristic of the measured apparent gravity vector is substantially unchanging during said specified time interval, adjusting said at least one adjustable parameter to adjust one or more of the at least one offset and the at least one gain of the at least one axis of said multi-axis sensor, thereby setting one or more of the magnitude of the measured apparent gravity vector to a first specified value and the direction of the measured apparent gravity vector to a second specified value.
18. The method of claim 16 further including, prior to said measuring of said apparent gravity vector, placing said multi- axis sensor in a predetermined reference position.
20. A method of monitoring an externally applied acceleration of a user, comprising the steps of: measuring an acceleration vector representing an acceleration of said user, the measured acceleration vector having a vertical component in a vertical direction, and at least one horizontal component in at least one horizontal direction; monitoring, during a specified time interval, one or more of a magnitude of said vertical component and at least one magnitude of said at least one horizontal component of the measured acceleration vector; and in the event one or more of the magnitude of said vertical component and the at least one magnitude of said at least one horizontal component of the measured acceleration vector remain in excess of at least one specified threshold value during said specified time interval, generating at least one signal to indicate a presence of said externally applied acceleration.
22. A method of calibrating an accelerometer, said accelerometer including a multi-axis sensor, said method comprising the steps of: measuring, by said multi-axis sensor, an acceleration vector with reference to a plurality of axes to provide a plurality of analog sensor signals, each of said plurality of analog sensor signals representing a projection of said acceleration vector on a respective one of said plurality of axes, each of said plurality of analog sensor signals having an associated offset proportional to an adjustable parameter; for at least one of the respective analog sensor signals, adjusting said associated adjustable parameter to adjust the offset of the respective analog sensor signal; and converting each of the respective analog sensor signals to a corresponding digital sensor signal, wherein the digital sensor signal corresponding to at least one of the respective analog sensor signals has a dynamic range determined at least in part by the adjusted offset of the respective analog sensor signal.
24. A method of monitoring energy expended by a user during a specified activity, comprising the steps of: measuring, by a sensor mounted to said user, an acceleration vector representing an acceleration of said user, the measured acceleration vector having a plurality of horizontal components; determining a reference direction relative to said sensor, said reference direction lying in a horizontal plane; and calculating a plurality of amounts of energy expended by said user based at least upon said plurality of horizontal components of the measured acceleration vector, respectively, each of said plurality of amounts of energy corresponding to a sum of changes in kinetic energy of said user in a respective one of a plurality of horizontal directions lying in said horizontal plane, each of said plurality of horizontal directions being defined relative to said reference direction of said sensor.
25. The method of claim 24 wherein said sensor has an associated axis defining a reference orientation of said sensor, and wherein said determining of said reference direction of said sensor includes determining a direction of a horizontal projection of said axis of said sensor on said horizontal plane.
27. The method of claim 25 further including: mounting said sensor to said user in such a manner that a value of an alignment offset between said reference orientation of said sensor and at least one axis of said user is known,- and wherein said calculating of said plurality of amounts of energy expended by said user includes calculating, using the known value of said alignment offset, said plurality of amounts of energy expended by said user in a plurality of horizontal directions defined relative to said at least one axis of said user.
28. The method of claim 24: wherein said plurality of amounts of energy expended by said user in said plurality of horizontal directions defines a pattern of energy expenditure for said specified activity, wherein said method further includes determining, based upon the defined pattern of energy expenditure for said specified activity, an approximate direction of a horizontal projection of at least one axis of said user on said horizontal plane; and wherein said calculating of said plurality of amounts of energy expended by said user includes calculating said plurality of amounts of energy expended by said user in a plurality of horizontal directions defined relative to said at least one axis of said user.
29. The method of claim 28 further including determining an alignment offset between said approximate direction of said horizontal projection of said at least one axis of said user and said reference direction of said sensor.
PCT/US2008/003745 2007-05-18 2008-03-21 Newtonian physical activity monitor WO2008143738A1 (en)
US93074907P true 2007-05-18 2007-05-18
US60/930,749 2007-05-18
US93421607P true 2007-06-11 2007-06-11
US60/934,216 2007-06-11
US93611707P true 2007-06-18 2007-06-18
US60/936,117 2007-06-18
WO2008143738A1 true WO2008143738A1 (en) 2008-11-27
WO2008143738A8 WO2008143738A8 (en) 2009-04-09
PCT/US2008/003745 WO2008143738A1 (en) 2007-05-18 2008-03-21 Newtonian physical activity monitor
WO (1) WO2008143738A1 (en)
WO2017005130A1 (en) * 2015-07-09 2017-01-12 深圳市声禾科技有限公司 Method and device for measuring energy used by human body during exercise, and pedometer
EP2536334B1 (en) * 2010-02-15 2013-12-25 Cardiac Pacemakers, Inc. Methods for constructing posture calibration matrices
JP5515875B2 (en) * 2010-03-08 2014-06-11 セイコーエプソン株式会社 Fall detection device, fall detection method
ES2394842B1 (en) * 2010-10-22 2013-11-26 Universidad De Sevilla Portable and adaptive multimodal monitor for humans based on biomechanical-physiological avatar for the detection of physical risk events.
US9451916B2 (en) 2011-08-02 2016-09-27 Sway Medical Llc System and method for assessing postural sway and human motion
CA2868009A1 (en) * 2012-01-18 2013-07-25 Nike Innovate C.V. Activity and inactivity monitoring
JP5978687B2 (en) * 2012-03-22 2016-08-24 富士ゼロックス株式会社 Moving situation determination program and moving situation determination apparatus
CN104519787B (en) * 2012-06-04 2017-03-22 耐克创新有限合伙公司 Fitness training system having a plurality of sensor inputs using an energy consumption calculation
KR101774120B1 (en) 2012-06-05 2017-09-01 나이키 이노베이트 씨.브이. Multi-activity platform and interface
JP6273364B2 (en) * 2013-08-23 2018-01-31 ナイキ イノベイト シーブイ Energy consuming equipment
US9608087B2 (en) * 2015-04-08 2017-03-28 Globalfoundries Inc. Integrated circuits with spacer chamfering and methods of spacer chamfering
EP3206654A4 (en) * 2014-10-17 2018-06-06 Stryker Corporation Person support apparatuses with motion monitoring
DE102015219572A1 (en) * 2015-10-09 2017-04-13 Sivantos Pte. Ltd. Method for operating a hearing device and hearing device
EP3364150A4 (en) * 2015-10-14 2018-09-26 Alps Electric Co., Ltd. Wearable device, method for measuring orientation of same, and program
DE69628407T2 (en) 1995-09-12 2004-05-06 Omron Corp. Pedometer
CN1256752A (en) 1998-02-25 2000-06-14 皇家菲利浦电子有限公司 Method of and system for measuring performance during exercise activity, and athletic shoe for use in the system
WO2002102475A1 (en) 2001-06-07 2002-12-27 Rutgers, The State University Of New Jersey Method and apparatus for analyzing a golf stroke
US6816326B2 (en) 2002-07-12 2004-11-09 Schott Glas Optical system with compensated spatial dispersion
EP1571990B1 (en) 2002-12-10 2008-06-04 Philips Electronics N.V. Activity monitoring
WO2007008930A2 (en) 2005-07-13 2007-01-18 Ultimate Balance, Inc. Orientation and motion sensing in athletic training systems, physical rehabilitation and evaluation systems, and hand-held devices
2008-03-21 US US12/077,955 patent/US7634379B2/en active Active
2008-03-21 WO PCT/US2008/003745 patent/WO2008143738A1/en active Application Filing
US20080288200A1 (en) 2008-11-20
WO2008143738A8 (en) 2009-04-09
US7634379B2 (en) 2009-12-15
US9186096B2 (en) 2015-11-17 System and method for measuring balance and track motion in mammals
EP2387357B1 (en) 2016-04-20 Gait monitor
CN101969848B (en) 2012-11-28 An activity monitoring system insensitive to accelerations induced by external motion factors
Marschollek et al. 2008 A performance comparison of accelerometry-based step detection algorithms on a large, non-laboratory sample of healthy and mobility-impaired persons
Ref document number: 08727062