CN109864705A - The method, apparatus and computer equipment that pulse wave is filtered - Google Patents
The method, apparatus and computer equipment that pulse wave is filtered Download PDFInfo
- Publication number
- CN109864705A CN109864705A CN201910012502.9A CN201910012502A CN109864705A CN 109864705 A CN109864705 A CN 109864705A CN 201910012502 A CN201910012502 A CN 201910012502A CN 109864705 A CN109864705 A CN 109864705A
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- wave
- difference
- height
- extreme point
- 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.)
- Granted
Links
Abstract
Present applicant proposes method, apparatus and computer equipment that a kind of pair of pulse wave is filtered, wherein above-mentioned includes: to obtain collected pulse wave to the method that pulse wave is filtered；Extremum extracting is carried out to the pulse wave, obtains the extreme point of the pulse wave；Classify to the extreme point of the pulse wave, obtains the wave crest and trough of the pulse wave；To each wave crest of the pulse wave, the difference in height of first trough after each wave crest and the wave crest is calculated separately, and calculates the difference in height of second trough after each wave crest and the wave crest；It is averaged to the difference in height obtained is calculated, obtains the first mean value of difference in height, and first threshold is arranged according to the first mean value of the difference in height；The extreme point that the absolute value of amplitude in the extreme point of the pulse wave is greater than or equal to the first threshold is filtered out.The application may be implemented to filter out the unusual waveforms in pulse wave, improve to the wave crest of pulse wave and the accuracy of Trough Location.
Description
[technical field]
The method, apparatus that is filtered this application involves signal processing technology field more particularly to a kind of pair of pulse wave and
Computer equipment.
[background technique]
Pulse, breathing etc. are the important vital sign of human body, intensity, form and the frequency that these vital signs are presented
It is able to reflect out physical condition, the state of mind and the healthy water equality physiological and pathological information of human body.
Foundation of the physiological and pathological information of human body as clinical diagnosis and treatment is extracted from pulse wave, always by
The attention of outer medical field the reasons such as unskilled or physical fitness would generally produce during acquiring pulse wave because of user's operation
First portion unusual waveforms, such as: dicrotic wave and main wave are in same level or because abnormal high level or different occurs in external factor
Normal low value.
It is existing in the related technology, the general method using difference or sliding window mentions the Wave crest and wave trough of pulse wave
It takes, but above two method can not identify the above-mentioned unusual waveforms in pulse wave, it is even more impossible to remove above-mentioned unusual waveforms,
It is inaccurate so as to cause the positioning of wave crest and trough to pulse wave.
[summary of the invention]
The embodiment of the present application provides the method, apparatus and computer equipment that a kind of pair of pulse wave is filtered, to realize
Unusual waveforms in pulse wave are filtered out, are improved to the wave crest of pulse wave and the accuracy of Trough Location.
In a first aspect, the embodiment of the present application provides the method that a kind of pair of pulse wave is filtered, comprising: acquisition collects
Pulse wave；Extremum extracting is carried out to the pulse wave, obtains the extreme point of the pulse wave；To the extreme point of the pulse wave
Classify, obtains the wave crest and trough of the pulse wave；To each wave crest of the pulse wave, calculate separately each wave crest with
The difference in height of first trough after the wave crest, and calculate each wave crest and second trough after the wave crest
Difference in height；It is averaged to the difference in height obtained is calculated, obtains the first mean value of difference in height, and according to the first of the difference in height
First threshold is arranged in mean value；The absolute value of amplitude in the extreme point of the pulse wave is greater than or equal to the pole of the first threshold
Value point filters out.
Wherein in one possible implementation, the absolute value by amplitude in the extreme point of the pulse wave is greater than
Or after being filtered out equal to the extreme point of the first threshold, further includes: classify to the remaining extreme point of the pulse wave, obtain
Obtain the wave crest and trough in the pulse wave residue extreme point；For each wave crest in the remaining extreme point, calculate separately
The difference in height of first trough after each wave crest and the wave crest, and calculate after each wave crest and the wave crest the
The difference in height of two troughs；It is averaged to the difference in height obtained is calculated, obtains the second mean value of difference in height, and according to described the
Second threshold is arranged in two mean values；The absolute value of amplitude in the remaining extreme point of the pulse wave is less than to the pole of the second threshold
Value point filters out.
Wherein in one possible implementation, described that first threshold packet is arranged according to the first mean value of the difference in height
It includes: the first mean value of the difference in height is multiplied with pre-determined factor, using the product of acquisition as the first threshold.
Wherein in one possible implementation, described that extremum extracting is carried out to the pulse wave, obtain the pulse
Before the extreme point of wave, further includes: carry out noise reduction process to collected pulse wave by polynomial curve fitting.
Second aspect, the embodiment of the present application provide the device that a kind of pair of pulse wave is filtered, comprising: obtain module, use
In the collected pulse wave of acquisition；Detection module, the pulse wave for obtaining to the acquisition module carry out extremum extracting, obtain
The extreme point of the pulse wave；The extreme point of categorization module, the pulse wave for obtaining to the detection module is classified, and is obtained
Obtain the wave crest and trough of the pulse wave；Computing module calculates separately each wave crest for each wave crest to the pulse wave
With second trough after the difference in height of first trough after the wave crest, and each wave crest of calculating and the wave crest
Difference in height；Filter module obtains the first mean value of difference in height for being averaged to calculating the difference in height obtained, and according to
First threshold is arranged in first mean value of the difference in height；The absolute value of amplitude in the extreme point of the pulse wave is greater than or equal to
The extreme point of the first threshold filters out.
Wherein in one possible implementation, the categorization module is also used to the arteries and veins in the filter module
After the extreme point that the absolute value of amplitude is greater than or equal to the first threshold in the extreme point of wave of fighting filters out, to the pulse wave
Remaining extreme point classify, obtain the wave crest and trough in the pulse wave residue extreme point；The computing module, is also used
In for each wave crest in the remaining extreme point, first trough after each wave crest and the wave crest is calculated separately
Difference in height, and calculate the difference in height of second trough after each wave crest and the wave crest；To the difference in height for calculating acquisition
It is averaged, obtains the second mean value of difference in height, and second threshold is arranged according to second mean value；The filter module, also
The extreme point that absolute value for amplitude in the remaining extreme point by the pulse wave is less than the second threshold filters out.
Wherein in one possible implementation, the filter module, specifically for equal by the first of the difference in height
Value is multiplied with pre-determined factor, using the product of acquisition as the first threshold.
Wherein in one possible implementation, the device that pulse wave is filtered further include: noise reduction module,
For carrying out extremum extracting to the pulse wave in the detection module, before the extreme point for obtaining the pulse wave, by more
Item formula curve matching carries out noise reduction process to collected pulse wave.
The third aspect, the embodiment of the present application provide a kind of computer equipment, including memory, processor and are stored in described
It is real when the processor executes the computer program on memory and the computer program that can run on the processor
Now method as described above.
Fourth aspect, the embodiment of the present application provide a kind of non-transitorycomputer readable storage medium, are stored thereon with meter
Calculation machine program, the computer program realize method as described above when being executed by processor.
In above technical scheme, after obtaining collected pulse wave, extremum extracting is carried out to above-mentioned pulse wave, in acquisition
The extreme point for stating pulse wave classifies to the extreme point of above-mentioned pulse wave, the wave crest and trough of above-mentioned pulse wave is obtained, to upper
The each wave crest for stating pulse wave calculates separately the difference in height of first trough after each wave crest and above-mentioned wave crest, Yi Jiji
The difference in height of second trough after each wave crest and above-mentioned wave crest is calculated, is then averaged to the difference in height for calculating acquisition,
The first mean value of difference in height is obtained, and first threshold is arranged according to the first mean value of above-mentioned difference in height, by the pole of above-mentioned pulse wave
The extreme point that the absolute value of amplitude is greater than or equal to above-mentioned first threshold in value point filters out, so as to realize in pulse wave
Unusual waveforms are filtered out, and are improved to the wave crest of pulse wave and the accuracy of Trough Location.
[Detailed description of the invention]
Technical solution in ord to more clearly illustrate embodiments of the present application, below will be to needed in the embodiment attached
Figure is briefly described, it should be apparent that, the drawings in the following description are only some examples of the present application, for this field
For those of ordinary skill, without creative efforts, it can also be obtained according to these attached drawings other attached drawings.
Fig. 1 is the flow chart of method one embodiment that the application is filtered pulse wave；
Fig. 2 is the flow chart for method another embodiment that the application is filtered pulse wave；
Fig. 3 is the flow chart for the method further embodiment that the application is filtered pulse wave；
Fig. 4 is the structural schematic diagram of device one embodiment that the application is filtered pulse wave；
Fig. 5 is the structural schematic diagram of the application computer equipment one embodiment.
[specific embodiment]
In order to better understand the technical solution of the application, the embodiment of the present application is retouched in detail with reference to the accompanying drawing
It states.
It will be appreciated that described embodiments are only a part of embodiments of the present application, instead of all the embodiments.Base
Embodiment in the application, it is obtained by those of ordinary skill in the art without making creative efforts it is all its
Its embodiment, shall fall in the protection scope of this application.
The term used in the embodiment of the present application is only to be not intended to be limiting merely for for the purpose of describing particular embodiments
The application.In the embodiment of the present application and the "an" of singular used in the attached claims, " described " and "the"
It is also intended to including most forms, unless the context clearly indicates other meaning.
Fig. 1 is the flow chart of method one embodiment that the application is filtered pulse wave, as shown in Figure 1, above-mentioned right
The method that pulse wave is filtered may include:
Step 101, collected pulse wave is obtained.
Step 102, extremum extracting is carried out to above-mentioned pulse wave, obtains the extreme point of above-mentioned pulse wave.
In specific implementation, extreme value inspection can be carried out to above-mentioned pulse wave using the argrelmax algorithm in signal detection
It surveys, and when detecting the extreme point of above-mentioned pulse wave using argrelmax algorithm, using apart from the smallest principle.
It is of course also possible to use other extremum extracting methods detect the extreme point of above-mentioned pulse wave, the present embodiment is to this
It is not construed as limiting.
Step 103, classify to the extreme point of above-mentioned pulse wave, obtain the wave crest and trough of above-mentioned pulse wave.
Specifically, after the extreme point for obtaining above-mentioned pulse wave, can classify to the extreme point of above-mentioned pulse wave, it will
Wave crest and trough in the extreme point of above-mentioned pulse wave distinguish.
Step 104, to each wave crest of above-mentioned pulse wave, each wave crest and first after above-mentioned wave crest are calculated separately
The difference in height of trough, and calculate the difference in height of second trough after each wave crest and above-mentioned wave crest.
Specifically, the difference in height for calculating first trough after each wave crest and above-mentioned wave crest can be with are as follows: calculates each
The difference of the amplitude of first trough after the amplitude of wave crest and above-mentioned wave crest；Equally, each wave crest and above-mentioned wave crest are calculated
The difference in height of second trough later can be with are as follows: the amplitude for calculating each wave crest and second trough after above-mentioned wave crest
The difference of amplitude.
Step 105, it is averaged to the difference in height for calculating acquisition, obtains the first mean value of difference in height, and according to above-mentioned height
First threshold is arranged in first mean value of degree difference.
Specifically, to calculate that the difference in height obtained carries out averagely can be with are as follows: to the difference in height obtained is calculated, to carry out arithmetic flat
, or to the difference in height obtained is calculated it is weighted and averaged, the present embodiment is not construed as limiting this.
Specifically, first threshold is arranged according to the first mean value of above-mentioned difference in height can be with are as follows: by the first of above-mentioned difference in height
Mean value is multiplied with pre-determined factor, using the product of acquisition as above-mentioned first threshold.
Wherein, above-mentioned pre-determined factor can voluntarily be set according to system performance and/or realization demand etc. in specific implementation
Fixed, the present embodiment is not construed as limiting the size of above-mentioned pre-determined factor, for example, above-mentioned pre-determined factor can be 1.5.
Step 106, the absolute value of amplitude in the extreme point of above-mentioned pulse wave is greater than or equal to the pole of above-mentioned first threshold
Value point filters out.
In the present embodiment, after first threshold is set according to the first mean value, first threshold can use to above-mentioned pulse
Wave carries out first time filtering, the amplitude of the extreme point of above-mentioned pulse wave is compared with first threshold, by the absolute value of amplitude
Extreme point more than or equal to above-mentioned first threshold filters out, to filter out the abnormal high level in pulse wave and abnormal low value, effectively
Remove the unusual waveforms in pulse wave.
In the above-mentioned method being filtered to pulse wave, after obtaining collected pulse wave, above-mentioned pulse wave is carried out
Extremum extracting obtains the extreme point of above-mentioned pulse wave, then classifies to the extreme point of above-mentioned pulse wave, obtains above-mentioned pulse
The wave crest and trough of wave calculate separately each wave crest and first after above-mentioned wave crest to each wave crest of above-mentioned pulse wave
The difference in height of trough, and the difference in height of second trough after each wave crest and above-mentioned wave crest is calculated, then calculating is obtained
The difference in height obtained is averaged, and obtains the first mean value of difference in height, and the first threshold is arranged according to the first mean value of above-mentioned difference in height
Value filters out the extreme point that the absolute value of amplitude in the extreme point of above-mentioned pulse wave is greater than or equal to above-mentioned first threshold, thus
It may be implemented to filter out the unusual waveforms in pulse wave, improve to the wave crest of pulse wave and the accuracy of Trough Location.
Fig. 2 is the flow chart for method another embodiment that the application is filtered pulse wave, as shown in Fig. 2, this Shen
Please after step 106, can also include: in embodiment illustrated in fig. 1
Step 201, classify to the remaining extreme point of above-mentioned pulse wave, obtain in the remaining extreme point of above-mentioned pulse wave
Wave crest and trough.
Specifically, it after carrying out first time filtering, needs again to classify to the remaining extreme point of above-mentioned pulse wave, it will
Wave crest and trough in the remaining extreme point of above-mentioned pulse wave distinguish.
Step 202, for each wave crest in above-mentioned remaining extreme point, after calculating separately each wave crest and above-mentioned wave crest
First trough difference in height, and calculate the difference in height of second trough after each wave crest and above-mentioned wave crest.
Wherein, the mode of computed altitude difference may refer to the description of the application embodiment illustrated in fig. 1 step 104, herein not
It repeats again.
Step 203, it is averaged to calculating the difference in height obtained, obtains the second mean value of difference in height, and according to above-mentioned the
Second threshold is arranged in two mean values.
Equally, to calculate that the difference in height obtained carries out averagely can be with are as follows: carry out arithmetic average to the difference in height obtained is calculated,
Or be weighted and averaged to the difference in height obtained is calculated, the present embodiment is not construed as limiting this.
Wherein, the mode of second threshold and the side that first threshold is arranged according to the first mean value are set according to above-mentioned second mean value
Formula is identical, and details are not described herein.When second threshold is arranged, above-mentioned second mean value multiplied by pre-determined factor can be with the first mean value
The pre-determined factor multiplied is identical or different, is not limited thereto.
Step 204, the absolute value of amplitude in the remaining extreme point of above-mentioned pulse wave is less than to the extreme value of above-mentioned second threshold
Point filters out.
In the present embodiment, after second threshold is set according to the second mean value, second threshold can use to above-mentioned pulse
Wave carries out second and filters, and the amplitude of the remaining extreme point of above-mentioned pulse wave is compared with second threshold, by the exhausted of amplitude
The extreme point for being less than above-mentioned second threshold to value filters out.
Since when human heart is shunk, blood is from heart to outflow, if blood is in the process flowed outward from heart
In by resistance, can flow back toward the direction of heart, the blood at this moment to flow back bumps against with from heart to the blood of outflow, will form weight
Rich wave, the extreme point for winning wave again is generally all smaller, therefore by filtering out the extreme point less than second threshold, so that it may effectively filter
Except wave extreme point is won again, without considering baseline drift problem, so that it may the starting of the single waveform of effective position pulse wave
Point and terminal, eliminate the influence of dicrotic wave.
Fig. 3 is the flow chart for the method further embodiment that the application is filtered pulse wave, as shown in figure 3, this Shen
It please be in embodiment illustrated in fig. 1, before step 102, further includes:
Step 301, noise reduction process is carried out to collected pulse wave by polynomial curve fitting.
It, can be by polynomial curve fitting to collected after obtaining collected pulse wave in the present embodiment
Pulse wave carries out noise reduction process, filters out the small noise spot in collected pulse wave, then carries out pole to above-mentioned pulse wave again
Value detection, obtains the extreme point of above-mentioned pulse wave.
In specific implementation, generally using scheduled sample frequency (such as: 30Hz) pulse wave is sampled, is adopted in this way
The waveform for the pulse wave that the sampled point collected is constituted be it is jagged, then can be using in polynomial curve fitting algorithm
Three rank B-spline curves fitting algorithms, carry out curve fitting to collected sampled point, filter out not on the curve that fitting obtains
Sampled point, to obtain the smooth pulse wave of waveform, then the pulse wave smooth to waveform carries out extremum extracting again, obtains above-mentioned
The extreme point of pulse wave.
Polynomial curve fitting algorithm used in the present embodiment is three rank B-spline curves fitting algorithms, but this implementation
Example is not limited to that, can also use other polynomial curve fitting algorithms, and the present embodiment is not construed as limiting this.
Fig. 4 is the structural schematic diagram of device one embodiment that the application is filtered pulse wave, in the present embodiment
The method provided by the embodiments of the present application being filtered to pulse wave may be implemented to the device that pulse wave is filtered.Such as Fig. 4
Shown, above-mentioned may include: to obtain module 41, detection module 42, categorization module 43, calculating to the device that pulse wave is filtered
Module 44 and filter module 45；
Wherein, module 41 is obtained, for obtaining collected pulse wave；
Detection module 42 obtains above-mentioned pulse wave for carrying out extremum extracting to the pulse wave for obtaining the acquisition of module 41
Extreme point；In specific implementation, detection module 42 can using the argrelmax algorithm in signal detection to above-mentioned pulse wave into
Row extremum extracting, and when detecting the extreme point of above-mentioned pulse wave using argrelmax algorithm, using apart from the smallest original
Then.
It is of course also possible to use other extremum extracting methods detect the extreme point of above-mentioned pulse wave, the present embodiment is to this
It is not construed as limiting.
The extreme point of categorization module 43, the pulse wave for obtaining to detection module 42 is classified, and above-mentioned pulse is obtained
The wave crest and trough of wave；Specifically, after detection module 42 obtains the extreme point of above-mentioned pulse wave, categorization module 43 can be to upper
The extreme point for stating pulse wave is classified, by the extreme point of above-mentioned pulse wave wave crest and trough distinguish.
Computing module 44, for each wave crest to above-mentioned pulse wave, after calculating separately each wave crest and above-mentioned wave crest
First trough difference in height, and calculate the difference in height of second trough after each wave crest and above-mentioned wave crest；
Specifically, the difference in height that computing module 44 calculates first trough after each wave crest and above-mentioned wave crest can be with
Are as follows: calculate the difference of the amplitude of each wave crest and the amplitude of first trough after above-mentioned wave crest；Equally, computing module 44 is counted
The difference in height for calculating second trough after each wave crest and above-mentioned wave crest can be with are as follows: calculate each wave crest amplitude and above-mentioned wave
The difference of the amplitude of second trough after peak.
Filter module 45 obtains the first mean value of difference in height for being averaged to calculating the difference in height obtained, and according to
First threshold is arranged in first mean value of above-mentioned difference in height；The absolute value of amplitude in the extreme point of above-mentioned pulse wave is greater than or equal to
The extreme point of the first threshold filters out.
Specifically, to calculate that the difference in height obtained carries out averagely can be with are as follows: 45 pairs of filter module calculate the difference in height obtained
Arithmetic average is carried out, or is weighted and averaged to the difference in height obtained is calculated, the present embodiment is not construed as limiting this.
In the present embodiment, filter module 45, specifically for the first mean value of above-mentioned difference in height is multiplied with pre-determined factor, with
The product of acquisition is as above-mentioned first threshold.Wherein, above-mentioned pre-determined factor can in specific implementation, according to system performance and/
Or the sets itselfs such as realization demand, the present embodiment are not construed as limiting the size of above-mentioned pre-determined factor, for example, above-mentioned predetermined system
Number can be 1.5.
In the present embodiment, after first threshold is arranged according to the first mean value, filter module 45 can use first threshold
First time filtering is carried out to above-mentioned pulse wave, the amplitude of the extreme point of above-mentioned pulse wave is compared with first threshold, by width
The extreme point that the absolute value of value is greater than or equal to above-mentioned first threshold filters out, to filter out abnormal high level and the exception in pulse wave
Low value effectively removes the unusual waveforms in pulse wave.
Further, categorization module 43, be also used to filter module 45 by the extreme point of above-mentioned pulse wave amplitude it is exhausted
After being filtered out to the extreme point that value is greater than or equal to above-mentioned first threshold, classify to the remaining extreme point of above-mentioned pulse wave,
Obtain the wave crest and trough in above-mentioned pulse wave residue extreme point；Specifically, after carrying out first time filtering, categorization module 43
It needs again to classify to the remaining extreme point of above-mentioned pulse wave, by the wave crest and trough in the remaining extreme point of above-mentioned pulse wave
It distinguishes.
Computing module 44, is also used to for each wave crest in above-mentioned remaining extreme point, calculate separately each wave crest with it is upper
The difference in height of first trough after wave crest is stated, and calculates the height of second trough after each wave crest and above-mentioned wave crest
It is poor to spend；Be averaged to the difference in height obtained is calculated, obtain the second mean value of difference in height, and is arranged the according to above-mentioned second mean value
Two threshold values；Equally, 44 pairs of computing module calculate that the differences in height obtained carry out averagely can be with are as follows: carried out to the difference in height obtained is calculated
Arithmetic average, or be weighted and averaged to the difference in height obtained is calculated, the present embodiment is not construed as limiting this.
Wherein, the mode of second threshold and the side that first threshold is arranged according to the first mean value are set according to above-mentioned second mean value
Formula is identical, and details are not described herein.When second threshold is arranged, above-mentioned second mean value multiplied by pre-determined factor can be with the first mean value
The pre-determined factor multiplied is identical or different, is not limited thereto.
Filter module 45, the absolute value for being also used to amplitude in the remaining extreme point by above-mentioned pulse wave are less than above-mentioned second threshold
The extreme point of value filters out.
In the present embodiment, after second threshold is arranged according to the second mean value in computing module 44, second threshold can use
Second is carried out to above-mentioned pulse wave to filter, and the amplitude of the remaining extreme point of above-mentioned pulse wave is compared with second threshold,
The extreme point that the absolute value of amplitude is less than above-mentioned second threshold is filtered out.
Since when human heart is shunk, blood is from heart to outflow, if blood is in the process flowed outward from heart
In by resistance, can flow back toward the direction of heart, the blood at this moment to flow back bumps against with from heart to the blood of outflow, will form weight
Rich wave, the extreme point for winning wave again is generally all smaller, therefore filter module 45 passes through the extreme point filtered out less than second threshold, just
It can effectively filter out and win wave extreme point again, without considering baseline drift problem, so that it may which effective position pulse wave is single
The starting point and terminal of waveform, eliminate the influence of dicrotic wave.
Further, the above-mentioned device being filtered to pulse wave can also include: noise reduction module 46, in detection mould
Block 42 before the extreme point for obtaining above-mentioned pulse wave, passes through polynomial curve fitting pair to above-mentioned pulse wave progress extremum extracting
Collected pulse wave carries out noise reduction process.Specifically, after obtaining module 41 and obtaining collected pulse wave, noise reduction module
46 can carry out noise reduction process to collected pulse wave by polynomial curve fitting, filter out micro- in collected pulse wave
Then small noise spot carries out extremum extracting to above-mentioned pulse wave by detection module 43 again, obtains the extreme point of above-mentioned pulse wave.
In specific implementation, generally using scheduled sample frequency (such as: 30Hz) pulse wave is sampled, is adopted in this way
The waveform for the pulse wave that the sampled point collected is constituted be it is jagged, then noise reduction module 46 can be quasi- using polynomial curve
Three rank B-spline curves fitting algorithms in hop algorithm, carry out curve fitting to collected sampled point, filter out and do not obtain in fitting
Curve on sampled point, to obtain the smooth pulse wave of waveform, then the pulse wave smooth to waveform carries out extreme value inspection again
It surveys, obtains the extreme point of above-mentioned pulse wave.
Polynomial curve fitting algorithm used in the present embodiment is three rank B-spline curves fitting algorithms, but this implementation
Example is not limited to that, can also use other polynomial curve fitting algorithms, and the present embodiment is not construed as limiting this.
In the above-mentioned device being filtered to pulse wave, after acquisition module 41 obtains collected pulse wave, mould is detected
Block 42 carries out extremum extracting to above-mentioned pulse wave, obtains the extreme point of above-mentioned pulse wave, then categorization module 43 is to above-mentioned pulse
The extreme point of wave is classified, and the wave crest and trough of above-mentioned pulse wave are obtained, to each wave crest of above-mentioned pulse wave, computing module
44 calculate separately the difference in height of first trough after each wave crest and above-mentioned wave crest, and calculate each wave crest and above-mentioned wave
The difference in height of second trough after peak, filter module 45 filter above-mentioned pulse wave according to the difference in height obtained is calculated
Wave filters out the unusual waveforms in pulse wave so as to realize, improves to the wave crest of pulse wave and the standard of Trough Location
True property.
Fig. 5 is the structural schematic diagram of the application computer equipment one embodiment, and above-mentioned computer equipment may include depositing
Reservoir, processor and it is stored in the computer program that can be run on above-mentioned memory and on above-mentioned processor, above-mentioned processor
When executing above-mentioned computer program, the method provided by the embodiments of the present application being filtered to pulse wave may be implemented.
Wherein, above-mentioned computer equipment can be server, such as: Cloud Server, or electronic equipment, such as:
The intelligent electronic devices such as smart phone, smartwatch or tablet computer, specific form of the present embodiment to above-mentioned computer equipment
It is not construed as limiting.
Fig. 5 shows the block diagram for being suitable for the exemplary computer device 12 for being used to realize the application embodiment.Fig. 5 is shown
Computer equipment 12 be only an example, should not function to the embodiment of the present application and use scope bring any restrictions.
As shown in figure 5, computer equipment 12 is showed in the form of universal computing device.The component of computer equipment 12 can be with
Including but not limited to: one or more processor or processing unit 16, system storage 28 connect different system components
The bus 18 of (including system storage 28 and processing unit 16).
Bus 18 indicates one of a few class bus structures or a variety of, including memory bus or Memory Controller,
Peripheral bus, graphics acceleration port, processor or the local bus using any bus structures in a variety of bus structures.It lifts
For example, these architectures include but is not limited to industry standard architecture (Industry Standard
Architecture；Hereinafter referred to as: ISA) bus, microchannel architecture (Micro Channel Architecture；Below
Referred to as: MAC) bus, enhanced isa bus, Video Electronics Standards Association (Video Electronics Standards
Association；Hereinafter referred to as: VESA) local bus and peripheral component interconnection (Peripheral Component
Interconnection；Hereinafter referred to as: PCI) bus.
Computer equipment 12 typically comprises a variety of computer system readable media.These media can be it is any can be by
The usable medium that computer equipment 12 accesses, including volatile and non-volatile media, moveable and immovable medium.
System storage 28 may include the computer system readable media of form of volatile memory, such as arbitrary access
Memory (Random Access Memory；Hereinafter referred to as: RAM) 30 and/or cache memory 32.Computer equipment 12
It may further include other removable/nonremovable, volatile/non-volatile computer system storage mediums.Only conduct
Citing, storage system 34 can be used for reading and writing immovable, non-volatile magnetic media, and (Fig. 5 do not show, commonly referred to as " hard disk
Driver ").Although being not shown in Fig. 5, the magnetic for reading and writing to removable non-volatile magnetic disk (such as " floppy disk ") can be provided
Disk drive, and to removable anonvolatile optical disk (such as: compact disc read-only memory (Compact Disc Read Only
Memory；Hereinafter referred to as: CD-ROM), digital multi CD-ROM (Digital Video Disc Read Only
Memory；Hereinafter referred to as: DVD-ROM) or other optical mediums) read-write CD drive.In these cases, each driving
Device can be connected by one or more data media interfaces with bus 18.Memory 28 may include that at least one program produces
Product, the program product have one group of (for example, at least one) program module, and it is each that these program modules are configured to perform the application
The function of embodiment.
Program/utility 40 with one group of (at least one) program module 42 can store in such as memory 28
In, such program module 42 includes --- but being not limited to --- operating system, one or more application program, other programs
It may include the realization of network environment in module and program data, each of these examples or certain combination.Program mould
Block 42 usually executes function and/or method in embodiments described herein.
Computer equipment 12 can also be with one or more external equipments 14 (such as keyboard, sensing equipment, display 24
Deng) communication, can also be enabled a user to one or more equipment interact with the computer equipment 12 communicate, and/or with make
The computer equipment 12 any equipment (such as network interface card, the modulatedemodulate that can be communicated with one or more of the other calculating equipment
Adjust device etc.) communication.This communication can be carried out by input/output (I/O) interface 22.Also, computer equipment 12 may be used also
To pass through network adapter 20 and one or more network (such as local area network (Local Area Network；Hereinafter referred to as:
LAN), wide area network (Wide Area Network；Hereinafter referred to as: WAN) and/or public network, for example, internet) communication.Such as figure
Shown in 5, network adapter 20 is communicated by bus 18 with other modules of computer equipment 12.Although should be understood that in Fig. 5 not
It shows, other hardware and/or software module can be used in conjunction with computer equipment 12, including but not limited to: microcode, equipment are driven
Dynamic device, redundant processing unit, external disk drive array, RAID system, tape drive and data backup storage system etc..
Processing unit 16 by the program that is stored in system storage 28 of operation, thereby executing various function application and
Data processing, such as realize the method provided by the embodiments of the present application being filtered to pulse wave.
The embodiment of the present application also provides a kind of non-transitorycomputer readable storage medium, is stored thereon with computer journey
The side provided by the embodiments of the present application being filtered to pulse wave may be implemented in sequence, above-mentioned computer program when being executed by processor
Method.
Above-mentioned non-transitorycomputer readable storage medium can appointing using one or more computer-readable media
Meaning combination.Computer-readable medium can be computer-readable signal media or computer readable storage medium.Computer can
Reading storage medium for example may be-but not limited to-the system of electricity, magnetic, optical, electromagnetic, infrared ray or semiconductor, device
Or device, or any above combination.The more specific example (non exhaustive list) of computer readable storage medium includes:
Electrical connection, portable computer diskette, hard disk, random access memory (RAM), read-only storage with one or more conducting wires
Device (Read Only Memory；Hereinafter referred to as: ROM), erasable programmable read only memory (Erasable
Programmable Read Only Memory；Hereinafter referred to as: EPROM) or flash memory, optical fiber, portable compact disc are read-only deposits
Reservoir (CD-ROM), light storage device, magnetic memory device or above-mentioned any appropriate combination.In this document, computer
Readable storage medium storing program for executing can be any tangible medium for including or store program, which can be commanded execution system, device
Either device use or in connection.
Computer-readable signal media may include in a base band or as carrier wave a part propagate data-signal,
Wherein carry computer-readable program code.The data-signal of this propagation can take various forms, including --- but
It is not limited to --- electromagnetic signal, optical signal or above-mentioned any appropriate combination.Computer-readable signal media can also be
Any computer-readable medium other than computer readable storage medium, which can send, propagate or
Transmission is for by the use of instruction execution system, device or device or program in connection.
The program code for including on computer-readable medium can transmit with any suitable medium, including --- but it is unlimited
In --- wireless, electric wire, optical cable, RF etc. or above-mentioned any appropriate combination.
Can with one or more programming languages or combinations thereof come write for execute the application operation computer
Program code, described program design language include object oriented program language-such as Java, Smalltalk, C++,
It further include conventional procedural programming language-such as " C " language or similar programming language.Program code can be with
It fully executes, partly execute on the user computer on the user computer, being executed as an independent software package, portion
Divide and partially executes or executed on a remote computer or server completely on the remote computer on the user computer.?
It is related in the situation of remote computer, remote computer can pass through the network of any kind --- including local area network (Local
Area Network；Hereinafter referred to as: LAN) or wide area network (Wide Area Network；Hereinafter referred to as: WAN) it is connected to user
Computer, or, it may be connected to outer computer (such as being connected using ISP by internet).
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is contained at least one embodiment or example of the application.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present application, the meaning of " plurality " is at least two, such as two, three
It is a etc., unless otherwise specifically defined.
Any process described otherwise above or method description are construed as in flow chart or herein, and expression includes
It is one or more for realizing custom logic function or process the step of executable instruction code module, segment or portion
Point, and the range of the preferred embodiment of the application includes other realization, wherein can not press shown or discussed suitable
Sequence, including according to related function by it is basic simultaneously in the way of or in the opposite order, Lai Zhihang function, this should be by the application
Embodiment person of ordinary skill in the field understood.
Depending on context, word as used in this " if " can be construed to " ... when " or " when ...
When " or " in response to determination " or " in response to detection ".Similarly, depend on context, phrase " if it is determined that " or " if detection
(condition or event of statement) " can be construed to " when determining " or " in response to determination " or " when the detection (condition of statement
Or event) when " or " in response to detection (condition or event of statement) ".
It should be noted that terminal involved in the embodiment of the present application can include but is not limited to personal computer
(Personal Computer；Hereinafter referred to as: PC), personal digital assistant (Personal Digital Assistant；Below
Referred to as: PDA), radio hand-held equipment, tablet computer (Tablet Computer), mobile phone, MP3 player, MP4 player etc..
In several embodiments provided herein, it should be understood that disclosed systems, devices and methods, it can be with
It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit
It divides, only a kind of logical function partition, there may be another division manner in actual implementation, for example, multiple units or group
Part can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown
Or the mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, device or unit it is indirect
Coupling or communication connection can be electrical property, mechanical or other forms.
It, can also be in addition, each functional unit in each embodiment of the application can integrate in one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of hardware adds SFU software functional unit.
The above-mentioned integrated unit being realized in the form of SFU software functional unit can store and computer-readable deposit at one
In storage media.Above-mentioned SFU software functional unit is stored in a storage medium, including some instructions are used so that a computer
It is each that device (can be personal computer, server or network equipment etc.) or processor (Processor) execute the application
The part steps of embodiment the method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory (Read-
Only Memory；Hereinafter referred to as: ROM), random access memory (Random Access Memory；Hereinafter referred to as: RAM),
The various media that can store program code such as magnetic or disk.
The foregoing is merely the preferred embodiments of the application, not to limit the application, all essences in the application
Within mind and principle, any modification, equivalent substitution, improvement and etc. done be should be included within the scope of the application protection.
Claims (10)
1. the method that a kind of pair of pulse wave is filtered characterized by comprising
Obtain collected pulse wave；
Extremum extracting is carried out to the pulse wave, obtains the extreme point of the pulse wave；
Classify to the extreme point of the pulse wave, obtains the wave crest and trough of the pulse wave；
To each wave crest of the pulse wave, the height of first trough after each wave crest and the wave crest is calculated separately
Difference, and calculate the difference in height of second trough after each wave crest and the wave crest；
It is averaged to the difference in height obtained is calculated, obtains the first mean value of difference in height, and first according to the difference in height is equal
Value setting first threshold；
The extreme point that the absolute value of amplitude in the extreme point of the pulse wave is greater than or equal to the first threshold is filtered out.
2. the method according to claim 1, wherein it is described by the extreme point of the pulse wave amplitude it is absolute
After the extreme point that value is greater than or equal to the first threshold filters out, further includes:
Classify to the remaining extreme point of the pulse wave, obtains the wave crest and trough in the pulse wave residue extreme point；
For each wave crest in the remaining extreme point, first trough after each wave crest and the wave crest is calculated separately
Difference in height, and calculate the difference in height of second trough after each wave crest and the wave crest；
Be averaged to the difference in height obtained is calculated, obtain the second mean value of difference in height, and is arranged the according to second mean value
Two threshold values；
The extreme point that the absolute value of amplitude in the remaining extreme point of the pulse wave is less than the second threshold is filtered out.
3. the method according to claim 1, wherein described according to the first mean value of difference in height setting first
Threshold value includes:
First mean value of the difference in height is multiplied with pre-determined factor, using the product of acquisition as the first threshold.
4. method according to claim 1 to 3, which is characterized in that described to carry out extreme value inspection to the pulse wave
It surveys, before the extreme point for obtaining the pulse wave, further includes:
Noise reduction process is carried out to collected pulse wave by polynomial curve fitting.
5. the device that a kind of pair of pulse wave is filtered characterized by comprising
Module is obtained, for obtaining collected pulse wave；
Detection module, the pulse wave for obtaining to the acquisition module carry out extremum extracting, obtain the extreme value of the pulse wave
Point；
The extreme point of categorization module, the pulse wave for obtaining to the detection module is classified, and the pulse wave is obtained
Wave crest and trough；
Computing module calculates separately each wave crest and first after the wave crest for each wave crest to the pulse wave
The difference in height of a trough, and calculate the difference in height of second trough after each wave crest and the wave crest；
Filter module obtains the first mean value of difference in height, and according to the height for being averaged to the difference in height for calculating acquisition
First threshold is arranged in first mean value of degree difference；The absolute value of amplitude in the extreme point of the pulse wave is greater than or equal to described the
The extreme point of one threshold value filters out.
6. device according to claim 5, which is characterized in that
The categorization module, be also used to be greater than the absolute value of amplitude in the extreme point of the pulse wave in the filter module or
Equal to the first threshold extreme point filter out after, classify to the remaining extreme point of the pulse wave, obtain the arteries and veins
The wave crest and trough fought in wave residue extreme point；
The computing module, is also used to for each wave crest in the remaining extreme point, calculate separately each wave crest with it is described
The difference in height of first trough after wave crest, and calculate the height of second trough after each wave crest and the wave crest
Difference；It is averaged to the difference in height obtained is calculated, obtains the second mean value of difference in height, and according to second mean value setting second
Threshold value；
The filter module, the absolute value for being also used to amplitude in the remaining extreme point by the pulse wave are less than the second threshold
Extreme point filter out.
7. device according to claim 5, which is characterized in that
The filter module is made specifically for being multiplied the first mean value of the difference in height with pre-determined factor with the product of acquisition
For the first threshold.
8. according to device described in claim 5-7 any one, which is characterized in that further include:
Noise reduction module obtains the extreme value of the pulse wave for carrying out extremum extracting to the pulse wave in the detection module
Before point, noise reduction process is carried out to collected pulse wave by polynomial curve fitting.
9. a kind of computer equipment, which is characterized in that including memory, processor and be stored on the memory and can be in institute
The computer program run on processor is stated, when the processor executes the computer program, is realized as in claim 1-4
Any method.
10. a kind of non-transitorycomputer readable storage medium, is stored thereon with computer program, which is characterized in that the meter
The method as described in any in claim 1-4 is realized when calculation machine program is executed by processor.
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 true CN109864705A (en) | 2019-06-11 |
| CN109864705B 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) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110960203A (en) * | 2019-12-13 | 2020-04-07 | 心核心科技(北京)有限公司 | Cardiovascular characteristic parameter detection method and device |
| CN111551952A (en) * | 2020-05-08 | 2020-08-18 | 中国地质大学（武汉） | Extraction method of water depth measurement laser data, water depth measurement method and device |
| WO2021103067A1 (en) * | 2019-11-29 | 2021-06-03 | 苏州润迈德医疗科技有限公司 | Method and system for acquiring extreme points of diastolic pressure and systolic pressure according to invasive blood pressure |
| CN112932428A (en) * | 2021-02-05 | 2021-06-11 | 江西国科美信医疗科技有限公司 | Blood oxygen and heart rate monitoring system and portable wireless remote monitoring device thereof |
| CN113297943A (en) * | 2021-05-18 | 2021-08-24 | 北京远舢智能科技有限公司 | Equipment auxiliary control technology based on mixed reality |
| CN113545753A (en) * | 2021-07-21 | 2021-10-26 | 芯原微电子（成都）有限公司 | Waveform detection method and device, wearable device and computer-readable storage medium |
Citations (4)
| 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 |
| US20180132744A1 (en) * | 2015-04-20 | 2018-05-17 | Vita-Course Technologies Co., Ltd. | Systems and methods for physiological sign analysis |
-
2019
- 2019-01-07 CN CN201910012502.9A patent/CN109864705B/en active Active
Patent Citations (4)
| 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 |
| US20180132744A1 (en) * | 2015-04-20 | 2018-05-17 | Vita-Course Technologies Co., Ltd. | Systems and methods for physiological sign analysis |
Non-Patent Citations (1)
| Title |
|---|
| 丑永新: "动态脉搏信号检测与脉率变异性实时分析方法研究", 《中国博士学位论文全文数据库（电子期刊） 信息科技辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021103067A1 (en) * | 2019-11-29 | 2021-06-03 | 苏州润迈德医疗科技有限公司 | Method and system for acquiring extreme points of diastolic pressure and systolic pressure according to invasive blood pressure |
| CN110960203A (en) * | 2019-12-13 | 2020-04-07 | 心核心科技(北京)有限公司 | Cardiovascular characteristic parameter detection method and device |
| CN111551952A (en) * | 2020-05-08 | 2020-08-18 | 中国地质大学（武汉） | Extraction method of water depth measurement laser data, water depth measurement 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 |
| CN113297943A (en) * | 2021-05-18 | 2021-08-24 | 北京远舢智能科技有限公司 | Equipment auxiliary control technology based on mixed reality |
| CN113297943B (en) * | 2021-05-18 | 2022-12-02 | 北京远舢智能科技有限公司 | Equipment auxiliary control technology based on mixed reality |
| CN113545753A (en) * | 2021-07-21 | 2021-10-26 | 芯原微电子（成都）有限公司 | Waveform detection method and device, wearable device and computer-readable storage medium |
| CN113545753B (en) * | 2021-07-21 | 2024-04-05 | 芯原微电子（成都）有限公司 | Waveform detection method and device, wearable equipment and computer readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109864705B (en) | 2022-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109864705A (en) | The method, apparatus and computer equipment that pulse wave is filtered | |
| US10866637B2 (en) | Gesture classification apparatus and method using EMG signal | |
| Liu et al. | Spectrum analysis of EEG signals using CNN to model patient’s consciousness level based on anesthesiologists’ experience | |
| US10980429B2 (en) | Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time | |
| Fishman et al. | A method for analyzing temporal patterns of variability of a time series from Poincare plots | |
| CN110461215B (en) | Determining health signs using a portable device | |
| CN103251391A (en) | Determining arterial pulse transit time from video images | |
| EP3165158B1 (en) | Signal feature extracting apparatus | |
| 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 | |
| CN108403101A (en) | A kind of model training method and device, the method and device for measuring blood pressure | |
| CN110477865A (en) | A kind of epileptic attack detection device, terminal device and storage medium | |
| CN106659404A (en) | Continuous blood pressure measurement method, apparatus and device | |
| US11596764B2 (en) | Electronic device and method for providing information for stress relief by same | |
| CN108937880A (en) | Wearable device and pulse detection method for pulse detection | |
| KR101977745B1 (en) | Method and apparatus for classifying cardiac arrhythmia | |
| Signorini et al. | Antepartum fetal monitoring through a wearable system and a mobile application | |
| EP3025649B1 (en) | Biometric sound testing device and biometric sound testing method | |
| CN109077721A (en) | atrial fibrillation detection device and storage medium | |
| Kwon et al. | Sinabro: A smartphone-integrated opportunistic electrocardiogram monitoring system | |
| CN110200620A (en) | Multi-modal cardiac data detection method, device, storage medium and terminal device | |
| CN116687409A (en) | Emotion recognition method and system based on digital twin and deep learning | |
| CN116327133A (en) | Multi-physiological index detection method, device and related equipment | |
| JP5911840B2 (en) | Diagnostic data generation device and diagnostic device | |
| US20220211283A1 (en) | Methods for blood pressure calibration selection and modeling methods thereof | |
| CN108937925A (en) | A kind of method, apparatus, equipment and the storage medium of evoked brain potential signal extraction |
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 |