CN117643480A - Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom - Google Patents
Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom Download PDFInfo
- Publication number
- CN117643480A CN117643480A CN202410115721.0A CN202410115721A CN117643480A CN 117643480 A CN117643480 A CN 117643480A CN 202410115721 A CN202410115721 A CN 202410115721A CN 117643480 A CN117643480 A CN 117643480A
- Authority
- CN
- China
- Prior art keywords
- human body
- body sample
- depth
- frequency
- attenuation coefficient
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000004075 alteration Effects 0.000 title claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 27
- 238000001514 detection method Methods 0.000 claims abstract description 25
- 230000005540 biological transmission Effects 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 7
- 238000002592 echocardiography Methods 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000003902 lesion Effects 0.000 description 2
- 210000005228 liver tissue Anatomy 0.000 description 2
- 208000004930 Fatty Liver Diseases 0.000 description 1
- 206010019708 Hepatic steatosis Diseases 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 208000010706 fatty liver disease Diseases 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 231100000240 steatosis hepatitis Toxicity 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention relates to the technical field of ultrasonic detection, in particular to a method for estimating an ultrasonic attenuation coefficient by utilizing aberration compensation and a reference phantom, which respectively utilizes a single emission frequency and a plurality of different emission frequencies to carry out ultrasonic detection on a human body sample to be detected and the reference phantom, calculates the center frequency of an echo signal at the depth of the beginning of an ultrasonic attenuation coefficient region to be estimated, selects one reference phantom and emission frequency with the center frequency closest to the human body sample to be detected from the center frequency, and finally utilizes the selected reference phantom and the emission frequency to estimate the attenuation coefficient of the human body sample to be detected according to a reference phantom method by utilizing the power spectrum of the echo signal measured by the human body sample to be detected; the beneficial effects are that: the center frequency of the reference body model and the center frequency of the human body sample to be detected at the detection starting depth are compensated, so that the compensation of the aberration of signals brought by different acoustic paths of the human body sample to be detected and the reference body model is realized, and the accuracy of estimating the ultrasonic attenuation coefficient of the human body sample to be detected is effectively improved.
Description
Technical Field
The invention relates to the technical field of ultrasonic detection, in particular to a method for estimating an ultrasonic attenuation coefficient by utilizing aberration compensation and a reference phantom.
Background
In clinical diagnosis of fatty liver, it is generally determined whether or not a lesion has occurred in liver tissue and the degree of the lesion by estimating the ultrasonic attenuation coefficient of the liver tissue; the current method for estimating the ultrasonic attenuation coefficient is to estimate by a reference body model method. Specifically, the power spectrums of echo signals of a human body sample to be detected and a reference phantom with known attenuation coefficient are measured under the same transmitting and receiving conditions, and the ratio of the power spectrums is only related to the difference between depth and the attenuation coefficient due to the same transmitting and receiving conditions, specifically, the attenuation coefficient of the human body sample to be detected is calculated by the following formula:wherein, the method comprises the steps of, wherein,is the attenuation coefficient of the human body sample to be measured,for the power spectrum of the echo signal of the human body sample to be measured at a certain transmission frequency,for the power spectrum of the echo signal of the reference phantom measured at a certain transmit frequency, when the transmit frequency is fixed,andonly in relation to the measurement depth z,taking logarithm of the ratio of the power spectrum of the echo signal of the human body sample to be measured and the power spectrum of the echo signal of the reference phantom, and measuring in a measuring areaDifferential value of depth z. When the reference body model method estimates the ultrasonic attenuation coefficient, because the acoustic paths of the reference body model and the human body sample to be measured are completely the same, and the ultrasonic wave can pass through transition areas of different sound speeds and attenuation coefficients of skin, bones, blood vessels, fat and the like before reaching a measurement area to be actually measured and estimated, the acoustic paths of the reference body model and the human body sample to be measured are identical, and certain deviation and difference exist in the acoustic paths, the ultrasonic attenuation coefficient estimated according to the conventional reference body model method has larger error.
Disclosure of Invention
The invention aims to provide a method for estimating an ultrasonic attenuation coefficient by using aberration compensation and a reference phantom, which can effectively compensate the deviation of a human body sample to be detected and the reference phantom caused by different acoustic paths and improve the accuracy of ultrasonic attenuation coefficient estimation.
In order to achieve the above purpose, the invention adopts the following technical scheme: method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom, human body sample to be detected and N different attenuation coefficients alpha 1 、α 2 、α 3 、……α N Respectively calculating the central frequency of echo signals of a human body sample to be detected and the reference body model at depth d, wherein the depth d is the depth of the beginning of a measurement area of the human body sample to be detected, and when the central frequency of the echo signals at the depth d is calculated by carrying out ultrasonic detection on the human body sample to be detected, adopting a single transmitting frequency f 0 Measuring and calculating echo signals after transmission, carrying out ultrasonic detection on N reference phantoms, and calculating the central frequency of echo at depth d, wherein a plurality of different transmission frequencies f are adopted for each reference phantoms 1 、f 2 、f 3 、……f m Measuring and calculating echo signals after transmission, and selecting one reference phantom alpha closest to the central frequency of the echo signals of the human body sample to be measured at the depth d from the central frequencies of the echoes of the N reference phantoms at the depth d obtained by measuring and calculating a plurality of different transmission frequencies i Transmitting frequencyRate f j The method comprises the steps of carrying out a first treatment on the surface of the Finally according to the formulaCalculating the attenuation coefficient of the human body sample to be measured, wherein,is the attenuation coefficient of the human body sample to be measured,for the human body sample to be measured, the emission frequency is f 0 Is a power spectrum of the echo signalAnd a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at timeTaking the logarithm of the ratio of (c) and taking the derivative value of the depth in the measurement area,is the attenuation coefficient of the reference phantom.
Specifically, in the human body sample to be tested and N different attenuation coefficients alpha 1 、α 2 、α 3 、……α N When the central frequencies of echo signals of a human body sample to be detected and the reference body mould at the depth d are calculated respectively, firstly, ultrasonic detection of different emission frequencies is adopted for the reference body moulds with N different attenuation coefficients, echo signal data of each reference body mould at the different emission frequencies are stored, when the ultrasonic detection is actually carried out, ultrasonic detection is carried out on the human body sample to be detected, the central frequency of the echo signals at the depth d is calculated, and the central frequencies of the echo signals of each reference body mould at the depth d at the different emission frequencies are calculated from the stored echo signal data.
Specifically, when calculating the center frequency of the echo signals of the human body sample to be measured and the reference phantom at the depth d, a centroid method is adopted for calculation, and the following formula is specifically adopted for calculationWherein, the method comprises the steps of, wherein,as a result of the center frequency,for the frequency of the echo signal,is the spectrum of the echo signal at depth d.
Specifically, when the attenuation coefficient of the human body sample to be measured is calculated, the human body sample to be measured is first transmitted at the frequency f 0 Is a power spectrum of the echo signal of (a)And a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at timeAnd (5) performing noise reduction processing, and calculating an attenuation coefficient by using a formula.
The invention has the beneficial effects that: the method comprises the steps of respectively carrying out ultrasonic detection on a human body sample to be detected and a reference phantom by utilizing a single emission frequency and a plurality of different emission frequencies, calculating the center frequency of an echo signal at the beginning depth of an ultrasonic attenuation coefficient region to be estimated, selecting one reference phantom with the center frequency closest to the human body sample to be detected and the corresponding emission frequency from the center frequency, and finally estimating the attenuation coefficient of the human body sample to be detected by utilizing the selected reference phantom and the power spectrum of the echo signal measured by the corresponding emission frequency and the human body sample to be detected according to a reference phantom method, thereby realizing the compensation of the aberration of signals brought by different acoustic paths of the human body sample to be detected and the reference phantom and effectively improving the accuracy of estimating the ultrasonic attenuation coefficient of the human body sample to be detected.
Drawings
Fig. 1 is a schematic diagram of the structure of a reference phantom and a human body sample to be tested and the acoustic path during ultrasonic detection in the embodiment.
Detailed Description
Embodiment 1, a method for estimating an ultrasonic attenuation coefficient using aberration compensation and a reference phantom, a human sample to be measured and N different attenuation coefficients α 1 、α 2 、α 3 、……α N Respectively calculating the central frequency of echo signals of a human body sample to be detected and the reference body model at depth d, wherein the depth d is the depth of the beginning of a measurement area of the human body sample to be detected, and when the central frequency of the echo signals at the depth d is calculated by carrying out ultrasonic detection on the human body sample to be detected, adopting a single transmitting frequency f 0 Measuring and calculating echo signals after transmission, carrying out ultrasonic detection on N reference phantoms, and calculating the central frequency of echo at depth d, wherein a plurality of different transmission frequencies f are adopted for each reference phantoms 1 、f 2 、f 3 、……f m Measuring and calculating echo signals after transmission, and selecting one reference phantom alpha closest to the central frequency of the echo signals of the human body sample to be measured at the depth d from the central frequencies of the echoes of the N reference phantoms at the depth d obtained by measuring and calculating a plurality of different transmission frequencies i Transmission frequency f j The method comprises the steps of carrying out a first treatment on the surface of the Finally according to the formulaCalculating the attenuation coefficient of the human body sample to be measured, wherein,is the attenuation coefficient of the human body sample to be measured,for the human body sample to be measured, the emission frequency is f 0 Is a power spectrum of the echo signal of (a)And a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at timeTaking the logarithm of the ratio of (c) and taking the derivative value of the depth in the measurement area,is the attenuation coefficient of the reference phantom. It should be noted that the power spectrum is affectedAndcomprising the transmission frequencyAnd a detection depth z, the power spectrum being related only to the detection depth when the emission frequency is fixed. In this embodiment, by selecting a reference phantom with a center frequency closest to the human body sample to be measured at the depth d and a corresponding emission frequency, aberration of signals caused by different acoustic paths of the human body sample to be measured and the reference phantom can be compensated, so that accuracy of estimating the ultrasonic attenuation coefficient is improved.
In the estimation method of the embodiment, the reference body model and the emission frequency closest to the central frequency of the echo signal of the human body sample to be measured at the depth d can be obtained more accurately by adopting a plurality of different emission frequencies to measure and calculate the central frequencies of the echo signals of the reference body models with N different attenuation coefficients, so that the efficiency of estimating the ultrasonic attenuation coefficients of the human body sample to be measured is improved.
Specifically, in the human body sample to be tested and N different attenuation coefficients alpha 1 、α 2 、α 3 、……α N When the central frequencies of echo signals of a human body sample to be detected and the reference body mould at the depth d are calculated respectively, firstly, ultrasonic detection of different emission frequencies is adopted for the reference body moulds with N different attenuation coefficients, echo signal data of each reference body mould at the different emission frequencies are stored, when the ultrasonic detection is actually carried out, ultrasonic detection is carried out on the human body sample to be detected, the central frequency of the echo signals at the depth d is calculated, and the central frequencies of the echo signals of each reference body mould at the depth d at the different emission frequencies are calculated from the stored echo signal data. By first generating a plurality of different patterns for N reference phantomsThe radio frequency carries out ultrasonic detection and saves echo signal data, and when the attenuation coefficient of the human body sample to be detected is actually measured, the saved echo signal data can be directly called, so that the data acquisition time is shortened, and the detection calculation efficiency is improved.
Specifically, when calculating the center frequency of the echo signals of the human body sample to be measured and the reference phantom at the depth d, the center-of-mass method is adopted for calculation, and the following formula is specifically adopted for calculationWherein, the method comprises the steps of, wherein,as a result of the center frequency,for the frequency of the echo signal,is the spectrum of the echo signal at depth d. In addition, the center frequency of the echo signal may be calculated or estimated by a method conventionally used in the art, in addition to the centroid method.
Specifically, in this embodiment, when calculating the attenuation coefficient of the human body sample to be measured, the human body sample to be measured is first transmitted at a frequency f 0 Is a power spectrum of the echo signal of (a)And a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at timeAnd (5) performing noise reduction processing, and calculating an attenuation coefficient by using a formula. The accuracy of the estimation of the ultrasonic attenuation coefficient of the human body sample to be detected can be further improved by carrying out noise reduction treatment on the power spectrum of the echo signal; specifically, the noise reduction processing of the power spectrum can be performed by adopting a filtering mode, and other conventionally used noise reduction processing modes can also be adopted.
Of course, the above embodiments are only preferred embodiments of the present invention, and the scope of the present invention is not limited thereto, so that all equivalent modifications made in the principles of the present invention are included in the scope of the present invention.
Claims (4)
1. A method for estimating an ultrasonic attenuation coefficient using aberration compensation and a reference phantom, characterized by: for human body sample to be measured and N different attenuation coefficients alpha 1 、α 2 、α 3 、……α N Respectively calculating the central frequency of echo signals of a human body sample to be detected and the reference body model at depth d, wherein the depth d is the depth of the beginning of a measurement area of the human body sample to be detected, and when the central frequency of the echo signals at the depth d is calculated by carrying out ultrasonic detection on the human body sample to be detected, adopting a single transmitting frequency f 0 Measuring and calculating echo signals after transmission, carrying out ultrasonic detection on N reference phantoms, and calculating the central frequency of echo at depth d, wherein a plurality of different transmission frequencies f are adopted for each reference phantoms 1 、f 2 、f 3 、……f m Measuring and calculating echo signals after transmission, and selecting one reference phantom alpha closest to the central frequency of the echo signals of the human body sample to be measured at the depth d from the central frequencies of the echoes of the N reference phantoms at the depth d obtained by measuring and calculating a plurality of different transmission frequencies i Transmission frequency f j The method comprises the steps of carrying out a first treatment on the surface of the Finally according to the formulaCalculating attenuation coefficient of human body sample to be measured, wherein +_>For the attenuation coefficient of the human body sample to be measured, +.>For the human body sample to be measured, the emission frequency is f 0 Power spectrum of time echo signal->And a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at timeTaking the logarithm of the ratio of (2) and taking the differential value of the measured area to the depth,/>Is the attenuation coefficient of the reference phantom.
2. A method of estimating an ultrasound attenuation coefficient using aberration compensation and reference phantom as set forth in claim 1, wherein: in the human body sample to be tested and N different attenuation coefficients alpha 1 、α 2 、α 3 、……α N When the central frequencies of echo signals of a human body sample to be detected and the reference body mould at the depth d are calculated respectively, firstly, ultrasonic detection of different emission frequencies is adopted for the reference body moulds with N different attenuation coefficients, echo signal data of each reference body mould at the different emission frequencies are stored, when the ultrasonic detection is actually carried out, ultrasonic detection is carried out on the human body sample to be detected, the central frequency of the echo signals at the depth d is calculated, and the central frequencies of the echo signals of each reference body mould at the depth d at the different emission frequencies are calculated from the stored echo signal data.
3. A method of estimating an ultrasound attenuation coefficient using aberration compensation and reference phantom as set forth in claim 1, wherein: when calculating the center frequency of the echo signals of the human body sample to be measured and the reference phantom at the depth d, calculating by adopting a centroid method, and specifically adopting the following formula to calculateWherein->For the center frequency +.>For the frequency of the echo signal +.>Is the spectrum of the echo signal at depth d.
4. A method of estimating an ultrasound attenuation coefficient using aberration compensation and reference phantom as set forth in claim 1, wherein: when the attenuation coefficient of the human body sample to be measured is calculated, the human body sample to be measured is firstly transmitted at the frequency f 0 Is a power spectrum of the echo signal of (a)And a reference phantom at a transmit frequency f j Power spectrum of echo signal measured at the time +.>And (5) performing noise reduction processing, and calculating an attenuation coefficient by using a formula.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410115721.0A CN117643480B (en) | 2024-01-29 | Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410115721.0A CN117643480B (en) | 2024-01-29 | Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117643480A true CN117643480A (en) | 2024-03-05 |
CN117643480B CN117643480B (en) | 2024-07-02 |
Family
ID=
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104116524A (en) * | 2014-01-21 | 2014-10-29 | 深圳市一体医疗科技有限公司 | Ultrasonic attenuation coefficient compensation system and liver fat detection system |
US20200146656A1 (en) * | 2017-06-23 | 2020-05-14 | Mayo Foundation For Medical Education And Research | Methods for ultrasound system independent attenuation coefficient estimation |
CN111466951A (en) * | 2020-04-15 | 2020-07-31 | 深圳开立生物医疗科技股份有限公司 | Method and device for generating ultrasonic attenuation image, ultrasonic equipment and storage medium |
WO2023274763A1 (en) * | 2021-06-29 | 2023-01-05 | Koninklijke Philips N.V. | Ultrasonic quantification of acoustic attenuation coefficient in the presence of elevation aperture blockage |
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104116524A (en) * | 2014-01-21 | 2014-10-29 | 深圳市一体医疗科技有限公司 | Ultrasonic attenuation coefficient compensation system and liver fat detection system |
US20200146656A1 (en) * | 2017-06-23 | 2020-05-14 | Mayo Foundation For Medical Education And Research | Methods for ultrasound system independent attenuation coefficient estimation |
CN111466951A (en) * | 2020-04-15 | 2020-07-31 | 深圳开立生物医疗科技股份有限公司 | Method and device for generating ultrasonic attenuation image, ultrasonic equipment and storage medium |
WO2023274763A1 (en) * | 2021-06-29 | 2023-01-05 | Koninklijke Philips N.V. | Ultrasonic quantification of acoustic attenuation coefficient in the presence of elevation aperture blockage |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8523776B2 (en) | Ultrasonic doppler imaging apparatus and method with blood velocity waveform processing | |
US7798968B2 (en) | Automatic detection system and method of spectral Doppler blood flow velocity | |
CN101138507B (en) | Doppler bloodstream aural signal processing method and device thereof | |
US20070016045A1 (en) | Method and apparatus for extracting an envelope curve of a spectrogram | |
CN104997534B (en) | Ultrasonic attenuation coefficient imaging method based on belt restraining least square method | |
CN101766497B (en) | Method for processing signal of sound spectrogram image and system therefor | |
US4574635A (en) | Monitoring of frequency shift of ultrasound pulses in tissue | |
CN112515704B (en) | Blood vessel hardness measuring method based on ultrasound | |
CN107106125B (en) | System and method for measuring arterial parameters | |
Kiserud et al. | Validation of diameter measurements by ultrasound: intraobserver and interobserver variations assessed in vitro and in fetal sheep | |
Ikeshita et al. | Noninvasive measurement of transient change in viscoelasticity due to flow-mediated dilation using automated detection of arterial wall boundaries | |
Parker et al. | Comparison of techniques for in vivo attenuation measurements | |
CN117643480B (en) | Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom | |
JP2003126090A (en) | In vivo signal measurement device and ultrasonic diagnostic system | |
US20220280138A1 (en) | Ultrasonic imaging of acoustic attenuation coefficients with elevation compounding | |
CN117643480A (en) | Method for estimating ultrasonic attenuation coefficient by using aberration compensation and reference phantom | |
CN115844452B (en) | Pulse wave detection method, device and storage medium | |
Reuwer et al. | Characteristics of flow velocities in the umbilical arteries, assessed by Doppler ultrasound | |
CN111466951B (en) | Method and device for generating ultrasonic attenuation image, ultrasonic equipment and storage medium | |
Ibrahim et al. | Detection of boundaries of carotid arterial wall by analyzing ultrasonic RF signals | |
Kathpalia et al. | A robust Doppler spectral envelope detection technique for automated blood flow measurements | |
Strowitzki et al. | Ultrasonic radio-frequency spectrum analysis of normal brain tissue | |
Gao et al. | Design and Development of a Bladder Volume Determination Device Based on A-mode Ultrasound | |
JP2010068955A (en) | Ultrasonic diagnostic apparatus and data processing program of ultrasonic diagnostic apparatus | |
WO2023274763A1 (en) | Ultrasonic quantification of acoustic attenuation coefficient in the presence of elevation aperture blockage |
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 | ||
GR01 | Patent grant |