CN105653880B - A kind of heart sounding field source tracing method based on heart sound - Google Patents
A kind of heart sounding field source tracing method based on heart sound Download PDFInfo
- Publication number
- CN105653880B CN105653880B CN201610120681.4A CN201610120681A CN105653880B CN 105653880 B CN105653880 B CN 105653880B CN 201610120681 A CN201610120681 A CN 201610120681A CN 105653880 B CN105653880 B CN 105653880B
- Authority
- CN
- China
- Prior art keywords
- heart
- sound
- field
- heart sound
- speech
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Abstract
The invention discloses a kind of methods using cardiechema signals reconfiguring heart sounding field.Heart sound inverse problem is the heart sound information obtained from torso model surface, unfaithful intention speech field energy distribution map is described, then converter technique is applied, obtain heart speech field energy distribution map, the sounding action parameter of interior heart sound sonification model is calculated with this, and the complete information of entire heart sounding field is reconstructed, realize a kind of heart sounding field source tracing method based on heart sound.Theoretical research foundation can be provided using the influence etc. of the heart tissue structure, transmittance process of cardiechema signals data conversed analysis generation heart sound, the generation process to reproduce the heart sound by being traced to the source by heart.
Description
Technical field
The present invention relates to a kind of methods using cardiechema signals reconfiguring heart sounding field.
Background technology
Heart sound inverse problem is the heart sound information obtained from torso model surface, description unfaithful intention speech field energy distribution
Then figure applies converter technique, obtain heart speech field energy distribution map, and the sounding that interior heart sound sonification model is calculated with this is dynamic
Make parameter, and reconstruct the complete information of entire heart sounding field, realizes a kind of heart sounding field source tracing method based on heart sound.
Being traced to the source by heart can be using the influence of the heart tissue structure, transmittance process of cardiechema signals data conversed analysis generation heart sound
Deng the generation process to reproduce the heart sound provides the basis of theoretical research.
Invention content
In view of the above-mentioned fact existing in the prior art, the present invention according to unfaithful intention speech field (body surface heart sound field) is counter in order to push away
The distribution of heart speech field establishes the plane between a kind of outer heart sound energy profile and interior heart sound energy profile to flat field
Reconstruction relationship.A kind of heart sounding field source tracing method based on heart sound of the present invention, if one synthesized by heart sounding chamber
The interior cardiechema signals in period can be described as：
Human body heart sound again through thoracic cavity echo model V radiation patterns R and the be cascaded as H of multichannel heart sound transducer Model B (Z)
The outer cardiechema signals that Transmission system is obtained can be described as：
Wherein Gu (n) is normalization shock response and its product of gain coefficient of system.G control heart sound gives great volume
It is small.It is then promising：
Here H (Z) is equivalent to a time varing filter stable in short-term.Its parameter ajIt is the human body device by generation heart sound
Official is determined, actually slowly varying at any time, but one relatively stablize the short time in, it is believed that be stablize it is constant
's.
If heart sounding chamber can be equivalent to several points sound source, each point source is sent out the heart sound sound wave of single-frequency,
So in the unfaithful intention speech field at thoracic cavity surface any pointIt can be expressed as the superposition of acoustic pressure caused by the point：
Wherein, TrIt is the surface area of heart sound transducer, rTIt is the position of heart sounding chamber point sound source, r is rTTo (x, y, z)
Distance, c is the velocity of sound, and ds is the area of heart sounding chamber point sound source.Damages of the human body heart sound Transmission system H (Z) to heart sound energy
Appraxia propagation attenuation factor D indicates.
Heart speech field zsFace and the sound field transformation relation on the faces unfaithful intention sound field z, according to Helmholtz equations, outer heart sound
At t moment any point, sound pressure is represented by sound field：
According to the convolution theorem that Fu Shi changes, the Fu Shi variations of the convolution of above-mentioned two formula are equal to multiplying for two formula Fu Shi variations
Product can be obtained on the faces unfaithful intention sound field z and heart speech field z further according to two-dimentional Fu Shi variation relationssFu Shi variation relations on face,
And inverse transformation is carried out, it can get heart speech field zsInverse transformation relational expression on face and the faces unfaithful intention sound field z：
Wherein
In order to solve inverse transformation, setting is needed to solve domain, the good grid of subdivision carries out finite element solving.According to real heart
Size dimension, in zsThe length of side of heart speech field prediction model is designed as the rectangular of 150mm × 100mm on plane formula
Shape, die body geometry cellular are defined as the square of 10mm, and die body is divided into 15 × 10 grid.
To the cardiechema signals S of a cycleTDivide by frame, then seeks each frame time tkInterior corresponding heart speech field
In zsAverage sound intensity distribution on plane formula.
The inverse transformation equation of above formula is subjected to discretization, and considers the acoustic pressure of each heart sound transducer acquisition to each
The reconstruct of place's heart sound field all contributes, and direct back projection method for reconstructing can be used, have heart speech field in zsOn plane formula
Average sound intensity is distributed as：
Wherein,For average sound intensity of the heart speech field at grid r of heart sound kth frame,For heart sound
The unfaithful intention speech field of kth frame is in grid rdThe average sound intensity at place, b (rd) be multichannel heart sound transducer Model B receptance function,It is torso model to the transmission coefficient of cardiechema signals, includes the influence of thoracic cavity echo model V and radiation patterns R.To asking
The heart speech obtained can be used by force color mapping to describe its distribution situation.
Although heart sound inverse transformation spatial resolution is determined by the size of geometry cellular, practical heart sound transducer contacts thoracic cavity
Area be bigger, as long as in principle cellular composition solid size be more than true heart sound transducer size surface area
.
As a diameter of 20mm of heart sound transducer, distribution can be described heart heart speech with 16 grids by force, such as attached
Shown in Fig. 2.There is the average sound intensity in kth frame to be：
Advantageous effect：It, can be first by cardiechema signals outside multichannel for the outer cardiechema signals that multichannel heart sound acquisition device obtains
Blind source separating is carried out, the outer cardiechema signals in the roads Shi Ge are mutual indepedent, eliminate influencing each other between them, then heart sound recycled to believe
Number frequency domain inverse transformation formula obtain the mutually independent interior cardiechema signals of multichannel, on this basis reality in heart sound sonification model
Reconstruct.
Description of the drawings
The heart speech field faces zs of Fig. 1 embodiment of the present invention and the sound field transformation relation schematic diagram on the faces unfaithful intention sound field z.
The heart based on heart sound of Fig. 2 embodiment of the present invention is traced to the source 16 grid sound pressure plane distribution of model heart speech field
Schematic diagram.
Specific implementation mode
The present invention is described in further detail below in conjunction with the drawings and specific embodiments.
In order to according to unfaithful intention speech field (body surface heart sound field) the anti-distribution for pushing away heart speech field, establish a kind of outer heart sound energy
The reconstruction relationship of plane between distribution map and interior heart sound energy profile to flat field；If one synthesized by heart sounding chamber
The interior cardiechema signals in period can be described as：
Human body heart sound again through thoracic cavity echo model V radiation patterns R and the be cascaded as H of multichannel heart sound transducer Model B (Z)
The outer cardiechema signals that Transmission system is obtained can be described as：
Wherein Gu (n) is normalization shock response and its product of gain coefficient of system.G control heart sound gives great volume
It is small.It is then promising：
Here H (Z) is equivalent to a time varing filter stable in short-term.Its parameter ajIt is the human body device by generation heart sound
Official is determined, actually slowly varying at any time, but one relatively stablize the short time in, it is believed that be stablize it is constant
's.
If heart sounding chamber can be equivalent to several points sound source, each point source is sent out the heart sound sound wave of single-frequency,
So in the unfaithful intention speech field at thoracic cavity surface any pointIt can be expressed as the superposition of acoustic pressure caused by the point：
Wherein, TrIt is the surface area of heart sound transducer, rTIt is the position of heart sounding chamber point sound source, r is rTTo (x, y, z)
Distance, c is the velocity of sound, and ds is the area of heart sounding chamber point sound source.Damages of the human body heart sound Transmission system H (Z) to heart sound energy
Appraxia propagation attenuation factor D indicates.
Heart speech field z as shown in Fig. 1sFace and the sound field transformation relation on the faces unfaithful intention sound field z, according to
Helmholtz equations are represented by t moment any point sound pressure on unfaithful intention speech field：
According to the convolution theorem that Fu Shi changes, the Fu Shi variations of the convolution of above-mentioned two formula are equal to multiplying for two formula Fu Shi variations
Product can be obtained on the faces unfaithful intention sound field z and heart speech field z further according to two-dimentional Fu Shi variation relationssFu Shi variation relations on face,
And inverse transformation is carried out, it can get heart speech field zsInverse transformation relational expression on face and the faces unfaithful intention sound field z：
Wherein
In order to solve inverse transformation, setting is needed to solve domain, the good grid of subdivision carries out finite element solving.According to real heart
Size dimension, in zsThe length of side of heart speech field prediction model is designed as the rectangular of 150mm × 100mm on plane formula
Shape, die body geometry cellular are defined as the square of 10mm, and die body is divided into 15 × 10 grid.
To the cardiechema signals S of a cycleTDivide by frame, then seeks each frame time tkInterior corresponding heart speech field
In zsAverage sound intensity distribution on plane formula.
The inverse transformation equation of (8) formula is subjected to discretization, and considers the acoustic pressure of each heart sound transducer acquisition to every
The reconstruct of heart sound field all contributes at one, and direct back projection method for reconstructing can be used, and has heart speech field in zsOn plane formula
Average sound intensity be distributed as：
Wherein,For average sound intensity of the heart speech field at grid r of heart sound kth frame,For heart sound
The unfaithful intention speech field of kth frame is in grid rdThe average sound intensity at place, b (rd) be multichannel heart sound transducer Model B receptance function,It is torso model to the transmission coefficient of cardiechema signals, includes the influence of thoracic cavity echo model V and radiation patterns R.To asking
The heart speech obtained can be used by force color mapping to describe its distribution situation.
Although heart sound inverse transformation spatial resolution is determined by the size of geometry cellular, practical heart sound transducer contacts thoracic cavity
Area be bigger, as long as in principle cellular composition solid size be more than true heart sound transducer size surface area
.
As a diameter of 20mm of heart sound transducer, distribution can be described heart heart speech with 16 grids by force, such as attached
Shown in Fig. 2.There is the average sound intensity in kth frame to be：
Frequency domain can be transformed by (2) formula for single channel cardiechema signals, had：
S (ω, m)=s ' (ω, m) H (ω, m) (12)
Log s (ω, m)=logs ' (ω, m)+logH (ω, m) (13)
Therefore, cardiechema signals are in the inverse transformation formula of frequency domain:
Logs ' (ω, m)=log s (ω, m)-logH (ω, m) (14)
It, can be by the blind source of the advanced row of cardiechema signals outside multichannel point for the outer cardiechema signals that multichannel heart sound acquisition device obtains
From cardiechema signals are mutual indepedent outside the roads Shi Ge, eliminate influencing each other between them, then recycle cardiechema signals in frequency domain
Inverse transformation formula obtain the mutually independent interior cardiechema signals of multichannel, on this basis reality in heart sound sonification model reconstruct.
The implementation of the present invention is described in detail above in conjunction with the drawings and specific embodiments, but the present invention is not
Be limited to the above embodiment, technical field those of ordinary skill within the scope of knowledge, can also not depart from
It is made a variety of changes under the premise of present inventive concept.
Claims (4)
1. a kind of heart sounding field source tracing method based on heart sound, it is characterised in that：
If the interior cardiechema signals of a cycle synthesized by heart sounding chamber are described as：
Human body heart sound transmission through thoracic cavity echo model V radiation patterns R and the be cascaded as H of multichannel heart sound transducer Model B (Z) again
The outer cardiechema signals that system is obtained are described as：
Wherein Gu (n) is normalization shock response and its product of gain coefficient of system；G controls the volume of heart sound；Then
It is promising：
Here H (Z) is equivalent to a time varing filter stable in short-term, its parameter ajIt is the human organ institute by generation heart sound
It determines, it is actually slowly varying at any time, but stablize in the short time relatively at one, stablize constant；
If heart sounding chamber can be equivalent to several points sound source, each point source is sent out the heart sound sound wave of single-frequency, then
Unfaithful intention speech field at thoracic cavity surface any pointIt is expressed as the superposition of acoustic pressure caused by the point：
Wherein, TrIt is the surface area of heart sound transducer, rTIt is the position of heart sounding chamber point sound source, r is rTTo (x, y, z) away from
From c is the velocity of sound, and ds is the area of heart sounding chamber point sound source；Human body heart sound Transmission system H (Z) uses the loss of heart sound energy
Propagation attenuation factor D indicates；
Heart speech field zsSound field transformation relation on face and the faces unfaithful intention sound field z, according to Helmholtz equations, on unfaithful intention speech field
At t moment any point, sound pressure is represented by：
According to the convolution theorem that Fu Shi changes, the Fu Shi variations of the convolution of above-mentioned two formula are equal to the product of two formula Fu Shi variations,
Further according to two-dimentional Fu Shi variation relations, obtain on the faces unfaithful intention sound field z and heart speech field zsFu Shi variation relations on face, and carry out
Inverse transformation, can get the faces unfaithful intention sound field z on and heart speech field zsInverse transformation equation on face：
Wherein
2. the heart sounding field source tracing method according to claim 1 based on heart sound, it is characterised in that：In order to solve inversion
It changes, setting is needed to solve domain, the good grid of subdivision carries out finite element solving；According to the size dimension of real heart, speech at heart
Field zsThe length of side of heart speech field prediction model is designed as the rectangle of a 150mm × 100mm, die body geometry member on plane formula
Born of the same parents are defined as the square of 10mm, and die body is divided into 15 × 10 grid；
To the cardiechema signals S of a cycleTDivide by frame, then seeks each frame time tkInterior corresponding heart speech field is in zs
Average sound intensity distribution on plane formula.
3. the heart sounding field source tracing method according to claim 1 based on heart sound, it is characterised in that：By inverse transformation equation
Discretization is carried out, and considers that the acoustic pressure that each heart sound transducer obtains contributes the reconstruct of everywhere heart sound field,
Using direct back projection method for reconstructing, there is heart speech field in zsAverage sound intensity on plane formula is distributed as：
Wherein,For average sound intensity of the heart speech field at grid r of heart sound kth frame,For heart sound kth frame
Unfaithful intention speech field in grid rdThe average sound intensity at place, b (rd) be multichannel heart sound transducer Model B receptance function,For
Torso model includes the influence of thoracic cavity echo model V and radiation patterns R to the transmission coefficient of cardiechema signals.
4. the heart sounding field source tracing method according to claim 1 based on heart sound, it is characterised in that：Work as heart sound transducer
A diameter of 20mm when, heart heart speech is distributed by force to be described with 16 grids, has the average sound intensity in kth frame to be：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610120681.4A CN105653880B (en) | 2016-03-03 | 2016-03-03 | A kind of heart sounding field source tracing method based on heart sound |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610120681.4A CN105653880B (en) | 2016-03-03 | 2016-03-03 | A kind of heart sounding field source tracing method based on heart sound |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105653880A CN105653880A (en) | 2016-06-08 |
| CN105653880B true CN105653880B (en) | 2018-08-24 |
Family
ID=56492743
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610120681.4A Active CN105653880B (en) | 2016-03-03 | 2016-03-03 | A kind of heart sounding field source tracing method based on heart sound |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105653880B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108392220A (en) * | 2018-02-01 | 2018-08-14 | 南京邮电大学 | A method of obtaining cardiechema signals derived components |
| CN114159091B (en) * | 2021-12-17 | 2023-10-13 | 大连理工大学 | Heart sound propagation relation detection system based on wearable sensor array |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102542864A (en) * | 2011-11-09 | 2012-07-04 | 南京邮电大学 | Achieving method of heart sound signal generator |
| CN103054572A (en) * | 2012-12-18 | 2013-04-24 | 太原理工大学 | Blind source separation method on basis of late potential signals of single-lead cardiac ventricle |
| US9008762B2 (en) * | 2009-11-03 | 2015-04-14 | Vivaquant Llc | Method and apparatus for identifying cardiac risk |
| CN105266788A (en) * | 2015-11-11 | 2016-01-27 | 四川长虹电器股份有限公司 | Fetal heart sound separating method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008123432A1 (en) * | 2007-03-30 | 2008-10-16 | Kyoto University | Device and method for acquiring a field by measurement |
-
2016
- 2016-03-03 CN CN201610120681.4A patent/CN105653880B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9008762B2 (en) * | 2009-11-03 | 2015-04-14 | Vivaquant Llc | Method and apparatus for identifying cardiac risk |
| CN102542864A (en) * | 2011-11-09 | 2012-07-04 | 南京邮电大学 | Achieving method of heart sound signal generator |
| CN103054572A (en) * | 2012-12-18 | 2013-04-24 | 太原理工大学 | Blind source separation method on basis of late potential signals of single-lead cardiac ventricle |
| CN105266788A (en) * | 2015-11-11 | 2016-01-27 | 四川长虹电器股份有限公司 | Fetal heart sound separating method |
Non-Patent Citations (2)
| Title |
|---|
| 《一种心音信号盲源分离方法》;孙科学 等；;《电子测量与仪器学报》;20120630;第26卷(第6期);第498-502页； * |
| 《一种改进型综合去噪算法在心音信号预处理上的研究》;陈新华 等；;《南京邮电大学学报》;20101231;第30卷(第6期);第96-106页； * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105653880A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ziegelwanger et al. | Mesh2hrtf: Open-source software package for the numerical calculation of head-related transfer functions | |
| Algazi et al. | Approximating the head-related transfer function using simple geometric models of the head and torso | |
| CN106125047B (en) | A kind of sound localization method under air-flow environment | |
| CN102158637B (en) | Method for denoising space self-adaption threshold video based on Surfacelet transform domain | |
| CN104122585A (en) | Seismic forward modeling method based on elastic wave field vector decomposition and low-rank decomposition | |
| CN109883532A (en) | A kind of identification of sound source and sound field forecasting procedure | |
| CN103616071B (en) | Patch near field acoustic holography-sound quality objective parameter distributed in three dimensions method for visualizing | |
| Kitić et al. | Physics-driven inverse problems made tractable with cosparse regularization | |
| CN105653880B (en) | A kind of heart sounding field source tracing method based on heart sound | |
| CN101847256A (en) | Image denoising method based on adaptive shear wave | |
| Wang et al. | Prediction of sound absorption of a periodic groove structure with rectangular profile | |
| CN105120419B (en) | Method and system for enhancing effect of multichannel system | |
| Ziegelwanger et al. | The three-dimensional quasi-periodic boundary element method: Implementation, evaluation, and use cases | |
| CN108154557B (en) | Triangular surface tiling method for silent area based on home environment | |
| Brick et al. | Fast direct solution of 3-D scattering problems via nonuniform grid-based matrix compression | |
| Xin et al. | Dynamic response and acoustic radiation of double-leaf metallic panel partition under sound excitation | |
| Grialou et al. | Characterization of surface impedance of vibro-acoustic subdomains with experimental measurements | |
| CN104318064A (en) | Three-dimensional head-related impulse response data compressing method based on canonical multi-decomposition | |
| CN102880590A (en) | Method for constructing non-split complete matching layer of second-order fluctuation equation | |
| CN106874606A (en) | A kind of near field acoustic holography sound field Forecasting Methodology based on Rayleigh second integrals | |
| Gao et al. | Multizone sound reproduction with adaptive control of scattering effects | |
| Betlehem et al. | Sound field of a directional source in a reverberant room | |
| Murata et al. | Sparse sound field decomposition with parametric dictionary learning for super-resolution recording and reproduction | |
| Krijnen et al. | openPSTD: the open source implementation of the Pseudo Spectral Time-Domain method | |
| Wang et al. | A Nautilus bionic multi-information fusion compressed-sensing acoustic imaging device |
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 |