CN106037815A - 一种用于热凝固监测的超声回波统计参量成像***与方法 - Google Patents

一种用于热凝固监测的超声回波统计参量成像***与方法 Download PDF

Info

Publication number
CN106037815A
CN106037815A CN201610326626.0A CN201610326626A CN106037815A CN 106037815 A CN106037815 A CN 106037815A CN 201610326626 A CN201610326626 A CN 201610326626A CN 106037815 A CN106037815 A CN 106037815A
Authority
CN
China
Prior art keywords
ultrasonic echo
statistical parameter
ultrasonic
imaging
thermocoagulation
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
CN201610326626.0A
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610326626.0A priority Critical patent/CN106037815A/zh
Publication of CN106037815A publication Critical patent/CN106037815A/zh
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种用于热凝固监测的超声回波统计参量成像***与方法，其成像***包括主控设备、功率放大器、无创功率换能器、超声检测与成像换能器、全数字化超声设备、数据采集设备、样品、图像采集设备、指示灯、微创功率换能器、透光水槽和三维移动设备；本发明采用克服微泡覆盖和上部强回声区域遮挡的超声回波统计参量成像方法监测热凝固的形成过程，实现了热凝固起始的提早检测，增强了热凝固下部的显示，提高热消融过程中热凝固微泡比增强了热凝固监测对比度。

Description

一种用于热凝固监测的超声回波统计参量成像***与方法
【技术领域】
本发明属于超声检测和监控成像技术，涉及生物医学和信息领域，具体涉及一种用于热凝固监测的超声回波统计参量成像***与方法。
【背景技术】
以温度升高引起蛋白非可逆性热凝固为机制的热消融技术具有无创或微创、快速、副作用小等优势，国内外使用射频、微波、超声、激光等“绿色”物理能量实现的热消融技术在生物、医学、化学和工业等领域取得了令人鼓舞的成果，特别在国内使用高强度聚焦超声热消融技术***方面取得了开创性的成果。热凝固过程的检测与监控成像是实现高效和安全热消融的基础和保障。超声和磁共振成像是目前常见的监控成像方式。超声由于其成本低、实时性好、适用范围广和易于与诸多热消融技术结合等优势已成为目前使用较广的检测与监控成像方式。已有国内外的专利使用传统超声B模式图像对热消融过程进行监控成像，例如美国专利US6425867B1和中国专利CN1565671A等。在热消融过程中除了引起蛋白非可逆性热凝固外，还会出现空化汽化引起的微泡。B超图像实际得到的是组织背向散射信息，热消融引起的热凝固区域的背向散射信号与热消融前变化并不大，B超图像观测到的信息大部分是空化汽化引起的气泡的散射信息。因此，B超图像用于热凝固过程的检测与监控成像具有一定的局限性。许多学者开展了其他超声成像技术在热凝固过程的检测与监控成像的研究，包括超声温度估计与成像、超声组织定征参量监控成像、弹性成像等。
由于受消融过程中空化汽化微泡及强回声区域的影响，热消融过程中凝固起始与凝固底部的检测与监控成像具有一定的困难。超声回波统计参量成像方法依赖于背向散射回波包络的统计分布，受回波幅度的影响较小。
【发明内容】
针对上述检测与监控中的不足，本发明涉及一种用于热凝固监测的超声回波统计参量成像***与方法。该发明同步监测热消融过程中的超声背向散射信号和超声统计参量；并利用超声统计参量图像和B模式超声图像对目标区域进行同步动态成像监控。
为了实现上述任务，本发明采取如下的技术方案：
一种用于热凝固监测的超声回波统计参量成像***，包括主控设备；所述的主控设备上连接有功率放大器、数据采集设备、图像采集设备和三维移动设备；其中，所述的功率放大器的输出端连接有无创功率换能器和微创功率换能器；所述的数据采集设备的输入端连接有全数字化超声设备，全数字化超声设备的输入端连接有超声检测与成像换能器；所述的三维移动设备用于调节样品在透光水槽中的空间位置；所述的图像采集设备用于采集样品热消融动态图像；
所述的主控设备控制功率放大器驱动无创功率换能器和/或微创功率换能器对样品发射能量；全数字化超声设备通过超声检测与成像换能器采集超声回波射频数据，并通过数据采集设备输入到主控设备进行存贮与处理。
所述的图像采集设备分辨率为1280×1024的高清摄像机，其设置于样品侧面。
还包括受主控设备的时序控制、在图像采集设备采集到的动态图像中表征时间的指示灯。
所述的无创功率换能器为聚焦超声热消融辐射器；所述的微创功率换能器为微波热消融辐射器、射频热消融辐射器或激光热消融辐射器；
所述的透光水槽为透明的玻璃或有机玻璃水槽；所述的三维移动设备为微米级扫描精度且可编程的全数字化三维扫描移动设备。
一种用于热凝固监测的超声回波统计参量成像***的成像方法，包括以下步骤：
1)主控设备控制三维移动设备调节样品的空间位置，将样品放置于无创功率换能器和/或微创功率换能器的有效辐射区内；
2)主控设备控制功率放大器驱动无创功率换能器和/或微创功率换能器发射能量，作用于 样品；
3)将超声检测与成像换能器对准样品；
4)主控设备控制全数字化超声设备通过超声检测与成像换能器采集超声回波射频数据；
5)超声检测与成像换能器采集到的超声回波射频数据通过数据采集设备输入到主控设备并存贮在主控设备上；
6)主控设备上的超声回波射频数据按照超声回波统计参量算法得到样品中的超声回波统计参量及其动态变化过程。
根据权利要求7所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，还包括：
7)由步骤5)中得到的超声回波射频数据和步骤6)中得到的超声回波统计参量数据同步构建动态变化的超声回波统计参量图像和B模式超声图像；
8)由步骤5)中得到的超声回波射频数据和步骤6)中得到的超声回波统计参量数据，提取热凝固区域和微泡区域的回波幅值和超声回波统计参量，按照热凝固微泡比的计算方法得到热凝固微泡比值及其动态变化过程。
作为本发明方进一步改进，步骤6)具体为基于Nakagami统计模型统计超声回波射频数据：按照预定将射频数据划分为多个目标区域，选取组织中单个区域内的超声回波射频数据按照Nakagami参量的算法计算得到该区域内的超声回波统计分析的Nakagami参量；对每个区域的射频数据单独处理，同时对多个区域的超声回波统计参量及其动态变化进行检测。
作为本发明方进一步改进，采用基于Nakagami参量的超声回波统计参量图像对目标区域进行成像：超声回波统计参量成像首先确定选取的估计方法及滑动窗尺寸，使滑动窗以单像素点为步长在成像平面被进行纵向及横向滑动，每滑动一个像素点，利用选取的估计方法计算滑动窗被回波信号包络的统计参量值，并将该值赋给滑动窗的中心像素点，遍历整个成像区域以 得到超声回波统计参量图像。作为本发明方进一步改进，构建超声回波统计参量图像时首先对超声回波射频信号进行去噪处理，基于最大似然估计Nakagami参量成像的去噪方法为：
首先，将采集到的超声回波射频数据分别加入两次随机白噪声，得到有噪声的信号S1和S2，加入的白噪声幅值等于***白噪声幅值，***白噪声幅值由无回声区超声信号估计所得；
其次，以脉冲宽度作为采样窗宽宽计算S1和S2的局部相关系数，得到相关系数矩阵；
接着，以80％作为阈值，找出相关系数矩阵中小于80％的值，将它们对应的值全部设置为回波射频信号的最小值，得到去***白噪声后的射频信号；
最后，将判断为噪声的点置为回波信号中的最小值，将y小于0.25对应的统计参量值置为0，将y大于0.2的对应的统计参量值按照最大似然估计算法进行求解。
本发明与现有技术相比，具有下列优点：
本发明的监控***由主控设备及连接在主控设备上的功率放大器、数据采集设备、图像采集设备和三维移动设备构成，结构简单，功能完善，本发明检测多区域的超声统计参量，实现了不同方式热消融过程中超声统计参量同步检测、热凝固微泡比同步检测与损伤动态变化同步检测。
与现有热凝固检测与监控成像技术相比，本发明的方法对超声回波射频数据按照超声回波统计参量算法得到超声回波统计参量，通过分析超声回波统计参量的动态变化和超声统计参量成像，可以克服热消融过程中空化汽化微泡的覆盖和上部强回声区域的遮挡，将超声统计参量成像方法用于不同方式热消融过程中凝固的检测与监控成像，实现热凝固起始的提早检测，增强凝固下部显示，优化热凝固过程的检测与监控成像效果。
进一步，与现有热凝固检测与监控成像技术相比，本发明还提高了热凝固微泡比增强了图像的对比度。
【附图说明】
图1是本发明提出的用于热凝固监测的超声回波统计参量监控***的框图。
图2是热消融过程所中超声回波统计参量及回波幅值同步检测与成像***与方法控制的流程图。
图3是改进的基于最大似然估计的Nakagami参量图像去噪及成像算法流程图。
图4是热消融过程热凝固起始超声回波统计参量图像、B模式图像同步监测的结果。
图5是热消融过程中出现上部强回声区域遮挡时超声回波统计参量图像与B模式图像同步监测及统计参量与回波幅值动态变化结果。
图6是热消融过程中超声回波统计参量及B模式同步监测与图像热凝固微泡比动态变化的结果。
其中，1、为主控设备；2、功率放大器；3、无创功率换能器；4、超声检测与成像换能器；5、全数字化超声设备；6、数据采集设备；7、样品；8、图像采集设备；9、指示灯；10、微创功率换能器；11、透光水槽；12、三维移动设备。
【具体实施方式】
为了更清楚的理解本发明，以下结合附图及实施方式作进一步的详细说明：
一种用于热凝固监测的超声回波统计参量成像***。
参见图1，本发明的一种用于热凝固监测的超声回波统计参量成像***包括：主控设备1、功率放大器2、无创功率换能器3、超声检测与成像换能器4、全数字化超声设备5、数据采集设备6、样品7、图像采集设备8、指示灯9、微创功率换能器10、透光水槽11和三维移动设备12；在主控设备1上连接有功率放大器2、数据采集设备6、图像采集设备8、指示灯9和三维移动设备12；其中，主控设备1的输出端与功率放大器2的控制端相连，功率放大器2的输出端与无创功率换能器3和微创功率换能器10的输入端相连，无创功率换能器3和微创功率换能器10设置在透光水槽11的周向；主控设备1的输出端与数据采集设备6的控制端相 连，数据采集设备6的输入端与全数字化超声设备5的输出端相连，全数字化超声设备5的输入端与超声检测与成像换能器4的输出端相连；主控设备1的输出端与三维移动设备12的控制端相连，三维移动设备12用于调节样品7在透光水槽11的空间位置。
该***通过主控设备1控制功率放大器2驱动无创功率换能器3和/或微创功率换能器10发射能量；该***通过主控设备1控制全数字化超声设备5通过超声检测与成像换能器4采集超声回波射频数据，并通过数据采集设备6输入到主控设备1并存贮在主控设备1上；图像采集设备6进行采集，指示灯9受主控设备1控制，在采集到的动态图像中表征时间。
其中，主控设备1为时间序列和能量输出双重控制与具备数据存贮及处理的多功能主控设备，其中，主控设备1能够对功率放大器2能量输出的时间和幅度同步时空控制，并对全数字化超声设备5通过超声检测与成像换能器4采集的超声回波射频数据进行存贮与处理。
所述的无创功率换能器3为聚焦超声热消融辐射器；所述的微创功率换能器10为微波热消融辐射器；所述的微创功率换能器10为射频热消融辐射器；所述的微创功率换能器10为激光热消融辐射器；
所述的透光水槽11为透明的玻璃或有机玻璃水槽；所述的三维移动设备12为微米级扫描精度且能够编程的全数字化三维扫描移动设备。
所述的图像采集设备8分辨率为1280×1024的高清摄像机，其设置于样品7侧面。
一种克服微泡覆盖反映热凝固起始阶段的超声回波统计参量检测方法，具体步骤如下：
本发明检测热消融过程中凝固起始目标区域的超声回波统计参量及其动态变化过程，见图2。主控设备1控制三维移动设备12调节样品7使其处于无创功率换能器3和/或微创功率换能器10的有效辐射区内；同时将超声检测与成像换能器4对准换能器的有效辐射区内的目标区域；***通过预设热消融参数控制功率放大器2驱动无创功率换能器3和/或微创功率换能器10发射能量，作用于样品；***控制全数字化超声设备通过超声检测与成像换能器4采集 目标区域热消融凝固起始的超声回波射频数据。
***使用全数字化超声设备5采集超声回波射频数据，按照预定将射频数据划分为多个目标区域，对每个区域的射频数据单独处理，同时对多个区域的超声回波统计参量及其动态变化进行检测。本发明中使用Nakagami统计模型统计超声背向散射信号。选取组织中单个区域内的超声回波射频数据按照Nakagami参量的算法计算得到该区域内的超声回波统计分析的Nakagami参量。用于获取Nakgami参量的估计量的方法有很多，大致可以分为基于矩估计的参数估计方法和基于最大似然估计的参数估计方法。
一种克服微泡覆盖提早检测热凝固起始的超声回波统计参量成像方法，具体步骤如下：
本发明同时使用B模式图像和超声回波统计参量图像对热消融凝固起始进行检测与监控成像。主控设备1通过预设参数控制功率放大器2驱动无创功率换能器3和/或微创功率换能器10发射能量，作用于样品7；同时控制全数字化超声设备5通过超声检测与成像换能器4采集目标区域热消融凝固起始的超声回波射频数据。使用采集到的超声回波射频数据同步构建动态变化的超声回波统计参量图像和B模式超声图像。
常规B模式图像反应的目标区域的超声背向散射信息，反应了组织内部声阻抗差的差异，是目前较为成熟的超声成像方法。热消融过程中会产生空化汽化微泡，空化汽化微泡与损伤在B模式图像中均表现为高回声，从而导致热消融凝固起始检测具有一定的困难。Nakagmi统计模型中，不同散射子分布得到的Nakagami参量不同，Nakagmi参量图像具有分辨不同散射子浓度的能力，并且不受回波幅值的影响。
为了克服热消融凝固起始空化汽化微泡遮挡的影响，采用基于Nakagami参量的超声回波统计参量图像对目标区域进行成像。超声回波统计参量成像首先确定选取的估计方法及滑动窗尺寸，使滑动窗以单像素点为步长在成像平面被进行纵向及横向滑动，每滑动一个像素点，利用选取的估计方法计算滑动窗被回波信号包络的统计参量值，并将该值赋给滑动窗的中心像素 点，遍历整个成像区域以得到超声回波统计参量图像。由于热凝固组织和微泡的散射子浓度不同，因此可以通过超声回波统计参量图像对热凝固起始检测与监控成像。
一种克服上部强回声区域遮挡反映热凝固下部的超声回波统计参量检测方法，具体步骤如下：
本发明监测热消融过程中凝固底部目标区域的超声回波统计参量及其动态变化过程。主控设备1控制三维移动设备12调节样品7使其处于无创功率换能器3和/或微创功率换能器10的有效辐射区内；同时将超声检测与成像换能器4对准换能器的有效辐射区内的目标区域；***通过预设热消融参数控制功率放大器驱动无创功率换能器3和/或微创功率换能器7发射能量，作用于样品；***控制全数字化超声设备通过超声检测与成像换能器采集目标区域热消融凝固底部被遮挡时的超声回波射频数据。
随着热消融时间的增加，由于损伤上部分强散射系数组织及微泡的影响，下部超声严重衰减，使得损伤难以观察，超声回波统计参量图像可以有效的抑制超声散射回波中的上部强回声区域遮挡效应，可以识别B模式图像不能描述的较弱信号。由于无回声区的白噪声会影响超声回波统计参量图像检测热凝固区域，需要在回波散射波形不变的情况下取出白噪声。NCA去噪方法在基于矩估计的Nakagami参量成像中取得较好效果。NCA去噪方法直接将噪声点置为0，基于矩估计的超声回波统计参量成像算法中，超声回波统计参量通过利用回波射频数据包络R的k阶矩μk计算而得，令E(·)表示期望，则，
μk＝E(Rk)
因此矩估计的超声回波统计参量成像算法中将噪声点置为零可以实现较好的去噪效果，但是基于最大似然估计的统计参量成像算法中
y = ln ( μ 2 G )
G = ( Π i = 1 N R i 2 ) 1 / N
由上式可知，按照传统的NCA去噪方法，噪声点被置为0，会导致计算其他点的m值时，只要滑动窗包含该点，就会导致G＝0，同时导致y趋向于无穷大，导致计算出的m值无意义，本发明提出了一种改进的基于最大似然估计的Nakagami图像去噪方法。参见图3，首先，将采集到的超声回波射频数据分别加入两次随机白噪声，得到有噪声的信号S1和S2，加入的白噪声幅值等于***白噪声幅值，由无回声区超声信号估计所得；其次，以脉冲宽度作为采样窗宽计算S1和S2的局部相关系数，得到相关系数矩阵；接着，以80％作为阈值，找出相关系数矩阵中小于80％的值，将它们对应的值全部赋值为回波射频信号的最小值，得到去***白噪声后的射频信号。将判断为噪声的点置为回波信号中的最小值，这样会导致计算出来的y值远远小于其正常值，因此将y小于0.2对应的统计参量值置为0，将y大于0.2的对应的统计参量值按照最大似然估计算法进行求解。可以达到较好的去噪效果。
一种克服上部强回声区域遮挡增强热凝固下部显示的超声回波统计参量成像方法，具体步骤如下：
本发明同时使用超声回波统计参量图像和B模式图像对热消融过程中出现上部强回声区域遮挡的热凝固进行检测与监控成像。主控设备通过预设参数控制功率放大器驱动无创功率换能器3和/或微创功率换能器10发射能量，作用于样品7；同时控制全数字化超声设备5通过超声检测与成像换能器4采集目标区域热消融凝固起始的超声回波射频数据。
使用采集到的超声回波射频数据同步构建动态变化的超声回波统计参量图像和B模式超声图像。构建超声回波统计参量图像的时首先对超声回波射频信号进行去噪处理。在消融结束后，沿超声检测与成像换能器4采集切面剖开样品，并对切面进行拍照获取光学图像。将光学图像、B模式图像和超声回波统计参量图像进行对比分析，研究随着热消融时间增加出现上部强回声区域遮挡时光学图像、B模式图像和超声回波统计参量图像与损伤的对应关系。
一种提高热凝固微泡比增强热凝固监测对比度的超声回波统计参量成像方法，包括以下步骤：
本发明检测热消融过程中目标区域的热凝固微泡比(lesion-to-bubble ratio，LBR)及其动态变化过程。***通过主控设备控制三维移动设备调节样品使其处于无创功率换能器3和/或微创功率换能器10的有效辐射区内；同时将超声检测与成像换能4对准换能器的有效辐射区内的目标区域；***通过预设热消融参数控制功率放大器驱动无创功率换能器和/或微创功率换能器发射能量，作用于样品7；***控制全数字化超声设备通过超声检测与成像换能器采集目标区域热消融凝固底部被遮挡时的超声回波射频数据。
LBR是指凝固与微泡散射能力的比率，定义为：
L B R = 20 log 10 I l e s i o n I b u b b l e
其中，Ilesion和Ibubble分别为凝固和微泡目标区域的平均幅值。利用采集到的回波射频数据和计算得到超声回波统计参量数据。提取感兴趣区域(region of interest，ROI)的回波幅值和超声回波统计参量，根据公式计算得到图像的LBR值及其动态变化过程。
以下是具体的实施例，但并不局限与这些实施例。
实施例1
参见图4，热消融过程中凝固起始检测与超声监控成像。在热消融治疗早期，使用超声回波统计参量成像方法对热凝固起始进行检测与监控成像。
本实施例中使用聚焦超声进行热消融治疗。使用的声功率为85W，占空比为70％，总时间为4s，在热消融过程中，使用图像采集设备采集聚焦超声热消融治疗过程中透明仿生物组织体模中损伤的动态变化过程，如果4(a)所示。同时在热消融过程中是用全数字超声设备同步采集热消融治疗过程中目标区域的超声回波背向散射数据。使用这些采集到的超声回波射频数据同步构建超声B模式图像图4(b)和超声回波统计参量图像4(c)，将热消融凝固起 始的超声回波统计参量图像、B模式图像和光学图像进行对比分析，
实施例2
参见图5，热消融过程中出现上部强回声区域遮挡热时热凝固下部超声回波统计参量动态变化与监控成像。在热消融过程中，对样品多个区域的超声回波幅值与超声回波统计参量进行同步检测。同时同步构建超声回波统计参量图像与超声B模式图像。
本实施例中使用微波治疗仪进行热消融治疗。微波热消融参数为：消融功率为40W，消融时间为5min，在热消融过程中，使用全数字化超声设备控制超声检测与成像换能器采集热消融目标区域的超声回波射频数据。同步构建B模式图像图5(a)和超声回波统计参量图像图5(b)，同时选取感兴趣区域分别监测区域内超声回波幅值5(c)与超声回波统计参量5(d)的动态变化。
实施例3
参见图6，基于提高热凝固微泡比的超声回波统计参量的热凝固监控成像。在热消融过程中，分别提取凝固区域与微泡区域的超声回波数据与超声回波统计参量数据，计算ROI区域的平均幅值，按照公式计算得到LBR值及其动态变化，同时同步构建超声回波统计参量图像和B模式超声图像。
本实施例中使用聚焦超声进行热消融治疗，使用声功率为140W，连续辐射32s，在热消融治疗过程中，使用全数字化超声设备同步采集热消融过程中目标区域超声回波射频数据。使用这些射频数据同步构建热消融过程中B模式图像图6(a)和超声回波统计参量图像图6(b)。同时利用采集到的射频数据获取B模式图像和统计参量图像的LBR的动态变化过程图6(c)，具体的方法为是分别选取凝固区域和微泡区域为ROI区域，计算ROI区域内回波幅值和超声回波统计参量值的平均值。按照LBR值的计算方法得到LBR值及动态变化过程。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的 具体实施方式仅限于此，对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单的推演或替换，都应当视为属于本发明由所提交的权利要求书确定专利保护范围。

