CN106037792B - A kind of heart sound wave data processing and graphic display method and its auscultation system - Google Patents
A kind of heart sound wave data processing and graphic display method and its auscultation system Download PDFInfo
- Publication number
- CN106037792B CN106037792B CN201610570124.2A CN201610570124A CN106037792B CN 106037792 B CN106037792 B CN 106037792B CN 201610570124 A CN201610570124 A CN 201610570124A CN 106037792 B CN106037792 B CN 106037792B
- Authority
- CN
- China
- Prior art keywords
- heart sound
- peak
- data
- heart
- characteristic point
- 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.)
- Active
Links
- 238000012545 processing Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000002555 auscultation Methods 0.000 title claims abstract description 34
- 238000003672 processing method Methods 0.000 claims abstract description 10
- 238000012216 screening Methods 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000003745 diagnosis Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 2
- 230000037430 deletion Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 13
- 210000003128 head Anatomy 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 230000013707 sensory perception of sound Effects 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 208000020446 Cardiac disease Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 210000000613 ear canal Anatomy 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B7/00—Instruments for auscultation
- A61B7/02—Stethoscopes
- A61B7/04—Electric stethoscopes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
-
- 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
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/15—Biometric patterns based on physiological signals, e.g. heartbeat, blood flow
Abstract
The invention discloses a kind of heart sound wave data processings and graphic display method and its auscultation system, any parameter is inputted without user by heart sound wave data processing method, can accurately and automatically identify the feature and extract real-time heart rate of first heart sound and second heart sound;Pass through heart sound graphic display method, simultaneously by treated the similar electrocardiogram displaying of original caardiophonogram and data on same figure, share a time shaft, convenient for the comparison up and down of characteristic point, both the original fluctuation situation of the original caardiophonogram of observable, can also from the similar electrocardiogram after data processing clearly heart sound feature, realize intuitiveization of heart sound figure, clear caardiophonogram is obtained, makes to diagnose more acurrate;By electronic auscultation system, simultaneously extract real-time heart rate can be recorded in real time, shows clear caardiophonogram.
Description
Technical field
The present invention relates to a kind of processing of electronic auscultation device heart sound data and display technology, in particular to a kind of heart sound waveform numbers
According to processing and graphic display method and its auscultation system.
Background technique
Medical staff is often living by obtaining human viscera organ outside stethoscope by means of stethoscope in clinical position
Dynamic sound, medical staff can tentatively judge the cause of disease.
Heart sound is to reflect that heart is normal or the sound equipment of pathology, it is the opening and closing due to valve each during heartbeat
And vibration formation caused by cardiac muscle and blood motion.Cardiophony is the important evidence of Diagnosing Cardiac disease.Traditional
Method is using auscultation with stethoscope heart sound, and the foundation of diagnosis is mainly the experience of doctor.
However, old auscultation with stethoscope heart sound is often difficult to capture the sending of inside of human body internal organs although method is simple
Some faint but very important biological sound, cause doctor that can not make diagnosis in time, and diagnosis according to main basis
The experience of doctor, accuracy are poor.It is said from another angle, human ear is the comprehensive effect of the sound intensity and frequency to the perception of sound, because
And some pathological characteristics are difficult to capture, and accuracy is poor.
Meanwhile there is the discomfort of the external auditory canal as caused by pressure pipe stress problems, sound quality vulnerable to disadvantages such as interference in conventional stethoscope
End, is used for a long time conventional stethoscope, and the hearing of medical staff is affected.
This just needs a kind of novel electronic auscultation device to carry out quantitative, accurately analysis to auscultatory sound.Domestic and state at present
Exoelectron stethoscope product differing greatly in price, function, external electronic auscultation device properties of product are preferably at high price,
For home products is more external, there is also larger gaps for performance, and majority amplifies after only acquiring heart sound, although to a certain degree
On improve stethoscope, but heart sound waveform can't be well reflected.
Professional person both knows about in the industry, and the sound data processing to be carried out that electronic auscultation device obtains is very complicated, is not easy to obtain
To the physiology and pathological data model of heart, how extract real-time heart rate, and obtain interpretation clear, convenient for caardiophonogram
Visual caardiophonogram is all difficult to the technical issues of solving.
Meanwhile the electronic auscultation device instrument of a set of simple, cheap, novel audio-visual is built, and optimize heart sound
Data processing technique and algorithm, in real time record, extract real-time heart rate, obtain clear caardiophonogram, make to diagnose it is more acurrate,
Convenient for the use of the personnel such as medical care, meaning is also very great.
It is accordingly required in particular to a kind of heart sound wave data processing and graphic display method and its auscultation system, on solving
State existing problem.
Summary of the invention
The first purpose of this invention is to provide a kind of heart sound wave data processing and heart sound graphic display method, for
The heart sound data that electronic auscultation device obtains is converted into the similar figure of electrocardiogram, to identify first by the deficiencies in the prior art
The feature of heart sound and second heart sound and extract real-time heart rate, while by treated the similar electrocardiogram exhibition of original caardiophonogram and data
Show on same figure, convenient for the interpretation of caardiophonogram.
Second object of the present invention is to provide a kind of electronic auscultation system, in view of the deficiencies of the prior art, can be real-time
Simultaneously extract real-time heart rate is recorded, shows clear caardiophonogram, it is simple novel, convenient for the use of the personnel such as medical care.
Technical problem solved by the invention can be realized using following technical scheme:
In a first aspect, the present invention provides a kind of heart sound wave data processing method, which is characterized in that it includes following step
It is rapid:
(1) one section of heart sound data is chosen;
(2) data processing and conversion are carried out to heart sound data, picks up the peak dot of heart sound waveform;
(3) first heart sound and second heart sound are identified, and counts the average heart rate in the segment data.
In one embodiment of the invention, specific step is as follows for the processing of the heart sound data:
Step 1: choosing one section of heart sound data, and data are long to be not less than 3 seconds;
Step 2: data bandpass filtering: frequency filtering range is 20-100Hz;
Step 3: data positive value conversion: data take absolute value;
Step 4: data smoothing filtering;
Step 5: pick up the peak dot of heart sound waveform: empirical parameter can be used in peak-seeking parameter, and guarantee peak dot is effective heart sound
Characteristic point.
Step 6: calculate the characteristic point parameter in each heart sound period: position, peak height and it is peak-to-peak away from;
Step 7: identification first heart sound S1 and second heart sound S2: according to the time attribute of heart sound, to identify first heart sound
S1 and second heart sound S2;
Step 8: it counts the average heart rate in this section of heart sound data: the heart is calculated using the time between adjacent first heart sound S1
Rate and average heart rate.
Further, the length of the heart sound data is 3-10 seconds, and data volume is made to be unlikely to excessive and influence calculating speed.
Further, the characteristic point refers to first heart sound S1, second heart sound S2, third heart sound within a heart sound period
S3, fourth heart sound S4;The characteristic point parameter include peak position, peak height, it is peak-to-peak away from, caardiophonogram adjacent first heart sound S1 when
Between time interval be dS1S1, the time interval of first heart sound S1 and second heart sound S2 indicate in a heart sound period in caardiophonogram
For dS1S2.
Further, in above-mentioned 7th step, first heart sound S1 and second heart sound S2 is identified, the purpose is to handle heart sound data
The characteristic point parameter in each heart sound period that preceding 6 step obtains, is screened, rejects invalid peak dot, and finally automatically identifies the
One heart sound S1 and second heart sound S2, the recognition methods include the following steps:
Step 1: being screened according to peak number, screening technique 1) first calculate maximum effectively peak number, maximum effectively peak number
Calculation formula are as follows: maximum effectively peak number=time span (unit second)/60 ' 120 ' 4;2) if feature point number is no more than most
Big effectively peak number, then do not have to screening;3) if necessary to screen, then first to characteristic point parametric results with peak value from high to low into
Then row sequence is more than the feature point deletion of effective number tail portion;
Step 2: being screened according to peak value, screening technique 1) lift the practical maximum peak value of characteristic point;2) determination has
Imitate peak height minimum value, maximum peak value/4 of effective peak height minimum value=reality;3) peak height is less than to the feature of effective peak height minimum value
Point is rejected;
Step 3: peak-to-peak away from screening according to peak, screening technique is peak-to-peak away from the principle that must not drop below 0.2 second according to peak
Screened, when two it is adjacent it is peak-to-peak away from less than 0.2 second when, the small characteristic point of peak height is rejected, it is peak-to-peak then to search peak again again
Away from, until meet peak it is peak-to-peak away from requirement;
Step 4: identification first heart sound: first heart sound S1 and second heart sound S2 spacing have following time response:
DS1S2 < dS2S1, i.e., in a heart sound period, the time between first heart sound S1 to second heart sound S2 is less than second heart sound S2 and arrives
Time between first heart sound S1;According to the time response of first heart sound S1 and second heart sound S2 spacing, it is peak-to-peak to characteristic point peak away from
The selection result table carries out first heart sound identification.
Step 5: identifying second heart sound: recognition methods is the subsequent peak first heart sound S1 if not first heart sound S1,
It can be labeled as S2, while meeting dS1S2 > 0.2s.
Further, in above-mentioned 8th step, average heart rate in this section of heart sound data, calculation method is this section of heart sound
Time in data between all adjacent first heart sound S1 calculates heart rate, and counts the corresponding heart rate of its arithmetic mean of instantaneous value, calculates
Formula is r=60/ (pdS1S1), and wherein pdS1S1 is the arithmetic mean of instantaneous value in S1 period.
Second aspect, the present invention also provides a kind of heart sound graphic display methods, which is characterized in that its method is by the original heart
Treated that similar electrocardiogram is arranged above and below for sound figure and data, and shares a time shaft, the comparison up and down of characteristic point, at data
The comparable clearly mark heart sound S1 and heart sound S2 of similar electrocardiogram after reason, and original caardiophonogram can reflect heart sound S1 and heart sound
The original fluctuation situation of S2 merges the auxiliary diagnosis explained and be more advantageous to heart sound;Treated for the heart sound data similar to electrocardio
Figure refers to the datagraphic after data processing third step, it is mainly characterized by data after positive value conversion, no negative peak shape, class
It is similar to electrocardiogram.
The third aspect, the present invention also provides a kind of electronic auscultation systems, which is characterized in that it by auscultation head, pickup component,
Data record and processing component composition, the auscultation head are that stethoscopic sound obtains and conducting parts, the pickup component will
Electric signal is converted to by the sound that auscultation head transmits, the data record and processing component are heart sound waveform number through the invention
The terminal device of voice data record and processing, the auscultation head and pickup component are carried out according to processing and heart sound graphic display method
It is connected, the pickup component can be connected with data record and processing component by wiredly and/or wirelessly mode.
Further, the terminal device includes mobile phone, tablet computer and computer.
Heart sound wave data processing method of the invention inputs any parameter without user, can accurately and automatically identify
The feature of first heart sound and second heart sound and extract real-time heart rate.
Heart sound graphic display method of the invention, while treated that similar electrocardiogram is shown by original caardiophonogram and data
On same figure, a time shaft is shared, convenient for the comparison up and down of characteristic point, the both original fluctuation of the original caardiophonogram of observable
Situation, can also from the similar electrocardiogram after data processing clearly heart sound feature, realize intuitiveization of heart sound figure, obtain
Clear caardiophonogram is obtained, makes to diagnose more acurrate.
Electronic auscultation system of the invention, because it uses heart sound wave data processing of the invention and heart sound graphical display side
Method, therefore simultaneously extract real-time heart rate can be recorded in real time, show clear caardiophonogram.
The features of the present invention sees the detailed description of the drawings of the present case and following preferable embodiment and obtains clearly
Solution.
Detailed description of the invention
Fig. 1 is the structure principle chart of electronic auscultation device system of the invention;
Fig. 2 is the characteristic point schematic diagram in heart sound period of the invention;
Fig. 3 is heart sound data process flow diagram of the invention;
Fig. 4 is identification first heart sound S1 of the invention and second heart sound S2 flow chart;
Fig. 5 is sample heart sound initial data caardiophonogram of the invention;
Fig. 6 is caardiophonogram after sample heart sound data processing of the invention;
Fig. 7 is the peak dot figure of pickup heart sound waveform of the invention;
Fig. 8 is characteristic point parametric results table of the invention;
Fig. 9 is peak dot schematic diagram after characteristic point peak height screening of the invention;
Figure 10 is characteristic point peak height the selection result table of the invention;
Figure 11 is that characteristic point peak of the invention is peak-to-peak away from screening peak dot schematic diagram;
Figure 12 is that characteristic point peak of the invention is peak-to-peak away from the selection result table;
Figure 13 is heart sound characteristic point final result schematic diagram of the invention;
Figure 14 is heart sound characteristic point final result table of the invention;
Figure 15 is that heart sound data merges display figure to sample data of the invention before and after the processing.
Specific embodiment
In order to be easy to understand the technical means, the creative features, the aims and the efficiencies achieved by the present invention, tie below
Conjunction is specifically illustrating, and the present invention is further explained.
Embodiment
It is shown in Figure 1 which shows the schematic diagram of electronic auscultation system of the invention.As seen from the figure, the electronic auscultation
System mainly includes that auscultation head 101, pickup component 102, data record and processing component 103 are total to three parts composition.Wherein, described
Auscultation head 101 is the acquisition of stethoscopic sound and conducting parts;The pickup component 102 converts the sound transmitted by auscultation head
For electric signal;The data record and processing equipment 103 are that heart sound wave data processing through the invention and heart sound figure are aobvious
Show that method carries out the terminal device of voice data record and processing, terminal device includes mobile phone, tablet computer and computer;It is described
Auscultation head 101 is connected with pickup component 102;The pickup component 102 and data record and processing component 103 can be by wired
And/or wireless mode is connected.
It is shown in Figure 2 which shows the characteristic point schematic diagram in heart sound period.As seen from the figure, heart sound characteristic point includes the
One heart sound S1 and second heart sound S2, third heart sound S3 and fourth heart sound S4.
It is shown in Figure 3 which shows heart sound data process flow diagram of the invention.Data processing includes 8 steps.
It is shown in Figure 5, it is a heart sound primary data sample caardiophonogram, sample data length is 3 seconds.
It is shown in Figure 6, it is caardiophonogram after a sample heart sound data processing.
By taking Fig. 5 sample data as an example, heart sound data treatment process of the present invention is explained:
Step 1: choosing one section of heart sound data, as shown in figure 5, being sample initial data caardiophonogram, sample data is long here
Degree is 3 seconds, and the length that this method chooses one section of heart sound data is preferably selected between 3 to 10 seconds, and one side data volume was unlikely to
Influence calculating speed greatly.
Step 2: data bandpass filtering: frequency filtering range is 20-100Hz.
Step 3: data positive value conversion: data take absolute value.
Step 4: data smoothing filtering: filtering time width is 0.1 second.
It is handled by above four step data, caardiophonogram is as shown in Figure 6 after data processing.
Step 5: pick up the peak dot of heart sound waveform: empirical parameter can be used in peak-seeking parameter, and guarantee peak dot is effective heart sound
Characteristic point;As shown in fig. 7, showing in figure to pick up the peak dot figure of heart sound waveform, peak dot 18 are picked up altogether.
Step 6: calculate the characteristic point parameter in each heart sound period: position, peak height, it is peak-to-peak away from;Each point data is as shown in Figure 8
Characteristic point parametric results table.
Step 7: identification first heart sound S1 and second heart sound S2: being according to Fig.4, identification first heart sound S1 and second heart
Sound S2 process screened and identified to the characteristic point parameter in each heart sound period that step 6 obtains, the final sample data
4 S1 and 4 S2 are obtained, as shown in Figure 13 heart sound characteristic point final result schematic diagram;
Step 8: it counts the average heart rate in this section of heart sound data: the heart is calculated using the time between adjacent first heart sound S1
Rate and average heart rate.In this example sample data, using Figure 14 heart sound characteristic point final result table, sees in figure and share 4 S1 spies
Point is levied, average heart rate is 78 beats/min.
It is shown in Figure 4 which shows identification first heart sound S1 and second heart sound S2 flow chart, are heart sound numbers in the present invention
According to the flow chart of data processing figure of the step 7 " identifying first heart sound S1 and second heart sound S2 " in process flow.The data processing
Process includes 5 steps.It is every to being obtained to step 6 (using Fig. 5 as heart sound primary data sample) by taking this data processing as an example
The characteristic point parameter in a heart sound period is shown in Fig. 8 characteristic point parametric results table, is screened and identified, steps are as follows:
Step 1: being screened according to peak number, this example sample data is 3 seconds, maximum effectively peak number=time span (unit
Second)/60 ' 120 ' 4=24, fact characteristic point parametric results number is 18, does not have to delete redundant character point;
Step 2: screened according to peak value, in the characteristic point parametric results table (Fig. 8) of this example sample data, maximum peak height
Value=12431, maximum peak value/4=12431/4=3108 of effective peak height minimum value=reality, is less than effective peak height most for peak height
Small value rejects the characteristic point that the characteristic point that peak height is less than effective peak height minimum value is rejected, and characteristic point data is shown in Figure 10 spy after screening
Sign point peak height the selection result table;Fig. 9 show peak dot schematic diagram after the screening of characteristic point peak height.
Step 3: peak-to-peak away from screening according to peak, in this example sample data, to result of the characteristic point after peak height screens into
Row peak is peak-to-peak away from screening, No. 3 points in Figure 10 characteristic point peak height the selection result table is rejected, characteristic point data is shown in after being screened
Figure 12 characteristic point peak is peak-to-peak away from the selection result table;It is peak-to-peak away from peak dot schematic diagram after screening that Figure 11 show characteristic point peak.
Step 4: identification first heart sound: in this example sample data, using characteristic point peak it is peak-to-peak away from the selection result table (see figure
12) first heart sound identification, is carried out, result is Figure 14 heart sound characteristic point final result table, and it is final that Figure 13 show heart sound characteristic point
Result schematic diagram is seen in figure and shares 4 S1 characteristic points.
Step 5: identifying second heart sound: recognition methods is the subsequent peak S1 if not S1, can be labeled as S2, together
When meet dS1S2 > 0.2s.In this example sample data, after first heart sound identification, according to the principle of mark S2, S2 point is marked, most
Termination fruit is Figure 14 heart sound characteristic point final result table, and Figure 13 show heart sound characteristic point final result schematic diagram, sees in figure
Share 4 S2 characteristic points.
It is shown in Figure 4 which shows heart sound data merges display figure to sample data before and after the processing;Figure middle and upper part is divided into original
Beginning datagram is schemed after lower partial data processing, and the upper following figure shares time shaft, facilitates comparative analysis.
The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the above embodiments and description only describe this
The principle of invention, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes
Change and improvement all fall within the protetion scope of the claimed invention, the claimed scope of the invention by appended claims and its
Equivalent thereof.
Claims (9)
1. a kind of heart sound wave data processing method, which is characterized in that it includes the following steps:
(1) one section of heart sound data is chosen;
(2) data processing and conversion are carried out to heart sound Wave data, picks up the peak dot of heart sound waveform;
(3) first heart sound and second heart sound are identified, and counts the average heart rate in the segment data;
Identification first heart sound S1 and second heart sound S2 method include the following steps:
Step 1: being screened according to peak number, screening technique 1) maximum effectively peak number is first calculated, maximum effectively peak number calculates
Formula are as follows: maximum effectively peak number=time span/60 × 120 × 4;2) if feature point number is no more than maximum effectively peak
Number does not have to screening then;3) if necessary to screen, then first characteristic point parametric results are ranked up from high to low with peak value, so
It is afterwards more than the feature point deletion of effective number tail portion;
Step 2: being screened according to peak value, screening technique 1) extract the practical maximum peak value of characteristic point;2) effective peak is determined
High minimum value, maximum peak value/4 of effective peak height minimum value=reality;3) characteristic point that peak height is less than effective peak height minimum value is picked
It removes;
Step 3: peak-to-peak away from screening according to peak, screening technique is to carry out according to peak is peak-to-peak away from the principle that must not drop below 0.2 second
Screening, when two it is adjacent it is peak-to-peak away from less than 0.2 second when, the small characteristic point of peak height is rejected, it is peak-to-peak away from directly then to search peak again again
To meet peak it is peak-to-peak away from requirement;
Step 4: identification first heart sound: the spacing of first heart sound S1 and second heart sound S2 are with following time response: dS1S2 <
DS2S1, i.e., in a heart sound period, the time dS1S2 between first heart sound S1 to second heart sound S2 is less than second heart sound S2 and arrives
Time dS2S1 between first heart sound S1;According to the time response of first heart sound S1 and second heart sound S2 spacing, to characteristic point peak peak
Spacing the selection result table carries out first heart sound identification;
Step 5: identifying second heart sound: recognition methods is the subsequent peak first heart sound S1 if not first heart sound S1, that is, is marked
It is denoted as S2, while meeting dS1S2 > 0.2s.
2. heart sound wave data processing method as described in claim 1, which is characterized in that the processing of the heart sound Wave data
Specific step is as follows:
Step 1: choosing one section of heart sound data, and data are long to be not less than 3 seconds;
Step 2: data bandpass filtering: frequency filtering range is 20-100Hz;
Step 3: data positive value conversion: data take absolute value;
Step 4: data smoothing filtering;
Step 5: picking up the peak dot of heart sound waveform: peak-seeking parameter use experience parameter, and guarantee peak dot is effective heart sound characteristic point;
Step 6: calculate the characteristic point parameter in each heart sound period: position, peak height and it is peak-to-peak away from;
Step 7: identification first heart sound S1 and second heart sound S2: according to the time attribute of heart sound, come identify first heart sound S1 and
Second heart sound S2;
Step 8: count the average heart rate in this section of heart sound data: using the time between adjacent first heart sound S1 calculate heart rate and
Average heart rate.
3. heart sound wave data processing method as described in claim 1, which is characterized in that the length of the heart sound data is 3-
10 seconds, data volume is made to be unlikely to excessive and influence calculating speed.
4. heart sound wave data processing method as claimed in claim 2, which is characterized in that the characteristic point refers in a heart
First heart sound S1, second heart sound S2, third heart sound S3, fourth heart sound S4 in the sound period;The characteristic point parameter includes peak position
It sets, peak height, the peak-to-peak time interval away from adjacent first heart sound S1 in, caardiophonogram are expressed as dS1S1, a heart sound week in caardiophonogram
The time interval of first heart sound S1 and second heart sound S2 are expressed as dS1S2 in phase.
5. heart sound wave data processing method as claimed in claim 2, which is characterized in that in above-mentioned step 7, identify first heart
Sound S1 and second heart sound S2, the purpose is to handle heart sound data the characteristic point parameter in each heart sound period that preceding 6 step obtains, into
Row screening, rejects invalid peak dot, and finally automatically identify first heart sound S1 and second heart sound S2.
6. heart sound wave data processing method as claimed in claim 2, which is characterized in that in above-mentioned step 8, this section of heart
Average heart rate in sound data, calculation method are to calculate the time in this section of heart sound data between all adjacent first heart sound S1
Heart rate, and the corresponding heart rate of its arithmetic mean of instantaneous value is counted, calculation formula is r=60/ (pdS1S1), and wherein pdS1S1 is the S1 period
Arithmetic mean of instantaneous value.
7. a kind of heart sound graphic display method, which is characterized in that its method is that treated is similar by original caardiophonogram and data
Electrocardiogram is arranged above and below, and shares a time shaft, the comparison up and down of characteristic point, and the similar electrocardiogram after data processing is comparable
Heart sound S1 and heart sound S2 is clearly indicated, and original caardiophonogram can reflect the original fluctuation situation of heart sound S1 and heart sound S2, merge
Explain the auxiliary diagnosis for being more advantageous to heart sound;Treated that similar electrocardiogram refers to data such as claim 2 for the heart sound data
Described in step 3 treated datagraphic, it is characterized in that data, after positive value conversion, no negative peak shape is similar to electrocardio
Figure.
8. a kind of electronic auscultation system, which is characterized in that it is made of auscultation head, pickup component, data record and processing component,
The auscultation head is that stethoscopic sound obtains and conducting parts, the pickup component are converted to the sound transmitted by auscultation head
Electric signal, the data record and processing component are to pass through heart sound wave data processing method described in claim 1 and right
It is required that heart sound graphic display method described in 7 carries out the terminal device of voice data record and processing, the auscultation head and pickup
Component is connected, and the pickup component can be connected with data record and processing component by wiredly and/or wirelessly mode.
9. electronic auscultation system as claimed in claim 8, which is characterized in that the terminal device includes mobile phone and computer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610570124.2A CN106037792B (en) | 2016-07-20 | 2016-07-20 | A kind of heart sound wave data processing and graphic display method and its auscultation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610570124.2A CN106037792B (en) | 2016-07-20 | 2016-07-20 | A kind of heart sound wave data processing and graphic display method and its auscultation system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106037792A CN106037792A (en) | 2016-10-26 |
CN106037792B true CN106037792B (en) | 2019-11-26 |
Family
ID=57188470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610570124.2A Active CN106037792B (en) | 2016-07-20 | 2016-07-20 | A kind of heart sound wave data processing and graphic display method and its auscultation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106037792B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106551705B (en) * | 2016-11-09 | 2019-11-08 | 广州贝护佳医疗科技有限公司 | Processing, playback method and the system of fetal rhythm audio data and fetal rhythm spectrum data |
TWI644651B (en) * | 2017-11-17 | 2018-12-21 | 國立成功大學 | Signal synchronization device, as well as stethoscope, auscultation information output system and symptom diagnosis system capable of signal synchronization |
CN109009059B (en) * | 2018-09-11 | 2021-03-30 | 江苏鹿得医疗电子股份有限公司 | Heart rate calculation method based on heart sounds |
CN109998514B (en) * | 2019-04-25 | 2021-06-29 | 四川长虹电器股份有限公司 | Heart rate calculation method for irregular heart sound signals |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1850007A (en) * | 2006-05-16 | 2006-10-25 | 清华大学深圳研究生院 | Heart disease automatic classification system based on heart sound analysis and heart sound segmentation method |
CN101083938A (en) * | 2004-12-22 | 2007-12-05 | 大日本住友制药株式会社 | Electrocardiogram waveform correction display and electrocardiogram waveform correction display method |
CN101357076A (en) * | 2008-08-08 | 2009-02-04 | 重庆大学 | Remote heart reservation index analysis system and method |
CN102068244A (en) * | 2011-01-13 | 2011-05-25 | 南京邮电大学 | Storing, processing and displaying method for heart sound in intelligent mobile phone visual heart sound examination device |
CN102283670A (en) * | 2011-07-08 | 2011-12-21 | 西华大学 | Method and device for quantitatively analyzing heart sound signals |
CN102302373A (en) * | 2011-06-30 | 2012-01-04 | 重庆大学 | Method for detecting heart sound signal feature points based on multiplexing multi-resolution wavelet transformation |
CN104287772A (en) * | 2014-10-16 | 2015-01-21 | 任瑞 | Digital stethoscope |
-
2016
- 2016-07-20 CN CN201610570124.2A patent/CN106037792B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101083938A (en) * | 2004-12-22 | 2007-12-05 | 大日本住友制药株式会社 | Electrocardiogram waveform correction display and electrocardiogram waveform correction display method |
CN1850007A (en) * | 2006-05-16 | 2006-10-25 | 清华大学深圳研究生院 | Heart disease automatic classification system based on heart sound analysis and heart sound segmentation method |
CN101357076A (en) * | 2008-08-08 | 2009-02-04 | 重庆大学 | Remote heart reservation index analysis system and method |
CN102068244A (en) * | 2011-01-13 | 2011-05-25 | 南京邮电大学 | Storing, processing and displaying method for heart sound in intelligent mobile phone visual heart sound examination device |
CN102302373A (en) * | 2011-06-30 | 2012-01-04 | 重庆大学 | Method for detecting heart sound signal feature points based on multiplexing multi-resolution wavelet transformation |
CN102283670A (en) * | 2011-07-08 | 2011-12-21 | 西华大学 | Method and device for quantitatively analyzing heart sound signals |
CN104287772A (en) * | 2014-10-16 | 2015-01-21 | 任瑞 | Digital stethoscope |
Also Published As
Publication number | Publication date |
---|---|
CN106037792A (en) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10117635B2 (en) | Electronic acoustic stethoscope with ECG | |
US10362997B2 (en) | System and method of extraction, identification, marking and display of heart valve signals | |
CN106037792B (en) | A kind of heart sound wave data processing and graphic display method and its auscultation system | |
CN108647565A (en) | A kind of data preprocessing method classified to electrocardiosignal based on deep learning model | |
CN109276272A (en) | A kind of multifunctional intellectual stethoscope | |
CN107106118A (en) | The method for detecting dicrotic notch | |
Chan et al. | Respiratory rate estimation using u-net-based cascaded framework from electrocardiogram and seismocardiogram signals | |
RU2346653C2 (en) | Cardiovascular sound analysis method and system | |
CN102657525A (en) | Nondestructive detection system and method of cardiac function | |
CN111419208A (en) | Acceleration sensor-based unbound real-time heart rate monitoring method and system | |
EP3399905A1 (en) | System and method of extraction, identification, making and display of the heart valve signals | |
CN107049357B (en) | Heart sound signal analysis method, device and system | |
Luo | Portable Bluetooth visual electrical stethoscope research | |
Ghassemian et al. | Early detection of pediatric heart disease by automated spectral analysis of phonocardiogram | |
CN102920475A (en) | Gastrointestinal sound monitor system | |
Behbahani | A hybrid algorithm for heart sounds segmentation based on phonocardiogram | |
CN109431493A (en) | Wearable body surface potential acquisition device and method based on range segment separating weighting algorithm | |
CN206491813U (en) | The electronic auscultation device that a kind of novel visual is listened | |
CN212853501U (en) | Heart sound signal acquisition and analysis system | |
CN109346193A (en) | One kind being based on stethoscopic speech recognition equipment and its method | |
CN205107727U (en) | Visible stethoscope | |
JP7332723B2 (en) | System for detecting QRS complexes in an electrocardiogram (ECG) signal | |
CN107361786A (en) | The medical electronic stethoscope of more collection points | |
Monika et al. | Embedded Stethoscope for Real Time Diagnosis of Cardiovascular Diseases | |
CN206453785U (en) | A kind of smart electronicses stethoscope |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |