CN107273827A - A kind of R wave of electrocardiosignal detection method and device - Google Patents
A kind of R wave of electrocardiosignal detection method and device Download PDFInfo
- Publication number
- CN107273827A CN107273827A CN201710398845.4A CN201710398845A CN107273827A CN 107273827 A CN107273827 A CN 107273827A CN 201710398845 A CN201710398845 A CN 201710398845A CN 107273827 A CN107273827 A CN 107273827A
- Authority
- CN
- China
- Prior art keywords
- ripples
- order difference
- electrocardiosignal
- measured
- targeted customer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/14—Classification; Matching by matching peak patterns
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The present invention provides a kind of R wave of electrocardiosignal detection method and device.This method includes：The electrocardiosignal of the targeted customer of collection is pre-processed；First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the smooth signal of first-order difference of the electrocardiosignal is obtained；For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, the minimum value of the maximum for the smooth signal of first-order difference for being more than default rising edge differential threshold and the smooth signal of first-order difference less than default trailing edge differential threshold is obtained；The width of the R ripples to be measured of the targeted customer is obtained according to the maximum of the smooth signal of the first-order difference and the minimum value of the smooth signal of the first-order difference；The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer.The present invention carries out first-order difference extreme value dynamic threshold to R ripples bilateral and judged, finds the R crests position that first-order difference is zero, quickly and accurately realizes the R ripples detection of electrocardiosignal.
Description
Technical field
The present invention relates to technical field of data processing, and in particular to a kind of R wave of electrocardiosignal detection method and device.
Background technology
Containing abundant human body in daily lifes such as tranquillization, motions in electrocardiosignal (Electrocardiograph, ECG)
Physiological data in activity.Electrocardiosignal QRS complex (Q, R, S represent three home positions of ecg wave form respectively) detection is in motion
Occupy critical role in many life sign monitor systems of type continuous non-invasive.
The noise being ordinarily incorporated into electrocardiosignal has baseline drift, direct current offset, industrial frequency noise, ECG electrode piece and human body
The impedance capacitive reactance change that contact is produced, and the motor electrical noise that limb motion is produced.Wherein, motor electrical noise is with human body
Motion mode and intensity and change, band ratio baseline drift is slightly higher, in 1~5Hz or so, is maximum on ECG signal detection influence
Noise.
The R wave detecting methods of existing electrocardiosignal are by the identification to each characteristic point of ecg wave form.Existing feature
Point detecting method mainly has difference threshold algorithm, template matching method, Wavelet Transform and neural network etc..Difference threshold algorithm algorithm
Simply, processing speed is fast, it is easy to Project Realization, but often missing inspection false retrieval R ripples.Template matching method principle is simple, but to high frequency
Noise and the removal effect of baseline drift are poor.Wavelet Transform has good time frequency localization characteristic, and detection is accurate but right
The description of primary signal different scale needs substantial amounts of memory space, calculates time-consuming, adds hardware cost, is not easy to engineering real
It is existing.Neural network can be realized differentiates effect well, but the training time is longer, actually more difficult application.
The content of the invention
The embodiment of the present invention provides a kind of R wave of electrocardiosignal detection method and device, for how quickly and accurately to solve
Realize the problem of R ripples of electrocardiosignal are detected.
The embodiments of the invention provide a kind of R wave of electrocardiosignal detection method, including：
The electrocardiosignal of the targeted customer of collection is pre-processed；
First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the electrocardiosignal is obtained
The smooth signal of first-order difference；
For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, obtain poor more than default rising edge
The maximum of the smooth signal of first-order difference of point threshold value and the smooth signal of first-order difference less than default trailing edge differential threshold
Minimum value；
Institute is obtained according to the maximum of the smooth signal of the first-order difference and the minimum value of the smooth signal of the first-order difference
State the width of the R ripples to be measured of targeted customer；
The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer.
Alternatively, methods described also includes：
The default trailing edge differential threshold is updated according to the maximum of the smooth signal of the first-order difference.
Alternatively, default trailing edge differential threshold is updated according to equation below：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；
Wherein, spk=max (spk, max [y (i)])；Thr represents constant threshold value coefficient,
0≤thr≤1；Max [y (i)] represents the maximum of the smooth signal of first-order difference；A and b are constant, 0≤a≤1,0
≤b≤1。
Alternatively, the crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer, including：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than
Predetermined threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
The embodiment of the present invention provides a kind of R wave of electrocardiosignal detection means, including：
Pretreatment unit, the electrocardiosignal for the targeted customer to collection is pre-processed；
The smooth signal acquiring unit of first-order difference, for carrying out first-order difference processing peace to pretreated electrocardiosignal
Equal smooth treatment, obtains the smooth signal of first-order difference of the electrocardiosignal；
Max min acquiring unit, the first-order difference for the electrocardiosignal for each sampling period is smooth
Signal, obtains the maximum for the smooth signal of first-order difference for being more than default rising edge differential threshold and less than default trailing edge difference
The minimum value of the smooth signal of first-order difference of threshold value；
R wave widths acquiring unit to be measured, for the maximum according to the smooth signal of the first-order difference and a jump
Divide the width of the R ripples to be measured of the minimum value acquisition targeted customer of smooth signal；
Crest location acquiring unit, the width for the R ripples to be measured according to the targeted customer determines the R ripples to be measured
Crest location.
Alternatively, in addition to：
Trailing edge differential threshold updating block, for the maximum according to the smooth signal of the first-order difference to described default
Trailing edge differential threshold is updated.
Alternatively, the trailing edge differential threshold updating block is further used for：
Default trailing edge differential threshold is updated according to equation below：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；
Wherein, spk=max (spk, max [y (i)])；Thr represents constant threshold value coefficient,
0≤thr≤1；Max [y (i)] represents the maximum of the smooth signal of first-order difference；A and b are constant, 0≤a≤1,0
≤b≤1。
Alternatively, the crest location acquiring unit is further used for：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than
Predetermined threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
The embodiment of the present invention provides a kind of electronic equipment, it is characterised in that including：Processor, memory and bus；Its
In,
Processor and memory complete mutual communication by bus；
Processor is used to call the programmed instruction in memory, to perform above-mentioned R wave of electrocardiosignal detection method.
The embodiment of the present invention provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage
Medium storing computer is instructed, and the computer instruction makes the computer perform above-mentioned R wave of electrocardiosignal detection method.
R wave of electrocardiosignal detection method and device provided in an embodiment of the present invention, believe the electrocardio of the targeted customer of collection
Number pre-processed；First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the electrocardio is obtained
The smooth signal of first-order difference of signal；For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, obtain
One more than the maximum of the smooth signal of first-order difference of default rising edge differential threshold and less than default trailing edge differential threshold
The minimum value of the smooth signal of order difference；According to the maximum of the smooth signal of the first-order difference and the smooth signal of the first-order difference
Minimum value obtain the targeted customer R ripples to be measured width；
The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer.The embodiment of the present invention
First-order difference extreme value dynamic threshold is carried out to R ripples bilateral to judge, finds the R crests position that first-order difference is zero, it is quickly and accurately real
The R ripples detection of existing electrocardiosignal.
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 do simply to introduce, it should be apparent that, drawings in the following description are the present invention
Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can also basis
These accompanying drawings obtain other accompanying drawings.
Fig. 1 is the schematic flow sheet of the R wave of electrocardiosignal detection method of one embodiment of the invention；
Fig. 2 is the structural representation of many life sign monitor systems of sports type continuous non-invasive of one embodiment of the invention；
Fig. 3 is the structural representation of the R wave of electrocardiosignal detection means of one embodiment of the invention；
Fig. 4 is the entity structure schematic diagram of the electronic equipment of one embodiment of the invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, clear, complete description is carried out to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is
A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Fig. 1 is the schematic flow sheet of the R wave of electrocardiosignal detection method of one embodiment of the invention.As shown in figure 1, the reality
Applying the method for example includes：
S11：The electrocardiosignal of the targeted customer of collection is pre-processed；
It should be noted that being put down to the preprocessing process of the electrocardiosignal of targeted customer including removing baseline drift, M point
Equal smooth treatment, pretreated electrocardiosignal is x (i), and i is sample data sequence number.
In actual applications, can use has limit for length's unit impulse response (Finite Impulse Response, FIR) filter
Ripple device removes baseline drift, it would however also be possible to employ other modes remove baseline drift, the invention is not limited in this regard.
It should be noted that the number M of the point of average smooth treatment, need to integrate R wave widths TR(the width T of R ripplesRFor R ripples
Time interval between rising edge and trailing edge) and sample rate fsTo determine.Usual M point at least needs to cover TR, i.e.,
M≥fs·TR
For example, fs=500Hz, TR=0.030 second, then M >=15.
S12：First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the electrocardio letter is obtained
Number the smooth signal of first-order difference；
First-order difference is asked to x (i)：
D (i)=x (i+1)-x (i)
N number of point is done to d (i) again average smooth, obtain y (i).
Averagely smooth data point number N is to first-order difference d (i), human limb trunk motion frequency f, R ripple need to be integrated
Width TRWith sample rate fsIt is determined that.Usual N number of point at least needs to cover limb motion cycle 1/f, i.e.,
N≥fs/f
For example：fs=500Hz, f=10Hz, then N >=50.
S13：For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, obtain and be more than default rising
Smoothly believe along the maximum of the smooth signal of the first-order difference of differential threshold and less than the first-order difference of default trailing edge differential threshold
Number minimum value；
In actual applications, in order to improve the accuracy that R ripples are detected, rising edge differential threshold (spkthr) is preset and pre-
If trailing edge differential threshold (- spkthr), takes 10≤spk≤100 and 0.1≤thr≤1.0.
S14：Obtained according to the minimum value of the maximum of the smooth signal of the first-order difference and the smooth signal of the first-order difference
Take the width of the R ripples to be measured of the targeted customer；
S15：The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer；
It should be noted that being sought between the maximum of rising edge and the minimum value of trailing edge of the smooth signal of first-order difference
Look for the position that first-order difference is 0, the crest location of R ripples as to be measured.
R wave of electrocardiosignal detection method provided in an embodiment of the present invention, first-order difference extreme value dynamic threshold is carried out to R ripples bilateral
Value judges, finds the R crests position that first-order difference is zero, quickly and accurately realizes the R ripples detection of electrocardiosignal.
In a kind of optional embodiment of the embodiment of the present invention, similar with the method in Fig. 1, methods described is also wrapped
Include：
The default trailing edge differential threshold is updated according to the maximum of the smooth signal of the first-order difference.
Specifically, default trailing edge differential threshold is updated according to equation below：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；
Wherein, spk=max (spk, max [y (i)])；Thr represents constant threshold value coefficient,
0≤thr≤1；Ma [y (i)] represents the maximum of the smooth signal of first-order difference；A and b are constants, 0≤a≤1,0≤
b≤1。
Further, the crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer, is wrapped
Include：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
It should be noted that the quiescent condition R wave widths of targeted customer take targeted customer to measure acquisition under quiescent condition
Multiple R ripples rising edge and trailing edge between time interval average value；Quiescent condition R wave widths are that have with measurement object
The dynamic value of pass.
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than
Predetermined threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
It will be appreciated that the embodiment of the present invention is examined after the judgement of first-order difference extreme value dynamic threshold is carried out to R ripples bilateral
The width of R ripples to be measured meets after the quiescent condition R wave widths of targeted customer, finds the crest location for the R ripples that first-order difference is zero.
Specifically, in sample frequency fsIn length, the R ripple rising edges for being more than given default rising edge differential threshold (spkthr) are found
Y (i) maximum：y(is)=max [y (i)].If it is found, setting new threshold value
Spk=aspk+bmax [y (i)]
Wherein a and b are undetermined constant, 0≤a≤1,0≤b≤1.When taking a=0.4, b=0.6, good R ripples can be obtained
Detection results.
Be further continued for find y (i) be less than trailing edge differential threshold (- spkthr) R ripple trailing edges y (i) it is minimum
Value：y(if)=min [y (i)].
The width T of R ripplesRFor the time interval between R ripples rising edge and trailing edge, i.e.,
TR=(if-is)fs
Wherein, ifRepresent the trailing edge of R ripples, isRepresent the rising edge of R ripples.
Under normal circumstances, the width T of R ripplesR≈ 30 milliseconds=0.030 second.
If have found R ripple trailing edges, the time interval (i between the R ripple bilateral difference extreme values for calculating and obtaining is examinedf-is)
fs.If (if-is)fsIn (TR± e%) in the range of, then it is assumed that R ripples are in isAnd ifBetween.Here e is default error constant.
Afterwards, then in isAnd ifBetween find first-order difference y (iR)=0, so that it is determined that R ripples position iR。
In actual applications, when continually looking for R ripples, using refractory period condition simplified operation, i.e., in refractory period without R ripples
Detection, can save the time, improve the efficiency of R ripples detection positioning.
When the musculus pectoralis abdominis that human body occurs 10 times per second or so is acutely shaken, the motor electrical noise fluctuation of generation, width
Similar to the width of R ripples, steepness absolute value is less than R ripples rising edge and the absolute value of trailing edge steepness.Although this kind of noise is more
It is rare, but it has the continuous repeated of certain time length, the present invention is proposed：It is by additional frequency
Wave filter remove.
The musculus pectoralis abdominis of high frequency is trembled again, and signal smooth treatment can be removed.
In actual applications, the R wave of electrocardiosignal detection method of the embodiment of the present invention can be applied to motion as shown in Figure 2
In many life sign monitor systems of type continuous non-invasive.The many life sign monitor systems of sports type continuous non-invasive, are after applied to hand
After the various invasive and noninvasive life sign monitor system that art room, hospital, family use when being monitored static person and rest,
It is of future generation to carry out health and the core technology of physiological parameter monitoring in long-range tranquillization and motion.
As shown in Fig. 2 many life sign monitor systems of sports type continuous non-invasive include monitoring host machine part 21 and wearing portion
Divide 22；Monitor host machine part 21 include the peripheral data such as system host, keyboard, mouse, display screen, printer, network, figure and
The input such as energy, output and communication apparatus；Wearing portion 22 includes electrocardiogram acquisition electrode 221, photoplethysmographic and gathered
Device 222 and monitor 223.Monitor 223 gathers electrocardiogram acquisition electrode 221 and photoplethysmographic collector 222
Electrocardiosignal and pulse signal are sent to monitoring host machine part 21；Monitor the method pair according to embodiments of the present invention of host machine part 21
Electrocardiosignal is handled；And many vital signs of the sports type continuous non-invasive obtained after processing cardioelectric signals are monitored into supplemental characteristic
Send to monitor 223；The many vital sign prisons of sports type continuous non-invasive after the processing cardioelectric signals that 223 pairs of monitor is received
Survey supplemental characteristic and carry out numeral, figure shows.In actual applications, monitoring host machine part 21 and monitor 223 can be gone here and there by USB
Mouth, bluetooth, netting twine and WiFi are communicated, can also shared processor, memory and bus.
The software for the R wave of electrocardiosignal detection method that the present invention is provided is applied, can also be loaded into and meet hardware configuration
It is required that PC computers, the webserver, notebook computer, iPad, tablet personal computer, single-chip microcomputer, the intelligence of iOS and Android operation system
On energy mobile phone, above-mentioned smart machine is set to turn into many life sign monitor system host machine parts of sports type continuous non-invasive.
Tables 1 and 2 shows the testing result of the R wave of electrocardiosignal detection method using the embodiment of the present invention, such as He of table 1
Shown in table 2, the degree of accuracy that the method for the embodiment of the present invention is detected to R wave of electrocardiosignal is very high.
The testing result of the true man's R wave of electrocardiosignal of table 1
| Manner | R ripples sum | Missing inspection number | Flase drop number | Mistake sum | Accuracy rate |
| It is quiet | 2645 | 0 | 0 | 0 | 100% |
| It is careful | 3489 | 0 | 0 | 0 | 100% |
| Slow upper extremity exercise | 2677 | 0 | 0 | 0 | 100% |
| Do gymnastics | 2311 | 0 | 0 | 0 | 100% |
| Run | 3123 | 0 | 1 | 1 | 99.97% |
| Fast upper extremity exercise | 2155 | 0 | 1 | 1 | 99.95% |
| Fast lower extremity movement | 3688 | 0 | 1 | 1 | 99.97% |
| It is total | 20088 | 0 | 3 | 3 | 99.99% |
The testing result of the EGC analog device R wave of electrocardiosignal of table 2
Fig. 3 is the structural representation of the R wave of electrocardiosignal detection means of one embodiment of the invention.As shown in figure 3, the reality
Apply example device include pretreatment unit 31, the smooth signal acquiring unit 32 of first-order difference, max min acquiring unit 33,
R wave widths acquiring unit 34 and crest location acquiring unit 35 to be measured, specifically：
Pretreatment unit 31, the electrocardiosignal for the targeted customer to collection is pre-processed；
The smooth signal acquiring unit 32 of first-order difference, for pretreated electrocardiosignal carry out first-order difference processing and
Average smooth treatment, obtains the smooth signal of first-order difference of the electrocardiosignal；
Max min acquiring unit 33, the first-order difference for the electrocardiosignal for each sampling period is put down
Sliding signal, obtains the maximum of the smooth signal of first-order difference for being more than default rising edge differential threshold and poor less than default trailing edge
Divide the minimum value of the smooth signal of first-order difference of threshold value；
R wave widths acquiring unit 34 to be measured, for the maximum and the single order according to the smooth signal of the first-order difference
The minimum value of the smooth signal of difference obtains the width of the R ripples to be measured of the targeted customer；
Crest location acquiring unit 35, the width for the R ripples to be measured according to the targeted customer determines the R ripples to be measured
Crest location.
R wave of electrocardiosignal detection means provided in an embodiment of the present invention, first-order difference extreme value dynamic threshold is carried out to R ripples bilateral
Value judges, finds the R crests position that first-order difference is zero, quickly and accurately realizes the R ripples detection of electrocardiosignal.
In a kind of optional embodiment of the embodiment of the present invention, the device also includes：
Trailing edge differential threshold updating block, for the maximum according to the smooth signal of the first-order difference to described default
Trailing edge differential threshold is updated.
Trailing edge differential threshold updating block is further used for：
Default trailing edge differential threshold is updated according to equation below：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；Wherein, spk=max (spk,
max[y(i)])；Thr represents constant threshold value coefficient, 0≤thr≤1；Max [y (i)] represents the smooth signal of first-order difference most
Big value；A and b are constant, 0≤a≤1,0≤b≤1.
Crest location acquiring unit 35 is further used for：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than
Predetermined threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
The R wave of electrocardiosignal detection means of the embodiment of the present invention can be used for perform above method embodiment, its principle and
Technique effect is similar, and here is omitted.
Fig. 4 is the structural representation of the electronic equipment of one embodiment of the invention.
Reference picture 4, electronic equipment includes：Processor (processor) 41, memory (memory) 42 and bus (bus)
43；Wherein,
Processor 41 and memory 42 complete mutual communication by bus 43；
Processor 41 is used to call the programmed instruction in memory 42, to perform the heart that above-mentioned each method embodiment is provided
Electric signal R wave detecting methods.
In addition, the logical order in above-mentioned memory 42 is realized and is used as independent in the form by SFU software functional unit
Production marketing or in use, can be stored in a computer read/write memory medium.Understood based on such, the present invention
The part that is substantially contributed in other words to prior art of technical scheme or the part of the technical scheme can be with software
The form of product is embodied, and the computer software product is stored in a storage medium, including some instructions are to cause
One computer equipment (can be personal computer, notebook computer, server, tablet personal computer, smart mobile phone, single-chip microcomputer or
Person's network equipment etc.) perform all or part of step of each of the invention embodiment methods described.And foregoing storage medium bag
Include：USB flash disk, hard 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 be with the medium of store program codes.
The present embodiment provides a kind of computer program product, and the computer program product includes being stored in non-transient calculating
Computer program on machine readable storage medium storing program for executing, the computer program includes programmed instruction, when described program instruction is calculated
When machine is performed, computer is able to carry out the R wave of electrocardiosignal detection method that above-mentioned each method embodiment is provided.
The present embodiment provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage medium storing program for executing
Computer instruction is stored, the computer instruction makes the computer perform the electrocardiosignal that above-mentioned each method embodiment is provided
R wave detecting methods.
R wave of electrocardiosignal detection method and device provided in an embodiment of the present invention, believe the electrocardio of the targeted customer of collection
Number pre-processed；First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the electrocardio is obtained
The smooth signal of first-order difference of signal；For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, obtain
One more than the maximum of the smooth signal of first-order difference of default rising edge differential threshold and less than default trailing edge differential threshold
The minimum value of the smooth signal of order difference；According to the maximum of the smooth signal of the first-order difference and the smooth signal of the first-order difference
Minimum value obtain the targeted customer R ripples to be measured width；
The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer.The embodiment of the present invention
First-order difference extreme value dynamic threshold is carried out to R ripples bilateral to judge, finds the R crests position that first-order difference is zero, it is quickly and accurately real
The R ripples detection of existing electrocardiosignal.
It should be understood by those skilled in the art that, embodiments of the invention can be method, system or computer program production
Product.Therefore, the present invention can be using the implementation in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
The form of example.Moreover, of the invention can use can in one or more computers for wherein including computer usable program code
With the computer program product implemented on storage medium (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Form.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product
Figure and/or block diagram are described.Being interpreted as can be by each in computer program instructions implementation process figure and/or block diagram
Flow and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer journeys can be provided
Sequence instruction to all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices processor with
A machine is produced, the instruction for becoming the computing device by computer or other programmable data processing devices is produced
For realizing the function of being specified in one flow of flow chart or multiple flows and/or one square frame of block diagram or multiple square frames
Device.
It should be noted that term " comprising ", "comprising" or its any other variant are intended to the bag of nonexcludability
Contain, so that process, method, article or equipment including a series of key elements are not only including those key elements, but also including
Other key elements being not expressly set out, or also include for this process, method, article or the intrinsic key element of equipment.
In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that including the key element
Process, method, article or equipment in also there is other identical element.
In the specification of the present invention, numerous specific details are set forth.Although it is understood that, embodiments of the invention can
To be put into practice in the case of these no details.In some instances, known method, structure and skill is not been shown in detail
Art, so as not to obscure the understanding of this description.Similarly, it will be appreciated that disclose in order to simplify the present invention and helps to understand respectively
One or more of individual inventive aspect, above in the description of the exemplary embodiment of the present invention, each of the invention is special
Levy and be grouped together into sometimes in single embodiment, figure or descriptions thereof.However, should not be by the method solution of the disclosure
It is interpreted into the following intention of reflection：I.e. the present invention for required protection requirement is than the feature that is expressly recited in each claim more
Many features.More precisely, as the following claims reflect, inventive aspect is to be less than single reality disclosed above
Apply all features of example.Therefore, it then follows thus claims of embodiment are expressly incorporated in the embodiment,
Wherein each claim is in itself as the separate embodiments of the present invention.
Above example is merely to illustrate technical scheme, rather than its limitations；Although with reference to the foregoing embodiments
The present invention is described in detail, it will be understood by those within the art that：It still can be to foregoing each implementation
Technical scheme described in example is modified, or carries out equivalent substitution to which part technical characteristic；And these are changed or replaced
Change, the essence of appropriate technical solution is departed from the spirit and scope of various embodiments of the present invention technical scheme.
Claims (10)
1. a kind of R wave of electrocardiosignal detection method, it is characterised in that including：
The electrocardiosignal of the targeted customer of collection is pre-processed；
First-order difference processing and average smooth treatment are carried out to pretreated electrocardiosignal, the single order of the electrocardiosignal is obtained
The smooth signal of difference；
For the smooth signal of first-order difference of the electrocardiosignal in each sampling period, obtain and be more than default rising edge difference threshold
The minimum of the maximum of the smooth signal of first-order difference of value and the smooth signal of first-order difference less than default trailing edge differential threshold
Value；
The mesh is obtained according to the maximum of the smooth signal of the first-order difference and the minimum value of the smooth signal of the first-order difference
Mark the width of the R ripples to be measured of user；
The crest location of the R ripples to be measured is determined according to the width of the R ripples to be measured of the targeted customer.
2. according to the method described in claim 1, it is characterised in that methods described also includes：
The default trailing edge differential threshold is updated according to the maximum of the smooth signal of the first-order difference.
3. method according to claim 2, it is characterised in that carried out according to equation below to default trailing edge differential threshold
Update：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；
Wherein, spk=max (spk, max [y (i)])；Thr represents constant threshold value coefficient, 0≤thr≤1；Max [y (i)] table
Show the maximum of the smooth signal of first-order difference；A and b are constant, 0≤a≤1,0≤b≤1.
4. according to the method described in claim 1, it is characterised in that determined according to the width of the R ripples to be measured of the targeted customer
The crest location of the R ripples to be measured, including：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than default
Threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
5. a kind of R wave of electrocardiosignal detection means, it is characterised in that including：
Pretreatment unit, the electrocardiosignal for the targeted customer to collection is pre-processed；
The smooth signal acquiring unit of first-order difference, for carrying out first-order difference processing and average light to pretreated electrocardiosignal
Sliding processing, obtains the smooth signal of first-order difference of the electrocardiosignal；
Max min acquiring unit, the first-order difference for the electrocardiosignal for each sampling period is smoothly believed
Number, obtain the maximum for the smooth signal of first-order difference for being more than default rising edge differential threshold and less than default trailing edge difference threshold
The minimum value of the smooth signal of first-order difference of value；
R wave widths acquiring unit to be measured, puts down for the maximum and the first-order difference according to the smooth signal of the first-order difference
The minimum value of sliding signal obtains the width of the R ripples to be measured of the targeted customer；
Crest location acquiring unit, the width for the R ripples to be measured according to the targeted customer determines the crest of the R ripples to be measured
Position.
6. device according to claim 5, it is characterised in that also include：
Trailing edge differential threshold updating block, for the maximum according to the smooth signal of the first-order difference to the default decline
It is updated along differential threshold.
7. device according to claim 6, it is characterised in that the trailing edge differential threshold updating block is further used
In：
Default trailing edge differential threshold is updated according to equation below：
Rd=-spkthr
Wherein, RdRepresent trailing edge differential threshold；Spk=aspk+bmax [y (i)]；Wherein, spk=max (spk, max [y
(i)])；Thr represents constant threshold value coefficient, 0≤thr≤1；Max [y (i)] represents the maximum of the smooth signal of first-order difference；a
It is constant, 0≤a≤1,0≤b≤1 with b.
8. device according to claim 5, it is characterised in that the crest location acquiring unit is further used for：
The width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer are compared；
If the difference of the width of the R ripples to be measured of the targeted customer and the quiescent condition R wave widths of the targeted customer is less than default
Threshold value, then determine the crest location of the R ripples to be measured according to the width of the R ripples to be measured of the targeted customer.
9. a kind of electronic equipment, it is characterised in that including：Processor, memory and bus；Wherein,
Processor and memory complete mutual communication by bus；
Processor is used to call the programmed instruction in memory, and the electrocardiosignal R described in 1~4 any one is required with perform claim
Wave detecting method.
10. a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage medium storing program for executing storage computer refers to
Order, the computer instruction makes the computer perform claim require the R wave of electrocardiosignal detection method described in 1~4 any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710398845.4A CN107273827B (en) | 2017-05-31 | 2017-05-31 | Electrocardiosignal R wave detection method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710398845.4A CN107273827B (en) | 2017-05-31 | 2017-05-31 | Electrocardiosignal R wave detection method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107273827A true CN107273827A (en) | 2017-10-20 |
| CN107273827B CN107273827B (en) | 2020-01-07 |
Family
ID=60064940
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710398845.4A Active CN107273827B (en) | 2017-05-31 | 2017-05-31 | Electrocardiosignal R wave detection method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107273827B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107997759A (en) * | 2017-10-27 | 2018-05-08 | 北京康博众联电子科技有限公司 | The detection method of QRS complex, storage medium and computer equipment in electrocardiosignal |
| CN111046836A (en) * | 2019-12-24 | 2020-04-21 | 杭州电力设备制造有限公司 | Method, system, equipment and storage medium for filtering, denoising and analyzing partial discharge signal |
| CN111345812A (en) * | 2018-12-24 | 2020-06-30 | Zoll医疗公司 | R wave detection method and device, and computer readable storage medium |
| CN111938640A (en) * | 2020-08-05 | 2020-11-17 | 深圳大学 | Method and device for positioning point C of cardiac impedance differential signal and storage medium |
| CN114521901A (en) * | 2021-12-28 | 2022-05-24 | 宁波慈溪生物医学工程研究所 | Electrocardio feature extraction method and device and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101828917A (en) * | 2010-05-07 | 2010-09-15 | 深圳大学 | Method and system for extracting electrocardiosignal characteristic |
| US20120095358A1 (en) * | 2010-10-14 | 2012-04-19 | Denso Corporation | Living body state monitor apparatus |
| CN103549950A (en) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring |
| CN103654770A (en) * | 2013-12-03 | 2014-03-26 | 上海理工大学 | Method and device for real-time detection of moving electrocardiosignal QRS (magnetic resonance angiography) waves |
| WO2016035701A1 (en) * | 2014-09-04 | 2016-03-10 | 日本電信電話株式会社 | Heartbeat detection method and heartbeat detection device |
-
2017
- 2017-05-31 CN CN201710398845.4A patent/CN107273827B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101828917A (en) * | 2010-05-07 | 2010-09-15 | 深圳大学 | Method and system for extracting electrocardiosignal characteristic |
| US20120095358A1 (en) * | 2010-10-14 | 2012-04-19 | Denso Corporation | Living body state monitor apparatus |
| CN103549950A (en) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring |
| CN103654770A (en) * | 2013-12-03 | 2014-03-26 | 上海理工大学 | Method and device for real-time detection of moving electrocardiosignal QRS (magnetic resonance angiography) waves |
| WO2016035701A1 (en) * | 2014-09-04 | 2016-03-10 | 日本電信電話株式会社 | Heartbeat detection method and heartbeat detection device |
Non-Patent Citations (1)
| Title |
|---|
| 王蔷薇: ""面向智能终端的心电信号处理方法研究及实现"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107997759A (en) * | 2017-10-27 | 2018-05-08 | 北京康博众联电子科技有限公司 | The detection method of QRS complex, storage medium and computer equipment in electrocardiosignal |
| CN111345812A (en) * | 2018-12-24 | 2020-06-30 | Zoll医疗公司 | R wave detection method and device, and computer readable storage medium |
| CN111046836A (en) * | 2019-12-24 | 2020-04-21 | 杭州电力设备制造有限公司 | Method, system, equipment and storage medium for filtering, denoising and analyzing partial discharge signal |
| CN111938640A (en) * | 2020-08-05 | 2020-11-17 | 深圳大学 | Method and device for positioning point C of cardiac impedance differential signal and storage medium |
| CN111938640B (en) * | 2020-08-05 | 2023-02-07 | 深圳大学 | Method and device for positioning point C of cardiac impedance differential signal and storage medium |
| CN114521901A (en) * | 2021-12-28 | 2022-05-24 | 宁波慈溪生物医学工程研究所 | Electrocardio feature extraction method and device and storage medium |
| CN114521901B (en) * | 2021-12-28 | 2024-03-19 | 宁波慈溪生物医学工程研究所 | Electrocardiogram feature extraction method, device and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107273827B (en) | 2020-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107273827A (en) | A kind of R wave of electrocardiosignal detection method and device | |
| US10866637B2 (en) | Gesture classification apparatus and method using EMG signal | |
| Bashar et al. | Atrial fibrillation detection from wrist photoplethysmography signals using smartwatches | |
| Machado et al. | Human activity data discovery from triaxial accelerometer sensor: Non-supervised learning sensitivity to feature extraction parametrization | |
| CN105813558B (en) | System and method for determining sleep stages based on sleep cycle | |
| CN105451652B (en) | For determining the processing unit and processing method of the breath signal of object | |
| CN105578960B (en) | For handling processing unit, the processing method and system of physiological signal | |
| US11412972B2 (en) | Detection of atrial fibrillation | |
| Benson et al. | Classifying running speed conditions using a single wearable sensor: Optimal segmentation and feature extraction methods | |
| CN107595243A (en) | A kind of illness appraisal procedure and terminal device | |
| CN112165897A (en) | Photoplethysmogram data analysis and presentation | |
| CN107714024A (en) | For monitoring the method, system and device of cardiomotility | |
| CN105792740A (en) | Detection and calculation of heart rate recovery in non-clinical settings | |
| CN105637502B (en) | For determining the method and apparatus for being smoothed data point in data point stream | |
| CN105844075A (en) | Method for low-power-consumption, robust estimation of cardiovascular periodicity, contour analysis, and heart rate | |
| CN104523281A (en) | Movement monitoring method and system and movement monitoring clothes | |
| Panahi et al. | Application of fractional Fourier transform in feature extraction from ELECTROCARDIOGRAM and GALVANIC SKIN RESPONSE for emotion recognition | |
| EP3165158B1 (en) | Signal feature extracting apparatus | |
| CN106886675A (en) | Method and system for monitoring continuous biomedicine signals | |
| Pander et al. | An automatic saccadic eye movement detection in an optokinetic nystagmus signal | |
| KR101977745B1 (en) | Method and apparatus for classifying cardiac arrhythmia | |
| US20120277816A1 (en) | Adjusting neighborhood widths of candidate heart beats according to previous heart beat statistics | |
| JP6842306B2 (en) | Biological signal processing method, its device, and computer-readable recording medium | |
| Khairuddin et al. | A general framework for improving electrocardiography monitoring system with machine learning | |
| Mishra et al. | Performance Evaluation of Various Window Techniques for Noise Cancellation from ECG Signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A method and device for detecting R-wave of electrocardiogram signals Effective date of registration: 20230717 Granted publication date: 20200107 Pledgee: Jiangsu Jiangyin Rural Commercial Bank Co.,Ltd. high tech Zone sub branch Pledgor: JIANGSU STANDARD MEDICAL TECHNOLOGY Co.,Ltd. Registration number: Y2023980048651 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |