CN106175742A - A kind of heart sign acquisition methods and device - Google Patents

A kind of heart sign acquisition methods and device Download PDF

Info

Publication number
CN106175742A
CN106175742A CN201610571030.7A CN201610571030A CN106175742A CN 106175742 A CN106175742 A CN 106175742A CN 201610571030 A CN201610571030 A CN 201610571030A CN 106175742 A CN106175742 A CN 106175742A
Authority
CN
China
Prior art keywords
described
according
crest
ppg
data
Prior art date
Application number
CN201610571030.7A
Other languages
Chinese (zh)
Inventor
陈远翔
马斌
Original Assignee
北京心量科技有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 北京心量科技有限公司 filed Critical 北京心量科技有限公司
Priority to CN201610571030.7A priority Critical patent/CN106175742A/en
Publication of CN106175742A publication Critical patent/CN106175742A/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405Determining heart rate variability
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infra-red radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infra-red light
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/029Measuring or recording blood output from the heart, e.g. minute volume
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/026Measuring blood flow
    • A61B5/0295Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06KRECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K9/00Methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
    • G06K9/00496Recognising patterns in signals and combinations thereof
    • G06K9/00503Preprocessing, e.g. filtering
    • G06K9/0051Denoising
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06KRECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K9/00Methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
    • G06K9/00496Recognising patterns in signals and combinations thereof
    • G06K9/00503Preprocessing, e.g. filtering
    • G06K9/0051Denoising
    • G06K9/00516Denoising by applying a scale-space analysis, e.g. using wavelet analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06KRECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K9/00Methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
    • G06K9/00496Recognising patterns in signals and combinations thereof
    • G06K9/00523Feature extraction
    • G06K9/0053Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06KRECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K2209/00Indexing scheme relating to methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
    • G06K2209/05Recognition of patterns in medical or anatomical images
    • G06K2209/051Recognition of patterns in medical or anatomical images of internal organs

Abstract

The present invention relates to Computer Applied Technology field, in particular to a kind of heart sign acquisition methods and device.Obtain PPG data;Described PPG numerical value is done the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes high-frequency noise;The PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revises baseline drift;According to the described PPG data of correction baseline drift, obtain PPG waveform medium wave peak and the position of trough, and according to described crest and the position of described trough and the sample rate of described PPG data, obtain the time interval numerical value between crest;According to the time interval numerical value between described crest, draw RR interval data, and calculate heart sign according to described RR interval data.The method and device can allow user more fully obtain the heart sign numerical value of self.

Description

A kind of heart sign acquisition methods and device

Technical field

The present invention relates to Computer Applied Technology field, in particular to a kind of heart sign acquisition methods and dress Put.

Background technology

At present, except, in the detection that hospital, MEC etc. are proprietary, current user oneself is in the measurement one to heart sign As concentrate in the measurement to heart rate, and precision is low, it is impossible to obtained other heart sign, such as heart rate variability by heart rate Property, heart are concerned with index, heart amount output, spiritual degree at random and heart volume index etc..

Therefore, a kind of method that user can be allowed at any time to obtain more comprehensively heart sign becomes the most urgently to be resolved hurrily Problem.

Summary of the invention

In view of this, the purpose of the embodiment of the present invention is to provide a kind of heart sign acquisition methods and device, it is possible to User is allowed to obtain self more comprehensively heart sign.

First aspect, embodiments provides a kind of heart sign acquisition methods, including:

Obtain PPG data;

Described PPG numerical value is done the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes high-frequency noise;

The PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revises baseline drift;

According to revising the described PPG data of baseline drift, obtain PPG waveform medium wave peak and the position of trough, and according to Described crest and the position of described trough and the sample rate of described PPG data, obtain the time interval numerical value between crest;

According to the time interval numerical value between described crest, draw RR interval data, and according to described RR interval data meter Calculate heart sign.

In conjunction with first aspect, embodiments provide the first possible embodiment of first aspect, wherein:, The described described PPG data according to correction baseline drift, obtains PPG waveform medium wave peak and the position of trough, specifically includes:

Described PPG data is carried out second dervative process, obtains described crest and the position of described trough.

In conjunction with first aspect, embodiments provide the embodiment that the second of first aspect is possible, wherein: institute Before stating the time interval numerical value obtained between crest, also include:

According to the position at described crest place, obtain the height of described crest, and the height of described crest is carried out normalizing Change processes, and obtains normalized crest height numerical value;

Calculate the meansigma methods of all normalized described crest height numerical value;

Normalized described crest height numerical value is compared with described meansigma methods respectively, and according to the result of comparison, Ignore and exceed the Height Anomalies crest of described meansigma methods 50% with described meansigma methods difference.

In conjunction with first aspect, embodiments provide the third possible embodiment of first aspect, wherein: institute Before stating the time interval numerical value obtained between crest, also include:

According to the position at described crest place, calculate the time difference between two the most adjacent described crests;

Described time difference is compared with the threshold value preset;

If described time difference is less than described default threshold value, then ignore one of them among two described crests.

In conjunction with first aspect, embodiments provide the 4th kind of possible embodiment of first aspect, wherein: institute State and calculate heart sign according to described RR interval data, specifically include: calculate heart rate variability, heart phase according to RR interval data In dry index, heart amount output, spiritual degree at random, heart volume index at least one;

Wherein, when calculating heart rate variability rate according to RR interval data, specifically include:

According to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, draw HRV signal curve;

According to following formula calculating HRV Time-Domain Technique index:

NNVGR is the meansigma methods of phase between the most normal NN;

SDNN is the standard deviation of phase between the most normal NN;

RNSSD is the root-mean-square value that between whole adjacent NN, the phase is the poorest;

SDSD is the standard deviation of the difference of phase between whole adjacent NN;

NN50, for the heart rate more than 50ms of the difference of phase between whole adjacent NN;

PNN50 is the percentage of the difference of the phase heart rate more than 50ms between whole adjacent NN Ratio;

In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is the phase between i+1 RR Data.It it is the meansigma methods of the RR interval data of N number of heartbeat;Between NN, the phase is hole IBI;

RRi'=RRi-RRi-1

RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 ;

When according to RR interval data calculate heart be concerned with index time, specifically include:

Described RR interval data is done fast fourier transform, generates FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the general power in the range of the first predeterminated frequency;

According to described power spectrum, calculate peak value peak power in the range of the second predeterminated frequency;

By the ratio of described peak value peak power Yu described general power, calculate heart and be concerned with index;

When calculating the output of heart amount according to RR interval data, specifically include:

Described PPG waveform is done data segment cutting, obtains PPG Wave data section;

To PPG Wave data section each described: described PPG Wave data section to be done wavelet transformation, revise baseline drift; Calculate the meansigma methods of the PPG Wave data section medium wave peak height after revising baseline drift;Calculate RR crest in PPG Wave data section Ratio between height and the meansigma methods of crest height, and using ratio minima as waveform quality E, wherein 0 < E≤1;

Extract all PPG Wave data sections of E > 0.9;

From the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, obtain ripple Shape time T;