Claims (10)

1.一种用于热凝固监测的超声回波统计参量成像***，其特征在于：包括主控设备(1)；所述的主控设备(1)上连接有功率放大器(2)、数据采集设备(6)、图像采集设备(8)和三维移动设备(12)；其中，所述的功率放大器(2)的输出端连接有无创功率换能器(3)和微创功率换能器(10)；所述的数据采集设备(6)的输入端连接有全数字化超声设备(5)，全数字化超声设备(5)的输入端连接有超声检测与成像换能器(4)；所述的三维移动设备(12)用于调节样品(7)在透光水槽(11)中的空间位置；所述的图像采集设备(8)用于采集样品(7)热凝固动态图像；
所述的主控设备(1)控制功率放大器(2)驱动无创功率换能器(3)和/或微创功率换能器(10)对样品(7)辐射能量；全数字化超声设备(5)通过超声检测与成像换能器(4)采集超声回波射频数据，并通过数据采集设备(6)输入到主控设备(1)进行存贮与处理。
2.根据权利要求1所述的一种用于热凝固监测的超声回波统计参量成像***，其特征在于：所述的图像采集设备(8)分辨率为1280×1024的高清摄像机，其设置于样品(7)侧面。
3.根据权利要求1所述的一种用于热凝固监测的超声回波统计参量成像***，其特征在于：还包括受主控设备(1)的时序控制、在图像采集设备(8)采集到的动态图像中表征时间的指示灯(9)。
4.根据权利要求1所述的一种用于热凝固监测的超声回波统计参量成像***，其特征在于：所述的无创功率换能器(3)为聚焦超声热消融辐射器；所述的微创功率换能器(10)为微波热消融辐射器、射频热消融辐射器或激光热消融辐射器。
5.根据权利要求1所述的一种用于热凝固监测的超声回波统计参量成像***，其特征在于：所述的透光水槽(11)为透明的玻璃或有机玻璃水槽；所述的三维移动设备(12)为微米级扫描精度且可编程的全数字化三维扫描移动设备。
6.基于权利要求1所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，包括以下步骤：
1)主控设备(1)控制三维移动设备(12)调节样品(7)的空间位置，将样品(7)放置于无创功率换能器(3)和/或微创功率换能器(10)的有效辐射区内；
2)主控设备(1)控制功率放大器(2)驱动无创功率换能器(3)和/或微创功率换能器(10)发射能量，作用于样品(7)；
3)将超声检测与成像换能器(4)对准样品(7)；
4)主控设备(1)控制全数字化超声设备(5)通过超声检测与成像换能器(4)采集超声回波射频数据；
5)超声检测与成像换能器(4)采集到的超声回波射频数据通过数据采集设备(6)输入到主控设备(1)并存贮在主控设备(1)上；
6)主控设备(1)上的超声回波射频数据按照超声回波统计参量算法得到样品(7)中的超声回波统计参量及其动态变化过程。
7.根据权利要求6所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，还包括：
7)由步骤5)中得到的超声回波射频数据和步骤6)中得到的超声回波统计参量数据同步构建动态变化的超声回波统计参量图像和B模式超声图像；
8)由步骤5)中得到的超声回波射频数据和步骤6)中得到的超声回波统计参量数据，提取热凝固区域与微泡区域的回波幅值和超声回波统计参量，按照热凝固微泡比的计算方法得到热凝固微泡比值及其动态变化过程。
8.根据权利要求6或7所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，步骤6)具体为基于Nakagami统计模型统计超声回波射频数据：按照预定将射频数据划分为多个目标区域，选取组织中单个区域内的超声回波射频数据按照Nakagami参量的算法计算得到该区域内的超声回波统计分析的Nakagami参量；对每个区域的射频数据单独处理，同时对多个区域的超声回波统计参量及其动态变化进行检测。
9.根据权利要求8所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，采用基于Nakagami参量的超声回波统计参量图像对目标区域进行成像：超声回波统计参量成像首先确定选取的估计方法及滑动窗尺寸，使滑动窗以单像素点为步长在成像平面被进行纵向及横向滑动，每滑动一个像素点，利用选取的估计方法计算滑动窗被回波信号包络的统计参量值，并将该值赋给滑动窗的中心像素点，遍历整个成像区域以得到超声回波统计参量图像。
10.根据权利要求9所述的一种用于热凝固监测的超声回波统计参量成像***的成像方法，其特征在于，构建超声回波统计参量图像时首先对超声回波射频信号进行去噪处理，基于最大似然估计Nakagami参量成像的去噪方法为：
首先，将采集到的超声回波射频数据分别加入两次随机白噪声，得到有噪声的信号S1和S2，加入的白噪声幅值等于***白噪声幅值，***白噪声幅值由无回声区超声信号估计所得；
其次，以脉冲宽度作为采样窗宽计算S1和S2的局部相关系数，得到相关系数矩阵；
接着，以80％作为阈值，找出相关系数矩阵中小于80％的值，将它们对应的值全部设置为回波射频信号的最小值，得到去***白噪声后的射频信号；
最后，将判断为噪声的点置为回波信号中的最小值，将y小于0.2对应的统计参量值置为0，将y大于0.2的对应的统计参量值按照最大似然估计算法进行求解。
CN201610326626.0A 2016-05-17 2016-05-17 一种用于热凝固监测的超声回波统计参量成像***与方法 Pending CN106037815A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610326626.0A CN106037815A (zh) 2016-05-17 2016-05-17 一种用于热凝固监测的超声回波统计参量成像***与方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610326626.0A CN106037815A (zh) 2016-05-17 2016-05-17 一种用于热凝固监测的超声回波统计参量成像***与方法

