CN105054926A - Electrocardiosignal feature information extraction method and device - Google Patents
Electrocardiosignal feature information extraction method and device Download PDFInfo
- Publication number
- CN105054926A CN105054926A CN201510172501.2A CN201510172501A CN105054926A CN 105054926 A CN105054926 A CN 105054926A CN 201510172501 A CN201510172501 A CN 201510172501A CN 105054926 A CN105054926 A CN 105054926A
- Authority
- CN
- China
- Prior art keywords
- point
- ripple
- peak
- electrocardiosignal
- maximum
- 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.)
- Pending
Links
Abstract
The present invention discloses an electrocardiosignal feature information extraction method. The method comprises acquiring a to-be-processed electrocardiosignal and determinging the range of a QRS complex; within the range of the QRS complex, determining the position of an R peak and determining positions of a Q peak and an S peak according to the range of the QRS complex and the R peak; and extracting relative information of a T wave and a P wave between the QRS complex and an adjacent QRS complex, the determined QRS complex, the position of the R peak, the positions of the Q peak and the S peak and the extracted relative information of the T wave and the P wave being feature information of the to-be-processed electrocardiosignal, and the relative information including wave crest positions of the T wave and the P wave and starting points and stop points of the T wave and the P wave. An an electrocardiosignal feature information extraction device is also disclosed. The acquisition of the electrocardiosignal feature information is well-adapted, and a calculated amount of electrocardiosignal feature information extraction is reduced.
Description
Technical field
The present invention relates to medicine technology field, particularly relate to extracting method and the device of electrocardiosignal characteristic information.
Background technology
Heart disease is one of the highest disease of current sickness rate and fatality rate, and cardiopathic control and diagnosis have become the major issue of current medical circle.At present, ECG techniques is easy with its operational approach, to patient without advantages such as intrusion infringements, by the important means as heart disease diagnosis.
The diagnostic method of routine electrocardiogram is: gather a bit of ecg wave form, carries out interpreting blueprints diagnose by doctor.This diagnostic method, owing to only acquiring the ecg wave form of a bit of time, truly can not reflect the situation of heart.And traditional Holter ECG detecting technology (i.e. " ambulatory electrocardiogram " technology), although acquire patient's electrocardio-activity overall process of 24 hours, but still main dependence manually carries out off-line diagnosis, doctor often needs to read a large amount of ECG data, and this causes serious impact to the feasibility of diagnosis and accuracy.For the problems referred to above, need automatically to extract electrocardiosignal feature, need the data gone out according to Automatic Feature Extraction just can carry out judging and analyzing according to electrocardiogram automatic diagnosis heart condition.
At present, the method for electrocardiosignal feature extraction is broadly divided into three major types, and a class is morphologic detection technology, utilize the characteristic parameters such as slope to detect, advantage is that amount of calculation is little, real-time, but shortcoming is also very remarkable, and cannot tackle the signal intensity of more susceptible condition complexity, universality is poor; The modern signal processing technology of one class to be the method such as wavelet transformation, empirical mode decomposition be representative, detect after rationally converting electrocardiosignal, advantage is that universality is comparatively strong, but amount of calculation is large; Also have a class to be take neutral net, support vector machine (SVM) etc. as the machine learning algorithm of representative, to detecting after algorithm parameter reasonable training, advantage be universality more by force, but need mass data to train, and amount of calculation is large.
To sum up, the mode of existing electrocardiosignal feature extraction cannot be accomplished adaptable while, reduces the amount of calculation of electrocardiosignal.
Foregoing, only for auxiliary understanding technical scheme of the present invention, does not represent and admits that foregoing is prior art.
Summary of the invention
Main purpose of the present invention is the extracting method and the device that provide a kind of electrocardiosignal characteristic information, and the mode being intended to solve existing electrocardiosignal feature extraction cannot be accomplished adaptable while, reduces the problem of the amount of calculation of electrocardiosignal.
For achieving the above object, the extracting method of a kind of electrocardiosignal characteristic information provided by the invention, comprising:
Obtain pending electrocardiosignal, determine the scope of QRS wave group in described pending electrocardiosignal;
Within the scope of QRS wave group, determine R peak position, and determine the position at Q peak and S peak according to QRS wave group scope and determined R peak position;
Between described QRS wave group and adjacent QRS wave group, extract the crest relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.
Preferably, the relevant information step of described extraction T ripple and P ripple comprises:
In the both sides at T wave-wave peak and P wave-wave peak, extracted the starting point and ending point of T ripple, P ripple by default partial transformation method.
Preferably, describedly determine that the step of the scope of QRS wave group in described pending electrocardiosignal comprises:
Empirical mode decomposition is carried out to described pending electrocardiosignal, obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
According to described single order high frequency intrinsic mode functions and moving average filter function, determine the characteristic signal of described pending electrocardiosignal, determine the scope of described QRS wave group according to the scope of described characteristic signal.
Preferably, describedly determine R peak position, and determine that the step of the position at Q peak and S peak comprises according to QRS wave group scope and determined R peak position:
Obtain the maximum point in described QRS wave group;
When described maximum point is one, determine described R peak position according to described maximum point;
When described maximum point is multiple, obtain the maximum point that in maximum point preset range described in each, slope characteristics is maximum, determine described R peak position according to obtained maximum point;
Within the scope of QRS wave group starting point and R peak, obtain minimum point as the first minimum point;
When described first minimum point is one, determine described Q peak position according to described first minimum point;
When described first minimum point is multiple, obtain the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum, determine described Q peak position according to the minimum point that slope characteristics in the first minimum point is minimum;
Within the scope of QRS wave group terminal and R peak, obtain minimum point as the second minimum point;
When described second minimum point is one, determine described S peak position according to described second minimum point;
When described second minimum point is multiple, obtain the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum, determine described S peak position according to the minimum point that slope characteristics in the second minimum point is minimum.
Preferably, the step of the relevant information of described extraction T ripple and P ripple comprises:
Calculate the separation of T ripple and P ripple;
Between the starting point and described separation of described QRS wave group, obtain maximum point as the first maximum point, determine to extract the crest location of described T ripple according to the first maximum point that amplitude is maximum by the first maximum point that amplitude is maximum;
Between the terminal and described separation of described QRS wave group, obtain maximum point as the second maximum point, determine to extract the crest location of described P ripple according to the second maximum point that amplitude is maximum by the second maximum point that amplitude is maximum.
In addition, for achieving the above object, the present invention also provides a kind of extraction element of electrocardiosignal characteristic information, comprising:
Acquisition module, for obtaining pending electrocardiosignal;
Processing module, for determining the scope of QRS wave group in described pending electrocardiosignal; Also within the scope of QRS wave group, determine R peak position, and determine the position at Q peak and S peak according to QRS wave group scope and determined R peak position;
Extraction module, for between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.
Preferably, described extraction module, also for the both sides at T wave-wave peak and P wave-wave peak, extracts the starting point and ending point of T ripple, P ripple by default partial transformation method.
Preferably, described processing module comprises:
Resolving cell, for carrying out empirical mode decomposition to described pending electrocardiosignal, obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
Processing unit, for according to described single order high frequency intrinsic mode functions and moving average filter function, determines the characteristic signal of described pending electrocardiosignal, determines the scope of described QRS wave group according to the scope of described characteristic signal.
Preferably, described processing module also comprises:
Acquiring unit, for obtaining the maximum point in described QRS wave group;
Described processing unit, also for when described maximum point is one, determines described R peak position according to described maximum point;
Described acquiring unit, time also for being multiple at described maximum point, obtains the maximum point that in maximum point preset range described in each, slope characteristics is maximum;
Described processing unit, also for determining described R peak position according to obtained maximum point;
Described acquiring unit, also within the scope of QRS wave group starting point and R peak, obtains minimum point as the first minimum point;
Described processing unit, also for when described first minimum point is one, determines described Q peak position according to described first minimum point;
Described acquiring unit, also for when described first minimum point is multiple, obtains the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum;
Described processing unit, also determines described Q peak position for the minimum point minimum according to slope characteristics in the first minimum point;
Described acquiring unit, also within the scope of QRS wave group terminal and R peak, obtains minimum point as the second minimum point;
Described processing unit, also for when described second minimum point is one, determines described S peak position according to described second minimum point;
Described acquiring unit, also for when described second minimum point is multiple, obtains the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum;
Described processing unit, also determines described S peak position for the minimum point minimum according to slope characteristics in the second minimum point.
Preferably, described extraction module comprises:
Computing unit, for calculating the separation of T ripple and P ripple;
Determining unit, for obtaining maximum point as the first maximum point between the starting point and described separation of described QRS wave group, determines the first maximum point that amplitude is maximum;
Extraction unit, extracts the crest location of described T ripple for the first maximum point maximum according to amplitude;
Described determining unit, also for obtaining maximum point as the second maximum point between the terminal and described separation of described QRS wave group, determines the second maximum point that amplitude is maximum;
Described extraction unit, also determines the crest location of described P ripple for the second maximum point maximum according to amplitude.
The present invention by determining the scope of QRS wave group according to electrocardiosignal, and based on QRS wave group scope determine R peak position, Q peak and S peak position and extract the crest location of T ripple and P ripple, realize the extraction of electrocardiosignal characteristic information.Make the strong adaptability that the characteristic information of electrocardiosignal obtains, and decrease the amount of calculation of electrocardiosignal feature information extraction.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the preferred embodiment of the extracting method of electrocardiosignal characteristic information of the present invention;
Fig. 2 is the schematic flow sheet that the present invention determines scope one embodiment of QRS wave group in described pending electrocardiosignal;
Fig. 3 is that the present invention determines R peak position, and determines the schematic flow sheet of position one embodiment at Q peak and S peak according to QRS wave group scope and determined R peak position;
Fig. 4 is the schematic flow sheet that the present invention extracts crest location one embodiment of T ripple and P ripple;
Fig. 5 is the schematic flow sheet that the present invention extracts starting point and ending point one embodiment of T ripple and P ripple;
Fig. 6 is the extraction effect figure of electrocardiosignal characteristic information of the present invention;
Fig. 7 is the schematic flow sheet of extraction one embodiment of electrocardiosignal characteristic information of the present invention;
Fig. 8 is the high-level schematic functional block diagram of the preferred embodiment of the extraction element of electrocardiosignal characteristic information of the present invention;
Fig. 9 is the refinement high-level schematic functional block diagram of processing module one embodiment in Fig. 8;
Figure 10 is the refinement high-level schematic functional block diagram of another embodiment of processing module in Fig. 8;
Figure 11 is the refinement high-level schematic functional block diagram of extraction module one embodiment in Fig. 8.
The realization of the object of the invention, functional characteristics and advantage will in conjunction with the embodiments, are described further with reference to accompanying drawing.
Detailed description of the invention
Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
The primary solutions of the embodiment of the present invention is: obtain pending electrocardiosignal, determines the scope of QRS wave group in described pending electrocardiosignal; Within the scope of QRS wave group, determine R peak position, and determine the position at Q peak and S peak according to QRS wave group scope and determined R peak position; Between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.Extract the characteristic information of electrocardiosignal by the way, make the strong adaptability that the characteristic information of electrocardiosignal obtains, and decrease the amount of calculation of electrocardiosignal feature information extraction.
Due to mode or the bad adaptability of existing electrocardiosignal feature extraction, or it is large to calculate amount, cannot accomplish adaptable while, reduces the problem of the amount of calculation of electrocardiosignal.
Therefore, based on the problems referred to above, the invention provides a kind of method.
With reference to the schematic flow sheet that Fig. 1, Fig. 1 are the preferred embodiment of the extracting method of electrocardiosignal characteristic information of the present invention.
In one embodiment, the extracting method of described electrocardiosignal characteristic information comprises:
Step S10, obtains pending electrocardiosignal, determines the scope of QRS wave group in described pending electrocardiosignal;
In the present embodiment, the obtain manner of described pending electrocardiosignal comprises: obtaining from the electrocardiosignal that this locality is preserved needs electrocardiosignal to be processed, or the electrocardiosignal of Real-time Collection user to be detected.The electrocardiosignal of described Real-time Collection user to be detected can be gathered by electrocardiograph.
After getting pending electrocardiosignal, in order to ensure the accuracy of the electrocardiosignal gathered and reduce data volume:
Carry out pretreatment to described pending electrocardiosignal, described pretreatment includes but not limited to by digital filtering algorithm filtering Hz noise, baseline drift and other noise jamming etc.Utilize the scope of QRS wave group in empirical mode decomposition (EMD, ExperimentalModeDecomposition) algorithm determination electrocardiosignal.QRS wave group: the electric ignition process of reflection Ventricular, the width claiming QRS wave group is the QRS time limit, represents the time required for whole ventricular muscles ignition process.Normal person is the highest is no more than 0.10s.QRS wave group time lengthening, sees ventricular hypertrophy, intraventricular conductional block and preexcitation syndrome.
Concrete, with reference to figure 2, describedly determine that the step of the scope of QRS wave group in described pending electrocardiosignal comprises:
Step S11, carries out empirical mode decomposition to described pending electrocardiosignal, obtains its intrinsic mode functions, and from the intrinsic mode functions obtained, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
Step S12, according to described single order high frequency intrinsic mode functions and moving average filter function, determines the characteristic signal of described pending electrocardiosignal, determines the scope of described QRS wave group according to the scope of described characteristic signal.
In the present embodiment, described empirical mode decomposition is suitable for analyzing non-linear, non-stationary signal sequence, has very high signal to noise ratio.Empirical mode decomposition, it can make sophisticated signal be decomposed into limited intrinsic mode functions (IntrinsicModeFunction, be called for short IMF), each IMF component of decomposing out contain the local feature signal of the Different time scales of original signal.It carries out signal decomposition according to the time scale feature of data self, need not preset any basic function.
Empirical mode decomposition is carried out to electrocardiosignal and obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose single order high frequency intrinsic mode functions I1 (x) of described pending electrocardiosignal, calculate the characteristic signal of described pending electrocardiosignal:
Ftsg (x)=MovAvg (│ I
1(x) │, fs/K), wherein, K is configuration parameter, and for the smoothness of adjustment feature signal, function y=MovAvg (x, N) for moving average filter function, its transfer function is:
setting threshold value TH, the scope meeting ftsg (x) >TH characteristic signal correspond to the scope of QRS wave group.
According to described high frequency single order high frequency intrinsic mode functions and moving average filter function, determine the characteristic signal of described pending electrocardiosignal, determine the scope of described QRS wave group according to the scope of described characteristic signal.
Step S20, within the scope of QRS wave group, determines R peak position, and determines the position at Q peak and S peak according to QRS wave group scope and determined R peak position; Described R peak, Q peak and S peak are the proper noun of medical domain, do not repeat them here!
Concrete, with reference to figure 3, determine R peak position, and determine that the process of the position at Q peak and S peak can comprise according to QRS wave group scope and determined R peak position:
Step S21, obtains the maximum point in described QRS wave group;
Step S22, when described maximum point is one, determines described R peak position according to described maximum point;
Step S23, when described maximum point is multiple, obtains the maximum point that in maximum point preset range described in each, slope characteristics is maximum, determines described R peak position according to obtained maximum point;
Step S24, within the scope of QRS wave group starting point and R peak, obtains minimum point as the first minimum point;
Step S25, when described first minimum point is one, determines described Q peak position according to described first minimum point;
Step S26, when described first minimum point is multiple, obtains the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum, determines described Q peak position according to the minimum point that slope characteristics in the first minimum point is minimum;
Step S27, within the scope of QRS wave group terminal and R peak, obtains minimum point as the second minimum point;
Step S28, when described second minimum point is one, determines described S peak position according to described second minimum point;
Step S29, when described second minimum point is multiple, obtains the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum, determines described S peak position according to the minimum point that slope characteristics in the second minimum point is minimum.
In the present embodiment, step S21 to step S23 is the deterministic process of R peak position; Step S24 to step S26 is the deterministic process of Q peak position; Step S27 to step S29 is the deterministic process of S peak position.The deterministic process at above-mentioned each peak can be carried out simultaneously, or selects wherein one or more preferentially to carry out according to arranging.
The concrete process finding out each peak is: find out the maximum point in QRS wave group, be designated as R
peak(i), i=1,2 ...If R
peaki (), more than 1, compares slope characteristics R near maximum point
ft(i):
Rft (i)=2*E (R
peak(i))-E (R
peak(i)-M)-E (R
peak(i)+M), i=1,2 ..., M is the integer between 2 ~ 6.Slope characteristics R near R peak
fti () should be maximum, can determine R peak position thus.Then, in QRS wave group starting point within the scope of R peak, find minimum point, be designated as Q
peak(i), i=1,2 ..., according to slope calculations feature R
fti method that () is identical, slope calculations feature Q
ft(i), and determine Q peak position thus.Finally, at R peak in QRS wave group ending range, adopting the mode identical with detecting Q peak, determining S peak position.
Step S30, between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.Described T ripple and P ripple are the proper noun of medical domain, do not repeat them here!
Concrete, with reference to figure 4, the process of the crest location of described extraction T ripple and P ripple can comprise:
Step S31, calculates the separation of T ripple and P ripple;
Step S32, obtains maximum point as the first maximum point between the starting point and described separation of described QRS wave group, determines to extract the crest location of described T ripple according to the first maximum point that amplitude is maximum by the first maximum point that amplitude is maximum;
Step S33, obtains maximum point as the second maximum point between the terminal and described separation of described QRS wave group, determines to extract the crest location of described P ripple according to the second maximum point that amplitude is maximum by the second maximum point that amplitude is maximum.
First, the separation TP of T ripple and P ripple is determined
sep(i) position:
TP
sep(i)=3/5 (QRS
off(i)-QRS
on(i))+QRS
on(i), QRS
on(i) and QRS
offi () represents starting point and the final position of QRS wave group respectively.Then, at QRS
on(i)-TP
sepfind maximum point in (i) scope, the crest location of what amplitude was the highest be T ripple.Equally, at TP
sep(i)-QRS
offfind maximum point in (i) scope, the crest location of what amplitude was the highest be P ripple.
Concrete, with reference to figure 5, the process of the starting point and ending point of described extraction T ripple and P ripple can comprise:
Step S34, in the both sides at T wave-wave peak and P wave-wave peak, extracts the starting point and ending point of T ripple, P ripple by default partial transformation method.
In the present embodiment, described default partial transformation method is for improving partial transformation method.To extract the starting point T of T ripple
oni () is example, the described process by presetting partial transformation method extraction T ripple starting point is as follows: first calculate S peak and be connected straight line L with T peak
sTexpression formula.Then, calculating S peak and the peak-to-peak ecg signal data point of T are to L
sTvertical dimension, draw characteristic sequence D (i).Finally, find the minimum point of sequence D (i), corresponding to S peak and the maximum point of T peak-to-peak electrocardiosignal slope variation, be the starting point T of T ripple
onthe position at (i) place.The final position T of T ripple
offthe starting point of (i), P ripple and final position P
on(i), P
offi () is determined all in the same way.With reference to figure 6, be the design sketch of the extracting method one embodiment feature information extraction of above electrocardiosignal characteristic information, R peak marks with " △ ", and Q peak, S peak, T peak and P peak mark with " * ", T
on(i) and P
oni () marks with " ", T
off(i) and P
offi () marks with " ".
In order to better describe the process of the extraction of electrocardiosignal characteristic information of the present invention, with reference to figure 7, the flow process that described electrocardiosignal is extracted comprises: 1, ECG signal processing; 2, QRS wave group scope is determined; 3, R peak, Q peak, S peak position are determined; 4, T ripple and P wave-wave peak position are determined; 5, T ripple and P ripple starting point, terminating point are determined.The detailed process of 1 to 5 is as described in each embodiment above-mentioned.
The present embodiment by determining the scope of QRS wave group according to electrocardiosignal, and based on QRS wave group scope determine R peak position, Q peak and S peak position and extract the crest location of T ripple and P ripple, realize the extraction of electrocardiosignal characteristic information.Make the strong adaptability that the characteristic information of electrocardiosignal obtains, and decrease the amount of calculation of electrocardiosignal feature information extraction.
The executive agent of the extracting method of the electrocardiosignal characteristic information of above-described embodiment can be all ECG Signal Analysis equipment.Further, the extracting method of this electrocardiosignal characteristic information can be realized by the extraction procedure of the client electrocardiosignal characteristic information be arranged on ECG Signal Analysis equipment, wherein, this ECG Signal Analysis equipment can include but not limited to ECG signal sampling instrument, computer, PC etc.
The extracting method of corresponding above-mentioned electrocardiosignal characteristic information, the present invention further provides a kind of extraction element of electrocardiosignal characteristic information.
With reference to the high-level schematic functional block diagram that Fig. 8, Fig. 8 are the preferred embodiment of the extraction element of electrocardiosignal characteristic information of the present invention.
In one embodiment, the extraction element of described electrocardiosignal characteristic information comprises: acquisition module 10, processing module 20 and extraction module 30.
Described acquisition module 10, for obtaining pending electrocardiosignal;
Described processing module 20, for determining the scope of QRS wave group in described pending electrocardiosignal;
In the present embodiment, the obtain manner of described pending electrocardiosignal comprises: obtaining from the electrocardiosignal that this locality is preserved needs electrocardiosignal to be processed, or the electrocardiosignal of Real-time Collection user to be detected.The electrocardiosignal of described Real-time Collection user to be detected can be gathered by electrocardiograph.
After getting pending electrocardiosignal, in order to ensure the accuracy of the electrocardiosignal gathered and reduce data volume:
Carry out pretreatment to described pending electrocardiosignal, described pretreatment includes but not limited to by digital filtering algorithm filtering Hz noise, baseline drift and other noise jamming etc.Utilize the scope of QRS wave group in empirical mode decomposition (EMD, ExperimentalModeDecomposition) algorithm determination electrocardiosignal.QRS wave group: the electric ignition process of reflection Ventricular, the width claiming QRS wave group is the QRS time limit, represents the time required for whole ventricular muscles ignition process.Normal person is the highest is no more than 0.10s.QRS wave group time lengthening, sees ventricular hypertrophy, intraventricular conductional block and preexcitation syndrome.
Concrete, with reference to figure 9, described processing module 20 comprises resolving cell 21 and processing unit 22,
Described resolving cell 21, for carrying out empirical mode decomposition to described pending electrocardiosignal, obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
Described processing unit 22, for according to described single order high frequency intrinsic mode functions and moving average filter function, determines the characteristic signal of described pending electrocardiosignal, determines the scope of described QRS wave group according to the scope of described characteristic signal.
In the present embodiment, described empirical mode decomposition is suitable for analyzing non-linear, non-stationary signal sequence, has very high signal to noise ratio.Empirical mode decomposition, it can make sophisticated signal be decomposed into limited intrinsic mode functions (IntrinsicModeFunction, be called for short IMF), each IMF component of decomposing out contain the local feature signal of the Different time scales of original signal.It carries out signal decomposition according to the time scale feature of data self, need not preset any basic function.
Empirical mode decomposition is carried out to electrocardiosignal and obtains its intrinsic mode functions (IMF, IntrinsicModeFunction), and from obtained intrinsic mode functions, choose single order high frequency intrinsic mode functions I1 (x) of described pending electrocardiosignal, calculate the characteristic signal of described pending electrocardiosignal:
Ftsg (x)=MovAvg (│ I
1(x) │, fs/K), wherein, K is configuration parameter, and for the smoothness of adjustment feature signal, function y=MovAvg (x, N) for moving average filter function, its transfer function is:
setting threshold value TH, the scope meeting ftsg (x) >TH characteristic signal correspond to the scope of QRS wave group.
According to described single order high frequency intrinsic mode functions and moving average filter function, determine the characteristic signal of described pending electrocardiosignal, determine the scope of described QRS wave group according to the scope of described characteristic signal.
Described processing module 20, also within the scope of QRS wave group, determines R peak position, and determines the position at Q peak and S peak according to QRS wave group scope and determined R peak position;
Concrete, with reference to Figure 10, described processing module 20 also comprises acquiring unit 23,
Described acquiring unit 23, for obtaining the maximum point in described QRS wave group;
Described processing unit 22, also for when described maximum point is one, determines described R peak position according to described maximum point;
Described acquiring unit 23, time also for being multiple at described maximum point, obtains the maximum point that in maximum point preset range described in each, slope characteristics is maximum;
Described processing unit 22, also for determining described R peak position according to obtained maximum point;
Described acquiring unit 23, also within the scope of QRS wave group starting point and R peak, obtains minimum point as the first minimum point;
Described processing unit 22, also for when described first minimum point is one, determines described Q peak position according to described minimum point;
Described acquiring unit 23, also for when described first minimum point is multiple, obtains the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum;
Described processing unit 22, also determines described Q peak position for the minimum point minimum according to slope characteristics in the first minimum point;
Described acquiring unit 23, also within the scope of QRS wave group terminal and R peak, obtains minimum point as the second minimum point;
Described processing unit 22, also for when described second minimum point is one, determines described S peak position according to described minimum point;
Described acquiring unit 23, also for when described second minimum point is multiple, obtains the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum;
Described processing unit 22, also determines described S peak position for the minimum point minimum according to slope characteristics in the second minimum point.
In the present embodiment, the concrete process finding out each peak is: find out the maximum point in QRS wave group, be designated as R
peak(i), i=1,2 ...If R
peaki (), more than 1, compares slope characteristics R near maximum point
ft(i):
R
ft(i)=2 × E (R
peak(i)-E (R
peak(i)-M)-E (R
peak(i)+M), i=1,2 ..., M is the integer between 2 ~ 6.Slope characteristics R near R peak
fti () should be maximum, can determine R peak position thus.Then, in QRS wave group starting point within the scope of R peak, find minimum point, be designated as Q
peak(i), i=1,2 ..., according to slope calculations feature R
fti method that () is identical, slope calculations feature Q
ft(i), and determine Q peak position thus.Finally, at R peak in QRS wave group ending range, adopting the mode identical with detecting Q peak, determining S peak position.
Described extraction module 30, for between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.
With reference to Figure 11, described extraction module 30 comprises: computing unit 31, determining unit 32 and extraction unit 33,
Described computing unit 31, for calculating the separation of T ripple and P ripple;
Described determining unit 32, for obtaining maximum point as the first maximum point between the starting point and described separation of described QRS wave group, determines the first maximum point that amplitude is maximum;
Described extraction unit 33, extracts the crest location of described T ripple for the first maximum point maximum according to amplitude;
Described determining unit 32, also for obtaining maximum point as the second maximum point between the terminal and described separation of described QRS wave group, determines the second maximum point that amplitude is maximum;
Described extraction unit 33, also extracts the crest location of described P ripple for the second maximum point maximum according to amplitude.
First, the separation TP of T ripple and P ripple is determined
sep(i) position:
TP
sep(i)=3/5 (QRS
off(i)-QRS
on(i))+QRS
on(i), QRS
on(i), QRS
off(i) and QRS
offi () represents starting point and the final position of QRS wave group respectively.Then, at QRS
on(i)-TP
sepfind maximum point in (i) scope, the crest location of what amplitude was the highest be T ripple.Equally, at TP
sep(i)-QRS
offfind maximum point in (i) scope, the crest location of what amplitude was the highest be P ripple.
The present embodiment by determining the scope of QRS wave group according to electrocardiosignal, and based on QRS wave group scope determine R peak position, Q peak and S peak position and extract the crest location of T ripple and P ripple, realize the extraction of electrocardiosignal characteristic information.Make the strong adaptability that the characteristic information of electrocardiosignal obtains, and decrease the amount of calculation of electrocardiosignal feature information extraction.
Further, described extraction unit 33, also for the both sides at T wave-wave peak and P wave-wave peak, extracts the starting point and ending point of T ripple, P ripple by default partial transformation method.
In the present embodiment, described default partial transformation method is for improving partial transformation method.To extract the starting point T of T ripple
oni () is example, the described process by presetting partial transformation method extraction T ripple starting point is as follows: first calculate S peak and be connected straight line L with T peak
sTexpression formula.Then, calculating S peak and the peak-to-peak ecg signal data point of T are to L
sTvertical dimension, draw characteristic sequence D (i).Finally, find the minimum point of sequence D (i), corresponding to S peak and the maximum point of T peak-to-peak electrocardiosignal slope variation, be the starting point T of T ripple
onthe position at (i) place.The final position T of T ripple
offthe starting point of (i), P ripple and final position P
on(i), P
offi () is determined all in the same way.Reference
fig. 6, be the effect of the extracting method one embodiment feature information extraction of above electrocardiosignal characteristic information
figure, R peak marks with " △ ", and Q peak, S peak, T peak and P peak mark with " * ", T
on(i) and P
oni () is used
mark, T
off(i) and P
offi () is used
mark.
In order to better describe the process of the extraction of electrocardiosignal characteristic information of the present invention, with reference to figure 7, the flow process that described electrocardiosignal is extracted comprises: 1, ECG signal processing; 2, QRS wave group scope is determined; 3, R peak, Q peak, S peak position are determined; 4, T ripple and P wave-wave peak position are determined; 5, T ripple and P ripple starting point, terminating point are determined.The detailed process of 1 to 5 is as described in each embodiment above-mentioned.
The present embodiment is by extracting the starting point and ending point of T ripple and P ripple according to default localized variation method, the characteristic information of the electrocardiosignal improved further, ensureing the adaptable colleague of electrocardiosignal feature extraction, reducing amount of calculation, the analysis efficiency of electrocardiosignal is improved.
These are only the preferred embodiments of the present invention; not thereby the scope of the claims of the present invention is limited; every utilize description of the present invention and accompanying drawing content to do equivalent structure or equivalent flow process conversion; or be directly or indirectly used in other relevant technical fields, be all in like manner included in scope of patent protection of the present invention.
Claims (10)
1. an extracting method for electrocardiosignal characteristic information, is characterized in that, comprising:
Obtain pending electrocardiosignal, determine the scope of QRS wave group in described pending electrocardiosignal;
Within the scope of QRS wave group, determine R peak position, and determine the position at Q peak and S peak according to QRS wave group scope and determined R peak position;
Between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.
2. the extracting method of electrocardiosignal characteristic information as claimed in claim 1, is characterized in that, the step of the relevant information of described extraction T ripple and P ripple comprises:
In the both sides at T wave-wave peak and P wave-wave peak, extracted the starting point and ending point of T ripple, P ripple by default partial transformation method.
3. the extracting method of electrocardiosignal characteristic information as claimed in claim 1 or 2, is characterized in that, describedly determines that the step of the scope of QRS wave group in described pending electrocardiosignal comprises:
Empirical mode decomposition is carried out to described pending electrocardiosignal, obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
According to described single order high frequency intrinsic mode functions and moving average filter function, determine the characteristic signal of described pending electrocardiosignal, determine the scope of described QRS wave group according to the scope of described characteristic signal.
4. the extracting method of electrocardiosignal characteristic information as claimed in claim 1 or 2, is characterized in that, describedly determines R peak position, and determines that the step of the position at Q peak and S peak comprises according to QRS wave group scope and determined R peak position:
Obtain the maximum point in described QRS wave group;
When described maximum point is one, determine described R peak position according to described maximum point;
When described maximum point is multiple, obtain the maximum point that in maximum point preset range described in each, slope characteristics is maximum, determine described R peak position according to obtained maximum point;
Within the scope of QRS wave group starting point and R peak, obtain minimum point as the first minimum point;
When described first minimum point is one, determine described Q peak position according to described first minimum point;
When described first minimum point is multiple, obtain the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum, determine described Q peak position according to the minimum point that slope characteristics in the first minimum point is minimum;
Within the scope of QRS wave group terminal and R peak, obtain minimum point as the second minimum point;
When described second minimum point is one, determine described S peak position according to described second minimum point;
When described second minimum point is multiple, obtain the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum, determine described S peak position according to the minimum point that slope characteristics in the second minimum point is minimum.
5. the extracting method of electrocardiosignal characteristic information as claimed in claim 1 or 2, is characterized in that, the step of the relevant information of described extraction T ripple and P ripple comprises:
Calculate the separation of T ripple and P ripple;
Between the starting point and described separation of described QRS wave group, obtain maximum point as the first maximum point, determine to extract the crest location of described T ripple according to the first maximum point that amplitude is maximum by the first maximum point that amplitude is maximum;
Between the terminal and described separation of described QRS wave group, obtain maximum point as the second maximum point, determine to extract the crest location of described P ripple according to the second maximum point that amplitude is maximum by the second maximum point that amplitude is maximum.
6. an extraction element for electrocardiosignal characteristic information, is characterized in that, comprising:
Acquisition module, for obtaining pending electrocardiosignal;
Processing module, for determining the scope of QRS wave group in described pending electrocardiosignal; Also within the scope of QRS wave group, determine R peak position, and determine the position at Q peak and S peak according to QRS wave group scope and determined R peak position;
Extraction module, for between described QRS wave group and adjacent QRS wave group, extract the relevant information of T ripple and P ripple, the scope of determined QRS wave group, R peak position, Q peak and the position at S peak and the relevant information of the T ripple extracted and P ripple are the characteristic information of described pending electrocardiosignal, and described relevant information comprises the crest location of T ripple and P ripple and the starting point and ending point of T ripple and P ripple.
7. the extraction element of electrocardiosignal characteristic information as claimed in claim 6, is characterized in that, described extraction module, also for the both sides at T wave-wave peak and P wave-wave peak, is extracted the starting point and ending point of T ripple, P ripple by default partial transformation method.
8. the extraction element of electrocardiosignal characteristic information as claimed in claims 6 or 7, it is characterized in that, described processing module comprises:
Resolving cell, for carrying out empirical mode decomposition to described pending electrocardiosignal, obtains its intrinsic mode functions, and from obtained intrinsic mode functions, choose the single order high frequency intrinsic mode functions of described pending electrocardiosignal;
Processing unit, for according to described single order high frequency intrinsic mode functions and moving average filter function, determines the characteristic signal of described pending electrocardiosignal, determines the scope of described QRS wave group according to the scope of described characteristic signal.
9. the extraction element of electrocardiosignal characteristic information as claimed in claim 8, it is characterized in that, described processing module also comprises:
Acquiring unit, for obtaining the maximum point in described QRS wave group;
Described processing unit, also for when described maximum point is one, determines described R peak position according to described maximum point;
Described acquiring unit, time also for being multiple at described maximum point, obtains the maximum point that in maximum point preset range described in each, slope characteristics is maximum;
Described processing unit, also for determining described R peak position according to obtained maximum point;
Described acquiring unit, also within the scope of QRS wave group starting point and R peak, obtains minimum point as the first minimum point;
Described processing unit, also for when described first minimum point is one, determines described Q peak position according to described first minimum point;
Described acquiring unit, also for when described first minimum point is multiple, obtains the minimum point that in the first minimum point preset range described in each, slope characteristics is minimum;
Described processing unit, also determines described Q peak position for the minimum point minimum according to slope characteristics in the first minimum point;
Described acquiring unit, also within the scope of QRS wave group terminal and R peak, obtains minimum point as the second minimum point;
Described processing unit, also for when described second minimum point is one, determines described S peak position according to described second minimum point;
Described acquiring unit, also for when described second minimum point is multiple, obtains the minimum point that in the second minimum point preset range described in each, slope characteristics is minimum;
Described processing unit, also determines described S peak position for the minimum point minimum according to slope characteristics in the second minimum point.
10. the extraction element of electrocardiosignal characteristic information as claimed in claims 6 or 7, it is characterized in that, described extraction module comprises:
Computing unit, for calculating the separation of T ripple and P ripple;
Determining unit, for obtaining maximum point as the first maximum point between the starting point and described separation of described QRS wave group, determines the first maximum point that amplitude is maximum;
Extraction unit, extracts the crest location of described T ripple for the first maximum point maximum according to amplitude;
Described determining unit, also for obtaining maximum point as the second maximum point between the terminal and described separation of described QRS wave group, determines the second maximum point that amplitude is maximum;
Described extraction unit, also determines the crest location of described P ripple for the second maximum point maximum according to amplitude.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510172501.2A CN105054926A (en) | 2015-04-13 | 2015-04-13 | Electrocardiosignal feature information extraction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510172501.2A CN105054926A (en) | 2015-04-13 | 2015-04-13 | Electrocardiosignal feature information extraction method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105054926A true CN105054926A (en) | 2015-11-18 |
Family
ID=54484630
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510172501.2A Pending CN105054926A (en) | 2015-04-13 | 2015-04-13 | Electrocardiosignal feature information extraction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105054926A (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105877739A (en) * | 2016-02-25 | 2016-08-24 | 姜坤 | Clinical examination method of electrocardio intelligent analyzing system |
| CN105997055A (en) * | 2016-07-11 | 2016-10-12 | 吉林大学 | Automatic classification method, system and device of electrocardiosignal ST band |
| CN106344004A (en) * | 2016-09-28 | 2017-01-25 | 清华大学 | Electrocardiosignal feature point detecting method and device |
| CN106725450A (en) * | 2017-02-20 | 2017-05-31 | 武汉中旗生物医疗电子有限公司 | QRS complex detection method |
| CN106974644A (en) * | 2017-04-25 | 2017-07-25 | 深圳开立生物医疗科技股份有限公司 | A kind of electrocardiogram R wave detection method and system |
| CN107212881A (en) * | 2017-05-26 | 2017-09-29 | 广东工业大学 | A kind of T ripples electrical alternations detection method |
| CN107997759A (en) * | 2017-10-27 | 2018-05-08 | 北京康博众联电子科技有限公司 | The detection method of QRS complex, storage medium and computer equipment in electrocardiosignal |
| CN109171711A (en) * | 2018-11-14 | 2019-01-11 | 中国计量大学 | A kind of quick P wave detecting method based on extremum method |
| CN109875550A (en) * | 2019-04-02 | 2019-06-14 | 东北大学 | A kind of sequences of ventricular depolarization critical point detection method |
| CN110897627A (en) * | 2018-09-14 | 2020-03-24 | 杭州脉流科技有限公司 | Electrocardiogram signal feature extraction method, device, equipment, system and storage medium |
| CN111248895A (en) * | 2018-11-30 | 2020-06-09 | 无锡祥生医疗科技股份有限公司 | Electrocardiosignal characteristic detection method and system |
| CN111543977A (en) * | 2020-05-09 | 2020-08-18 | 益体康(北京)科技有限公司 | Multi-cascade artificial intelligence vagina discharge method based on 12-lead resting electrocardiogram |
| CN113116358A (en) * | 2019-12-30 | 2021-07-16 | 华为技术有限公司 | Display method and device of electrocardiogram, terminal equipment and storage medium |
| CN113786202A (en) * | 2021-08-17 | 2021-12-14 | 广东省科学院健康医学研究所 | Electrocardiogram characteristic starting point and end point detection method, system, device and medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4896677A (en) * | 1987-12-26 | 1990-01-30 | Fukuda Denshi Kabushiki Kaisha | Electrocardiographic waveform display apparatus, and method of expressing electrocardiographic waveforms |
| WO2001091627A2 (en) * | 2000-05-30 | 2001-12-06 | Vladimir Shusterman | System and device for multi-scale analysis and representation of electrocardiographic data |
| CN101088456A (en) * | 2007-07-13 | 2007-12-19 | 天津大学 | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution |
| CN101766484A (en) * | 2010-01-18 | 2010-07-07 | 董军 | Method and equipment for identification and classification of electrocardiogram |
| CN102085095A (en) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | Method, system and electrocardioscanner for detecting ST segment in electrocardiogram |
| CN102626310A (en) * | 2012-04-23 | 2012-08-08 | 天津工业大学 | Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving |
| CN104382590A (en) * | 2014-12-11 | 2015-03-04 | 赖大坤 | Automatic shockable rhythm identification and classification method combined with electrocardio time-frequency domain feature analysis |
-
2015
- 2015-04-13 CN CN201510172501.2A patent/CN105054926A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4896677A (en) * | 1987-12-26 | 1990-01-30 | Fukuda Denshi Kabushiki Kaisha | Electrocardiographic waveform display apparatus, and method of expressing electrocardiographic waveforms |
| WO2001091627A2 (en) * | 2000-05-30 | 2001-12-06 | Vladimir Shusterman | System and device for multi-scale analysis and representation of electrocardiographic data |
| CN101088456A (en) * | 2007-07-13 | 2007-12-19 | 天津大学 | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution |
| CN102085095A (en) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | Method, system and electrocardioscanner for detecting ST segment in electrocardiogram |
| CN101766484A (en) * | 2010-01-18 | 2010-07-07 | 董军 | Method and equipment for identification and classification of electrocardiogram |
| CN102626310A (en) * | 2012-04-23 | 2012-08-08 | 天津工业大学 | Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving |
| CN104382590A (en) * | 2014-12-11 | 2015-03-04 | 赖大坤 | Automatic shockable rhythm identification and classification method combined with electrocardio time-frequency domain feature analysis |
Non-Patent Citations (3)
| Title |
|---|
| 李雪飞: "心电自动分析系统的研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
| 覃爱娜等: "一种改进的QRS波检测方法", 《中南大学学报》 * |
| 黄家洺: "基于物联网的移动医疗监护系统的设计与实现", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105877739A (en) * | 2016-02-25 | 2016-08-24 | 姜坤 | Clinical examination method of electrocardio intelligent analyzing system |
| CN105997055A (en) * | 2016-07-11 | 2016-10-12 | 吉林大学 | Automatic classification method, system and device of electrocardiosignal ST band |
| CN106344004A (en) * | 2016-09-28 | 2017-01-25 | 清华大学 | Electrocardiosignal feature point detecting method and device |
| CN106725450B (en) * | 2017-02-20 | 2020-03-27 | 武汉中旗生物医疗电子有限公司 | QRS complex detection method |
| CN106725450A (en) * | 2017-02-20 | 2017-05-31 | 武汉中旗生物医疗电子有限公司 | QRS complex detection method |
| CN106974644A (en) * | 2017-04-25 | 2017-07-25 | 深圳开立生物医疗科技股份有限公司 | A kind of electrocardiogram R wave detection method and system |
| CN107212881B (en) * | 2017-05-26 | 2020-08-11 | 广东工业大学 | T-wave electricity alternative detection method |
| CN107212881A (en) * | 2017-05-26 | 2017-09-29 | 广东工业大学 | A kind of T ripples electrical alternations detection method |
| CN107997759A (en) * | 2017-10-27 | 2018-05-08 | 北京康博众联电子科技有限公司 | The detection method of QRS complex, storage medium and computer equipment in electrocardiosignal |
| CN110897627A (en) * | 2018-09-14 | 2020-03-24 | 杭州脉流科技有限公司 | Electrocardiogram signal feature extraction method, device, equipment, system and storage medium |
| CN109171711A (en) * | 2018-11-14 | 2019-01-11 | 中国计量大学 | A kind of quick P wave detecting method based on extremum method |
| CN111248895A (en) * | 2018-11-30 | 2020-06-09 | 无锡祥生医疗科技股份有限公司 | Electrocardiosignal characteristic detection method and system |
| CN109875550A (en) * | 2019-04-02 | 2019-06-14 | 东北大学 | A kind of sequences of ventricular depolarization critical point detection method |
| CN113116358A (en) * | 2019-12-30 | 2021-07-16 | 华为技术有限公司 | Display method and device of electrocardiogram, terminal equipment and storage medium |
| CN111543977A (en) * | 2020-05-09 | 2020-08-18 | 益体康(北京)科技有限公司 | Multi-cascade artificial intelligence vagina discharge method based on 12-lead resting electrocardiogram |
| CN111543977B (en) * | 2020-05-09 | 2023-04-07 | 益体康(北京)科技有限公司 | Multi-cascade artificial intelligence vagina discharge method based on 12-lead resting electrocardiogram |
| CN113786202A (en) * | 2021-08-17 | 2021-12-14 | 广东省科学院健康医学研究所 | Electrocardiogram characteristic starting point and end point detection method, system, device and medium |
| CN113786202B (en) * | 2021-08-17 | 2023-12-26 | 广东省科学院健康医学研究所 | Electrocardiogram characteristic starting point and ending point detection method, system, device and medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105054926A (en) | Electrocardiosignal feature information extraction method and device | |
| CN101449973B (en) | Judgment index generation method and device for cardiac interference signal identification | |
| CN107041743B (en) | Real-time R wave detection method for electrocardiosignals | |
| Zhu et al. | An R-peak detection method based on peaks of Shannon energy envelope | |
| Karvounis et al. | Fetal heart rate extraction from composite maternal ECG using complex continuous wavelet transform | |
| Zou et al. | An ultra-low power QRS complex detection algorithm based on down-sampling wavelet transform | |
| CN103584854B (en) | Extraction method of electrocardiosignal R waves | |
| Zhang et al. | QRS detection based on morphological filter and energy envelope for applications in body sensor networks | |
| CN102626310A (en) | Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving | |
| CN108888259A (en) | A kind of real-time QRS wave detection method of electrocardiosignal | |
| Manikandan et al. | Straightforward and robust QRS detection algorithm for wearable cardiac monitor | |
| CN102258368B (en) | Time-domain sparsity linear aliasing blind separation model discrimination method in fetal electrocardiogram detection | |
| CN104305992B (en) | A kind of interactive fast automatic extracting method of Fetal ECG | |
| CN104783787A (en) | J-wave detecting method based on neural network | |
| CN105286860A (en) | Motor imagery brain electrical signal recognition method based on dual-tree complex wavelet energy difference | |
| Boichat et al. | Wavelet-based ECG delineation on a wearable embedded sensor platform | |
| CN103829944A (en) | TTI signal processing method based on pattern recognition | |
| CN107320096B (en) | Electrocardio R wave positioning method | |
| CN108836305B (en) | A kind of ECG feature extracting method of fusion Butterworth filtering and wavelet transformation | |
| Manikandan et al. | Robust detection of premature ventricular contractions using sparse signal decomposition and temporal features | |
| Liu et al. | A novel R-peak detection method combining energy and wavelet transform in electrocardiogram signal | |
| Narsimha et al. | Denoising and QRS detection of ECG signals using Empirical Mode Decomposition | |
| CN208511016U (en) | A kind of vital sign comprehensive detection analysis system | |
| KR101048763B1 (en) | Apparauts and method for detecting signal | |
| Dessì et al. | Identification of fetal QRS complexes in low density non-invasive biopotential recordings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151118 |
|
| RJ01 | Rejection of invention patent application after publication |