CN106963368A - QRS complex detection method and device - Google Patents
QRS complex detection method and device Download PDFInfo
- Publication number
- CN106963368A CN106963368A CN201710067025.7A CN201710067025A CN106963368A CN 106963368 A CN106963368 A CN 106963368A CN 201710067025 A CN201710067025 A CN 201710067025A CN 106963368 A CN106963368 A CN 106963368A
- Authority
- CN
- China
- Prior art keywords
- qrs complex
- new
- electrocardiosignal
- threshold value
- slope1
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/366—Detecting abnormal QRS complex, e.g. widening
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Abstract
The present invention discloses a kind of QRS complex detection method and device, and wherein QRS complex detection method comprises the following steps：Electrocardiosignal is obtained, the signal of acquisition is pre-processed；Core one section of beginning calculating electrocardiosignal initial slope threshold value of electric signal；Electrocardiosignal is detected according to initial slope threshold value, the position of QRS complex is obtained；Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.The QRS complex detection method and device of the proposition of the present invention, upgraded in time QRS complex judgment threshold using QRS complex, compared with prior art, substantially increase judging nicety rate, amount of calculation has been greatly reduced simultaneously, anti-dry RR immunity is strong, improve treatment effeciency, adaptivity is stronger.
Description
Technical field
The present invention relates to signal processing technology in Medical Devices, more particularly to a kind of QRS complex detection method and device.
Background technology
Cardiovascular disease incidence rate more and more higher, it has also become constitute one of maximum harm disease to human life at present.
And electrocardiogram is detection, diagnoses and prevent the main effective ways of this kind of disease.Traditional PC patient monitors are expensive, volume
Huge, inconvenience is mobile and is concentrated mainly on hospital, and can not real-time monitoring patient the state of an illness, brought very very much not to patient and doctor
Just.For these drawbacks, also with developing rapidly for embedded and network communications technology, traditional mobile electrocardioscanner is overcome
Arise at the historic moment, it has low cost, small volume, reliability height, simple operation and other advantages, it is adaptable to personal, family, middle-size and small-size doctor
Institute and community medical therapy unit.Mobile ecg telemonitoring equipment can quickly carry out detection and heart disease in real time to ambulatory ecg signal
Automatic diagnosis, inevitable so as to which a variety of services such as urgent rescue, Disease Warning Mechanism, medical consultation and guidance can be provided the user, choosing
A kind of accurate heart rhythm detection algorithms being quickly applicable are selected, it is significant.
Detection QRS complex algorithm mainly had following several in recent years：
（1）Threshold detection method based on slope, amplitude and width characteristics value：The algorithm is simpler, and processing speed is fast.It is not enough：
Accuracy rate needs further raising.
（2）Detection method based on wavelet transformation：This method is used for the detection of electrocardiosignal, has obtained preferable effect, small
There is wave conversion good time-frequency domain to change and localization analysis ability, but its deficiency is：Amount of calculation is than larger, processing high-volume
Data age rate is not high；
（3）Detection method based on template matches：This method principle is electrocardiosignal and the waveform template pointwise prestored
Compare, when pending signal is coupled with template, correlation is maximum, it is not enough：This method easily by high-frequency noise and baseline drift and
The influence of individual difference, shows poor in terms of adaptive ability；
（4）Detection method based on neutral net：This method is used as a kind of simulation human brain mechanism and work(developed rapidly in recent years
The novel computer and artificial intelligence technology of energy, mainly as self-adaptation nonlinear fallout predictor and in QRS complex detection
Vector quantization network is practised, there is faster detection speed and higher verification and measurement ratio, but its deficiency is：This method is by artificial network's
Learning training, and popularity to training sample and it is representative have certain requirement, training study will spend the more time, real
Border application difficult；
（5）Detection method based on mathematical morphology：This method extracts signal component（Crest, trough）Form be by tying
Constitutive element is determined, high-frequency noise is removed to electrocardiosignal first with parallel opening and closing, make and break computing, then electrocardiosignal is passed through
PVE methods carry out peak valley and extract computing to detect QRS complex, and the QRS complex detection architecture that this method is constituted has mathematics
On preciseness and accuracy, have extraordinary inhibition to low-frequency noises such as baseline drifts, but its deficiency is：This method pair
High-frequency noise is sensitive, and the efficiency of algorithm is confined to the length of structural element and the sample rate of signal etc..
The content of the invention
The main object of the present invention is to disclose a kind of QRS complex detection method and device, is mainly used in Medical Devices
In, to solve that current QRS complex detection method accuracy rate is low, computationally intensive, treatment effeciency is low, the adaptive technology such as poor asks
Topic.
In order to solve the above technical problems, one aspect of the present invention proposes a kind of QRS complex detection method, including following step
Suddenly：
Electrocardiosignal is obtained, the signal of acquisition is pre-processed；
Core one section of beginning calculating electrocardiosignal initial slope threshold value of electric signal；
Electrocardiosignal is detected according to initial slope threshold value, the position of QRS complex is obtained；
Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.
Preferably, in above-mentioned QRS complex detection method, one section of calculating electrocardiosignal of the beginning for electric signal of coring
Initial slope threshold value；Electrocardiosignal is detected according to initial slope threshold value, the position specific steps of QRS complex are obtained such as
Under：
10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, calculates the difference maximum value of each signal segment；
Remove one minimum differential absolute value of a maximum, the average for remembering remaining difference absolute value is avgSlope；
Calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*avgSlope；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
Preferably, in above-mentioned QRS complex detection method, after n QRS complex of acquisition, according to the QRS of acquisition
Complex wave updates slope threshold value；Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained
Comprise the following steps that：
N QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
The average for calculating difference maximum value is avgSlope_New；
It is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New=2/ to calculate the slope threshold value after updating
9*avgSlope_New；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, work as slope1_
It is each before and after slope is slope1_New point when New >=SlopeThr1_New, slope2_New >=SlopeThr2_New
The spike point position of maximum, as QRS complex position are searched out in 0.6s electrocardiosignal.
Preferably, in above-mentioned QRS complex detection method, after taking n QRS complex, according to QRS complex more
New slope threshold value；Electrocardiosignal is detected according to the slope threshold value after renewal, obtaining the position of QRS complex also includes step
Suddenly：
The average of QRS complex peak value, the i.e. QRS complex for subsequent analysis are calculated from the QRS complex of n acquisition
The judgment threshold avgQRSPeak of peak value；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, and meet
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, in the point that slope is slope1_New
The spike point position of maximum is searched out in front and rear each 0.6s electrocardiosignal, then it is assumed that be tentative QRS complex position, and remember
Its lower peak value is QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then it is assumed that should
Position is the position of QRS complex.
Preferably, in above-mentioned QRS complex detection method, step is also included after acquisition QRS complex：
Whether judge QRS complex is interference signal section；
Reject the QRS complex for being judged as interference signal section.
Meanwhile, invention additionally discloses a kind of QRS complex detection means, including：
Acquisition module, for obtaining electrocardiosignal, is pre-processed to the signal of acquisition；
Computing module, one section of calculating electrocardiosignal initial slope threshold value for the beginning for electric signal of coring；
Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
Judge module, for being detected according to initial slope threshold value to electrocardiosignal, obtains the position of QRS complex；
Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.
Preferably, in above-mentioned QRS complex detection means, the computing module is cored one section of beginning of electric signal
Calculate electrocardiosignal initial slope threshold value；The judge module is detected according to initial slope threshold value to electrocardiosignal, is obtained
The position of QRS complex is comprised the following steps that：
10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, calculates the difference maximum value of each signal segment；
Remove one minimum differential absolute value of a maximum, the average for remembering remaining difference absolute value is avgSlope；
Calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*avgSlope；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
Preferably, in above-mentioned QRS complex detection means, the computing module is obtained after n QRS complex,
Slope threshold value is updated according to the QRS complex of acquisition；The judge module enters according to the slope threshold value after renewal to electrocardiosignal
Row detection, the position for obtaining QRS complex is comprised the following steps that：
N QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
The average for calculating difference maximum value is avgSlope_New；
It is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New=2/ to calculate the slope threshold value after updating
9*avgSlope_New；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, work as slope1_
It is each before and after slope is slope1_New point when New >=SlopeThr1_New, slope2_New >=SlopeThr2_New
The spike point position of maximum, as QRS complex position are searched out in 0.6s electrocardiosignal.
Preferably, in above-mentioned QRS complex detection means, the computing module is additionally operable to answer from the QRS of n acquisition
Calculate the average of QRS complex peak value in multiplex, i.e., for subsequent analysis QRS complex peak value judgment threshold
avgQRSPeak；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, and meet
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, in the point that slope is slope1_New
The spike point position of maximum is searched out in front and rear each 0.6s electrocardiosignal, then it is assumed that be tentative QRS complex position, and remember
Its lower peak value is QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then it is assumed that should
Position is the position of QRS complex.
Preferably, in above-mentioned QRS complex detection means, the QRS complex detection means also includes screening mould
Block, the screening module is used for after obtaining QRS complex, judges whether QRS complex is interference signal section and rejects judgement
For the QRS complex of interference signal section.
The QRS complex detection method and device of the proposition of the present invention, are upgraded in time QRS complex using QRS complex
Judgment threshold, is compared with prior art, substantially increases judging nicety rate, while being greatly reduced amount of calculation, having resisted and do
Immunity is strong, improve treatment effeciency, and adaptivity is stronger.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Structure according to these accompanying drawings obtains other accompanying drawings.
Fig. 1 is the schematic flow sheet of one embodiment in QRS complex detection method of the present invention；
Fig. 2 is the schematic flow sheet in QRS complex detection method of the present invention according to initial slope threshold test QRS complex；
Fig. 3 is combined to obtain renewal slope threshold value detection QRS after n QRS complex in QRS complex detection method of the present invention
The schematic flow sheet of ripple；
Fig. 4 is the structural representation of one embodiment in QRS complex detection means of the present invention；
Fig. 5 is the structural representation of another embodiment in QRS complex detection means of the present invention.
The realization, functional characteristics and advantage of the object of the invention will be described further referring to the drawings in conjunction with the embodiments.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the present invention, rather than whole embodiments.Base
Embodiment in the present invention, those of ordinary skill in the art obtained under the premise of creative work is not made it is all its
His embodiment, belongs to the scope of protection of the invention.
Indicated it is to be appreciated that institute is directional in the embodiment of the present invention（Such as up, down, left, right, before and after ...）Only use
In explanation in a certain particular pose（As shown in drawings）Under relative position relation, motion conditions between each part etc., if should
When particular pose changes, then directionality indicates also correspondingly therewith to change.
In addition, the description for being related to " first ", " second " etc. in the present invention is only used for describing purpose, and it is not intended that referring to
Show or imply its relative importance or the implicit quantity for indicating indicated technical characteristic.Thus, " first ", " are defined
At least one this feature can be expressed or be implicitly included to two " feature.In addition, the technical scheme between each embodiment can
To be combined with each other, but must can be implemented as basis with those of ordinary skill in the art, when the combination outlet of technical scheme
It is conflicting or when can not realize it will be understood that the combination of this technical scheme is not present, also not in the protection model of application claims
Within enclosing.
First, the invention discloses a kind of QRS complex detection method, reference picture 1, QRS complex disclosed by the invention
Detection method comprises the following steps：
11st, electrocardiosignal is obtained, the signal of acquisition is pre-processed；
12nd, core one section of beginning calculating electrocardiosignal initial slope threshold value of electric signal；
13rd, electrocardiosignal is detected according to initial slope threshold value, obtains the position of QRS complex；
14th, obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
15th, electrocardiosignal is detected according to the slope threshold value after renewal, obtains the position of QRS complex.
Obtain after electrocardiosignal, it usually needs the signal of acquisition is pre-processed, to reduce signal noise, improve detection
In accuracy rate, the embodiment, the step of being pre-processed to signal includes but is not limited to 50Hz notch filters and 1-30Hz bands
Pass filter.
In reference picture 2, the technical scheme of the invention, one section of calculating electrocardiosignal initial slope of the beginning for electric signal of coring
Threshold value；Electrocardiosignal is detected according to initial slope threshold value, the position for obtaining QRS complex is comprised the following steps that：
121st, 10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, and the difference for calculating each signal segment is maximum definitely
Value；
122nd, one minimum differential absolute value of a maximum is removed, the average for remembering remaining difference absolute value is avgSlope；
123rd, calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*
avgSlope；
131st, obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
In above-described embodiment, by the electrocardiosignal initial segment 10s to acquisition and 5 2s signal segment is divided into, meter
Calculate the difference maximum value of each signal segment, Difference Calculation formula：
slope (i) = ecg (i+2) – ecg (i)
In the embodiment, by calculating the initial slope threshold value of electrocardiosignal, and by continuously connecting a slope on electrocardiosignal
Judged, when meeting default threshold condition, i.e. slope1 >=SlopeThr1, slope2 >=SlopeThr2 is in slope
The spike point position of maximum, as QRS complex position are searched out before and after slope1 point in each 0.6s electrocardiosignal.
In reference picture 3, the technical scheme of the invention, after n QRS complex is obtained using the above method, according to obtaining
The QRS complex taken updates slope threshold value；Electrocardiosignal is detected according to the slope threshold value after renewal, QRS is obtained and is combined
The position of ripple is comprised the following steps that：
141st, n QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
142nd, the average for calculating difference maximum value is avgSlope_New；
143rd, it is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New to calculate the slope threshold value after updating
= 2/9*avgSlope_New；
151st, obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, when
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, in the point that slope is slope1_New
The spike point position of maximum, as QRS complex position are searched out in front and rear each 0.6s electrocardiosignal.
In the embodiment, 8 QRS complexes are generally obtained, the difference that each QRS complex is calculated afterwards is maximum definitely
Value, calculation formula is as follows：
slope(i) = ecg(i+2) – ecg (i)
In the embodiment, slope threshold value is constantly updated using the QRS complex of acquisition so that follow-up QRS complex detection is accurate
True rate is higher.
In above-described embodiment, the absolute refractory period of a period of time is set after QRS complex is detected, usually 200ms,
During this period without the detection of QRS complex, the refractory period is set according to the electrocardiosignal cycle, in refractory period, generally not
It is possible that QRS complex, the amount of calculation of technical solution of the present invention can be substantially reduced by the technical method, detection is improved
Efficiency.
As another embodiment of the present invention, take after n QRS complex, slope threshold value is updated according to QRS complex；
Electrocardiosignal is detected according to the slope threshold value after renewal, obtaining the position of QRS complex also includes step：
The judgment threshold avgRR of phase between the RR of phase average, i.e. subsequent analysis between RR is calculated from the QRS complex of n acquisition；
The QRS complex peak value of the average, i.e. subsequent analysis of QRS complex peak value is calculated from the QRS complex of n acquisition
Judgment threshold avgQRSPeak；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, and meet
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, it is slope1_New's to write down slope
Point peak value is QRSPeak, works as RR>200ms, and QRSPeak>It is slope1_New's in slope during avgQRSPeak/2
The spike point position of maximum, as QRS complex position are searched out before and after point in each 0.6s electrocardiosignal.
The embodiment also add phase and QRS complex peak between RR while test point progress slope threshold value judgement is treated
The judgment threshold of value, is meeting RR>200ms, and QRSPeak>It is slope1_New's in slope during avgQRSPeak/2
The spike point position of maximum is searched out before and after point in each 0.6s electrocardiosignal, technical solution of the present invention detection QRS is greatly improved
The accuracy rate of complex wave.
In said one or multiple embodiments, also include step after acquisition QRS complex：
Whether judge QRS complex is interference signal section；
Reject the QRS complex for being judged as interference signal section.
The embodiment eliminates the part for belonging to interference signal in QRS complex, directly enhances the detection of QRS complex
Accuracy rate.
In technical scheme, after a QRS complex is detected, if the phase is longer between current RR, RR is met
>Then it is probably to occur in that missing inspection during 1.5*avgRR, takes an interim segment signal between current RR, QRS will be re-started to it and answered
Multiplex is detected, its corresponding slope judgment threshold is all reduced into half, when the slope of continuous two points meets slope threshold conditions
When, position of the spike point position of maximum for tentative QRS complex is searched out before and after the point in an each 0.6s segment signal,
And its peak value is write down for QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then
It is the physical location of QRS complex to think the position, and judges whether the signal segment is interference signal section simultaneously, conversely, then not
It is the position of QRS complex.
The QRS complex detection method of the proposition of the present invention, threshold is judged using the QRS complex QRS complex that upgrades in time
Value, compared with prior art, substantially increase judging nicety rate, while be greatly reduced amount of calculation, strong interference immunity,
Treatment effeciency is improved, adaptivity is stronger.
The present invention also provides a kind of QRS complex detection means that QRS complex detection is carried out using the above method
Reference picture 4, invention additionally discloses a kind of QRS complex detection means, including：
Acquisition module 21, for obtaining electrocardiosignal, is pre-processed to the signal of acquisition；
Computing module 22, one section of calculating electrocardiosignal initial slope threshold value for the beginning for electric signal of coring；
Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
Judge module 23, for being detected according to initial slope threshold value to electrocardiosignal, obtains the position of QRS complex；
Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.
Obtain after electrocardiosignal, acquisition module 21 usually requires to pre-process the signal of acquisition, is made an uproar with reducing signal
Sound, is improved in Detection accuracy, the embodiment, the step of being pre-processed to signal includes but is not limited to 50Hz notch filters
With 1-30Hz bandpass filterings.
Preferably, in above-mentioned QRS complex detection means, computing module 22 is cored one section of beginning meter of electric signal
Calculate electrocardiosignal initial slope threshold value；Judge module 23 is detected according to initial slope threshold value to electrocardiosignal, is obtained QRS and is answered
The position of multiplex is comprised the following steps that：
10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, calculates the difference maximum value of each signal segment；
Remove one minimum differential absolute value of a maximum, the average for remembering remaining difference absolute value is avgSlope；
Calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*avgSlope；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
In above-described embodiment, by the electrocardiosignal initial segment 10s to acquisition and 5 2s signal segment is divided into, meter
Calculate the difference maximum value of each signal segment, Difference Calculation formula：
slope (i) = ecg (i+2) – ecg (i)
In the embodiment, by calculating the initial slope threshold value of electrocardiosignal, and by continuously connecting a slope on electrocardiosignal
Judged, when meeting default threshold condition, i.e. slope1 >=SlopeThr1, slope2 >=SlopeThr2 is in slope
The spike point position of maximum, as QRS complex position are searched out before and after slope1 point in each 0.6s electrocardiosignal.
Preferably, in above-mentioned QRS complex detection means, computing module 22 is obtained after n QRS complex, root
Slope threshold value is updated according to the QRS complex of acquisition；Judge module 23 is examined according to the slope threshold value after renewal to electrocardiosignal
Survey, the position for obtaining QRS complex is comprised the following steps that：
N QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
The average for calculating difference maximum value is avgSlope_New；
It is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New=2/ to calculate the slope threshold value after updating
9*avgSlope_New；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, work as slope1_
It is each before and after slope is slope1_New point when New >=SlopeThr1_New, slope2_New >=SlopeThr2_New
The spike point position of maximum, as QRS complex position are searched out in 0.6s electrocardiosignal.
In the embodiment, 8 QRS complexes are generally obtained, the difference that each QRS complex is calculated afterwards is maximum definitely
Value, calculation formula is as follows：
slope(i) = ecg(i+2) – ecg (i)
In the embodiment, slope threshold value is constantly updated using the QRS complex of acquisition so that follow-up QRS complex detection is accurate
True rate is higher.
Preferably, in above-mentioned QRS complex detection means, computing module 22 is additionally operable to be combined from the QRS of n acquisition
The judgment threshold avgRR of phase between the RR of phase average, i.e. subsequent analysis between RR is calculated in ripple；
The QRS complex peak value of the average, i.e. subsequent analysis of QRS complex peak value is calculated from the QRS complex of n acquisition
Judgment threshold avgQRSPeak；
Judge module 23 be used for obtain electrocardiosignal on continuous 2 points and calculate respectively its slope for slope1_New,
Slope2_New, and slope1_New >=SlopeThr1_New is met, during slope2_New >=SlopeThr2_New, in slope
The spike point position of maximum is searched out in electrocardiosignal for each 0.6s before and after slope1_New point, then it is assumed that be tentative
QRS complex position, and its peak value is write down for QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>
AvgQRSPeak/2, then it is assumed that the position is the position of QRS complex.
Reference picture 5, in above-mentioned QRS complex detection means, QRS complex detection means also includes screening module
24, screening module is used for after obtaining QRS complex, judges whether QRS complex is interference signal section and rejecting is judged as doing
Disturb the QRS complex of signal segment.
Screening module 24 in the embodiment eliminates the part for belonging to interference signal in QRS complex, directly enhances
The accuracy rate of the detection of QRS complex.
In technical scheme, after a QRS complex is detected, if the phase is longer between current RR, RR is met>
Then it is probably to occur in that missing inspection during 1.5*avgRR, takes an interim segment signal between current RR, QRS will be re-started to it and be combined
Ripple is detected, its corresponding slope judgment threshold is all reduced into half, when the slope of continuous two points meets slope threshold conditions,
Position of the spike point position of maximum for tentative QRS complex is searched out before and after the point in an each 0.6s segment signal, and
Its peak value is write down for QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then it is assumed that the position
Put be QRS complex physical location, and judge simultaneously the signal segment whether be interference signal section, conversely, be not then QRS be combined
The position of ripple.
Each functional unit in each embodiment of the invention can be integrated in a processing unit or each
Unit is individually physically present, can also two or more units it is integrated in a unit.Above-mentioned integrated unit both may be used
To be realized in the form of hardware, it would however also be possible to employ the form of SFU software functional unit is realized.
If the integrated unit is realized using in the form of SFU software functional unit and as independent production marketing or used
When, it can be stored in a computer read/write memory medium.Understood based on such, technical scheme is substantially
The part contributed in other words to prior art or all or part of the technical scheme can be in the form of software products
Embody, the computer software product is stored in a storage medium, including some instructions are to cause a computer
Equipment (can be personal computer, server, or network equipment etc.) performs the complete of each embodiment methods described of the invention
Portion or part steps.And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can store journey
The medium of sequence code.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the scope of the invention, it is every to utilize
The equivalent structure transformation that description of the invention and accompanying drawing content are made, or directly or indirectly it is used in other related technology necks
Domain, is included within the scope of the present invention.
Claims (10)
1. a kind of QRS complex detection method, it is characterised in that comprise the following steps：
Electrocardiosignal is obtained, the signal of acquisition is pre-processed；
Core one section of beginning calculating electrocardiosignal initial slope threshold value of electric signal；
Electrocardiosignal is detected according to initial slope threshold value, the position of QRS complex is obtained；
Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.
2. QRS complex detection method as claimed in claim 1, it is characterised in that one section of meter of the beginning for electric signal of coring
Calculate electrocardiosignal initial slope threshold value；Electrocardiosignal is detected according to initial slope threshold value, the position of QRS complex is obtained
Comprise the following steps that：
10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, calculates the difference maximum value of each signal segment；
Remove one minimum differential absolute value of a maximum, the average for remembering remaining difference absolute value is avgSlope；
Calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*avgSlope；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
3. QRS complex detection method as claimed in claim 2, it is characterised in that after n QRS complex of acquisition, according to
The QRS complex of acquisition updates slope threshold value；Electrocardiosignal is detected according to the slope threshold value after renewal, QRS is obtained and answers
The position of multiplex is comprised the following steps that：
N QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
The average for calculating difference maximum value is avgSlope_New；
It is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New=2/ to calculate the slope threshold value after updating
9*avgSlope_New；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, work as slope1_
It is each before and after slope is slope1_New point when New >=SlopeThr1_New, slope2_New >=SlopeThr2_New
The spike point position of maximum, as QRS complex position are searched out in 0.6s electrocardiosignal.
4. QRS complex detection method as claimed in claim 3, it is characterised in that after taking n QRS complex, according to
QRS complex updates slope threshold value；Electrocardiosignal is detected according to the slope threshold value after renewal, QRS complex is obtained
Position also includes step：
The average of QRS complex peak value, the i.e. QRS complex for subsequent analysis are calculated from the QRS complex of n acquisition
The judgment threshold avgQRSPeak of peak value；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, and meet
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, in the point that slope is slope1_New
The spike point position of maximum is searched out in front and rear each 0.6s electrocardiosignal, then it is assumed that be tentative QRS complex position, and remember
Its lower peak value is QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then it is assumed that should
Position is the position of QRS complex.
5. QRS complex detection method as claimed in claim 2, it is characterised in that also include step after acquisition QRS complex
Suddenly：
Whether judge QRS complex is interference signal section；
Reject the QRS complex for being judged as interference signal section.
6. a kind of QRS complex detection means, it is characterised in that including：
Acquisition module, for obtaining electrocardiosignal, is pre-processed to the signal of acquisition；
Computing module, one section of calculating electrocardiosignal initial slope threshold value for the beginning for electric signal of coring；
Obtain after n QRS complex, slope threshold value is updated according to the QRS complex of acquisition；
Judge module, for being detected according to initial slope threshold value to electrocardiosignal, obtains the position of QRS complex；
Electrocardiosignal is detected according to the slope threshold value after renewal, the position of QRS complex is obtained.
7. QRS complex detection means as claimed in claim 6, it is characterised in that the computing module is cored electric signal
The one section of calculating electrocardiosignal initial slope threshold value started；The judge module is carried out according to initial slope threshold value to electrocardiosignal
Detection, the position for obtaining QRS complex is comprised the following steps that：
10s that electric signal of coring starts simultaneously is divided into 5 2s signal segments, calculates the difference maximum value of each signal segment；
Remove one minimum differential absolute value of a maximum, the average for remembering remaining difference absolute value is avgSlope；
Calculating initial slope threshold value is：SlopeThr1=2/7*avgSlope, SlopeThr2=2/11*avgSlope；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1, slope2, when slope1 >=
When SlopeThr1, slope2 >=SlopeThr2, searched out before and after slope is slope1 point in each 0.6s electrocardiosignal
Maximum spike point position, as QRS complex position.
8. QRS complex detection means as claimed in claim 7, it is characterised in that the computing module obtains n QRS and answered
After multiplex, slope threshold value is updated according to the QRS complex of acquisition；The judge module is according to the slope threshold value after renewal to the heart
Electric signal is detected that the position for obtaining QRS complex is comprised the following steps that：
N QRS complex is obtained, the difference maximum value of each QRS complex is calculated；
The average for calculating difference maximum value is avgSlope_New；
It is SlopeThr1_New=2/5*avgSlope_New, SlopeThr2_New=2/ to calculate the slope threshold value after updating
9*avgSlope_New；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, work as slope1_
It is each before and after slope is slope1_New point when New >=SlopeThr1_New, slope2_New >=SlopeThr2_New
The spike point position of maximum, as QRS complex position are searched out in 0.6s electrocardiosignal.
9. the QRS complex detection means as described in right 8, it is characterised in that the computing module is additionally operable to what is obtained from n
Calculate the average of QRS complex peak value in QRS complex, i.e., for subsequent analysis QRS complex peak value judgment threshold
avgQRSPeak；
Obtain on electrocardiosignal at continuous 2 points and calculate its slope respectively for slope1_New, slope2_New, and meet
When slope1_New >=SlopeThr1_New, slope2_New >=SlopeThr2_New, in the point that slope is slope1_New
The spike point position of maximum is searched out in front and rear each 0.6s electrocardiosignal, then it is assumed that be tentative QRS complex position, and remember
Its lower peak value is QRSPeak, if the point is met：Current RR>200ms, and QRSPeak>AvgQRSPeak/2, then it is assumed that should
Position is the position of QRS complex.
10. QRS complex detection means as claimed in claim 7, it is characterised in that the QRS complex detection means is also
Including screening module, the screening module is used for after obtaining QRS complex, and whether judge QRS complex is interference signal section
And reject the QRS complex for being judged as interference signal section.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710067025.7A CN106963368A (en) | 2017-02-07 | 2017-02-07 | QRS complex detection method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710067025.7A CN106963368A (en) | 2017-02-07 | 2017-02-07 | QRS complex detection method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106963368A true CN106963368A (en) | 2017-07-21 |
Family
ID=59335218
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710067025.7A Pending CN106963368A (en) | 2017-02-07 | 2017-02-07 | QRS complex detection method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106963368A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114027853A (en) * | 2021-12-16 | 2022-02-11 | 安徽心之声医疗科技有限公司 | QRS complex detection method, device, medium and equipment based on feature template matching |
| CN116407132A (en) * | 2023-04-17 | 2023-07-11 | 曦成半导体技术(上海)有限公司 | R wave positioning method and Q, S, T wave positioning method based on R wave |
| TWI810447B (en) * | 2020-03-23 | 2023-08-01 | 大陸商疆域康健創新醫療科技成都有限公司 | Method for detecting wave complex, ecg detection device and readable storage medium |
-
2017
- 2017-02-07 CN CN201710067025.7A patent/CN106963368A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI810447B (en) * | 2020-03-23 | 2023-08-01 | 大陸商疆域康健創新醫療科技成都有限公司 | Method for detecting wave complex, ecg detection device and readable storage medium |
| CN114027853A (en) * | 2021-12-16 | 2022-02-11 | 安徽心之声医疗科技有限公司 | QRS complex detection method, device, medium and equipment based on feature template matching |
| CN114027853B (en) * | 2021-12-16 | 2022-09-27 | 安徽心之声医疗科技有限公司 | QRS complex detection method, device, medium and equipment based on feature template matching |
| CN116407132A (en) * | 2023-04-17 | 2023-07-11 | 曦成半导体技术(上海)有限公司 | R wave positioning method and Q, S, T wave positioning method based on R wave |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gupta et al. | QRS complex detection using STFT, chaos analysis, and PCA in standard and real-time ECG databases | |
| Yang et al. | Using PPG signals and wearable devices for atrial fibrillation screening | |
| KR102008196B1 (en) | Apparatus for generating of Potassium level determination model using ECG and method thereof | |
| Orhan | Real-time CHF detection from ECG signals using a novel discretization method | |
| CN109009073A (en) | Atrial fibrillation detection device and storage medium | |
| CN114052744B (en) | Electrocardiosignal classification method based on impulse neural network | |
| CN104840186A (en) | Evaluation method of autonomic nervous function of patient suffering from CHF (Congestive Heart-Failure) | |
| CN103970975A (en) | Electrocardio data processing method and electrocardio data processing system | |
| CN105147252A (en) | Heart disease recognition and assessment method | |
| CN106963368A (en) | QRS complex detection method and device | |
| CN109044348A (en) | atrial fibrillation detection device and storage medium | |
| CN110226928A (en) | Atrial fibrillation is early with room and atrial fibrillation is with the poor recognition methods passed and device | |
| CN106874872A (en) | Industrial frequency noise filtering device and method | |
| RU2016146176A (en) | METHOD AND SYSTEM OF NON-INVASIVE SCREENING ASSESSMENT OF PHYSIOLOGICAL PARAMETERS AND PATHOLOGIES | |
| CN107530015A (en) | A kind of vital sign analysis method and system | |
| US20230153391A1 (en) | Electrocardiogram data classification method, computer device, and storage medium | |
| CN109948396A (en) | A kind of beat classification method, beat classification device and electronic equipment | |
| Kumari et al. | Performance evaluation of neural networks and adaptive neuro fuzzy inference system for classification of Cardiac Arrhythmia | |
| CN110327032A (en) | It is a kind of singly to lead the accurate recognizer of electrocardiosignal PQRST wave joint | |
| Mahesh et al. | ECG arrhythmia classification based on logistic model tree | |
| Yang et al. | ECG events detection and classification using wavelet and neural networks | |
| CN109674474B (en) | Sleep apnea recognition method, device and computer readable medium | |
| Xiong et al. | Enhancing the detection of atrial fibrillation from wearable sensors with neural style transfer and convolutional recurrent networks | |
| CN111839494A (en) | Heart rate monitoring method and system | |
| KR20140097039A (en) | Method and apparatus for classifying cardiac arrhythmia using an auto associative neural network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170721 |
|
| WD01 | Invention patent application deemed withdrawn after publication |