CN108451498A - A kind of heart sound analysis system and method for fetal congenital heart disease - Google Patents
A kind of heart sound analysis system and method for fetal congenital heart disease Download PDFInfo
- Publication number
- CN108451498A CN108451498A CN201710090859.XA CN201710090859A CN108451498A CN 108451498 A CN108451498 A CN 108451498A CN 201710090859 A CN201710090859 A CN 201710090859A CN 108451498 A CN108451498 A CN 108451498A
- Authority
- CN
- China
- Prior art keywords
- heart sound
- fetal
- data
- heart
- frequency
- 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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4343—Pregnancy and labour monitoring, e.g. for labour onset detection
- A61B5/4362—Assessing foetal parameters
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- 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
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/02—Foetus
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Acoustics & Sound (AREA)
- Cardiology (AREA)
- Mathematical Physics (AREA)
- Pediatric Medicine (AREA)
- Pregnancy & Childbirth (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention discloses a kind of heart sound analysis method and system of fetal congenital heart disease.This approach includes the following steps:S1:Obtain fetal heart sound digital signal;S2:The heart sound digital signal is pre-processed, digital audio and video signals are obtained;S3:To carrying out frequency domain processing by the pretreated digital audio and video signals, frequency distribution characteristic, amplitude versus frequency characte and the power Spectral Estimation of fetal heart sound are obtained, pathology is obtained and judges data;S4:The pathology is judged into data, is shown by heart rate, heart sound curve or fetal rhythm state diagram.The heart sound analysis method and system of the fetal congenital heart disease of the present invention, can be provided for medical worker it is accurate, noninvasive, inexpensively, inexpensive, devoid of risk basis for estimation;Neonatal survival rate and the general level of the health are helped to improve, the universal of congenital heart disease pregnancy period screening is conducive to.
Description
Technical field
The present invention relates to medical instruments and digital signal processing technique field, more particularly to a kind of fetal congenital heart disease
Heart sound analysis system and method.
Background technology
Congenital heart disease (abbreviation congenital heart disease) be seriously affect fetus, baby's normal development heart organic disease,
Congenital heart disease refers in embryo development procedure, due to dissecting knot caused by the formation obstacle or dysplasia of heart and big blood vessel
Structure is abnormal, leads to the channel of blood flow that should be automatically closed after fetal birth, fails the situation being closed or channel of blood flow stenotic, from
And lead to the blood supply insufficiency of corporal parts, the normal development of fetus is seriously affected, or even die young.Congenital heart disease is congenital abnormal
Most common one kind in shape, accounts for about the 28% of various congenital abnormalities, and accounts for about the 0.4%~1% of life baby.
At present in fetal congenital heart disease diagnostic field, currently in use is more complex, relatively large expensive instrument;It is logical
It crosses when checking in 18-20 weeks pregnancy period, carries out outside four-dimensional ultrasound abdomen, the inspection of birth canal or enteral, and carry out the inspection of nuclear magnetic resonance
It looks into, it can be seen that whether fetus suffers from congenital heart disease, but the sensitivity of this kind of inspection is low, of high cost.But this kind of detection is adopted
Instrument exist in use higher risk, it is complicated for operation, be not easy the shortcomings of portable, use cost is excessively high, for universal
Property screening be difficult carry out.
Invention content
It is an object of the invention to for presently, there are defect, a kind of heart sound analysis of fetal congenital heart disease is provided
Method and system;Method through the invention obtains fetal heart sound signal to be had on this basis by Digital Signal Processing
Fetal heart sound since noise caused by heart organic disease is separated, is then determined the heart by effect ground according to the characteristic of noise
Dirty anatomic variation;This system and method for carrying out Digital Signal Processing for fetal heart sound signal, can be medicine work
Author provides accurate, noninvasive, cheap, inexpensive, devoid of risk basis for estimation;Be conducive to the screening of congenital heart disease pregnancy period
It is universal.
The present invention provides a kind of heart sound analysis method of fetal congenital heart disease, and this approach includes the following steps:
S1:Obtain fetal heart sound digital signal;
S2:The heart sound digital signal is pre-processed, digital audio and video signals are obtained;
S3:To carrying out frequency domain processing by the pretreated digital audio and video signals, the frequency point of fetal heart sound is obtained
Cloth characteristic, amplitude versus frequency characte and power Spectral Estimation obtain pathology and judge data;
S4:The pathology is judged into data, is shown by heart rate, heart sound curve or fetal rhythm state diagram.
The present invention also provides a kind of heart sound analysis system of fetal congenital heart disease, the system comprises:For obtaining
The fetal heart sound data receipt unit of the heart sound digital signal of fetus;For being pre-processed to fetal heart sound data receipt unit
Data pre-processing unit;For carrying out frequency domain processing to the cardiechema signals after the data prediction, the fetus heart is obtained
Frequency distribution characteristic, amplitude versus frequency characte and the power Spectral Estimation of sound obtain the data processing unit that pathology judges data;For inciting somebody to action
The pathology judges data, the display unit shown by heart rate, heart sound curve or fetal rhythm state diagram;And for whole
The main control unit that the control of a system and data calculate;Wherein, the fetal heart sound data receipt unit, data prediction list
Member, data processing unit and display unit are sequentially connected electrically, the main control unit respectively with fetal heart sound data receiver list
Member, data pre-processing unit, data processing unit, display unit electrical connection.
Beneficial effects of the present invention:Method through the invention obtains fetal heart sound signal, passes through number on this basis
Signal Pretreatment technology is then right effectively by fetal heart sound since noise caused by heart organic disease is separated
Frequency domain processing is carried out by the pretreated digital audio and video signals, frequency distribution characteristic, the amplitude-frequency for obtaining fetal heart sound are special
Property and power Spectral Estimation, obtain pathology and judge data;This system for carrying out Digital Signal Processing for fetal heart sound signal
And method, accurate, noninvasive, cheap, inexpensive, devoid of risk basis for estimation can be provided for medical worker;Have
Help improve neonatal survival rate and the general level of the health, is conducive to the universal of congenital heart disease pregnancy period screening.
Description of the drawings
Present invention will be further explained below with reference to the attached drawings and examples, in attached drawing:
Fig. 1 is a kind of flow chart of the heart sound analysis method of fetal congenital heart disease of the preferred embodiment of the present invention;
Fig. 2 is a kind of function module of the heart sound analysis system of fetal congenital heart disease of the preferred embodiment of the present invention
Figure;
Fig. 3 is the best finite impulse response filtering of the heart sound analysis method of fetal congenital heart disease through the invention
The design drawing of device;
Fig. 4 is the minimum Mean Square Error Linear FIR of the heart sound analysis method of fetal congenital heart disease through the invention
The design drawing of filter;
Fig. 5 is the heart sound analysis method application least mean-square error of the fetal congenital heart disease of the embodiment of the present invention 1
Cardiechema signals figure before linear FIR filtering;
Fig. 6 is the heart sound analysis method application least mean-square error of the fetal congenital heart disease of the embodiment of the present invention 1
The filtered cardiechema signals figures of linear FIR;
Fig. 7 is the normal fetus heart that the heart sound analysis method of the fetal congenital heart disease of the embodiment of the present invention 2 obtains
The distribution map of sound;
Fig. 8 is the doubtful atrioventricular valve that the heart sound analysis method of the fetal congenital heart disease of the embodiment of the present invention 2 obtains
The distribution map of the fetal heart sound of insufficiency;
Fig. 9 is the heart sound analysis method of the fetal congenital heart disease of the embodiment of the present invention 3 before carrying out frequency domain processing
Fetal heart sound distribution map;
Figure 10 is the heart sound analysis method of the fetal congenital heart disease of the embodiment of the present invention 3 after carrying out frequency domain processing
Fetal heart sound distribution map.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, not
For limiting the present invention.
Refering to fig. 1 and Fig. 2, the preferred embodiment of the present invention provide a kind of heart sound analysis system of fetal congenital heart disease
System, the analysis system includes sequentially connected fetal heart sound data receipt unit 1, data pre-processing unit 2, data processing list
Member 3 and display unit 4.Wherein, fetal heart sound data receipt unit 1 is placed at pregnant woman's lower abdomen, the heart sound for acquiring fetus
Signal.Data pre-processing unit 2 is amplified the fetal heart sound signal that fetal heart sound data receipt unit 1 is got, filters
With A/D data conversions, the digital audio and video signals of fetal heart sound are then sent to data processing unit 3 by being wirelessly transferred.Institute
It includes amplifier, filter and A/D data conversions to state data pre-processing unit 2, is respectively used to fetal heart sound data receiver list
The 1 fetal heart sound signal that gets of member is amplified, filters and A/D data conversions;Preferably, the filter is filtered using FIR
Wave device 21 carries out minimum Mean Square Error Linear FIR filtering to fetal heart sound signal.Data processing unit 3, wireless receiving is arrived
Digital audio and video signals are handled, including carry out FFT processing, work(by FFT processing modules 31 and power Spectral Estimation module 32
The processing such as rate Power estimation obtain pathology and judge data;Finally handling result is sent to display unit 4.The analysis system is also wrapped
Include main control unit 5, the main control list, 5 respectively with fetal heart sound data receipt unit 1, data pre-processing unit 2, data
Processing unit 3, display unit 4 are electrically connected, and for controlling the effect of these modules, and will pass through what power Spectral Estimation module 32 obtained
Pathology judges data, is shown in display unit 4 by heart rate, heart sound curve or fetal rhythm state diagram.The main control unit 5
It is calculated for realizing the control of whole system and data, and with analysis and drawing function, the equipment with these functions
For the present invention, such as can be digital signal processor commonly used in the art, host computer etc., the software journey of installation
Sequence, data operation program and mapping program can be realized by those skilled in the art according to the prior art.
Referring to Fig. 1, the preferred embodiment of the present invention provides a kind of heart sound analysis method of fetal congenital heart disease, packet
Include following steps:
S1:Obtain fetal heart sound digital signal;In the present embodiment, it using fetal heart sound data receipt unit 1, places pregnant
At woman's lower abdomen, the heart sound digital signal for acquiring fetus.
Preferably, heart sound digital signal is acquired at the lower abdomen of pregnant woman by infrared photoelectric sensing technology.Specifically, pregnant
Heart sound digital signal sensors are placed at the lower abdomen of woman, sending infrared light by heart sound digital signal sensors is used as detection light,
And receive to detect the reflected light of light, cardiechema signals are obtained by reflected light, to acquire the cardiechema signals of fetus.The infrared electro
Infrared photoelectric sensing technology known to those skilled in the art can be used in sensing technology, certainly other to acquire cardiechema signals
The prior art can be used for the present invention.
S2:The heart sound digital signal got is pre-processed, digital audio and video signals are obtained;Because of fetal heart sound
Acquisition when usually contain the interference from parent and the noise of measuring system itself, these are to correctly carrying out fetal heart sound point
Analysis can bring certain influence, thus need to pre-process the heart sound digital signal got.In this step, right
Heart sound digital signal is pre-processed, include the heart sound digital signal of fetus is amplified, FIR filtering, at A/D data conversions
Reason, obtains digital audio and video signals.
Thus, it is first link that the present invention is handled to carry out purification to the fetal heart sound of acquisition;Using adaptive
Linear filter.Minimum Mean Square Error Linear FIR filtering is carried out to the digital fetal heart sound signal acquired, is desirably to obtain preferably
True cardiechema signals, that is, purify.For heart sound frequency characteristic, using FIR low-pass filtering, cardiechema signals sample frequency 2~
3KHz, 0.8~1KHz of filter cutoff frequency.
Specifically, there was only sensor all the way when practical heart sound obtains, so being used as line by postponing digital cardiechema signals
Property prediction desired value time samples value.X (n) is the sampling time sample value of heart sound.N is the sequence of digital fetal heart sound signal
Length.In view of M=5, M are known as linear mean squared error estimator exponent number here for the real-time and convenience of calculation of system.
Y (n)=0≤n of x (n)≤N (1)
xk(n)=x (n+1-k) 1≤k≤M (2)
Because in linear prediction, the sampling of expected response and heart sound is the time samples value of same signal, so formula (1) is
The sampling time sequence of heart sound, the i.e. expected response of linear prediction;Formula (2) is the sampling time sequence after delay.
The present invention uses linear mean squared error estimator, such as (3) formula, is that one kind property grasped work is strong, is easy to mathematical computations
Method.
By to coefficient ck(n) determination so that (5) formula is minimum, and (4) formula ε (n) is error signal.
ε (n)=y (n)-Y (n) (4)
P (n)=E | ε (n) |2} (5)
By positive definite matrix theory, c can be calculatedk(n) value.
Above-mentioned result is applied in best FIR (finite impulse response filter) designs, as shown in Figure 3.
The output form of FIR filter is determined by following formula (6).
Here coefficient hk(n) impulse response of FIR filter is indicated.Equation group is solved by the above positive definite matrix theory, is obtained
To best coefficient ck(n).Also, h herek(n)=ck(n), minimum Mean Square Error Linear FIR filter, as shown in Figure 4.
Certainly other prior arts that can cardiechema signals be carried out with purified treatment can be used for the present invention.But this hair
The sef-adapting filter of bright use is different from conventional filter in advance it is to be understood that the spectrum distribution of signal, and it is only necessary to know that adopting
The sample value of sample signal.This adaptive-filtering is the sample value according to sampled signal, before the optimization of sample value statistical property
The filter put.
Using the above filtered fetal heart sound data, the frequecy characteristic of itself of fetal heart sound is remained, and by parent
Or the random disturbances that system itself generates are purified.Best letter under least mean-square error is provided for lower surface analysis
Number.
S3:To carrying out frequency domain processing by the pretreated digital audio and video signals, the frequency point of fetal heart sound is obtained
Cloth characteristic, amplitude versus frequency characte and power Spectral Estimation obtain pathology and judge data;
Specifically, in this step S3, data processing is carried out to digital audio and video signals, including carry out to digital audio and video signals
FFT processing, power Spectral Estimation obtain pathology and judge data.
Wherein, mathematical tools of the FFT (Fast Fourier Transform (FFT)) as time-frequency convert, its main function are exactly to solve letter
Number information that can not be obtained in the time domain.To solve the practical problem in various engineerings.
It is filtered to obtain the time samples value f (n) after time-domain filtering by above-mentioned FIR, FFT formulas (7) is passed through to f (n)
Frequency domain processing is carried out, uses 512 FFT transform of decimation in time here.
Wherein, f (n) is the time samples value after time-domain filtering, and e is exponential function, n=0,1,2 ..., 255;M=0,1,
2 ..., N-1, N=256;It is handled by FFT complex operations, because FFT operation results are mirror symmetries, obtains fetal heart sound letter
Number 256 domain complex handling result, i.e. R (n), I (n);P (n)=R is carried out to R (n), I (n)2(n)+I2(n) it calculates, n
=0,1,2 ... .255, obtain corresponding power spectrum.When carrying out FFT processing in order to reduce, time domain data blocks generation frequency leakage
Effect and resolution loss can carry out windowing process according to actual conditions to time domain fetal heart sound signal.It is commonly used in the art
Windowing process method be quarter window, Hanning window, hamming window, Gaussian window etc..
Why the present invention considers 512 FFT transform, is since it is considered that the component characteristic based on heart sound.Ordinary circumstance
Lower heart sound is made of first heart sound, second heart sound, third heart sound and fourth heart sound, in time prolonging in their having times
Late, they are the fairly close family's audio signals of frequency in frequency.If frequency resolution is too big, it would be possible that will produce
Spectral overlay causes the error in analysis.Since the frequency of sampling fetal heart sound is 2KHz, handled using 512 points, frequency discrimination
Rate is(Δ t=0.5ms is the time interval of sampling, and N=512, Δ F are frequency resolutions), Δ F=3.9Hz.
It is acceptable when as heart sound analysis by frequency resolution as experiment.In addition, in subsequent power Spectral Estimation
In, also make sequence length N sufficiently large as far as possible.
For there is the heart of fetus of anatomic variation, cause to flow through the blood of heart because narrow, insufficiency, defect etc. because
Element cannot be flowed by normal direction, and improper vibration is generated when flowing through these positions so as to cause blood, these vibrations are exactly to produce
The reason of sick rationality fetal heart sound.Pathologic fetal heart sound using fetal heart sound auscultation or carries out time-domain analysis in the pregnancy period to it
Processing, it is difficult to extract apparent feature, makes correct judgement.This case is handled by FFT turns the time-domain signal of fetal heart sound
After being changed to frequency-region signal, to the frequency distribution and amplitude versus frequency characte of its frequency-region signal be further analyzed and statistical significance under feature
Extraction is able to detect that blood stream through vibrating generated tiny dynamic change when heart in this way, and it is tool that these vibrations, which generate,
There is randomness, assurance in the time domain is very difficult, and can eaily carry out frequency spectrum, energy by frequency domain data
The comparison of the parameters such as amount and the statistical properties estimation.This is the analysis foundation that heart of fetus whether there is organic disease, and
Can be that fetus is made in the normal room of embryonic development period heart itself or ventricular septal defect and reliably differentiates foundation.
It is by frequency domain power during power Spectral Estimation is used as the estimation theory analyzed in fetal heart sound in this case
Spectral sequence P (n) is handled as a random sequence.Here inventor still assumes that the power spectrum sequence P (n) of cardiechema signals is one
A linear stationary random process.I.e. (8) formula is equal to 0:
According to a limited heart sound power spectrum sample value sequence P (n), (0≤n≤N-1) calculates stochastic variable { PnIt is flat
Equal estimated value.To assess the power spectrum data of fetal heart sound due to different cardioanatomy variations, system of the power spectrum in frequency domain
Meter learns feature.
If N is sufficiently large, the average value of power spectrum sample value sequence is shown in (9) formula.
The variance of power spectrum sample value sequence is shown in (10) formula.
εP=E [(P (n)-aP)2] (10)
The auto-covariance (auto-correlation function) of power spectrum sample value sequence is shown in (11) formula.
rxx(m)=[(P (n)-aP)(P*(n+m)-aP *)] (11)
These parameters can be used as our estimated service lives, for assess fetal heart sound frequency domain correlated characteristic.Especially
Analysis (for third heart sound and the more difficult judgement of fourth heart sound) to the correlated characteristic of first heart sound and second heart sound, to judging tire
The congenital organic disease of youngster's heart has very important meaning.Under normal circumstances, because the cardiac muscular tissue of fetus compares
Softness, so the heart sound frequency for thus vibrating generation is also relatively grown up low, but the heart rate of fetus is relatively grown up much higher, this is the fetus heart
The important feature of sound.And said relative to first heart sound and second heart sound because first heart sound is ventricular contraction sound, at this moment along with
Bicuspid valve, the closure of tricuspid valve, opening of aortic valve and rapid blood flow to aorta, generated heart sound frequency is higher,
Energy is larger, but the duration is shorter.Second heart sound is the diastole sound of ventricle, at this moment along with the opening of bicuspid valve, tricuspid valve
With the contraction in atrium because blood flow velocity is relatively low, generate heart sound frequency is relatively low, energy is also smaller, but duration
It is longer.
According to experiment, for the fetus with native heart organic disease, according to frequency of its cardiechema signals on frequency domain
Spectrum and energy-distributing feature and evaluation operators with statistical significance are obtained as a result, being obtained under normal circumstances with above-mentioned
Result it is then entirely different.And also differed with the variation of evaluation operators on frequency domain caused by different cardioanatomy variations,
Just for our screenings, there may be the types of congenital heart disease and congenital heart disease to provide basis for estimation with fetus is differentiated for this.
In practical application (9), (10), (11) parameter, because what they represented is feature under statistical significance, they are anti-
The statistics feature of fetal congenital heart disease different anatomic venereal disease change is reflected.aPAverage value reflects the average level of energy, for example,
In application example b, aPWant much higher in the case value with doubtful atrioventricular insufficiency (Fig. 8), it means that heart sound is deposited
In the higher noise of energy, and εPVariance yields is also bigger, this explanation compares concentration in high frequency section Energy distribution, in conjunction with
Fig. 8 is analyzed, doubtful atrioventricular insufficiency.
The present invention will be described in further detail by embodiment below.
Embodiment 1
Refering to Fig. 5 and Fig. 6, compared using the front and back cardiechema signals of minimum Mean Square Error Linear FIR filtering, Fig. 5 is filtering
Preceding cardiechema signals, Fig. 6 are filtered cardiechema signals;Fig. 6 is can be seen that from the comparison of Fig. 5 and Fig. 6 to carry out cardiechema signals
Purified treatment remains the frequecy characteristic of itself of fetal heart sound, provides the optimum signal under least mean-square error.
Embodiment 2
Refering to Fig. 7 and Fig. 8, fetal heart sound the power spectral energies of frequency domain distribution, to passing through the pretreated number
Word audio signal carries out frequency domain processing, obtains frequency distribution characteristic, amplitude versus frequency characte and the power Spectral Estimation of fetal heart sound, Fig. 7
For the distribution map of normal fetus heart sound;Fig. 8 is the distribution map of the fetal heart sound of doubtful atrioventricular insufficiency.It can thus be appreciated that passing through
Digital audio and video signals carry out frequency domain processing, obtain frequency distribution characteristic, amplitude versus frequency characte and the power Spectral Estimation of fetal heart sound, can
For normal and abnormal fetal heart sound respectively;Be conducive to the screening of congenital heart disease pregnancy period.
Embodiment 3
Refering to Fig. 9 and Figure 10, respectively for fetal heart sound in the comparison diagram of time domain and frequency domain, detection object is doubtful with room
The narrow fetus of interventricular septum.Fig. 9 is the filtered time domain waveform of fetal heart sound, if being difficult merely to carry with auscultation or time-domain analysis
For more information, the foundation that makes accurate judgment.Figure 10 is the spectrogram of the corresponding frequency domain of time domain waveform, the frequency of fetal heart sound
Spectrum is reflected when blood flows through room channel, causes velocity of blood flow to accelerate to hit with ventricle wall and atrioventricular valve because channel is narrow
It hits to form vortex, to reflect abundant harmonic component on frequency spectrum;It absolutely proves and digital audio and video signals is carried out at frequency domain
Reason, is more advantageous to the pregnancy period screening of congenital heart disease.
In conclusion method through the invention obtains fetal heart sound signal, located in advance by digital signal on this basis
Reason technology, effectively by fetal heart sound since noise caused by heart organic disease is separated, then to described in process
Pretreated digital audio and video signals carry out frequency domain processing, obtain frequency distribution characteristic, amplitude versus frequency characte and the work(of fetal heart sound
Rate Power estimation obtains pathology and judges data;This system and method for carrying out Digital Signal Processing for fetal heart sound signal, energy
Enough accurate, noninvasive, cheap, inexpensive, devoid of risk basis for estimation is provided for medical worker;It helps to improve new
The survival rate and the general level of the health of raw youngster is conducive to the universal of congenital heart disease pregnancy period screening.
It should be understood that for those of ordinary skills, it can be modified or changed according to the above description,
All these modifications and variations should all belong to the protection domain of appended claims of the present invention.
Claims (10)
1. a kind of heart sound analysis method of fetal congenital heart disease, which is characterized in that this approach includes the following steps:
S1:Obtain fetal heart sound digital signal;
S2:The heart sound digital signal is pre-processed, digital audio and video signals are obtained;
S3:To carrying out frequency domain processing by the pretreated digital audio and video signals, the frequency distribution for obtaining fetal heart sound is special
Property, amplitude versus frequency characte and power Spectral Estimation, obtain pathology and judge data;
S4:The pathology is judged into data, is shown by heart rate, heart sound curve or fetal rhythm state diagram.
2. the heart sound analysis method of fetal congenital heart disease according to claim 1, which is characterized in that in step S3
In, after the frequency distribution characteristic, amplitude versus frequency characte and the power Spectral Estimation that obtain fetal heart sound, by extracting each stage heart sound feature
And power spectrum, it obtains pathology and judges data.
3. the heart sound analysis method of fetal congenital heart disease according to claim 1, which is characterized in that in step S2
In, heart sound digital signal is pre-processed, include the heart sound digital signal of fetus is amplified, FIR filtering, A/D data
Conversion process obtains digital audio and video signals.
4. the heart sound analysis method of fetal congenital heart disease according to claim 1 or 2, which is characterized in that in step
In S3, data processing is carried out to digital audio and video signals, including FFT processing, power Spectral Estimation are carried out to digital audio and video signals, obtained
Pathology judges data.
5. the heart sound analysis method of fetal congenital heart disease according to claim 1 or 2, which is characterized in that described
In step S2, minimum Mean Square Error Linear FIR filtering is carried out to the fetal heart sound digital signal of acquisition, for fetal heart sound number
The frequency characteristic of signal, using low-pass filtering, cardiechema signals 2~3KHz of sample frequency, 0.8~1KHz of filter cutoff frequency.
6. the heart sound analysis method of fetal congenital heart disease according to claim 1, which is characterized in that in the step
In S3, real-time time extraction method FFT is carried out to digital audio and video signals, FFT complex operation processing is carried out by formula (7):
Wherein, f (n) is the time samples value after time-domain filtering, and e is exponential function, n=0,1,2 ..., 255;M=0,1,2 ...,
N-1, N=256;It is handled by FFT complex operations, obtains the handling result of 256 domain complex of fetal heart sound signal, i.e. R
(n)、I(n);
P (n)=R2(n)+I2(n), wherein the .255 of n=0,1,2 ... obtains corresponding power spectrum.
7. according to claim 1 or the heart sound analysis method of the fetal congenital heart disease, which is characterized in that in step S3
In, after carrying out FFT processing to digital audio and video signals, power Spectral Estimation is carried out by formula (8):
According to limited heart sound power spectrum sample value sequence P (n), wherein (0≤n≤N-1) calculates stochastic variable { PnBe averaged and estimate
Evaluation;
The average value of power spectrum sample value sequence is shown in (9) formula:
The variance of power spectrum sample value sequence is shown in (10) formula:
εP=E [(P (n)-aP)2] (10);
The auto-covariance of power spectrum sample value sequence is shown in (11) formula:
rxx(m)=[(P (n)-aP)(P*(n+m)-aP *)] (11)
Wherein, formula (9), (10), (11) parameter aP、εPAnd rxx(m) data are judged as pathology, for assessing fetal heart sound
The correlated characteristic of frequency domain.
8. a kind of heart sound analysis system of fetal congenital heart disease, which is characterized in that the system comprises:
Fetal heart sound data receipt unit for the heart sound digital signal for obtaining fetus;
For carrying out pretreated data pre-processing unit to fetal heart sound data receipt unit;
For carrying out frequency domain processing to the cardiechema signals after the data prediction, the frequency distribution for obtaining fetal heart sound is special
Property, amplitude versus frequency characte and power Spectral Estimation, obtain the data processing unit that pathology judges data;
For the pathology to be judged data, the display unit that is shown by heart rate, heart sound curve or fetal rhythm state diagram;
And the main control unit of the control and data calculating for whole system;
Wherein, the fetal heart sound data receipt unit, data pre-processing unit, data processing unit and display unit be successively
Electrical connection, the main control unit respectively with fetal heart sound data receipt unit, data pre-processing unit, data processing unit,
Display unit is electrically connected.
9. the heart sound analysis system of fetal congenital heart disease according to claim 8, which is characterized in that the data are pre-
Processing unit includes sequentially connected amplifier, filter and A/D data conversions;The filter is FIR filter.
10. the heart sound analysis system of fetal congenital heart disease according to claim 8, which is characterized in that the data
Processing unit includes sequentially connected FFT processing modules and power Spectral Estimation module.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710090859.XA CN108451498A (en) | 2017-02-20 | 2017-02-20 | A kind of heart sound analysis system and method for fetal congenital heart disease |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710090859.XA CN108451498A (en) | 2017-02-20 | 2017-02-20 | A kind of heart sound analysis system and method for fetal congenital heart disease |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108451498A true CN108451498A (en) | 2018-08-28 |
Family
ID=63221645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710090859.XA Pending CN108451498A (en) | 2017-02-20 | 2017-02-20 | A kind of heart sound analysis system and method for fetal congenital heart disease |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108451498A (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2675047Y (en) * | 2004-02-20 | 2005-02-02 | 肖青 | Nonlinear fetus heart rate tester |
US20090143650A1 (en) * | 2007-11-28 | 2009-06-04 | Regents Of The University Of Minnesota | Miniaturized, dermal-adhesive-based device for position-independent, non-invasive fetal monitoring |
CN101530331A (en) * | 2008-09-11 | 2009-09-16 | 兰州理工大学 | Heart sound signal detection and analysis system based on Labview development platform |
CN101790346A (en) * | 2007-07-24 | 2010-07-28 | 皇家飞利浦电子股份有限公司 | Method of monitoring a fetal heart rate |
CN101862201A (en) * | 2010-06-23 | 2010-10-20 | 南京邮电大学 | Method for representing cardiophonogram in intelligent cardiophone of mobile phone |
CN104274209A (en) * | 2014-08-22 | 2015-01-14 | 广东睿超电子科技有限公司 | Novel fetus-voice meter based on mobile intelligent terminal |
JP2015229030A (en) * | 2014-06-05 | 2015-12-21 | 国立大学法人福井大学 | Heartbeat state analysis device |
-
2017
- 2017-02-20 CN CN201710090859.XA patent/CN108451498A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2675047Y (en) * | 2004-02-20 | 2005-02-02 | 肖青 | Nonlinear fetus heart rate tester |
CN101790346A (en) * | 2007-07-24 | 2010-07-28 | 皇家飞利浦电子股份有限公司 | Method of monitoring a fetal heart rate |
US20090143650A1 (en) * | 2007-11-28 | 2009-06-04 | Regents Of The University Of Minnesota | Miniaturized, dermal-adhesive-based device for position-independent, non-invasive fetal monitoring |
CN101530331A (en) * | 2008-09-11 | 2009-09-16 | 兰州理工大学 | Heart sound signal detection and analysis system based on Labview development platform |
CN101862201A (en) * | 2010-06-23 | 2010-10-20 | 南京邮电大学 | Method for representing cardiophonogram in intelligent cardiophone of mobile phone |
JP2015229030A (en) * | 2014-06-05 | 2015-12-21 | 国立大学法人福井大学 | Heartbeat state analysis device |
CN104274209A (en) * | 2014-08-22 | 2015-01-14 | 广东睿超电子科技有限公司 | Novel fetus-voice meter based on mobile intelligent terminal |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5036857A (en) | Noninvasive diagnostic system for coronary artery disease | |
US5109863A (en) | Noninvasive diagnostic system for coronary artery disease | |
US5685317A (en) | Apparatus for measuring cardiac signals, using acoustic and ecg signals | |
RU2449730C2 (en) | Multiparameter classification of cardiovascular tones | |
CN103313662B (en) | System, the stethoscope of the risk of instruction coronary artery disease | |
US5727561A (en) | Method and apparatus for non-invasive detection and analysis of turbulent flow in a patient's blood vessels | |
CN103610460B (en) | A kind of Fetal ECG method for extracting signal based on self adaptation FLANN wave filter | |
JPH11506970A (en) | Method and apparatus for characterizing gastrointestinal sounds | |
CN104414632B (en) | Signal processing apparatus and signal processing method | |
CN104367344A (en) | Method and device for detecting instant heart rate of fetus on basis of Shanon envelope | |
CN203555751U (en) | Abdomen bowel sound analyzer system based on time-frequency analysis | |
TWI572327B (en) | Apparatus, computer program product and computer readable medium using audio signal for detection and determination of narrowing condition of fluid pipe | |
CN104873184B (en) | A kind of method and system for measuring heart rate | |
CN110313903B (en) | Pulse wave frequency domain characteristic parameter extraction method and device | |
CN107249466B (en) | Risk indication for coronary artery disease | |
CN109778485B (en) | Control method and system for clothes processing device | |
CN107392172B (en) | fetal heart sound processing method and processing device thereof | |
CN106175835B (en) | A kind of method of superonic spectrum complete cardiac cycle label | |
CN108451498A (en) | A kind of heart sound analysis system and method for fetal congenital heart disease | |
JP6103591B2 (en) | Auscultation heart sound signal processing method, auscultation heart sound signal processing apparatus, and program for processing auscultation heart sound signal | |
KR20130015368A (en) | Apparatus for diagnosing cardiac disorder by combining hidden markov model and murmur log-likelyhood | |
Hellicar et al. | An algorithm for the automatic analysis of signals from an oyster heart rate sensor | |
CN110881958A (en) | Pulse signal non-physiological signal removing method for traditional Chinese medicine pulse diagnosis instrument | |
CN109884483A (en) | Insulating tube type busbar shelf depreciation acoustics on-line monitoring method and device | |
Ari et al. | On a robust algorithm for heart sound segmentation |
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: 20180828 |
|
RJ01 | Rejection of invention patent application after publication |