According to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein:

Obtain crest Ps, trough Pd, according to following formula calculating K value:

K = P m - P d P s - P d ;

And according to following formula, calculate heart amount output CO:

When calculating spiritual degree at random according to RR interval data, specifically include:

RR interval data is done fast fourier transform, it is thus achieved that FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the power P between the 3rd predeterminated frequency scope;

According to described power spectrum, calculate general power A between the 4th predeterminated frequency scope;

By the ratio of P Yu A, calculate spiritual degree at random;

When calculating heart volume index according to RR interval data, specifically include:

Described RR interval data is done fast fourier transform, it is thus achieved that FFT data;

Described FFT data is converted to power spectrum;

According to described power spectrum, calculate power P V between the 5th predeterminated frequency scope;

According to described power spectrum, calculate power P L between the 6th predeterminated frequency scope;

According to described power spectrum, calculate power P H between the 7th predeterminated frequency scope;

Heart rate HR is calculated according to RR interval data;

According to following formula calculating heart volume index HA:

H A = P H P V + P L 2 × H R .

In conjunction with first aspect, embodiments provide the 5th kind of possible embodiment of first aspect, wherein: institute State acquisition PPG data, specifically include:

Obtain the detection video of predetermined time period;

Following method is used to obtain the PPG numerical value of each two field picture in described detection video:

Obtain the red channel numerical value of each pixel in image, and the described red channel numerical value of all pixels is entered Row is cumulative, forms the PPG numerical value of described image;

Wherein, described detection video includes: finger detection video, cervical region detection video, forehead detection video, wrist detection Video, instep detect at least one in video.

Second aspect, the embodiment of the present invention also provides for a kind of heart sign acquisition device, including:

PPG data acquisition module, is used for obtaining PPG data;

Filtration module, for described PPG numerical value does the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes height Frequently noise;And the PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revise baseline drift;

Crest and wave trough position determine module, according to the described PPG data of correction baseline drift, obtain in PPG waveform Crest and the position of trough, and according to described crest and the position of described trough and the sample rate of described PPG data, obtain Take the time interval numerical value between crest;

Heart sign computing module, for according to the time interval numerical value between described crest, draws RR interval data, and Heart sign is calculated according to described RR interval data.

In conjunction with second aspect, embodiments provide the first possible embodiment of second aspect, wherein: also Including:

Crest height acquisition module, for the position according to described crest place, obtains the height of described crest, and to institute The height stating crest is normalized, and obtains normalized crest height numerical value;

Average wave peak heights Numerical Simulation Module, for calculating the average of all normalized described crest height numerical value Value;

Height Anomalies crest ignores module, for being entered with described meansigma methods respectively by normalized described crest height numerical value Row comparison, and according to the result of comparison, ignores and exceedes the Height Anomalies crest of described meansigma methods 50% with described meansigma methods difference.

In conjunction with second aspect, embodiments provide the embodiment that the second of second aspect is possible, wherein:

Also include:

Peak time difference acquisition module, for the position according to described crest place, calculates two the most adjacent described ripples Time difference between peak;

Comparing module, for comparing described time difference with the threshold value preset;

Time anomaly crest ignores module, for when described time difference is less than described default threshold value, ignoring two institutes State one of them among crest.

In conjunction with second aspect, embodiments provide the third possible embodiment of second aspect, wherein: institute State heart sign computing module to include: heart rate variability rate calculating sub module, heart be concerned with Index for Calculation submodule, the heart amount output meter In operator module, spiritual degree calculating sub module at random and heart volume index calculating sub module at least one;

Wherein, described heart rate variability rate calculating sub module is used for:

According to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, draw HRV signal curve;

According to following formula calculating HRV Time-Domain Technique index:

NNVGR is the meansigma methods of phase between the most normal NN;

SDNN is the standard deviation of phase between the most normal NN;

RNSSD is the root-mean-square value that between whole adjacent NN, the phase is the poorest;

SDSD is the standard deviation of the difference of phase between whole adjacent NN;

NN50, for the heart rate more than 50ms of the difference of phase between whole adjacent NN;

PNN50 is the percentage of the difference of the phase heart rate more than 50ms between whole adjacent NN Ratio;

In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is the phase between i+1 RR Data.It it is the meansigma methods of the RR interval data of N number of heartbeat;Between NN, the phase is hole IBI;

RRi'=RRi-RRi-1

RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 ;

The described heart Index for Calculation submodule that is concerned with is used for:

Described RR interval data is done fast fourier transform, generates FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the general power in the range of the first predeterminated frequency;

According to described power spectrum, calculate peak value peak power in the range of the second predeterminated frequency;

By the ratio of described peak value peak power Yu described general power, calculate heart and be concerned with index;

Described heart amount output calculating sub module is used for:

Described PPG waveform is done data segment cutting, obtains PPG Wave data section;

To PPG Wave data section each described: described PPG Wave data section to be done wavelet transformation, revise baseline drift; Calculate the meansigma methods of the PPG Wave data section medium wave peak height after revising baseline drift;Calculate RR crest in PPG Wave data section Ratio between height and the meansigma methods of crest height, and using ratio minima as waveform quality E, wherein 0 < E≤1;

Extract all PPG Wave data sections of E > 0.9;

From the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, obtain ripple Shape time T;

According to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein:

Obtain crest Ps, trough Pd, according to following formula calculating K value:

K = P m - P d P s - P d ;

And according to following formula, calculate heart amount output CO:

Described spirit degree at random calculating sub module is used for:

RR interval data is done fast fourier transform, it is thus achieved that FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the power P between the 3rd predeterminated frequency scope;

According to described power spectrum, calculate general power A between the 4th predeterminated frequency scope;

By the ratio of P Yu A, calculate spiritual degree at random;

Described heart volume index calculating sub module is used for: described RR interval data is done fast fourier transform, it is thus achieved that FFT number According to;

Described FFT data is converted to power spectrum;

According to described power spectrum, calculate power P V between the 5th predeterminated frequency scope;

According to described power spectrum, calculate power P L between the 6th predeterminated frequency scope;

According to described power spectrum, calculate power P H between the 7th predeterminated frequency scope;

Heart rate HR is calculated according to RR interval data;

According to following formula calculating heart volume index HA:

H A = P H P V + P L 2 × H R .

Heart sign that the embodiment of the present invention is provided obtains crime and device, will first obtain PPG data, and PPG data Obtain easily, after obtaining PPG data, first PPG data done the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, Remove high-frequency noise, then the PPG data removing high-frequency noise is carried out discrete wavelet Filtering Processing, revise baseline drift.So After, according to PPG data, obtain PPG waveform medium wave peak and the position of trough, and according to the sample rate of PPG data, obtain crest Between time interval numerical value, finally according to this time interval numerical value, calculate between RR the phase, and finally calculate heart according to the phase between RR Sign.In this process, computational accuracy is high, and error is little, and user can be allowed to obtain more comprehensively heart sign at any time Numerical value.

For making the above-mentioned purpose of the present invention, feature and advantage to become apparent, preferred embodiment cited below particularly, and coordinate Appended accompanying drawing, is described in detail below.

