CN105147327A - Generation method for multi-frequency shear wave - Google Patents
Generation method for multi-frequency shear wave Download PDFInfo
- Publication number
- CN105147327A CN105147327A CN201510323019.4A CN201510323019A CN105147327A CN 105147327 A CN105147327 A CN 105147327A CN 201510323019 A CN201510323019 A CN 201510323019A CN 105147327 A CN105147327 A CN 105147327A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- code
- tissue
- shearing wave
- 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.)
- Granted
Links
Abstract
The invention provides a generation method for multi-frequency shear waves. The method comprises: obtaining a frequency range [fg, fh] of shear waves in an interested tissue; selecting a first signal and a second signal, performing nonlinear computing on binary codes which represent the first signal and the second signal, to obtain a code excitation signal, and through the code excitation signal, driving a low-frequency vibrator to excite to generate multi-frequency shear waves in a to-be-detected tissue. In an elasticity detection technology applying a low-frequency vibrator excitation tissue, shear waves in multiple-frequency are excited and generated in the tissue, so as to solve an elastic coefficient and a viscosity coefficient of the tissue at the same time.
Description
[technical field]
The present invention relates to a kind of Medical Image Processing, particularly relate to a kind of production method of multi-frequency shearing wave.
[background technology]
The viscoelasticity information of biological tissue and the biological characteristics of tissue have and contact closely, and the change of viscoelasticity information can the Pathologic changes of characterising biological tissue well, for clinical diagnosis provides important reference frame.Ultrasound elastography imaging technology can carry out qualitative, quantitative measurement to the elasticity of tissue, viscosity.This technology encourages interested tissue regions generation deformation by certain method or produces vibration, and the information such as the deformation organized with medical ultrasonic detection or vibration, carry out qualitative or qualitative assessment to mechanical attributes such as the elastic modelling quantity organized or modulus of shearing.In ultrasonic tissue viscoelasticity detection technique, first needing to make tissue produce deformation or vibration by certain means, then detecting organizing deformation or vibration with medical ultrasonic.The method of four kinds of conventional excitation tissues has: natural excitation method, quasistatic motivational techniques, low-frequency vibration motivational techniques and acoustic radiation force motivational techniques.Natural motion such as the contraction of heartbeat, blood vessel or diastole etc. that namely natural excitation method is organized can drive surrounding tissue to move; Quasistatic motivational techniques utilize ultrasonic probe slowly to apply pressure directly to soft tissue, thus make tissue generation deformation; Low-frequency vibration motivational techniques, produce vibration with low-frequency vibrator at body surface, cause organization internal to vibrate; Acoustic radiation force motivational techniques, cause tissue vibration with focus supersonic at focus place, produce the shearing wave propagated in tissue.
The elastic characteristic of tissue, viscous characteristics and tissue be interior shears wave frequency, shear velocity of wave propagation becomes non-linear relation.When tissue vibration frequency is lower, organize main manifestations to go out elastic characteristic, viscosity can be ignored; When tissue vibration frequency is higher, main manifestations is organized to go out viscous characteristics.So, need to measure the shearing velocity of wave propagation obtaining multiple frequency content in tissue simultaneously and could solve the coefficient of elasticity and viscosity that obtain tissue simultaneously.
In the elastomeric check technology of existing application low-frequency vibrator excitation tissue vibration, what drive low-frequency vibrator is the signal of a single lower frequency.When shearing wave frequency is lower, affect the elastic characteristic mainly organized of shearing wave spread speed, so in the elastomeric check technology of existing application low-frequency vibrator excitation tissue, the elasticity number of tissue can only be obtained, have ignored the viscous characteristics of tissue, the information provided is not comprehensive.
[summary of the invention]
The production method that the technical problem to be solved in the present invention is to overcome multi-frequency shearing wave in prior art can only obtain the elasticity number of tissue, there is provided a kind of low-frequency vibrator excitation tissue of applying to produce the method including the shearing wave of multiple frequency content, thus coefficient of elasticity and the viscosity of tissue to be measured can be obtained simultaneously.
For solving the problems of the technologies described above, the present invention adopts following technical proposals:
A production method for multi-frequency shearing wave, comprises the steps:
Step S110: the frequency range [f obtaining shearing wave in tissue of interest
g, f
h];
Step S120: at described frequency range [f
g, f
h] in choose two signals, be designated as the first signal and secondary signal respectively, described first signal and secondary signal frequency are designated as f respectively
1, f
2;
Step S130: to representing that the binary code of described first signal and secondary signal carries out nonlinear operation and obtains encoded excitation signal, wherein, represent that the chip number of the binary code of described first signal and secondary signal is expressed as N1
code, N2
code;
Step S140: described encoded excitation signal drives low-frequency vibrator to produce multi-frequency shearing wave at underexcitation of organizing to be detected.
In certain embodiments, in step S130, described nonlinear operation is position XOR.
In certain embodiments, also comprise the steps: after completing steps S130
Change the chip number N1 of described first signal
code, secondary signal chip number N2
code, encoded excitation signal pulse number at least one, obtain and there is the encoded excitation signal of different spectral characteristic.
Adopt technique scheme, beneficial effect of the present invention is:
The production method of multi-frequency shearing wave provided by the invention, by obtaining the frequency range [f of shearing wave in tissue of interest
g, f
h]; Choose the first signal and secondary signal, and to representing that the binary code of described first signal and secondary signal carries out nonlinear operation and obtains encoded excitation signal, make in the elastomeric check technology of application low-frequency vibrator excitation tissue, the shearing wave that can underexcitation organized to produce multiple frequency, to try to achieve coefficient of elasticity and the viscosity of tissue simultaneously.
In addition, the production method of multi-frequency shearing wave provided by the invention, can participate in the first signal of computing and the chip number N1 of secondary signal by adjustment
code, N2
codevalue and the transmitted pulse number of pumping signal, the spectrum distribution adjusting the encoded excitation signal obtained neatly obtains the coded signal with different frequency scope and different harmonic value, thus can when applying, according to system needs, excitation produces the shearing wave with different frequency composition, to estimate tissue viscoelasticity's value.
[accompanying drawing explanation]
The flow chart of steps of the production method of a kind of multi-frequency shearing wave that Fig. 1 provides for one embodiment of the invention.
Fig. 2 is the time domain waveform and the amplitude spectrum waveform that participate in the time domain waveform of the first signal of nonlinear operation and the binary code of secondary signal, amplitude spectrum waveform and coded signal.
The chip number N2 of Fig. 3 for a change secondary signal
codethe frequency spectrum that obtains of numerical value, and the frequency spectrum of the coded signal obtained.
The chip number N1 of Fig. 4 for a change the first signal
codeafter, the amplitude spectrogram of the first signal.
Figure 5 shows that under the condition that secondary signal is identical, the first signal chip number N1
codedifference, the spectrum distribution of the coded signal obtained.
[detailed description of the invention]
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and the specific embodiments, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
Refer to Fig. 1, a kind of flow chart of steps applying the production method of the multi-frequency shearing wave of low-frequency vibrator excitation that Fig. 1 provides for the embodiment of the present invention, comprise the steps:
Step S110: the frequency range [f obtaining shearing wave in tissue of interest
g, f
h];
Step S120: at described frequency range [f
g, f
h] in choose two signals, be designated as the first signal and secondary signal respectively, described first signal and secondary signal frequency are designated as f respectively
1, f
2;
Preferably, [f within the scope of said frequencies
g, f
h] the frequency values f of the first selected signal and secondary signal
1, f
2lay respectively near low frequency, high frequency boundary, i.e. f
1>=f
g, f
2≤ f
h.
Step S130: to representing that the binary code of described first signal and secondary signal carries out nonlinear operation and obtains encoded excitation signal, wherein, represent that the chip number of the binary code of described first signal and secondary signal is expressed as N1
code, N2
code;
Preferably, in step S130, described nonlinear operation is position XOR.
Be appreciated that, described first signal and secondary signal are have certain length in time domain, it not endless, obvious spectral leakage can be there is in low frequency signal, contain multiple frequency content, the encoded excitation signal that described like this first signal and secondary signal obtain after carrying out nonlinear operation also contains multiple frequency content in interested frequency range.
Be appreciated that, by adopting above-mentioned nonlinear operation, when obtaining encoded excitation signal, the adjustment of parameter more flexibly, directly, and the frequency spectrum of the coded signal obtained each harmonic energy within the scope of frequency-of-interest increases exponentially with frequency to be increased, under the condition that gross energy is identical, the coded signal high frequency components energy obtained is higher, and the highdensity plywood wave energy that excitation tissue produces is also higher.
Preferably, in certain embodiments, also comprise the steps: after completing steps S130
Change the chip number N1 of described first signal
code, secondary signal chip number N2
code, encoded excitation signal pulse number at least one, obtain and there is the encoded excitation signal of different spectral characteristic.
Be appreciated that and work as described N1
codequantity, described N2
codequantity constant time, change the quantity of the driving pulse of described encoded excitation signal, be equivalent to the fundamental frequency changing driving pulse, described N1 can be made
codewith described N2
codefrequency shift, the corresponding change of coded frequency range signal can be caused equally.
Step S140: described encoded excitation signal drives low-frequency vibrator to produce multi-frequency shearing wave at underexcitation of organizing to be detected.
Be appreciated that and use encoded excitation signal to drive low-frequency vibrator, make in the elastomeric check technology of application low-frequency vibrator excitation tissue, the shearing wave that can underexcitation organized to produce multiple frequency, tries to achieve coefficient of elasticity and the viscosity of tissue simultaneously.
In addition, be with above-described embodiment difference: to representing that the binary code of described first signal and secondary signal adopts climbing method computing to obtain encoded excitation signal, the spectral shape of pumping signal as required arranges a window function, with this window function, binary signal is screened, finally obtain the coding waveforms meeting demand, the step due to other is consistent with above-described embodiment here to be repeated no more.
The production method of multi-frequency shearing wave provided by the invention, by obtaining the frequency range [f of shearing wave in tissue of interest
g, f
h]; Choose the first signal and secondary signal, and to representing that the binary code of described first signal and secondary signal carries out nonlinear operation and obtains encoded excitation signal, in the elastomeric check technology of application low-frequency vibrator excitation tissue, apply this encoded excitation signal and drive low-frequency vibrator excitation tissue, the shearing wave that can underexcitation organized to produce multiple frequency, to try to achieve coefficient of elasticity and the viscosity of tissue simultaneously; In addition, owing to can be participated in the N1 of computing by adjustment
codeand N2
codechip number and the driving pulse number of encoded excitation signal, adjust the spectrum distribution of the encoded excitation signal obtained flexibly, acquisition has different frequency scope, thus interested frequency range and harmonic value can be determined when applying, can derive and draw the value of relevant parameter, final acquisition meets the encoded excitation signal of different demand.
Set forth the present invention further by the following examples, these embodiments are only presented for purposes of illustration, do not limit the scope of the invention.
embodiment
Suppose that relative parameters setting is as follows, the frequency range of interested tissue vibration is set to [1,1000] Hz, the driving pulse number N launched in a sense cycle
shoots=512, the chip number participating in the first signal of nonlinear operation and the binary code of secondary signal is respectively N1
code=4, N2
code=128, the frequency of two corresponding signals is respectively 28.9Hz, 926Hz.To representing that the binary code step-by-step XOR of these two signals obtains encoded excitation signal.
Refer to time domain waveform and the amplitude spectrum waveform of time domain waveform, amplitude spectrum waveform and the coded signal that Fig. 2 is participation the first signal of nonlinear operation and the binary code of secondary signal, wherein, a, b is expressed as time domain waveform and the amplitude spectrum waveform of the binary code of the first signal, c, d are expressed as time domain waveform and the amplitude spectrum waveform of the binary code of secondary signal, and e, f are time domain waveform and the amplitude spectrum waveform of encoded excitation signal.
Can obtain from Fig. 2, the spectrum distribution of the encoded excitation signal that nonlinear operation obtains, gross energy in interested frequency range accounts for the major part of whole pumping signal gross energy, and interested [1,1000], in Hz frequency range, the energy of each harmonic components increases along with the increase exponentially of frequency.
Adjustment participates in the chip number in two signals of nonlinear operation, i.e. N1
code, N2
codenumerical value, change the frequency values of two signals, interested frequency range can be changed and number of harmonic wave within the scope of this.Under other parameter the same terms, change second signal N2
codechip number, make N
shoots=512, N1
code=4, N2
code=256, the highest energy of the coded signal that nonlinear operation obtains is got near 1852Hz, interested frequency range is become [1,2000] Hz.As shown in Figure 3, wherein (a), (b) is respectively and changes N2
codethe frequency spectrum that chip number obtains, and the frequency spectrum of the coded signal obtained.
If change N1
codemiddle chip number, then can change the number of harmonic wave within the scope of frequency-of-interest, when, N
shoots=512, N1
code=16, N2
codewhen=128, obtain first signal N1
codewith the frequency spectrum of encoded excitation signal as shown in Figure 4.
As shown in Figure 5, wherein, scheming (a) with figure (b) is at second signal N2
codechip number is identical, under launching the identical condition of driving pulse number, and N1
codechip number is different, and what obtain is code-excited in [1,1000] interested Hz frequency range, and the number of harmonic wave is different.
And when the first signal, the chip number N1 of secondary signal
code, N2
codenumerical value constant, change the number of driving pulse of launching, be equivalent to the fundamental frequency changing driving pulse, can N1 be made
code, N2
codethe frequency shift of (first, second signal), can cause the corresponding change of coded frequency range signal equally.
The above, it is only preferred embodiment of the present invention, not any pro forma restriction is done to the present invention, although the present invention discloses as above with preferred embodiment, but and be not used to limit the present invention, any those skilled in the art, do not departing within the scope of technical solution of the present invention, make a little change when the technology contents of above-mentioned announcement can be utilized or be modified to the Equivalent embodiments of equivalent variations, in every case be do not depart from technical solution of the present invention content, according to any simple modification that technical spirit of the present invention is done above embodiment, equivalent variations and modification, all still belong in the scope of technical solution of the present invention.
Claims (3)
1. a production method for multi-frequency shearing wave, is characterized in that, comprises the steps:
Step S110: the frequency range [f obtaining shearing wave in tissue of interest
g, f
h];
Step S120: at described frequency range [f
g, f
h] in choose two signals, be designated as the first signal and secondary signal respectively, described first signal and secondary signal frequency are designated as f respectively
1, f
2;
Step S130: to representing that the binary code of described first signal and secondary signal carries out nonlinear operation and obtains encoded excitation signal, wherein, represent that the chip number of the binary code of described first signal and secondary signal is expressed as N1
code, N2
code;
Step S140: described encoded excitation signal drives low-frequency vibrator to produce multi-frequency shearing wave at underexcitation of organizing to be detected.
2. the production method of multi-frequency shearing wave according to claim 1, is characterized in that, in step S130, described nonlinear operation is position XOR.
3. the production method of multi-frequency shearing wave according to claim 1, is characterized in that, also comprises the steps: after completing steps S130
Change the chip number N1 of described first signal
code, secondary signal chip number N2
code, encoded excitation signal pulse number at least one, obtain and there is the encoded excitation signal of different spectral characteristic.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510323019.4A CN105147327B (en) | 2015-06-12 | 2015-06-12 | A kind of production method of multi-frequency shearing wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510323019.4A CN105147327B (en) | 2015-06-12 | 2015-06-12 | A kind of production method of multi-frequency shearing wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105147327A true CN105147327A (en) | 2015-12-16 |
| CN105147327B CN105147327B (en) | 2018-03-09 |
Family
ID=54788639
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510323019.4A Expired - Fee Related CN105147327B (en) | 2015-06-12 | 2015-06-12 | A kind of production method of multi-frequency shearing wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105147327B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017121170A1 (en) * | 2016-01-11 | 2017-07-20 | 无锡海斯凯尔医学技术有限公司 | Tissue parameter detection method and system |
| CN107928703A (en) * | 2016-10-12 | 2018-04-20 | 甘肃农业大学 | A kind of biological information imaging device and method |
| CN108982674A (en) * | 2017-06-05 | 2018-12-11 | 郑翊 | Wave property is sheared using single driving pulse and with different in width and multiple driving pulses of shape measurement tissue |
| CN109996498A (en) * | 2016-11-16 | 2019-07-09 | 罗切斯特大学 | The reverberation shearing wave field estimation of subject characteristic |
| CN112638275A (en) * | 2018-06-27 | 2021-04-09 | 皇家飞利浦有限公司 | Shear wave detection of anatomical viscosity and associated devices, systems, and methods |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101803933A (en) * | 2009-02-13 | 2010-08-18 | 深圳市一体医疗科技有限公司 | Liver fibrosis detection device |
| WO2013091630A2 (en) * | 2011-12-21 | 2013-06-27 | Charité-Universitätsmedizin Berlin | Method for examining human or animal tissue |
-
2015
- 2015-06-12 CN CN201510323019.4A patent/CN105147327B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101803933A (en) * | 2009-02-13 | 2010-08-18 | 深圳市一体医疗科技有限公司 | Liver fibrosis detection device |
| WO2013091630A2 (en) * | 2011-12-21 | 2013-06-27 | Charité-Universitätsmedizin Berlin | Method for examining human or animal tissue |
Non-Patent Citations (3)
| Title |
|---|
| 曾斯宁 等: "超声弹性成像编码激励产生声辐射力方法的研究", 《2010年超声医学和医学超声论坛》 * |
| 林浩铭 等: "基于正交频率脉冲编码激励的软组织粘弹性定量测量", 《中国生物医学工程学报》 * |
| 陈思平 等: "超声编码激励产生声辐射力的数值模拟", 《深圳大学学报理工版》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017121170A1 (en) * | 2016-01-11 | 2017-07-20 | 无锡海斯凯尔医学技术有限公司 | Tissue parameter detection method and system |
| CN107928703A (en) * | 2016-10-12 | 2018-04-20 | 甘肃农业大学 | A kind of biological information imaging device and method |
| CN109996498A (en) * | 2016-11-16 | 2019-07-09 | 罗切斯特大学 | The reverberation shearing wave field estimation of subject characteristic |
| CN109996498B (en) * | 2016-11-16 | 2022-11-29 | 罗切斯特大学 | Reverberation shear wave field estimation of subject characteristics |
| CN108982674A (en) * | 2017-06-05 | 2018-12-11 | 郑翊 | Wave property is sheared using single driving pulse and with different in width and multiple driving pulses of shape measurement tissue |
| CN112638275A (en) * | 2018-06-27 | 2021-04-09 | 皇家飞利浦有限公司 | Shear wave detection of anatomical viscosity and associated devices, systems, and methods |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105147327B (en) | 2018-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105212968B (en) | Elastomeric check method and apparatus | |
| RU2688299C1 (en) | Method and device for determining viscoelastic medium parameter | |
| CN102641137B (en) | Measurement of viscoelastic property using amplitude-phase modulation | |
| CN105147327A (en) | Generation method for multi-frequency shear wave | |
| Shriki | Ultrasound physics | |
| CN104622502B (en) | Diagnostic ultrasound equipment and method | |
| US5474070A (en) | Method and apparatus for elastographic measurement and imaging | |
| CN105266851B (en) | The viscosity parameter detection method and equipment of viscoelastic medium | |
| CN109077754B (en) | Method and equipment for measuring tissue mechanical characteristic parameters | |
| CN104622509A (en) | Ultrasonic diagnostic apparatus and elastic evaluation method | |
| Urban | Production of acoustic radiation force using ultrasound: methods and applications | |
| CN109717899A (en) | Estimated in ultrasound medical imaging according to the tissue viscoelasticity of shear rate | |
| CN103006274A (en) | Method and system for ultrasonic detection of cornea viscoelasticity | |
| WO2017121170A1 (en) | Tissue parameter detection method and system | |
| CN104825195A (en) | Shear wave viscoelasticity imaging method and system | |
| CN102813533A (en) | Ultrasonic bone evaluation device based on acoustic radiation force and shear wave parameter detection method | |
| CN103267800B (en) | The ultrasonic fast imaging of cavitation microvesicle high s/n ratio and dynamic dimension distribution estimation method | |
| Ashir et al. | Skeletal Muscle assessment using quantitative ultrasound: A narrative review | |
| CN105636520A (en) | Ultrasound diagnosis device, ultrasound diagnosis method, and ultrasound diagnosis program | |
| CN105078512B (en) | A kind of method and device generating composite ultraphonic blood-stream image | |
| CN103142216B (en) | Method for calculating multilayer medium sound velocity based on photoacoustic imaging technology | |
| CN109498057A (en) | A kind of ultrasonic contrast imaging method, system, control equipment and storage medium | |
| JP5325670B2 (en) | Ultrasonic diagnostic equipment | |
| CN106659470A (en) | Ultrasonic diagnostic device | |
| EP3505927A1 (en) | Method and system for measuring a characteristic loop sensitivity for an acoustic transducer |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180309 Termination date: 20190612 |