CN109864705B - Method and device for filtering pulse wave and computer equipment - Google Patents
Method and device for filtering pulse wave and computer equipment Download PDFInfo
- Publication number
- CN109864705B CN109864705B CN201910012502.9A CN201910012502A CN109864705B CN 109864705 B CN109864705 B CN 109864705B CN 201910012502 A CN201910012502 A CN 201910012502A CN 109864705 B CN109864705 B CN 109864705B
- Authority
- CN
- China
- Prior art keywords
- peak
- pulse wave
- height difference
- extreme
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The application provides a method, a device and computer equipment for filtering pulse waves, wherein the method for filtering the pulse waves comprises the following steps: acquiring the acquired pulse wave; carrying out extreme value detection on the pulse wave to obtain an extreme value point of the pulse wave; classifying the extreme points of the pulse waves to obtain peaks and troughs of the pulse waves; for each peak of the pulse wave, respectively calculating the height difference between each peak and a first trough after the peak, and calculating the height difference between each peak and a second trough after the peak; averaging the height differences obtained by calculation to obtain a first average value of the height differences, and setting a first threshold value according to the first average value of the height differences; and filtering the extreme points of the pulse wave, wherein the absolute value of the amplitude value of the extreme points is greater than or equal to the first threshold value. This application can realize carrying out the filtering to the unusual wave form in the pulse ripples, improves the accuracy to the crest and the trough location of pulse ripples.
Description
[ technical field ] A method for producing a semiconductor device
The present application relates to the field of signal processing technologies, and in particular, to a method and an apparatus for filtering a pulse wave, and a computer device.
[ background of the invention ]
Pulse, respiration and the like are important vital signs of a human body, and the intensity, the form and the frequency of the vital signs can reflect the physiological and pathological information of the human body, such as the physical state, the mental state, the health level and the like.
The extraction of physiological and pathological information of human body from pulse wave as the basis of clinical diagnosis and treatment has been regarded by the medical field, and some abnormal waveforms are usually generated during the pulse wave collection process due to unskilled user operation or physical quality, such as: the dicrotic wave is at the same level as the main wave, or has an abnormally high value or an abnormally low value due to external factors.
In the prior art, a difference value or a sliding window method is generally adopted to extract peaks and troughs of pulse waves, but the two methods cannot identify the abnormal waveforms in the pulse waves and even remove the abnormal waveforms, so that the peaks and the troughs of the pulse waves are not accurately positioned.
[ summary of the invention ]
The embodiment of the application provides a method, a device and computer equipment for filtering pulse waves, so that abnormal waveforms in the pulse waves are filtered, and the accuracy of positioning wave crests and wave troughs of the pulse waves is improved.
In a first aspect, an embodiment of the present application provides a method for filtering a pulse wave, including: acquiring the acquired pulse wave; carrying out extreme value detection on the pulse wave to obtain an extreme value point of the pulse wave; classifying the extreme points of the pulse waves to obtain peaks and troughs of the pulse waves; for each peak of the pulse wave, respectively calculating the height difference between each peak and a first trough after the peak, and calculating the height difference between each peak and a second trough after the peak; averaging the height differences obtained by calculation to obtain a first average value of the height differences, and setting a first threshold value according to the first average value of the height differences; and filtering the extreme points of the pulse wave, wherein the absolute value of the amplitude value of the extreme points is greater than or equal to the first threshold value.
In a possible implementation manner, after filtering the extreme points of the pulse wave whose absolute value of the amplitude is greater than or equal to the first threshold, the method further includes: classifying the residual extreme points of the pulse waves to obtain wave crests and wave troughs in the residual extreme points of the pulse waves; for each peak in the remaining extreme points, respectively calculating the height difference between each peak and a first trough after the peak, and calculating the height difference between each peak and a second trough after the peak; averaging the height difference obtained by calculation to obtain a second average value of the height difference, and setting a second threshold value according to the second average value; and filtering the extreme points of which the absolute value of the amplitude values in the residual extreme points of the pulse wave is smaller than the second threshold value.
In a possible implementation manner, the setting a first threshold according to the first mean value of the height differences includes: multiplying the first mean value of the height difference by a predetermined coefficient to obtain a product as the first threshold value.
In a possible implementation manner, before performing extremum detection on the pulse wave and obtaining an extremum point of the pulse wave, the method further includes: and carrying out noise reduction processing on the acquired pulse wave through polynomial curve fitting.
In a second aspect, an embodiment of the present application provides an apparatus for filtering a pulse wave, including: the acquisition module is used for acquiring the acquired pulse waves; the detection module is used for carrying out extreme value detection on the pulse waves acquired by the acquisition module to acquire extreme value points of the pulse waves; the classification module is used for classifying the extreme points of the pulse waves obtained by the detection module to obtain the wave crests and the wave troughs of the pulse waves; the calculation module is used for calculating the height difference between each peak and a first trough behind the peak and the height difference between each peak and a second trough behind the peak respectively for each peak of the pulse waves; the filtering module is used for averaging the height difference obtained by calculation to obtain a first average value of the height difference and setting a first threshold value according to the first average value of the height difference; and filtering the extreme points of the pulse wave, wherein the absolute value of the amplitude value of the extreme points is greater than or equal to the first threshold value.
In a possible implementation manner, the classification module is further configured to classify remaining extreme points of the pulse wave after the filtering module filters the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold from the extreme points of the pulse wave, so as to obtain peaks and troughs in the remaining extreme points of the pulse wave; the calculation module is further configured to calculate, for each peak in the remaining extreme points, a height difference between each peak and a first trough after the peak, and a height difference between each peak and a second trough after the peak, respectively; averaging the height difference obtained by calculation to obtain a second average value of the height difference, and setting a second threshold value according to the second average value; the filtering module is further configured to filter the extreme point, where an absolute value of the amplitude in the remaining extreme points of the pulse wave is smaller than the second threshold.
In a possible implementation manner, the filtering module is specifically configured to multiply the first mean value of the height difference by a predetermined coefficient, and use the obtained product as the first threshold.
In a possible implementation manner, the apparatus for filtering a pulse wave further includes: and the noise reduction module is used for carrying out noise reduction processing on the collected pulse waves through polynomial curve fitting before the detection module carries out extreme value detection on the pulse waves and obtains extreme values of the pulse waves.
In a third aspect, an embodiment of the present application provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor executes the computer program to implement the method described above.
In a fourth aspect, embodiments of the present application provide a non-transitory computer-readable storage medium having stored thereon a computer program, which when executed by a processor, implements the method as described above.
In the above technical solution, after the acquired pulse wave is acquired, the extreme value detection is performed on the pulse wave to obtain the extreme value point of the pulse wave, classifying the extreme points of the pulse wave to obtain peaks and troughs of the pulse wave, calculating a height difference between each peak and a first trough after the peak and a height difference between each peak and a second trough after the peak for each peak of the pulse wave, then averaging the height differences obtained by calculation to obtain a first average value of the height differences, setting a first threshold value according to the first average value of the height differences, filtering the extreme points of which the absolute value of the amplitude value is greater than or equal to the first threshold value in the extreme points of the pulse wave, therefore, the abnormal waveform in the pulse wave can be filtered, and the accuracy of positioning the wave crest and the wave trough of the pulse wave is improved.
[ description of the drawings ]
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a flowchart illustrating an embodiment of a method for filtering a pulse wave according to the present application;
FIG. 2 is a flowchart illustrating another embodiment of a method for filtering a pulse wave according to the present application;
FIG. 3 is a flowchart illustrating a method for filtering a pulse wave according to another embodiment of the present invention;
FIG. 4 is a schematic structural diagram of an embodiment of an apparatus for filtering a pulse wave according to the present application;
FIG. 5 is a schematic structural diagram of an embodiment of a computer apparatus according to the present application.
[ detailed description ] embodiments
For better understanding of the technical solutions of the present application, the following detailed descriptions of the embodiments of the present application are provided with reference to the accompanying drawings.
It should be understood that the embodiments described are only a few embodiments of the present application, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The terminology used in the embodiments of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in the examples of this application and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
Fig. 1 is a flowchart illustrating an embodiment of a method for filtering a pulse wave according to the present application, and as shown in fig. 1, the method for filtering a pulse wave may include:
In a specific implementation, the argrelmax algorithm in signal detection may be used to perform extreme value detection on the pulse wave, and when the argrelmax algorithm is used to detect the extreme point of the pulse wave, the principle of minimum distance is used.
Of course, other extreme value detection methods may be used to detect the extreme value point of the pulse wave, which is not limited in this embodiment.
Specifically, after obtaining the extreme points of the pulse wave, the extreme points of the pulse wave may be classified to distinguish peaks and troughs of the extreme points of the pulse wave.
And 104, respectively calculating the height difference between each peak and a first trough after the peak and the height difference between each peak and a second trough after the peak for each peak of the pulse wave.
Specifically, calculating the height difference between each peak and the first valley after the peak may be: calculating the difference between the amplitude of each peak and the amplitude of the first trough after the peak; similarly, the height difference between each peak and the second trough after the peak may be calculated as: the difference between the amplitude of each peak and the amplitude of the second trough following said peak is calculated.
And 105, averaging the calculated height differences to obtain a first average value of the height differences, and setting a first threshold according to the first average value of the height differences.
Specifically, the calculated height difference may be averaged as follows: the calculated height difference is subjected to arithmetic averaging or weighted averaging, which is not limited in this embodiment.
Specifically, setting the first threshold according to the first average of the height differences may be: multiplying the first average value of the height difference by a predetermined coefficient to obtain a product as the first threshold value.
The predetermined coefficient may be set according to system performance and/or implementation requirements during specific implementation, and the size of the predetermined coefficient is not limited in this embodiment, for example, the predetermined coefficient may be 1.5.
In this embodiment, after the first threshold is set according to the first mean value, the first threshold may be used to perform first filtering on the pulse wave, compare the amplitude of the extreme point of the pulse wave with the first threshold, and filter the extreme point whose absolute value is greater than or equal to the first threshold, so as to filter an abnormally high value and an abnormally low value in the pulse wave, and effectively remove an abnormal waveform in the pulse wave.
In the method for filtering the pulse wave, after the acquired pulse wave is acquired, the extreme value detection is carried out on the pulse wave to acquire the extreme value point of the pulse wave, then classifying the extreme points of the pulse wave to obtain the wave crest and the wave trough of the pulse wave, calculating a height difference between each peak and a first valley after the peak and a height difference between each peak and a second valley after the peak, respectively, for each peak of the pulse wave, then averaging the height differences obtained by calculation to obtain a first average value of the height differences, setting a first threshold value according to the first average value of the height differences, filtering the extreme points of which the absolute value of the amplitude value is greater than or equal to the first threshold value in the extreme points of the pulse wave, therefore, the abnormal waveform in the pulse wave can be filtered, and the accuracy of positioning the wave crest and the wave trough of the pulse wave is improved.
Fig. 2 is a flowchart illustrating another embodiment of the method for filtering a pulse wave according to the present application, as shown in fig. 2, in the embodiment shown in fig. 1, after step 106, the method may further include:
Specifically, after the first filtering, the remaining extreme points of the pulse wave need to be classified, and the peaks and troughs in the remaining extreme points of the pulse wave need to be distinguished.
The method for calculating the height difference may refer to the description of step 104 in the embodiment shown in fig. 1 of the present application, and is not described herein again.
Also, the calculated height difference may be averaged as: the calculated height difference is subjected to arithmetic averaging or weighted averaging, which is not limited in this embodiment.
The manner of setting the second threshold according to the second average value is the same as the manner of setting the first threshold according to the first average value, and is not described herein again. When the second threshold is set, the predetermined coefficient multiplied by the second average value may be the same as or different from the predetermined coefficient multiplied by the first average value, and is not limited herein.
In this embodiment, after the second threshold is set according to the second average value, the second threshold may be used to perform a second filtering on the pulse wave, compare the amplitude of the remaining extreme point of the pulse wave with the second threshold, and filter out the extreme point whose absolute value of the amplitude is smaller than the second threshold.
Because blood flows outwards from the heart when the heart of a person contracts, if the blood is subjected to resistance in the process of flowing outwards from the heart, the blood flows back towards the heart, the backflow blood collides with the blood flowing outwards from the heart to form a dicrotic wave, and extreme points of the dicrotic wave are generally small, so that the extreme points of the dicrotic wave can be effectively filtered by filtering the extreme points smaller than the second threshold, the initial point and the final point of a single pulse wave can be effectively positioned without considering the problem of baseline drift, and the influence of the dicrotic wave is eliminated.
Fig. 3 is a flowchart illustrating a method for filtering a pulse wave according to another embodiment of the present invention, as shown in fig. 3, in the embodiment of the present invention shown in fig. 1, before step 102, the method further includes:
In this embodiment, after acquiring the acquired pulse wave, noise reduction processing may be performed on the acquired pulse wave through polynomial curve fitting, so as to filter out minute noise points in the acquired pulse wave, and then extremum detection may be performed on the pulse wave to obtain an extremum point of the pulse wave.
In specific implementation, a predetermined sampling frequency (for example, 30Hz) is generally adopted to sample the pulse wave, the waveform of the pulse wave formed by the acquired sampling points is jagged, then a third-order B-spline curve fitting algorithm in a polynomial curve fitting algorithm can be adopted to perform curve fitting on the acquired sampling points, the sampling points which are not on the curve obtained by fitting are filtered out, so that the pulse wave with a smooth waveform is obtained, and then extreme value detection is performed on the pulse wave with a smooth waveform, so that the extreme value point of the pulse wave is obtained.
The polynomial curve fitting algorithm used in this embodiment is a third-order B-spline curve fitting algorithm, but this embodiment is not limited thereto, and other polynomial curve fitting algorithms may also be used, and this embodiment is not limited thereto.
Fig. 4 is a schematic structural diagram of an embodiment of the apparatus for filtering a pulse wave according to the present application, where the apparatus for filtering a pulse wave in the present embodiment can implement the method for filtering a pulse wave according to the present application. As shown in fig. 4, the apparatus for filtering pulse waves may include: an acquisition module 41, a detection module 42, a classification module 43, a calculation module 44 and a filtering module 45;
the acquisition module 41 is configured to acquire the acquired pulse wave;
a detection module 42, configured to perform extremum detection on the pulse wave acquired by the acquisition module 41 to obtain an extremum point of the pulse wave; in a specific implementation, the detecting module 42 may perform extreme value detection on the pulse wave by using an argrelmax algorithm in signal detection, and use a principle of minimum distance when detecting the extreme point of the pulse wave by using the argrelmax algorithm.
Of course, other extreme value detection methods may be used to detect the extreme value point of the pulse wave, which is not limited in this embodiment.
A classification module 43, configured to classify the extreme points of the pulse wave obtained by the detection module 42, so as to obtain peaks and troughs of the pulse wave; specifically, after the detecting module 42 obtains the extreme points of the pulse wave, the classifying module 43 may classify the extreme points of the pulse wave to distinguish peaks and troughs of the extreme points of the pulse wave.
A calculating module 44, configured to calculate, for each peak of the pulse wave, a height difference between each peak and a first trough after the peak, and a height difference between each peak and a second trough after the peak;
specifically, the calculation module 44 may calculate the height difference between each peak and the first valley after the peak as follows: calculating the difference between the amplitude of each peak and the amplitude of the first trough after the peak; similarly, the calculation module 44 may calculate the height difference between each peak and the second trough after the peak as: the difference between the amplitude of each peak and the amplitude of the second trough following said peak is calculated.
The filtering module 45 is configured to average the calculated height differences to obtain a first average value of the height differences, and set a first threshold according to the first average value of the height differences; and filtering the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold value in the extreme points of the pulse waves.
Specifically, the calculated height difference may be averaged as follows: the filtering module 45 performs an arithmetic average on the calculated height difference, or performs a weighted average on the calculated height difference, which is not limited in this embodiment.
In this embodiment, the filtering module 45 is specifically configured to multiply the first average value of the height difference by a predetermined coefficient, so as to obtain a product, which is used as the first threshold. The predetermined coefficient may be set according to system performance and/or implementation requirements during specific implementation, and the size of the predetermined coefficient is not limited in this embodiment, for example, the predetermined coefficient may be 1.5.
In this embodiment, after the first threshold is set according to the first mean value, the filtering module 45 may perform first filtering on the pulse wave by using the first threshold, compare the amplitude of the extreme point of the pulse wave with the first threshold, and filter the extreme point of which the absolute value is greater than or equal to the first threshold, so as to filter the abnormal high value and the abnormal low value in the pulse wave, and effectively remove the abnormal waveform in the pulse wave.
Further, the classification module 43 is further configured to, after the filtering module 45 filters the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold from the extreme points of the pulse wave, classify the remaining extreme points of the pulse wave to obtain peaks and troughs in the remaining extreme points of the pulse wave; specifically, after the first filtering, the classification module 43 needs to classify the remaining extreme points of the pulse wave, and distinguish the peak from the trough in the remaining extreme points of the pulse wave.
The calculating module 44 is further configured to calculate, for each peak in the remaining extreme points, a height difference between each peak and a first trough after the peak, and a height difference between each peak and a second trough after the peak, respectively; averaging the height differences obtained by calculation to obtain a second average value of the height differences, and setting a second threshold value according to the second average value; likewise, the calculation module 44 may average the calculated height difference as follows: the calculated height difference is subjected to arithmetic averaging or weighted averaging, which is not limited in this embodiment.
The manner of setting the second threshold according to the second average value is the same as the manner of setting the first threshold according to the first average value, and is not described herein again. When the second threshold is set, the predetermined coefficient multiplied by the second average value may be the same as or different from the predetermined coefficient multiplied by the first average value, and is not limited herein.
The filtering module 45 is further configured to filter the extreme points where the absolute value of the amplitude in the remaining extreme points of the pulse wave is smaller than the second threshold.
In this embodiment, after the calculating module 44 sets the second threshold according to the second mean value, the second threshold may be used to perform a second filtering on the pulse wave, compare the amplitude of the remaining extreme points of the pulse wave with the second threshold, and filter out the extreme points whose absolute values are smaller than the second threshold.
Because blood flows outwards from the heart when the heart of a person contracts, if the blood is subjected to resistance in the process of flowing outwards from the heart, the blood flows back towards the heart, at the moment, the returned blood collides with the blood flowing outwards from the heart to form a repeating pulse wave, and the extreme point of the repeating pulse wave is generally smaller, so that the filtering module 45 can effectively filter the extreme point of the repeating pulse wave by filtering the extreme point smaller than the second threshold value, the initial point and the final point of a single pulse wave can be effectively positioned without considering the baseline drift problem, and the influence of the repeating pulse wave is eliminated.
Further, the apparatus for filtering a pulse wave may further include: and a noise reduction module 46, configured to perform noise reduction processing on the acquired pulse wave through polynomial curve fitting before the detection module 42 performs extremum detection on the pulse wave to obtain an extremum point of the pulse wave. Specifically, after the obtaining module 41 obtains the collected pulse wave, the noise reduction module 46 may perform noise reduction processing on the collected pulse wave through polynomial curve fitting, filter out tiny noise points in the collected pulse wave, and then perform extremum detection on the pulse wave by the detection module 43 to obtain an extremum point of the pulse wave.
In a specific implementation, a predetermined sampling frequency (e.g., 30Hz) is generally adopted to sample the pulse wave, so that the waveform of the pulse wave formed by the collected sampling points is jagged, then the noise reduction module 46 may adopt a third-order B-spline curve fitting algorithm in a polynomial curve fitting algorithm to perform curve fitting on the collected sampling points, filter out the sampling points that are not on the curve obtained by fitting, thereby obtaining the pulse wave with a smooth waveform, and then perform extreme value detection on the pulse wave with a smooth waveform, so as to obtain the extreme value points of the pulse wave.
The polynomial curve fitting algorithm used in this embodiment is a third-order B-spline curve fitting algorithm, but this embodiment is not limited thereto, and other polynomial curve fitting algorithms may also be used, and this embodiment is not limited thereto.
In the above apparatus for filtering a pulse wave, after the obtaining module 41 obtains the collected pulse wave, the detecting module 42 performs extremum detection on the pulse wave to obtain an extremum point of the pulse wave, then the classifying module 43 classifies the extremum point of the pulse wave to obtain a peak and a trough of the pulse wave, for each peak of the pulse wave, the calculating module 44 calculates a height difference between each peak and a first trough after the peak, and calculates a height difference between each peak and a second trough after the peak, and the filtering module 45 filters the pulse wave according to the calculated height difference, so as to filter an abnormal waveform in the pulse wave and improve accuracy of positioning the peak and the trough of the pulse wave.
Fig. 5 is a schematic structural diagram of an embodiment of a computer device according to the present application, where the computer device may include a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the method for filtering a pulse wave according to the embodiment of the present application may be implemented.
The computer device may be a server, for example: the cloud server may also be an electronic device, for example: the present embodiment does not limit the specific form of the computer device.
FIG. 5 illustrates a block diagram of an exemplary computer device 12 suitable for use in implementing embodiments of the present application. The computer device 12 shown in fig. 5 is only an example and should not bring any limitation to the function and scope of use of the embodiments of the present application.
As shown in FIG. 5, computer device 12 is in the form of a general purpose computing device. The components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that couples various system components including the system memory 28 and the processing unit 16.
The system Memory 28 may include computer system readable media in the form of volatile Memory, such as Random Access Memory (RAM) 30 and/or cache Memory 32. Computer device 12 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, storage system 34 may be used to read from and write to non-removable, nonvolatile magnetic media (not shown in FIG. 5, and commonly referred to as a "hard drive"). Although not shown in FIG. 5, a magnetic disk drive for reading from and writing to a removable nonvolatile magnetic disk (e.g., a "floppy disk") and an optical disk drive for reading from or writing to a removable nonvolatile optical disk (e.g., a Compact disk Read Only Memory (CD-ROM), a Digital versatile disk Read Only Memory (DVD-ROM), or other optical media) may be provided. In these cases, each drive may be connected to bus 18 by one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to carry out the functions of embodiments of the application.
A program/utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28, such program modules 42 including, but not limited to, an operating system, one or more application programs, other program modules, and program data, each of which examples or some combination thereof may comprise an implementation of a network environment. Program modules 42 generally perform the functions and/or methodologies of the embodiments described herein.
The processing unit 16 executes programs stored in the system memory 28 to perform various functional applications and data processing, such as implementing the method for filtering pulse waves provided by the embodiments of the present application.
Embodiments of the present application further provide a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, can implement the method for filtering a pulse wave provided in the embodiments of the present application.
The non-transitory computer readable storage medium described above may take any combination of one or more computer readable media. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a Read Only Memory (ROM), an Erasable Programmable Read Only Memory (EPROM), a flash Memory, an optical fiber, a portable compact disc Read Only Memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
Computer program code for carrying out operations for aspects of the present application may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + +, and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of Network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
In the description herein, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present application, "plurality" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present application.
The word "if" as used herein may be interpreted as "at … …" or "when … …" or "in response to a determination" or "in response to a detection", depending on the context. Similarly, the phrases "if determined" or "if detected (a stated condition or event)" may be interpreted as "when determined" or "in response to a determination" or "when detected (a stated condition or event)" or "in response to a detection (a stated condition or event)", depending on the context.
It should be noted that the terminal according to the embodiments of the present application may include, but is not limited to, a Personal Computer (Personal Computer; hereinafter, referred to as PC), a Personal Digital Assistant (Personal Digital Assistant; hereinafter, referred to as PDA), a wireless handheld device, a Tablet Computer (Tablet Computer), a mobile phone, an MP3 player, an MP4 player, and the like.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions in actual implementation, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, or in a form of hardware plus a software functional unit.
The integrated unit implemented in the form of a software functional unit may be stored in a computer readable storage medium. The software functional unit is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) or a Processor (Processor) to execute some steps of the methods according to the embodiments of the present application. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above description is only a preferred embodiment of the present application and should not be taken as limiting the present application, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (8)
1. A method of filtering a pulse wave, comprising:
acquiring the acquired pulse wave;
carrying out extreme value detection on the pulse wave to obtain an extreme value point of the pulse wave;
classifying the extreme points of the pulse waves to obtain peaks and troughs of the pulse waves;
for each peak of the pulse wave, respectively calculating the height difference between each peak and a first trough after the peak, and calculating the height difference between each peak and a second trough after the peak;
averaging the height differences obtained by calculation to obtain a first average value of the height differences, and setting a first threshold value according to the first average value of the height differences;
filtering the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold value in the extreme points of the pulse wave;
classifying the residual extreme points of the pulse waves to obtain wave crests and wave troughs in the residual extreme points of the pulse waves;
for each peak in the remaining extreme points, respectively calculating the height difference between each peak and a first trough after the peak, and calculating the height difference between each peak and a second trough after the peak;
averaging the height difference obtained by calculation to obtain a second average value of the height difference, and setting a second threshold value according to the second average value;
and filtering the extreme points of which the absolute value of the amplitude values in the residual extreme points of the pulse wave is smaller than the second threshold value.
2. The method of claim 1, wherein setting the first threshold value according to the first mean of the height differences comprises:
multiplying the first mean value of the height difference by a predetermined coefficient to obtain a product as the first threshold value.
3. The method according to any one of claims 1-2, wherein before the extreme value detection of the pulse wave to obtain the extreme value point of the pulse wave, the method further comprises:
and carrying out noise reduction processing on the acquired pulse wave through polynomial curve fitting.
4. An apparatus for filtering a pulse wave, comprising:
the acquisition module is used for acquiring the acquired pulse waves;
the detection module is used for carrying out extreme value detection on the pulse waves acquired by the acquisition module to acquire extreme value points of the pulse waves;
the classification module is used for classifying the extreme points of the pulse waves obtained by the detection module to obtain the wave crests and the wave troughs of the pulse waves;
the calculation module is used for calculating the height difference between each peak and a first trough behind the peak and the height difference between each peak and a second trough behind the peak respectively for each peak of the pulse waves;
the filtering module is used for averaging the height difference obtained by calculation to obtain a first average value of the height difference and setting a first threshold value according to the first average value of the height difference; filtering the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold value in the extreme points of the pulse wave;
the classification module is further configured to classify the remaining extreme points of the pulse wave after the filtering module filters the extreme points of which the absolute value of the amplitude is greater than or equal to the first threshold from the extreme points of the pulse wave, so as to obtain peaks and troughs of the remaining extreme points of the pulse wave;
the calculation module is further configured to calculate, for each peak in the remaining extreme points, a height difference between each peak and a first trough after the peak, and a height difference between each peak and a second trough after the peak, respectively; averaging the height difference obtained by calculation to obtain a second average value of the height difference, and setting a second threshold value according to the second average value;
the filtering module is further configured to filter the extreme point, where an absolute value of the amplitude in the remaining extreme points of the pulse wave is smaller than the second threshold.
5. The apparatus of claim 4,
the filtering module is specifically configured to multiply the first mean value of the height difference by a predetermined coefficient, and use an obtained product as the first threshold.
6. The apparatus of any one of claims 4-5, further comprising:
and the noise reduction module is used for carrying out noise reduction processing on the collected pulse waves through polynomial curve fitting before the detection module carries out extreme value detection on the pulse waves and obtains extreme values of the pulse waves.
7. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method according to any one of claims 1-3 when executing the computer program.
8. A non-transitory computer-readable storage medium having stored thereon a computer program, wherein the computer program, when executed by a processor, implements the method of any one of claims 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910012502.9A CN109864705B (en) | 2019-01-07 | 2019-01-07 | Method and device for filtering pulse wave and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910012502.9A CN109864705B (en) | 2019-01-07 | 2019-01-07 | Method and device for filtering pulse wave and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109864705A CN109864705A (en) | 2019-06-11 |
| CN109864705B true CN109864705B (en) | 2022-08-12 |
Family
ID=66917495
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910012502.9A Active CN109864705B (en) | 2019-01-07 | 2019-01-07 | Method and device for filtering pulse wave and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109864705B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112155531A (en) * | 2019-11-29 | 2021-01-01 | 苏州润迈德医疗科技有限公司 | Method and system for obtaining diastolic pressure and systolic pressure extreme points according to invasive blood pressure |
| CN110960203B (en) * | 2019-12-13 | 2022-11-25 | 心核心科技(北京)有限公司 | Cardiovascular characteristic parameter detection method and device |
| CN111551952B (en) * | 2020-05-08 | 2022-10-18 | 中国地质大学（武汉） | Extraction method of water depth measurement laser data, water depth measurement method and device |
| CN112932428A (en) * | 2021-02-05 | 2021-06-11 | 江西国科美信医疗科技有限公司 | Blood oxygen and heart rate monitoring system and portable wireless remote monitoring device thereof |
| CN113297943B (en) * | 2021-05-18 | 2022-12-02 | 北京远舢智能科技有限公司 | Equipment auxiliary control technology based on mixed reality |
| CN113545753B (en) * | 2021-07-21 | 2024-04-05 | 芯原微电子（成都）有限公司 | Waveform detection method and device, wearable equipment and computer readable storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1813628A (en) * | 2006-01-18 | 2006-08-09 | 哈尔滨工业大学 | Pulse condition detecting device capable of automatically regulating pulse-feeling pressure and pulse condition characteristic extracting method |
| CN102626307A (en) * | 2012-01-16 | 2012-08-08 | 兰州理工大学 | Dynamic pulse signal real-time detection system and detection method |
| CN103027667A (en) * | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | Characteristic parameter extraction of pulse wave |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016168979A1 (en) * | 2015-04-20 | 2016-10-27 | 深圳市长桑技术有限公司 | Vital sign analysis method and system |
-
2019
- 2019-01-07 CN CN201910012502.9A patent/CN109864705B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1813628A (en) * | 2006-01-18 | 2006-08-09 | 哈尔滨工业大学 | Pulse condition detecting device capable of automatically regulating pulse-feeling pressure and pulse condition characteristic extracting method |
| CN103027667A (en) * | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | Characteristic parameter extraction of pulse wave |
| CN102626307A (en) * | 2012-01-16 | 2012-08-08 | 兰州理工大学 | Dynamic pulse signal real-time detection system and detection method |
Non-Patent Citations (1)
| Title |
|---|
| 动态脉搏信号检测与脉率变异性实时分析方法研究;丑永新;《中国博士学位论文全文数据库（电子期刊） 信息科技辑》;20161215(第12期);第I136-2页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109864705A (en) | 2019-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109864705B (en) | Method and device for filtering pulse wave and computer equipment | |
| JP6811773B2 (en) | How to Quantify PhotoPretismogram (PPG) Signal Quality | |
| CN110461215B (en) | Determining health signs using a portable device | |
| CN109452935B (en) | Non-invasive method and system for estimating blood pressure from a vascular plethysmogram using statistical post-processing | |
| US9380956B2 (en) | Method and apparatus for classifying cardiac arrhythmia | |
| CN108742697B (en) | Heart sound signal classification method and terminal equipment | |
| Paula Rocha et al. | Beyond long memory in heart rate variability: an approach based on fractionally integrated autoregressive moving average time series models with conditional heteroscedasticity | |
| KR101977745B1 (en) | Method and apparatus for classifying cardiac arrhythmia | |
| CN107273827A (en) | A kind of R wave of electrocardiosignal detection method and device | |
| CN112006673A (en) | Human body heart rate detection method and system, storage medium, computer equipment and terminal | |
| CN106821336A (en) | Sleep monitor method and system | |
| Guven et al. | Computationally efficient real‐time interpolation algorithm for non‐uniform sampled biosignals | |
| CN111079560B (en) | Tumble monitoring method and device and terminal equipment | |
| CN116327133A (en) | Multi-physiological index detection method, device and related equipment | |
| US20220211283A1 (en) | Methods for blood pressure calibration selection and modeling methods thereof | |
| CN116269297A (en) | Vital sign detection method and device based on radar, radar and storage medium | |
| CN113616217B (en) | Method and device for generating baseline drift curve | |
| JP6706996B2 (en) | Biological signal processing device, abnormality determination method and program | |
| CN114724693A (en) | Method and device for detecting abnormal diagnosis and treatment behaviors, electronic equipment and storage medium | |
| CN112598033B (en) | Physiological signal processing method, device, equipment and storage medium | |
| CN116776087A (en) | Heart rate detection method and related equipment | |
| CN113040757B (en) | Head posture monitoring method and device, head intelligent wearable device and storage medium | |
| CN110664392B (en) | Blood flow velocity calculation method and device, terminal equipment and storage medium | |
| CN113743704A (en) | Event monitoring method and device, computer equipment and storage medium | |
| US20190163837A1 (en) | Digital data filtering method, apparatus, and terminal device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |