CN110514742A - 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法 - Google Patents

基于仿体和频域cbe的聚焦式超声波消融的温度成像方法 Download PDF

Info

Publication number
CN110514742A
CN110514742A CN201910817973.7A CN201910817973A CN110514742A CN 110514742 A CN110514742 A CN 110514742A CN 201910817973 A CN201910817973 A CN 201910817973A CN 110514742 A CN110514742 A CN 110514742A
Authority
CN
China
Prior art keywords
cbe
frequency domain
imitative body
ultrasonic
window
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910817973.7A
Other languages
English (en)
Inventor
杨昆
李锵
关欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201910817973.7A priority Critical patent/CN110514742A/zh
Publication of CN110514742A publication Critical patent/CN110514742A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654—Imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
    • G01N31/12—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using combustion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration by non-spatial domain filtering
    • G06T5/70—
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids

Abstract

本发明涉及一种基于仿体和频域CBE的聚焦式超声波消融的温度成像方法，包括下列步骤：制作仿体；对仿体进行HIFU灼烧实验；超声影像系统置于仿体上方以获取超声逆散射信号，即RF data；通过滑动窗，将影像切割为数个小窗；利用傅里叶变换将窗内的空间域信号转换至频率域；截取用于CBE成像频率范围即3±0.75MHz内的信号，并找出范围内的最大振幅；将最大振幅取平方，则获得频域范围内最大逆散射能量；根据各个窗的最大逆散射能量的变化值，得到频域CBE值，进行成像。

Description