Accompanying drawing explanation

In order to be illustrated more clearly that the technical scheme of the embodiment of the present invention, below by embodiment required use attached Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, and it is right to be therefore not construed as The restriction of scope, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to according to this A little accompanying drawings obtain other relevant accompanying drawings.

Fig. 1 shows the flow chart of a kind of heart sign acquisition methods that the embodiment of the present invention provided;

Fig. 2 shows in the heart sign acquisition methods that the embodiment of the present invention is provided, and obtains the concrete grammar of PPG data Flow chart;

Fig. 3 shows the flow chart of the another kind of heart sign acquisition methods that the embodiment of the present invention provided;

Fig. 4 shows the flow chart of the another kind of heart sign acquisition methods that the embodiment of the present invention provided;

Fig. 5 shows the flow chart of the another kind of heart sign acquisition methods that the embodiment of the present invention provided;

Fig. 6 shows the structural representation of a kind of heart sign acquisition device that the embodiment of the present invention provided;

Fig. 7 shows the structural representation of the another kind of heart sign acquisition device that the embodiment of the present invention provided.

Detailed description of the invention

For making the purpose of the embodiment of the present invention, technical scheme and advantage clearer, below in conjunction with the embodiment of the present invention Middle accompanying drawing, is clearly and completely described the technical scheme in the embodiment of the present invention, it is clear that described embodiment is only It is a part of embodiment of the present invention rather than whole embodiments.Generally real with the present invention illustrated described in accompanying drawing herein The assembly executing example can be arranged with various different configurations and design.Therefore, below to the present invention's provided in the accompanying drawings The detailed description of embodiment is not intended to limit the scope of claimed invention, but is merely representative of the selected reality of the present invention Execute example.Based on embodiments of the invention, the institute that those skilled in the art are obtained on the premise of not making creative work There are other embodiments, broadly fall into the scope of protection of the invention.

Current heart sign acquisition methods carries out detecting generally by the detection equipment of hospital or health check-up mechanism, Making user the heart sign of self cannot be carried out real-time monitoring, based on this, a kind of heart sign that the application provides obtains Access method, device and system, can allow user be obtained from the dirty sign of body and mind in real time.

For ease of the present embodiment is understood, first to a kind of heart sign acquisition side disclosed in the embodiment of the present invention Method describes in detail.

Shown in Figure 1, the heart sign acquisition methods that the embodiment of the present invention is provided includes:

S101: obtain PPG data;

Implementing when, the acquisition methods of PPG data has multiple, for example, it is possible to directly survey from proprietary PPG data PPG data is directly obtained, additionally it is possible to use other mode to obtain PPG data on amount equipment.

Seeing described in Fig. 2, the embodiment of the present invention provides a kind of concrete grammar obtaining PPG data, including:

S201: obtain the detection video of predetermined time period.

Implementing when, the when of flowing in the blood vessel due to the blood of human body, blood vessel, especially tremulous pulse, meeting Along with the expansion of occurrence law of beating and the contraction of heart, and due to the expansion of blood vessel and contraction, cause covering on blood vessel The skin of side also there will be corresponding fluctuations.And in the case of light stability, the shadow on skin can be along with skin Rising and falling and produce certain change, this change is difficult to observe at human eye, but after being taken into video, can pass through Video is processed, quantify after obtain, accordingly, it would be desirable to obtain detection video.This detection video is that the detection for human body regards Frequently, and in order to amplify skin to greatest extent with the expansion of blood vessel and to shrink produced change, acquired video one As jump more visibility point, such as finger, cervical region, forehead, instep, wrist etc., therefore, this detection for body surface blood vessel Video may include that finger detection video, cervical region detection video, forehead detection video, wrist detection video, instep detection video In at least one.And for the detection video of different health different parts, its particular location obtained also is different.Example As, if detection video is finger detects video, then with the finger detection video of the joint finger pulp loss side near palm For most preferably;If cervical region detection video, then the cervical region detection video above carotid artery is optimal.In addition, it is necessary to note Be, it is also possible to obtain the detection video at other positions of health.

It addition, obtain detection video when, it is possible to use mobile phone camera obtains, it is also possible to by other Image acquisition equipment obtains.Such as, using mobile phone camera to carry out Image Acquisition when, finger is attached to photographic head Camera lens on, and shooting when open flash lamp, captured image is a red image by eye-observation.And Along with blood vessel contraction and expansion when, detection video in red channel numerical value can be As time goes at different images The change that middle appearance is different.

And in order to ensure the precision of acquired heart amount sign, be unlikely to again follow-up, detection video to be analyzed simultaneously During devote a tremendous amount of time calculating therefore, acquired image resolution ratio can have according to concrete practical situation The setting of body.If it is intended to obtain the PPG data that precision is higher, then the image resolution ratio obtained will be high accordingly.Need If it should be noted that the resolution of image is bigger in acquired detection video, then need to intercept a part wherein and come It is analyzed calculating.Usually, it addition, PPG data is obtaining when, there is sample rate.For example, it is possible to adopting detection video Sample rate is as the sample rate of PPG data.Such as, obtaining finger detection video when, keep finger static, taken the photograph by mobile phone As head gathers 640 × 480 resolution, the finger detection video of 30HZ sample rate.

Detection video, obtaining when, is the detection video obtaining predetermined time period, in the shooting of detection video Cheng Zhong, then captured video time length can intercept one section of bat with this time span preset of heavy rain from this video Take the photograph the detection video of relatively stable predetermined time period.

S202: use following method to obtain the PPG numerical value of each two field picture in described detection video: to obtain in image each The red channel numerical value of pixel, and the described red channel numerical value of all pixels is added up, form described image PPG numerical value.

Implementing when, the when of there is fluctuations due to human body skin with the expansion of blood vessel and contraction, institute The shadow change caused is that comparison is small, and therefore, detection video needs to keep certain duration, and (this duration can be according to reality Detection case carry out concrete setting), and each two field picture in the detection video of this certain time length is all processed, then This change is amplified.Specifically, after obtaining detection video, RGB can be carried out for each two field picture in detection video Color Channel is analyzed.Carrying out Color Channel analysis when, can be by the red channel numerical analysis of pixel each on image Out, then the red channel numerical value of pixels all in this image is added up, form the PPG of this two field picture (PhotoPlethysmoGraphy, photoplethysmographic) numerical value.With the time as transverse axis, with the size of PPG numerical value as the longitudinal axis Setting up coordinate system, the line of all points representing PPG numerical value under this coordinate system is a fluctuation with crest and trough Figure.

With resolution for 640 × 480, as a example by the finger that sample rate is 30HZ, predetermined time period is 10 seconds detects video: Resolution be 640 × 480 image have 307200 pixels, for each two field picture, obtaining, this two field picture is all of During the red channel numerical value of pixel, it is possible to obtain 307200 red channel numerical value;Then, to these 307300 red channels Numerical value adds up, the final PPG numerical value obtaining this two field picture.And the sample rate of this detection video is 30HZ, be equivalent to per second Clock gathers 30 images, within 10 seconds, can gather 300 images, the final PPG numerical value that PPG numerical value is 300 images obtained, I.e. have 300 PPG numerical value altogether.

