Patent Application
20250050170
Provided are systems and methods for measurement of maximal oxygen consumption of a subject.
Cardiorespiratory fitness represents an important aspect of a subject's physical condition and is a key indicator of a medical condition of the subject. The measurement of cardiorespiratory fitness is useful in many different situations, such as improving performance in elite athletes, monitoring a cardiac patient undergoing a rehabilitation program, monitoring cardiovascular system of deconditioning astronauts while in space, etc. The gold standard measurement for cardiorespiratory fitness assessment is the maximal oxygen uptake rate one can withstand, and more commonly known as the VO2 max measurement. Obtaining an accurate VO2 max measurement requires a maximal effort ergospirometry test. The test involves complex gas exchangers and heavy machinery, requires expertise and so is costly, and not always possible in all environments (e.g., sports centres, developing countries, local clinics, space vehicles).
An aim to provide an affordable, light-weight, and operationally more efficient solution to the problem of cardiorespiratory fitness measurement, that offers a high level of measurement accuracy.
Provided herein is a system (100) for determining a cardiorespiratory fitness measurement of a subject (50), comprising a processing unit (160), the processing unit (160) configured to:
The cardiorespiratory fitness measurement of the subject may be determined from the SPMKE by:
The transformation of the SPMKE into the LVO2F, may comprise a step of converting each PMKE of the SPMKE into a scalar value representative of the non-exercise period (N), wherein the scalar value is an order statistic.
The conversion of each PMKE of the SPMKE into a scalar value may be performed:
The same order statistic may be determined for each non-exercise period (N), at least some of these order statistics (OS) form a set of order statistics, SOS, for the exercise session (E), and the SOS is transformed into the LVO2F of the subject.
The transformation of the SOS into the LVO2F may be performed:
The order statistic (OS) may be either a median or a quantile, and the order statistics (OS) within each set of order statistics (SOS) are either all medians or all quantiles.
The exercise session (E) may bring the subject to a level of exhaustion, preferably to a level of full exhaustion.
The time window (W) may be for a period of w seconds, and the value of w is greater than a duration of 1 cardiac cycle of the subject.
The value of w may be the same during each non-exercise period (N) and during the exercise session (E).
The cardiorespiratory fitness measurement of the subject (50) may be a VO2 max measurement of the subject (50).
One of the MSM (111) may be a seismocardiograph (114), SCG, and optionally a further MSM (111) may be a ballistocardiograph (112), BCG.
Further provided is a computer-implemented method for determining a cardiorespiratory fitness measurement of a subject (50), comprising:
The computer-implemented method may comprise a limitation as defined herein, in particular above.
Further provided is a computer program or computer program product having instructions which when executed by a computing device or system cause the computing device or system to perform the computer-implemented method as defined herein.
Before the present system and method of the invention are described, it is to be understood that this invention is not limited to particular systems and methods or combinations described, since such systems and methods and combinations may, of course, vary. It is also to be understood that the terminology used herein is not intended to be limiting, since the scope of the present invention will be limited only by the appended claims.
As used herein, the singular forms “a”, “an”, and “the” include both singular and plural referents unless the context clearly dictates otherwise.
The terms “comprising”, “comprises” and “comprised of” as used herein are synonymous with “including”, “includes” or “containing”, “contains”, and are inclusive or open-ended and do not exclude additional, non-recited members, elements or method steps. It will be appreciated that the terms “comprising”, “comprises” and “comprised of” as used herein comprise the terms “consisting of”, “consists” and “consists of”.
The recitation of numerical ranges by endpoints includes all numbers and fractions subsumed within the respective ranges, as well as the recited endpoints.
The term “about” or “approximately” as used herein when referring to a measurable value such as a parameter, an amount, a temporal duration, and the like, is meant to encompass variations of +/−10% or less, preferably +/−5% or less, more preferably +/−1% or less, and still more preferably +/−0.1% or less of and from the specified value, insofar such variations are appropriate to perform in the disclosed invention. It is to be understood that the value to which the modifier “about” or “approximately” refers is itself also specifically, and preferably, disclosed.
Whereas the terms “one or more” or “at least one”, such as one or more or at least one member(s) of a group of members, is clear per se, by means of further exemplification, the term encompasses inter alia a reference to any one of said members, or to any two or more of said members, such as, e.g., any ≥3, ≥4, ≥5, ≥6 or ≥7 etc. of said members, and up to all said members.
The term “each” means an individual value or group of values within a set. It does not necessarily mean each and every value within the set.
All references cited in the present specification are hereby incorporated by reference in their entirety. In particular, the teachings of all references herein specifically referred to are incorporated by reference.
Unless otherwise defined, all terms used in disclosing the invention, including technical and scientific terms, have the meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. By means of further guidance, term definitions are included to better appreciate the teaching of the present invention.
In the following passages, different aspects of the invention are defined in more detail. Each aspect so defined may be combined with any other aspect or aspects unless clearly indicated to the contrary. In particular, any feature indicated as being preferred or advantageous may be combined with any other feature or features indicated as being preferred or advantageous.
Reference throughout this specification to “one embodiment” or “an embodiment” means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, appearances of the phrases “in one embodiment” or “in an embodiment” in various places throughout this specification are not necessarily all referring to the same embodiment, but may. Furthermore, the particular features, structures or characteristics may be combined in any suitable manner, as would be apparent to a person skilled in the art from this disclosure, in one or more embodiments. Furthermore, while some embodiments described herein include some but not other features included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention, and form different embodiments, as would be understood by those in the art. For example, in the appended claims, any of the claimed embodiments can be used in any combination.
In the present description of the invention, reference is made to the accompanying drawings that form a part hereof, and in which are shown by way of illustration only of specific embodiments in which the invention may be practiced. Parenthesized or emboldened reference numerals affixed to respective elements merely exemplify the elements by way of example, with which it is not intended to limit the respective elements. Unless otherwise indicated, all figures and drawings in this document are not to scale and are chosen for the purpose of illustrating different embodiments of the invention. In particular the dimensions of the various components are depicted in illustrative terms only, and no relationship between the dimensions of the various components should be inferred from the drawings, unless so indicated.
It is to be understood that other embodiments may be utilised and structural or logical changes may be made without departing from the scope of the present invention. The following detailed description, therefore, is not to be taken in a limiting sense, and the scope of the present invention is defined by the appended claims.
Provided is a system (100) and method for determining a cardiorespiratory fitness measurement, in particular, a maximal oxygen consumption, VO2 max, of a subject (50). An exemplary system (100) is shown in
The method or system comprising a processing unit (160) configured to receive a stream of signals (120) outputted by a sensor unit (110). The sensor unit (110) comprises one or more motion sensing modules (MSMs) (111), each motion sensing module (MSM) being configured for detection of bodily vibrations induced by cardiovascular activity within the body. The sensor unit (110) may comprise a motion sensing module (MSM) (111) that is a seismocardiograph (SCG) (114), and optionally a motion sensing module (MSM) (111) that is a ballistocardiograph (BCG) (112), and optionally one or more electrocardiograph, ECG, electrodes.
The method or processing unit (160) determines, from the signal stream (120), a kinetic energy data stream, KEDS, for each non-exercise period (N) of an exercise session (E) for gradually increasing the oxygen consumption, VO2, of the subject, wherein the exercise session (E) comprises a plurality of alternating exertion segments (S) and non-exercise periods (N).
The method or processing unit (160) is further configured to, for each non-exercise period (N):
The method or processing unit (160) is further configured to obtain a set of PMKEs (SPMKE), wherein a SPMKE contains a plurality of PMKEs (periods of mean kinetic energy) obtained during the course of the exercise session (E). The cardiorespiratory fitness measurement of the subject can be determined from the set of PMKEs (SPMKE), as the SPMKE is indicative of the oxygen consumption rate (VO2) over time for the exercise session (E). The steps are shown in
The cardiorespiratory fitness measurement of a subject is a measure of a subject's physical condition. In particular, it is a measure of the maximal rate of oxygen consumption rate (by volume) under maximal exertion by the subject. In other words, the more oxygen the subject is able to consume, the greater amount of oxygen arrives in the muscle cells, allowing the subject to withstand longer and/or more intensive exercise. Maximal oxygen uptake can be detected if an increase in power (for example the resistance on a bike exercise) is not met by an increase in oxygen consumption rate. This means that the body is no longer capable of increasing its oxygen consumption and enters the anaerobic phase (energy consumption in the absence of oxygen). VO2 max is one measurement of cardiorespiratory fitness.
An exercise session (E) undertaken by the subject contains a plurality of alternating exertion segments (S) and non-exercise periods (N). The exercise session (E) gradually increases the VO2 of the subject. An example of an exercise session (E) is shown in
The exertion segment (S) may involve using a stationary exercise apparatus such as a bicycle, elliptical trainer, treadmill.
The intensity of the exertion segment (S) is set to increase the heart and respiration rates. The intensity of the exertion segment (S) may be between 50 and 300 W, i.e., the subject expends this amount of energy during the exertion segment (S). The intensity of the exertion segment (S) may be the same for each exertion segment (S) throughout the exercise session (E). The intensity of an exertion segment (S) may gradually increase as the exercise session (E) progresses.
The length of an exertion segment (S) may be 30 to 180 seconds, more preferably 60 to 180 seconds. The length of an exertion segment (S) may gradually increase as the exercise session (E) progresses. Preferably, the length of the exertion segment (S) is the same for each exertion segment (S) throughout the exercise session (E).
The first exertion segment (S) may be a warm-up session designed to gently raise the heart and respiration rates. The intensity of the warm-up session may be lower than in subsequent exertion segments (S). The length of the warm-up session may be longer than in subsequent exertion segments (S).
The number of exertion segments (S) within the exercise session (E) may vary from subject to subject and, as mentioned elsewhere, is determined by a level of exhaustion of the subject. The subject may continue to exercise until they could no longer continue (full exhaustion), or until they reach a pre-determined level of exhaustion, e.g., 16 to 19 on the Borg scale (rating of perceived exertion (RPE)).
The non-exercise period (N) is a period of non-activity. During the non-exercise period (N) the subject may sit, stand up, or lie down. The non-exercise period (N) typically follows an exertion segment (S); in such case it is a rest or recovery period.
The length of a non-exercise period (N) may be 20 to 60 seconds, more preferably 30 to 45 seconds. Preferably, the length of the non-exercise period (N) is the same for each non-exercise period (N) throughout the exercise session (E).
The period prior to the first exertion segment (S) is a baseline period in which a baseline measurement is taken of the subject; this may or may not be used as non-exercise period (N) for determining the cardiorespiratory fitness, in particular, VO2 max.
The period after the first exertion segment (warm-up session) may or may not be used as non-exercise period (N) for determining the cardiorespiratory fitness measurement, in particular, the VO2 max. Preferably it is used as non-exercise period (N) for determining the cardiorespiratory fitness measurement, in particular, the VO2 max.
The sensor unit (110) comprises one or more motion sensing modules (MSMs) (111). It may further optionally comprise one or more electrocardiograph (ECG) electrodes.
The sensor unit (110) outputs a signal stream (120) from one of more movement sensors in the one or more motion sensing modules (MSMs) (111). The ECG output may or may not be part of the signal stream (120).
A motion sensing module (MSM) (111) detects bodily vibrations induced by cardiovascular activity within the body. The bodily vibrations originate from the heart, and are transmitted through the body via tissue structures, including the vasculature and the skeleton. The bodily vibrations hence demonstrate a certain pattern that repeats with each heartbeat. The bodily vibrations are detected by mechanical transference to the motion sensing module (MSM) (111). A motion sensing module (MSM) (111) is discrete. By discrete, it is meant that it occupies or detects movements from a discrete region of the body. Where there are multiple motion sensing modules (MSMs), they may be disposed at different discrete locations of the body. Typically, a motion sensing module (MSM) (111) comprises an exterior (rigid) housing that receives the bodily vibrations, which are transmitted to one or more movement sensors disposed in fixed relation to the housing.
Non-limiting examples of motion sensing module (MSM) (111) include the seismocardiograph (SCG) (114) and the ballistocardiograph (BCG) (112) described elsewhere herein.
The motion sensing module (MSM) (111) outputs one or more channel outputs (e.g., 111a-a to f or 111b-a to f), each channel output reporting a different movement type (linear or rotation) and/or different direction (x, y, z).
The motion sensing module (MSM) (111) is configured for placement on the subject at a discrete location. The placement is typically on the skin, or over a material (e.g., garment) in contact with the skin.
A motion sensing module (MSM) (111) comprises one or more movement sensors for detection of body movements. The one or more movement sensors may be contained within the motion sensing module (MSM) (111) housing. A movement sensor detects linear or rotational movement. Each movement sensor may generate one or more channel outputs. The sampling rate may be between 40 Hz and 1000 Hz. The motion sensing module (MSM) may comprise a movement sensor that is an accelerometer for measurement of linear acceleration. The accelerometer may be a one-, two-, or three-dimensional accelerometer for measurement of linear acceleration in one, two, or three different directions (1, 2 or 3 linear degrees of freedom (1, 2, or 3DOF)). The motion sensing module (MSM) (111) may generate one, two, or three channel outputs depending on the number of degrees of freedom being measured. Preferably, the motion sensing module (MSM) (111) comprises a three-dimensional accelerometer for measurement of linear acceleration in three different directions (3 linear degrees of freedom (3DOF)); the three-dimensional accelerometer preferably generates three channel outputs (
Preferably, the motion sensing module (MSM) (111) comprises a three-dimensional gyroscope for measurement of rotational movement around three different axes (3 linear degrees of freedom (3DOF)); the three-dimensional gyroscope preferably generates three channel outputs (
For motion sensing module (MSM) (111) measurements, an x-axis may be a lateral axis (transverse axis: from the subject's left to the subject's right); a y-axis may be a longitudinal axis (caudocranial, oriented from feet to head), a z-axis may be an antero-posterior axis (ventrodorsal) (see
Preferably, the motion sensing module (MSM) (111) comprises at least two movement sensors, one being a one-, two-, or three- (preferably three) dimensional accelerometer and another being a one-, two-, or three- (preferably three) dimensional gyroscope. Preferably, the motion sensing module (MSM) (111) generates two or more channel outputs, one, two, or three (preferably three) from the accelerometer (
The motion sensing module (MSM) may further comprise a means for attachment to the body of the subject, e.g., a band or strap (e.g., for attachment to chest, head or waist, wrist etc), an adhesive, a garment (e.g., hat, sweat band, vest, glove, shorts, sock, shoe), or a wearable (e.g., harness, watch, glasses, ring, hat). It may be incorporated into a fitness tracker (e.g., chest band, watch, ankle strap).
The motion sensing module (MSM) may include an input for one or more electrocardiograph (ECG) electrodes for measurement of ECG signals. The motion sensing module (MSM) may include one or more interface modules for communication, for example, wireless or wired data transfer of the one or more channel outputs. Data transfer may be between the motion sensing module (MSM), and a smartphone, tablet, or other computing device. Data transfer may be between one motion sensing module (MSM) and one or more further motion sensing modules (MSM).
One example of a motion sensing module (MSM)(111) is a seismocardiograph (SCG)(114). The seismocardiograph (SCG) is configured for placement on the subject on the thorax, for instance, over the sternum, upper sternum (manubrium), or sternum body (corpus sternal).
The seismocardiograph (SCG) comprises one or more movement sensors for detection of bodily movements. A movement sensor detects linear or rotational movement. Each movement sensor may generate one or more channel outputs. The sampling rate may be between 40 Hz and 1000 Hz.
The seismocardiograph (SCG) may comprise a movement sensor that is an accelerometer for measurement of linear acceleration. The accelerometer may be a one-, two-, or three-dimensional accelerometer for measurement of linear acceleration in one, two, or three different directions (1, 2, or 3 linear degrees of freedom (1, 2, or 3DOF)). The seismocardiograph may generate one, two, or three channel outputs depending on the number of degrees of freedom being measured. Preferably, the seismocardiograph comprises a three-dimensional accelerometer for measurement of linear acceleration in three different directions (3 linear degrees of freedom (3DOF)); the three-dimensional accelerometer preferably generates output three channel outputs (
The seismocardiograph (SCG) may alternatively or additionally comprise a movement sensor that is a gyroscope for measurement of rotational movement. The gyroscope may be a one-, two-, or three-dimensional gyroscope for measurement of rotational movement in one, two, or three different directions (1, 2 or 3 linear degrees of freedom (1, 2, or 3DOF)). The seismocardiograph may generate one, two, or three channel outputs depending on the number of degrees of freedom being measured. Preferably, the seismocardiograph comprises a three-dimensional gyroscope for measurement of rotational movement in three different axes (3 linear degrees of freedom (3DOF)); the three-dimensional gyroscope preferably generates three channel outputs (
For seismocardiograph measurements, an x-axis may be a lateral axis (transverse axis: from the subject's right to the subject's left); a y-axis may be a longitudinal axis (caudocranial, oriented from feet to head); a z-axis may be an antero-posterior axis (dorsoventral, from back to stomach) (see
Preferably, the seismocardiograph (SCG) comprises at least two movement sensors, one being a one-, two-, or three- (preferably three) dimensional accelerometer and another being a one-, two-, or three- (preferably three) dimensional gyroscope. Preferably, the seismocardiograph generates two or more channel outputs, one, two, or three (preferably three) from the accelerometer (
The seismocardiograph may further comprise a means for attachment to the body of the subject, e.g., mounting to a bodily-worn belt, adhesive gel pad. The seismocardiograph may include an input for one or more electrocardiograph (ECG) electrodes for measurement of ECG signals. The seismocardiograph may include one or more interface modules for communication, for example, wireless or wired data transfer of the one or more channel outputs. Data transfer may be between the seismocardiograph and a smartphone, tablet, or other computing device. Data transfer may be between the seismocardiograph and the ballistocardiograph.
One example of a motion sensing module (MSM)(111) is a ballistocardiocardiograph (BCG)(112). The ballistocardiograph (BCG) is configured for placement on the subject at or around the centre of mass (or gravity) the subject.
The ballistocardiograph (BCG) comprises one or more movement sensors for detection of body movements. A movement sensor detects linear or rotational movement. Each movement sensor may generate one or more channel outputs. The sampling rate may be between 40 Hz and 1000 Hz. The ballistocardiograph (BCG) may comprise a movement sensor that is an accelerometer for measurement of linear acceleration. The accelerometer may be a one-, two-, or three-dimensional accelerometer for measurement of linear acceleration in one, two, or three different directions (1, 2 or 3 linear degrees of freedom (1, 2, or 3DOF)). The ballistocardiograph may generate one, two, or three channel outputs depending on the number of degrees of freedom being measured. Preferably, the ballistocardiograph comprises a three-dimensional accelerometer for measurement of linear acceleration in three different directions (3 linear degrees of freedom (3DOF)); the three-dimensional accelerometer preferably generates three channel outputs (
The ballistocardiograph (BCG) may alternatively or additionally comprise a movement sensor that is a gyroscope for measurement of rotational movement. The gyroscope may be a one-, two-, or three-dimensional gyroscope for measurement of rotational movement in one, two, or three different directions (1, 2 or 3 linear degrees of freedom (1, 2, or 3DOF)). The ballistocardiograph may generate one, two, or three channel outputs depending on the number of degrees of freedom being measured. Preferably, the ballistocardiograph comprises a three-dimensional gyroscope for measurement of rotational movement around three different axes (3 linear degrees of freedom (3DOF)); the three-dimensional gyroscope preferably generates three channel outputs (
For ballistocardiograph (BCG) measurements, an x-axis may be a lateral axis (transverse axis: from the subject's left to the subject's right); a y-axis may be a longitudinal axis (caudocranial, oriented from feet to head), a z-axis may be an antero-posterior axis (ventrodorsal, from stomach to back) (see
Preferably, the ballistocardiograph (BCG) comprises at least two movement sensors, one being a one-, two-, or three- (preferably three) dimensional accelerometer and another being a one-, two-, or three- (preferably three) dimensional gyroscope. Preferably, the ballistocardiograph generates two or more channel outputs, one, two, or three (preferably three) from the accelerometer (
The ballistocardiograph (BCG) may further comprise a means for attachment to the body of the subject, e.g., mounting to a bodily-worn belt, adhesive gel pad. The ballistocardiograph may include an input for one or more electrocardiograph (ECG) electrodes for measurement of ECG signals. The ballistocardiograph may include one or more interface modules for communication, for example, wireless or wired transfer of the one or more channel outputs. Data transfer may be between the ballistocardiograph and a smartphone, tablet, or other computing device.
The sensor unit (110) may be configured to generate a stream of signals (120) from the seismocardiograph (SCG) and optionally from the ballistocardiograph (BCG). The sensor unit (100) signal stream (120) may contain:
The sensor unit (100) signal stream (120) may contain channel outputs (115b) from the y-linear axis accelerometer axis of seismocardiograph (SCG). The sensor unit (100) signal stream (120) may contain channel outputs (115a-c) from the three accelerometer axes of ballistocardiograph (BCG) (all 3 linear (x, y, z) axes).
The sensor unit (100) signal stream (120) may contain channel outputs (113b) from the y-linear axis accelerometer axis of ballistocardiograph (BCG). The sensor unit (100) signal stream (120) may contain channel outputs (113a-c) from the three accelerometer axes of ballistocardiograph (BCG) (all 3 linear (x, y, z) axes).
It has been found that cardiorespiratory fitness measurement, in particular, the VO2 max value can be determined from as little as one seismocardiograph (SCG) axis of movement (e.g., y-linear axis, one channel output (115b)), or one ballistocardiograph (BCG) axis of movement (e.g., y-linear axis, one channel output (113b)), however, the more axes of movement that can be taken into account the more optimised the cardiorespiratory fitness measurement, in particular, the VO2 max value.
The sensor unit (110) may comprise (only) one motion sensing module (MSM)(111), for convenience to the subject. For ease of wearing, the sensor unit (110) may comprise (only) the seismocardiograph, SCG (114). In such case, the seismocardiograph (SCG) comprises at least one movement sensor; the movement sensor may be a one-, two-, or three-(preferably three) dimensional accelerometer, or a one-, two-, or three- (preferably three) dimensional gyroscope.
The sensor unit (110) may comprise both the ballistocardiograph, BCG, and the seismocardiograph, SCG. In such case, the ballistocardiograph (BCG) comprises at least one movement sensor being a one-, two-, or three- (preferably three) dimensional accelerometer, or a one-, two-, or three- (preferably three) dimensional gyroscope. The seismocardiograph (SCG) may comprise at least one movement sensor being a one-, two-, or three- (preferably three) dimensional accelerometer or a one-, two-, or three- (preferably three) dimensional gyroscope. The sensor unit (110) may generate two or more channel outputs:
One example of a sensor unit (110) comprising both SCG and the BCG is a “Kino” created by HeartKinetics; Gosselies, Belgium. It is made of two small modules that independently collect the SCG and BCG signals of the subjects together with a 1-lead ECG (one lead contains of electrodes with opposite polarity).
The processing unit (160) comprises circuitry configured performing the method of the invention. Typically, the circuitry comprises a processor and a memory. The processing unit (160) may be implemented in a computing device such as a desktop PC, laptop, dedicated programmable controller, as a collection of connected computing devices. The processing unit may be provided in part or entirely by a processor and a memory disposed within a housing of the sensor unit (110), or ballistocardiograph (112) or seismocardiograph (114), or within a housing of a separate unit. The processing unit (160) may be configured for performing one of more of the methods, or parts thereof, as described herein.
The processing unit (160) contain or be part of a standard computer system such as an Intel Architecture IA-32 based computer system 2, and may implement programming instructions of one or more software modules stored on non-volatile (e.g., hard disk or solid-state drive) storage associated with the corresponding computer system. However, it will be apparent that at least some of the steps of any of the described processes could alternatively be implemented, either in part or in its entirety, as one or more dedicated hardware components, such as gate configuration data for one or more field programmable gate arrays (FPGAs), or as application-specific integrated circuits (ASICs), for example.
A kinetic energy data stream (KEDS) is determined for each non-exercise period (N) from the signal stream (120). The kinetic energy data stream (KEDS) is determined from the sensor unit (100) signal stream (120) on one or more output channels. The determination of KEDS for each non-exercise period (N) is shown schematically in
The kinetic energy data stream (KEDS) may be determined from one or more of:
The sensor unit (100) signal stream (120) may be transformed into a kinetic energy data stream (KEDS) according to the methods described in Hossein et al “Accurate Detection of Dobutamine-induced Haemodynamic Changes by Kino-Cardiography: A Randomised Double-Blind Placebo- Controlled Validation Study”, Nature Scientific Reports (2019) 9:10479 (https://doi.org/10.1038/s41598-019-46823-3).
Linear kinetic energy (KLin) at time (t) for a single linear axis ((ax) which may be x, y, or z) may be determined according to the equation:
K
Lin(ax)(t)=0.5 [m·vax(t)2], where
Where there are multiple channel outputs for multiple linear axis (x, y, and z), the channel outputs are summed, viz:
Linear kinetic energy (KLin) at time (t) for a multiple linear axis (x, y, and z) may be determined according to the equation:
Where two of the three axes are used, the axis not being used is removed from the equation.
Rotational kinetic energy (KRot) at time (t) around a single linear axis ((ax) which may be x, y, or z) may be determined according to the equation:
K
Rot(ax)(t)=0.5 [I(ax)·W(ax)(t)2], where:
Where there are multiple channel outputs for rotations around multiple linear axis (x, y, and z), the channel outputs are summed, viz:
Rotational kinetic energy (KRot) at time (t) around multiple linear axes (x, y, and z) may be determined according to the equation:
Where two of the three axes are used, the axis not being is used is removed from the equation.
The kinetic energy data stream (KEDS) may be determined from one or more parameters obtained from the sensor unit (110) that is correlated to the linear kinetic energy and/or to the rotational kinetic energy (i.e. one or more correlated parameters).
An example of a correlated parameter is amplitude (also known as magnitude) of one or more channel outputs (e.g., 111a-a to f or 111b-a to f) of the motion sensing module (MSM) (111) of the sensor unit (110). In particular, an example of a correlated parameter is the amplitude expressed in acceleration, velocity, displacements, or force. Another example is signal energy of one or more channel outputs (e.g., 111a-a to f or 111b-a to f) of the motion sensing module (MSM) (111) of the sensor unit (110). In particular, an example of a correlated parameter is the signal energy expressed in acceleration, velocity, displacements or force.
The correlated parameter that is amplitude is a norm of a signal from one or more channel outputs (e.g., 111a-a to f or 111b-a to f) of the motion sensing module (MSM) (111) of the sensor unit (110). This signal may be an acceleration vector, a velocity vector, a displacement vector or a force vector. The units of acceleration vector are preferably m/s2. The units of the velocity vector are preferably m/s. The units of the displacement vector are preferably m. The units of the force vector are preferably N (Newton).
The amplitude at time (t) for a single linear axis ((ax) which may be x, y, or z) for a VoM may be determined according to the equation:
A
Lin(ax)(t)=sqrt [VoMax(t)2], where
Where there are multiple channel outputs for multiple linear axis (x, y, and z) for a VoM, the channel outputs are summed, viz:
The amplitude at time (t) for multiple linear axes (x, y, and z) for a VoM may be determined according to the equation:
A
Lin(x,y,z)(t)=sqrt[((VoMx(t)2)+(VoMy(t)2)+(VoMz(t)2))], where
Where two of the three axes are used, the axis not being is used is removed from the equation.
The amplitude at time (t) around a single linear axis ((ax) which may be x, y, or z) for a VoM may be determined according to the equation:
A
Rot(ax)(t)=sqrt [VOMax(t)2], where
Where there are multiple channel outputs for multiple linear axis (x, y, and z) for a VoM, the channel outputs are summed, viz:
The amplitude at time (t) around multiple linear axes (x, y, and z) for a VoM may be determined according to the equation:
A
Rot(x, y, z)(t)=sqrt[((VoMx(t)2)+(VoMy(t)2)+(VoMz(t)2))], where
Where two of the three axes are used, the axis not being is used is removed from the equation.
The correlated parameter that is signal energy is an area under a squared amplitude (or magnitude) of a signal from one or more channel outputs (e.g., 111a-a to f or 111b-a to f) of the motion sensing module (MSM) (111) of the sensor unit (110). This signal may be an acceleration vector, a velocity vector, a displacement vector or a force vector. The units of acceleration vector are preferably m/s2. The units of the velocity vector are preferably m/s. The units of the displacement vector are preferably m. The units of the force vector are preferably N (Newton).
The signal energy at time (t) for a single linear axis ((ax) which may be x, y, or z) for a VoM may be determined according to the equation:
E
Lin(ax)(t)=integral(ALin(ax)(t)2),
The signal energy at time (t) for multiple linear axes (x, y, and z) of a VoM may be determined according to the equation:
E
Lin(x, y, z)(t)=integral(ALin(x, y, z)(t)2), where
Where two of the three axes are used, the axis not being is used is removed from the equation.
The signal energy at time (t) around a single linear axis ((ax) which may be x, y, or z) for a VoM may be determined according to the equation:
E
Rot(ax)(t)=integral(ARot(ax)(t)2),
The signal energy at time (t) around multiple linear axes (x, y, and z) of a VoM may be determined according to the equation:
E
Rot(x, y, z)(t)=Integral(ARot(x, y, z)(t)2), where
Where two of the three axes are used, the axis not being is used is removed from the equation.
A mean kinetic energy (MKE) that is a time-averaged kinetic energy value from a time window (W) of the kinetic energy data stream is obtained. The time window (W) is for a period of w seconds. The value of w is less than the duration of the non-exercise period (N). Typically, w is greater than a duration of 1 cardiac cycle of the subject. The value of w may be in a range of 2 to 10 seconds. The value of w is preferably the same during each non-exercise period (N). The value of w is preferably the same during each non-exercise period (N) and during the exercise session (E).
A mean kinetic energy (MKE) may be determined from a time window (W) of the kinetic energy data stream (KEDS) by any method to compute the mean of discrete kinetic energy values (DKEV) over a given time window. For instance, the mean kinetic energy (MKE) may be determined from an integration of the discrete kinetic energy values (DKEV) measured within the time window (W) of the kinetic energy data stream (KEDS). For instance, the mean kinetic energy (MKE) may be determined from a time average of the discrete kinetic energy values (DKEV) measured within the time window (W) of the kinetic energy data stream (KEDS).
A period of mean kinetic energies (PMKEs) is determined within each non-exercise period (N), by determining each mean kinetic energy (MKE) of the multiple at a different starting time points within the non-exercise period (N). Some periods of mean kinetic energies (PMKEs) may overlap but have different starting time points. Preferably, a first mean kinetic energy (MKE) is determined at the start of the non-exercise period (N), a last mean kinetic energy (MKE) is determined at the end of the non-exercise period (N), and intermediate mean kinetic energies (MKEs) are determined at regular intervals between the first and last MKEs; this plurality of MKEs constitute the period of mean kinetic energies PMKEs. The time difference between consecutive time windows (W) and hence mean kinetic energy (MKE) may be constant or different within a non-exercise period (N). For instance, consecutive time windows (W) may differ in starting time of 0.5 to 2 seconds.
An example of a plurality of mean kinetic energies (MKEs) (MKEt0, MKEt1, MKEt2, MKEtp) determined from different time windows (Wt0, Wt1, Wt2, Wtp) during a non-exercise period (N) is shown in
A plurality of PMKEs (periods of mean kinetic energy) obtained during the course of the exercise session (E), one PMKE per non-exercise period, is known as set of PMKEs (SPMKE). The determination of PMKEs for each non-exercise period (N) within the exercise session (E) (the SPMKE) is shown schematically in
The cardiorespiratory fitness measurement of the subject can be determined from the set of PMKEs (SPMKE). The SPMKE is indicative of the oxygen consumption rate (VO2) over time for the exercise session (E). From the oxygen consumption rate (VO2) over time, the VO2 max value can be determined, which as mentioned elsewhere is one measurement of cardiorespiratory fitness.
In particular, the SPMKE may be transformed into a linear VO2 function (LVO2F) for the subject, which VO2 function (LVO2F) is representative of the oxygen consumption rate (VO2) as a function of time of the subject during the exercise period. The determination of LVO2F from the SPMKE within the exercise session (E) is shown schematically in
In order to transform the SPMKE into the LVO2F, each PMKE of the set of PMKEs (SPMKE) is converted into a scalar value (single value) representative of the non-exercise period (N), wherein the scalar value is an order statistic. Thus, a set of order statistic (SOS) for the exercise session (E) is obtained, and from the evolution over time of the order statistics (OSs) of the SOS, the linear VO2 function (LVO2F) for the subject can be determined. The determination of LVO2F from the SPMKE within the exercise session (E) using a set of order statistic (SOS) is shown schematically in
An order statistic (OS) is determined for each non-exercise period (N). An order statistic (OS) may be a median, or another quantile (e.g., 1st and 3rd quartile).
To generate the order statistic (OS) for each non-exercise period (N), each PMKE of the set of PMKEs (SPMKE) is converted into a scalar value (single value). This conversion may be performed:
In the case of transforming indirectly the PMKE into an order statistic (OS), the period of mean kinetic energies (PMKEs) within a non-exercise period (N) is transformed into an intermediate linear function (ILF); this may be achieved using a scaling law. For instance, each mean kinetic energy (MKE) from the period of mean kinetic energies (PMKEs) may be plotted on a graph (y-axis: MKE, x-axis: time), to which a curve can be fitted. The fitted curve is transformed into an intermediate straight line (the intermediate linear function) by the scaling law. The intermediate linear function may be expressed graphically or as a mathematical formula. The scaling law may be any function that linearises a curve. In particular, the scaling law may use a cubic root of the mean kinetic energy (MKE). In particular, the scaling law may use the heart rate (HR) of the subject. One example of a scaling law comprises the equation HR*CubicRoot(MKE/HR).
In the case of transforming directly the PMKE into an order statistic, each mean kinetic energy (MKE) from the period of mean kinetic energies (PMKEs) within a non-exercise period (N) is arranged in value order (e.g., largest to the smallest or vice versa). From the value-ordered mean kinetic energies (MKEs), the median (mid-value), or quantile(s) can be determined.
The same order statistic (median or quantile) is determined for each non-exercise period (N), and at least some of these order statistics (OS) form a set of order statistics (SOS) for the exercise session (E). For example, the set of order statistics (SOS) may contain only median order statistics (and not a mix of median and quantile order statistics).
The set of order statistics (SOS) preferably contains only order statistics (OS) showing a linear evolution over time by the subject. The order statistics (OS) for inclusion in the set of order statistics (SOS) may be determined by including the first three order statistics (OS) related to the three first non-exercise periods (N) (excluding the very first measurement taken before the first exercise period), then progressively adding the following order statistics and evaluating if the new set of order statistics (SOS) has a better linear correlation coefficient (e.g., Pearson Correlation Coefficient) than the previous one. The set of order statistics (SOS) is considered complete whenever the addition of a new order statistics results in a decrease of the linear correlation coefficient. Typically, the last (by time) order statistics (OS) are not included in the set of order statistics (SOS) because the subjects reach a point of non-linearity. Advantageously, measuring cardiorespiratory fitness, in particular, the VO2 max based on a submaximal effort is more comfortable (less exhausting) for the subject. The cardiorespiratory fitness measurement, in particular, the VO2 max can be determined from the set of order statistics (SOS). As elaborated elsewhere herein, cardiorespiratory fitness measurement, in particular, the VO2 max can be determined from a set of order statistics (SOS) by transformation into a VO2 function (VO2F) for the subject In order to transform the SOS into the LVO2F, a proportionality function is employed. The proportionality function is applied to the evolution over time of the order statistics (OSs) of the SOS to arrive at the linear VO2 function (LVO2F) for the subject using a proportionality function. The proportionality function is described in more detail below.
Once the linear VO2 function (LVO2F) is obtained for the subject, the value on the linear VO2F curve at a point that corresponds to the actual (or estimated) maximal load is determined. This can be based on observation of a plateau in the order statistics (OS). Alternatively, it can be based on the point of the linear VO2 function (LVO2F) corresponding to the actual (or estimated) maximal heart rate.
This estimated value may be corrected by a constant (‘d’). The constant (d) corrects for recovery by the subject during the non-exercise period that would not arise if measurement has been made during an exercise session (S).
This transformation of the SOS into the LVO2F by the proportionality function may be performed
In the case of transforming directly the SOS into the LVO2F, the proportionality function is applied to the SOS. For instance, the OS values of the SOS may be plotted on a graph (y-axis: OS (t), x-axis: time) as a first straight line, and the SOS is transformed into the LVO2F by the proportionality function. The VO2 function (VO2F) may be expressed graphically or as a mathematical formula. In
In the case of transforming indirectly the SOS into the LVO2F, the set of order statistics (SOS) may be first transformed into an intermediate linear function (ILF) for the exercise session (E); this may be achieved using a scaling law. For instance, the set of order statistics (SOS) may be plotted on a graph (y-axis: order statistic value, x-axis: time), to which a curve can be fitted. The fitted curve is transformed into a first straight line (the linear function) by the scaling law. The linear function may be expressed graphically or as a mathematical formula.
Then, the intermediate linear function (ILF) is transformed into the VO2 function (VO2F) for the subject; this may be achieved using the proportionality function. The proportionality function is described in more detail below. The VO2 function (VO2F) is representative of the oxygen consumption of the subject as a function of time during the exercise session (E). For instance, the intermediate linear function (ILF) may be plotted on a graph (y-axis: LF (t), x-axis: time) as a first straight line, and the LF is transformed into a second straight line (VO2 function) by the proportionality function. The intermediate linear function (ILF) is transformed into the VO2 function (second straight line) by the proportionality function. The VO2 function (VO2F) may be expressed graphically or as a mathematical formula. In
As mentioned earlier, the set of order statistics (SOS) or intermediate linear function (LF) is transformed into a VO2 function (VO2F) for the subject; this may be achieved using a proportionality function.
The proportionality function may include a step of multiplication by a constant (Q). The constant (Q) may be fixed for a subject; it may be different from one subject to another.
One example of a proportionality function comprises the equation Q(LF). One example of a function for determining the LVO2F is
The value of the constant Q may be determined by calibrating the subject using an ergospirometry measurement. After this first measurement, the value of Q and additional VO2 max measurements may be consulted on the same subject. Alternatively, the value of Q based on parameters from the subjects (height, weight, body composition, cross-section area of the aorta, blood composition). The value of d, as mentioned elsewhere herein, corrects for recovery by the subject during the non-exercise period that would not arise if measurement has been made during an exercise session (S).
When it is mentioned that a kinetic energy data stream (KEDS) is determined for each non-exercise period (N) from the signal stream (120), it is understood that one or more kinetic energy data streams (KEDS) may be generated, depending on whether one or multiple channel outputs are being used.
For instance, where only one channel output (e.g., BCG Lin(y), 113b) of the signal stream (120) is being used, only one kinetic energy data stream (KEDS) may be determined from the one channel (e.g., BCG Lin(y), 113b). Where multiple (two or more) channel outputs (e.g., BCG Lin(x), 113a; BCG Lin(y), 113b; BCG Lin(z), 113c) are being utilised, only one kinetic energy data stream (KEDS) may be determined by combining the channel outputs. For instance, one kinetic energy data stream (KEDS) that is a single BCG linear kinetic energy (KLin(BCG)) may be determined from a combination of BCG channel outputs BCG Lin(x), 113a; BCG Lin(y), 113b; BCG Lin(z), 113c. Other combinations of channel outputs are foreseen: combining KLin (BCG) with: KLin (SCG), and/or with KRot (BCG), and/or with KRot (SCG). It is foreseen that one kinetic energy data stream (KEDS) may be determined by combining any two of more channel outputs (113a-f; 115a-f)). It is described elsewhere exemplary equations for generating a KEDS from one channel output or multiple channel outputs.
The subsequent steps of determining period of mean kinetic energy (PMKEs), order statistics, transformations using a scaling law and a proportionality function to arrive at the VO2 function are hence performed using only one channel output (e.g. BCG Lin(y), 113b) from the sensor unit (100) signal stream (120).
Multiple kinetic energy data stream (KEDSs) may be determined from multiple (two or more) channel outputs (e.g., BCG Lin(x), 113a; BCG Lin(y), 113b; BCG Lin(z), 113c). For example, one KEDS may be determined from each of the multiple channel outputs (e.g., one KEDS for each of BCG Lin(x), 113a; BCG Lin(y), 113b; BCG Lin(z), 113c).
Where multiple KEDS are determined, the subsequent steps of determining periods of mean kinetic energies (PMKEs), order statistics, transformations using a scaling law and a proportionality function to arrive at the VO2 function may be performed for each of the different channel outputs (113a-f; 115a-f), to arrive at a plurality of VO2Fs, one per channel output (113a-f; 115a-f). The steps after generating the plurality of VO2Fs may combine channel output transformations, for instance, by summing or averaging, to arrive to a single function or value from which the cardiorespiratory fitness measurement, in particular, the VO2 max value can ultimately be calculated.
For instance, a plurality of VO2Fs (each generated from a different channel output (113a-f; 115a-f)) may be combined, e.g., by averaging (weighted or non-weighted), to arrive at a single VO2F. The cardiorespiratory fitness measurement, in particular, the VO2 max value may be determined from the combined VO2F.
For instance, a plurality of the cardiorespiratory fitness measurements, in particular, of VO2 max values (each generated from a different channel output (113a-f; 115a-f)) may be combined, e.g., by averaging (weighted or non-weighted), to arrive at a single cardiorespiratory fitness measurement, in particular, a single VO2 max value.
The processing steps of the present method are performed in vitro.
The system (100) may include the sensor unit (110).
Further provided is a computer-implemented method for determining a cardiorespiratory fitness measurement, in particular, a maximal oxygen consumption, VO2 max, of a subject comprising the steps carried out by the processing unit (160) as described elsewhere herein.
Further provided is a computer program or computer program product having instructions which when executed by a computing device or system cause the computing device or system to perform the computer-implemented method.
Further provided is a computer readable medium having stored thereon the computer program or computer program product.
Further provided is a data stream which is representative of the computer program or computer program product.
16 healthy subjects were recruited (8 males and 8 females). The demographics of the test subjects are presented in Table 1.
These subjects performed a standard VO2 max protocol using bicycle ergospirometry (ESt), as well as an exercise session (E1) as described herein including alternating exertion segments (S1) and non-exercise periods (N1). These two protocols are described in Table 2 and include Kino records in parallel with traditional analysis of gas exchanges.
The results of VO2 max given by ergospirometry (see Table 3) show that the test protocol (E1) gives results similar to the standard protocol (ESt). A Wilcoxon matched test gives indeed a p-value of 1.00, strongly supporting the equivalence of the two protocols.
To estimate the value of VO2 max based on the Kino records, we choose to focus on a kinetic energy data stream (KED) computed based on the linear axes of BCG:
During each non-exercise period (N), the mean kinetic energy (MKE) was computed by integration of the KEDS signal over a window W of 5 seconds and normalisation by the length of this window. This mean kinetic energy (MKE) was computed by shifting the window every 1 second during a given non-exercise period (N) and resulted in a period mean kinetic energy (PMKE) for each non-exercise period (N).
The PMKE for each non-exercise period (N) was transformed into a linear function (LF) using the following scaling law:
LF=HR·CubicRoot(MKE/HR)
where HR is the heart rate of the subject during the associated time-window W.
One LF was determined for each PMKE and hence for each non-exercise period (N1). An LF of a non-exercise period (N1) is populated by a block of LF values. From the block of LF values, an order statistic (OS) was extracted, given by the median of the block of LF values within the LF. For the exercise session (E1) a set of order statistics (SOS) was obtained, the SOS containing a plurality of order statics (medians), each order statics obtained from a block of LF values. An example of a set of order statistics (SOS) for an exercise session (E1) is shown in
In
For each of the SOS points that are in the linear section, the ratio between the median value of measured VO2 and OS was determined. The median of these ratios gave the proportionality coefficient Q. The SOS curve was then transformed in a VO2F line, by multiplying each value by Q and fitting a line.
The value of the VO2F line at the maximal load reached by the subject was evaluated, which gave a preliminary evaluation of VO2 max.
The value of VO2 max was further optimised to correct for a slight underestimate, because it is based on measurements performed during non-exercise periods (N), during which the subject had some time to recover (depending on the duration of the non-exercise period and the level of fitness of the subject). A correction factor ‘d’ for each subject was determined, based on ergospirometry data. Alternatively, it was found that this correction factor was positively correlated to the value of the estimated VO2 max, and may thus be approximated using a linear relationship.
The results of the VO2 max estimation with a known maximal workload and correction factor ‘d’ are given in Table 4. The estimates are very accurate, with a relative error that is mostly under 10%. A scatter plot of the estimated vs. measured value of VO2 max is also given in
We also verified the effect of missing input information, by estimating the maximum workload based on measured maximum heart rate and based on estimated maximum heart rate (220—age). Each approximation leads to a decrease in the accuracy of the estimates, however the classification of the subjects according to their global fitness level remains very good (results not presented).
Additional experiments (results not presented) show that the proportionality coefficient Q is constant for a given subject.
| Number | Date | Country | Kind |
|---|---|---|---|
| 22156055.0 | Feb 2022 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2023/053302 | 2/10/2023 | WO |