基于仿体和频域CBE的聚焦式超声波消融的温度成像方法
技术领域
本发明属于医学图像研究领域，涉及一种高能聚焦式超声波热疗过程中的超声温度成像方法。
背景技术
近些年来，癌症愈演愈烈，成为夺走人类生命健康的一大杀手。根据美国癌症协会(American Cancer Society,ACS)的纪录，现今癌症治疗的主要方式分别为：手术切除、化学治疗、放射治疗、标靶治疗、免疫治疗。除上述五种治疗方法外，肿瘤热治疗在近二十年也逐渐受到医学界的注目。
高强度聚焦超声(HIFU)是一种广泛使用的无创热疗工具。为控制HIFU的治疗质量，需要对组织中的灼烧点进行定位[1-2]。超声成像能够监测加热区域；特别地，逆散射能量变化(CBE)参数成像可以使组织中的热信息可视化。逆散射能量变化法为近代测量温度影像的主流，具有图像处理快速、可直接转换能量信号等优势。
华盛顿大学Arthur学者提出超声波逆散射能量变化温度影像主要和逆散射子特性有关[3]，而在后续研究中发现，组织内散射子的CBE值与温度变化几乎呈线性关系[4]。因此，可通过取得不同温度下CBE值以得出组织温度，并可进一步转换为温度影像。相较于其他超声波影像成像方式，CBE演算法于37℃到50℃之间几乎呈单调函数的关系[5]。推算出单一散射子背向散射能量变化与温度变化关系的归一化模型，如(1-1)式：
其中，α(T)为衰减系数，η(T)为背向散射系数。若散射子的尺寸小于超声波波长，则模型可简化为式(1-2)
CBE＝η(T)/η(37) (1-2)
CBE影像是在像素点到像素点的基础上进行的成像，对散斑运动和超声背散射信号特征的变化是敏感的。集成逆散射能量变化(ICBE)影像算法[6]在传统CBE算法的基础上运用滑窗技术，这种窗口到窗口的方法可以实现温度分布的可视化，且无需执行额外的运动补偿或相位校正处理。
然而，在HIFU热点定位过程中，温度上升较小，导致热点处的CBE值较小。因此，在空间域中构建的常规CBE成像虽可以呈现HIFU热区，但容易受到噪声的影响。在这种情况下，空间域CBE图像亮暗区对比度较低，导致HIFU加热点可视化成程度低，严重影响到CBE影像对HIFU灼烧点的术中监控。
参考文献：
[1]Wu F,Wang Z B,Chen W Z,et al.Extracorporeal high intensity focusedultrasound ablation in the treatment of patients with large hepatocellularcarcinoma[J].Annals of surgical oncology,2004,11(12):1061.
[2]Illing R O,Kennedy J E,Wu F,et al.The safety and feasibility ofextracorporeal high-intensity focused ultrasound(HIFU)for the treatment ofliver and kidney tumours in a Western population[J].British journal ofcancer,2005,93(8):890.
[3]Arthur R M,Broadstone S R.Imaging via inversion of ellipsoidalprojections of solutions to the linear acoustic wave equation(medicalultrasound)[J].IEEE transactions on medical imaging,1989,8(1):89-95.
[4]Arthur R M,Straube W L,Trobaugh J W,et al.Non-invasive estimationof hyperthermia temperatures with ultrasound[J].International journal ofhyperthermia,2005,21(6):589-600.
[5]W L S.Theoretical estimation of the temperature dependence ofbackscattered ultrasonic power for noninvasive thermometry.Ultrasound inmedicine&biology1994；9.
[6]An approach for the visualization of temperature distribution intissues according to changes in ultrasonic backscattered energy
发明内容
本发明的目的是提供一种基于HIFU灼烧仿体模拟实验的超声波温度成像方法，为克服空间域CBE影像算法在HIFU监控过程中热区可视化水平低的问题提供仿真基础，技术方案如下:
一种基于仿体和频域CBE的聚焦式超声波消融的温度成像方法，包括下列步骤：
1)制作仿体；
2)对仿体进行HIFU灼烧实验；
3)超声影像系统置于仿体上方以获取超声逆散射信号，即RF data；
4)通过滑动窗，将影像切割为数个小窗。
5)利用傅里叶变换将窗内的空间域信号转换至频率域；
6)截取用于CBE成像频率范围即3±0.75MHz内的信号，并找出范围内的最大振幅；
7)将最大振幅取平方，则获得频域范围内最大逆散射能量；
8)根据各个窗的最大逆散射能量的变化值，得到频域CBE值，进行成像。
实验成果表明，与空间域CBE成像相比，频域CBE成像将影像亮暗区对比度由8.3提升至12.06，可以提高HIFU热区的可视化水平。为进一步实现HIFU灼烧热点的监控奠定了基础。
附图说明
附图1为本发明实验使用的石墨仿体图。
附图2为本发明使用的实验装置图。
附图3为本发明的算法流程图。
附图4为空间域及频域超声背散射能量变化影像对比度的比较图。
具体实施方法
本发明基于空间域集成背散射能量变化方法发明而来，利用频域信号进行CBE成像，有效的抑制了CBE影像背景处的噪声，提升了亮暗区的对比度，改善了热区的可视化程度。
在进行频域CBE测算过程中，将背散射信号经过傅里叶变换转换至频率域，通过截取特定频域范围内的信号可有效的消除噪声的影响。且计算出的频域CBE值也与温度变化值有较高的相关性。
具体步骤如下：
1)利用Terason3000超声影像系统获取超声逆散射信号，即RF data。
2)通过滑动窗，将影像切割为数个小窗。
3)利用傅里叶变换将窗内的空间域信号转换至频率域。
4)截取用于CBE成像频率范围内的信号，并找出该范围内的最大振幅。
5)将最大振幅取平方，则可获得频域范围内最大逆散射能量。
6)将各个窗的最大逆散射能量使用常规逆散射能量变化的计算方式，可以得到频域CBE值进行成像。
实施例如下：
本发明通过HIFU灼烧体外仿体实验进行验证。仿体实体图如图1所示。将2g琼脂粉溶解在200ml水中而制备成琼脂液，再将10g平均直径小于20μm的石墨粉加入到琼脂液中。煮沸冷却后制成散射子浓度为128scatterers/mm3的仿体，大小约为5×6×7cm3。
高能聚焦式超声波灼烧仿体的实验装置如图2所示，此装置由三个系统组成：HIFU系统，超声波监测系统和温度系统。在HIFU系统中，连续波首先由信号发生器产生，然后由功率放大器放大。最后，将放大的信号发送到HIFU探头。HIFU探头的频率为2.12MHz，焦距为5.5厘米。放大的信号也将被发送到功率计，以监控HIFU的输出功率。超声波监测系统是超声波扫描仪，使用弧形探头进行扫描。探头的中心频率为3MHz，脉冲长度为0.7mm。若要使用平面波复合成像则需更换相应的成像系统。在温度系统中，将热电偶温度计插入仿体中以收取温度数据。为了保证每次实验的起始温度相同，恒温器须将水箱温度控制在33℃。
HIFU分别以5W，10W，15W的功率对仿体进行灼烧。因为需要比较超声图像和模拟温度之间的关系，所以需要获得超声信号和温度信息。超声信息通过Terason3000系统获取，温度信息有热电偶温度计插入仿体获取。然而，当热电偶刺入仿体时，由于热电偶由金属制成，因此超声灰度图像上将出现阴影伪影。当RF数据用于计算时，所获得的图像和信号强度可能具有偏差。为了解决不能同时获得温度信息和射频信号的问题，本实验将进行两次,分别获得超声信号和温度信息。
在第一次实验中，在没有热电偶的情况下获得超声RF数据。第二次插入热电偶。为避免伪温升的发生，热电偶应略微偏移3mm，以便可以离开HIFU的聚焦区域，然后记录仿体的温度变化。在每次实验之前，将仿体冷却至33℃的背景温度。
对获取的超声数据信息进行频域CBE的测算，并计算影像热区与背景区的对比度。之后需将超声信息与温度信息进行相关性的测算，以确保所提出的频域CBE算法对温度监测的准确性。将取得的结果与空间域CBE影像进行比较，验证频域CBE提升HIFU灼烧区可视化水平的潜力。