Publications (1)

Publication Number Publication Date
CN106037815A true CN106037815A (zh) 2016-10-26

Family

ID=57177025

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610326626.0A Pending CN106037815A (zh) 2016-05-17 2016-05-17 一种用于热凝固监测的超声回波统计参量成像***与方法

Country Status (1)

Country Link
CN (1) CN106037815A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108553763A (zh) * 2018-01-19 2018-09-21 北京工业大学 一种基于超声回波去相关成像技术的微波热疗监测方法
CN109171998A (zh) * 2018-10-22 2019-01-11 西安交通大学 基于超声深度学习的热消融区域识别监测成像方法与***
CN112566694A (zh) * 2018-05-25 2021-03-26 卡尔迪亚韦弗公司 包括用于使空化泡成像的装置的超声处理设备
CN112966212A (zh) * 2021-02-07 2021-06-15 西安交通大学 基于超声回波背向散射能量变化的多参量实时监控成像***
CN113117258A (zh) * 2019-12-30 2021-07-16 重庆融海超声医学工程研究中心有限公司 组织凝固性坏死的检测装置及方法
CN113907792A (zh) * 2021-08-31 2022-01-11 西安交通大学 非线性解相关融合超声多参量的热凝固检测与成像方法
CN114305668A (zh) * 2021-12-22 2022-04-12 聚融医疗科技(杭州)有限公司 基于解调域参量成像的超声热消融多参数监测方法及***
CN114587565A (zh) * 2022-03-01 2022-06-07 河南中医药大学 一种射频消融中温度控制方法及***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1847824A (zh) * 2006-02-27 2006-10-18 西安交通大学 基于超高速摄影技术的超声场中微泡行为分析***和方法
US20080114242A1 (en) * 2006-11-10 2008-05-15 National Taiwan University Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue
CN101623203A (zh) * 2009-08-07 2010-01-13 西安交通大学 一种瞬态物理过程中多模式和多参量同步检测成像监控***及监控方法
CN103948402A (zh) * 2014-05-13 2014-07-30 中国科学院深圳先进技术研究院 肿瘤超声成像特征提取方法和***
CN104398271A (zh) * 2014-11-14 2015-03-11 西安交通大学 血管与斑块的三维力学及组织特性成像检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1847824A (zh) * 2006-02-27 2006-10-18 西安交通大学 基于超高速摄影技术的超声场中微泡行为分析***和方法
US20080114242A1 (en) * 2006-11-10 2008-05-15 National Taiwan University Ultrasonic imaging technique for differentiating the distribution of scatterers within a tissue
CN101623203A (zh) * 2009-08-07 2010-01-13 西安交通大学 一种瞬态物理过程中多模式和多参量同步检测成像监控***及监控方法
CN103948402A (zh) * 2014-05-13 2014-07-30 中国科学院深圳先进技术研究院 肿瘤超声成像特征提取方法和***
CN104398271A (zh) * 2014-11-14 2015-03-11 西安交通大学 血管与斑块的三维力学及组织特性成像检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SIYUAN ZHANG 等: "Enhanced Lesion-to-Bubble Ratio on Ultrasonic Nakagami Imaging for Monitoring of High-Intensity Focused Ultrasound", 《J ULTRASOUND MED 2014》 *
SIYUAN ZHANG 等: "Feasibility of using Nakagami distribution in evaluating the formation of ultrasound-induced thermal lesions", 《2012 ACOUSTICAL SOCIETY OF AMERICA》 *
宋晓阳: "《医学超声图像的轮廓波方法研究及其在相控HIFU治疗***中的应用》", 《上海交通大学博士学位论文》 *
张思远等: ""高强度聚焦超声热消融中参量监控成像研究"", 《中国超声医学工程学会超声诊疗、生物效应、仪器工程、重庆超声医学工程学会学术会议论文集》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108553763A (zh) * 2018-01-19 2018-09-21 北京工业大学 一种基于超声回波去相关成像技术的微波热疗监测方法
CN112566694A (zh) * 2018-05-25 2021-03-26 卡尔迪亚韦弗公司 包括用于使空化泡成像的装置的超声处理设备
CN109171998A (zh) * 2018-10-22 2019-01-11 西安交通大学 基于超声深度学习的热消融区域识别监测成像方法与***
US11058473B2 (en) 2018-10-22 2021-07-13 Xi'an Jiaotong University Recognition, monitoring, and imaging method and system for thermal ablation region based on ultrasonic deep learning
CN113117258B (zh) * 2019-12-30 2023-04-28 重庆融海超声医学工程研究中心有限公司 组织凝固性坏死的检测装置
CN113117258A (zh) * 2019-12-30 2021-07-16 重庆融海超声医学工程研究中心有限公司 组织凝固性坏死的检测装置及方法
CN112966212A (zh) * 2021-02-07 2021-06-15 西安交通大学 基于超声回波背向散射能量变化的多参量实时监控成像***
CN113907792A (zh) * 2021-08-31 2022-01-11 西安交通大学 非线性解相关融合超声多参量的热凝固检测与成像方法
CN113907792B (zh) * 2021-08-31 2022-08-05 西安交通大学 非线性解相关融合超声多参量的热凝固检测与成像方法
CN114305668A (zh) * 2021-12-22 2022-04-12 聚融医疗科技(杭州)有限公司 基于解调域参量成像的超声热消融多参数监测方法及***
CN114305668B (zh) * 2021-12-22 2023-10-31 聚融医疗科技(杭州)有限公司 基于解调域参量成像的超声热消融多参数监测方法及***
CN114587565A (zh) * 2022-03-01 2022-06-07 河南中医药大学 一种射频消融中温度控制方法及***
CN114587565B (zh) * 2022-03-01 2023-09-29 河南中医药大学 一种射频消融中温度控制方法及***

Similar Documents

Publication Publication Date Title
CN106037815A (zh) 一种用于热凝固监测的超声回波统计参量成像***与方法
US11058399B2 (en) Bubble-induced color doppler feedback during histotripsy
US20070049827A1 (en) Clinical feedback of ablation efficacy during ablation procedure
DE69634714T2 (de) Therapeutisches Ultraschallgerät
JP5150267B2 (ja) 肝臓の障害を検出する超音波診断イメージングシステム
DE69333866T2 (de) Gerät zur medizinischen Ultraschallbehandlung unter Verwendung von Computertomographie
DE19548988C2 (de) Ultraschall-Diagnosesystem
JP4682149B2 (ja) 血流及び潅流パラメータを同時に表示するための超音波イメージングシステムおよび方法
US20110208061A1 (en) Ultrasonic lesion identification using temporal parametric contrast images
CN101623203B (zh) 一种瞬态物理过程中多模式和多参量同步检测成像监控***及监控方法
US5419332A (en) Mapping of flow parameters
US20130231565A1 (en) Ultrasound imaging system and method with automatic adjustment and/or multiple sample volumes
Mast et al. Ultrasound monitoring of in vitro radio frequency ablation by echo decorrelation imaging
CN104203115B (zh) 超声成像的方法和***
US20170011508A1 (en) Three-dimensional cavitation quantitative imaging method for microsecond-resolution cavitation spatial-temporal distribution
CN103123721B (zh) 一种实时减少图像中伪影的方法以及装置
CN107260217A (zh) 用于脑部聚焦超声空化实时监控的三维无源成像方法及***
Miller et al. Bubble-induced color Doppler feedback for histotripsy tissue fractionation
CN106264722A (zh) 窗口调制结合高斯多项式拟合监测射频消融设备及方法
KR20090081340A (ko) 의료 진단 초음파 이미지에 대한 조영제 파괴 효과성 결정
WO2007088477A2 (en) Improvements in methods for diagnosis and treatment of vessel occlusion
Chen et al. In vivo 3D power Doppler imaging using continuous translation and ultrafast ultrasound
Gu et al. Parametric perfusion imaging based on low-cost ultrasound platform
Rahimian et al. An acoustic backscatter-based method for estimating attenuation towards monitoring lesion formation in high intensity focused ultrasound
CN103220982B (zh) 用于检测对象的性质的装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161026

RJ01 Rejection of invention patent application after publication