In this process, image is generally by different Color Channels, and the form of image is different, and the color of pie graph picture is led to Road is also the most different, and owing to the blood of human body is red, in acquired detection video, redness is the most prominent, therefore use R, The image that G, B triple channel is constituted is analyzed, same be analyzed for each two field picture in detection video time Wait, be the red channel numerical value directly obtaining each pixel of image.

S203: using all images PPG numerical value as described PPG data.

PPG data is Wave data, and forms PPG waveform.

S102: described PPG numerical value does the low-pass filtering treatment being less than or equal to 4HZ by frequency, removes high-frequency noise.

S103: the PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revises baseline drift.

Implementing when, owing to people are using image acquisition equipment shot detection video when, it is difficult to protect The position that card is taken is totally stationary motionless, is likely to owing to human body is conscious or cannot control in the middle of the process of shooting The shake of system is moved, and causes there will be during some Image Acquisition in video bigger deviation, therefore, at acquired PPG numerical value In be likely to there is bigger error, affect final result of calculation, therefore, in embodiments of the present invention, according to PPG Before numerical value obtains the position of crest and trough, in addition it is also necessary to be filtered PPG numerical value processing.Specifically it is being filtered process When, to first use FIR filter that PPG numerical value does the low-pass filtering treatment being less than or equal to 4HZ by frequency, remove high frequency Noise, then re-uses 6 layer scattering wavelet filters and the PPG fluctuation data eliminating high-frequency noise is carried out discrete wavelet filtering Process, revise baseline drift, the accurate PPG data of final acquisition.

S104: according to the described PPG data of correction baseline drift, obtains PPG waveform medium wave peak and the position of trough, and According to described crest and the position of described trough and the sample rate of described PPG data, obtain the time interval between crest Numerical value;

Implementing when, the when of obtaining the position of the crest in PPG waveform and trough, it is by all images PPG numerical value utilize 2 order derivatives, find the position of crest and trough, concrete, be by composition PPG data multiple PPG numbers In value, the PPG numerical value of consecutive frame carries out additive operation, and the result subtracted each other is done normalized, such as, if PPG numerical value There are 10, respectively ppg1, ppg2, ppg3 ..., ppg9, ppg10, carrying out additive operation when, ppg2-ppg2, Ppg2-ppg3 ... ppg9-ppg8, ppg10-ppg9.The like;Adjacent PPG numerical value is carried out additive operation it After, operation result is normalized, the first derivative numerical value obtained includes :+1,0 and-1, this mediation number numerical tabular Having reached blood when trickle in the blood vessel, the blood in blood vessel is in increase trend and still reduces trend, the most all of PPG The PPG waveform that numerical value is formed is in propradation and is in decline state.Then, then by obtained first derivative numerical value In adjacent numerical value carry out additive operation, and the result subtracted each other is done again normalized, finally can obtain second dervative Numerical value also includes :+1,0 and-1.And now ,+1 represents crest ,-1 represents trough, 0 PPG represented between crest and trough Numerical value.In the middle of it should be noted that the first derivative numerical value obtained by the process to PPG numerical value and second dervative numerical value are only Calculating process.Furthermore it is also possible to go to judge the position of crest and trough according to overall increase, the trend of reduction of PPG numerical value, And the judgement of crest and wave trough position is not limited solely to the one provided in above-described embodiment.Such as can also be with relative method PPG numerical value is contrasted respectively, thus judges the position of crest and trough.

It addition, crest and the position of trough, corresponding to the position at place on the time shaft of PPG numerical value,.And according to sampling Rate and the position at PPG numerical value place, it is possible to calculate the time difference when obtaining of the PPG numerical value corresponding to adjacent crest, Obtain the time interval numerical value between crest, the time interval numerical value between adjacent peaks, i.e. illustrate adjacent twice jumping of heart Time difference between Dong.

S105: according to the time interval numerical value between described crest, draws RR interval data, and according to issue between described RR According to calculating heart sign.

After RR interval data obtains, just can calculate heart sign according to RR interval data.

The heart features acquisition methods that the embodiment of the present invention is provided, will first obtain PPG data, and PPG data obtains and holds Easily, after obtaining PPG data, first PPG data is done the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes height Frequently noise, then carries out discrete wavelet Filtering Processing to the PPG data removing high-frequency noise, revises baseline drift.Then, according to PPG data, obtains PPG waveform medium wave peak and the position of trough, and according to the sample rate of PPG data, obtains between crest Time interval numerical value, finally according to this time interval numerical value, calculates the phase between RR, and finally calculates heart sign according to the phase between RR. In this process, computational accuracy is high, and error is little, and user can be allowed to obtain more comprehensively heart sign numerical value at any time.

Shown in Figure 3, the embodiment of the present invention thus provide heart sign acquisition methods, in above-mentioned several embodiments On the basis of, before the time interval numerical value between described acquisition crest, also include:

S301: according to described crest and the position at trough place, obtains the height of described crest, and to described crest Highly it is normalized, obtains normalized crest height numerical value.

Implementing when, with the time as transverse axis, in the coordinate system that PPG numerical value is set up by the longitudinal axis, the height of crest Degree is: crest is done vertical line, and the line of this crest both sides trough forms intersection point with this vertical line, and this vertical line is positioned at crest and intersection point Between distance.

S302: calculate the meansigma methods of all normalized described crest height numerical value;

S303: normalized described crest height numerical value is compared with described meansigma methods respectively, and according to comparison As a result, ignore and exceed the Height Anomalies crest of described meansigma methods 50% with described meansigma methods difference.

Implementing when, after the height of crest is normalized, the PPG that can be more smoothed The curve of cyclical fluctuations, is more beneficial for follow-up calculating.Calculate the meansigma methods of all normalized crest height numerical value, and return with each One changes crest height numerical value compares with this meansigma methods, if normalization crest height numerical value is more than or equal to the 50% of meansigma methods Time, it is meant that this crest is probably abnormal crest, therefore to be ignored by this crest.Will PPG data corresponding to this crest not It is re-used as crest to treat.

Shown in Figure 4, in the heart sign acquisition methods that the embodiment of the present invention is provided, between described acquisition crest Before time interval numerical value, also include:

S401: according to the position at described crest place, calculate the time difference between two the most adjacent described crests;

S402: described time difference is compared with the threshold value preset;

S403: if described time difference is less than described default threshold value, then the camber ignoring two described crests is less Crest.

Implementing when, after the crest that have ignored Height Anomalies, can be according to the position at remaining crest place Put, calculate the time difference between two the most adjacent crests, described time difference is contrasted with the threshold value preset, if the time Difference is less than the threshold value preset, then it is assumed that wherein having crest is interference crest, and non-cardiac normally beat in the case of produce, because of This ignores in two crests.Which is specifically ignored, concrete judgement can be carried out according to actual situation, such as, Having three adjacent peaks, wherein, the time difference between middle crest and other two crests is respectively less than the threshold value preset, So it is ignored as being positioned at the crest of centre.

