The present invention relates to the field of determining vital signs of a subject, in particular, to a processing apparatus, processing method and system for processing a physiological signal.
Vital signs of a subject are powerful indicators in determining a medical condition and fitness of a subject. Vital signs include but are not limited to a respiratory rate (RR) and a heart rate (HR) or pulse rate. The vital signs are descriptive of an underlying physiological process such as heart-beats or a respiratory movement. A physiological signal descriptive of the underlying physiological process can be measured and evaluated.
For example, respiration monitoring can be based on different principles: The measurement of respiratory effort, for example, thorax impedance plethysmography, respiratory inductance plethysmography, often referred to as respiration band or respiband, photoplethysmography or the measurement of respiratory effect, for example, sound recording, temperature sensing, carbon dioxide sensing.
Exemplary sensors for measuring the physiological signal for determining the heart rate comprise a photoplethysmograph sensor, an ECG (electrocardiography) sensor, or a heart sound sensor.
US 2004/0039420 A1 further discloses a pulse detection apparatus that uses physiological signal data obtained from an accelerometer placed on a patient's body to detect the presence of a cardiac pulse. The accelerometer is adapted to sense movement due to a cardiac pulse and produce accelerometer signal data in response thereto as the physiological signal. Processing circuitry analyzes the accelerometer signal data for a feature indicative of a cardiac pulse and determines whether a cardiac pulse is present in the patient based on the feature. The apparatus can be used in a defibrillator for delivering defibrillation pulses to a patient to stimulate a cardiac activity.
The feature indicative of a cardiac pulse can be a temporal energy feature, for example, a peak-to-peak amplitude, a peak-peak amplitude of a derivative, or an energy of the accelerometer signal. Alternatively, the feature can be a spectral energy feature such as the energy or frequency of a peak in the energy spectrum of the signal. In yet another aspect, the feature may be obtained by comparing the accelerometer signal data with a previously-identified pattern known to predict the presence of a cardiac pulse. This is also referred to as pattern matching, wherein the measured physiological signal is compared with a known pattern.
US 2006/0224074 A1 discloses a heartbeat measuring apparatus using pattern matching. Instead of comparing waveform data with a known pattern, a section of waveform data corresponding to one heartbeat is used as a template. A waveform similarity such as correlation between said template and the measured waveform is calculated and used for subsequent analysis.
US 2010/0312075 A1 further discloses signal processing techniques for aiding the interpretation of respiration signals. The peaks of a signal can be rescaled to yield a normalized waveform.
Since the vital signs are commonly derived from such characteristic features of the physiological signal, there is a need to determine whether the features, for example signal peaks, have been identified correctly. Current approaches to classify whether the peak-detections were done correctly evaluate, e.g., a heart-rate variability or perform a spectral analysis. A prerequisite for applying these techniques is an at least weak-stationary signal.
It is an object of the present invention to provide a processing apparatus, processing method and system for more accurately determining a quality of signal features, in particular signal features of a non-stationary physiological signal. It is a further object of the present invention to provide a processing apparatus, processing method and system for distinguishing different contributions in a physiological signal.
In a first aspect of the present invention a processing apparatus for processing a signal is presented that is configured to perform the steps of:
obtaining the physiological signal containing at least two wave cycles originating from a physiological process,
obtaining a feature signal descriptive of occurrences of a feature in the physiological signal,
determining an average waveform of the physiological signal around said occurrences of the feature in the physiological signal,
determining a model signal comprising amplitude-scaled instances of the average waveform of the physiological signal placed at the occurrences of the feature in the physiological signal
In a further aspect of the present invention a system for processing a physiological signal is presented comprising:
a sensor for measuring a physiological signal, and
a processing unit comprising a first processing apparatus for determining a feature signal descriptive of occurrences of a feature in the physiological signal, and a second processing apparatus for processing a physiological signal, wherein the second processing apparatus is a processing apparatus as described above configured to determine a model signal,
wherein the processing unit is further configured to determine the vital sign of the subject based on the physiological signal taking into account the model signal.
In yet further aspects of the present invention there are provided a corresponding processing method, a computer program which comprises program code means for causing a computer to perform the steps of the method disclosed herein when said computer program is carried out on a computer as well as a non-transitory computer-readable recording medium that stores therein a computer program product, which, when executed by a processor, causes the method disclosed herein to be performed.
Preferred embodiments of the invention are defined in the dependent claims. It shall be understood that the claimed method, system, computer program and recording or storage medium have similar and/or identical preferred embodiments as the claimed processing apparatus and as defined in the dependent claims.
Determining a vital sign such as a heart rate or respiratory rate from a physiological signal is far from trivial, in particular if the physiological signal is obtained via an unobtrusive sensor. The measured physiological signal often comprises signal distortions such as drifts, noise or movement artifacts. Furthermore, the measured physiological signal may comprise information about more than one underlying physiological process such as information about cardiac activity and respiration. For example, a constant heart rate and a constant respiration rate show as distinct peaks in a spectrum of the physiological signal. They can be separated from each other by filtering in frequency-domain. Correspondingly, drifts, noise and artifacts can be suppressed by filtering the spectrum.
While these conventional approaches work well in many applications, there are several limitations: Evaluating a physiological signal by spectral analysis assumes an at least weak stationary signal. However, in case of arrhythmias, i.e., variations, of the underlying physiological process or processes, such approaches do not work anymore. Furthermore, information about the curve shape or morphology is lost. For example, a heart rate is extracted from the spectrum but the shape of the physiological signal corresponding to individual heart beats is distorted by the filtering process.
Furthermore, if the spectra of two underlying physiological phenomena overlap, filtering the physiological signal to extract a first signal, corresponding to a first physiological phenomenon, may impact a spectrum of a second signal, corresponding to a second physiological phenomenon. The signal shape or morphology of the first and second signal should be preserved.
In a further aspect, vital signs can be derived from characteristic features of the physiological signal. However, one has to make sure the features are extracted correctly. This can be indicated by a quality metric. Alternatively, the quality metric can be used for classification of the physiological signal as being good or bad. For example, if the signal feature is a signal peak, a classifier output indicates whether peak-detections from a peak-detector contain errors, e.g., false-positives or false-negatives. Current approaches for determining a quality metric based on heart-rate variability or spectral analysis assume an at least weak stationary signal. In case of arrhythmias, i.e., large variations in peak-to-peak distances, such approaches do not work anymore.
Furthermore, the inventors have found that a shape and spectrum of the physiological signal largely depends on the signal source from which the physiological signal is obtained. For example, a cardiac pulse in an ECG signal differs considerably from new measurement techniques, for example wherein an accelerometer is used for measuring a physiological signal. In particular with new measurement techniques, the curve shape corresponding to a specific physiological phenomenon such as a cardiac pulse can vary significantly from subject to subject.
Furthermore, the waveform of a physiological signal corresponding to a physiologic phenomenon can change depending on the posture of the subject. Both the curve shape and amplitude can be affected. Thus, a plurality of different waveforms can represent one and the same physiological phenomenon which renders a generic pattern comparison difficult.
However, the inventors have found that despite the variation between patients and postures, for one patient, in a given posture, the waveforms corresponding to a physiological phenomenon are very similar. This also holds true for the case of arrhythmias.
Therefore, according to the present invention, a processing apparatus is proposed that takes the physiological signal as well as a feature signal descriptive of occurrences of a feature in the physiological signal as inputs and determines an average waveform around the occurrences. The physiological signal contains at least two wave cycles originating from a physiological process to allow averaging. The average waveform can be invariant to timing and amplitude variations. Based on this average waveform a model signal is determined that comprises amplitude-scaled instances of the average waveform of the physiological signal at occurrences of the feature in the physiological signal. In other words, a copy of the average waveform of the physiological signal around occurrences of the feature in the physiological signal is placed at occurrences of the feature. The amplitudes of each copy of the average waveform can be scaled depending on an amplitude of the physiological signal at or around the occurrence of the feature where the copy is placed. The plurality of copies of the average waveform makes up the model signal.
An advantage of this approach is that it is not necessary that the waveform of physiological phenomena in the physiological signal is known a priori. Furthermore, the processing apparatus can be used for evaluating a non-stationary physiological signal. Thus, also non-periodic physiological signals can be processed. Moreover, the shape of the model signal is not manipulated by frequency-selective filters.
Referring to the system for determining a vital sign of a subject, the processing unit can be configured to determine the vital sign of the subject based on the physiological signal taking into account the model signal. In a first aspect, the model signal is taken into account in that the model signal is used to determine a quality metric of the physiological signal and/or the feature signal. For example, the quality metric can be provided along with the vital sign. In a second aspect, the model signal is taken into account in that the vital sign is determined from the model signal. The model signal in turn is based on the physiological signal. Alternatively, the vital sign is determined based on a difference of physiological signal and model signal as will be described in detail further below.
In a preferred embodiment of the processing apparatus, the physiological signal is a movement signal descriptive of a physiological movement. The processing apparatus can be configured for obtaining a movement signal descriptive of a physiological movement, for example, from an accelerometer or a gyroscope. The proposed processing apparatus is particular advantageous for processing movement signals. A problem of movement signals is that the waveform largely depends on a placement of the sensor and may change significantly depending on the posture of the subject. However, since the processing apparatus according to the present invention determines an average waveform of the physiological signal around the occurrences of the feature in the physiological signal, the shape does not have to be known beforehand.
In a further embodiment, the feature of the physiological signal is a peak of the physiological signal. A peak in this context refers to a local maximum or minimum of the physiological signal. In this embodiment, the feature signal indicates the locations of peaks in the physiological signal. Alternatively, a feature refers to a part of the waveform comprising a predetermined energy, slope, inflection points, local extrema or combinations of the above. In general, any characteristic of the physiological signal can be evaluated as a feature representing an underlying physiological phenomenon.
The resulting model signal can directly be evaluated or used for further processing. For example, the model signal can serve as a “clean” physiological signal. In an alternative embodiment, the processing apparatus is further configured to perform the step of comparing the model signal with the physiological signal.
In an embodiment, the comparing step further comprises the computation of a difference signal between the model signal and the physiological signal. For example, the model signal, corresponding to a first physiological phenomenon, can be removed from the physiological signal which gives a residual physiological signal. An advantage of this embodiment is that the waveform of further contributions to the physiological signal remains untouched. This is in contrast to other filtering techniques (like low-pass filtering, median filtering or Savitsky-Golay filtering) where the residual physiological signal is modified by the filtering process.
In a further refinement, the comparing step further comprises the computation of an energy of the difference signal descriptive of the quality of the feature signal or the physiological signal. The term energy in this context refers to a value descriptive of the difference signal. Optionally, the energy can be normalized by, for example, a value of the physiological signal. For example, the energy represents the value of an integral over a segment of the difference signal. A high energy results from a large difference between physiological signal and model signal and thus represents a large mismatch. A large mismatch can be interpreted as a poor signal quality whereas a small mismatch can be interpreted as a good signal quality. Alternatively, the energy can be compared with a threshold to determine a digital value, for example, indicating 1 if the difference is below a threshold value for high signal quality and 0 if the difference exceeds the threshold. In other words, the energy of the difference signal is inversely proportional to the quality of the feature signal or the physiological signal. The energy value can thus be used for determining a quality metric descriptive of a quality of the feature signal and/or the model signal and/or the physiological signal.
In an embodiment, the step of determining the average waveform comprises the step of determining an average relative change between a value of the physiological signal at occurrences of the feature and a value of a neighboring sample of the physiological signal. Thereby, the analysis is independent of predetermined patterns descriptive of physiological features but purely relies on the actually measured physiological signal. For example, at each occurrence of the feature (as specified by the feature signal) the value of the physiological signal is taken and a value of a neighboring sample is subtracted therefrom. The result can be divided by a value of the physiological signal, in particular the value of the physiological signal at the occurrence of the feature or the value of the neighboring sample, to obtain a relative change. For example, for each peak of the physiological signal the value of the sample immediately preceding the respective peak is subtracted and divided by the respective peak value. Thereby, a relative change for this particular neighbor is determined. The relative changes for the neighbors, wherein the neighbors refer to the neighbors immediately preceding the respective peaks, are averaged to determine an average relative change of the immediately preceding sample with respect to a peak.
In other words, the step of determining the average waveform comprises the scaling of the difference between the value of the physiological signal at an occurrence of the feature and the value of a neighboring sample. Thereby the differences for determining the average relative change can be weighted, for example, by the value of the physiological signal at an occurrence of the feature or an inverse thereof. In the latter case, the impact of samples with high values of the average waveform is reduced.
In a further refinement of this embodiment, the step of determining the average waveform comprises repeating the step of determining an average relative change for different neighboring samples to a right and/or a left side, i.e., preceding and succeeding samples, of occurrences of the feature. Thus, not only a relative change for the immediate neighbor is determined but relative changes for samples spaced further apart with respect to the occurrence of the feature to a left and/or right side of the feature are determined and thereby construct the average waveform. In other words, in this embodiment, the average waveform represents the average relative change at predetermined distances from the occurrences of the feature in the physiological signal.
In a further refinement of this embodiment, said repeating is performed until reaching a stop criterion, in particular until a relative change has been determined for all samples of the physiological signal between at least one pair of succeeding occurrences of the features having the smallest temporal separation.
In an embodiment, the step of determining the model signal further comprises scaling an instance of the average waveform depending on a value of the physiological signal at the occurrence of the feature in the physiological signal. Not only the temporal separation of features in the physiological signal can vary, but also the amplitude. Thus, the instance of the average waveform at a particular occurrence of the feature in the physiological signal can be scaled to match an amplitude of the physiological signal. For example, the instance of the average waveform can be scaled such that a value of the average waveform at the occurrence of the feature corresponds to a value of the physiological signal at the occurrence of the feature.
In an embodiment, the step of determining the model signal further comprises the step of interpolating the model signal between instances of the average waveform. As indicated above, there can be a stop criterion for determining the average waveform which is based on the smallest temporal separation of two occurrences of the physiological signal. Thus, there can be a gap between two instances of the average waveform in the model signal which can be filled by interpolating values in-between. Any type of interpolation can be used in particular linear cubic and spline interpolation.
In an embodiment, the processing apparatus is further configured to evaluate a variation of the temporal separation of occurrences of a feature in the feature signal. Advantageously, the proposed determination of a quality metric based on an energy of the difference signal as disclosed herein, can be combined with alternative approaches for determining a quality of the feature signal and/or the physiological signal. However, even if a large variation of the temporal separation is detected, for example due to a non-life threatening arrhythmia, the signal can still be classified as good if the quality metric based on the energy of the difference signal indicates so.
In yet another embodiment, the processing apparatus is further configured to evaluate the physiological signal and/or feature signal in time-domain. Advantageously, blocks, i.e., segments of the physiological signal and corresponding feature signal, of limited length are evaluated. Advantageously, blocks of a duration of 1 to 8 seconds are evaluated. This is an advantage compared to frequency domain processing, wherein typically longer windows of, for example, 30 seconds are evaluated for meaningful results in frequency domain.
These and other aspects of the invention will be apparent from and elucidated with reference to the embodiment(s) described hereinafter. In the following drawings
The device 1 for determining a vital sign of a subject 100 according to an aspect of the present invention enables an unobtrusive measurement of, for example, a movement signal descriptive of a respiratory movement and/or a movement due to cardiac pulses. The signals form the basis for determining a respiratory rate and a heart rate (pulse rate) of the subject. As can be seen in
Conventional respiration monitoring devices, for example, a nasal canula are inconvenient for the patient. Alternative conventional respiration monitoring devices such as a chest belt for respiratory inductance plethysmography require more effort in handling. For respiratory inductance plethysmography (RIP) one or two RIP-bands have to be wrapped around the chest and/or abdomen of the subject 100. Especially for a heavy patient, the application of these bands can be cumbersome for a nurse. Nonetheless, respiratory inductance plethysmography is one of the most common techniques used for determining a respiration rate today. Further alternative conventional respiration monitoring devices include devices for thorax impedance plethysmography that are used in intensive care units. In thorax impedance plethysmography, the impedance is measured via two electrodes, which can be inconvenient for the patient. Correspondingly, for conventional heart rate measurements either ECG electrodes or a finger clip for photoplethysmography are applied to the patient.
The sensor 5 may, for example, comprise an accelerometer which is adapted to generate a movement signal indicative of the acceleration along one or more spatial axes. In this embodiment, the accelerometer is a tri-axial accelerometer adapted to generate a movement signal that comprises three accelerometer signals indicative of the acceleration along three orthogonal spatial axes. For example, three-axial accelerometers named Bosch BMA355, ST Microelectronics LIS3DSH, ST Microelectronics LIS344ALH or Kionix KXM52 can be used. Alternatively, the sensor 5 may comprise a gyroscope which is adapted to generate a movement signal indicative of an angular rotation. For example, the gyroscope ST Microelectronics LYPR540AH can be used to acquire angular rates for three orthogonal axes from which angular rotations can easily be derived.
The movement signal 6 is provided to the processing unit 7 for determining a vital sign of the subject. An exemplary embodiment of a processing unit 7 will be described further below with reference to
The interface 10 can be a wireless interface for providing a vital sign signal of the subject. Advantageously, the device 1 is battery powered device and the interface 10 is a wireless interface such that no cables are required. For example, the interface 10 can provide the vital sign to a patient monitor for displaying the vital sign on a screen. Alternatively, the interface 10 can be a human machine interface (HMI) for showing the vital sign and/or the physiological signal and/or any quantity derived thereof on a display of the device 1.
The memory 11 can store the physiological signal and/or the feature signal derived therefrom and/or the determined vital signs that have been determined by the processing unit 7. Alternatively, the memory 11 can be a non-transitory storage medium containing instructions for execution by the processing unit 7 or components thereof, wherein the instructions cause the processing unit 7 to perform at least some of the steps of the flowchart shown in
In an alternative embodiment, the sensor 5 and the processing unit 7 or parts thereof are not implemented in one device. For example, a minimum configuration of sensor 5 and interface 10 can be implemented as a device 1 that is worn by the patient. The physiological signal 6 is thus transmitted via the interface 10 to a processing unit 7 at a remote location. For example, the processing unit 7 can be implemented as a part of a patient monitor or part of a health care infrastructure such as a hospital IT system. In a further embodiment, a smartphone or similar device serves as the processing unit 7 and is configured to perform the steps for determining a vital sign of the subject and/or to perform the steps for processing a physiological signal according to an aspect of the present invention.
The second physiological signal 6b has been obtained using an alternative measurement technique, i.e. respiratory inductance plethysmography, which is also often referred to as respiration band or respiband. The second feature signal 12b, indicated by crosses, denotes peaks of the second physiological signal 6b. It should be noted that the respiband physiological signal 6b shows a very clean signal. Thus, the features of the feature signal 12b can be derived from the physiological signal 6b with high accuracy and the feature signal 6b has an inherently high quality. This respiband physiological signal 6b serves for comparison and for evaluating the performance of the accelerometer-based physiological signal 6a.
The dashed vertical lines in
Both physiological signals 6c and 6d are obtained from the same subject. The time delay between the two signals is indicated by vertical dashed lines, wherein the accelerometer-based physiological signal 6c is slightly delayed with respect to the ECG-based physiological signal 6d.
It should be noted that the subject in this example has a non-steady, irregular heart rate that shows arrhythmias. Just after second 2400 and after second 2410, there are wide intervals between two heart beats (inter-beat-interval, IBI). In this example, a classifier or quality metric that is solely based on the peak-to-peak variances, i.e. a classifier that evaluates how the separation of neighboring peaks fluctuates may classify the signal as being bad. Furthermore, a classifier that evaluates spectral features may fail and classify the signal as being bad because arrhythmias cause distortions to the spectrum of the measured physiological signal. For example, multiple peaks may be present in the spectrum or the spectrum around a peak may be significantly broadened.
A classifier is usually determined for segments of the physiological signal. If the classifier indicates a poor signal quality, the entire segment is discarded and not evaluated. As a consequence a reduced number of segments of the physiological signal is evaluated and the coverage of the physiological signal is lowered.
In
In this embodiment, the pre-processing apparatus 13 receives the output signals acc_x, acc_y, acc_z of a tri-axial accelerometer and combines these signals into a physiological signal 6 that is provided to the first processing apparatus 8 and the second processing apparatus 9. The pre-processing apparatus 13 is an optional element that can be part of a sensor 5 (as shown in
In the embodiment described herein, the first processing apparatus 8 for determining a feature signal descriptive of occurrences of the feature in the physiological signal 6 is a peak-detector. However, it should be noted that the invention is not limited in this respect and that other features of the signal such as extrema of the wave, extrema of a gradient of the wave, i.e. a slope of the signal, inflection points, a peak energy and the like can also be considered as features. The peak-detector (first processing apparatus 8) provides a feature signal 12 to the second processing apparatus 9. The features signal 12 is descriptive of occurrences of peaks in the physiological signal. However, the peak-detection is not necessarily perfect. Sometimes a peak is missed (FN: false negative) and sometimes peaks are wrongly detected (FP: false positive). Thus, the second processing apparatus 9 in this embodiment is configured to classify whether the peak-detection is accurate.
The second processing apparatus 9 in this embodiment can also be referred to as a classifier unit. In an embodiment, the classifier unit computes a quality metric descriptive of a quality of the feature signal 12 based on a window, i.e. a section, of the physiological signal 6 of, for example 8 seconds length, and classifies the physiological signal 6 and/or the feature signal 12 as either being good or bad. The classifier result is provided at an output 14 of the classifier unit 9.
Optionally, the classifier unit 9 also provides a feedback path 15 to the peak-detector 8. Thereby, the peak-detector 8 receives information whether a segment of the physiological signal 6 is good and should be used for determining a vital sign signal 16 or whether this particular segment should be discarded and not used for determining a vital sign of a subject. In this exemplary embodiment, the vital sign of the subject is a heart rate. However, the invention is not limited in this respect and could also be used for different physiological signals and vital signs, for example, for processing a respiratory signal as the physiological signal 6 and for determining a respiratory rate as the vital sign.
The processing of a physiological signal 6 by the processing apparatus 9 is described in more detail with reference to the flowchart shown in
In a second step S12, a feature signal descriptive of occurrences of a feature in the physiological signal is obtained. For example, the feature signal can indicate the time or locations such as the samples of the physiological signal, where a feature such as a peak occurs. For example, the feature signal can be obtained from a peak-detector 8 as shown in
In a third step S13 an average waveform of the physiological signal around said occurrences of the feature in the physiological signal is determined. For example, the peaks of the physiological signal are taken as a starting point and for each peak, the difference of its neighboring values with respect to the value of the peak are determined. Based on these differences, for each distance with respect to a peak, an average difference is determined. Optionally, the difference between a peak and neighboring value is weighted, for example, weighted by the height of the corresponding peak or the inverse thereof to obtain a relative change.
In a next step S14, a model signal comprising instances of the average waveform of the physiological signal at the occurrences of the feature in the physiological signal is determined. In other words, for each peak a copy of the average waveform is instantiated and arranged such that the peak of the average waveform corresponds to the respective peak of the physiological signal. Optionally, the height of the average waveform is scaled by the height of the respective peak such that the model signal comprises amplitude-scaled instances of the average waveform.
In a last optional step S15, a quality metric descriptive of a quality of the feature signal is determined based on a comparison of the physiological signal and the model signal. For example, a difference signal between the physiological signal and the model signal is computed. Optionally, the difference signal is normalized, for example by peak values or an average value of the physiological signal. An energy of the difference signal is computed which is descriptive of the quality of the feature signal and/or the physiological signal. A small difference and thus small energy of the difference signal indicates a good matching and high signal quality whereas a high difference indicates poor matching and low signal quality. The output 14 of this step is a classifier output based on the quality metric. This quality metric is computed as, for example directly an energy value or normalized energy value of the difference signal. Alternatively, the computed difference is compared with a threshold to determine a digital output, for example 1 if the signal quality is good and the respective segment of the physiological signal should be used for determining a vital sign, or 0 which means that the segment should be discarded. In an embodiment, a classifier determines whether a particular excerpt belongs to a specific, distinct class. The output of the classifier can be a label or identifier of one of these classes
It should be noted that the steps of determining an average waveform of the physiological signal around occurrences of the feature in the physiological signal S13 and a step S14 of determining a reference signal can also be performed in a combined step or also in an iterative manner. For example, not the entire average waveform is determined before determining the model signal. Instead, the average waveform and the model signal can also be determined step by step on a sample-by-sample basis.
An exemplary embodiment of determining a model signal based on averaging the physiological signal around occurrences of the feature in the physiological signal is illustrated with reference to
Referring to
In a second step S22, for every neighbor value of the peaks, a relative change with respect to the peak_value is determined by calculating change(p)=(peak_value(p)−neighbor_value(p))*peak_value_inverse(p), where p indicates a peak-number index p=[0 . . . P−1]. A neighbor_value as used herein refers to a sample that is located at a left or right side with respect to a peak. For example, the first neighbor to the left is the first sample to the left in
In step S23, an average relative change of a particular neighbor value is determined for all peaks, for example by calculating Avg_change=sum(change(p))/P. For example, for the first left neighbor for each peak the relative change of the first left neighbor to the respective peak is calculated. Then the sum of all these relative changes is calculated and the result divided by the number of peaks.
In a next step S24, the auxiliary vector m_est(p) is updated based on the value of the peaks peak_value(p) and the determined average relative change Avg_change by calculating m_est(p)=(1−Avg_change(p))*peak_value(p). The multiplication by the peak value gives an amplitude scaled value of the average relative change. This auxiliary vector is used for updating the reference signal in that, for example, the first left neighbor of the first peak is assigned the value of the first element of the auxiliary vector m_est(p), the first left neighbor of the second peak is assigned the second value of the auxiliary vector m_est(p), and so on.
The same calculations can be performed for further left neighbors and further right neighbors. This is indicated by a loop 17 in
In a last optional step S26, sample values of the reference signal that have not been calculated in the afore-mentioned process can be interpolated to close gaps. An exemplary gap, indicated by X2, occurs between peaks P8 and P9.
The procedure according to this aspect of the present invention is illustrated in with reference to the graphs shown in
Finally, as described with reference to step S26 in
Optionally, a classifier output is computed based on a quality metric. This quality metric is descriptive of a quality of the feature signal, calculated based on a comparison of the physiological signal and the reference signal. The quality metric can be computed as an energy of a difference signal between the model signal 18 and the physiological signal 6c. To evaluate how well the model signal matches the physiological signal, the quality metric can be computed, for example, by energy=sum(abs(physiological_signal−model_signal))/sum(abs(physiological_signal)). If the model signal is not computed for all samples of the physiological_signal, the computation of the quality metric or energy can be limited to those samples where a model signal is available. The normalization by the absolute values of the physiological signal is optional but can help in comparing segments with different signal magnitude. For the example shown in
The situation for a correct detection, i.e. wherein the peak P2 at sample-index 230, instead of peak P1 at sample-index 200, is detected, is shown in
It should be noted that in case of large movement artifacts, the value of the quality metric or energy of the difference signal can exhibit even more extreme values, for example 10 or 100. However, the value of the quality metric is not limited in this respect and depends on the desired normalization scaling and amplitude of the evaluated physiological and reference signal.
A more detailed example of how a quality metric can be determined including details of an exemplary calculation of average waveform and model signal is described with reference to
The proposed signal processing is based on the finding that each wave cycle in the physiological_signal has a similar shape, however, not necessarily a certain (i.e. regular) temporal separation and/or amplitude of the wave cycles. For example, amplitudes of a heart-wave can depend on the breathing.
The basic steps of determining the model signal in this example can be summarized as (1) finding the peaks as starting points, (2) for each (left and/or right) neighboring sample: compute the average relative drop/increase in level compared to the initial-peaks, and compute the values of the model signal based on the average drop/increase and initial-peaks, and (3) a (linear) interpolation of further values of the model signal, when any two neighboring samples coincide.
For describing the further algorithm hereafter, the physiological_signal is denoted as a vector xpP[wP] with wP being the index of the 8 second window. The length of x equals Nx. In the example in
The waveform values at the detected peaks-locations are denoted by:
(xp)j(x)(p)
where j=0, . . . , Np−1 is the index of the peak. In the example of
The next step in determining the average waveform and determining the model is to look for neighboring samples (toward the left and right with respect to the initial peaks) and to compute the average drop or increase compared to the initial peak-levels.
The neighboring locations with respect to the peak-locations can be defined as:
p
Δk
=p+Δ
k.
where k=1 defines the scanning of neighboring samples to the left, whereas k=2 defines the scanning of neighboring samples to the right, thus
Δk<0 for: k=1
Δk>0 for: k=2.
For all neighboring locations with respect to. the peak-locations which are inside the 8 second window, the average level-drop/increase is computed as:
where N equals the number of neighboring locations that are inside the region [0, . . . , Nx−1] computed as:
N=#{∀
j=0 . . . Np-1:0≦└(pj┘<Nk}.
Thereby, the quantity ψ for the neighbors of the respective peaks gives the average waveform of the physiological signal around the occurrences of the peaks. In the given example, the drops/increases are scaled by the values (xp)j of the peaks.
The model signal, for comparison with the physiological signal x, is denoted by m. The length of m is equal to the length of x and is denoted by Nm=Nx.
The values of the model signal are computed as:
∀j:0≦└(p)j+Δk┘<Nx:(m)j=(xp)j·(1−ψΔk).
In
An exemplary calculation of the model signal involves iteratively decreasing values of Δ1 and increases the values of Δ2, i.e. this gives a sequence of (k, Δk)=(1,−1), (2, 1), (1,−2), (2, 2), . . . etc. After a limited number of iterations, obviously two indices will coincide and the computation of the reference waveform can be stopped:
∃j:0≦j<N
In this non-limiting example, the stopping occurs when:
The term └(p)j+1−(p)j┘ essentially gives the interval between two neighboring peaks. The values of Δk with k=1, 2 at the moment that the condition of the previous equation occurs will be denoted by Δk* with k=1, 2. These values denote the offsets (toward the left and right) that represent the boundaries of the model for each wave (centered around the occurrences of peak-detections). For the given example, the values of Δ1 and Δ2 are respectively −36 and 36. The dashed curve in
Alternatively, a different stopping criterion can be employed for left and/or right side. For example, the determination of values of the model signal towards the left-side (Δ1<0) ends earlier than toward the right-side (Δ2>0). For example, the stopping criterion to the left-side can be defined as:
∀j:(xp)j·(1−ψΔ1)<0.
The two stopping criteria described above can also be combined or alternative stopping criteria can be used.
After the aforementioned determination of values of the model signal m, there can still be unmodelled values in m. In other words, not all values of m are filled. This can also be seen in
The unmodeled values in-between the peak-detections of the segment can be filled by means of, for example, linear interpolation. For each pair subsequent peak-detections (Np−1 pairs), a vector g having Np−1 elements can be computed. The elements of the vector indicate the lengths of the respective gaps in samples. For the i'th pair the i'th element of the vector can be computed by
∀0≦i<(N
where Δ1* and Δ2* are defined as being the offsets (toward the left and right respectively) that represent the samples at boundaries of the model signal for each scaled average waveform (centered around the peak-detections). The elements in the vector g for the example shown in
The gaps in the model m are then filled up by computing:
Optionally, a quality metric descriptive of a quality of the feature signal can be determined by computing a normalized magnitude difference between the original data x and the model m. Since the model signal may not have been determined for all samples of the physiological signal, the calculation of the quality metric may be limited to those values where the model signal exists. For example if samples at a boundary of the segment have not been modeled, the determination of the quality metric may be limited to the range between a leftmost sample and a rightmost sample of the model signal, hence:
The quality metric as calculated above can be directly used for indicating a quality of the feature signal or can be used in a classifier for evaluating the physiological signal. For example, a classifier-tree can be used:
Advantageously, the energy of the physiological signal should be below a certain bound (bound1) in order to avoid classifications in case of heavy movements measured by the accelerometer, e.g. periodic walking signals.
The peak-to-peak variance can be used for the classification in case of regular signals, i.e. signals without arrhythmia, whereas the classifier according to an aspect of the present invention can be used in case of arrhythmias.
A further advantageous application of the processing apparatus 9 according to an aspect of the present invention is described with reference to
The system shown in
The processing apparatus 9 obtains the physiological signal 6, containing at least two wave cycles originating from cardiac activity (ECG signal) and respiration (EMG signal), and the feature signal 12 from the feature detector 8. The model signal 18 is provided at an output of the processing apparatus 9. In this embodiment, the feature signal 12 indicates occurrences of cardiac pulses. Therefore, the resulting model signal 18 comprises instances of the average waveform of the physiological signal 6 at the occurrences of the cardiac pules in the physiological signal. Thus, the model signal 18 represents a model of the ECG signal.
The model signal 18 can be subtracted from the physiological signal 6 with the adder 19, wherein the model signal is provided at an inverted input of the adder 19. The output signal is a difference signal or residual signal 20 which is effectively free from the ECG signal. The difference signal can be evaluated, for example to determine the respiration rate or the respiratory effort. The model signal 18 can also be used for further processing or for example be presented to the subject or physician as a clean ECG signal.
The ECG waves due to cardiac pulses are cancelled in the residual signal 20. This can also be seen from the spectrum shown in the lower plot of
In summary, the present invention proposes a processing apparatus, processing method and device for processing a physiological signal, based on an obtained physiological signal and an obtained feature signal descriptive of occurrences of a feature in the physiological signal, by determining an average waveform of the physiological signal around said occurrences of the feature in the physiological signal and by determining a model signal comprising amplitude-scaled instances of the average waveform of the physiological signal at the occurrences of the feature in the physiological signal. Thus, the proposed processing apparatus provides a model signal which can, for example, be used for determining a quality of the feature signal (which forms the basis for determining a vital sign of the subject) even in case of arrhythmias and also in case of fluctuating amplitudes. Furthermore, the model signal can serve for distinguishing different contributions to a physiological signal.
While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive; the invention is not limited to the disclosed embodiments. Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims.
In the claims, the word “comprising” does not exclude other elements or steps, and the indefinite article “a” or “an” does not exclude a plurality. A single element or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.
A computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.
Any reference signs in the claims should not be construed as limiting the scope.
Number | Date | Country | Kind |
---|---|---|---|
13186423.3 | Sep 2013 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2014/069768 | 9/17/2014 | WO | 00 |