This invention relates to the field of plethysmography, and, in particular, to a system and method of extracting features from a plethysmographic waveform to determine various characteristics of a human subject.
Water is vital in the life of all living things. It plays a central role in such functions as nutrient transport, waste removal, maintenance of cell volume, and thermal regulation. Assessment of hydration provides insight into basic biological processes and therefore is an important clinical parameter in determining health status and treatment of disease. A loss of only 2-3% of hydration by body weight is enough to lead to loss of physical performance, a loss of 6% can have serious repercussions, including heat exhaustion, and a loss of 8% and greater can lead to delirium, temporary deafness, cessation of urination and, at greater levels, death. Both dehydration and over hydration can have tragic consequences. Thus, maintaining a normal hydration state is a primary concern in all aspects of medicine. Although there are no quantitative measurements of hydration currently being used, qualitative parameters such as skin turgor, urine production, and the presence of tears are commonly used to determine a patient's hydration status. Although the laboratory setting lends itself to more quantitative assessments, accurate measurements of hydration are often difficult and time consuming. A more readily attainable, accurate, and quantitative measure of hydration would provide scientists with a powerful tool in studying health and disease.
The amount of water in the blood is perhaps the most relevant indicator of the overall hydration status of the body, as blood is delivered to almost every tissue in the body. Water is also contained in the cells of the body, referred to as intracellular fluid, which helps the cells maintain their shape, and facilitates the myriad processes that a living organism needs. Water is also present outside the cells, referred to as the extracellular reservoir. The remainder of the water in a human body is contained in the blood plasma. The water in the different compartments of the body's reservoirs is constantly shifting back and forth in response to changing physiological and environmental conditions, with blood generally being the vehicle of redistribution. Blood is also the gateway to many of the body's water loss mechanisms such as the kidney, breath, and sweating. Therefore, measuring the water content of the blood gives a good relative measurement of the hydration status of the individual.
It is well known in the art that the pumping of blood by the heart causes a pressure wave to form in the arteries and arterioles in the subcutaneous tissue, known as a plethysmographic waveform. This waveform can be detected by illuminating the skin with the light from a light-emitting diode (LED) and then measuring the amount of light either transmitted or reflected to a photodiode. When the plethysmographic waveform is detected in this manner, it is referred to as a photoplethysmogram.
Such plethysmographic waveforms are typically collected and used by a pulse oximetry device to detect oxygen saturation in the blood of a human subject. Pulse oximeters typically use LEDs transmitting light through the skin in two wavelengths, typically in the red at about 660 nm and in the infrared at about 940 nm. Waveforms at each wavelength can be collected using a transmissive method, where the light from each LED is passed through a translucent portion of the subject's body and detected by a photodiode on the other side. The absorption of light at these wavelengths differs significantly between blood loaded with oxygen and blood lacking oxygen. Oxygenated blood absorbs more infrared light and allows more red light to pass through, while deoxygenated blood allows more infrared light to pass through and absorbs more red light. The detected light from each wavelength can be compared to derive a plethysmographic waveform, an example of which is shown in
The present invention relates to using feature extraction from a plethysmographic waveform to determine the hydration status of an individual. The plethysmographic waveform contains information about a host of physiological properties, including the relative amount of water in the blood. This invention relates to a method for measuring the amount of water in the blood by algorithmic feature extraction of the plethysmographic waveform. A method of determining the subject's pulse rate is also detailed in this invention.
The sensing of a plethysmographic waveform is based on the principal that absorption of a specific wavelength of light energy in tissue varies with the amount of oxygenated blood in vessels such as arteries or arterioles. When the heart beats, the volume of blood increases and travels as a pressure wave through the circulatory system. When this pressure wave passes the sensor, more of the light energy is absorbed, and after the pressure wave has passed less light energy is absorbed. This time variant signal can be detected using a sensor consisting of a light source combined with a photodetector. Amplifying the signal gives an electrical representation of the plethysmographic waveform.
This invention relates to a sensor system capable of recording and analyzing the plethysmographic signal by means of either transmitted or reflected light energy, either separately or in combination. This system uses a single monochromatic light source, for example an LED, between the wavelengths of 400 nm and 2000 nm and a photo-detector, an example being a PIN photo diode. Light is projected into oxygenated tissue, and is transmitted or reflected (or both) into a photo-detector, producing a modulated electrical signal which is recorded as a photoplethysmographic waveform.
Another aspect this invention also relates to the use of a motion sensor, such as an accelerometer or gyroscope, to compensate for distortion in plethysmographic sensor signals induced by motion of the subject.
In accordance with the above, a system and algorithmic methods to measure the relative water content of the blood includes a sensor and a system capable of measuring, recording and performing analysis on a plethysmographic waveform such that relevant features can be extracted to measure the relative amount of water in the blood. This measurement allows for the determination of the relative hydration status of an individual, as well as their pulse rate.
In one embodiment of the invention, a system for the non-invasive measurement and recording of a photoplethysmographic waveform includes a light source with a narrow spectral bandwidth in combination with a photodetector capable of sensing light energy within the spectral bandwidth of the light source. The light energy is projected into tissue with oxygen carrying blood vessels, and either transmitted, reflected or both back through tissue to the photodetector. The light energy reaching the photodetector contains information of the complete plethysmographic waveform. A means for analyzing the waveform to extract biometric information such as hydration status and heart rate is in communication with the photodetector for receiving information about the plethysmographic waveform.
In a second embodiment of the invention, the system is also provided with a motion sensor, such as an accelerometer or gyroscope, that measures the relative motion of the sensing environment, which tends to distort the plethysmographic waveform. This motion information can subsequently be subtracted from the distorted plethysmographic waveform, thereby allowing for the recovery of the true plethysmographic waveform from the distorted waveform.
This invention relates to using feature extraction from a plethysmographic waveform to determine the hydration status of an individual. The plethysmographic wave is the pressure wave present in the arteries and arterioles of the body that is generated by the beating of the heart. This waveform contains information about a host of physiological properties, including the relative amount of water in the blood. It has been found that an inference can be made as to water content of the entire body based on the water content in the blood plasma.
This invention comprises the use of one or more of six different algorithmic methods to extract the salient features from the plethysmographic waveform for the measurement of the amount of water in the blood. These methods are:
Each of these methods is a mathematical calculation which can be used to analyze the salient features of the plethysmographic waveform to determine the overall hydration status of the human body. In each case, these mathematical operations provide a number from which an inference can be made as to the amount of water in the blood.
Before the plethysmographic waveform can be analyzed, the raw data needs to be processed to remove noise and baseline wander. The data may also be inverted to make the waveform conform to the normal plethysmographic standard, but this for convenience and is not an absolute requirement for the processing of the data. The plethysmographic waveform can be processed in the inverted or non-inverted form.
Noise in the waveform can occur due to either electrical or mechanical factors, such as 60 Hz power line noise or a mechanical vibration induced into the patient. Baseline wander is the slowly changing DC bias variation caused by subtle changes in the DC bias of the sensor amplifiers, which causes a DC offset in the waveform.
Because the frequencies of interest in the plethysmographic waveform related to this invention are below 15 Hz, proper selection of sample frequency and low pass filtering of the signal are usually sufficient to remove the electrical noise from the signal. Baseline wander is removed by subtracting the long term running average of the signal.
Once the signal has been processed such that it is symmetrically periodic around zero, with all DC bias removed, the waveform is ready for feature extraction to determine the amount of water in the blood.
Those skilled in the art will understand that practical methods of measuring and recording a biometric waveform, such as a plethysmogram, involve analog to digital sampling of the waveform which has the effect of breaking the waveform into discrete samples representative of the original, inherently analog waveform. This discrete data is now ready for analysis by a microcontroller or a microprocessor.
The first method of feature extraction related to this algorithm is to obtain a frequency spectrum of the waveform through the use of Fourier transform analysis. A Fourier transform converts time or space information to frequency. Of particular usefulness, is the method known as the fast Fourier transform (FFT), as it can be used to process a discrete dataset in a computationally fast and efficient manner. This is the preferred embodiment of the invention. Equation 1 shows the mathematical relation between the Fourier transform and the function f(x) that represents the waveform.
f(ξ)=∫−∞∞f(χ)e−2πixtdx, for any real number ζ (1)
An FFT is an algorithm to compute the discrete Fourier transform (DFT) and it's inverse. The DFT is defined by the formula shown in Equation 2.
An FFT computes the DFT and produces the same result as evaluating the DFT definition directly; the most important difference is that an FFT is much faster. The most common method is the Cooley-Tukey FFT algorithm, but others include the Prime-factor FFT algorithm, Bruun's FFT algorithm, Rader's FFT algorithm, and Bluestein's FFT algorithm. Algorithmically, and as related to this invention, the magnitude of the maxima of the frequency spectrum returned by a Fourier transform, specifically the first (primary) harmonic of a plethysmographic waveform, is inversely proportional to the hydration status of the individual the waveform is taken from.
Additionally, the spectrum produced by the Fourier transform can be used to measure the pulse rate of an individual. The first harmonic of the spectrum is the fundamental frequency and therefore represents the heart rate. To determine heart rate, the first harmonic maxima can be multiplied by sixty with the resulting number being the heart rate in beats per minute, as shown in Equation 3.
HR=f1st(60) (3)
Root Mean Square
The first step in the feature extraction related to this algorithm is to determine the Root Mean Square (RMS) of the plethysmographic waveform. In mathematics, the root mean square, also known as the quadratic mean, is a statistical measure of the magnitude of a varying quantity. The RMS of a waveform can be calculated for either a continuous or a discrete signal. Equation 4 shows the mathematical relation for continuous case where the waveform function RMS is between times T1 and T2, and Equation 5 shows it for the discrete case with n set of values.
In either case, the result returned is a real number that can be used to analyze a waveform or set of data. Algorithmically, as related to this invention, the RMS value of the plethysmographic waveform is inversely proportional to the hydration status of the individual from which the waveform is taken.
The derivative of an equation or a waveform gives the slope of the line that is tangential to the function or waveform at that specific point.
Algorithmically, as related to this invention, the magnitude of the derivative of the plethysmographic waveform, meaning the slope of the line at any given point in the waveform, is inversely proportional to the hydration status of the individual from which the waveform is taken. The maximum derivative values can be used alone or in combination with the other maximum derivative values to assess hydration.
The crest factor is the ratio of the systolic peak amplitude to the RMS of the waveform, as shown in Equation 7. The magnitude of the crest factor measurement is inversely proportional to the hydration status of the individual.
The peak to average power (PAPR) ratio is the square of the crest factor and is shown in Equation 8. The magnitude of the PAPR measurement is inversely proportional to the hydration status of the individual.
A single wavelength optical biometric sensor system such as described and used to measure a plethysmographic waveform is prone to noise and signal distortion due to physical movement of the subject. The time variance and frequencies of such movement are on the same order as that of biological signals, with the plethysmographic pressure wave being a preferred example. Therefore, in addition to the sensing and recording of a plethysmographic waveform, the present invention comprises a motion sensor in direct proximity of the apparatus which senses the transmissive or reflective light. The sensed motion information can be used to compensate for motion distortion of the plethysmographic waveform, allowing for the recovery of the biometric plethysmographic signal.
One of three preferred methods of extracting the hydration information from the plethysmographic waveform is then utilized. It should be noted that not all three algorithms will be programmed into the system. This portion of
At 922, the derivative method, discussed above, is utilized. At 922, the first derivative of the plethysmographic waveform is derived, showing the slope of the curve. At 923, the second derivative is calculated, showing the rate of change of the slope curve calculated in 922. The rate of change of the derivative of the plethysmographic waveform, meaning the second derivative of the line at any given point in the waveform, is inversely proportional to the hydration status of the individual from which the waveform has been obtained.
At 924, the FFT method is used. The FFT method is the preferred embodiment of the invention. At 924, the FFT of the plethysmographic waveform is calculated and, at 925, the pulse rate is determined the based on the frequency of the first harmonic of the frequency spectrum, obtained using a Fourier Transform or another algorithm.
The hydration status of the individual is based on the magnitude of the peak frequency (i.e., the primary harmonic). Generally, the lower the magnitude, the more likely it is that the subject is de-hydrated. The exact magnitude to distinguish between a properly-hydrated state and a state of dehydration varies for each individual. Therefore, the system of the present invention will need to be calibrated on an individual basis. The FFT method for determining hydration status has been demonstrated through hydration testing and a small clinical study. The magnitude of the FFT primary harmonic correlates to the amount of blood volume, of which water is the largest component.
At 924, the last preferred method of interpreting the plethysmographic waveform is utilized. The RMS of the plethysmographic waveform is calculated. Generally, the RMS value of the plethysmographic waveform is inversely proportional to the hydration status of the individual. As with the other methods, the system of the present invention will need to be calibrated for each individual.
Regarding the processing of the plethysmographic waveform to remove motion artifacts,
However,
A high level block diagram of the hardware components of the invention is shown in
The foregoing invention has been explained in terms of function and purpose, and, as would be realized by skill in the art, may be implemented using many methods. For example, the portion 2106 of the invention that manipulates the plethysmographic waveform and the motion waveform may be implemented as software running in a general purpose computer, as hard circuitry, as a microprocessor running firmware, as a functionally programmed ASIC, or by using any other method known in the art. In any case, the sensor portion of the invention 2102, namely the LED(s), photodetector(s) and motion detector need to be configured to be in communication with the analysis components of the implementation. This can be accomplished via A-to-D converter 2104, which transforms the raw analog signal produced by the photodiodes and motion detector into a digital signal for processing by microcomputer 2106.
In one implementation of the invention, the device is entirely self-contained, including a means for communicating the results of the hydration analysis to the subject. Other possible implementations of the invention could communicate raw sensor data, a digitized version of the raw sensor data, or completely analyzed results to another device, for example, a mobile device such as a tablet or smartphone via, for example, WiFi or Bluetooth. The mobile device could provide the mathematical analysis of the raw sensor data or could be used merely to display results.
In addition, there will need to be a method of communicating the hydration status of the subject. This could be a simple as an LED indicating a dehydrated status, or more complicated, such as an LCD readout providing a numerical indication of the hydration status, or through the use of a “fuel gauge”, where a graphical indicator gives an indication of hydration status. In cases where the hydration status is displayed by a mobile computing device, more complicated feedback could be provided, for example, a graph of the subject's hydration status throughout the day or over longer periods. The system may also provide suggestions as to how the subject may improve his or her hydration status.
This application claims the benefit of the following co-pending U.S. Provisional Applications: U.S. Provisional No. 61/974,063, filed Apr. 2, 2014, U.S. Provisional No. 62/000,588, filed May 20, 2014 and U.S. Provisional No. 62/007,671, filed Jun. 4, 2014.
Number | Date | Country | |
---|---|---|---|
61974063 | Apr 2014 | US | |
62000588 | May 2014 | US | |
62007671 | Jun 2014 | US |