CN107898443A - Dicrotic pulse wave detecting method, device and computer-readable storage medium - Google Patents
Dicrotic pulse wave detecting method, device and computer-readable storage medium Download PDFInfo
- Publication number
- CN107898443A CN107898443A CN201711172968.2A CN201711172968A CN107898443A CN 107898443 A CN107898443 A CN 107898443A CN 201711172968 A CN201711172968 A CN 201711172968A CN 107898443 A CN107898443 A CN 107898443A
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- dicrotic
- signal
- wave
- interest
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The present invention relates to processing of biomedical signals field, specifically, there is provided a kind of dicrotic pulse wave detecting method, device and computer-readable storage medium.This method positions the dicrotic wave position of pulse wave signal by detecting the Local Extremum of pulse wave signal wavelet coefficient on some scale of Mexican hat wavelet transformations.The detection interval of dicrotic wave position is determined using the main ripple wave crest and cycle of each pulse wave cadence signal, can not only reduce calculation amount, moreover it is possible to reduces the influence that other waveforms detect dicrotic wave.The dicrotic wave position of pulse wave signal still keeps extreme point position on its each rank scale based on Mexican hat wavelet transformations, utilize the characteristic, the dicrotic wave trough of pulse wave signal and crest location are positioned by directly detecting local minimum (or minimum value) and the maximum (or maximum) of pulse wave signal wavelet coefficient on some scale of Mexican hat wavelet transformations, it not only can effectively reduce computation complexity, moreover it is possible to effectively improve dicrotic wave Feature point recognition accuracy.
Description
Technical field
The present invention relates to processing of biomedical signals field, in particular to a kind of dicrotic pulse wave detecting method, device and
Computer-readable storage medium.
Background technology
Dicrotic wave is characteristic point position important in pulse wave signal, and dicrotic wave can react the function shape of aorta petal
The states such as condition, blood vessel elasticity and blood flow.The dicrotic wave position of accurate detection pulse signal, for cardiovascular system health shape
The assessment of condition plays the role of important.Compared to the main ripple of pulse wave signal, dicrotic wave is fainter, and is easily subject to external environment condition
With the interference of internal waveform, this brings certain difficulty to dicrotic wave detection.Existing method is to obvious dicrotic wave detection result
Preferably, it is but undesirable to unconspicuous dicrotic wave detection result.
At present, the dicrotic wave detection of pulse wave signal mainly uses relative position method, difference zero-crossing method and wavelet analysis
Modulus maximum.
The main thought of relative position method is:Pulse wave signal pre-processes, and eliminates noise jamming;Accurate detection pulse wave letter
Number main crest value point position;Using main crest value point as reference point, on the right of it in certain time domain scale (0.255~0.605)
Pulse wave ripple signal in find minimum and maximum position as dicrotic wave trough and crest location.This method is using relatively
Position relationship shortens dicrotic wave detection interval, reduces the interference that other waveforms detect dicrotic wave, but the method is not to
Obvious dicrotic wave detection result is undesirable.
The method being combined using relative position method and difference zero crossing detects the main of pulse wave signal dicrotic wave position
Thought is:Synchronous electrocardiosignal and pulse wave signal pretreatment, for eliminating noise jamming;Accurate detection electrocardiosignal T ripples position
Put;Two difference zero crossings are found in front and rear a certain range in the pulse wave signal of T ripple position correspondences, respectively as dicrotic wave
Trough and crest location.This method shortens dicrotic wave detection interval also with relative position relation, but existing for the method
Shortcoming is that process is cumbersome, it is necessary to measure electrocardiosignal;In addition, when the dicrotic wave unobvious of pulse wave, the corresponding difference of pulse wave
Zero crossing may be not present in sub-signal, this prevents difference zero-crossing method from accurately detecting dicrotic wave position.
Using the modulus maximum detection pulse wave signal dicrotic wave position of wavelet analysis, its main thought is:Pulse wave
Signal Pretreatment, for eliminating noise jamming;Wavelet transformation is carried out to pulse wave signal using bior6.8 small echos, in scale j=
4 small echo, which becomes, carries out dicrotic wave wave crest positioning in modulus maximum sequence.This method protrudes the energy of local feature using wavelet transformation
Power, can effectively highlight position of the dicrotic wave in pulse wave, but after wavelet transformation, there are multiple mould poles on wavelet coefficient
Big value pair, accurately filters out the corresponding modulus maximum of dicrotic wave to relatively difficult.
The content of the invention
The purpose of the embodiment of the present invention is to provide a kind of dicrotic pulse wave detecting method, device and computer-readable storage medium, with
Detect the position of dicrotic wave.
To achieve the above object, the technical solution that the embodiment of the present invention uses is as described below:
In a first aspect, an embodiment of the present invention provides a kind of dicrotic pulse wave detecting method, the described method includes:Obtain pulse wave
Signal;Detect the main ripple wave crest of the pulse wave signal;According to the position of the main ripple wave crest, calculate in the pulse wave signal
The cycle of each cadence signal aroused in interest;Determine the first detection time domain model of each cadence signal aroused in interest respectively according to each cycle
Enclose;Wavelet transformation is carried out to the pulse wave signal in the cycle by Mexican-hat small echos, obtains wavelet coefficient;Detection
Amplitude smallest point of the wavelet coefficient in the described first detection time domain scale, corresponds to the arteries and veins by the amplitude smallest point
Trough of the position fought in ripple signal as the dicrotic wave of the cycle corresponding cadence signal aroused in interest;According to the dicrotic wave
Trough determines the second detection time domain scale;Detect amplitude maximum of the wavelet coefficient in the described second detection time domain scale
Point, corresponds to the position in the pulse wave signal as the cycle corresponding cadence signal aroused in interest using the amplitude maximum point
Dicrotic wave wave crest.
Further, the method further includes:Determined to correct time domain model according to the trough of the dicrotic wave and crest location
Enclose;Calculus of differences is carried out to the pulse wave signal in the amendment time domain scale, differential signal is obtained, when the differential signal
When more than zero curve, be incremented by first zero crossing of direction finding towards the time in the starting point of the amendment time domain scale, will be described
First zero crossing that point is incremented by direction towards the time is corresponding as the cycle corresponding to the position in the pulse wave signal
The trough of the dicrotic wave of cadence signal aroused in interest, successively decreases first mistake of direction finding towards the time in the terminal of the amendment time domain scale
Zero point, using the terminal towards the time successively decrease direction first zero crossing correspond to the pulse wave signal in position as institute
State the wave crest of the dicrotic wave of cycle corresponding cadence signal aroused in interest.
Further, further included after the step of acquisition pulse wave signal:Denoising is carried out to the pulse wave signal.
Further, it is described that pulse wave signal progress denoising is included:Baseline drift is removed by wavelet filtering;It is logical
The integral coefficient LP filter for crossing default cutoff frequency removes high-frequency noise.
Further, the cycle for calculating each cadence signal aroused in interest in the pulse wave signal passes through formula Ti=
Pi+1-PiRealize, the first detection time domain scale is:[Pi+0.2Ti,Pi+0.5Ti], wherein, Ti(i=1,2 ..., N-1) be
The cycle of cadence signal aroused in interest, Pi+1When (i=1,2 ..., N) is corresponding for the main ripple wave crest of (i+1) a cadence signal aroused in interest
Between, Pi(i=1,2 ..., N) is the main ripple wave crest corresponding time of i-th cadence signal aroused in interest, wherein, i represent i-th it is aroused in interest
Cadence signal, N represent the sum of cadence signal aroused in interest in whole pulse wave signal.
Further, the second detection time domain scale is:[Di,Pi+0.5Ti], wherein, Di(i=1,2 ..., N) is the
The trough corresponding time of the dicrotic wave of i cadence signals aroused in interest.
Further, the amendment time domain scale is:[Di-0.02Ti,Fi+0.05i], wherein, Fi(i=1,2 ..., N) be
The wave crest corresponding time of the dicrotic wave of i-th of cadence signal aroused in interest.
It is further, described that wavelet transformation is carried out to the pulse wave signal in the cycle by Mexican-hat small echos,
The step of obtaining wavelet coefficient includes:4 scales are carried out to the pulse wave signal in the cycle by Mexican-hat small echos
One Dimension Continuous Wavelet Transform, obtains wavelet coefficient.
Further, when the differential signal is below zero curve, keep the trough of the dicrotic wave and wave crest constant.
Second aspect, the embodiment of the present invention additionally provide a kind of computer-readable storage medium, it is characterised in that the storage medium
Including one group of instruction, when executed, cause at least one processor to perform first aspect of the embodiment of the present invention and provide
The method.
The third aspect, the embodiment of the present invention additionally provide a kind of dicrotic wave detection device, and described device includes:Obtain mould
Block, for obtaining pulse wave signal;First detection module, for detecting the main ripple wave crest of the pulse wave signal;Computing module,
For the position according to the main ripple wave crest, the cycle of each cadence signal aroused in interest in the pulse wave signal is calculated;First model
Determining module is enclosed, for determining the first detection time domain scale of each cadence signal aroused in interest respectively according to each cycle;Become mold changing
Block, for carrying out wavelet transformation to the pulse wave signal in the cycle by Mexican-hat small echos, obtains wavelet coefficient;
Second detection module, for detecting amplitude smallest point of the wavelet coefficient in the described first detection time domain scale, by described in
Amplitude smallest point corresponds to dicrotic wave of the position as the cycle corresponding cadence signal aroused in interest in the pulse wave signal
Trough;Second range determination module, for determining the second detection time domain scale according to the trough of the dicrotic wave;3rd detection
Module, for detecting amplitude maximum point of the wavelet coefficient in the described second detection time domain scale, by the amplitude maximum
Point corresponds to wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in the pulse wave signal.
Further, described device further includes the 3rd range determination module, for the trough and ripple according to the dicrotic wave
Peak determines to correct time domain scale;Correcting module, for carrying out difference fortune to the pulse wave signal in the amendment time domain scale
Calculate, obtain differential signal, when the differential signal is more than zero curve, be incremented by the starting point of the amendment time domain scale towards the time
First zero crossing of direction finding, first zero crossing that the starting point is incremented by direction towards the time are believed corresponding to the pulse wave
Trough of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in number, at the end of the amendment time domain scale
Point successively decreased first zero crossing of direction finding towards the time, and the terminal is corresponded to towards successively decrease first zero crossing in direction of time
Wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in the pulse wave signal.
Further, described device further includes denoising module, for carrying out denoising to the pulse wave signal.
Dicrotic pulse wave detecting method and computer-readable storage medium provided in an embodiment of the present invention, are believed by directly detecting pulse wave
The Local Extremum of wavelet coefficient positions the dicrotic pulse of pulse wave signal number on some scale of Mexican-hat wavelet transformations
Ripple position.Each pulse wave cadence signal dicrotic wave position is determined using the main ripple wave crest and cycle of each pulse wave cadence signal
The detection interval put, can not only reduce calculation amount, moreover it is possible to reduce the influence that other waveforms detect dicrotic wave.Pulse wave signal
Dicrotic wave position still keeps extreme point position on its each rank scale based on Mexican-hat wavelet transformations, using the characteristic,
By the local minimum for directly detecting pulse wave signal wavelet coefficient on some scale of Mexican-hat wavelet transformations
(or minimum value) and maximum (or maximum) position the dicrotic wave trough of pulse wave signal and crest location, not only can be effective
Reduce computation complexity, moreover it is possible to effectively improve dicrotic wave Feature point recognition accuracy.In addition, utilize first-order difference signal zero-crossing
To correct the pulse wave signal dicrotic wave trough and crest location that are detected based on Mexican-hat wavelet transformations, can reduce
Because of the influence that wavelet transformation time shift detects dicrotic pulse wave characteristic point, the accuracy of dicrotic wave characteristic point detection is further improved.This
The dicrotic pulse wave detecting method that inventive embodiments provide can be the cardiovascular system health status based on pulse wave signal characteristic parameter
Assessment provides technical foundation.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, therefore be not construed as pair
The restriction of scope, for those of ordinary skill in the art, without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 is the oscillogram of typical pulse wave signal.
Fig. 2 is the flow chart of dicrotic pulse wave detecting method provided in an embodiment of the present invention.
Fig. 3 is first-order difference signal corresponding when the dicrotic wave of pulse wave signal is obvious and small based on Mexican-hat
Wavelet coefficient under 4 scales of wave conversion.
First-order difference signal corresponding when being the dicrotic wave unobvious of pulse wave signal Fig. 4 and based on Mexican-hat
Wavelet coefficient under 4 scales of wavelet transformation.
Fig. 5 is the flow chart of another embodiment of dicrotic pulse wave detecting method provided in an embodiment of the present invention.
Fig. 6 is the dicrotic pulse wave detecting method inspection provided according to embodiments of the present invention in pulse wave signal dicrotic wave unobvious
The main ripple wave crest for the pulse wave signal measured, dicrotic wave trough, the oscillogram of dicrotic wave crest location.
Fig. 7 is the dicrotic pulse wave detecting method detection provided according to embodiments of the present invention when pulse wave signal dicrotic wave is obvious
The main ripple wave crest of the pulse wave signal gone out, dicrotic wave trough, the oscillogram of dicrotic wave crest location.
Fig. 8 is the arteries and veins that the dicrotic pulse wave detecting method provided according to embodiments of the present invention when there is obvious tidal wave to disturb detects
Fight the main ripple wave crest of ripple signal, dicrotic wave trough, the oscillogram of dicrotic wave crest location.
Fig. 9 is the functional module framework schematic diagram of pulse wave detection device provided in an embodiment of the present invention.
Icon:901- acquisition modules;902- first detection modules;903- computing modules;The first range determination modules of 904-;
905- conversion modules;The second detection modules of 906-;The second range determination modules of 907-;The 3rd detection modules of 908-;909- the 3rd
Range determination module;910- correcting modules;911- denoising modules.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical solution in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
Part of the embodiment of the present invention, instead of all the embodiments.The present invention implementation being usually described and illustrated herein in the accompanying drawings
The component of example can be arranged and designed with a variety of configurations.
Therefore, below the detailed description of the embodiment of the present invention to providing in the accompanying drawings be not intended to limit it is claimed
The scope of the present invention, but be merely representative of the present invention selected embodiment.Based on the embodiments of the present invention, this area is common
Technical staff's all other embodiments obtained without making creative work, belong to the model that the present invention protects
Enclose.
It should be noted that:Similar label and letter represents similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined, then it further need not be defined and explained in subsequent attached drawing in a attached drawing.
Pulse wave signal contains abundant physiology and pathological information, is the important physiological parameter of human body, pulse wave signal
The diagnosis accurately detected for disease and prevention of characteristic point have directive function.Typical pulse wave signal waveform such as Fig. 1 institutes
Show, the rising edge of pulse wave signal includes a main ripple wave crest (P points), and trailing edge includes a dicrotic wave.Dicrotic wave is pulse
Important characteristic point position in ripple signal, it is due in the ventricular ejection later stage that it, which is produced, and the intra-arterial blood that backflows impinges upon suddenly
The aorta petal of closing and be shot back, aortic pressure is slightly increased again, the secondary ascending wave of formation.Dicrotic wave can react
The physiological status such as function status, blood vessel elasticity and the blood flow of aorta petal, its trough (D points as shown in Figure 1) are that heart is received
Contracting and the separation of diastole, key reaction sustainer static pressure emptying time, and its wave crest (F points as shown in Figure 1) can react blood
The size of pipe peripheral resistance.Therefore, the dicrotic wave of pulse wave signal, the assessment for cardiovascular system health status are accurately detected
Have great significance.
For the deficiencies in the prior art, the embodiment of the present invention is faint in order to solve the dicrotic wave of pulse wave, easily
A kind of the problem of being subject to various interference and causing dicrotic wave detection difficult, it is proposed that weight based on Mexican-hat wavelet transformations
Fight wave detecting method, to improve the dicrotic wave accuracy in detection to pulse wave signal, for based on pulse wave signal characteristic parameter
Cardiovascular system health Evaluation provides technical foundation.
Fig. 2 is refer to, is the flow chart of dicrotic pulse wave detecting method provided in an embodiment of the present invention, the dicrotic pulse wave detecting method
Comprise the following steps:
Step S101, obtains pulse wave signal.
The pulse wave signal is the pulse wave signal of the subject gathered by pulse wave signal collecting device, it is preferable that
In the present embodiment, signal sampling rate is more than 200Hz, and the pulse wave signal gathered can be pressure pulse wave signal, can also
It is photoplethysmographic signal.
Step S102, detects the main ripple wave crest of the pulse wave signal.
The interference that the main ripple wave crest of pulse wave signal is subject to is smaller, is easily detected, existing to be based on wavelet method or difference
Detection of the method to the main ripple wave crest of pulse wave signal all has higher accuracy rate.Due to the use of wavelet method or calculus of finite differences detection arteries and veins
The main crest value point of ripple signal of fighting be to those skilled in the art it is known, therefore the present embodiment this is not done it is superfluous
State.It should be noted that before being detected to pulse wave, denoising can be carried out to pulse wave signal, specifically, going
Processing of making an uproar mainly includes the removal to baseline drift and high-frequency noise.Pulse wave signal frequency is mainly in 0~20Hz, and 99%
Energy concentrate in 0~10Hz, the presence of noise can interfere the detection of pulse wave signal characteristic point.In the present embodiment
In, wavelet filtering can be used to remove baseline drift, use cutoff frequency to remove high frequency for the integral coefficient LP filter of 40Hz
Noise.
Step S103, according to the position of main ripple wave crest, calculates the week of each cadence signal aroused in interest in the pulse wave signal
Phase.
The cycle of each cadence signal aroused in interest can pass through formula Ti=Pi+1-PiRealize, Ti(i=1,2 ..., N-1) it is the heart
The cycle of dynamic cadence signal, Pi+1(i=1,2 ..., N) is the main ripple wave crest corresponding time of (i+1) a cadence signal aroused in interest,
Pi(i=1,2 ..., N is the main ripple wave crest corresponding time of i-th of cadence signal aroused in interest.
Step S104, the first detection time domain scale of each cadence signal aroused in interest is determined according to each cycle respectively.
The first detection time domain scale is used for the detection time scope for determining each cadence signal aroused in interest, in order to determine weight
Fight the trough of ripple, in the present embodiment, determine that rule of thumb formula is calculated for the mode of the first detection time domain.
Step S105, carries out wavelet transformation to the pulse wave signal in the cycle by Mexican-hat small echos, obtains
Wavelet coefficient.
Mexican-hat small echos are the second dervatives of Gaussian function, and inventor passes through experimental study, finds pulse wave signal
Dicrotic wave trough and wave crest still kept on each rank scale of its wavelet transformation minimum (or minimum value) and maximum (or
Maximum), and corresponding minimum and maximum position are unique.Further, since wavelet transformation protrudes local feature
Ability, in dicrotic wave position during unobvious, it still corresponds to the extreme value place on each rank scale of wavelet transformation.For example, it please join
According to Fig. 3 and Fig. 4, Fig. 3 is first-order difference signal corresponding when the dicrotic wave of pulse wave signal is obvious and based on Mexican-hat
Wavelet coefficient under 4 scales of wavelet transformation, first-order difference letter corresponding when being the dicrotic wave unobvious of pulse wave signal Fig. 4
Number and 4 scales based on Mexican-hat wavelet transformations under wavelet coefficient.It is from figure 3, it can be seen that obvious in dicrotic wave
When, dicrotic wave wave trough position corresponds to the minimum point position of differential signal zero crossing, wavelet coefficient, and dicrotic wave crest location corresponds to
The maximum point position of differential signal zero crossing, wavelet coefficient.Figure 4, it is seen that in dicrotic wave unobvious, dicrotic wave
Zero crossing is not present on the differential signal of position correspondence, but dicrotic wave trough and crest location remain to the minimum of corresponding wavelet coefficient
Value point and maximum point position.Therefore, using Mexican-hat Wavelet Properties, by the extreme point for directly detecting wavelet transformation
To detect dicrotic wave position, the accuracy of detection of pulse wave signal dicrotic wave position is greatly improved, while reduces detection algorithm
Complexity.
It should be noted that Mexican-hat wavelet transformations, the small echo being not limited under 4 scales are carried out to pulse wave signal
Decompose, suitable decomposition scale can be selected according to experimental result, such as can also be the wavelet decomposition under 6 scales, in this implementation
In example, the One Dimension Continuous Wavelet Transform of 4 scales is carried out to the pulse wave signal in the cycle by Mexican-hat small echos,
Obtain wavelet coefficient s_cwti(n), n is wavelet coefficient s_cwti(n) sampling number included in.
Step S106, amplitude smallest point of the detection wavelet coefficient in the described first detection time domain scale are minimum by amplitude
Point corresponds to trough of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in the pulse wave signal.
In the present embodiment, the first detection time domain scale is:[Pi+0.2Ti,Pi+0.5Ti], it is necessary to explanation, the scope
For empirical value, in practical applications, can appropriate adjustment determine suitable scope.
Step S107, the second detection time domain scale is determined according to the trough of dicrotic wave.
The second detection time domain scale is used for the wave crest for determining dicrotic wave.The second detection time domain scale is used to determine dicrotic pulse
The detection time scope of each cadence signal aroused in interest of ripple, in order to determine the wave crest of dicrotic wave, in the present embodiment, determines second
Rule of thumb formula is calculated the mode of detection time domain.
Step S108, amplitude maximum point of the detection wavelet coefficient in the second detection time domain scale, by amplitude maximum point pair
Wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in pulse wave signal described in Ying Yu.
In the present embodiment, the second detection time domain scale is:[Di,Pi+0.5Ti], wherein, Di(i=1,2 ..., N) is the
The trough corresponding time of the dicrotic wave of i cadence signals aroused in interest.It should be noted that the scope is empirical value, actually should
In, can appropriate adjustment determine suitable scope.
Pass through above step, it is possible to achieve the accurate positionin to dicrotic wave signal, is believed using each pulse wave beat aroused in interest
Number main ripple wave crest and cycle determine the detection interval of each pulse wave cadence signal dicrotic wave position, can not only reduce calculating
Amount, moreover it is possible to reduce the influence that other waveforms detect dicrotic wave.In addition, the dicrotic wave position of pulse wave signal is based at it
Extreme point position is still kept on each rank scale of Mexican-hat wavelet transformations, using the characteristic, by directly detecting pulse wave
The local minimum (or minimum value) and maximum of signal wavelet coefficient on some scale of Mexican-hat wavelet transformations
(or maximum) positions the dicrotic wave trough of pulse wave signal and crest location, not only can effectively reduce computation complexity, also
Dicrotic wave Feature point recognition accuracy can be effectively improved.
Inventor passes through experimental study, finds the true wave trough position of dicrotic wave and the local minimum position of wavelet coefficient
It is sufficiently close to or a slightly biased left side, the local maximum position of dicrotic wave true peaks position and wavelet coefficient is sufficiently close to or slightly biased
The right side, when with error, it is necessary to be modified to the position of dicrotic wave.Fig. 5 is refer to, it is provided in an embodiment of the present invention
Dicrotic pulse wave detecting method further includes:
Step S109, determines to correct time domain scale according to the trough of dicrotic wave and crest location.
Correct time domain scale to be used to the dicrotic wave wave crest primarily determined that and the position of trough are tested and corrected, this is repaiied
Positive time domain scale can be rule of thumb derived from.In the present embodiment, correcting time domain scale is:[Di-0.02Ti,Fi+
0.05Ti], wherein, Fi(i=1,2 ..., N) for the dicrotic wave of i-th of cadence signal aroused in interest the wave crest corresponding time.Need
Bright, which is empirical value, in practical applications, can appropriate adjustment determine suitable scope.
Step S110, carries out calculus of differences to the pulse wave signal in the range of the positive time domain of repairing, obtains differential signal.
To in [Di-0.02Ti,Fi+0.05Ti] pulse wave signal in time domain scale carries out first-order difference computing, obtain one
Jump sub-signal s_difi(n), n represents the first-order difference signal s_difi(n) sampling number included in.
Whether step S111, judge differential signal more than zero curve.If more than zero curve, step S112 is performed, if
Not more than zero curve, then keep the trough of the dicrotic wave and wave crest constant, end step.
Step S112, the starting point of the positive time domain scale of repairing are incremented by first zero crossing of direction finding towards the time, described will rise
First zero crossing that point is incremented by direction towards the time is corresponding aroused in interest as the cycle corresponding to the position in pulse wave signal
The trough of the dicrotic wave of cadence signal.
Step S113, the terminal of the positive time domain scale of repairing successively decreased first zero crossing of direction finding towards the time, by the end
Point towards the time successively decrease direction first zero crossing correspond to pulse wave signal in position it is corresponding aroused in interest as the cycle
The wave crest of the dicrotic wave of cadence signal.
If s_difi(n) below zero curve, that is, s_difi(n) maximum is both less than equal to 0, then dicrotic wave
Position remains unchanged, that is, i-th of cadence signal dicrotic wave wave trough position D aroused in interest after correctingi, i-th of cadence signal weight after amendment
Fight ripple crest location Fi.Otherwise, the positive time domain scale [D of repairingi-0.02Ti,Fi+0.05Ti] in, starting point is turned right, and (time is incremented by side
To) find s_difi(n) first zero crossing, the position of the zero crossing correspondence in pulse wave signal is corresponding as the cycle
The trough of the revised dicrotic wave of cadence signal aroused in interest;Terminal turn left (time successively decrease direction) find s_difi(n) first mistake
Zero point, using the position of the zero crossing correspondence in pulse wave signal as the cycle corresponding revised dicrotic pulse of cadence signal aroused in interest
The wave crest of ripple.
Fig. 6~Fig. 8 is refer to, wherein, Fig. 6 is shown to be implemented according to the present invention in pulse wave signal dicrotic wave unobvious
Main ripple wave crest, dicrotic wave trough, the dicrotic wave crest location for the pulse wave signal that the dicrotic pulse wave detecting method that example provides detects.
As can be seen that by being detected to extreme point on corresponding wavelet coefficient, it can accurately detect dicrotic wave position in pulse wave signal, say
The bright embodiment of the present invention has dicrotic wave preferable detection result in dicrotic wave unobvious.
Fig. 7 shows the dicrotic pulse wave detecting method provided according to embodiments of the present invention when pulse wave signal dicrotic wave is obvious
Main ripple wave crest, dicrotic wave trough, the dicrotic wave crest location of the pulse wave signal detected.It is it can be seen from the figure that of the invention
Embodiment is high to dicrotic wave Detection accuracy when dicrotic wave is obvious.
Fig. 8 shows that the dicrotic pulse wave detecting method provided according to embodiments of the present invention when there is obvious tidal wave to disturb detects
Pulse wave signal main ripple wave crest, dicrotic wave trough, dicrotic wave crest location.The tidal wave of pulse wave signal is main ripple and dicrotic pulse
Wave crest between ripple.Under normal circumstances, tidal wave is fainter, but when tidal wave is obvious, it may be to the detection of dicrotic wave
Interfere.By setting suitable time domain scale, method provided in an embodiment of the present invention is to there are remained to during the interference of obvious tidal wave
Accurate detection dicrotic wave position.
The embodiment of the present invention also provides a kind of computer-readable storage medium, which includes one group of instruction, when execution institute
When stating instruction, at least one processor is caused to perform the above-mentioned dicrotic pulse wave detecting method of the embodiment of the present invention.
Fig. 9 is refer to, the embodiment of the present invention additionally provides a kind of dicrotic wave detection device 900, the dicrotic wave detection device
900 can be applied to terminal device, such as PC (personal computer, PC), tablet computer, personal digital assistant
(personal digital assistant, PDA), mobile internet surfing equipment (mobile Internet device, MID) etc..
The dicrotic wave detection device 900 includes acquisition module 901, first detection module 902, computing module 903, the first scope and determines mould
Block 904, conversion module 905, the second detection module 906, the second range determination module 907, the 3rd detection module 908, the 3rd model
Enclose determining module 909, correcting module 910 and denoising module 911.
Wherein, acquisition module 901 is used to obtain pulse wave signal.In the present embodiment, acquisition module 901 can be used for performing
Step S101.
First detection module 902 is used for the main ripple wave crest for detecting the pulse wave signal.In the present embodiment, the first detection
Module 902 can be used for performing step S102.
Computing module 903 is used for the position according to the main ripple wave crest, calculates each section aroused in interest in the pulse wave signal
Clap the cycle of signal.In the present embodiment, computing module 903 can be used for performing step S103.
First range determination module 904 is used for the first detection for determining each cadence signal aroused in interest respectively according to each cycle
Time domain scale.In the present embodiment, the first range determination module 904 can be used for performing step S104.
Conversion module 905 is used to carry out small echo change to the pulse wave signal in the cycle by Mexican-hat small echos
Change, obtain wavelet coefficient.In the present embodiment, conversion module 905 can be used for performing step S105.
Second detection module 906 is used to detect amplitude minimum of the wavelet coefficient in the described first detection time domain scale
Point, corresponds to the position in the pulse wave signal as the cycle corresponding cadence signal aroused in interest using the amplitude smallest point
Dicrotic wave trough.In the present embodiment, the second detection module 906 can be used for performing step S106.
Second range determination module 907 is used to determine the second detection time domain scale according to the trough of the dicrotic wave.At this
In embodiment, the second range determination module 907 can be used for performing step S107.
3rd detection module 908 is used to detect amplitude maximum of the wavelet coefficient in the described second detection time domain scale
Point, corresponds to the position in the pulse wave signal as the cycle corresponding cadence signal aroused in interest using the amplitude maximum point
Dicrotic wave wave crest.In the present embodiment, the 3rd detection module 908 can be used for performing step S108.
3rd range determination module 909 is used to be determined to correct time domain scale according to the trough and wave crest of the dicrotic wave.
In the present embodiment, the 3rd range determination module 909 can be used for performing step S109.
Correcting module 910 is used to carry out calculus of differences to the pulse wave signal in the amendment time domain scale, obtains difference
Sub-signal, when the differential signal is more than zero curve, is incremented by direction finding towards the time in the starting point of the amendment time domain scale
First zero crossing, first zero crossing that the starting point is incremented by direction towards the time correspond to the position in the pulse wave signal
Put the trough of the dicrotic wave as the cycle corresponding cadence signal aroused in interest, the amendment time domain scale terminal towards the time
Successively decrease first zero crossing of direction finding, and the terminal is corresponded to the pulse towards successively decrease first zero crossing in direction of time
Wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in ripple signal.In the present embodiment, correct
Module 910 can be used for performing step S111~step S113.
Denoising module 911 is used to carry out denoising to the pulse wave signal.
Since the principle of each step being described in preceding method embodiment, dicrotic pulse provided in this embodiment
The modules of ripple detection device 900 can be used for performing each step, and the principle of modules is not repeated herein.
In conclusion an embodiment of the present invention provides a kind of dicrotic pulse wave detecting method, device and computer-readable storage medium, sheet
The dicrotic pulse wave detecting method that embodiment provides is by directly detecting some ruler of pulse wave signal in Mexican-hat wavelet transformations
The Local Extremum of wavelet coefficient positions the dicrotic wave position of pulse wave signal on degree.Utilize each pulse wave cadence signal
Main ripple wave crest and cycle determine the detection interval of each pulse wave cadence signal dicrotic wave position, can not only reduce calculation amount,
The influence that other waveforms detect dicrotic wave can also be reduced.The dicrotic wave position of pulse wave signal is based on Mexican-hat at it
Extreme point position is still kept on each rank scale of wavelet transformation, using the characteristic, is existed by directly detecting pulse wave signal
The local minimum (or minimum value) of wavelet coefficient and maximum on some scale of Mexican-hat wavelet transformations (or it is maximum
Value) the dicrotic wave trough of pulse wave signal and crest location are positioned, it not only can effectively reduce computation complexity, moreover it is possible to effectively carry
High dicrotic wave Feature point recognition accuracy.In addition, corrected using first-order difference signal zero-crossing based on Mexican-hat small echos
The pulse wave signal dicrotic wave trough and crest location that change detection comes out, can reduce because wavelet transformation time shift is to dicrotic pulse wave characteristic
The influence of point detection, further improves the accuracy of dicrotic wave characteristic point detection.Dicrotic wave detection provided in an embodiment of the present invention
Method can provide technical foundation for the cardiovascular system health Evaluation based on pulse wave signal characteristic parameter.
In embodiment provided herein, it should be understood that disclosed method, can also be by another way
Realize.Embodiments described above is only schematical, for example, the flow chart in attached drawing show according to the present invention it is more
Architectural framework in the cards, function and the operation of the method and computer program product of a embodiment.At this point, flow chart
Or each square frame in block diagram can represent a part for a module, program segment or code, the module, program segment or code
A part include one or more be used for realization as defined in logic function executable instructions.It should also be noted that in some works
Function in the implementation of replacement, to be marked in square frame can also be with different from the order marked in attached drawing generation.Example
Such as, two continuous square frames can essentially perform substantially in parallel, they can also be performed in the opposite order sometimes, this according to
Depending on involved function.It is also noted that each square frame and block diagram and/or flow chart in block diagram and/or flow chart
In square frame combination, function or the dedicated hardware based system of action can be realized as defined in execution, Huo Zheke
Realized with combination with specialized hardware and computer instruction.
Based on such understanding, part that technical scheme substantially in other words contributes the prior art or
The part of person's technical solution can be embodied in the form of software product, which is stored in a storage
In medium, including some instructions are used so that a computer equipment (can be that personal computer, server, or network are set
It is standby etc.) perform all or part of step of each embodiment the method for the present invention.And foregoing storage medium includes:USB flash disk, shifting
Dynamic hard disk, read-only storage (ROM, Read-Only Memory), random access memory (RAM, Random Access
Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.It should be noted that herein, such as the
One and second or the like relational terms are used merely to distinguish one entity or operation from another entity or operation,
Without necessarily requiring or implying there are any actual relationship or order between these entities or operation.Moreover, art
Language " comprising ", "comprising" or any other variant thereof is intended to cover non-exclusive inclusion so that including it is a series of will
Process, method, article or the equipment of element not only include those key elements, but also including other elements that are not explicitly listed,
Either further include as elements inherent to such a process, method, article, or device.In the absence of more restrictions,
The key element limited by sentence "including a ...", it is not excluded that in the process including the key element, method, article or equipment
In also there are other identical element.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the invention, for the skill of this area
For art personnel, the invention may be variously modified and varied.Within the spirit and principles of the invention, that is made any repaiies
Change, equivalent substitution, improvement etc., should all be included in the protection scope of the present invention.It should be noted that:Similar label and letter exists
Similar terms is represented in following attached drawing, therefore, once being defined in a certain Xiang Yi attached drawing, is then not required in subsequent attached drawing
It is further defined and is explained.
Claims (13)
- A kind of 1. dicrotic pulse wave detecting method, it is characterised in that the described method includes:Obtain pulse wave signal;Detect the main ripple wave crest of the pulse wave signal;According to the position of the main ripple wave crest, the cycle of each cadence signal aroused in interest in the pulse wave signal is calculated;Determine the first detection time domain scale of each cadence signal aroused in interest respectively according to each cycle;Wavelet transformation is carried out to the pulse wave signal in the cycle by Mexican-hat small echos, obtains wavelet coefficient;Amplitude smallest point of the wavelet coefficient in the described first detection time domain scale is detected, the amplitude smallest point is corresponded to Trough of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in the pulse wave signal;Second detection time domain scale is determined according to the trough of the dicrotic wave;Amplitude maximum point of the wavelet coefficient in the described second detection time domain scale is detected, the amplitude maximum point is corresponded to Wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest in the pulse wave signal.
- 2. dicrotic pulse wave detecting method according to claim 1, it is characterised in that the method further includes:Determined to correct time domain scale according to the trough of the dicrotic wave and crest location;Calculus of differences is carried out to the pulse wave signal in the amendment time domain scale, obtains differential signal, when the difference is believed When number more than zero curve, be incremented by first zero crossing of direction finding towards the time in the starting point of the amendment time domain scale, described in general First zero crossing that starting point is incremented by direction towards the time is corresponded to corresponding to the position in the pulse wave signal as the cycle Cadence signal aroused in interest dicrotic wave trough, successively decrease direction finding first towards the time in the terminal of the amendment time domain scale Zero crossing, using the terminal towards the time successively decrease direction first zero crossing correspond to the pulse wave signal in position as The wave crest of the dicrotic wave of cycle corresponding cadence signal aroused in interest.
- 3. dicrotic pulse wave detecting method according to claim 1, it is characterised in that it is described acquisition pulse wave signal the step of it After further include:Denoising is carried out to the pulse wave signal.
- 4. dicrotic pulse wave detecting method according to claim 3, it is characterised in that described to be gone to the pulse wave signal Make an uproar including:Baseline drift is removed by wavelet filtering;High-frequency noise is removed by the integral coefficient LP filter for presetting cutoff frequency.
- 5. dicrotic pulse wave detecting method according to claim 2, it is characterised in that described to calculate in the pulse wave signal often The cycle of a cadence signal aroused in interest passes through formula Ti=Pi+1-PiRealize, the first detection time domain scale is:[Pi+0.2Ti,Pi+ 0.5Ti], wherein, Ti(i=1,2 ..., N-1) be cadence signal aroused in interest cycle, Pi+1(i=1,2 ..., it is (i+1) a heart The main ripple wave crest corresponding time of dynamic cadence signal, Pi(i=1,2 ..., N) for i-th of cadence signal aroused in interest main ripple wave crest pair The time answered, wherein, i represents i-th of cadence signal aroused in interest, and N represents the sum of cadence signal aroused in interest in whole pulse wave signal.
- 6. dicrotic pulse wave detecting method according to claim 5, it is characterised in that described second, which detects time domain scale, is:[Di, Pi+0.5Ti], wherein, Di(i=1,2 ..., N) for the dicrotic wave of i-th of cadence signal aroused in interest the trough corresponding time.
- 7. dicrotic pulse wave detecting method according to claim 6, it is characterised in that it is described amendment time domain scale be:[Di- 0.02Ti,Fi+0.05i], wherein, Fi(i=1,2 ..., N) for the wave crest of the dicrotic wave of i-th of cadence signal aroused in interest it is corresponding when Between.
- 8. according to claim 1-7 any one of them dicrotic pulse wave detecting methods, it is characterised in that described to pass through Mexican- Hat small echos carry out wavelet transformation to the pulse wave signal in the cycle, and the step of obtaining wavelet coefficient includes:The One Dimension Continuous Wavelet Transform of 4 scales is carried out to the pulse wave signal in the cycle by Mexican-hat small echos, is obtained To wavelet coefficient.
- 9. according to claim 5-7 any one of them dicrotic pulse wave detecting methods, it is characterised in that when the differential signal is zero When below line, keep the trough of the dicrotic wave and wave crest constant.
- A kind of 10. computer-readable storage medium, it is characterised in that the storage medium includes one group of instruction, when executed, At least one processor perform claim is caused to require the method described in any claims of 1-9.
- 11. a kind of dicrotic wave detection device, it is characterised in that described device includes:Acquisition module, for obtaining pulse wave signal;First detection module, for detecting the main ripple wave crest of the pulse wave signal;Computing module, for the position according to the main ripple wave crest, calculates each cadence signal aroused in interest in the pulse wave signal Cycle;First range determination module, for determining the first detection time domain model of each cadence signal aroused in interest respectively according to each cycle Enclose;Conversion module, for carrying out wavelet transformation to the pulse wave signal in the cycle by Mexican-hat small echos, obtains Wavelet coefficient;Second detection module, will for detecting amplitude smallest point of the wavelet coefficient in the described first detection time domain scale The amplitude smallest point corresponds to weight of the position as the cycle corresponding cadence signal aroused in interest in the pulse wave signal Fight the trough of ripple;Second range determination module, for determining the second detection time domain scale according to the trough of the dicrotic wave;3rd detection module, will for detecting amplitude maximum point of the wavelet coefficient in the described second detection time domain scale The amplitude maximum point corresponds to weight of the position as the cycle corresponding cadence signal aroused in interest in the pulse wave signal Fight the wave crest of ripple.
- 12. dicrotic wave detection device according to claim 11, it is characterised in that it is true that described device further includes the 3rd scope Cover half block, determines to correct time domain scale for the trough according to the dicrotic wave and wave crest;Correcting module, for carrying out calculus of differences to the pulse wave signal in the amendment time domain scale, obtains differential signal, When the differential signal is more than zero curve, it is incremented by first mistake of direction finding towards the time in the starting point of the amendment time domain scale Zero point, first zero crossing that the starting point is incremented by direction towards the time correspond to the position in the pulse wave signal as institute The trough of the dicrotic wave of cycle corresponding cadence signal aroused in interest is stated, is successively decreased direction towards the time in the terminal of the amendment time domain scale Find first zero crossing, by the terminal towards the time successively decrease direction first zero crossing correspond to the pulse wave signal in Wave crest of the position as the dicrotic wave of the cycle corresponding cadence signal aroused in interest.
- 13. dicrotic wave detection device according to claim 11, it is characterised in that described device further includes denoising module, For carrying out denoising to the pulse wave signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711172968.2A CN107898443B (en) | 2017-11-21 | 2017-11-21 | Method and device for detecting counterpulsation wave and computer storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711172968.2A CN107898443B (en) | 2017-11-21 | 2017-11-21 | Method and device for detecting counterpulsation wave and computer storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107898443A true CN107898443A (en) | 2018-04-13 |
CN107898443B CN107898443B (en) | 2020-11-24 |
Family
ID=61847151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711172968.2A Active CN107898443B (en) | 2017-11-21 | 2017-11-21 | Method and device for detecting counterpulsation wave and computer storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107898443B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110477889A (en) * | 2018-10-08 | 2019-11-22 | 合云力健康科技(北京)有限公司 | The method and apparatus for identifying radial artery wave shape reflection point |
CN110801210A (en) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | Pulse wave signal filtering method and device, readable medium and electronic equipment |
CN111513689A (en) * | 2020-05-27 | 2020-08-11 | 福州数据技术研究院有限公司 | Method for identifying aorta opening characteristic points in pulse waves |
CN113836995A (en) * | 2021-08-09 | 2021-12-24 | 新绎健康科技有限公司 | Age group identification method and device |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999026529A1 (en) * | 1997-11-20 | 1999-06-03 | Seiko Epson Corporation | Pulse wave diagnostic apparatus, blood pressure monitor, pulse wave shape monitor and pharmacologic effect monitor |
CN1792319A (en) * | 2005-11-03 | 2006-06-28 | 浙江大学 | Automatic testing method for traditional Chinese medical pulse manifestation characteristics parameter |
CN102293639A (en) * | 2011-06-30 | 2011-12-28 | 芜湖圣美孚科技有限公司 | Pulse condition signal time domain feature extraction method |
CN102440768A (en) * | 2010-10-13 | 2012-05-09 | 兰州理工大学 | Method for extracting pulse waveform feature points |
CN102551687B (en) * | 2012-01-18 | 2014-02-05 | 华韵之 | Extraction method of pulse signal feature points based on second-generation wavelets |
CN105286815A (en) * | 2015-11-02 | 2016-02-03 | 重庆大学 | Pulse wave signal feature point detection method based on waveform time domain features |
CN106974631A (en) * | 2017-03-17 | 2017-07-25 | 深圳大学 | A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal |
-
2017
- 2017-11-21 CN CN201711172968.2A patent/CN107898443B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999026529A1 (en) * | 1997-11-20 | 1999-06-03 | Seiko Epson Corporation | Pulse wave diagnostic apparatus, blood pressure monitor, pulse wave shape monitor and pharmacologic effect monitor |
CN1792319A (en) * | 2005-11-03 | 2006-06-28 | 浙江大学 | Automatic testing method for traditional Chinese medical pulse manifestation characteristics parameter |
CN102440768A (en) * | 2010-10-13 | 2012-05-09 | 兰州理工大学 | Method for extracting pulse waveform feature points |
CN102293639A (en) * | 2011-06-30 | 2011-12-28 | 芜湖圣美孚科技有限公司 | Pulse condition signal time domain feature extraction method |
CN102551687B (en) * | 2012-01-18 | 2014-02-05 | 华韵之 | Extraction method of pulse signal feature points based on second-generation wavelets |
CN105286815A (en) * | 2015-11-02 | 2016-02-03 | 重庆大学 | Pulse wave signal feature point detection method based on waveform time domain features |
CN106974631A (en) * | 2017-03-17 | 2017-07-25 | 深圳大学 | A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110477889A (en) * | 2018-10-08 | 2019-11-22 | 合云力健康科技(北京)有限公司 | The method and apparatus for identifying radial artery wave shape reflection point |
CN110477889B (en) * | 2018-10-08 | 2021-11-02 | 合肥伊阳健康科技有限公司 | Method and device for identifying reflection point of radial artery pressure waveform |
CN110801210A (en) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | Pulse wave signal filtering method and device, readable medium and electronic equipment |
CN111513689A (en) * | 2020-05-27 | 2020-08-11 | 福州数据技术研究院有限公司 | Method for identifying aorta opening characteristic points in pulse waves |
CN113836995A (en) * | 2021-08-09 | 2021-12-24 | 新绎健康科技有限公司 | Age group identification method and device |
CN113836995B (en) * | 2021-08-09 | 2024-04-02 | 新绎健康科技有限公司 | Age group identification method and device |
Also Published As
Publication number | Publication date |
---|---|
CN107898443B (en) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107898443A (en) | Dicrotic pulse wave detecting method, device and computer-readable storage medium | |
US10478129B2 (en) | Methods for improving response time, robustness and user comfort in continuous estimation of biophysiological rates | |
Sayadi et al. | Model-based fiducial points extraction for baseline wandered electrocardiograms | |
CN103405227B (en) | Double-layer morphological filter based electrocardiosignal preprocessing method | |
CN106974631A (en) | A kind of blood pressure measuring method and device based on pulse waveform and electrocardiosignal | |
CN101919704B (en) | Heart sound signal positioning and segmenting method | |
CN107708531A (en) | The method that physiological parameter is determined from physiological data | |
CN111528821A (en) | Method for identifying characteristic points of counterpulsation waves in pulse waves | |
Rezk et al. | An algebraic derivative-based method for R wave detection | |
WO2022237222A1 (en) | Pulse wave acquisition method and system based on signal-to-noise ratio improvement technology, and storage medium | |
CN104173030A (en) | Pulse wave starting point real-time detection method resisting waveform change interference and application thereof | |
Xiaomeng | A NEW real-time ECG R-wave detection algorithm | |
CN106137184B (en) | Electrocardiosignal QRS complex detection method based on wavelet transformation | |
Fedotov et al. | Effective QRS-detector based on Hilbert transform and adaptive thresholding | |
Yang et al. | Removal of pulse waveform baseline drift using cubic spline interpolation | |
KR101704491B1 (en) | Heart rate measufing apparatus and method, recording medium for performing the method | |
Tan et al. | EMD-based electrocardiogram delineation for a wearable low-power ECG monitoring device | |
KR20110094103A (en) | Method of deriving central aortic systolic pressure values and method for analysing an arterial dataset to derive the same | |
Lin et al. | A characteristic filtering method for pulse wave signal quality assessment | |
CN110269642B (en) | Doppler heart rate estimation method based on fractional Fourier transform and wavelet transform | |
CN109620198B (en) | Cardiovascular index detection and model training method and device | |
Uchaipichat et al. | Development of QRS detection using short-time fourier transform based technique | |
CN110151156B (en) | Fetal heart rate extraction method and system based on in-window automatic weighted average algorithm | |
CN104352234A (en) | Detection method for peak singular point of physiologic electric signal | |
CN111616695A (en) | Heart rate acquisition method, device, system and medium |
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 |