CN109497984A - A kind of processing method and relevant device of heart rate signal - Google Patents
A kind of processing method and relevant device of heart rate signal Download PDFInfo
- Publication number
- CN109497984A CN109497984A CN201811612483.5A CN201811612483A CN109497984A CN 109497984 A CN109497984 A CN 109497984A CN 201811612483 A CN201811612483 A CN 201811612483A CN 109497984 A CN109497984 A CN 109497984A
- Authority
- CN
- China
- Prior art keywords
- signal
- ppg
- jump
- heart rate
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 52
- 230000035772 mutation Effects 0.000 claims abstract description 29
- 238000012216 screening Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000001914 filtration Methods 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 32
- 230000001133 acceleration Effects 0.000 claims description 22
- 230000003313 weakening effect Effects 0.000 claims description 12
- 230000003044 adaptive effect Effects 0.000 claims description 11
- 238000009434 installation Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 26
- 230000033001 locomotion Effects 0.000 description 16
- 210000004369 blood Anatomy 0.000 description 8
- 239000008280 blood Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 210000001519 tissue Anatomy 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 238000013186 photoplethysmography Methods 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 230000003447 ipsilateral effect Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The embodiment of the present application discloses a kind of processing method of heart rate signal, for realizing the optimization to heart rate value detection process.The embodiment of the present application method includes: to obtain PPG signal；Screening identification is carried out to the PPG signal, determines the jump signal in the PPG signal；Mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal；Heart rate value is calculated according to the echo signal.Wherein, mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal and then heart rate value is calculated according to the echo signal, to reduce influence of the jump signal to heart rate value calculating process in PPG signal, the optimization to heart rate value detection process is realized.
Description
Technical field
This application involves intelligent wearable device field more particularly to the processing methods and relevant device of a kind of heart rate signal.
Background technique
Exercise is promotive of health, strong, and it is desirable to pass through the physical fitness of movement lifts oneself.The form of training,
Sporter can select the suitable movement of oneself according to the physical function state of oneself.Movement can monitor, and pass through prison
Heart rate measuring can guarantee that sporter will not cause movement effects unknown because of exercise intensity is too low to control the rhythm of motion exercise
It is aobvious, it is also possible to prevent that sporter's exercise intensity is excessive to lead to injury gained in sports.
In the prior art, the common heart rate detection mode of intelligent wearable device is detected in the tissue by photoelectric method
Volumetric blood variation obtains pulse signal.Intelligent wearable device control illumination is mapped to tested tissue, decays by tissue resorption
After effect, from light source is ipsilateral or heteropleural is received by photelectric receiver, some structural constituents such as muscle, bone etc. in tested tissue
The effect of light is considered as constant, and volumetric blood can change with the contraction of heart with diastole in arteries,
Cause blood to change the uptake of light, cause output intensity that cyclically-varying is presented, the variation of final light intensity is converted to
Electric signal is photoplethysmography (photoplethysmography, PPG) signal.To the PPG signal processing fortune of acquisition
It calculates, the value of heart rate and blood oxygen saturation can be acquired.
However, using photoelectric sphyg wave carry out heart rate detection when, since signal is weaker, be highly prone to it is various bounce or not
The interference of stabilization signal leads to the unstability of overall signal since movement generates the influence of jump, often normal heart rate letter
It number is completely covered, to affect the normal calculating of heart rate.
Summary of the invention
The embodiment of the present application provides a kind of processing method of heart rate signal, for realizing to the excellent of heart rate value detection process
Change.
The embodiment of the present application first aspect provides a kind of processing method of heart rate signal, comprising:
Obtain PPG signal；
Screening identification is carried out to the PPG signal, determines the jump signal in the PPG signal；
Mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal；
Heart rate value is calculated according to the echo signal.
Optionally, the jump signal packet that screening identification is carried out to the PPG signal, determines in the PPG signal
It includes:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described
The presupposition multiple of second mean value.
Optionally, it carries out mutation to the jump signal in the PPG signal to handle, obtaining echo signal includes:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than described
Jump signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
Optionally, the acquisition PPG signal includes:
Obtain PPG original signal；
The PPG original signal is normalized, second signal is obtained；
The second signal is filtered, the PPG signal is obtained.
Optionally, the second signal is filtered, obtaining the PPG signal includes:
High-pass filtering and/or low-pass filtering are carried out to the second signal, obtain the PPG signal；
And/or
Acceleration signal is obtained, there are incidence relations with the PPG original signal for the acceleration signal；
Adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
The embodiment of the present application second aspect provides a kind of optimization system of heart rate signal detection, comprising:
Acquiring unit, for obtaining PPG signal；
Determination unit determines the jump signal in the PPG signal for carrying out screening identification to the PPG signal；
Processing unit is handled for carrying out mutation to the jump signal in the PPG signal, obtains target letter
Number；
Computing unit, for calculating heart rate value according to the echo signal.
Optionally, the determination unit is specifically used for:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described
The presupposition multiple of second mean value.
Optionally, the processing unit is specifically used for:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than described
Jump signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
Optionally, the acquiring unit is specifically used for:
Obtain PPG original signal；
The PPG original signal is normalized, second signal is obtained；
The second signal is filtered, the PPG signal is obtained.
Optionally, the acquiring unit is specifically used for:
High-pass filtering and/or low-pass filtering are carried out to the second signal, obtain the PPG signal；
And/or
Acceleration signal is obtained, there are incidence relations with the PPG original signal for the acceleration signal；
Adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
The embodiment of the present application third aspect provides a kind of computer installation, comprising:
Processor, memory, input-output equipment and bus；
The processor, memory, input-output equipment are connected with the bus respectively；
The processor is for executing such as the described in any item methods of preceding method embodiment.
The embodiment of the present application fourth aspect provides a kind of computer readable storage medium, is stored thereon with computer journey
Sequence, it is characterised in that: the side as described in any one of preceding method embodiment is realized when the computer program is executed by processor
The step of method.
As can be seen from the above technical solutions, the embodiment of the present application, which has the advantage that, obtains PPG signal；To the PPG
Signal carries out screening identification, determines the jump signal in the PPG signal；To the jump signal in the PPG signal
It carries out mutation to handle, obtains echo signal；Heart rate value is calculated according to the echo signal.Wherein, in the PPG signal
Mutation is carried out to the jump signal to handle, and obtains echo signal and then heart rate value is calculated according to the echo signal,
To reduce influence of the jump signal to heart rate value calculating process in PPG signal, the optimization to heart rate value detection process is realized.
Detailed description of the invention
Fig. 1 is an a kind of schematic diagram of the processing method embodiment of heart rate signal in the embodiment of the present application；
Fig. 2 is a kind of another schematic diagram of the processing method embodiment of heart rate signal in the embodiment of the present application；
Fig. 3 is a kind of another schematic diagram of the processing method embodiment of heart rate signal in the embodiment of the present application；
Fig. 4 is a kind of another schematic diagram of the processing method embodiment of heart rate signal in the embodiment of the present application；
Fig. 5 is a kind of another schematic diagram of the processing method embodiment of heart rate signal in the embodiment of the present application；
Fig. 6 is an a kind of schematic diagram of the processing system embodiment of heart rate signal in the embodiment of the present application；
Fig. 7 is an a kind of schematic diagram of computer installation embodiment in the embodiment of the present application.
Specific embodiment
The embodiment of the present application provides a kind of processing method of heart rate signal, for realizing to the excellent of heart rate value detection process
Change.
In order to make those skilled in the art more fully understand application scheme, below in conjunction in the embodiment of the present application
Attached drawing, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described embodiment is only
The embodiment of the application a part, instead of all the embodiments.Based on the embodiment in the application, ordinary skill people
Member's every other embodiment obtained without making creative work, all should belong to the model of the application protection
It encloses.
The description and claims of this application and term " first ", " second ", " third ", " in above-mentioned attached drawing
Four " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way
Data be interchangeable under appropriate circumstances, so that the embodiments described herein can be in addition to illustrating herein or describing
Sequence other than appearance is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that covering is non-exclusive
Include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to clearly arrange
Those of out step or unit, but may include be not clearly listed or it is solid for these process, methods, product or equipment
The other step or units having.
Exercise is promotive of health, strong, and it is desirable to pass through the physical fitness of movement lifts oneself.The form of training,
Sporter can select the suitable movement of oneself according to the physical function state of oneself.Movement can monitor, and pass through prison
Heart rate measuring can guarantee that sporter will not cause movement effects unknown because of exercise intensity is too low to control the rhythm of motion exercise
It is aobvious, it is also possible to prevent that sporter's exercise intensity is excessive to lead to injury gained in sports.In the prior art, the common heart rate of intelligent wearable device
Detection mode is to detect volumetric blood variation in the tissue by photoelectric method to obtain pulse signal.Intelligent wearable device control
Illumination processed is mapped to tested tissue, after tissue resorption attenuation, from light source is ipsilateral or heteropleural is received by photelectric receiver,
Some structural constituents such as muscle, bone etc. are considered as the effect of light constant in tested tissue, and blood in arteries
Volume can change with the contraction of heart with diastole, cause blood to change the uptake of light, cause output intensity
Cyclically-varying is presented, it is photoplethysmography that the variation of final light intensity, which is converted to electric signal,
(photoplethysmography, PPG) signal.PPG signal processing operation to acquisition can acquire heart rate and blood oxygen saturation
The value of degree.
In existing frequently-used wearing product, common rate calculation be usually can from when domain lookup signal peak value,
After time-domain signal can also being converted into frequency-region signal, effective heart rate peak in frequency domain is then found out as last output.
Usually carrying out signal processing from frequency domain has preferable stability, smaller by noise jamming.But due to the limit of FFT itself
System, is to carry out calculating analysis to the signal of a segment length, therefore will appear the inadequate situation of signal resolution power.Such as when in signal
There are some jump signals, then will lead to entire frequency domain and get muddled, i.e., effective heart rate peak is completely covered.Now
Common processing method has wavelet transform process, but wavelet transformation computation complexity is high, has using reduction time-domain signal length
Method, but will lead to the case where overall signal is unstable and frequency domain resolving power reduces in this way.
However, using photoelectric sphyg wave carry out heart rate detection when, since signal is weaker, be highly prone to it is various bounce or not
The interference of stabilization signal leads to the unstability of overall signal since movement generates the influence of jump, often normal heart rate letter
It number is completely covered, to affect the normal calculating of heart rate.Situation more than being directed to, the present invention provides a kind of hearts rate
The processing method and its relevant device of signal for realizing quick identification mutation, and restore the method for original signal, not only
The complexity of algorithm is reduced, and solves the problems, such as spectral resolution.
In order to make it easy to understand, the detailed process in the embodiment of the present application is described below, referring to Fig. 1, the application
A kind of one embodiment of the processing method of heart rate signal includes: in embodiment
101, PPG signal is obtained；
In the present embodiment, terminal obtains PPG signal, wherein and the terminal can be the wearable devices such as bracelet, wrist-watch,
It can be other intelligent terminals such as mobile phone, plate, herein without limitation.
It is directly acquired pair in addition, the process for obtaining PPG signal can be by the PPG signal acquisition sensor that terminal carries
User is acquired to obtain, and is also possible to acquire by way of receiving the PGG signal that other terminals collect, may be used also
By be it is other in the form of obtain PPG signal, herein without limitation.
102, screening identification is carried out to the PPG signal, determines the jump signal in the PPG signal；
In the present embodiment, terminal carries out screening identification to the PPG signal that step 101 acquires, and determines the PPG signal
In jump signal.
Specifically, which carries out screening identification to the PPG signal that step 101 acquires, the mistake of screening identification
Journey can be the mean value for calculating the PPG signal first, then set a threshold value, traverse all PPG signals, will exceed the threshold
The PPG signal value of value is determined as jump signal；In addition it is also possible to be segmented first to the PPG signal, then segmentation asks equal
Value, the new number section of PPG that given threshold will exceed the threshold value later are determined as jump signal, can also be that other similar screenings are known
Other process, details are not described herein again.
103, mutation is carried out to the jump signal in the PPG signal to handle, obtain echo signal；
In the present embodiment, terminal carries out mutation to the jump signal in PPG signal and handles, and obtains echo signal.
Specifically, terminal can be believed in PPG signal for the mutation after step 102 determines jump signal
It number is handled, to obtain echo signal.Wherein, after this goes mutation processing to can be zero setting processing, delete jump signal
Distortion compensation is carried out, simulation recovery is carried out or other similar processing modes carries out mutation and handle, obtains target letter
Number.
104, heart rate value is calculated according to the echo signal.
In the present embodiment, obtained echo signal is handled according to step 103 and calculates heart rate value.
It specifically, that is, can further root after the PPG signal after going mutation to handle after obtaining echo signal
Heart rate value is calculated according to the echo signal.Such as can be calculated by time-domain analysis, that is, calculate PPG signal in certain time
Wave crest number, it is assumed that the time of continuous sampling 5 seconds, wave crest number in 5s is N, then heart rate is exactly N*12；It is another
It is that the bigger letter of opposite variation is found out in frequency domain signal diagrams by carrying out the characteristics of FFT transform obtains frequency domain to PPG signal
It number is exactly the AC signal of blood flow transformation, the corresponding frequency of the signal is f so heart rate HeartRate HR=f x60；
In addition to this it is possible to heart rate value be calculated by similar methods using the echo signal, specifically herein without limitation.
In the present embodiment, PPG signal is obtained；Screening identification is carried out to the PPG signal, is determined in the PPG signal
Jump signal；Mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal；According to described
Echo signal calculates heart rate value.Wherein, mutation is carried out to the jump signal in the PPG signal to handle, obtain mesh
Mark and then calculates heart rate value according to the echo signal at signal, to reduce jump signal in PPG signal to heart rate value
The optimization to heart rate value detection process is realized in the influence of calculating process.
In the embodiment of the present application, based on embodiment described in Fig. 1, screening knowledge specifically is carried out to the PPG signal in step 102
Not, during determining the jump signal in the PPG signal, the method that can take mean value for segmentation carries out screening and identifies
Out.Below referring to Fig. 2, based on embodiment described in Fig. 1, a kind of processing method of heart rate signal of the embodiment of the present application another
In embodiment, specifically may include: in step 102
201, the PPG signal is divided into N number of signal segment, N is the integer greater than 1；
In the present embodiment, the PPG signal that step 101 acquires is divided into N number of signal segment by terminal, and N is whole greater than 1
Number, herein, can be segmented according to precision, and segmentation is thinner, and resolving power is higher, and corresponding calculation amount will be bigger.
202, the average value for calculating separately N number of signal segment obtains N number of first mean value；
In the present embodiment, the respective mean value of N number of signal segment of step 201 division is calculated separately, N number of first mean value is obtained.
203, N number of first mean of mean is calculated, the second mean value is obtained；
In the present embodiment, terminal averages again to the respective mean value of N number of signal segment, i.e., carries out to N number of first mean value
It averages, obtains N number of first mean of mean i.e. second mean value.
204, determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is
The presupposition multiple of second mean value.
In the present embodiment, specific terminal can preset a presupposition multiple, to use the presupposition multiple and the second mean value
Solution obtains preset value, after this, then determines that in N number of signal segment, the signal segment greater than preset value is jump signal, i.e.,
The corresponding signal segment of numerical value found out in N number of first mean value more than the preset value is jump signal.
In the embodiment of the present application, based on embodiment described in Fig. 1 or Fig. 2, specifically to institute in the PPG signal described in step 103
State jump signal and carry out mutation and handle, obtain echo signal mode can there are many, such as zero setting processing, lose frame compensation etc.
Weakening treatment is carried out to jump signal.Below referring to Fig. 3, based on embodiment described in Fig. 1 or Fig. 2, the embodiment of the present application is a kind of
In another embodiment of the processing method of heart rate signal, specifically may include: in step 103
301, Weakening treatment is carried out to the jump signal, obtains the first signal；
In the present embodiment, terminal carries out Weakening treatment to the jump signal, obtains the first signal, first signal
Absolute value is less than the jump signal.
Specifically, terminal can carry out Weakening treatment to the jump signal, the first signal be obtained, so that the first signal is exhausted
The jump signal is less than to value, wherein can handle by zero setting, jump signal is deleted into laggard line distortion compensation, carry out
Simulation recovery or other similar processing modes carry out mutation processing and obtain after Weakening treatment to the jump signal
To the first signal.Herein it is preferably zero setting processing, i.e., jump signal is assigned a value of 0, it is possible to reduce the complexity of calculating, realization pair
Scheme advanced optimizes.
302, echo signal is obtained using jump signal described in first signal substituting in the PPG signal.
In the present embodiment, terminal, using jump signal described in first signal substituting, obtains mesh in the PPG signal
It marks signal and obtains echo signal herein by i.e. the first signal substituting of the jump signal after the Weakening treatment jump signal, it can
The influence for greatly eliminating jump signal generation shake guarantees the stability of signal to improve the signal-to-noise ratio of signal.
In the embodiment of the present application, based on embodiment described in Fig. 1 to Fig. 3, the mistake of PPG signal is specifically obtained in a step 101
Cheng Zhong can obtain PPG original signal in advance, then carry out preparatory optimization processing to the PPG original signal and obtain PPG letter again
Number.Below referring to Fig. 4, based on embodiment described in Fig. 1, Fig. 2 or Fig. 3, a kind of processing side of heart rate signal of the embodiment of the present application
In another embodiment of method, specific steps 101 may include:
401, PPG original signal is obtained；
In the present embodiment, terminal obtains PPG original signal, i.e., without the original signal for being pre-filtered processing, with step
101 is similar, can be the PPG original signal that terminal itself collects herein, is also possible to receive the PPG that other terminals are sent
Original signal, herein without limitation.
402, the PPG original signal is normalized, obtains second signal；
In the present embodiment, the PPG original signal is normalized in terminal, obtains second signal.Specifically,
PPG original signal is normalized, normalized is conducive to the control of dynamic range of signals.Wherein, normalization is
A kind of mode of simplified calculating, i.e., the expression formula that will have dimension are converted into nondimensional expression formula, become scale by transformation,
Handled within the scope of PPG original signal is mapped to 0~1, it is more convenient and quick, it should be grouped into Digital Signal Processing scope it
It is interior.
403, the second signal is filtered, obtains the PPG signal.
In the present embodiment, terminal is filtered the second signal, obtains the PPG signal, filtering herein
Processing specifically includes the either other adaptive filtering processes of high-pass filtering, low-pass filtering, to remove to a certain extent
The influence of interference signal, obtained PPG signal are PPG signal used in subsequent step, are realized to the further excellent of scheme
Change.
In the embodiment of the present application, based on embodiment described in Fig. 4, specifically the second signal is filtered in step 403
Wave processing, obtaining the PPG signal can be there are many filtering mode, such as high-pass filtering, low-pass filtering, adaptive-filtering etc..
The process will be specifically described below, referring to Fig. 5, based on embodiment described in Fig. 4, a kind of heart rate of the embodiment of the present application
In another embodiment of the processing method of signal, specific steps 403 may include:
501, high-pass filtering and/or low-pass filtering are carried out to the second signal, obtains the PPG signal；
In the present embodiment, terminal can carry out high-pass filtering and/or low-pass filtering to second signal, obtain the PPG letter
Number, i.e., high-pass filtering and/or low-pass filtering treatment are carried out to second signal.
Specifically, since PPG original signal is highly prone to the influence dimmed in movement or actual moving process, baseline
Variation is very big, therefore signal, that is, second signal after normalization is passed through high pass filter, processes,；In addition, low-pass filter can
To filter out the radio-frequency component in signal, to guarantee the flatness of signal.Illustratively, since the range of heart rate arrives for 0.5Hz
Between 4Hz, therefore, the cutoff frequency of high-pass filter can be set as 0.5Hz, according to the range of heart rate, by low-pass filter
Cutoff frequency can be set as 4Hz.
In the actual application of scheme, in addition to the corresponding high-pass filtering of step 501 and/or low-pass filtering treatment, also
There can be adaptive filtering processing, be described below by step 502 to step 503, it is worth mentioning at this point that, high pass is filtered herein
Wave and/or low-pass filtering treatment and the processing of subsequent adaptive filtering be also and/or relationship.
502, acceleration signal is obtained；
In the present embodiment, terminal obtains acceleration signal, and the acceleration signal exists with the PPG original signal to be associated with
Relationship.
503, adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
In the present embodiment, terminal carries out adaptive filtering to the second signal according to the acceleration signal, obtains institute
State PPG signal.
Specifically, since in PPG signal, there is motion artifacts ingredients, therefore the influence of motion artifacts is excluded, it can
By the way that motion artifacts ingredient acquired in acceleration transducer is filtered as reference signal using sef-adapting filter
It removes, to ensure that the stability of signal.Wherein, it is many not to be present in signal processing, control, image procossing etc. for adaptive-filtering
Same domain, it is a kind of more targeted filtering method of intelligence, commonly used in denoising, research object be have it is uncertain
System or information process.Here " uncertainty " refers to that the mathematical model of studied processing information process and its environment is not
It is completely specified, therefore, adaptive filtering operation is carried out to the second signal using acceleration signal, after obtaining PPG signal,
The further promotion to signal stabilization may be implemented, increase the accuracy that subsequent heart rate value calculates.
The embodiment of the present application is illustrated from the angle of method above, below from the angle of virtual bench to the application
Embodiment is introduced.
Referring to Fig. 6, a kind of one embodiment of the optimization system of heart rate signal detection includes: in the embodiment of the present application
Acquiring unit 601, for obtaining PPG signal；
Determination unit 602 determines the mutation letter in the PPG signal for carrying out screening identification to the PPG signal
Number；
Processing unit 603 handles for carrying out mutation to the jump signal in the PPG signal, obtains target
Signal；
Computing unit 604, for calculating heart rate value according to the echo signal.
In the present embodiment, acquiring unit 601, for obtaining PPG signal；Determination unit 602, for the PPG signal
Screening identification is carried out, determines the jump signal in the PPG signal；Processing unit 603, for right in the PPG signal
The jump signal carries out mutation and handles, and obtains echo signal；Computing unit 604, for being calculated according to the echo signal
Heart rate value out.Wherein, mutation is carried out to the jump signal in the PPG signal to handle, obtain echo signal and then
Heart rate value is calculated according to the echo signal, to reduce the jump signal in PPG signal to the shadow of heart rate value calculating process
It rings, realizes the optimization to heart rate value detection process.
As a preferred embodiment, the determination unit 602 is specifically used for:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described
The presupposition multiple of second mean value.
As a preferred embodiment, the processing unit 603 is specifically used for:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than described
Jump signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
As a preferred embodiment, the acquiring unit 601 is specifically used for:
Obtain PPG original signal；
The PPG original signal is normalized, second signal is obtained；
The second signal is filtered, the PPG signal is obtained.
As a preferred embodiment, the acquiring unit 601 is specifically used for:
High-pass filtering and/or low-pass filtering are carried out to the second signal, obtain the PPG signal；
And/or
Acceleration signal is obtained, there are incidence relations with the PPG original signal for the acceleration signal；
Adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
The embodiment of the present application is described from modular device above, below from the angle of hardware device to this Shen
Please the computer installation in embodiment be described, referring to Fig. 7, in the embodiment of the present application one of computer installation it is specific real
Applying example includes:
The device 700 can generate bigger difference because configuration or performance are different, may include one or more
Central processing unit (central processing units, CPU) 701 (for example, one or more processors) and storage
Device 705 is stored with one or more application program or data in the memory 705.
Wherein, memory 705 can be volatile storage or persistent storage.The program for being stored in memory 705 can wrap
One or more modules are included, each module may include to the series of instructions operation in server.Further, in
Central processor 701 can be set to communicate with memory 705, and a series of fingers in memory 705 are executed on intelligent terminal 700
Enable operation.
The device 700 can also include one or more power supplys 702, one or more wired or wireless networks
Interface 703, one or more input/output interfaces 704, and/or, one or more operating systems, such as
Windows ServerTM, Mac OS XTM, UnixTM, LinuxTM, FreeBSDTM etc..
The processor 701 is specifically used for executing following steps:
Obtain PPG signal；
Screening identification is carried out to the PPG signal, determines the jump signal in the PPG signal；
Mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal；
Heart rate value is calculated according to the echo signal.
It is described that screening identification is carried out to the PPG signal as a preferred embodiment, it determines in the PPG signal
Jump signal includes:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described
The presupposition multiple of second mean value.
As a preferred embodiment, mutation is carried out to the jump signal in the PPG signal and is handled, target is obtained
Signal includes:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than described
Jump signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
As a preferred embodiment, the acquisition PPG signal includes:
Obtain PPG original signal；
The PPG original signal is normalized, second signal is obtained；
The second signal is filtered, the PPG signal is obtained.
As a preferred embodiment, the second signal is filtered, obtaining the PPG signal includes:
High-pass filtering and/or low-pass filtering are carried out to the second signal, obtain the PPG signal；
And/or
Acceleration signal is obtained, there are incidence relations with the PPG original signal for the acceleration signal；
Adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
It is understood that the size of the serial number of above steps is not meant in the various embodiments of the application
Execution sequence it is successive, the execution of each step sequence should be determined by its function and internal logic, without coping with the embodiment of the present application
Implementation process constitute any restriction.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided herein, it should be understood that disclosed system, device and method can be with
It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit
It divides, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components
It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or
The mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, the indirect coupling of device or unit
It closes or communicates to connect, can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, each functional unit in each embodiment of the application can integrate in one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product
When, it can store in a computer readable storage medium.Based on this understanding, the technical solution of the application is substantially
The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words
It embodies, which is stored in a storage medium, including some instructions are used so that a computer
Equipment (can be personal computer, server or the network equipment etc.) executes the complete of each embodiment the method for the application
Portion or part steps.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only
Memory), random access memory (RAM, RandomAccess Memory), magnetic or disk etc. are various can store journey
The medium of sequence code.
The above, above embodiments are only to illustrate the technical solution of the application, rather than its limitations；Although referring to before
Embodiment is stated the application is described in detail, those skilled in the art should understand that: it still can be to preceding
Technical solution documented by each embodiment is stated to modify or equivalent replacement of some of the technical features；And these
It modifies or replaces, the spirit and scope of each embodiment technical solution of the application that it does not separate the essence of the corresponding technical solution.
Claims (10)
1. a kind of processing method of heart rate signal characterized by comprising
Obtain PPG signal；
Screening identification is carried out to the PPG signal, determines the jump signal in the PPG signal；
Mutation is carried out to the jump signal in the PPG signal to handle, and obtains echo signal；
Heart rate value is calculated according to the echo signal.
2. being determined the method according to claim 1, wherein described carry out screening identification to the PPG signal
Jump signal in the PPG signal includes:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described second
The presupposition multiple of mean value.
3. the method according to claim 1, wherein being gone in the PPG signal to the jump signal
Mutation processing, obtaining echo signal includes:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than the mutation
Signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
4. method according to any one of claims 1 to 3, which is characterized in that the acquisition PPG signal includes:
Obtain PPG original signal；
The PPG original signal is normalized, second signal is obtained；
The second signal is filtered, the PPG signal is obtained.
5. according to the method described in claim 4, obtaining described it is characterized in that, be filtered to the second signal
PPG signal includes:
High-pass filtering and/or low-pass filtering are carried out to the second signal, obtain the PPG signal；
And/or
Acceleration signal is obtained, there are incidence relations with the PPG original signal for the acceleration signal；
Adaptive filtering is carried out to the second signal according to the acceleration signal, obtains the PPG signal.
6. a kind of optimization system of heart rate signal detection characterized by comprising
Acquiring unit, for obtaining PPG signal；
Determination unit determines the jump signal in the PPG signal for carrying out screening identification to the PPG signal；
Processing unit handles for carrying out mutation to the jump signal in the PPG signal, obtains echo signal；
Computing unit, for calculating heart rate value according to the echo signal.
7. system according to claim 6, which is characterized in that the determination unit is specifically used for:
The PPG signal is divided into N number of signal segment, N is the integer greater than 1；
The average value for calculating separately N number of signal segment obtains N number of first mean value；
N number of first mean of mean is calculated, the second mean value is obtained；
Determine that the signal segment greater than preset value is jump signal in N number of signal segment, the preset value is described second
The presupposition multiple of mean value.
8. system according to claim 7, which is characterized in that the processing unit is specifically used for:
Weakening treatment is carried out to the jump signal, obtains the first signal, the absolute value of first signal is less than the mutation
Signal；
Using jump signal described in first signal substituting in the PPG signal, echo signal is obtained.
9. a kind of computer installation characterized by comprising
Processor, memory, input-output equipment and bus；
The processor, memory, input-output equipment are connected with the bus respectively；
The processor is for executing such as method described in any one of claim 1 to 5.
10. a kind of computer readable storage medium, is stored thereon with computer program, it is characterised in that: the computer program
It realizes when being executed by processor such as the step of any one of claims 1 to 5 the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811612483.5A CN109497984A (en) | 2018-12-27 | 2018-12-27 | A kind of processing method and relevant device of heart rate signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811612483.5A CN109497984A (en) | 2018-12-27 | 2018-12-27 | A kind of processing method and relevant device of heart rate signal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109497984A true CN109497984A (en) | 2019-03-22 |
Family
ID=65755397
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811612483.5A Pending CN109497984A (en) | 2018-12-27 | 2018-12-27 | A kind of processing method and relevant device of heart rate signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109497984A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112568886A (en) * | 2020-11-03 | 2021-03-30 | 中国科学院深圳先进技术研究院 | Detection method of heart rhythm, electronic device and computer readable storage medium |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110201946A1 (en) * | 2004-07-19 | 2011-08-18 | Turcott Robert G | Reducing data acquisition, power and processing for hemodynamic signal sampling |
| US8571642B2 (en) * | 2010-09-14 | 2013-10-29 | Pacesetter, Inc. | Pre-ejection interval (PEI) monitoring devices, systems and methods |
| CN106175742A (en) * | 2016-07-19 | 2016-12-07 | 北京心量科技有限公司 | A kind of heart sign acquisition methods and device |
| CN106933483A (en) * | 2017-02-28 | 2017-07-07 | 清华大学 | A kind of touch interactive mode that can perceive user's impression |
| CN107205640A (en) * | 2014-12-23 | 2017-09-26 | 日东电工株式会社 | Apparatus and method for removing the pseudomorphism in physiological measurement |
| CN107874750A (en) * | 2017-11-28 | 2018-04-06 | 华南理工大学 | Pulse frequency variability and the psychological pressure monitoring method and device of sleep quality fusion |
| CN108697329A (en) * | 2017-06-29 | 2018-10-23 | 华为技术有限公司 | The detection method and wearable device of wearable device |
| CN108828266A (en) * | 2018-08-30 | 2018-11-16 | 衡阳市衡山科学城科技创新研究院有限公司 | A kind of median filter method and device for accelerometer calibration |
-
2018
- 2018-12-27 CN CN201811612483.5A patent/CN109497984A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110201946A1 (en) * | 2004-07-19 | 2011-08-18 | Turcott Robert G | Reducing data acquisition, power and processing for hemodynamic signal sampling |
| US8571642B2 (en) * | 2010-09-14 | 2013-10-29 | Pacesetter, Inc. | Pre-ejection interval (PEI) monitoring devices, systems and methods |
| CN107205640A (en) * | 2014-12-23 | 2017-09-26 | 日东电工株式会社 | Apparatus and method for removing the pseudomorphism in physiological measurement |
| CN106175742A (en) * | 2016-07-19 | 2016-12-07 | 北京心量科技有限公司 | A kind of heart sign acquisition methods and device |
| CN106933483A (en) * | 2017-02-28 | 2017-07-07 | 清华大学 | A kind of touch interactive mode that can perceive user's impression |
| CN108697329A (en) * | 2017-06-29 | 2018-10-23 | 华为技术有限公司 | The detection method and wearable device of wearable device |
| CN107874750A (en) * | 2017-11-28 | 2018-04-06 | 华南理工大学 | Pulse frequency variability and the psychological pressure monitoring method and device of sleep quality fusion |
| CN108828266A (en) * | 2018-08-30 | 2018-11-16 | 衡阳市衡山科学城科技创新研究院有限公司 | A kind of median filter method and device for accelerometer calibration |
Non-Patent Citations (1)
| Title |
|---|
| 徐菁,朱敏: "基于三轴加速度信号的心率修正算法", 《计算机应用与软件》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112568886A (en) * | 2020-11-03 | 2021-03-30 | 中国科学院深圳先进技术研究院 | Detection method of heart rhythm, electronic device and computer readable storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107949321B (en) | Temporal interference removal and improved heart rate measurement tracking mechanism | |
| JP6646081B2 (en) | System and method for heart rate estimation | |
| Schäck et al. | Computationally efficient heart rate estimation during physical exercise using photoplethysmographic signals | |
| US20190059822A1 (en) | Heart rate detection method and apparatus, and electronic terminal thereof | |
| CN109222948B (en) | Method, device, electronic equipment and storage medium for eliminating motion interference noise | |
| US10285650B2 (en) | Heart monitoring device and method | |
| CN107708538B (en) | Steady heart rate estimation | |
| CN107273827B (en) | Electrocardiosignal R wave detection method and device | |
| CN108903929B (en) | Heart rate detection correction method, device, storage medium and system | |
| CN104706336B (en) | A kind of photo-electric pulse signal measuring method, device and measuring apparatus | |
| CN113349752A (en) | Wearable device real-time heart rate monitoring method based on sensing fusion | |
| JP6933220B2 (en) | Biometric information processing device, biometric information processing method and information processing device | |
| US20210100463A1 (en) | Motion detection and cancellation using ambient light | |
| CN109864713A (en) | Rhythm of the heart method based on multi-channel parallel filtering and spectral peak weighting selection algorithm | |
| CN109875543A (en) | For the heart rate estimation method and device under a variety of fitness exercise conditions of wearable heart rate monitor apparatus | |
| TWI505816B (en) | Detecting method and apparatus for blood oxygen saturation | |
| US9729693B1 (en) | Determining measurement confidence for data collected from sensors of a wearable device | |
| CN112957055A (en) | EEMD-PCA-based method and device for removing motion artifacts in EEG signal | |
| CN109620195B (en) | Multi-wavelength signal fusion heart rate detection method and system based on wearable equipment | |
| CN109497984A (en) | A kind of processing method and relevant device of heart rate signal | |
| CN114929101A (en) | System and method for physiological measurement based on optical data | |
| CN109620198B (en) | Cardiovascular index detection and model training method and device | |
| CN112292069B (en) | Method and device for screening atrial fibrillation | |
| CN115118788B (en) | Time sequence data compression method and device, wearable intelligent equipment and storage medium | |
| CN111643053A (en) | Method and system for reducing motion artifacts in pulse wave signals |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190628 Address after: 518057 Yuanzheng Industrial Park, 4012 North Fifth Avenue, Bantian Street, Longgang District, Shenzhen City, Guangdong Province Applicant after: LAUNCH TECH Co.,Ltd. Address before: 518000 Room 201, building A, No. 1, Qian Wan Road, Qianhai Shenzhen Hong Kong cooperation zone, Shenzhen, Guangdong (Shenzhen Qianhai business secretary Co., Ltd.) Applicant before: HESVIT HEALTH TECHNOLOGY Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190322 |
|
| RJ01 | Rejection of invention patent application after publication |