CN107669334B - Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy - Google Patents
Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy Download PDFInfo
- Publication number
- CN107669334B CN107669334B CN201711049911.3A CN201711049911A CN107669334B CN 107669334 B CN107669334 B CN 107669334B CN 201711049911 A CN201711049911 A CN 201711049911A CN 107669334 B CN107669334 B CN 107669334B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- energy
- temperature
- backscatter
- envelope
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 15
- 238000007674 radiofrequency ablation Methods 0.000 title claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 239000000523 sample Substances 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 19
- 238000002679 ablation Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 5
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B18/04—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
- A61B18/12—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating by passing a current through the tissue to be heated, e.g. high-frequency current
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B2018/00636—Sensing and controlling the application of energy
- A61B2018/00773—Sensed parameters
- A61B2018/00791—Temperature
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Life Sciences & Earth Sciences (AREA)
- Otolaryngology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Plasma & Fusion (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention relates to an ultrasonic radio frequency ablation temperature imaging method based on ultrasonic back scattering energy, which comprises the following steps: obtaining a detection ultrasonic backscattering signal, namely initial data, by using an ultrasonic probe; performing band-pass filtering on the initial data to reduce noise; obtaining an envelope of the initial data; through sliding windows, when calculating in each window, obtaining the square of the envelope in the step 3) to obtain the energy at the moment, dividing the energy by the square of the envelope at the reference temperature, only keeping positive backscatter energy at the moment, and removing negative backscatter energy to obtain a matrix of the positive ultrasonic backscatter energy; interpolating to obtain the size of the original image; and combining the images at different moments to remove part of noise.
Description
Technical Field
The invention belongs to the field of medical image research, and relates to an ultrasonic radiofrequency ablation temperature imaging device.
Background
In today's society, tumors are one of the important causes for taking human life. Since the cause of the disease cannot be clarified yet and the prevention is carried out, the treatment of tumor is a major concern in the medical field. The traditional treatment mode mainly comprises surgical excision and has great physical trauma to patients. For tumors with small volume, ablation can be performed by using a high-energy focusing ultrasonic burning method, a radio frequency ablation method and the like, and the method has the greatest advantage of being minimally invasive, so that the audience range can be wide. The temperature condition in the burning process can be observed through nuclear magnetic resonance, ultrasonic waves, a thermocouple and the like, but the ultrasonic waves have the characteristics of harmlessness to a human body, low cost, simplicity and convenience in operation and the like, so that the method has great research value.
The existing methods for temperature monitoring by using ultrasonic waves include echo time shifting (echoing), backscatter energy (CBE) and the like. The echo time shift method estimates the temperature by means of tissue expansion and sound velocity change caused by the temperature rise of tissues in a burning area, but parameters such as a tissue expansion coefficient and the like need to be measured in advance, and the temperature rise is up to 40 ℃. The temperature of the back scattering energy is monitored by the change of the back scattering energy in the heating process, and the algorithm is simple and convenient and is suitable for real-time monitoring.
Straube and Arthur at washington university found that scatterers were divided into positive and negative, the backscattering energy value of positive scatterers increased with increasing temperature and showed positive values during ablation, and the backscattering energy value of negative scatterers decreased with increasing temperature and showed negative values during ablation. On the basis of sign and Reid, a normalized model [1-5] of the relation between the single scatterer backscatter energy change and the temperature change is calculated, and the normalized model is expressed as a formula (1-1):
wherein α (T) is the attenuation coefficient, η (T) is the backscattering coefficient, if the size of the scatterer is smaller than the wavelength of the ultrasonic wave, the model can be simplified to the formula (1-2):
CBE＝η(T)/η(37) (1-2)
tsui P H et al found that the back scattering energy and temperature relationship can be obtained by canceling the displacement compensation of the original data, and that the sensitivity is higher and the time is saved [6 ].
Based on the above theory, Xiashi et al propose a method of integrating ultrasonic backscattering energy (ICBE) and a sliding window to obtain a temperature distribution image of a heating region, and the integrated backscattering energy inverts negative backscattering energy in an ablation process and displays the temperature distribution together with positive backscattering energy. Integrated backscatter energy values within each window of the sliding window:
when the comprehensive ultrasonic backscattering energy method is used for real tissue radio frequency ablation, an artifact exists below a burning needle, the correlation between the comprehensive ultrasonic backscattering energy value and a temperature value is influenced, and a tissue temperature distribution image is also seriously influenced.
Reference documents:
[1]Straube W L and Arthur R M,Theoretical Estimation of theTemperature Dependence of Backscattered Ultrasonic Power for NoninvasiveThermometry,Ultrasound in Medicine&Biology,1994,20(9):915～922.
[2]Arthur R M,Trobaugh J W,Straube W L,et al.,Temperature Dependenceof Ultrasonic Backscattered Energy in Images Compensated for Tissue Motion,2003IEEE Symposium on Ultrasonics,2003,1:990～993.
[3]Arthur R M,Trobaugh J W,Straube W L,et al.,Temperature Dependenceof Ultrasonic Backscattered Energy in Motion Compensated Images,Ultrasonics,Ferroelectrics,and Frequency Control,2005,52(10):1644～1652.
[4]Arthur R M,Straube W L,Starman J D,et al.,Noninvasive TemperatureEstimation Based on the Energy of Backscattered Ultrasound,Medical physics,2003,30(6):1021～1029.
[5]Arthur R M,Straube W L,Trobaugh J W,et al.,Non-Invasive Estimationof Hyperthermia Temperatures with Ultrasound,International journal ofhyperthermia,2005,21(6):589～600.
[6]Po-Hsiang Tsui,Yu-Ting Chien,et al,Using Ultrasound CBE ImagingWithout Echo Shift Compensation for Temperature Estimation,Ultrasonics,2012,52:925–935.
disclosure of Invention
The invention aims to provide an ultrasonic radio frequency ablation temperature imaging device which can effectively remove the artifact and ensure that a temperature distribution image is more accurate and can better accord with a temperature value change curve, and the technical scheme is as follows:
1. an ultrasonic radio frequency ablation temperature imaging device based on ultrasonic back scattering energy comprises an ultrasonic probe, a band-pass filtering module, a temperature imaging module and a display module, wherein,
an ultrasonic probe for detecting an ultrasonic backscatter signal, i.e., initial data;
the band-pass filtering module is used for performing band-pass filtering on the initial data to reduce noise;
the temperature imaging module is used for completing the following functions:
(1) obtaining an envelope of the initial data;
(2) through sliding windows, when calculating in each window, obtaining the energy at the moment by the square of the envelope obtained in the step (1), dividing the energy by the square of the envelope at the reference temperature, only keeping positive backscatter energy at the moment, and removing negative backscatter energy to obtain a matrix of the positive ultrasonic backscatter energy;
(3) interpolating to obtain the size of the original image;
and the display module is used for combining the original images subjected to interpolation at different moments by the temperature imaging module to remove part of noise and displaying the combined original images.
The invention has the following beneficial effects:
1. when the comprehensive ultrasonic backscatter energy is used for radiofrequency ablation, severe artifacts can be produced below the burning needle. By applying the positive ultrasonic wave back scattering energy, the artifact can be effectively removed, so that the temperature distribution image is more accurate, and the method can be better applied to clinic.
2. The fitting effect of the positive ultrasonic backscattering energy value and the temperature is better than that of the comprehensive ultrasonic backscattering energy, so that the change of the temperature value in the ablation process can be better known by using the positive backscattering energy.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a graph comparing a forward ultrasonic backscatter energy method, a comprehensive ultrasonic backscatter energy method, with a temperature curve, respectively.
Detailed Description
The ultrasonic radio frequency ablation temperature imaging method based on the ultrasonic back scattering energy improves the artifact problem of a temperature distribution image displayed in real time in the radio frequency ablation process by using a Positive ultrasonic back scattering energy (PCBE) method, and the value of the Positive ultrasonic back scattering energy is more in line with the curve of a temperature value.
During the temperature monitoring process, only the positive backscatter energy part is reserved, and the negative backscatter energy part is removed, namely, the positive backscatter energy is only utilized to display a temperature distribution image and calculate the backscatter energy value, so that an artifact below a burning needle can be effectively improved during the real-time temperature monitoring of radio frequency ablation, and the positive ultrasonic backscatter energy value is more in line with a temperature curve than the comprehensive ultrasonic backscatter energy value.
The method comprises the following specific steps:
1) obtaining a detection ultrasonic backscattering signal, namely initial data, by using an ultrasonic probe;
2) band-pass filtering of initial data to reduce noise
3) Determining an envelope of initial data
4) Through sliding windows, when calculating in each window, obtaining the square of the envelope in the step 3) to obtain the energy at the moment, dividing the energy by the square of the envelope at the reference temperature, only keeping positive backscattering energy at the moment, removing negative backscattering energy, and obtaining a matrix of the positive ultrasonic backscattering energy
5) Interpolated imaging original image size
6) Combining the images at different moments to remove part of noise and display
7) And (4) selecting the interested area, and averaging the part of the matrix obtained in the step 4) in the interested area to obtain the forward and backward scattering energy value in the interested area. Subsequent calculation of ablation area and the like can also be performed
The examples are as follows:
the pig internal spine with proper size is placed in an acrylic box, a radio frequency burning needle is inserted into the internal spine through a small hole in the box, the burning needle is provided with water circulation, and the temperature is monitored by a thermocouple. And (4) finding a section of the needle tip by using an ultrasonic probe, and burning in a 50W mode. One initial data was recorded every two seconds with a 12min duration of heating.
And processing the obtained data: filtering each piece of data based on the center frequency of the ultrasonic probe for experiment to remove noise outside the bandwidth; performing Hilbert transform on the data without the noise to obtain signal envelope; dividing the image into a plurality of windows with the same size to obtain a positive ultrasonic backscattering energy value in each window; interpolating the obtained matrix to obtain the size of the original image; compounding data at different moments to eliminate partial noise; and displaying the temperature distribution image of the forward ultrasonic wave back scattering energy, and calculating the average value of the forward ultrasonic wave back scattering energy in the region of interest on the basis of the temperature distribution image and comparing the average value with the temperature value recorded by the thermocouple.
The obtained results can be compared with an integrated ultrasonic backscatter energy (ICBE) method, and the advantages of the positive ultrasonic backscatter energy (PCBE) in radio frequency ablation temperature monitoring can be verified.
Claims (1)
1. An ultrasonic radio frequency ablation temperature imaging device based on ultrasonic back scattering energy comprises an ultrasonic probe, a band-pass filtering module, a temperature imaging module and a display module, wherein,
an ultrasonic probe for detecting an ultrasonic backscatter signal, i.e., initial data;
the band-pass filtering module is used for performing band-pass filtering on the initial data to reduce noise;
the temperature imaging module is used for completing the following functions:
(1) obtaining an envelope of the initial data;
(2) through sliding windows, when calculating in each window, obtaining the energy of the window at the moment by the square of the envelope obtained in the step (1), dividing the energy by the square of the envelope at the reference temperature, only keeping positive backscatter energy at the moment, and removing negative backscatter energy to obtain a matrix of the positive ultrasonic backscatter energy;
(3) interpolating to obtain the size of the original image;
and the display module is used for combining the original images subjected to interpolation at different moments by the temperature imaging module to remove part of noise and displaying the combined original images.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711049911.3A CN107669334B (en) | 2017-10-31 | 2017-10-31 | Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711049911.3A CN107669334B (en) | 2017-10-31 | 2017-10-31 | Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107669334A CN107669334A (en) | 2018-02-09 |
| CN107669334B true CN107669334B (en) | 2020-06-16 |
Family
ID=61142508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711049911.3A Expired - Fee Related CN107669334B (en) | 2017-10-31 | 2017-10-31 | Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107669334B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109394263B (en) * | 2018-09-25 | 2021-06-18 | 北京工业大学 | Ultrasonic scatterer diameter multi-scale imaging method based on backscattering coefficient |
| CN110681077A (en) * | 2019-09-25 | 2020-01-14 | 天津大学 | HIFU hotspot imaging method based on phantom and frequency domain single-phase CBE |
| CN112966212B (en) * | 2021-02-07 | 2022-08-16 | 西安交通大学 | Multi-parameter real-time monitoring imaging system based on ultrasonic echo backscattering energy change |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103512960A (en) * | 2013-09-27 | 2014-01-15 | 中国科学院声学研究所 | Ultrasound array imaging method |
| CN104155362A (en) * | 2014-08-14 | 2014-11-19 | 东南大学 | Ultrasound imaging technology based method for detecting flow pattern of gas-liquid two-phase flow in micro channel |
| CN104330475A (en) * | 2014-10-23 | 2015-02-04 | 陕西师范大学 | Metal anti-fake identification method based on ultrasonic backscattering attenuation coefficient spectrum |
| CN106264722A (en) * | 2016-08-31 | 2017-01-04 | 天津大学 | Window modulation combines Gauss polynomial matching monitoring radio frequency ablation device and method |
| CN106373103A (en) * | 2016-09-08 | 2017-02-01 | 飞依诺科技(苏州)有限公司 | Ultrasonic data compounding method and apparatus |
| CN107212903A (en) * | 2016-03-22 | 2017-09-29 | 美国西门子医疗解决公司 | Relative backscattering coefficient in medical diagnostic ultrasound |
| WO2017176101A1 (en) * | 2016-04-05 | 2017-10-12 | Частное Учреждение "Назарбаев Университет Рисеч Энд Инновейшн Систем" | Method for distributed temperature sensing during thermal ablation of a tumour using a fibre optic temperature sensor with a linearly chirped bragg grating |
-
2017
- 2017-10-31 CN CN201711049911.3A patent/CN107669334B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103512960A (en) * | 2013-09-27 | 2014-01-15 | 中国科学院声学研究所 | Ultrasound array imaging method |
| CN104155362A (en) * | 2014-08-14 | 2014-11-19 | 东南大学 | Ultrasound imaging technology based method for detecting flow pattern of gas-liquid two-phase flow in micro channel |
| CN104330475A (en) * | 2014-10-23 | 2015-02-04 | 陕西师范大学 | Metal anti-fake identification method based on ultrasonic backscattering attenuation coefficient spectrum |
| CN107212903A (en) * | 2016-03-22 | 2017-09-29 | 美国西门子医疗解决公司 | Relative backscattering coefficient in medical diagnostic ultrasound |
| WO2017176101A1 (en) * | 2016-04-05 | 2017-10-12 | Частное Учреждение "Назарбаев Университет Рисеч Энд Инновейшн Систем" | Method for distributed temperature sensing during thermal ablation of a tumour using a fibre optic temperature sensor with a linearly chirped bragg grating |
| CN106264722A (en) * | 2016-08-31 | 2017-01-04 | 天津大学 | Window modulation combines Gauss polynomial matching monitoring radio frequency ablation device and method |
| CN106373103A (en) * | 2016-09-08 | 2017-02-01 | 飞依诺科技(苏州)有限公司 | Ultrasonic data compounding method and apparatus |
Non-Patent Citations (2)
| Title |
|---|
| 射频消融治疗中超声温度影像与弹性成像的可靠性研究;李锵;《万方学位论文》;20140731;全文 * |
| 超声波温度成像与弹性成像技术研究;夏静静;《万方学位论文》;20170801;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107669334A (en) | 2018-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Miller et al. | Fundamental limitations of noninvasive temperature imaging by means of ultrasound echo strain estimation | |
| Fahey et al. | Acoustic radiation force impulse imaging of the abdomen: demonstration of feasibility and utility | |
| Subramanian et al. | In vivo thermal ablation monitoring using ultrasound echo decorrelation imaging | |
| Miller et al. | Imaging of temperature-induced echo strain: preliminary in vitro study to assess feasibility for guiding focused ultrasound surgery | |
| EP1803403B1 (en) | Ultrasound diagnostic system of detecting a lesion | |
| US9643030B2 (en) | High-intensity focused ultrasound thermal ablation apparatus having integrated temperature estimation and elastography for thermal lesion determination and the method thereof | |
| CN107669334B (en) | Ultrasonic radio frequency ablation temperature imaging device based on ultrasonic backscattering energy | |
| Zhou et al. | Ultrasonic evaluation of microwave-induced thermal lesions based on wavelet analysis of mean scatterer spacing | |
| Dahl et al. | Coherence beamforming and its applications to the difficult-to-image patient | |
| Civale et al. | Calibration of ultrasound backscatter temperature imaging for high-intensity focused ultrasound treatment planning | |
| Zhang et al. | A fast tissue stiffness-dependent elastography for HIFU-induced lesions inspection | |
| Gudur et al. | High-frequency rapid B-mode ultrasound imaging for real-time monitoring of lesion formation and gas body activity during high-intensity focused ultrasound ablation | |
| Yang et al. | Post-procedure evaluation of microwave ablations of hepatocellular carcinomas using electrode displacement elastography | |
| Zhang et al. | Ultrasound single-phase CBE imaging for monitoring radiofrequency ablation | |
| CN106264607B (en) | Real-time ultrasound wave temperature imaging method and apparatus based on time signal offset | |
| Li et al. | Ultrasonic Nakagami visualization of HIFU-induced thermal lesions | |
| Anand et al. | Monitoring formation of high intensity focused ultrasound (HIFU) induced lesions using backscattered ultrasound | |
| Liu et al. | Instantaneous frequency-based ultrasonic temperature estimation during focused ultrasound thermal therapy | |
| Yang et al. | Delineation of post-procedure ablation regions with electrode displacement elastography with a comparison to acoustic radiation force impulse imaging | |
| Qiu et al. | A scanning-mode 2D shear wave imaging (s2D-SWI) system for ultrasound elastography | |
| Alvarenga et al. | Influence of temperature variations on the entropy and correlation of the Grey-Level Co-occurrence Matrix from B-Mode images | |
| Athamnah et al. | Diagnostice breast elastography estimation from doppler imaging using central difference (CD) and least-squares (LS) algorithms | |
| Shi et al. | Phase aberration in shear wave dispersion ultrasound vibrometry | |
| Yang et al. | Ultrasound monitoring of temperature and coagulation change during tumor treatment with microwave ablation | |
| Gargir et al. | Assessment of coded excitation implementation for estimating heat-induced speed of sound changes |
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 | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200616 |