CN106361283A - Heart sound signal optimization method - Google Patents
Heart sound signal optimization method Download PDFInfo
- Publication number
- CN106361283A CN106361283A CN201610804749.0A CN201610804749A CN106361283A CN 106361283 A CN106361283 A CN 106361283A CN 201610804749 A CN201610804749 A CN 201610804749A CN 106361283 A CN106361283 A CN 106361283A
- Authority
- CN
- China
- Prior art keywords
- cardiechema signals
- optimization method
- frequency
- carried out
- heart sound
- 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/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/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
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Abstract
The invention discloses a heart sound signal optimization method, which comprises the following steps of 1, performing morphologic processing on a heart sound signal; 2, performing frequency domain filtering on the heart sound signal subjected to the morphologic processing. The method provided by the invention has the advantages that the method of combining the morphology and the frequency filtering is used; signal burs are removed through the morphological calculation; then, through the frequency filtering, signals with too far frequency deviation are removed; the heart sound signal can be better optimized.
Description
Technical field
The present invention relates to medical treatment ＆ health technical field is and in particular to a kind of cardiechema signals optimization method.
Background technology
The cardiechema signals medically using are a kind of amplitude signal, for simultaneous interference environment, right and wrong
Often faint signal, detects that in the case of this low signal-to-noise ratio faint cardiechema signals need noise jamming is suppressed.
When being mixed with noise in signal, usual way is that the signal of collection is transformed to frequency domain, according to useful signal and
Residing frequency range is different in a frequency domain for noise, by frequency filter, noise is removed, then carries out the reconstruct of signal and recover former
Signal.But useful signal and noise are often overlapping on frequency domain in practice.
Prior art is generally processed to cardiechema signals on frequency domain, using methods such as frequency filter, wavelet filterings,
Reach the purpose of abating noises.But the frequency range of noise and useful signal is overlapping, if retaining narrower frequency band, going
Also greatly cut down signal itself while making an uproar, caused cardiechema signals distortion, and according to a wider frequency band, rise not again
To the effect removing noise.
Content of the invention
Instant invention overcomes the deficiencies in the prior art, provide a kind of cardiechema signals optimization method.
Cardiechema signals are a kind of quasi-periodic signals, and the waveform in each cycle has a typical case and similar morphological characteristic, adopts
Can be very good to lift heart sound quality with the method that morphology is combined with frequency domain filtering.
The present invention is directed to the wave character of heart sound, carries out Morphological scale-space to cardiechema signals first, can preferably remove
Interference noise.The method carrying out frequency domain filtering again, can avoid Morphological scale-space to introduce new high-frequency noise.
For solving above-mentioned technical problem, the present invention employs the following technical solutions:
A kind of cardiechema signals optimization method, described method comprises the following steps:
Step one, Morphological scale-space is carried out to cardiechema signals；
Step 2, frequency domain filtering is carried out to the cardiechema signals after Morphological scale-space.
Further technical scheme is that described step one includes:
Step a, cardiechema signals are carried out with morphology opening operation process；
Step b, closing operation of mathematical morphology process is carried out to cardiechema signals；
Step c, take the average of step a and step b operation result as the cardiechema signals after Morphological scale-space.
Further technical scheme is that described step a includes:
Erosion operation is carried out to original cardiechema signals, then carries out dilation operation,
Wherein ,
Wherein, s is cardiechema signals, and f is structural element.
Further technical scheme is that described step b includes:
Dilation operation is carried out to original cardiechema signals, then carries out erosion operation,
Wherein ,
Wherein, s is cardiechema signals, and f is structural element.
Further technical scheme is that described step c includes y=[s f+s ο f]/2.
Further technical scheme is that to carry out frequency domain filtering described in described step 2 be to be carried out using low pass filter
Frequency domain filtering.
Further technical scheme is that described step 2 also includes: the cut-off frequency of setting wave filter is 600hz.
Compared with prior art, one of beneficial effect of the embodiment of the present invention is: the inventive method adopts morphology and frequency
The method that rate filtering combines, removes signal burr by morphology operations, then removes frequency departure too far through frequency filtering
Signal, cardiechema signals can be optimized well.
Specific embodiment
All features disclosed in this specification, or disclosed all methods or during step, except mutually exclusive
Feature and/or step beyond, all can combine by any way.
Any feature disclosed in this specification (including any accessory claim, summary), unless specifically stated otherwise,
Replaced by other alternative features equivalent or that there is similar purpose.I.e., unless specifically stated otherwise, each feature is a series of
One of equivalent or similar characteristics example.
With reference to embodiment, the specific embodiment of the present invention is described in detail.
According to one embodiment of present invention, the present embodiment discloses a kind of cardiechema signals optimization method, and the present embodiment is with one
As a example cardiechema signals:
Cardiechema signals are carried out Morphological scale-space by the first step
Choose flat-structure element f, morphology opening operation and closed operation are carried out respectively to cardiechema signals, takes both averages
As the cardiechema signals after Morphological scale-space.I.e. y=[s f+s ο f]/2, wherein
Specifically, in the present embodiment, the first step carries out Morphological scale-space to cardiechema signals and comprises the following steps:
Cardiechema signals are carried out morphology opening operation process by step a
Erosion operation is carried out to original cardiechema signals, then carries out dilation operation,Wherein ,
Wherein, s is cardiechema signals, and f is structural element；
Cardiechema signals are carried out closing operation of mathematical morphology process by step b
Dilation operation is carried out to original cardiechema signals, then carries out erosion operation,Wherein ,
Wherein, s is cardiechema signals, and f is structural element；
Step c, takes the result as cardiechema signals Morphological scale-space for the average of a and b operation result, and y=[s f+s ο f]/
2.
Cardiechema signals after Morphological scale-space are carried out frequency domain filtering by second step, and the present embodiment chooses Butterworth low pass
Wave filter, setting cut-off frequency is 600hz, and passband maximum attenuation is 3db, and minimum attenuation in stop band is 18db.
The method that the present embodiment is combined using morphology and frequency filtering, removes signal burr by morphology operations,
Remove frequency departure signal too far through frequency filtering again, cardiechema signals can be optimized well.
" embodiment ", " another embodiment ", " embodiment " of being spoken of in this manual etc., refers to combine
The specific features of this embodiment description, structure or feature are included at least one embodiment of the application generality description.
Multiple local appearance statement of the same race in the description is not necessarily to refer to same embodiment.Furthermore, it is understood that combining arbitrary
When individual embodiment describes specific features, structure or feature, to be advocated be to realize with reference to other embodiment this
Feature, structure or feature also fall within the scope of the present invention.
Although reference be made herein to invention has been described for the multiple explanatory embodiments invented, however, it is to be understood that this
Skilled person can be designed that a lot of other modifications and embodiment, and these modifications and embodiment will fall in the application
Within disclosed spirit and spirit.More specifically, in the range of disclosure claim, can be to theme group
Close the building block of layout and/or layout carries out multiple modifications and improvement.Except the modification that building block and/or layout are carried out
Outer with improving, to those skilled in the art, other purposes also will be apparent from.
Claims (7)
1. a kind of cardiechema signals optimization method it is characterised in that: described method comprises the following steps:
Step one, Morphological scale-space is carried out to cardiechema signals；
Step 2, frequency domain filtering is carried out to the cardiechema signals after Morphological scale-space.
2. cardiechema signals optimization method according to claim 1 is it is characterised in that described step one includes:
Step a, cardiechema signals are carried out with morphology opening operation process；
Step b, closing operation of mathematical morphology process is carried out to cardiechema signals；
Step c, take the average of step a and step b operation result as the cardiechema signals after Morphological scale-space.
3. cardiechema signals optimization method according to claim 2 is it is characterised in that described step a includes: original heart sound is believed
Number carry out erosion operation, then carry out dilation operation,Wherein ,
Wherein, s is cardiechema signals, and f is structural element.
4. cardiechema signals optimization method according to claim 3 is it is characterised in that described step b includes: to the original heart
Message number carries out dilation operation, then carries out erosion operation,Wherein ,
Wherein, s is cardiechema signals, and f is structural element.
5. cardiechema signals optimization method according to claim 4 is it is characterised in that described step c includes
6. cardiechema signals optimization method according to claim 1 is it is characterised in that enter line frequency described in described step 2
Domain filtering is to carry out frequency domain filtering using low pass filter.
7. cardiechema signals optimization method according to claim 6 is it is characterised in that described step 2 also includes: setting filter
The cut-off frequency of ripple device is 600hz.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610804749.0A CN106361283A (en) | 2016-09-06 | 2016-09-06 | Heart sound signal optimization method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610804749.0A CN106361283A (en) | 2016-09-06 | 2016-09-06 | Heart sound signal optimization method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106361283A true CN106361283A (en) | 2017-02-01 |
Family
ID=57900390
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610804749.0A Pending CN106361283A (en) | 2016-09-06 | 2016-09-06 | Heart sound signal optimization method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106361283A (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103190901A (en) * | 2013-04-01 | 2013-07-10 | 天津工业大学 | R wave detection algorithm based on extremum field mean mode decomposition and improved Hilbert enveloping |
| CN103405227A (en) * | 2013-08-02 | 2013-11-27 | 重庆邮电大学 | Double-layer morphological filter based electrocardiosignal preprocessing method |
| CN104182625A (en) * | 2014-08-15 | 2014-12-03 | 重庆邮电大学 | Electrocardiosignal denoising method based on morphology and EMD (empirical mode decomposition) wavelet threshold value |
| CN104382590A (en) * | 2014-12-11 | 2015-03-04 | 赖大坤 | Automatic shockable rhythm identification and classification method combined with electrocardio time-frequency domain feature analysis |
| WO2015169724A1 (en) * | 2014-05-09 | 2015-11-12 | Universiteit Gent | Detection of pulmonary vein isolation |
| CN105662391A (en) * | 2016-01-27 | 2016-06-15 | 东北大学 | Feature extraction and classification system and method for gastric magnetic signals |
| CN105740845A (en) * | 2016-03-02 | 2016-07-06 | 深圳竹信科技有限公司 | Method and system for filtering baseline drift based on single layer morphology |
| CN105899268A (en) * | 2015-06-23 | 2016-08-24 | 中国科学院深圳先进技术研究院 | GPU-based parallel electrocardiosignal analyzing method |
-
2016
- 2016-09-06 CN CN201610804749.0A patent/CN106361283A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103190901A (en) * | 2013-04-01 | 2013-07-10 | 天津工业大学 | R wave detection algorithm based on extremum field mean mode decomposition and improved Hilbert enveloping |
| CN103405227A (en) * | 2013-08-02 | 2013-11-27 | 重庆邮电大学 | Double-layer morphological filter based electrocardiosignal preprocessing method |
| WO2015169724A1 (en) * | 2014-05-09 | 2015-11-12 | Universiteit Gent | Detection of pulmonary vein isolation |
| CN104182625A (en) * | 2014-08-15 | 2014-12-03 | 重庆邮电大学 | Electrocardiosignal denoising method based on morphology and EMD (empirical mode decomposition) wavelet threshold value |
| CN104382590A (en) * | 2014-12-11 | 2015-03-04 | 赖大坤 | Automatic shockable rhythm identification and classification method combined with electrocardio time-frequency domain feature analysis |
| CN105899268A (en) * | 2015-06-23 | 2016-08-24 | 中国科学院深圳先进技术研究院 | GPU-based parallel electrocardiosignal analyzing method |
| CN105662391A (en) * | 2016-01-27 | 2016-06-15 | 东北大学 | Feature extraction and classification system and method for gastric magnetic signals |
| CN105740845A (en) * | 2016-03-02 | 2016-07-06 | 深圳竹信科技有限公司 | Method and system for filtering baseline drift based on single layer morphology |
Non-Patent Citations (2)
| Title |
|---|
| 张晓燕等: "《数字视频处理及应用》", 31 January 2014, 西安电子科技大学出版社 * |
| 王金亮: "心电信号的预处理及R波检测的研究", 《中国优秀硕士学位论文全文数据库信息科技I辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108158573B (en) | Electrocardiosignal noise reduction method based on adaptive threshold wavelet transformation | |
| CN108985179B (en) | Electric energy quality signal denoising method based on improved wavelet threshold function | |
| CN110659621A (en) | Combined noise reduction method based on variational modal decomposition and permutation entropy | |
| CN107184203A (en) | Electrocardiosignal Feature point recognition method based on adaptive set empirical mode decomposition | |
| CN107361745B (en) | Supervised sleep electroencephalogram and electrooculogram mixed signal stage interpretation method | |
| CN108051857B (en) | A kind of local similarity wavelet threshold denoising method based on empirical mode decomposition | |
| CN107464226A (en) | A kind of image de-noising method based on improvement two-dimensional empirical mode decomposition algorithm | |
| CN106570843A (en) | Adaptive wavelet threshold function image noise suppression method | |
| CN107361764B (en) | Method for rapidly extracting electrocardiosignal characteristic waveform R wave | |
| CN111160090B (en) | BCG signal noise reduction method and system | |
| CN104835127A (en) | Adaptive smooth filtering method | |
| CN105030232A (en) | Baseline drift correction method for electrocardiosignal | |
| Janušauskas et al. | Ensemble empirical mode decomposition based feature enhancement of cardio signals | |
| CN107689814A (en) | A kind of MSK band spread receivers narrow-band interference rejection method | |
| Liu et al. | ECG signal denoising based on morphological filtering | |
| Rakshit et al. | An improved EMD based ECG denoising method using adaptive switching mean filter | |
| CN104732498B (en) | A kind of thresholded image denoising method based on non-downsampling Contourlet conversion | |
| CN113576495B (en) | Motor imagery assessment method combined with EEG data quality | |
| CN110010144A (en) | Voice signals enhancement method and device | |
| CN106361283A (en) | Heart sound signal optimization method | |
| Bauernfeind et al. | Physiological noise removal from fNIRS signals | |
| CN103915102B (en) | Method for noise abatement of LFM underwater sound multi-path signals | |
| CN105719658B (en) | Wavelet packet voice de-noising method based on new threshold function table and adaptive threshold | |
| Heydari et al. | Adaptive wavelet technique for EEG de-noising | |
| US11660062B2 (en) | Method and system for recognizing crackles |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170201 |