CN103690156A - Heart rate acquiring method and electrocardiosignal processing method - Google Patents
Heart rate acquiring method and electrocardiosignal processing method Download PDFInfo
- Publication number
- CN103690156A CN103690156A CN201310598475.0A CN201310598475A CN103690156A CN 103690156 A CN103690156 A CN 103690156A CN 201310598475 A CN201310598475 A CN 201310598475A CN 103690156 A CN103690156 A CN 103690156A
- Authority
- CN
- China
- Prior art keywords
- heart rate
- ripple
- moment
- time
- value
- 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
Images
Abstract
The invention relates to a heart rate acquiring method and an electrocardiosignal processing method. The heart rate acquiring method comprises the following steps of recording moments when an R wave appears on an electrocardiosignal waveform; on the basis of the moment TRi when the R wave appears for the ith times and the moment TRi-1 when the R wave appears for the i-1th times, acquiring an instantaneous heart rate at the moment TRi. The electrocardiosignal processing method comprises the following steps of on the basis of an input electrocardiosignal, acquiring the heart rate at the moment TRi by utilizing the heart rate acquiring method; processing the heart rate at the moment TRi. According to the technical scheme adopted by the invention, the heart rate can be more accurately acquired and real variation of the heart rate is reflected.
Description
Technical field
Technical solution of the present invention relates to the information processing of electrocardiosignal, particularly the processing method of a kind of heart rate acquisition methods and electrocardiosignal.
Background technology
Electrocardiosignal, also be that ECG(ECG is the abbreviation of electrocardiogram) signal is that heart is in each cardiac cycle, by the excitement in succession of pacemaker, atrium, ventricle, be accompanied by bioelectric variation, by electrocardiograph, from body surface, draw the figure of the potential change of various ways.When individual heartbeat, electric signal streams is crossed heart, and it can be acquired by the appropriate device being connected on individual skin.Typical and complete electrocardiosignal comprises and is considered to the complete waveform that consists of P ripple, Q ripple, R ripple, S ripple and T ripple, and wherein, Q ripple, R ripple and S ripple are become formation " QRS wave group " sometimes jointly.Generally speaking, P ripple is produced by the depolarization in atrium, and QRS wave group is produced by the depolarization of ventricle, and T ripple has the polarization again of ventricular muscle to produce.The meaning that detects described electrocardiosignal is, can to diseases such as various arrhythmia, ventricular atrial hypertrophy, myocardial infarction, arrhythmia, myocardial ischemia, check by detected electrocardiosignal.
Conventionally, R ripple is the most significant ripple in whole electrocardiosignal waveform, has the advantages that waveform is precipitous, amplitude is large, width is narrow, be easy to identification.Therefore, individual heart rate can be detected and be determined by the R ripple number to ECG signal in the unit interval.
Heart rate refers to the frequency of heartbeat, can be used for weighing heart number of times of beating per minute.Detect the minor variations of heart rate for people's health monitoring, exercise intensity judgement,, there is very important effect the aspects such as evaluate application of the early warning of physical fatigue degree and energy expenditure.
A kind of heart rate acquisition methods of the prior art can be with reference to figure 1, by gathering the number of beating of heart in a minute, obtain heart rate value, comprise: input ecg wave form, at 1 minute clocking internal and the ecg wave form of input in to 1 minute detect R ripple, and calculate R ripple number according to the R ripple detecting, after 1 minutes arrives, the R ripple number detecting in 1 minute of usining is exported as heart rate value.The method is very accurately concerning the measurement of heart rate, and can accomplish zero error.But because this method adopted in the set time, detect the R ripple number collecting, can lose the difference between adjacent two R ripples, can't detect the minor variations in each heartbeat interval, cannot get instantaneous heart rate value.
Summary of the invention
The technical problem that technical solution of the present invention solves is how to obtain more exactly heart rate, and reflect the real change of heart rate.
Based on above-mentioned technical problem, technical solution of the present invention provides a kind of heart rate acquisition methods, comprising:
Record the moment that on electrocardiosignal waveform, R ripple occurs;
Moment T based on the i time R ripple appearance
rimoment T with the i-1 time R ripple appearance
ri-1, obtain T constantly
rion instantaneous heart rate.
Optionally, describedly record the moment that R ripple on electrocardiosignal waveform occurs and comprise:
Described electrocardiosignal waveform is sampled to obtain sampled point, and the corresponding described electrocardiosignal wave-shape amplitude of described sampled point is worth time dependent relation;
Detect R ripple, range value is greater than to the sampled point of amplitude threshold as the sampled point of described R ripple;
The sampled point of searching for described R ripple is to obtain time corresponding to maximum point in the sampled point of described R ripple.
Optionally, the sampled point of the described R ripple of described search comprised with time corresponding to maximum point obtaining in the sampled point of described R ripple:
Calculate the sampled point slope of described R ripple, the sampled point slope of described R ripple along described R ripple gradually from the occasion of becoming negative value;
Identification slope is from the sampled point of the R ripple that is negative value on the occasion of saltus step;
Gather the time corresponding to sampled point of the R ripple of identifying.
Optionally, described heart rate acquisition methods also comprises: upgrade described amplitude threshold, the described amplitude threshold of described renewal comprises:
Obtain the amplitude equalizing value of the maximum point of the several R ripples that are consecutively detected;
Obtained amplitude equalizing value is multiplied by threshold value coefficient, to obtain new amplitude threshold.
Optionally, the amplitude threshold upgrading comprises the amplitude threshold that upgraded and initial amplitude threshold, and described heart rate acquisition methods also comprises:
Obtain the amplitude peak value in the some sampled points that are consecutively detected;
Obtained amplitude peak value is multiplied by described threshold value coefficient, to obtain described initial amplitude threshold.
Optionally, based on following formula, obtain T constantly
rion instantaneous heart rate:
F=60/T; Or,
F=1/T；
Wherein, T is occurred the i time R ripple and the interval between the i-1 time R ripple, and F is T constantly
rion instantaneous heart rate.
Based on above-mentioned technical problem, technical solution of the present invention also provides a kind of processing method of electrocardiosignal, comprising:
Electrocardiosignal based on input, utilizes method as above to obtain T constantly
ritime heart rate;
Process described moment T
ritime heart rate.
Optionally, the processing method of described electrocardiosignal also comprises: as described moment T
ritime heart rate when being less than the first heart rate value or being greater than the second heart rate value, obtain T constantly
ritime electrocardiosignal waveform wave characteristic long-pending;
The described moment T of described processing
riheart rate comprise: more described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform is long-pending, to transmit described moment T
riheart rate.
Optionally, described in, obtain T constantly
ritime electrocardiosignal waveform long-pending comprising of wave characteristic:
Extract T constantly
ri-1to moment T
ribetween ecg wave form, described moment T
ri-1to moment T
ribetween ecg wave form be divided into N section wavy curve, wherein N is greater than 2 natural number;
Calculate respectively the slope of described N section wavy curve, described moment T
ritime electrocardiosignal waveform wave characteristic amass the absolute value into the slope product of described N section wavy curve.
Optionally, the long-pending scope of living in of the wave characteristic of described R ripple comprises the first fluctuation range and the second fluctuation range that is greater than described the second fluctuation threshold value between the first fluctuation threshold value and the second fluctuation threshold value, and described the first fluctuation threshold value is less than the second fluctuation threshold value.
Optionally, described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform is long-pending to process described moment T
riheart rate comprise:
At described moment T
riheart rate while being less than the first heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic long-pending in described the first fluctuation range, transmit described moment T
ritime heart rate, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly greater than described the second fluctuation threshold value, abandon described moment T
ritime heart rate;
At described moment T
riheart rate while being greater than the second heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic long-pending in described the first fluctuation range or the second fluctuation range, transmit described moment T
ritime heart rate, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly less than described the first fluctuation threshold value, abandon described moment T
ritime heart rate.
The beneficial effect of technical solution of the present invention at least comprises:
A kind of heart rate acquisition methods is provided, and the heart real time image data (being electrocardiosignal waveform) that the method can be based on input, by obtaining the interval time between R ripple, extracts heart rate value exactly; The above-mentioned heart rate value of technical solution of the present invention output is instantaneous heart rate value but not average heart rate value, from the secondary R wave detecting, a R ripple often detected and just export in real time current instantaneous heart rate value, the small changes in heart rate of convenient detection, to identify the clinical problems such as heart rate variability; Compared to prior art, obtain the mode of heart rate, as long as detecting R ripple, technical solution of the present invention just can provide instantaneous heart rate accurately, not only from obtaining speed, be improved, can also be convenient to the minor variations of interior heart rate between detection period, improved on the whole the mode that heart rate obtains.
In the possibility of heart rate acquisition methods, also provide a kind of mode of obtaining more exactly the time of R ripple appearance: technical solution of the present invention is used amplitude threshold values to obtain the time that the sampled point on R ripple occurs to obtain R ripple, the time again time corresponding to maximum point of detecting in the sampled point of R ripple being occurred as this R ripple, make heart rate detection result more accurate and be detection in real time; In addition, corresponding each R ripple time, can obtain the temporal instantaneous heart rate of this R ripple, can make rate testing result embody the minor variations of heart rate, have more actual clinical meaning.
In possibility further, the amplitude threshold values using is dynamic, renewable.Described amplitude threshold values can be according to the fluctuation of detected R wave amplitude and real-time update, reduces the false drop rate of R ripple, improves the precision that instantaneous heart rate value detects, and controls the detection error of interval time between the R ripple that this programme obtains in several milliseconds.In addition, the amplitude threshold values of technical solution of the present invention only upgrades according to detected R crest value, is guaranteeing, under the prerequisite of data validity, further to reduce the memory space of system data, thereby is improving data processing rate.
The heart rate acquisition methods providing based on technical solution of the present invention, also provides a kind of processing method of electrocardiosignal.In possibility, a kind of scheme of the abnormal cardiac rate detecting being returned to inspection is also provided, the heart rate abnormal data relating to detecting is examined: technical solution of the present invention considers that the generation of abnormal cardiac rate value may be due to small magnitude R ripple, arrhythmia, the reasons such as amplitude T ripple flase drop cause, thereby when instantaneous heart rate exceeds the first heart rate value of setting and the scope between the second heart rate value, just utilize the Wave data of preserving, based on above-mentioned Wave data, calculating the now wave characteristic of Wave data amasss, the reason of utilizing the long-pending scope of wave characteristic wave characteristic long-pending and R ripple to produce recently distinguishing abnormal cardiac rate value.Utilize above-mentioned scheme of abnormal cardiac rate being returned to inspection, can accurately detect heart rate abnormality (comprising: small magnitude R ripple, rapid heart rate, heart rate are crossed slow or irregular heart rate), prevent the misinterpretation of the abnormal cardiac rate data such as or irregular heart beating abnormal to heart rate.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of a kind of heart rate acquisition methods of prior art;
The schematic flow sheet of a kind of heart rate acquisition methods that Fig. 2 provides for technical solution of the present invention embodiment;
Fig. 3 goes out the method flow schematic diagram of now for a kind of R ripple that technical solution of the present invention embodiment provides;
The method flow schematic diagram of the time that maximum point in a kind of sampled point that obtains described R ripple that Fig. 4 provides for technical solution of the present invention embodiment is corresponding;
The waveform schematic diagram of the electrocardiosignal of a kind of input that Fig. 5 provides for technical solution of the present invention embodiment;
The schematic flow sheet of the processing method of a kind of electrocardiosignal that Fig. 6 provides for technical solution of the present invention embodiment;
The schematic flow sheet of the processing method of a kind of concrete application electrocardiosignal that Fig. 7 provides for technical solution of the present invention embodiment;
The method flow schematic diagram of the abnormal cardiac rate that a kind of screening that Fig. 8 provides for technical solution of the present invention embodiment detects;
The waveform schematic diagram of the electrocardiosignal of a kind of appearance that Fig. 9 provides for technical solution of the present invention embodiment R ripple by a small margin;
A kind of appearance that Figure 10 provides for technical solution of the present invention embodiment is the waveform schematic diagram of the electrocardiosignal of T ripple significantly;
A kind of schematic flow sheet of describing intuitively the processing method of described electrocardiosignal that Figure 11 provides for technical solution of the present invention embodiment.
The specific embodiment
For object of the present invention, feature and effect can be become apparent more, below in conjunction with accompanying drawing, the specific embodiment of the present invention is elaborated.
Set forth in the following description a lot of details so that fully understand the present invention, but the present invention can also adopt other to be different from mode described here, implement, so the present invention has not been subject to the restriction of following public specific embodiment.
The present embodiment provides a kind of heart rate acquisition methods, as shown in Figure 2, comprising:
Step S100, records the moment that on electrocardiosignal waveform, R ripple occurs.
In this step, described electrocardiosignal waveform generally based on be the ecg signal data of instant input, because R ripple is that in whole electrocardiosignal waveform, the most significant ripple and its are periodically comparatively stable compared to other electrocardio ripples, technical solution of the present invention has utilized the waveform characteristic of described R ripple to carry out heart rate detection.
Consider R waveform span in time, its waveform has occupied regular hour section, and the moment that described R ripple occurs can be any time of R waveform in the shared time period.But in this step, also implied the real-time R ripple of electrocardiosignal waveform has been detected, consideration R ripple detects usings its range value as detecting foundation, thereby, can select range value on R ripple to there is compared to other electrocardio ripples the moment that moment that in obvious amplitude, the point of superiority is corresponding occurs as R ripple.
For performing step S100, a kind of R of detection ripple goes out the scheme of now, as shown in Figure 3, comprising:
Step S200, samples to obtain sampled point to described electrocardiosignal waveform.
Step S201, detects R ripple, and range value is greater than to the sampled point of amplitude threshold as the sampled point of described R ripple.
Step S202, the sampled point of searching for described R ripple is to obtain time corresponding to maximum point in the sampled point of described R ripple.
In the process of step S200～S202, described sampled point, comprises the maximum point in sampled point, and corresponding is that described electrocardiosignal wave-shape amplitude is worth time dependent relation, and the information of a sampled point has comprised the range value of electrocardiosignal on special time.Because R ripple has the advantages that waveform is precipitous, amplitude is large, width is narrow, be easy to identification, the present embodiment has been set amplitude threshold and has been identified the sampled point that range value is greater than amplitude threshold, these sampled points that recognize are the sampled point of R ripple, identify thus the R ripple on electrocardiosignal waveform.
Especially, in step S202, maximum point in the sampled point of described R ripple is actual is the peak point of described R ripple, the time occurring for obtaining R crest value point, can pass through the slope detection of R ripple sampled point, and obtain the maximum point of R ripple.R ripple only has a maximum point, according to its amplitude characteristic, when surpassing described amplitude threshold values the slope of R ripple sampled point be first on the occasion of, after R wave crest point, slope becomes negative value, so, as long as detect in real time the positive and negative trip point of slope of input R ripple, just can Real-time Collection to time of R ripple appearance.The method can collect accurately R crest value point and R crest value is put the corresponding moment, compared to obtaining on R ripple, spend moment that the sampled point of amplitude threshold is corresponding as R ripple mode constantly, the R ripple obtaining goes out now can be more accurate and amount of calculation is little, and the error that the R ripple obtaining goes out now can be reduced to several milliseconds by tens milliseconds.
Based on above-mentioned analysis, Fig. 4 has provided the mode of a kind of concrete execution step S202, comprising:
Step S300, calculates the sampled point slope of described R ripple.
Can to each point slope on R ripple, calculate by the slope detection module of R ripple, from R waveform characteristic, the sampled point slope of described R ripple along described R ripple gradually from the occasion of becoming negative value.
Step S301, identification slope is from the sampled point of the R ripple that is negative value on the occasion of saltus step.
Step S302, the time corresponding to sampled point of the R ripple that collection is identified.
What need to further describe is, in above-mentioned discussion, related amplitude threshold dynamically updates in the process of implementation, the time that renewal process can be based on input ecg signal Wave data, upgrade with the fixed cycle, the needs that also can detect according to R ripple aperiodically upgrade.
Wherein, the initial value of described amplitude threshold is determined by detecting the amplitude maximum of the electrocardiosignal sampled point of inputting in a Fixed Time Interval, in this Fixed Time Interval, the amplitude maximum of gained sampled point is multiplied by a threshold coefficient, obtains described amplitude threshold.
Also can said process be described with following formula (1):
C
0=A
max×α ----（1）
Wherein, C
0for the initial value of described amplitude threshold, A
maxfor the amplitude maximum of gained sampled point in this Fixed Time Interval, α is described threshold coefficient.
Under normal conditions, described threshold coefficient can be got the real number between 0～1, and the real number between described 0～1 comprises gets 1 but do not comprise and get 0, and in the present embodiment, described threshold coefficient gets 0.9.
And the mode that dynamically updates described amplitude threshold can be based on step S200～S202 R ripple testing process and the detection of R crest value, the maximum point of the R ripple that can obtain according to step S202 obtains the peak amplitude of R ripple.When continuous acquisition is to setting after the R ripple of number, by calculating the average of the R crest value amplitude of above-mentioned setting number, and the average of above-mentioned peak amplitude is multiplied by described threshold coefficient upgrades amplitude threshold values.
Also can said process be described with following formula (2):
Wherein, C
ifor the amplitude threshold after upgrading, A
i1to A
inbe respectively the 1st to n the R crest value amplitude collecting, n by the setting number of collection R ripple, α is described threshold coefficient.
Suppose to calculate amplitude threshold according to the R ripple number being consecutively detected in 10s, the peak amplitude data that only need to store R ripple in 10s for upgrading described amplitude threshold in the present embodiment, above-mentioned data volume is generally tens data, shared interior poke approximately only accounts for one of percentage, also very good aspect arithmetic speed.
Continuation is with reference to figure 2, and the heart rate acquisition methods that the present embodiment provides, also comprises:
Step S101, the moment T based on the i time R ripple appearance
rimoment T with the i-1 time R ripple appearance
ri-1, obtain T constantly
rion instantaneous heart rate.
Instantaneous heart rate in step S101 is the inverse at each heart time interval, i.e. the inverse of interval between two adjacent R ripples of electrocardiogram.If each heartbeat has small variation in interval, utilize instantaneous heart rate all can detect.Utilize the Real-time Collection of instantaneous heart rate, detect the minor variations of each heart beating instantaneous heart rate for people's health monitoring, exercise intensity judgement,, there is very important effect the aspects such as evaluate application of the early warning of physical fatigue degree and energy expenditure.
Can formula (3) or formula (4) based on following obtain T constantly
rion instantaneous heart rate:
F=F1=60/T----（3）
F=F2=1/T----（4）
Wherein, T is occurred the i time R ripple and the interval between the i-1 time R ripple, and F is T constantly
rion instantaneous heart rate.Work as F=F1, described instantaneous heart rate F is the number of beating of heart in a minute, and the unit of F is " inferior per minute "; Work as F=F2, described instantaneous heart rate F is the number of beating of heart beating in the second, and the unit of F is " inferior per second ".For clinical, obtain in the designed scheme of heart rate, generally select formula (3) output heart rate value.
Technical solution of the present invention provides a kind of heart rate of the adjacent R ripple for detecting to obtain scheme, realizes the real-time detection of heart rate, if having small variation in each heartbeat interval, utilizes this programme all can detect.The extraction of the collection of the present embodiment ecg wave form and R ripple is carried out simultaneously, as long as guarantee to occur R ripple, and time that just can its appearance of immediate record, the time that the R ripple by twice record occurs, can calculate real-time heart rate value, realize real-time the obtaining of heart rate.
In the application examples first of the present embodiment, can detect the electrocardiosignal waveform as shown in Figure 5 of input, the ecg signal data of inputting in Fig. 5 is the waveform of the electrocardiosignal of input in 1 minute.If the heart rate of the prior art that employing background technology is recorded obtains scheme, can only arrive heart rate value is that 76(also has 76 R ripples in 1 minute).But the scheme that adopts the present embodiment can detect the moment T of the 1st R ripple appearance in the electrocardiosignal oscillogram shown in Fig. 5
r1moment T to the 75th R ripple appearance
r75corresponding instantaneous heart rate value (here, the moment that R ripple occurs is also the moment that R wave-wave peak occurs), wherein, during thering are 75 R ripples on this electrocardiosignal waveform.The moment T that continuation occurs with reference to 5, the 1 R ripples of figure
r1moment T to the 75th R ripple appearance
r75corresponding instantaneous heart rate value is respectively: 60, 71, 68, 70, 48, 58, 71, 79, 76, 87, 107, 86, 86, 120, 80, 126, 83, 113, 77, 83, 76, 82, 63, 76, 88, 122, 122, 53, 62, 60, 62, 77, 69, 83, 70, 122, 116, 61, 39, 60, 62, 70, 74, 83, 72, 89, 77, 69, 71, 70, 75, 122, 135, 137, 81, 81, 104, 57, 78, 67, 82, 75, 76, 63, 57, 73, 60, 75, 83, 120, 70, 61, 81, 78, 73, such as, the moment T of the 1st R ripple appearance
r1corresponding instantaneous heart rate value is the moment T of 60, the 75 R ripple appearance
r75corresponding instantaneous heart rate value is 73.
The result of obtaining from above-mentioned heart rate, the scheme that prior art heart rate obtains can be lost the abnormal information of a large amount of hearts rate, can not reflect real changes in heart rate.The heart rate of the present embodiment obtains scheme can just export current heart rate value when a R ripple often being detected, can detect small changes in heart rate, is beneficial to and identifies exactly abnormal heart rate information.
Heart rate acquisition methods based on the present embodiment, the present embodiment also provides a kind of processing method of electrocardiosignal, as shown in Figure 6, comprising:
Step S400, the electrocardiosignal based on input is obtained T constantly
ritime heart rate.
Can utilize the heart rate acquisition methods of the present embodiment to obtain T constantly according to the electrocardiosignal of described input
riinstant heart rate.
Step S401, processes described moment T
ritime heart rate.
In step S401, the described moment T of described processing
ritime heart rate be intended to the heart rate value that gets of application, described application can comprise the abnormal cardiac rate that detects in heart rate, abnormal cardiac rate that screening detects, heart rate that transmission detects and the abnormal cardiac rate the detecting step of aspect such as report to the police is carried out.
Described step S401 can comprise the step of carrying out above-mentioned application function, also can comprise the part steps of carrying out in above-mentioned application function step, and other steps in above-mentioned application function step are performed between step S400 and step S401.In the present embodiment, the described moment T of described processing
ritime heart rate comprise the abnormal cardiac rate of carrying out screening and detecting, the heart rate that transmission detects and the step that the abnormal cardiac rate detecting is reported to the police, abnormal cardiac rate in described detection heart rate is executed between step S400 and step S401, can be with reference to figure 7, the application examples of the processing method of the electrocardiosignal first of the present embodiment can comprise the steps:
Step S500, can be with step S400;
Step S501, detects the abnormal cardiac rate in heart rate.
Step S502, the abnormal cardiac rate that screening detects.
Step S503, based on the selection result, the heart rate that transmission detects is also reported to the police to abnormal cardiac rate.
Wherein, step S502～step S503 is the specific embodiment of step S400.
Be different from prior art, step S501～step S503 has realized the inspection of returning of abnormal cardiac rate value has been processed.Concrete:
In step S501, abnormal cardiac rate in described detection heart rate can be directly by the real-time heart rate of heart rate value range detection of setting: the lower limit heart rate value of setting in heart rate value scope is the first heart rate value, upper limit heart rate value is the second heart rate value, such as, normal human's heart rate can be with situation fluctuations such as motions, when motion, the maximum heart rate of human body can reach left and right per minute 150 times, and therefore, it is 150 that the second heart rate value can be set, the first heart rate value is 40, and described heart rate value scope is 40～150.Certainly, also can the first heart rate value and second heart rate value of individuation be set separately according to individual physique, to obtain normal human's heart rate value scope.Usually, the natural number between described the first heart rate value desirable 40～50, the natural number between described the second heart rate value desirable 150～160.As described moment T
ritime heart rate when being less than the first heart rate value or being greater than the second heart rate value, can think described moment T
ritime heart rate be abnormal cardiac rate; As described moment T
ritime heart rate between the first heart rate value and the second heart rate value, can think described moment T
ritime heart rate be normal cardiac rate.
Step S502 is mainly used in confirming the generation of abnormal heart rate value, search the reason that abnormal cardiac rate produces: consider the generation of abnormal cardiac rate may be owing to detecting the reason such as small magnitude R ripple, R ripple are undetected, arrhythmia, amplitude T ripple flase drop cause, wherein, undetected and the amplitude T ripple flase drop of R ripple is not human heart rate's abnormal attribute, can be rejected the abnormal cardiac rate therefore producing; So, when the heart rate value detecting is not between the first heart rate value and the second heart rate value, just utilize the electrocardiosignal Wave data of preserving to confirm the reason that abnormal cardiac rate produces.Step S502 in the present embodiment has utilized the long-pending contrast of heart rate volatility characteristic wave characteristic long-pending and R ripple, distinguishes the reason that abnormal cardiac rate value produces, thereby realizes the screening process of abnormal cardiac rate.The process that realizes the abnormal cardiac rate that described screening detects can as shown in Figure 8, comprise the steps:
Step S600, as described moment T
ritime heart rate be abnormal cardiac rate, obtain constantly T
ritime electrocardiosignal waveform wave characteristic long-pending.
The long-pending absolute value that refers to the product of the n section slope of curve that front and back are adjacent of described wave characteristic, wherein, n is more than or equal to 2 natural number.When the heart rate value detecting is abnormal cardiac rate, the time being occurred by abnormal cardiac rate, can in the array of stored waveform, navigate to the time point that abnormal cardiac rate occurs, the Wave data of the R ripple of extraction time point appearance, obtains T constantly
ritime electrocardiosignal waveform wave characteristic long-pending.The described T constantly that obtains
ritime electrocardiosignal waveform long-pending comprising of wave characteristic:
Extract T constantly
ri-1to moment T
ribetween ecg wave form, described moment T
ri-1to moment T
ribetween ecg wave form be divided into N section wavy curve, wherein N is greater than 2 natural number;
Calculate respectively the slope of described N section wavy curve, described moment T
ritime electrocardiosignal waveform wave characteristic amass the absolute value into the slope product of described N section wavy curve.
Can the long-pending acquisition process of above-mentioned wave characteristic be described based on following formula (5):
V=k
1×k
2×...×k
N----（5）
Wherein, V is moment T
ritime electrocardiosignal waveform wave characteristic long-pending, k
1to k
nbe respectively T constantly
ri-1to moment T
ribetween the slope of the 1st section of wavy curve to N section.
Continuation is with reference to figure 8, and the process of the abnormal cardiac rate that described screening detects also comprises:
Step S601, more described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform is long-pending, to obtain the selection result.
The present embodiment will be by the described moment T of judgement
riwhether the wave characteristic of electrocardiosignal waveform is long-pending to be that the wave characteristic of R ripple is long-pending produce constantly (i.e. moment T as abnormal cardiac rate
ri) surge detection benchmark, for abnormal cardiac rate value, return inspection.
From the wave character of electrocardiosignal, analyze, its waveform after filtering after, only have the undulatory property of QRS wave group maximum; The wave characteristic of QRS wave group is long-pending is generally 70～120, and the wave characteristic of other waveforms is long-pending generally in 10, and the wave characteristic of QRS wave group is long-pending several, and to decuple the wave characteristics of other wave band long-pending.But there is small magnitude R ripple (the crest value of this R ripple is less) in R ripple, its wave characteristic amasss more than 15, with other waveforms and possess that the wave characteristic of general QRS wave group is long-pending produces obvious gap.By above-mentioned analysis, can set the long-pending scope of living in of wave characteristic of described R ripple.
In the present embodiment, the long-pending scope of living in of the wave characteristic of described R ripple comprises the first fluctuation range and the second fluctuation range that is greater than described the second fluctuation threshold value between the first fluctuation threshold value and the second fluctuation threshold value, and described the first fluctuation threshold value is less than the second fluctuation threshold value.Wherein, described the first fluctuation range is R ripple scope of living in by a small margin, and described the second fluctuation range is general QRS wave group scope of living in.Can utilize described the first fluctuation range detection small magnitude R ripple, utilize the general QRS wave group of described the second fluctuation range detection.Can set the span of above-mentioned fluctuation threshold value.The span of described the first fluctuation threshold value is the real number between 10～20, and the span of described the second fluctuation threshold value is the real number between 60～70.It is that 15, the second fluctuation threshold values are 60 that the present embodiment is got the first fluctuation threshold value, and described the second fluctuation range is 70～120.
By step S601, as described moment T
ritime heart rate be abnormal cardiac rate, described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform amasss and comprises following comparison procedure:
At described moment T
riheart rate while being less than the first heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic long-pending in described the first fluctuation range, show T constantly
ritime electrocardiosignal R ripple undulatory property less, amplitude is also less, constantly T
ritime electrocardiosignal there is R ripple by a small margin.Its undulatory property feature can be with reference to figure 9, in Fig. 9, and T constantly
ritime electrocardiosignal there is R ripple u1 by a small margin.It should be noted that, it is significant to clinical practice that R ripple is by a small margin carried out detecting separately.
At described moment T
riheart rate while being less than the first heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly greater than described the second fluctuation threshold value, show T constantly
ritime there is the undetected phenomenon of R ripple.R ripple is undetected refer to original for normal R ripple (being that its wave characteristic is long-pending in the R of the second fluctuation range ripple) adopt the present embodiment amplitude threshold carry out R ripple while detecting because of the too high and undetected situation of amplitude threshold, undetected R waveguide causes the abnormal cardiac rate that now obtained and has error.The screening process of the present embodiment can detect above-mentioned undetected phenomenon.
At described moment T
riheart rate while being greater than the second heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly less than described the first fluctuation threshold value, show flase drop T ripple.T ripple flase drop refers to T constantly
ritime electrocardiosignal waveform because it has the significantly T ripple of value, make in heart rate acquisition process, significantly the T ripple of value (mainly referring to that its peak value is larger) has surpassed described amplitude threshold and has been used as R ripple flase drop; T ripple flase drop also causes the abnormal cardiac rate that now obtained to have error.During flase drop T ripple, the undulatory property feature of electrocardiosignal can be with reference to Figure 10, in Figure 10, and T constantly
ritime electrocardiosignal there is significantly T ripple u2, because the peak amplitude of T ripple u2 is larger, when described amplitude threshold is less than the peak value of T ripple u2, T ripple u2 is used as to R ripple and detects, thereby affect the accuracy of heart rate detection, also caused the abnormal cardiac rate of non-physiologic (non-pathologic).
Certainly, for other abnormal cardiac rates, screening ascribes the other reasonses such as arrhythmia to.
Based on above-mentioned the selection result, the transferring content of step S503 is relevant with described the selection result:
For non-abnormal cardiac rate, be also normal cardiac rate, can directly to described heart rate value, transmit;
For abnormal cardiac rate, its transmission means is as follows:
The less abnormal cardiac rate of heart rate value that cause R ripple by a small margin causes, transmits corresponding heart rate value, and can attach the warning of transmission R ripple by a small margin.
Because of the less abnormal cardiac rate of the undetected heart rate value causing of R ripple, corresponding heart rate value is abandoned to processing, get rid of to the transmission of corresponding heart rate value.
Because of the excessive abnormal cardiac rate of heart rate value that significantly the T ripple of value causes, corresponding heart rate value is abandoned to processing, get rid of the transmission to the corresponding heart rate value of institute, but can attach the warning of transmission T ripple flase drop.
The excessive abnormal cardiac rate of heart rate value because the other reasonses such as arrhythmia cause, transmits corresponding heart rate value, can attach and send rapid heart rate warning.
The little heart rate value of transmission abnormality, sends the alarm of little R ripple
Can also, with flow chart as shown in figure 11, above-mentioned processing procedure be described intuitively.
With reference to Figure 11, based on described input ecg signal, obtain heart rate value process can refer step S100～S101, described input ecg signal is held to form Wave data, can, according to the time of input ecg signal, sequentially form the waveform array that stores described electrocardiosignal; The abnormal process of described preliminary judgement heart rate can refer step S501; For normal cardiac rate value, can directly transmit described heart rate value; For abnormal data, Wave data that can be based on maintained electrocardiosignal, the abnormal reason of judgement heart rate, to return inspection and to screen to described abnormal cardiac rate, selectively heart rate value transmitted or to abandon, and selectively report to the police based on described abnormal cardiac rate; Said process can be carried out according to step S502～step S503.
Although the present invention with preferred embodiment openly as above; but it is not for limiting the present invention; any those skilled in the art without departing from the spirit and scope of the present invention; can utilize method and the technology contents of above-mentioned announcement to make possible change and modification to technical solution of the present invention; therefore; every content that does not depart from technical solution of the present invention; any simple modification, equivalent variations and the modification above embodiment done according to technical spirit of the present invention, all belong to the protection domain of technical solution of the present invention.
Claims (14)
1. a heart rate acquisition methods, is characterized in that, comprising:
Record the moment that on electrocardiosignal waveform, R ripple occurs;
Moment T based on the i time R ripple appearance
rimoment T with the i-1 time R ripple appearance
ri-1, obtain T constantly
rion instantaneous heart rate.
2. heart rate acquisition methods as claimed in claim 1, is characterized in that, the described moment of recording R ripple appearance on electrocardiosignal waveform comprises:
Described electrocardiosignal waveform is sampled to obtain sampled point, and the corresponding described electrocardiosignal wave-shape amplitude of described sampled point is worth time dependent relation;
Detect R ripple, range value is greater than to the sampled point of amplitude threshold as the sampled point of described R ripple;
The sampled point of searching for described R ripple is to obtain time corresponding to maximum point in the sampled point of described R ripple.
3. heart rate acquisition methods as claimed in claim 2, is characterized in that, the sampled point of the described R ripple of described search comprised with time corresponding to maximum point obtaining in the sampled point of described R ripple:
Calculate the sampled point slope of described R ripple, the sampled point slope of described R ripple along described R ripple gradually from the occasion of becoming negative value;
Identification slope is from the sampled point of the R ripple that is negative value on the occasion of saltus step;
Gather the time corresponding to sampled point of the R ripple of identifying.
4. heart rate acquisition methods as claimed in claim 2, is characterized in that, also comprises: upgrade described amplitude threshold, the described amplitude threshold of described renewal comprises:
Obtain the amplitude equalizing value of the maximum point of the several R ripples that are consecutively detected;
Obtained amplitude equalizing value is multiplied by threshold value coefficient, to obtain new amplitude threshold.
5. heart rate acquisition methods as claimed in claim 4, is characterized in that, the amplitude threshold upgrading comprises the amplitude threshold that upgraded and initial amplitude threshold, and described heart rate acquisition methods also comprises:
Obtain the amplitude peak value in the some sampled points that are consecutively detected;
Obtained amplitude peak value is multiplied by described threshold value coefficient, to obtain described initial amplitude threshold.
6. heart rate acquisition methods as claimed in claim 5, is characterized in that, the span of described threshold value coefficient is between 0～1, to be not equal to 0 real number.
7. heart rate acquisition methods as claimed in claim 1, is characterized in that, based on following formula, obtains T constantly
rion instantaneous heart rate:
F=60/T; Or,
F=1/T；
Wherein, T is occurred the i time R ripple and the interval between the i-1 time R ripple, and F is T constantly
rion instantaneous heart rate.
8. a processing method for electrocardiosignal, is characterized in that, comprising:
Electrocardiosignal based on input, utilizes the method as described in claim 1～7 any one to obtain T constantly
ritime heart rate;
Process described moment T
ritime heart rate.
9. the processing method of electrocardiosignal as claimed in claim 8, is characterized in that, also comprises: as described moment T
ritime heart rate when being less than the first heart rate value or being greater than the second heart rate value, obtain T constantly
ritime electrocardiosignal waveform wave characteristic long-pending;
The described moment T of described processing
riheart rate comprise: more described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform is long-pending, to transmit described moment T
riheart rate.
10. the processing method of electrocardiosignal as claimed in claim 9, is characterized in that, the described moment T of described processing
ritime heart rate also comprise: as described moment T
ritime heart rate between the first heart rate value and the second heart rate value, transmit described moment T
riheart rate.
The processing method of 11. electrocardiosignaies as claimed in claim 9, is characterized in that, described in obtain constantly T
ritime electrocardiosignal waveform long-pending comprising of wave characteristic:
Extract T constantly
ri-1to moment T
ribetween ecg wave form, described moment T
ri-1to moment T
ribetween ecg wave form be divided into N section wavy curve, wherein N is greater than 2 natural number;
Calculate respectively the slope of described N section wavy curve, described moment T
ritime electrocardiosignal waveform wave characteristic amass the absolute value into the slope product of described N section wavy curve.
The processing method of 12. electrocardiosignaies as claimed in claim 9, it is characterized in that, the long-pending scope of living in of wave characteristic of described R ripple comprises the first fluctuation range and the second fluctuation range that is greater than described the second fluctuation threshold value between the first fluctuation threshold value and the second fluctuation threshold value, and described the first fluctuation threshold value is less than the second fluctuation threshold value.
The processing method of 13. electrocardiosignaies as claimed in claim 12, is characterized in that, the span of described the first fluctuation threshold value is the real number between 10～20, and the span of described the second fluctuation threshold value is the real number between 60～70.
The processing method of 14. electrocardiosignaies as claimed in claim 12, is characterized in that, described moment T
rithe wave characteristic wave characteristic long-pending and R ripple of electrocardiosignal waveform is long-pending to process described moment T
riheart rate comprise:
At described moment T
riheart rate while being less than the first heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic long-pending in described the first fluctuation range, transmit described moment T
ritime heart rate, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly greater than described the second fluctuation threshold value, abandon described moment T
ritime heart rate;
At described moment T
riheart rate while being greater than the second heart rate value, if described moment T
ritime electrocardiosignal waveform wave characteristic long-pending in described the first fluctuation range or the second fluctuation range, transmit described moment T
ritime heart rate, if described moment T
ritime electrocardiosignal waveform wave characteristic be long-pendingly less than described the first fluctuation threshold value, abandon described moment T
ritime heart rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310598475.0A CN103690156B (en) | 2013-11-22 | 2013-11-22 | The processing method of a kind of heart rate acquisition methods and electrocardiosignal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310598475.0A CN103690156B (en) | 2013-11-22 | 2013-11-22 | The processing method of a kind of heart rate acquisition methods and electrocardiosignal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103690156A true CN103690156A (en) | 2014-04-02 |
| CN103690156B CN103690156B (en) | 2016-01-27 |
Family
ID=50351942
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310598475.0A Active CN103690156B (en) | 2013-11-22 | 2013-11-22 | The processing method of a kind of heart rate acquisition methods and electrocardiosignal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103690156B (en) |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104367344A (en) * | 2014-10-10 | 2015-02-25 | 广东工业大学 | Method and device for detecting instant heart rate of fetus on basis of Shanon envelope |
| CN105125199A (en) * | 2015-07-29 | 2015-12-09 | 北京医康世纪科技有限公司 | Heart rate detection method and heart rate detection device |
| CN105212919A (en) * | 2014-06-11 | 2016-01-06 | 成功大学 | Physiological signal analysis method and system |
| CN105825071A (en) * | 2016-05-13 | 2016-08-03 | 北京精密机电控制设备研究所 | Sensorless heart rate estimation device and estimation method thereof |
| CN105997042A (en) * | 2016-04-28 | 2016-10-12 | 江苏盖睿健康科技有限公司 | Heart rate acquisition method capable of preventing false identification of tall peaked T waves |
| CN106308791A (en) * | 2016-08-29 | 2017-01-11 | 郑州大学第附属医院 | Auxiliary data long distance transmission method and device for electrocardiosignal recording system |
| CN106383808A (en) * | 2016-09-18 | 2017-02-08 | 时瑞科技（深圳）有限公司 | Universal heart rate and electrocardiogram quick calculation system and method |
| CN106529425A (en) * | 2016-10-20 | 2017-03-22 | 包磊 | R-wave extraction method and system for ECG signals |
| CN106604679A (en) * | 2014-09-09 | 2017-04-26 | 日本电信电话株式会社 | Heartbeat detecting method and heartbeat detecting device |
| CN107041740A (en) * | 2016-02-05 | 2017-08-15 | 南京国雅信息科技有限公司 | Animal heart rate monitoring system and the Heart Rate States recognition methods based on neutral net |
| CN107233093A (en) * | 2017-05-23 | 2017-10-10 | 安徽华米信息科技有限公司 | Detection method, device and the electronic equipment of R ripples |
| CN107847146A (en) * | 2015-07-21 | 2018-03-27 | 皇家飞利浦有限公司 | The method and system of automatic mark activity in ECG data |
| CN108154085A (en) * | 2017-12-06 | 2018-06-12 | 北京顺源开华科技有限公司 | The method, apparatus and electronic equipment of identification are carried out based on electrocardiogram (ECG) data |
| TWI635852B (en) * | 2017-06-09 | 2018-09-21 | 國立臺灣科技大學 | Signal processing method for plural sensors detecting pulse in a palm |
| CN109497994A (en) * | 2019-01-16 | 2019-03-22 | 上海掌门科技有限公司 | A kind of detection method of R wave of electrocardiosignal, equipment and computer-readable medium |
| CN109497995A (en) * | 2019-01-16 | 2019-03-22 | 上海掌门科技有限公司 | Detection method, equipment and the computer-readable medium of R wave of electrocardiosignal |
| CN110232132A (en) * | 2019-06-18 | 2019-09-13 | 北京天泽智云科技有限公司 | Time series data processing method and processing device |
| CN110857975A (en) * | 2018-08-22 | 2020-03-03 | 英飞凌科技股份有限公司 | Radar range accuracy improvement method |
| WO2020133536A1 (en) * | 2018-12-29 | 2020-07-02 | 深圳迈瑞生物医疗电子股份有限公司 | Sleep state determining method and apparatus |
| CN111714110A (en) * | 2020-05-19 | 2020-09-29 | 成都云卫康医疗科技有限公司 | Real-time heart rate calculation method based on PPG waveform |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3599033A (en) * | 1969-09-22 | 1971-08-10 | Hewlett Packard Co | System for simultaneous display of ecg and heart rate |
| JP2002000575A (en) * | 2000-06-23 | 2002-01-08 | Matsushita Electric Ind Co Ltd | Heart rate sensor and method of calculating heart rate |
| CN1552283A (en) * | 2003-05-27 | 2004-12-08 | 郭博昭 | Heart rate variance analyzing method and device |
| CN100372499C (en) * | 2004-02-26 | 2008-03-05 | 阮炯 | Method and apparatus for detecting and analysing heart rate variation Co mode degree index |
| CN201248704Y (en) * | 2008-06-24 | 2009-06-03 | 王培勇 | System for real-time detection of exercise heart rate and video data superposition |
| CN103222864B (en) * | 2013-04-07 | 2015-03-25 | 广东工业大学 | Self-adaption electrocardiograph (ECG) detection method and monitoring system thereof |
-
2013
- 2013-11-22 CN CN201310598475.0A patent/CN103690156B/en active Active
Cited By (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105212919A (en) * | 2014-06-11 | 2016-01-06 | 成功大学 | Physiological signal analysis method and system |
| CN105212919B (en) * | 2014-06-11 | 2018-10-19 | 成功大学 | physiological signal analysis system |
| US10750969B2 (en) | 2014-09-09 | 2020-08-25 | Nippon Telegraph And Telephone Corporation | Heartbeat detection method and heartbeat detection device |
| CN106604679A (en) * | 2014-09-09 | 2017-04-26 | 日本电信电话株式会社 | Heartbeat detecting method and heartbeat detecting device |
| CN104367344A (en) * | 2014-10-10 | 2015-02-25 | 广东工业大学 | Method and device for detecting instant heart rate of fetus on basis of Shanon envelope |
| CN107847146B (en) * | 2015-07-21 | 2021-04-27 | 皇家飞利浦有限公司 | Method and system for automatically tagging activity on ECG data |
| CN107847146A (en) * | 2015-07-21 | 2018-03-27 | 皇家飞利浦有限公司 | The method and system of automatic mark activity in ECG data |
| CN105125199B (en) * | 2015-07-29 | 2017-12-26 | 北京医康世纪科技有限公司 | A kind of heart rate detection method and device |
| CN105125199A (en) * | 2015-07-29 | 2015-12-09 | 北京医康世纪科技有限公司 | Heart rate detection method and heart rate detection device |
| CN107041740A (en) * | 2016-02-05 | 2017-08-15 | 南京国雅信息科技有限公司 | Animal heart rate monitoring system and the Heart Rate States recognition methods based on neutral net |
| CN105997042A (en) * | 2016-04-28 | 2016-10-12 | 江苏盖睿健康科技有限公司 | Heart rate acquisition method capable of preventing false identification of tall peaked T waves |
| CN105825071B (en) * | 2016-05-13 | 2019-06-25 | 北京精密机电控制设备研究所 | A kind of no sensor heart rate estimation device and its evaluation method |
| CN105825071A (en) * | 2016-05-13 | 2016-08-03 | 北京精密机电控制设备研究所 | Sensorless heart rate estimation device and estimation method thereof |
| CN106308791A (en) * | 2016-08-29 | 2017-01-11 | 郑州大学第附属医院 | Auxiliary data long distance transmission method and device for electrocardiosignal recording system |
| CN106383808A (en) * | 2016-09-18 | 2017-02-08 | 时瑞科技（深圳）有限公司 | Universal heart rate and electrocardiogram quick calculation system and method |
| CN106383808B (en) * | 2016-09-18 | 2019-08-02 | 时瑞科技（深圳）有限公司 | The processing system and method for heart rate electrocardiosignal |
| CN106529425A (en) * | 2016-10-20 | 2017-03-22 | 包磊 | R-wave extraction method and system for ECG signals |
| CN106529425B (en) * | 2016-10-20 | 2019-11-29 | 深圳市善行医疗科技有限公司 | A kind of the R wave extracting method and system of electrocardiosignal |
| CN107233093A (en) * | 2017-05-23 | 2017-10-10 | 安徽华米信息科技有限公司 | Detection method, device and the electronic equipment of R ripples |
| CN107233093B (en) * | 2017-05-23 | 2020-05-12 | 安徽华米智能科技有限公司 | R wave detection method and device and electronic equipment |
| TWI635852B (en) * | 2017-06-09 | 2018-09-21 | 國立臺灣科技大學 | Signal processing method for plural sensors detecting pulse in a palm |
| CN108154085A (en) * | 2017-12-06 | 2018-06-12 | 北京顺源开华科技有限公司 | The method, apparatus and electronic equipment of identification are carried out based on electrocardiogram (ECG) data |
| CN108154085B (en) * | 2017-12-06 | 2022-02-18 | 北京顺源开华科技有限公司 | Method and device for identity recognition based on electrocardiogram data and electronic equipment |
| CN110857975A (en) * | 2018-08-22 | 2020-03-03 | 英飞凌科技股份有限公司 | Radar range accuracy improvement method |
| WO2020133536A1 (en) * | 2018-12-29 | 2020-07-02 | 深圳迈瑞生物医疗电子股份有限公司 | Sleep state determining method and apparatus |
| CN109497995A (en) * | 2019-01-16 | 2019-03-22 | 上海掌门科技有限公司 | Detection method, equipment and the computer-readable medium of R wave of electrocardiosignal |
| CN109497995B (en) * | 2019-01-16 | 2021-08-24 | 上海掌门科技有限公司 | Electrocardiosignal R wave detection method, electrocardiosignal R wave detection equipment and computer readable medium |
| CN109497994B (en) * | 2019-01-16 | 2021-08-24 | 上海掌门科技有限公司 | Electrocardiosignal R wave detection method, equipment and computer readable medium |
| CN109497994A (en) * | 2019-01-16 | 2019-03-22 | 上海掌门科技有限公司 | A kind of detection method of R wave of electrocardiosignal, equipment and computer-readable medium |
| CN110232132A (en) * | 2019-06-18 | 2019-09-13 | 北京天泽智云科技有限公司 | Time series data processing method and processing device |
| CN110232132B (en) * | 2019-06-18 | 2020-11-06 | 北京天泽智云科技有限公司 | Time series data processing method and device |
| CN111714110A (en) * | 2020-05-19 | 2020-09-29 | 成都云卫康医疗科技有限公司 | Real-time heart rate calculation method based on PPG waveform |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103690156B (en) | 2016-01-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103690156B (en) | The processing method of a kind of heart rate acquisition methods and electrocardiosignal | |
| US11147500B2 (en) | Electrocardiogram processing system for delineation and classification | |
| CN101449973B (en) | Judgment index generation method and device for cardiac interference signal identification | |
| US8543194B2 (en) | System and method of detecting abnormal movement of a physical object | |
| CN101856225B (en) | Method for detecting R wave crest of electrocardiosignal | |
| CN107837082A (en) | Electrocardiogram automatic analysis method and device based on artificial intelligence self study | |
| CN101828918B (en) | Electrocardiosignal R peak detection method based on waveform characteristic matching | |
| US20150080752A1 (en) | Device for automatic mapping of complex fractionated atrial electrogram | |
| CN109411041A (en) | Ecg information processing method and electro cardio signal workstation system | |
| CN101065058A (en) | Monitoring physiological activity using partial state space reconstruction | |
| CN105796096A (en) | Heart rate variability analysis method, heart rate variability analysis system and terminal | |
| CN104586381A (en) | Electrocardiograph monitoring system based on Internet of Things | |
| CN102579039B (en) | Method for detecting TWA (T wave alternans) in electrocardiogram | |
| CN109411042A (en) | Ecg information processing method and electro cardio signal workstation | |
| CN102379692A (en) | System and method for remote electrocardiographic monitoring for area-end networks | |
| Ng et al. | Iterative method to detect atrial activations and measure cycle length from electrograms during atrial fibrillation | |
| CN103970975A (en) | Electrocardio data processing method and electrocardio data processing system | |
| CN109288515B (en) | Periodicity monitoring method and device based on premature beat signal in wearable electrocardiosignal | |
| CN110226919A (en) | Electrocardiosignal type detection method, apparatus, computer equipment and storage medium | |
| CN105411579B (en) | A kind of electrocardiogram R wave detection method and device | |
| Gupta et al. | A derivative-based approach for QT-segment feature extraction in digitized ECG record | |
| CN109431487B (en) | Method for identifying typical data of ventricular heart beat in electrocardiogram dynamic real-time analysis data | |
| US11672464B2 (en) | Electrocardiogram processing system for delineation and classification | |
| CN106974644A (en) | A kind of electrocardiogram R wave detection method and system | |
| CN105212921B (en) | A kind of lie detecting method estimated based on electrocardio |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: Room 2070, building 2, Lane 1800, Xinyang Road, Lingang New Area, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai, 200120 Patentee after: NEUSOFT XIKANG HEALTHCARE TECHNOLOGY Co.,Ltd. Address before: Room 12b01, 03, 05, 07, 09, No. 9, Beisihuan West Road, Haidian District, Beijing Patentee before: NEUSOFT XIKANG HEALTHCARE TECHNOLOGY Co.,Ltd. |