The threshold value preset can carry out concrete setting according to the sample rate of actual detection video, and sample rate is the highest, in advance If threshold value the least.

Shown in Figure 5, in the acquisition methods of the heart sign that the embodiment of the present invention is provided, described according to described RR between Issue, according to calculating heart sign, specifically includes:

S501: calculate heart rate variability according to RR interval data, heart is concerned with index, heart amount output, spiritual journey at random Degree, in heart volume index at least one;

It addition, the embodiment of the present invention also provide for by RR interval data calculate heart sign concrete grammar:

One, by RR interval data calculating HRV (heart rate variability):

1, according to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, HRV signal curve is drawn.

2, according to following formula calculating HRV Time-Domain Technique index:

All meansigma methodss (unit: ms) of phase between normal sinus heartbeat (NN);

The standard deviation (unit: ms) of phase between the most normal NN;

The root-mean-square value (unit: ms) of the difference of phase between whole adjacent NN

The standard deviation (unit: ms) of the difference of phase between whole adjacent NN;

The difference of the phase heart rate more than 50ms between NN50 whole adjacent NN;

The percentage ratio of the difference of the phase heart rate more than 50ms between whole adjacent NN.

In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is the phase between i+1 RR Data.It it is the meansigma methods of the RR interval data of N number of heartbeat;

RRi'=RRi-RRi-1

RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 .

Two: calculate MIND (heart be concerned with index) by RR interval data:

1, described RR interval data is done fast fourier transform, generate FFT data;

2, FFT data is converted to power spectrum;

3, according to described power spectrum, calculate the general power in the range of the first predeterminated frequency and (such as, calculate 0.003HZ extremely 0.4HZ general power);

4, according to described power spectrum, in the range of calculating the second predeterminated frequency, peak value peak power (such as calculates 0.04HZ extremely Maximum peak power between 0.26HZ, peak value window 0.03HZ.)

5, by the ratio of the peak value peak power in the range of predeterminated frequency Yu general power, MIND is drawn.

Three: by RR interval data calculating CO (cardiac output):

1, described PPG waveform is done data segment cutting, obtain PPG Wave data section, for example, it is possible to according to 5 seconds windows Mouthful, within 2 seconds, once go forward one by one and do data segment cutting.

2, PPG Wave data section each described is done as follows:

1., to described PPG Wave data section do wavelet transformation, revise baseline drift;

2. the meansigma methods of the PPG Wave data section medium wave peak height AVE after revising baseline drift, is calculated;

3., the ratio between RR crest height and AVE in PPG Wave data section is calculated, and using ratio minima as ripple Form quality amount E (0 < E≤1);(little value is divided by big value);

3, all PPG Wave data sections of E > 0.9 are extracted;

4, from the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, obtain Waveform time T;

5, according to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein,

Obtain crest Ps, trough Pd, according to following formula, calculate K value:

K = P m - P d P s - P d ;

6, according to following formula, calculating heart amount output CO:

Four: according to RR interval data calculating VIK (spiritual degree at random):

1, RR interval data is done fast fourier transform, it is thus achieved that FFT data;

2, FFT data is converted to power spectrum;

3, according to power spectrum, calculate the power P between the 3rd predeterminated frequency scope and (such as, calculate between 0HZ to 0.04HZ Power);

4, according to power spectrum, calculate general power A between the 4th predeterminated frequency scope (such as, calculate 0HZ to 0.4HZ it Between general power);

5, by the ratio of P Yu A, VIK is drawn.

Five, according to RR interval data, calculating HA (heart volume index):

1, described RR interval data is done fast fourier transform, it is thus achieved that FFT data;

2, described FFT data is converted to power spectrum;

3, according to described power spectrum, power P V calculated between the 5th predeterminated frequency scope (such as calculates 0.0033HZ extremely Power P V between 0.04HZ);

4, according to described power spectrum, power P L calculated between the 6th predeterminated frequency scope (such as calculates 0.04HZ extremely Power P L between 0.15HZ);

5, according to described power spectrum, power P H calculated between the 7th predeterminated frequency scope (such as calculates 0.15HZ extremely Power P H between 0.4HZ);

Wherein, the 5th predeterminated frequency, the 6th predeterminated frequency and the 7th predeterminated frequency are continuous print frequencies.

6, heart rate HR is calculated according to RR interval data;

7, according to following formula calculating HA:

H A = P H P V + P L 2 × H R .

Further embodiment of this invention also provides for a kind of heart sign acquisition device, shown in Figure 6, embodiment of the present invention institute The heart sign acquisition device provided includes:

PPG data acquisition module, is used for obtaining PPG data;

Filtration module, for described PPG numerical value does the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes height Frequently noise;And the PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revise baseline drift;

Crest and wave trough position determine module, according to the described PPG data of correction baseline drift, obtain in PPG waveform Crest and the position of trough, and according to described crest and the position of described trough and the sample rate of described PPG data, obtain Take the time interval numerical value between crest;

Heart sign computing module, for according to the time interval numerical value between described crest, draws RR interval data, and Heart sign is calculated according to described RR interval data.

In the present embodiment, PPG data acquisition module, filtration module, crest and wave trough position determine module and heart body Levy concrete function and the interactive mode of computing module, can be found in the record of embodiment corresponding for Fig. 1, do not repeat them here.

The heart sign acquisition device that the embodiment of the present invention is provided, will first obtain PPG data, and PPG data obtains and holds Easily, after obtaining PPG data, first PPG data is done the cut-off frequency low-pass filtering treatment less than Deng drum 4HZ, removes height Frequently noise, then carries out discrete wavelet Filtering Processing to the PPG data removing high-frequency noise, revises baseline drift.Then, according to PPG data, obtains PPG waveform medium wave peak and the position of trough, and according to the sample rate of PPG data, obtains between crest Time interval numerical value, finally according to this time interval numerical value, calculates the phase between RR, and finally calculates heart sign according to the phase between RR. In this process, computational accuracy is high, and error is little, and user can be allowed to obtain more comprehensively heart sign numerical value at any time.

Shown in Figure 7, further embodiment of this invention also provides for another kind of heart sign acquisition device, also includes:

Crest height acquisition module, for the position according to described crest place, obtains the height of described crest, and to institute The height stating crest is normalized, and obtains normalized crest height numerical value;

Average wave peak heights Numerical Simulation Module, for calculating the average of all normalized described crest height numerical value Value;

Height Anomalies crest ignores module, for being entered with described meansigma methods respectively by normalized described crest height numerical value Row comparison, and according to the result of comparison, ignores and exceedes the Height Anomalies crest of described meansigma methods 50% with described meansigma methods difference.

Also include:

Peak time difference acquisition module, for the position according to described crest place, calculates two the most adjacent described ripples Time difference between peak;

Comparing module, for comparing described time difference with the threshold value preset;

Time anomaly crest ignores module, for when described time difference is less than described default threshold value, ignoring two institutes State one of them among crest.

