CN109589130A - Heart vitality index calculation method - Google Patents
Heart vitality index calculation method Download PDFInfo
- Publication number
- CN109589130A CN109589130A CN201811433965.4A CN201811433965A CN109589130A CN 109589130 A CN109589130 A CN 109589130A CN 201811433965 A CN201811433965 A CN 201811433965A CN 109589130 A CN109589130 A CN 109589130A
- Authority
- CN
- China
- Prior art keywords
- heart
- heart sound
- movement
- score
- tranquillization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B7/00—Instruments for auscultation
- A61B7/006—Detecting skeletal, cartilage or muscle noise
Abstract
The present invention proposes a kind of heart vitality index calculation method, belongs to heart vitality index calculating field.Technical solution of the present invention main points are as follows: acquired under quiescent condition respectively with post exercise cardiechema signals, and carry out pre-process after resampling and auto-correlation segmentation；It determines first heart sound amplitude maximum S1 and second heart sound amplitude maximum S2, finds out SMovement/STranquillizationRelative value, and according to paragraph number by the relative value take it is average be used as first characteristic parameter；The S1/S2 under quiescent condition is calculated, and the S1/S2 is taken by second characteristic parameter of average conduct according to paragraph number；The ratio for calculating each paragraph diastole time limit and time limit in systole phase under quiescent condition takes average as third characteristic value；Characteristic parameter each in test sample is scaled to the score of heart vitality index corresponding with criteria scores respectively, is denoted as ScoreS respectivelyMovement/STranquillizationRelative value, ScoreS1/S2 and ScoreD/S；By the score weight multiplied by each score and addition respectively, heart vigor score Score to the end is obtained.
Description
Technical field
The present invention relates to heart vitality index calculating fields, in particular to based on the index calculating side of heart sound resampling signal
The field of method.
Background technique
Cardiac parameter facilitates diagnosis of cardiovascular diseases, observes the state of an illness, infers the pathogenetic pathology of disease, selects treatment side
Method, estimating prognosis etc..S1/S2 can effectively reflect the ability of movement up-regulation myocardial contractive power, SMovement/STranquillizationRelative value can reflect
Employ the ability of cardiac reserve, wherein S1 is first heart sound amplitude maximum, and S2 is second heart sound amplitude maximum, SMovement
For S1/S2, S after movementTranquillizationFor S1/S2 under quiescent condition, this 2 indexs can under effectively evaluating motion state heart function.
D/S is related to ejection fraction, shows that it evaluates that cardiac function is relatively reliable, it can heart function under effectively evaluating rest state
Energy.Cardiechema signals evaluation heart activity level is based on this and generates.
Summary of the invention
The object of the present invention is to provide a kind of heart vitality index calculation methods, can evaluate heart according to cardiechema signals
Vigor.
The present invention solves its technical problem, the technical solution adopted is that: heart vitality index calculation method, including following step
It is rapid:
A, respectively acquire quiescent condition under cardiechema signals and post exercise cardiechema signals, by the cardiechema signals after acquisition into
It is pre-processed after row resampling；
B, pretreated cardiechema signals are subjected to auto-correlation segmentation, then enter step C, D and E simultaneously；
C, the signal after being segmented according to heart sound determines first heart sound amplitude maximum and second heart sound amplitude maximum, the
One heart sound amplitude maximum is denoted as S1, and second heart sound amplitude maximum is denoted as S2, finds out SMovement/STranquillizationRelative value, SMovementFor movement
S1/S2 afterwards, STranquillizationIt is S1/S2 under quiescent condition, the ratio note of first heart sound amplitude maximum and second heart sound amplitude maximum
For S1/S2, and the relative value is averaged according to paragraph number, the result obtained is as SMovement/STranquillizationRelative value end value, should
End value is as first characteristic parameter；
D, the S1/S2 under quiescent condition is calculated, and is averaged the S1/S2 according to paragraph number, calculated each paragraph
S1/S2 average value as second characteristic parameter；
E, the ratio in each paragraph diastole time limit and time limit in systole phase under quiescent condition is calculated, and is averaged, is made
For third characteristic value, the ratio in the diastole time limit and time limit in systole phase is denoted as D/S；
F, sportsman's criteria scores are set, by fisrt feature parameter, second feature parameter and third feature in test sample
Parameter is scaled the score of heart vitality index corresponding with criteria scores respectively, is denoted as ScoreS respectivelyMovement/STranquillizationIt is opposite
Value, ScoreS1/S2 and ScoreD/S；
G, by ScoreSMovement/STranquillizationRelative value, ScoreS1/S2 and ScoreD/S are respectively multiplied by the weight of each score and phase
Add, obtains heart vigor score Score to the end namely heart vitality index.
Further, the run duration is the time of running in 10 minutes in step A.
Particularly, in step A, 5 resamplings, sample frequency 2205Hz are carried out to collected cardiechema signals.
It is further, in step A, after the cardiechema signals by after acquisition carry out resampling, the pretreatment of progress
Include:
A1, Butterworth low pass wave is carried out to the cardiechema signals after resampling；
A2, filtered cardiechema signals are denoised.
Particularly, in step A2, when being filtered to the signal after resampling, setting passband maximum attenuation is 3db, resistance
Band minimal attenuation is 18db, in step A3, is denoised using dmey wavelet transformation to filtered cardiechema signals.
It is further in step B, pretreated cardiechema signals to be subjected to auto-correlation segmentation, are specifically included as follows
Step:
B1, the amplitude equalizing value for calculating the heart sound data after denoising；
B2, setup parameter: after resampling, the minimal point s_min=750 in a heart sound period, the maximum in heart sound period
Count s_max=2500；
B3, it is divided by obtain the segments of the cardiechema signals by the maximum length of cardiechema signals length and a heart sound period
M；
B4, the successively combination of two from start to end by this M sections of heart sound data, have extra segmentation if M is odd number, more
Remaining segmentation is cast out；
B5, the starting point for finding first first heart sound in each combined segment；
Heart sound period where the first heart sound starting point found in B6, calculating combined segment and the next heart sound in this section
The auto-correlation coefficient in period, and save；
The heart sound period in B7, all combined segments of selection where maximum auto-correlation coefficient, as segmentation result.
Particularly, in step B5, using continuous 100 points be less than heart sound data amplitude mean value as judgement first heart sound
The condition of starting point.
It is further that in step F, the criteria scores of the sportsman set is 100 points.
Particularly, in step G, the calculation formula of the heart vigor score Score are as follows:
Score=ScoreSMovement/STranquillizationRelative value * coefSMovement/STranquillizationRelative value+ScoreD/S*coefD/S
+ScoreS1/S2*coefS1/S2
Wherein, coefSMovement/STranquillizationRelative value be ScoreSMovement/STranquillizationRelative value weight, coefD/S ScoreD/
The weight of S, coefS1/S2 are the weight of ScoreS1/S2.
It is further that in step G, the weight of each score can carry out free setting, the ScoreSMovement/STranquillization's
The weight that the weight that the weight of relative value is 0.5, ScoreD/S is 0.2, ScoreS1/S2 is 0.3.
The invention has the advantages that by above-mentioned heart vitality index calculation method, firstly, acquiring quiescent condition respectively
Under cardiechema signals and post exercise cardiechema signals, will cardiechema signals after acquisition carry out resampling after pre-process；Secondly,
Pretreated cardiechema signals are subjected to auto-correlation segmentation；Then, the signal after being segmented according to heart sound determines first heart sound amplitude
Maximum value and second heart sound amplitude maximum, first heart sound amplitude maximum are denoted as S1, and second heart sound amplitude maximum is denoted as S2,
Find out SMovement/STranquillizationRelative value, SMovementFor S1/S2, S after movementTranquillizationIt is S1/S2 under quiescent condition, first heart sound amplitude maximum
The ratio of value and second heart sound amplitude maximum is denoted as S1/S2, and is averaged the relative value according to paragraph number, the result obtained
As SMovement/STranquillizationRelative value end value, the end value is as first characteristic parameter；Meanwhile it calculating under quiescent condition
S1/S2, and be averaged the S1/S2 according to paragraph number, the S1/S2 average value of calculated each paragraph is as second feature
Parameter；Meanwhile the ratio in each paragraph diastole time limit and time limit in systole phase under quiescent condition is calculated, and be averaged, make
For third characteristic value, the ratio in the diastole time limit and time limit in systole phase is denoted as D/S；Subsequently, sportsman is set
Criteria scores, by fisrt feature parameter, second feature parameter in test sample and third feature parameter is scaled respectively and standard
The score of the corresponding heart vitality index of score, is denoted as ScoreS respectivelyMovement/STranquillizationRelative value, ScoreS1/S2 and
ScoreD/S；Finally, by ScoreSMovement/STranquillizationRelative value, ScoreS1/S2 and ScoreD/S are respectively multiplied by the weight of each score
And be added, obtain heart vigor score Score to the end namely heart vitality index.
The present invention evaluates the method for heart vigor using cardiechema signals with noninvasive, sensibility and specificity is high, easy
It is easy, while convenient for being transplanted to all kinds of software platforms, it is widely used.
Specific embodiment
Technical solution of the present invention is described below in detail.
Heart vitality index calculation method of the present invention when concrete application, is divided into following steps:
1) the tranquillization cardiechema signals of acquisition and post exercise cardiechema signals are subjected to resampling, the phases such as filtering, denoising respectively
Close pretreatment:
A1. resampling is carried out to the heart sound data received；
A2. Butterworth low pass wave is carried out to the signal after resampling；
A3. filtered cardiechema signals are denoised.
Wherein, step A1 carries out 5 resamplings, sample frequency 2205Hz to the cardiechema signals received；Step A2 counterweight
Signal after sampling is filtered, and setting passband maximum attenuation is 3db, minimum attenuation in stop band 18db；Step A3 uses small echo
Transformation denoises filtered cardiechema signals, using dmey small echo.
2) pretreated signal is subjected to auto-correlation segmentation:
B1. the amplitude equalizing value of the heart sound data after denoising is calculated；
B2. setup parameter: after resampling, the minimal point s_min=750 in a heart sound period, the maximum in heart sound period
Count s_max=2500；
B3. it is divided by obtain the segments of the cardiechema signals by the maximum length of cardiechema signals length and a heart sound period
M；
B4. the successively combination of two from start to end by this M sections of heart sound data, has extra segmentation if M is odd number, more
Remaining segmentation is cast out；
B5. the starting point of first first heart sound in each combined segment is found；
B6. the heart sound period where the first heart sound starting point found in calculating combined segment and the next heart sound in this section
The auto-correlation coefficient in period, and save；
B7. the heart sound period in all combined segments where maximum auto-correlation coefficient is selected, as segmentation result.
Wherein, step B5 avoids starting point from being placed exactly at first heart sound, and the amplitude of heart sound data is less than with continuous 100 points
Mean value is as the condition for determining first heart sound starting point.
3) multiple paragraph signals after being segmented according to heart sound determine that first heart sound and second heart sound are most in each paragraph
Amplitude finds out S1/S2, wherein S1, S2 are respectively first heart sound amplitude maximum and second heart sound amplitude maximum, SMovement/
STranquillizationRelative value and D/S, and taken according to paragraph number these values and average be used as final argument value.
4) cardiac reserve of cardiechema signals is assessed:
The average value of the performance of athletic heart is set to the cardiechema signals of " most standard ".It will be in all sportsman's samples
S1/S2、SMovement/STranquillizationRelative value and D/S average value as full marks desired value, parameters are scored at equal with sportsman
It is worth relevant value, has
ScoreSMovement/STranquillizationRelative value=f (StandardSMovement/STranquillizationRelative value)
ScoreD/S=g (Standard S1/S2)
ScoreS1/S2=h (StandardS1/S2)
Wherein, each parameter differs bigger with sportsman's desired value, and score is lower.
5) by ScoreS obtained in 4)Movement/STranquillizationRelative value, ScoreD/S, ScoreS1/S2 are respective multiplied by them
Weight coefficient is simultaneously added, and obtains final score, i.e. heart vitality index:
Score=ScoreSMovement/STranquillizationRelative value * coefSMovement/STranquillizationRelative value+ScoreD/S*coefD/S
+ScoreS1/S2*coefS1/S2
In above-mentioned step 4), one of parameter (ScoreS is takenMovement/STranquillizationRelative value) illustrate, such as by big
This feature value that amount sportsman's data sample obtains is 2, and corresponding is 100 points, and the value that test sample obtains is 1.8, that
The score of test sample is exactly 90, linear to calculate.Step 4) is to calculate extracted three features and movement in test sample
The value that member's full marks Comparison of standards obtains, the value for exactly setting these parameters of sportsman is all 100 points, we are according to current sample
Originally the value measured estimates score, for example practical score is respectively 80,70,60.So step 5) be exactly according to these values,
Multiplied by a weight coefficient, an end value is obtained, wherein weight coefficient is chosen according to the experience of many experiments, without spy
Other algorithm, the weight of each score can carry out free setting, in this example, ScoreSMovement/STranquillizationRelative value weight be 0.5,
The weight that the weight of ScoreD/S is 0.2, ScoreS1/S2 is 0.3.So according to above example final score score=
80*0.5+70*0.2+60*0.3=72 points.
The index score result has following characteristics: full marks are 100 points, and the heart of full marks scorekeeper and sportsman are very
It is similar, it is a standard of healthy heart sound；Reflect heart vigor condition by single score, is easily understood.Training can be with
The heart vigor for improving people, so the invention can supervise people to take exercise.
Claims (10)
1. heart vitality index calculation method, which comprises the following steps:
A, the cardiechema signals and post exercise cardiechema signals under quiescent condition are acquired respectively, and the cardiechema signals after acquisition are subjected to weight
It is pre-processed after sampling；
B, pretreated cardiechema signals are subjected to auto-correlation segmentation, then enter step C, D and E simultaneously；
C, the signal after being segmented according to heart sound determines first heart sound amplitude maximum and second heart sound amplitude maximum, first heart
Sound amplitude maximum is denoted as S1, and second heart sound amplitude maximum is denoted as S2, finds out SMovement/STranquillizationRelative value, SMovementAfter movement
S1/S2, STranquillizationIt is S1/S2 under quiescent condition, the ratio of first heart sound amplitude maximum and second heart sound amplitude maximum is denoted as
S1/S2, and be averaged the relative value according to paragraph number, the result obtained is as SMovement/STranquillizationRelative value end value, this is most
Final value is as first characteristic parameter；
D, the S1/S2 under quiescent condition is calculated, and is averaged the S1/S2 according to paragraph number, the S1/ of calculated each paragraph
S2 average value is as second characteristic parameter；
E, the ratio for calculating each paragraph diastole time limit and time limit in systole phase under quiescent condition, and is averaged, as the
The ratio in three characteristic values, the diastole time limit and time limit in systole phase is denoted as D/S；
F, sportsman's criteria scores are set, by fisrt feature parameter, second feature parameter and third feature parameter in test sample
It is scaled the score of heart vitality index corresponding with criteria scores respectively, is denoted as ScoreS respectivelyMovement/STranquillizationRelative value,
ScoreS1/S2 and ScoreD/S；
G, by ScoreSMovement/STranquillizationRelative value, ScoreS1/S2 and the ScoreD/S weight multiplied by each score and addition respectively,
Obtain heart vigor score Score namely heart vitality index to the end.
2. heart vitality index calculation method according to claim 1, which is characterized in that in step A, the run duration
For the time of running in 10 minutes.
3. heart vitality index calculation method according to claim 1, which is characterized in that in step A, to the collected heart
Sound signal carries out 5 resamplings, sample frequency 2205Hz.
4. heart vitality index calculation method according to claim 1, which is characterized in that in step A, it is described will acquisition after
Cardiechema signals carry out resampling after, the pretreatment of progress includes:
A1, Butterworth low pass wave is carried out to the cardiechema signals after resampling；
A2, filtered cardiechema signals are denoised.
5. heart vitality index calculation method according to claim 1, which is characterized in that in step A2, after resampling
Signal when being filtered, setting passband maximum attenuation is 3db, small using dmey in minimum attenuation in stop band 18db, step A3
Wave conversion denoises filtered cardiechema signals.
6. heart vitality index calculation method according to claim 1 or 4 or 5, which is characterized in that in step B, will locate in advance
Cardiechema signals after reason carry out auto-correlation segmentation, specifically comprise the following steps:
B1, the amplitude equalizing value for calculating the heart sound data after denoising；
B2, setup parameter: after resampling, the minimal point s_min=750 in a heart sound period, the maximum number of points in heart sound period
S_max=2500；
B3, it is divided by obtain the segments M of the cardiechema signals by the maximum length of cardiechema signals length and a heart sound period；
B4, the successively combination of two from start to end by this M sections of heart sound data, have extra segmentation if M is odd number, extra
Segmentation is cast out；
B5, the starting point for finding first first heart sound in each combined segment；
Heart sound period where the first heart sound starting point found in B6, calculating combined segment and the next heart sound period in this section
Auto-correlation coefficient, and save；
The heart sound period in B7, all combined segments of selection where maximum auto-correlation coefficient, as segmentation result.
7. heart vitality index calculation method according to claim 6, which is characterized in that in step B5, with continuous 100
Point is less than the amplitude mean value of heart sound data as the condition for determining first heart sound starting point.
8. heart vitality index calculation method according to claim 1, which is characterized in that in step F, the fortune of the setting
The criteria scores of mobilization are 100 points.
9. heart vitality index calculation method according to claim 1, which is characterized in that in step G, the heart vigor
The calculation formula of score Score are as follows:
Score=ScoreSMovement/STranquillizationRelative value * coefSMovement/STranquillizationRelative value+ScoreD/S*coefD/S
+ScoreS1/S2*coefS1/S2
Wherein, coefSMovement/STranquillizationRelative value be ScoreSMovement/STranquillizationRelative value weight, coefD/S is ScoreD/S
Weight, coefS1/S2 are the weight of ScoreS1/S2.
10. according to claim 1 or heart vitality index calculation method described in 9, which is characterized in that in step G, each score
Weight can carry out free setting, the ScoreSMovement/STranquillizationRelative value weight be 0.5, ScoreD/S weight be 0.2,
The weight of ScoreS1/S2 is 0.3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811433965.4A CN109589130A (en) | 2018-11-28 | 2018-11-28 | Heart vitality index calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811433965.4A CN109589130A (en) | 2018-11-28 | 2018-11-28 | Heart vitality index calculation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109589130A true CN109589130A (en) | 2019-04-09 |
Family
ID=65960512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811433965.4A Pending CN109589130A (en) | 2018-11-28 | 2018-11-28 | Heart vitality index calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109589130A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006123774A1 (en) * | 2005-05-19 | 2006-11-23 | Matsushita Electric Works, Ltd. | Sensor assembly for sleep apnea syndrome examination and sleep apnea syndrome examination instrument |
| US20070265533A1 (en) * | 2006-05-12 | 2007-11-15 | Bao Tran | Cuffless blood pressure monitoring appliance |
| US20110245651A1 (en) * | 2010-03-31 | 2011-10-06 | Fujifilm Corporation | Medical image playback device and method, as well as program |
| CN107822597A (en) * | 2017-10-17 | 2018-03-23 | 四川长虹电器股份有限公司 | Based on the calculating of intelligent platform heart vitality index and image conversion display system and method |
| CN108682433A (en) * | 2018-06-01 | 2018-10-19 | 四川长虹电器股份有限公司 | The heart sound kind identification method of first-order difference coefficient based on MFCC |
-
2018
- 2018-11-28 CN CN201811433965.4A patent/CN109589130A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006123774A1 (en) * | 2005-05-19 | 2006-11-23 | Matsushita Electric Works, Ltd. | Sensor assembly for sleep apnea syndrome examination and sleep apnea syndrome examination instrument |
| US20070265533A1 (en) * | 2006-05-12 | 2007-11-15 | Bao Tran | Cuffless blood pressure monitoring appliance |
| US20110245651A1 (en) * | 2010-03-31 | 2011-10-06 | Fujifilm Corporation | Medical image playback device and method, as well as program |
| CN107822597A (en) * | 2017-10-17 | 2018-03-23 | 四川长虹电器股份有限公司 | Based on the calculating of intelligent platform heart vitality index and image conversion display system and method |
| CN108682433A (en) * | 2018-06-01 | 2018-10-19 | 四川长虹电器股份有限公司 | The heart sound kind identification method of first-order difference coefficient based on MFCC |
Non-Patent Citations (2)
| Title |
|---|
| 张敏: "运动心力监测评价心脏储备的临床研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
| 胡彦: "负重越野训练对战士心功能的影响", 《中国应用生理学杂志》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106108889B (en) | Electrocardiogram classification method based on deep learning algorithm | |
| D'Andrea et al. | Right heart structural and functional remodeling in athletes | |
| CN108577883A (en) | A kind of Screening for coronary artery disease device, screening system and signal characteristic extracting methods | |
| CN109222949A (en) | Heart rate detection method and heartbeat detection device | |
| CA2565192A1 (en) | Apparatus and method for analysis of high frequency qrs complexes | |
| CN102509267B (en) | Retrospective off-line gating method for intravascular ultrasound image sequence | |
| EP1304074A3 (en) | Method and apparatus for determining the left-ventricular ejection time tlev of a heart of a subject | |
| CN110037668B (en) | System for judging age, health state and malignant arrhythmia identification by combining pulse signal time-space domain with model | |
| US9451901B2 (en) | Method and device for evaluation of myocardial damages based on the current density variations | |
| CN109381181A (en) | The end-point detecting method of electrocardiosignal signature waveform | |
| D Andrea et al. | Effects of different training protocols on left ventricular myocardial function in competitive athletes: a Doppler tissue imaging study | |
| CN108670291A (en) | The heart sound kind identification method of improved MFCC is combined based on EMD | |
| George et al. | Echocardiographic examination of cardiac structure and function in elite cross trained male and female Alpine skiers. | |
| CN109864713A (en) | Rhythm of the heart method based on multi-channel parallel filtering and spectral peak weighting selection algorithm | |
| US20110270059A1 (en) | Signal processing for pulse oximetry | |
| CN109330582A (en) | Heart rate and its characteristic index detection method based on ECG Signal Analysis | |
| Lin et al. | Investigation on pulse wave forward peak detection and its applications in cardiovascular health | |
| CN109589130A (en) | Heart vitality index calculation method | |
| CN102334986A (en) | Method for identifying dicrotic notch point in arterial tension signal | |
| KR101273652B1 (en) | A method to detect myocardial ischemia using comparison analysis of heart rate variability time and frequency domain features | |
| Li et al. | An approach of heartbeat segmentation in seismocardiogram by matched-filtering | |
| CN109431535A (en) | A kind of cardiac reserve index calculation method | |
| Wu et al. | Main pulmonary arterial distensibility different presentation between chronic pulmonary hypertension and acute pulmonary embolism | |
| Goncalvesova et al. | Changes of portal flow in heart failure patients with liver congestion | |
| CN104706373B (en) | Heart vital index calculating method based on heart sounds |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190409 |
|
| RJ01 | Rejection of invention patent application after publication |