Claims (1)

1.一种基于仿体和频域CBE的聚焦式超声波消融的温度成像方法，包括下列步骤：
1)制作仿体；
2)对仿体进行HIFU灼烧实验；
3)超声影像系统置于仿体上方以获取超声逆散射信号，即RF data；
4)通过滑动窗，将影像切割为数个小窗。
5)利用傅里叶变换将窗内的空间域信号转换至频率域；
6)截取用于CBE成像频率范围即3±0.75MHz内的信号，并找出范围内的最大振幅；
7)将最大振幅取平方，则获得频域范围内最大逆散射能量；
8)根据各个窗的最大逆散射能量的变化值，得到频域CBE值，进行成像。
CN201910817973.7A 2019-08-30 2019-08-30 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法 Pending CN110514742A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910817973.7A CN110514742A (zh) 2019-08-30 2019-08-30 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910817973.7A CN110514742A (zh) 2019-08-30 2019-08-30 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法

Publications (1)

Publication Number Publication Date
CN110514742A true CN110514742A (zh) 2019-11-29

Family

ID=68628669

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910817973.7A Pending CN110514742A (zh) 2019-08-30 2019-08-30 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法

Country Status (1)

Country Link
CN (1) CN110514742A (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107569256A (zh) * 2017-09-25 2018-01-12 南京广慈医疗科技有限公司 基于热膨胀和门控算法测量生物组织温度变化的超声方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107569256A (zh) * 2017-09-25 2018-01-12 南京广慈医疗科技有限公司 基于热膨胀和门控算法测量生物组织温度变化的超声方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XIA JINGJING等: "An Approach for the Visualization of Temperature Distribution in Tissues According to Changes in Ultrasonic Backscattered Energy", 《COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE》 *
黄振维: "超音波热影像及声学干扰特征用于海扶刀热点定位：可行性研究", 《台湾博硕士论文知识加值系统》 *

Similar Documents

Publication Publication Date Title
Miller et al. Fundamental limitations of noninvasive temperature imaging by means of ultrasound echo strain estimation
Miller et al. Imaging of temperature-induced echo strain: preliminary in vitro study to assess feasibility for guiding focused ultrasound surgery
Lewis et al. Thermometry and ablation monitoring with ultrasound
Varghese et al. Elastographic imaging of thermal lesions in the liver in vivo following radiofrequency ablation: preliminary results
US8016757B2 (en) Non-invasive temperature estimation technique for HIFU therapy monitoring using backscattered ultrasound
Anand et al. Three-dimensional spatial and temporal temperature imaging in gel phantoms using backscattered ultrasound
Arthur et al. Temperature dependence of ultrasonic backscattered energy in motion compensated images
Wang et al. Monitoring radiofrequency ablation with ultrasound Nakagami imaging
JP7417740B2 (ja) 超音波手技における適応的単一気泡ベースの自動焦点および出力調節
Civale et al. Calibration of ultrasound backscatter temperature imaging for high-intensity focused ultrasound treatment planning
Azhari Feasibility study of ultrasonic computed tomography–guided high-intensity focused ultrasound
Zhong et al. Monitoring imaging of lesions induced by high intensity focused ultrasound based on differential ultrasonic attenuation and integrated backscatter estimation
Barrere et al. Imaging of thermal effects during high-intensity ultrasound treatment in liver by passive elastography: a preliminary feasibility in vitro study
CN106264607B (zh) 基于时间信号偏移的实时超声波温度成像方法及设备
Karaböce et al. Investigation of different TMMs in high intensity focused ultrasound applications
Liu et al. Noninvasive estimation of temperature elevations in biological tissues using acoustic nonlinearity parameter imaging
Grutman et al. Dense speed-of-sound shift imaging for ultrasonic thermometry
Miller et al. Ultrasonic temperature imaging for guiding focused ultrasound surgery: Effect of angle between imaging beam and therapy beam
Hou et al. High-intensity focused ultrasound monitoring using harmonic motion imaging for focused ultrasound (HMIFU) under boiling or slow denaturation conditions
Yeshurun et al. Non-invasive measurement of thermal diffusivity using high-intensity focused ultrasound and through-transmission ultrasonic imaging
CN110514742A (zh) 基于仿体和频域cbe的聚焦式超声波消融的温度成像方法
Chen et al. Real‐time tissue change monitoring on the Sonablate® 500 during high intensity focused ultrasound (HIFU) treatment of prostate cancer
Takeuchi et al. Statistical analysis of ultrasonic scattered echoes enables the non-invasive measurement of temperature elevations inside tumor tissue during oncological hyperthermia
Liu et al. A unified approach to combine temperature estimation and elastography for thermal lesion determination in focused ultrasound thermal therapy
Huang et al. Ultrasound thermal mapping based on a hybrid method combining cross-correlation and zero-crossing tracking

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191129

RJ01 Rejection of invention patent application after publication