In the present embodiment, crest height acquisition module, average wave peak heights Numerical Simulation Module, Height Anomalies crest are ignored Module, peak time difference acquisition module, comparing module and time anomaly crest ignore concrete function and the interactive mode of module, can See the record of embodiment corresponding for 3-Fig. 4, do not repeat them here.

Shown in Figure 7, further embodiment of this invention also provides for another kind of heart sign acquisition device, also includes:

Described heart sign computing module includes: heart rate variability rate calculating sub module, heart be concerned with Index for Calculation submodule, Heart amount exports at least one in calculating sub module, spiritual degree calculating sub module at random and heart volume index calculating sub module;

Wherein, described heart rate variability rate calculating sub module is used for:

According to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, draw HRV signal curve;

According to following formula calculating HRV Time-Domain Technique index:

NNVGR is the meansigma methods of phase between the most normal NN;

SDNN is the standard deviation of phase between the most normal NN;

RNSSD is the root-mean-square value that between whole adjacent NN, the phase is the poorest;

SDSD is the standard deviation of the difference of phase between whole adjacent NN;

NN50, for the heart rate more than 50ms of the difference of phase between whole adjacent NN;

PNN50 is the percentage of the difference of the phase heart rate more than 50ms between whole adjacent NN Ratio;

In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is the phase between i+1 RR Data.It it is the meansigma methods of the RR interval data of N number of heartbeat;Between NN, the phase is hole IBI;

RRi'=RRi-RRi-1

RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 ;

The described heart Index for Calculation submodule that is concerned with is used for:

Described RR interval data is done fast fourier transform, generates FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the general power in the range of the first predeterminated frequency;

According to described power spectrum, calculate peak value peak power in the range of the second predeterminated frequency;

By the ratio of described peak value peak power Yu described general power, calculate heart and be concerned with index;

Described heart amount output calculating sub module is used for:

Described PPG waveform is done data segment cutting, obtains PPG Wave data section;

To PPG Wave data section each described: described PPG Wave data section to be done wavelet transformation, revise baseline drift; Calculate the meansigma methods of the PPG Wave data section medium wave peak height after revising baseline drift;Calculate RR crest in PPG Wave data section Ratio between height and the meansigma methods of crest height, and using ratio minima as waveform quality E, wherein 0 < E≤1;

Extract all PPG Wave data sections of E > 0.9;

From the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, obtain ripple Shape time T;

According to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein:

Obtain crest Ps, trough Pd, according to following formula calculating K value:

K = P m - P d P s - P d ;

And according to following formula, calculate heart amount output CO:

Described spirit degree at random calculating sub module is used for:

RR interval data is done fast fourier transform, it is thus achieved that FFT data;

FFT data is converted to power spectrum;

According to described power spectrum, calculate the power P between the 3rd predeterminated frequency scope;

According to described power spectrum, calculate general power A between the 4th predeterminated frequency scope;

By the ratio of P Yu A, calculate spiritual degree at random;

Described heart volume index calculating sub module is used for: described RR interval data is done fast fourier transform, it is thus achieved that FFT number According to;

Described FFT data is converted to power spectrum;

According to described power spectrum, calculate power P V between the 5th predeterminated frequency scope;

According to described power spectrum, calculate power P L between the 6th predeterminated frequency scope;

According to described power spectrum, calculate power P H between the 7th predeterminated frequency scope;

Heart rate HR is calculated according to RR interval data;

According to following formula calculating heart volume index HA:

H A = P H P V + P L 2 × H R .

Wherein, be concerned with Index for Calculation submodule, the heart amount output of calculating of heart rate variability rate calculating sub module, heart calculates submodule Block, spiritual degree calculating sub module at random and heart volume index calculating sub module concrete function and interactive mode, can be found in 5 corresponding The record of embodiment, do not repeat them here.

Heart sign acquisition methods that the embodiment of the present invention is provided and the computer program of device, including storing The computer-readable recording medium of program code, the instruction that described program code includes can be used for performing in previous methods embodiment Described method, implements and can be found in embodiment of the method, does not repeats them here.

Those skilled in the art is it can be understood that arrive, for convenience and simplicity of description, and the system of foregoing description With the specific works process of device, it is referred to the corresponding process in preceding method embodiment, does not repeats them here.

