CN100536765C - Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution - Google Patents
Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution Download PDFInfo
- Publication number
- CN100536765C CN100536765C CNB2007100580219A CN200710058021A CN100536765C CN 100536765 C CN100536765 C CN 100536765C CN B2007100580219 A CNB2007100580219 A CN B2007100580219A CN 200710058021 A CN200710058021 A CN 200710058021A CN 100536765 C CN100536765 C CN 100536765C
- Authority
- CN
- China
- Prior art keywords
- ripple
- point
- wave
- qrs
- peak
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The cardioelectric data characteristic extracting process includes the following steps: 1.preprocessing cardioelectric data; 2. extracting QRS wave group characteristic through wavelet transformation to extract the sub-frequency band of the QRS wave group from the preprocessed cardioelectric signal and the subsequent evolutive wavelet Wiener deconvolution process to extract the position of the characteristic point of QRS wave group; and 3. extracting the characteristics of P wave and T wave through substituting the time section in the QRS complex wave with base line, the subsequent wavelet transformation to extract sub-frequency band of P wave and T wave from the cardioelectric signal with the QRS wave group eliminated, and the final evolutive wavelet Wiener deconvolution process until extracting accurate characteristic point position of P wave and T wave. The present invention lays foundation for the characteristic detection.
Description
Technical field
The present invention relates to a kind of feature extracting method of electrocardiogram (ECG) data, belong to the biomedical engineering technology field.
Background technology
The feature extraction of electrocardiosignal is one of powerful measure of detecting of heart disease.The algorithm that is generally used for the electrocardiosignal feature extraction has following several: length conversion and energy change of variable, hidden Markov model, artificial neural network and wavelet transformation.But said method all has certain defective: length conversion and energy change of variable can't accurately detect the characteristic information of abnormal Q RS wave group; Hidden Markov model is powerless to the detection of accidental abnormal electrocardiogram ordinary wave shape and isolated P ripple; Artificial Neural Network lacks effectiveness to the detection of abnormal electrocardiogram waveform equally; Small wave converting method is better than said method to the detection of ill electrocardiosignal, but the detection that the ecg wave form unstable state is changed, particularly under the low signal-to-noise ratio situation, and still can not be satisfactory.
In recent years, in order to solve the defective of said method, many scholars have carried out going deep into extensive studies.People such as Messadeg are used for the tagsort of electrocardiosignal by hidden Markov model and wavelet transformation technique are combined; Yu and he work together and adopt independent component analysis and the bonded method of artificial neural network to realize the extraction and the classification of ecg characteristics; Utilized Liu adaptive wavelets transform and fuzzy logic principle realize location, electrocardiosignal R peak, or the like.As can be seen, above-mentioned scholar no longer is simple a kind of theory or the method for adopting, but by remedying mutually between the several different methods, realizes optimizing coupling.But then, though these new trials have solved the subproblem that original ecg characteristics extracts, since perfect not enough, much need improved place so exist.
Summary of the invention
Purport of the present invention is at the on-line monitoring simultaneously of multidigit needs of patients in the remote electrocardiogram monitor, the problem of carrying out the heart disease detection promptly and accurately, a kind of electrocardiogram (ECG) data feature extracting method of efficiently and accurately is proposed, to overcome the defective of prior art medium-long range cardiac monitoring cardiac disease detection.The present invention can improve the efficient and the accuracy of heart disease identification greatly, thereby reduces doctor's operating pressure, makes therefore that also more oversensitive dirty disease patient can accept to guard accurately and effectively and treat.
For this reason, the present invention adopts following technical scheme: a kind of cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution comprises the following steps:
(1) according to the sample frequency of electrocardiosignal, selects the morphological structure operator, utilize the baseline drift in the Mathematical Morphology Method removal original electrocardiographicdigital signal;
(2) initialization Wiener filter, scale parameter
, energy parameter
And iteration step length (β
p), get base and float electrocardiosignal after the filtering, carry out the wavelet transformation under the 1st yardstick, obtain the useful signal of low frequency part, and the noise signal of HFS, obtain the power spectrum of useful signal and noise signal and crosspower spectrum between the two, construct the non-causal Wiener filter, as the initialization Wei Na operator vector that deconvolutes with this
(3) QRS complex wave feature strengthens: adopt the quadrature discrete small echo to decompose pretreated electrocardiosignal, according to the frequency range of electrocardiosignal and QRS complex wave, carry out twice or twice above wavelet transformation to obtain the feature sub-band at QRS wave group place; With the Wei Na operator vector that deconvolutes
Act on the sub-band at QRS wave group place, then calculate the gross energy error
If error amount less than setting threshold, stops computing; Otherwise, utilize steepest descent gradient based method, upgrade the Wei Na operator vector that deconvolutes
Energy parameter
And scale parameter
Continue computing up to finding the satisfied Wei Na operator vector that deconvolutes
(4) P ripple and T wave characteristic strengthen: utilize the positional information of the QRS complex wave that step (3) obtains, the QRS wave group is substituted with baseline; Adopt the quadrature discrete small echo to decompose the electrocardiosignal that removes the QRS wave group, according to gained electrocardiosignal and P ripple, T wave frequency scope, carry out twice or twice above wavelet transformation to obtain the feature sub-band at P ripple and T ripple place; With the Wei Na operator vector that deconvolutes
Act on the feature sub-band at P ripple and T ripple place, then calculate the gross energy error
If error amount less than setting threshold, stops computing; Otherwise, utilize steepest descent gradient based method, upgrade the Wei Na operator vector that deconvolutes
Energy parameter
And scale parameter
Continue computing up to finding the satisfied Wei Na operator vector that deconvolutes
(5) detect QRS wave group feature locations point: get the QRS sub-band signal after feature that step (3) obtains strengthens, at first detect the R peak position,, find out first valley point and peak dot before and after the R peak then according to the R peak position; Then according to following criterion location feature point: if there is valley point of a peak dot at the R peak, peak dot is a QRS wave group starting point so, and the valley point is a Q ripple flex point; If have only a valley point before the R peak, this wave group is the RS wave group so, and this valley point is a RS wave group starting point; Valley point behind the R peak is a S ripple flex point, and peak dot is the wave group terminating point;
(6) detect P, T wave characteristic location point: get P ripple, T marble band signal after feature that step (4) obtains strengthens, detect the peak point position of P ripple and T ripple, then find out first preceding valley point of P ripple and T crest value point respectively, be the starting point of P ripple and T ripple; Find out first valley point behind P ripple and the T crest value point, be the terminating point of P ripple and T ripple.
In above-mentioned cardioelectric characteristic extracting process, the formula that steepest descent gradient based method adopts can be:
The essential characteristics of technical scheme of the present invention is: 1. electrocardiogram (ECG) data pretreatment: select the baseline drift in the appropriate mathematical morphological operator removal electrocardiosignal; 2. the feature extraction of QRS wave group: extract through the sub-band under the QRS wave group in the electrocardiosignal after the pretreatment by small wave converting method,, extract the characteristic point position of QRS wave group again to this sub-band signal utilization evolution type Wei Na method of deconvoluting.3. the feature extraction of P ripple and T ripple: the characteristic point position information of the QRS wave group that obtains according to the first step replaces time period at QRS complex wave place, to avoid the frequency domain interference when P ripple and T wave characteristic are extracted with baseline.Then adopt wavelet transformation, the sub-band of the electrocardiosignal behind the removal QRS wave group being made P ripple, T ripple extracts, and the method for at last the sub-band utilization evolution type Wei Na that obtains being deconvoluted is up to extracting P ripple, T wave characteristic point position accurately.
The present invention proposes based on the evolutive wavelet wiener deconvolution new technique and realize the new approaches that ecg characteristics extracts: utilize the localized advantage of wavelet transformation time-frequency, signature waveform is divided into each sub-frequency bands, then utilize the evolution type Wei Na new technique that deconvolutes that the feature sub-band signal is carried out iterative processing, up to finding satisfied characteristic point position information.Adopt this method, can so just can provide strong basis outstanding tangible peak dot or the valley point of changing into of inapparent key position point in the original electrocardiographicdigital data for feature detection.
Description of drawings
The characteristic information of Fig. 1 ecg wave form and definition.
Fig. 2 is based on the electrocardio baseline drift filtering sketch map of Mathematical Morphology Method, and Fig. 2 (a) is the original electrocardiographicdigital signal; Fig. 2 (b) is the electrocardiosignal after the filtering baseline drift.
The feature of Fig. 3 QRS wave group strengthens result schematic diagram, and Fig. 3 (a) is the electrocardiosignal before QRS wave group feature strengthens; Fig. 3 (b) is the electrocardiosignal after QRS wave group feature strengthens.
The feature of Fig. 4 P ripple and T ripple strengthens result schematic diagram, and Fig. 4 (a) is the electrocardiosignal of QRS wave group after shifting out; Fig. 4 (b) is the electrocardiosignal after P ripple and T wave characteristic strengthen.
The specific embodiment
Before the present invention was further described, existing knowledge such as the character of first principal character to ecg wave form, Mathematical Morphology theory, electrocardiosignal wavelet transform, evolutive wavelet wiener deconvolution technology were done introduction.
The principal character of ecg wave form
Electrocardiosignal is one of topmost physiological signal, and it has reflected the activity of heart, and mainly contain 3 parts and form: the P ripple is represented first deflection of atrium depolarization; The QRS complex wave is produced by sequences of ventricular depolarization; Ventricular bipolar has produced the T ripple.Because these ripples have special shape in time domain and frequency domain, so, by observation, can find the abnormal conditions of heart to the ECG signal.Fig. 1 has provided Electrocardiographic fundamental component, and the definition of every ecg characteristics.Purpose of the present invention will be found out a kind of characteristic detection method accurately and effectively exactly, extracts the characteristic information shown in publishing picture.
Mathematical Morphology theory
Mathematical morphology filter is a kind of nonlinear signal conversion, is used to change the local geometric characteristic of signal.Its theoretical basis comes from the basic operation method of set theory in the graphical analysis, i.e. Mathematical Morphology Method, this theory are proposed by Matheron and Serra the earliest.In this theory, all signals all are seen as the set that has a dimensioning in the Euclidean geometry space, and morphological method just can be counted as a kind of set operation that comes picked up signal geometry quantitative information by the signal conversion.For the binary signal that is counted as gathering, it being corroded and expand with certain structural element is the most basic morphological operation.Suppose a given morphological structure element B, it can be defined as respectively expansion and the erosion operation of signal X:
For primary electrocardiosignal, dilation operation (formula 1) is to be used to remove the crest that frequency is higher than the electrocardio composition of baseline drift, erosion operation (formula 2) then is to be used to remove the trough that frequency is higher than the electrocardio composition of baseline drift, and expansion and corrosive resultant effect are exactly the baseline drift that obtains in the electrocardiosignal.Therefore, can be with the procedure definition of removing baseline drift:
The character of electrocardiosignal wavelet transform
Discrete wavelet coefficient X for discrete signal x (n)
j(n) can obtain by a twin-channel orthogonal filter group:
X
j(n)=h
j(n)*x(n),j=1,…,J (4)
Here J is the sum of change of scale; Utilize high pass filter h
1(n) and low pass filter h
0(n), can calculate filter coefficient h under each yardstick
j(n):
h
1(n)=h
0(2n)*h
0(2
2n)*h
0(2
3n)*…*h
0(2
J-1n)
h
2(n)=h
1(2n)*h
0(2
2n)*h
0(2
3n)*…*h
0(2
J-1n)
h
3(n)=h
1(2n) * h
0(2
2N) * h
0(2
3N) * ... * h
0(2
J-2N) (5) so, public affairs
.
.
.
h
J(n)=h
1(2n)
Wavelet coefficient in the formula (4) under each yardstick can be expressed as with the form of matrix and vector:
X
j=H
jx,j=1,…,J (6)
The H here
jBe convolution type matrix, x is a signal vector.
For given electrocardiosignal x, on time domain, the size of its power spectral density is directly proportional with the inverse of frequency (1/f), so the variance of electrocardio wavelet coefficient can be expressed as:
Here σ
2Be and signal variance
And analysis filter h
0(n) Xiang Guan constant, γ is a scale parameter, and 0≤γ≤2.
Taken the logarithm in formula (7) both sides, the logarithm value that can obtain wavelet coefficient is a linear change, that is:
Here α=log
2(σ
2) be defined as energy parameter.
The evolutive wavelet wiener deconvolution technology
Dimension is received deconvolution techniques and is normally used in the analysis of spectrum of optical signalling.The present invention is deconvoluted method through improving with traditional Wei Na, forms the distinctive evolution type dimension of wavelet field and receives deconvolution techniques, makes full use of wavelet transformation in the spatial time frequency resolution characteristic of Besov with this, realizes the abundant enhancing of characteristic information in the electrocardiosignal.The core concept of evolutive wavelet wiener deconvolution technology is described below.
Arbitrary signal y (t) in the real world can regard the result of useful signal x (t) and interfering signal g (t) convolution as, can be expressed as:
The result of convolution has caused signal waveform to broaden, and this has just brought overlapping each other of adjacent waveform, thereby causes the fuzzy of signal characteristic information.
In order to obtain useful signal x (t), need find a function to offset the effect of interference function g (t), just find out a function
Satisfy following relation (δ (t) is a uni-impulse function)
So can obtain formula (11) by formula (9) and (10).
For the convenience that further describes, with continuous signal x (t), the g (t) of above-mentioned appearance,
, y (t) and δ (t) be rewritten as corresponding discrete form x (n), g (n),
, y (n) and δ (n).Like this, just can be with the form representation formula (10) of matrix and vector.
Here
It is the operator that deconvolutes
Vector representation, among the present invention, in order to improve algorithm effects and efficient, will
Be designed to the Wei Na operator that deconvolutes.G is the matrix notation of g (n), and δ=[0 ... 010 ... 0]
T
Corresponding, formula (11) can be rewritten as:
Here
Be the vector representation of x (n), and Y is the matrix notation of y (n).
Thus, utilize formula (6) and formula (13) can obtain wavelet coefficient under each yardstick:
In order to obtain the optimum Wei Na operator that deconvolutes
The evolution generative process in defined two error criterions, can obtain first error criterion by formula (12)
Be defined as follows:
Then, utilize the character (formula 8) of the wavelet coefficient of 1/f signal can obtain another error criterion corresponding to each yardstick
J=1 ..., J
Like this, the error criterion by application of formula (15) and formula (16) can obtain the gross energy error criterion under each yardstick:
By finding the solution the minima of gross energy error, just can obtain the optimum Wei Na operator that deconvolutes
The gross energy error
Minima can find the solution by gradient base minimization technique.In order to improve algorithm the convergence speed, adopted steepest descent gradient base algorithm, that is: among the present invention
Cardioelectric characteristic extracting process of the present invention specifically describes as follows:
1) removes baseline drift.Adopt the baseline drift in the Mathematical Morphology technology removal original electrocardiographicdigital signal.The selection of morphological structure operator should be decided according to the sample frequency of electrocardiosignal.The present invention is directed to the electrocardiosignal of 360Hz sample frequency, get the reef knot constitutive element of 125ms.
2) parameter initialization:
A) initialization Wiener filter
Get base and float electrocardiosignal after the filtering, carry out the wavelet transformation under the 1st yardstick, obtain the useful signal of low frequency part, and the noise signal of HFS, obtain the power spectrum of useful signal and noise signal and crosspower spectrum between the two, construct the non-causal Wiener filter with this, as the initialization wave filter
C) initialization iteration step length β
p: for reducing operand, the selection iteration step length is a constant, β
pSelection influence convergence of algorithm speed, the excessive algorithm that causes is not restrained, too small then algorithm the convergence speed is slack-off, through experiment in a large number, selects β
p=0.125 can guarantee algorithmic statement, can guarantee signal smoothing after the iteration again.
3) feature of QRS complex wave strengthens:
A) wavelet transform.Adopt the quadrature discrete small echo to decompose pretreated electrocardiosignal, according to the frequency range of electrocardiosignal and QRS complex wave, the wavelet transformation that carries out suitable number of times is to obtain the feature sub-band at QRS wave group place;
B) wavelet field evolution type Wei Na deconvolutes.With the Wei Na operator vector that deconvolutes
Act on the sub-band at QRS wave group place, then calculate the gross energy error
If error amount is less than setting threshold, i.e. ε≤ε
0, stop computing; Otherwise, upgrade the Wei Na operator vector that deconvolutes
Energy parameter
And scale parameter
Continue computing up to finding the satisfied Wei Na operator vector that deconvolutes
4) feature of P ripple and T ripple strengthens:
A) remove the QRS wave group.Utilize the positional information of the above-mentioned QRS complex wave that obtains, the QRS wave group is substituted with baseline;
B) wavelet transform.Adopt the quadrature discrete small echo to decompose the electrocardiosignal that removes the QRS wave group, according to gained electrocardiosignal and P ripple, T wave frequency scope, the wavelet transformation that carries out suitable number of times is to obtain the feature sub-band at P ripple and T ripple place;
C) wavelet field evolution type Wei Na deconvolutes.Adopt the algorithm in the step 3), up to obtaining the satisfied Wei Na operator vector that deconvolutes
5) feature locations point detects:
A) QRS wave group feature locations point detects: get the QRS sub-band signal after feature that step 3) obtains strengthens, at first detect the R peak position, then according to the R peak position, find out first valley point and peak dot before and after the R peak.Then according to following criterion location feature point: if there is valley point of a peak dot at the R peak, peak dot is a QRS wave group starting point so, and the valley point is a Q ripple flex point; If have only a valley point before the R peak, this wave group is the RS wave group so, and this valley point is a RS wave group starting point; Valley point behind the R peak is a S ripple flex point, and peak dot is the wave group terminating point.
B) P, T wave characteristic location point detect: get P ripple, T marble band signal after feature that step 4) obtains strengthens, detect the peak point position of P ripple and T ripple, then find out first preceding valley point of P ripple and T crest value point respectively, be the starting point of P ripple and T ripple; Find out first valley point behind P ripple and the T crest value point, be the terminating point of P ripple and T ripple.
Fig. 2 adopts Mathematical Morphology Method to carry out the result schematic diagram that base floats filtering, and the present invention is directed to sample frequency is the 360Hz electrocardiosignal, has adopted 128 milliseconds reef knot constitutive element.
Fig. 3 is the feature enhanced results sketch map of QRS complex wave, and Fig. 3 (a) is the figure that No. 116 is write down after the process base floats filtering among the MIT-BIH arrhythmia data base, the position of figure acceptance of the bid clear R ripple, P ripple and T ripple.Receive deconvolution techniques by adopting the evolution type dimension, the feature of QRS wave group has been strengthened (Fig. 3 (b)) greatly, the amplitude of P ripple and T ripple has then obviously been cut down: the R ripple becomes more sharp-pointed, the starting point of QRS wave group becomes first valley point before the R peak, and the terminal point of QRS wave group becomes tangible peak dot behind the R peak; P ripple and T ripple are eliminated substantially.
Fig. 4 is P ripple and T wave characteristic enhanced results sketch map.Fig. 4 (a) is that No. 116 record floats filtering through base among the arrhythmia data base, and the figure of QRS wave group after removing.Also provided the position of R ripple, P ripple and T ripple among the figure, receive deconvolution techniques by adopting the evolution type dimension, the feature of P ripple and T ripple has been strengthened (Fig. 4 (b)) greatly: P ripple and T ripple become more sharp-pointed, their starting point becomes first valley point before P crest value point or the T crest value point, and their terminating point becomes QRS ripple P crest value point or T crest value point first valley point afterwards.
Among the present invention, adopt the evolution type dimension to receive deconvolution techniques and carry out the data that QRS ripple and P wave characteristic extract and take from MIT-BIH arrhythmia data base, the result of feature extraction and the performance comparison of carrying out feature extraction with traditional wavelet thereof are as shown in Table 1 and Table 2.When the T wave characteristic is extracted, detection performance for important parameter relevant in the evaluate cardiac disease detection-QT interval (characterizing ventricle polarization and unpolarized total time) with T ripple terminal point, adopt the QT data base as Data Source, the testing result of QT interval is as shown in table 3.
Table 1 QRS complex wave feature extraction result
Table 2 P wave characteristic is extracted the result
Table 3 QT interval testing result
Claims (2)
1. the cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution is characterized in that, comprises the following steps:
(1) according to the sample frequency of electrocardiosignal, selects the morphological structure operator, utilize the baseline drift in the Mathematical Morphology Method removal original electrocardiographicdigital signal;
(2) initialization Wiener filter, scale parameter
Energy parameter
And iteration step length (β
p), get base and float electrocardiosignal after the filtering, the ground floor that carries out wavelet transformation decomposes, obtain the useful signal of low frequency part under the 1st yardstick, and the noise signal of HFS, obtain the power spectrum of useful signal and noise signal and crosspower spectrum between the two, construct the non-causal Wiener filter, as the initialization Wei Na operator vector that deconvolutes with this
(3) QRS complex wave feature strengthens: adopt the quadrature discrete small echo to decompose pretreated electrocardiosignal, according to the frequency range of electrocardiosignal and QRS complex wave, carry out two-layer or two-layer above wavelet transformation to obtain the feature sub-band at QRS wave group place; With the initialization Wei Na operator vector that deconvolutes
Act on the sub-band at QRS wave group place, then calculate the gross energy error
If error amount less than setting threshold, stops computing; Otherwise, utilize steepest descent gradient based method, upgrade the Wei Na operator vector that deconvolutes
Energy parameter
And scale parameter
Continue computing up to error amount less than setting threshold, promptly find the satisfied Wei Na operator vector that deconvolutes
(4) P ripple and T wave characteristic strengthen: utilize the positional information of the QRS complex wave that step (3) obtains, the QRS wave group is substituted with baseline; Adopt the quadrature discrete small echo to decompose the electrocardiosignal that removes the QRS wave group,, carry out two-layer or two-layer above wavelet transformation to obtain the feature sub-band at P ripple and T ripple place according to gained electrocardiosignal and P ripple, T wave frequency scope; With the initialization Wei Na operator vector that deconvolutes
Act on the feature sub-band at P ripple and T ripple place, then calculate the gross energy error
If error amount less than setting threshold, stops computing; Otherwise, utilize steepest descent gradient based method, upgrade the Wei Na operator vector that deconvolutes
Energy parameter
And scale parameter
Continue computing up to finding the satisfied Wei Na operator vector that deconvolutes
(5) detect QRS wave group feature locations point: get the QRS sub-band signal after feature that step (3) obtains strengthens, at first detect the R peak position,, find out before the R peak first valley point and peak dot behind first valley point and peak dot and R peak then according to the R peak position; Then according to following criterion location feature point: if valley point of a peak dot is arranged before the R peak, peak dot is a QRS wave group starting point so, and the valley point is a Q ripple flex point; If have only a valley point before the R peak, this wave group is the RS wave group so, and this valley point is a RS wave group starting point; Valley point behind the R peak is a S ripple flex point, and peak dot is the wave group terminating point;
(6) detect P, T wave characteristic location point: get P ripple, T marble band signal after feature that step (4) obtains strengthens, detect the peak point position of P ripple and T ripple, then find out first preceding valley point of P ripple and T crest value point respectively, be the starting point of P ripple and T ripple; Find out first valley point behind P ripple and the T crest value point, be the terminating point of P ripple and T ripple.
2. the cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution according to claim 1 is characterized in that,
The formula that described steepest descent gradient based method adopts is:
In the formula,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100580219A CN100536765C (en) | 2007-07-13 | 2007-07-13 | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100580219A CN100536765C (en) | 2007-07-13 | 2007-07-13 | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101088456A CN101088456A (en) | 2007-12-19 |
CN100536765C true CN100536765C (en) | 2009-09-09 |
Family
ID=38942075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100580219A Expired - Fee Related CN100536765C (en) | 2007-07-13 | 2007-07-13 | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100536765C (en) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101843480B (en) * | 2009-03-27 | 2013-04-24 | 华为技术有限公司 | Device for processing bioelectrical signals |
CN101828917B (en) * | 2010-05-07 | 2011-09-14 | 深圳大学 | Method and system for extracting electrocardiosignal characteristic |
CN101866423B (en) * | 2010-05-12 | 2012-07-25 | 重庆大学 | Method for dividing arterial pressure signal beat by beat |
CN102160787B (en) * | 2011-05-31 | 2013-06-26 | 广东工业大学 | Time-frequency-transformation-based blind extraction method of fetal electrocardiography |
CN103197017A (en) * | 2013-05-02 | 2013-07-10 | 云南烟草科学研究院 | Cigarette smoke colour spectrum baseline correction MPLS (mathematically proven learning system) method |
CN104783786B (en) * | 2014-01-20 | 2017-08-08 | 中国科学院苏州纳米技术与纳米仿生研究所 | Characteristics of electrocardiogram identifying system and method |
CN107072545B (en) * | 2014-09-02 | 2020-01-31 | 吉莫·博尔只斤 | Electrocardiogram data analysis method and system for rapid diagnosis |
CN104783780B (en) * | 2015-04-13 | 2017-11-24 | 深圳市飞马与星月科技研究有限公司 | ECG De method and device |
CN105054926A (en) * | 2015-04-13 | 2015-11-18 | 深圳市飞马与星月科技研究有限公司 | Electrocardiosignal feature information extraction method and device |
WO2016168979A1 (en) | 2015-04-20 | 2016-10-27 | 深圳市长桑技术有限公司 | Vital sign analysis method and system |
US9706938B2 (en) * | 2015-10-29 | 2017-07-18 | Medtronic Monitoring, Inc. | System and method to determine premature ventricular contraction (PVC) type and burden |
JP6646755B2 (en) * | 2016-03-01 | 2020-02-14 | セント・ジュード・メディカル,カーディオロジー・ディヴィジョン,インコーポレイテッド | Method and system for operating a system for mapping cardiac activity |
CN107766864B (en) * | 2016-08-23 | 2022-02-01 | 斑马智行网络(香港)有限公司 | Method and device for extracting features and method and device for object recognition |
CN107203692B (en) * | 2017-05-09 | 2020-05-05 | 哈尔滨工业大学(威海) | Electrocardio data digital signal processing method based on deep convolutional neural network |
CN107184198A (en) * | 2017-06-01 | 2017-09-22 | 广州城市职业学院 | A kind of electrocardiosignal classifying identification method |
CN107874755B (en) * | 2017-09-12 | 2019-09-03 | 广东工业大学 | A kind of P wave extracting method based on rectangular projection and independent component analysis |
CN109984740A (en) * | 2018-01-03 | 2019-07-09 | 四川锦江电子科技有限公司 | A kind of device and method positioning ecg characteristics in body surface ECG data |
CN108888263B (en) * | 2018-05-22 | 2020-10-30 | 郑州大学 | R wave detection method based on geometric form group characteristics |
CN108983279B (en) * | 2018-07-05 | 2021-12-21 | 南京航空航天大学 | Low-background spectrum resolving method based on sodium iodide detector |
CN109199355B (en) * | 2018-09-18 | 2021-09-28 | 深圳和而泰数据资源与云技术有限公司 | Heart rate information detection method and device and detection equipment |
CN110179455A (en) * | 2019-04-30 | 2019-08-30 | 厦门大学 | A kind of ECG pretreatment and waveforms detection method and apparatus based on sudden death early warning |
CN111134659B (en) * | 2020-01-17 | 2022-04-26 | 上海乐普云智科技股份有限公司 | Method and device for detecting P wave and T wave in electrocardiosignal |
-
2007
- 2007-07-13 CN CNB2007100580219A patent/CN100536765C/en not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
WAVELET DOMAIN WIENER FILTERING FORECGDENOISING USING IMPROVED SIGNAL ESTIMATE. N. Nikolaev, Z. Nikolov, A. Gotchev, K. Egiazarian.IEEE International Conference on 5-9 June 2000,Vol.6 . 2000 |
WAVELET DOMAIN WIENER FILTERING FORECGDENOISING USING IMPROVED SIGNAL ESTIMATE. N. Nikolaev, Z. Nikolov, A. Gotchev, K. Egiazarian.IEEE International Conference on 5-9 June 2000,Vol.6 . 2000 * |
小波域Wiener 滤波器信号的去噪方法. 李锵,滕建辅,李士心,肖志涛.天津大学学报,第36卷第3期. 2003 |
小波域Wiener 滤波器信号的去噪方法. 李锵,滕建辅,李士心,肖志涛.天津大学学报,第36卷第3期. 2003 * |
Also Published As
Publication number | Publication date |
---|---|
CN101088456A (en) | 2007-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100536765C (en) | Cardioelectric characteristic extracting process based on evolutive wavelet wiener deconvolution | |
CN108714026B (en) | Fine-grained electrocardiosignal classification method based on deep convolutional neural network and online decision fusion | |
CN110013247B (en) | Detection, differentiation and positioning method for electrocardiogram P wave and T wave | |
CN103405227B (en) | Double-layer morphological filter based electrocardiosignal preprocessing method | |
CN100418480C (en) | Heart disease automatic classification system based on heart sound analysis and heart sound segmentation method | |
Sabherwal et al. | Automatic detection of the R peaks in single-lead ECG signal | |
CN105902263B (en) | Assess the method and device of the RR interphase degree of irregularity of electrocardiographic wave | |
CN103083013A (en) | Electrocardio signal QRS complex wave detection method based on morphology and wavelet transform | |
CN102626310A (en) | Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving | |
CN108309284A (en) | Electrocardiogram T wave end point detection method and device | |
CN112932498B (en) | T waveform state classification system with generalization capability based on deep learning | |
CN108647584B (en) | Arrhythmia identification and classification method based on sparse representation and neural network | |
Lastre-Dominguez et al. | Denoising and features extraction of ecg signals in state space using unbiased fir smoothing | |
Tang et al. | ECG de-noising based on empirical mode decomposition | |
CN113951893B (en) | Multi-lead electrocardiosignal characteristic point extraction method combining deep learning and electrophysiological knowledge | |
Rai et al. | ECG signal classification using wavelet transform and Back Propagation Neural Network | |
CN109124620A (en) | A kind of atrial fibrillation detection method, device and equipment | |
CN111419213A (en) | ECG electrocardiosignal generation method based on deep learning | |
CN103190901A (en) | R wave detection algorithm based on extremum field mean mode decomposition and improved Hilbert enveloping | |
CN111643070B (en) | T wave starting point determining method and device, storage medium and electronic equipment | |
CN114469124A (en) | Method for identifying abnormal electrocardiosignals in motion process | |
Li et al. | A novel abnormal ECG beats detection method | |
Rai et al. | Classification of ECG waveforms for abnormalities detection using DWT and back propagation algorithm | |
CN113598785B (en) | Electrocardiosignal denoising method based on wavelet domain sparse characteristic | |
CN115177268A (en) | Atrial fibrillation heart beat detection method, device and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090909 Termination date: 20100713 |