If described function is using the form realization of SFU software functional unit and as independent production marketing or use, permissible It is stored in a computer read/write memory medium.Based on such understanding, technical scheme is the most in other words The part contributing prior art or the part of this technical scheme can embody with the form of software product, this meter Calculation machine software product is stored in a storage medium, including some instructions with so that a computer equipment (can be individual People's computer, server, or the network equipment etc.) perform all or part of step of method described in each embodiment of the present invention. And aforesaid storage medium includes: USB flash disk, portable hard drive, read only memory (ROM, Read-Only Memory), random access memory are deposited The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.

The above, the only detailed description of the invention of the present invention, but protection scope of the present invention is not limited thereto, and any Those familiar with the art, in the technical scope that the invention discloses, can readily occur in change or replace, should contain Cover within protection scope of the present invention.Therefore, protection scope of the present invention should described be as the criterion with scope of the claims.

Claims (10)

1. a heart sign acquisition methods, it is characterised in that including:
Obtain PPG data;
Described PPG numerical value is done the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes high-frequency noise;
The PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revises baseline drift;
According to the described PPG data of correction baseline drift, obtain PPG waveform medium wave peak and the position of trough, and according to described Crest and the position of described trough and the sample rate of described PPG data, obtain the time interval numerical value between crest;
According to the time interval numerical value between described crest, draw RR interval data, and calculate the heart according to described RR interval data Dirty sign.
Heart sign acquisition methods the most according to claim 1, it is characterised in that described according to the institute revising baseline drift State PPG data, obtain PPG waveform medium wave peak and the position of trough, specifically include:
Described PPG data is carried out second dervative process, obtains described crest and the position of described trough.
Heart sign acquisition methods the most according to claim 2, it is characterised in that between the time between described acquisition crest Before numerical value, also include:
According to the position at described crest place, obtain the height of described crest, and the height of described crest is normalized place Reason, obtains normalized crest height numerical value;
Calculate the meansigma methods of all normalized described crest height numerical value;
Normalized described crest height numerical value is compared with described meansigma methods respectively, and according to the result of comparison, ignores With the Height Anomalies crest that described meansigma methods difference exceedes described meansigma methods 50%.
Heart sign acquisition methods the most according to claim 3, it is characterised in that between the time between described acquisition crest Before numerical value, also include:
According to the position at described crest place, calculate the time difference between two the most adjacent described crests;
Described time difference is compared with the threshold value preset;
If described time difference is less than described default threshold value, then ignore one of them among two described crests.
5. according to the heart sign acquisition methods described in claim 1-4 any one, it is characterised in that described according to described RR Interval data calculates heart sign, specifically includes: calculate heart rate variability according to RR interval data, heart is concerned with index, heart amount In output, spiritual degree at random, heart volume index at least one;
Wherein, when calculating heart rate variability rate according to RR interval data, specifically include:
According to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, draw HRV signal curve;
According to following formula calculating HRV Time-Domain Technique index:
NNVGR is the meansigma methods of phase between the most normal NN;
SDNN is the standard deviation of phase between the most normal NN;
RNSSD is the root-mean-square value of the difference of phase between whole adjacent NN;
SDSD is the standard deviation of the difference of phase between whole adjacent NN;
NN50, for the heart rate more than 50ms of the difference of phase between whole adjacent NN;
PNN50 is the percentage ratio of the difference of the phase heart rate more than 50ms between whole adjacent NN;
In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is i+1 RR interval data.It it is the meansigma methods of the RR interval data of N number of heartbeat;Between NN, the phase is hole IBI;
RR′i=RRi-RRi-1
RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 ;
When according to RR interval data calculate heart be concerned with index time, specifically include:
Described RR interval data is done fast fourier transform, generates FFT data;
FFT data is converted to power spectrum;
According to described power spectrum, calculate the general power in the range of the first predeterminated frequency;
According to described power spectrum, calculate peak value peak power in the range of the second predeterminated frequency;
By the ratio of described peak value peak power Yu described general power, calculate heart and be concerned with index;
When calculating the output of heart amount according to RR interval data, specifically include:
Described PPG waveform is done data segment cutting, obtains PPG Wave data section;
To PPG Wave data section each described: described PPG Wave data section to be done wavelet transformation, revise baseline drift;Calculate Revise the meansigma methods of the PPG Wave data section medium wave peak height after baseline drift;Calculate RR crest height in PPG Wave data section And the ratio between the meansigma methods of crest height, and using ratio minima as waveform quality E, wherein 0 < E≤1;
Extract all PPG Wave data sections of E > 0.9;
From the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, when obtaining waveform Between T;
According to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein:
Obtain crest Ps, trough Pd, according to following formula calculating K value:
K = P m - P d P s - P d ;
And according to following formula, calculate heart amount output CO:
When calculating spiritual degree at random according to RR interval data, specifically include:
RR interval data is done fast fourier transform, it is thus achieved that FFT data;
FFT data is converted to power spectrum;
According to described power spectrum, calculate the power P between the 3rd predeterminated frequency scope;
According to described power spectrum, calculate general power A between the 4th predeterminated frequency scope;
By the ratio of P Yu A, calculate spiritual degree at random;
When calculating heart volume index according to RR interval data, specifically include:
Described RR interval data is done fast fourier transform, it is thus achieved that FFT data;
Described FFT data is converted to power spectrum;
According to described power spectrum, calculate power P V between the 5th predeterminated frequency scope;
According to described power spectrum, calculate power P L between the 6th predeterminated frequency scope;
According to described power spectrum, calculate power P H between the 7th predeterminated frequency scope;
Heart rate HR is calculated according to RR interval data;
According to following formula calculating heart volume index HA:
H A = P H P V + P L 2 × H R .
Heart sign acquisition methods the most according to claim 5, it is characterised in that described acquisition PPG data, specifically wraps Include:
Obtain the detection video of predetermined time period;
Following method is used to obtain the PPG numerical value of each two field picture in described detection video:
Obtain the red channel numerical value of each pixel in image, and the described red channel numerical value of all pixels is tired out Add, form the PPG numerical value of described image;
Using all images PPG numerical value as described PPG data;
Wherein, described detection video includes: finger detection video, cervical region detection video, forehead detection video, wrist detection regard Frequently, instep detection video at least one.
7. a heart sign acquisition device, it is characterised in that including:
PPG data acquisition module, is used for obtaining PPG data;
Filtration module, for described PPG numerical value does the cut-off frequency low-pass filtering treatment less than or equal to 4HZ, removes high frequency and makes an uproar Sound;And the PPG numerical value removing high-frequency noise is carried out discrete wavelet Filtering Processing, revise baseline drift;
Crest and wave trough position determine module, according to the described PPG data of correction baseline drift, obtain PPG waveform medium wave peak And the position of trough, and according to described crest and the position of described trough and the sample rate of described PPG data, obtain ripple Time interval numerical value between peak;
Heart sign computing module, for according to the time interval numerical value between described crest, draws RR interval data, and according to Described RR interval data calculates heart sign.
Heart sign acquisition device the most according to claim 7, it is characterised in that also include:
Crest height acquisition module, for the position according to described crest place, obtains the height of described crest, and to described ripple The height at peak is normalized, and obtains normalized crest height numerical value;
Average wave peak heights Numerical Simulation Module, for calculating the meansigma methods of all normalized described crest height numerical value;
Height Anomalies crest ignores module, for being compared with described meansigma methods respectively by normalized described crest height numerical value Right, and according to the result of comparison, ignore and exceed the Height Anomalies crest of described meansigma methods 50% with described meansigma methods difference.
Heart sign acquisition device the most according to claim 8, it is characterised in that also include:
Peak time difference acquisition module, for according to the position at described crest place, calculate two the most adjacent described crests it Between time difference;
Comparing module, for comparing described time difference with the threshold value preset;
Time anomaly crest ignores module, for when described time difference is less than described default threshold value, ignoring two described ripples One of them among peak.
10. according to the heart sign acquisition device described in claim 7-9 any one, it is characterised in that described heart sign Computing module includes: heart rate variability rate calculating sub module, heart be concerned with Index for Calculation submodule, the heart amount output calculating sub module, Spirit degree calculating sub module at random and heart volume index calculating sub module at least one;
Wherein, described heart rate variability rate calculating sub module is used for:
According to RR interval data, using the phase between RR as vertical coordinate, with i as abscissa, draw HRV signal curve;
According to following formula calculating HRV Time-Domain Technique index:
NNVGR is the meansigma methods of phase between the most normal NN;
SDNN is the standard deviation of phase between the most normal NN;
RNSSD is the root-mean-square value that between whole adjacent NN, the phase is the poorest;
SDSD is the standard deviation of the difference of phase between whole adjacent NN;
NN50, for the heart rate more than 50ms of the difference of phase between whole adjacent NN;
PNN50 is the percentage ratio of the difference of the phase heart rate more than 50ms between whole adjacent NN;
In above-mentioned formula, N is normal total heart beats, RRiIt is i-th RR interval data, RRi+1It it is i+1 RR interval data.It it is the meansigma methods of the RR interval data of N number of heartbeat;Between NN, the phase is hole IBI;
RR′i=RRi-RRi-1
RR ′ ‾ = ( RR i ‾ - RR i - 1 ) 2 ;
The described heart Index for Calculation submodule that is concerned with is used for:
Described RR interval data is done fast fourier transform, generates FFT data;
FFT data is converted to power spectrum;
According to described power spectrum, calculate the general power in the range of the first predeterminated frequency;
According to described power spectrum, calculate peak value peak power in the range of the second predeterminated frequency;
By the ratio of described peak value peak power Yu described general power, calculate heart and be concerned with index;
Described heart amount output calculating sub module is used for:
Described PPG waveform is done data segment cutting, obtains PPG Wave data section;
To PPG Wave data section each described: described PPG Wave data section to be done wavelet transformation, revise baseline drift;Calculate Revise the meansigma methods of the PPG Wave data section medium wave peak height after baseline drift;Calculate RR crest height in PPG Wave data section And the ratio between the meansigma methods of crest height, and using ratio minima as waveform quality E, wherein 0 < E≤1;
Extract all PPG Wave data sections of E > 0.9;
From the PPG Wave data section extracted, identify the position of crest and trough, and intercept complete crest, when obtaining waveform Between T;
According to described waveform time T, crest numerical value P (t) is done integral operation, calculate average value Pm, wherein:
Obtain crest Ps, trough Pd, according to following formula calculating K value:
K = P m - P d P s - P d ;
And according to following formula, calculate heart amount output CO:
Described spirit degree at random calculating sub module is used for:
RR interval data is done fast fourier transform, it is thus achieved that FFT data;
FFT data is converted to power spectrum;
According to described power spectrum, calculate the power P between the 3rd predeterminated frequency scope;
According to described power spectrum, calculate general power A between the 4th predeterminated frequency scope;
By the ratio of P Yu A, calculate spiritual degree at random;
Described heart volume index calculating sub module is used for:
Described RR interval data is done fast fourier transform, it is thus achieved that FFT data;
Described FFT data is converted to power spectrum;
According to described power spectrum, calculate power P V between the 5th predeterminated frequency scope;
According to described power spectrum, calculate power P L between the 6th predeterminated frequency scope;
According to described power spectrum, calculate power P H between the 7th predeterminated frequency scope;
Heart rate HR is calculated according to RR interval data;
According to following formula calculating heart volume index HA:
H A = P H P V + P L 2 × H R .
CN201610571030.7A 2016-07-19 2016-07-19 A kind of heart sign acquisition methods and device CN106175742A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610571030.7A CN106175742A (en) 2016-07-19 2016-07-19 A kind of heart sign acquisition methods and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610571030.7A CN106175742A (en) 2016-07-19 2016-07-19 A kind of heart sign acquisition methods and device

Publications (1)

Publication Number Publication Date
CN106175742A true CN106175742A (en) 2016-12-07

Family

ID=57494415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610571030.7A CN106175742A (en) 2016-07-19 2016-07-19 A kind of heart sign acquisition methods and device

Country Status (1)

Country Link
CN (1) CN106175742A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106377226A (en) * 2016-08-30 2017-02-08 苏州品诺维新医疗科技有限公司 Cardiovascular state detection device and determining method for slope information
CN107595271A (en) * 2017-07-27 2018-01-19 康美健康云服务有限公司 Method for measuring heart rate, electronic equipment, storage medium and system based on mobile phone
CN108451494A (en) * 2017-02-17 2018-08-28 祥明大学校产学协力团 The method and system of time domain cardiac parameters are detected using pupillary reaction

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012134395A1 (en) * 2011-03-25 2012-10-04 Nitto Denko Corporation A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system
WO2014089665A1 (en) * 2012-12-13 2014-06-19 Cnv Systems Ltd. System for measurement of cardiovascular health
CN104127194A (en) * 2014-07-14 2014-11-05 华南理工大学 Depression evaluating system and method based on heart rate variability analytical method
CN105433936A (en) * 2016-01-29 2016-03-30 北京心量科技有限公司 Method and device for obtaining cardiac output
CN105550530A (en) * 2016-01-29 2016-05-04 北京心量科技有限公司 Method and device for obtaining vikshipta degree
CN105595991A (en) * 2016-01-29 2016-05-25 北京心量科技有限公司 Method and device for measuring cardiac coherence index

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012134395A1 (en) * 2011-03-25 2012-10-04 Nitto Denko Corporation A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system
WO2014089665A1 (en) * 2012-12-13 2014-06-19 Cnv Systems Ltd. System for measurement of cardiovascular health
CN104127194A (en) * 2014-07-14 2014-11-05 华南理工大学 Depression evaluating system and method based on heart rate variability analytical method
CN105433936A (en) * 2016-01-29 2016-03-30 北京心量科技有限公司 Method and device for obtaining cardiac output
CN105550530A (en) * 2016-01-29 2016-05-04 北京心量科技有限公司 Method and device for obtaining vikshipta degree
CN105595991A (en) * 2016-01-29 2016-05-25 北京心量科技有限公司 Method and device for measuring cardiac coherence index

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106377226A (en) * 2016-08-30 2017-02-08 苏州品诺维新医疗科技有限公司 Cardiovascular state detection device and determining method for slope information
CN106377226B (en) * 2016-08-30 2019-05-03 苏州涵轩信息科技有限公司 A kind of determination method of cardiovascular status detection device and slope information
CN108451494A (en) * 2017-02-17 2018-08-28 祥明大学校产学协力团 The method and system of time domain cardiac parameters are detected using pupillary reaction
CN107595271A (en) * 2017-07-27 2018-01-19 康美健康云服务有限公司 Method for measuring heart rate, electronic equipment, storage medium and system based on mobile phone

Similar Documents

Publication Publication Date Title
Poh et al. Advancements in noncontact, multiparameter physiological measurements using a webcam
Tarvainen et al. An advanced detrending method with application to HRV analysis
JP2010510850A (en) Biomedical signal analysis method including improved automatic segmentation
US7239988B2 (en) Apparatus and method for efficient representation of periodic and nearly periodic signals for analysis
US20140303454A1 (en) Remote monitoring of vital signs
US20070179396A1 (en) Method and System for Detecting and Classifying Facial Muscle Movements
US7720306B2 (en) Systems and methods for displaying changes in biological responses to therapy
JP2004531309A (en) Side noise analysis support device
US20050004479A1 (en) Locating features in a photoplethysmograph signal
Wyart et al. Rhythmic fluctuations in evidence accumulation during decision making in the human brain
Bousefsaf et al. Continuous wavelet filtering on webcam photoplethysmographic signals to remotely assess the instantaneous heart rate
US8649562B2 (en) Method and system for processing a signal including at least a component representative of a periodic phenomenon in a living being
Sameni et al. A nonlinear Bayesian filtering framework for ECG denoising
US20110251493A1 (en) Method and system for measurement of physiological parameters
US9913588B2 (en) Method and system for screening of atrial fibrillation
US8152731B2 (en) Wavelet transform and pattern recognition method for heart sound analysis
JP2015521075A (en) Extraction of respiratory rate from cardiac signal
JP2007535392A (en) Apparatus and method for analysis of high frequency QRS complex
WO2003005879A3 (en) Method and apparatus for monitoring cardiac patients for t-wave alternans
DE19703403A1 (en) Robust time diversity filter for removing electromagnetic interference
Hsu Continuous EEG signal analysis for asynchronous BCI application
Zidelmal et al. QRS detection using S-Transform and Shannon energy
US9025826B2 (en) Formation of a time-varying signal representative of at least variations in a value based on pixel values
US9301712B2 (en) Method and apparatus for continuous measurement of motor symptoms in parkinson's disease and essential tremor with wearable sensors
US20030210047A1 (en) Filtering artifact from fMRI data using the Stockwell transform

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination