CN101623203A - 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 - Google Patents
一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 Download PDFInfo
- Publication number
- CN101623203A CN101623203A CN200910023529A CN200910023529A CN101623203A CN 101623203 A CN101623203 A CN 101623203A CN 200910023529 A CN200910023529 A CN 200910023529A CN 200910023529 A CN200910023529 A CN 200910023529A CN 101623203 A CN101623203 A CN 101623203A
- Authority
- CN
- China
- Prior art keywords
- master control
- sample
- data
- ultrasonic
- acquisition equipment
- 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
Images
Abstract
本发明公开了一种瞬态物理过程中多模式和多参量同步检测成像监控系统，在主控PC上连接有数据采集设备、指示灯、采集卡和温度采集设备；阵列换能器通过全数字化B超与数据采集设备连接；热电偶连接在温度采集设备上；检测换能器通过接收器与采集卡连接；样品附近设置有超声换能器、摄像机、指示灯、检测换能器和热电偶；超声换能器采集到样品7的射频数据通过数据采集设备输入到主控PC并保存在主控PC上。本发明将热消融治疗系统与包括声学、光学和温度测量在内的多种监控系统集成同步，实现了瞬态物理过程中多模式、多参量同步检测成像监控。
Description
技术领域
本发明属于检测与监控技术领域，涉及到一种适合瞬态物理过程中多模式、多参量同步检测和成像监控于一体的系统与方法。
背景技术
在现有的检测与监控技术中，超声由于其成本低，实时性好，易于结合等优点已成为目前使用较广泛的检测与监控方式。已有国内外的专利使用传统B模式超声图像对肿瘤热消融治疗过程进行监控成像，例如美国专利US6425867B1和中国专利CN1565671A等。B模式超声图像反映是治疗过程中声阻抗差的变化，不能反映肿瘤热消融治疗中复杂的瞬态物理的变化。美国的S.Vaezy和英国的G.R.ter Haar等人对比了肿瘤热消融治疗前后治疗区域超声组织定征参量的变化，
但目前对于肿瘤热消融治疗过程中瞬态物理的研究不足。现有的检测与监控模式和参量单一，特别是缺少适合瞬态物理过程中对多种参量变化的同步检测与监控技术。通过对瞬态过程中发生的温度升高与热膨胀、不可逆组织损伤与瞬态变化、凝固性坏死，超声组织定征参量改变、空化汽化微泡及其活动等瞬态物理进行研究和监控，可以深入了解瞬态物理机制和推进同步动态检测监控的研究。
发明内容
针对上述检测与监控中的不足，本发明涉及到了一种以全数字化超声设备为基础的适合瞬态物理过程中多模式、多参量同步检测和成像监控的系统与方法。该发明同步检测瞬态物理过程中空化汽化微泡及其活动、超声组织定征参量和温度的动态变化；并利用光学图像、B模式超声图像和超声组织定征参量及其减影图像对目标区域进行同步动态成像监控。
为了实现上述任务，本发明采取如下的技术解决方案：
一种瞬态物理过程中多模式和多参量同步检测成像监控系统，包括主控PC1、数据采集设备2、全数字化B超3、超声换能器4、摄像机5、指示灯6、样品7、采集卡8、接收器9、检测换能器10、温度采集设备11和热电偶12；在主控PC1上连接有数据采集设备2、指示灯6、采集卡8和温度采集设备11；阵列换能器4通过全数字化B超3与数据采集设备2连接；热电偶12连接在温度采集设备11上；检测换能器10通过接收器9与采集卡8连接；样品7附近设置有超声换能器4、摄像机5、指示灯6、检测换能器10和热电偶12；按同步时序，超声换能器4采集到样品7的射频数据通过数据采集设备2输入到主控PC1并保存在主控PC1上；按同步时序，检测换能器10采集到样品7的射频数据通过接收器9和采集卡8输入到主控PC1并保存在主控PC1上；按同步时序，热电偶12采集到样品7的温度数据通过温度采集设备11输入到主控PC1并保存在主控PC1上；按同步时序，摄像机5进行拍摄，指示灯6受主控PC1的控制，而从在光学图像中表征同步时序。
所述同步时序是指主控PC1控制数据采集设备2、指示灯6、采集卡8和温度采集设备11同时开始工作的时间点为初始时间，以初始时间为基点按照相同的时间步长采集数据的方式。
一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法：
(1)将热电偶12插入样品7；按同步时序，热电偶12采集到样品7的温度数据通过温度采集设备11输入到主控PC1并保存在主控PC1上；
(2)将超声换能器4对准样品7；按同步时序，超声换能器4采集到样品7的射频数据通过数据采集设备2输入到主控PC1并保存在主控PC1上；主控PC1上的射频数据按照对数普差法计算得到样品7的平均衰减系数，主控PC1上的射频数据按照背向散射普法计算得到样品7的平均背向散射积分；
(3)将步骤(1)中得到的温度数据和步骤(2)中得到的平均衰减系数及平均背向散射积分进行比较，得到样品7的温度与射频数据的对应关系，由对应关系判断样品7的组织定征参量。
一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法：
(1)将检测换能器10对准样品7；按同步时序，检测换能器10同步接收空化汽化微泡及其活动散射的超声回波射频信号，超声回波射频信号由宽带接收器放大，再经由采集卡8采集到保存在主控PC1中；主控PC1中的每段超声回波射频信号经FFT转换到频域，通过频域提取其次谐波和宽带噪声，根据每段信号对应的时间转换到时域，获得稳态空化随时间的变化参数和惯性空化随时间的变化参数；
(2)将超声换能器4对准样品7；按同步时序，超声换能器4采集到样品7的射频数据通过数据采集设备2输入到主控PC1并保存在主控PC1上；对主控PC1上的射频数据构建B超声图像和背向散射积分减影图像；
(3)将步骤(1)中得到的稳态空化随时间的变化参数与惯性空化随时间的变化参数和步骤(2)中得到的B超声图像和背向散射积分减影图像进行比较，得到样品7的空化与射频数据的对应关系，由对应关系判断样品7的组织定征参量。
一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法：
(1)将超声换能器4对准样品7；按同步时序，超声换能器4采集到样品7的射频数据通过数据采集设备2输入到主控PC1并保存在主控PC1上；
主控PC1上的射频数据按照对数普差法计算得到样品7的平均衰减系数，主控PC1上的射频数据按照背向散射普法计算得到样品7的平均背向散射积分；
(2)由步骤(1)中得到的平均衰减系数构建衰减成像和衰减减影成像；由步骤(1)中得到的平均背向散射积分构建平均背向散射积分成像和平均背向散射积分减影成像；
(3)根据步骤(2)的结果判断样品7的组织定征参量。
一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法：
(1)将检测换能器10对准样品7；按同步时序，检测换能器10同步接收空化汽化微泡及其活动散射的超声回波射频信号，超声回波射频信号由宽带接收器放大，再经由采集卡8采集到保存在主控PC1中；主控PC1中的每段超声回波射频信号经FFT转换到频域，通过频域提取其次谐波和宽带噪声，根据每段信号对应的时间转换到时域，获得稳态空化随时间的变化参数和惯性空化随时间的变化参数；
(2)将超声换能器4对准样品7；按同步时序，超声换能器4采集到样品7的射频数据通过数据采集设备2输入到主控PC1并保存在主控PC1上；对主控PC1上的射频数据构建B超声图像和背向散射积分减影图像；
(3)将热电偶12插入样品7；按同步时序，热电偶12采集到样品7的温度数据通过温度采集设备11输入到主控PC1并保存在主控PC1上；
(4)将摄像机5对准样品7、指示灯6设置在摄像机5旁边，主控PC1按同步时序控制摄像机5和指示灯6工作；
(3)将步骤(1)、步骤(2)、步骤(3)和步骤(4)中得到的结果进行比较，得到样品7的组织定征参量。
本发明的集空化汽化微泡及其活动检测、多区域超声组织定征参量瞬态变化检测、动态光学信息采集和多点温度检测于一体的瞬态物理过程中多模式、多参量同步检测和成像监控系统。
该检测监控系统可以与治疗设备同步工作，并控制全数字化超声设备采集超声回波射频数据，以及同步采集微泡活动散射的回波信号、温度数据和光学信息，用于对瞬态物理过程中多模式、多参量同步检测和成像监控。该系统分成4大部分，分别是微泡活动检测系统、射频数据采集系统、光学信息采集系统和温度监测系统。微泡活动检测系统由检测换能器、接收器和高速数据采集卡组成。检测换能器被动接收空化汽化微泡及其活动散射的超声回波射频信号，并由接收器放大，再经由高速数据采集卡采集保存到电脑硬盘中。射频数据采集系统包括全数字化超声设备和数据采集设备。全数字化超声设备通过前端超声换能器采集超声回波射频数据，数据采集设备将超声设备采集的射频数据传输和保存到电脑硬盘中。光学信息采集系统包括高清录像机和照相机，可以采集透明仿生物组织体模中损伤、空化汽化微泡及其活动的动态变化过程。温度监测系统由前端热电偶和后端数据采集设备组成。热电偶采集温度变化引起的电压变化信息，并通过数据采集设备转化为温度信息保存在电脑硬盘中。
适合瞬态物理过程中多区域超声组织定征参量与多点温度同步监测及损伤动态变化同步监控的方法。采用上述系统，同时监测瞬态物理过程中不同区域的超声组织定征参量和多点温度的同步变化，以及与损伤动态变化之间的同步对应关系。其操作步骤如下：
(1)将热电偶插入离体样品中，并调节测温点位于目标区域。
(2)将插有热电偶的离体样品放入监控系统中，同步监测多个位置的温度变化。
(3)使用全数字化超声设备同步采集目标区域超声回波射频数据。
(4)使用这些射频数据经过超声组织定征算法处理，同时对多个区域的超声组织定征参量的同步瞬态变化进行检测。
(5)将不同区域的温度和多种超声组织定征参量的监测与损伤动态变化进行同步对比分析，研究不同的损伤发展阶段下温度和超声组织定征参量的瞬态变化以及与损伤的相互影响。
一种适合瞬态物理过程的多模式空化汽化微泡及其活动同步监测与损伤动态变化同步监控的方法。
采用上述系统，同时使用B模式超声图，背向散射积分减影图，被动空化检测和光学图像对空化汽化微泡及其活动进行监测，并与损伤的动态变化同步对照。其操作步骤如下：
(1)将离体样品放入监控系统中。
(2)使用全数字化超声设备同步采集目标区域超声回波射频数据。
(3)使用这些射频数据同步构建动态变化的B模式超声图和背向散射积分减影图。
(4)使用空化检测换能器同步接收空化汽化微泡及其活动散射的超声回波射频信号，并由宽带接收器放大，再经由数据采集卡采集到保存在电脑硬盘中。
(5)空化声信号特征参数提取，即将所得时域波形转换到频域，提取次谐波和宽带噪声，分别表征稳态和惯性空化活动的信息。
(6)使用光学信息采集系统采集透明仿生物组织体模中空化汽化微泡及其活动的动态变化过程。
(7)利用B模式超声图，背向散射积分减影图，被动空化检测和光学图像对空化汽化微泡及其活动的监测与损伤动态变化进行同步对比分析，研究不同的损伤发展阶段下微泡的瞬态变化以及与损伤的相互影响。
一种适合瞬态物理过程的多模式、多参量同步动态监控成像方法。
采用上述系统，同时使用B模式超声图像，衰减参数和背向散射积分及其减影成像方法进行动态监控成像。其操作步骤如下：
(1)将离体样品放入监控系统中。
(2)使用全数字化超声设备同步采集目标区域超声回波射频数据。
(3)使用这些射频数据经过谱差法和回波功率谱法处理，同时对衰减和背向散射积分的同步瞬态变化进行检测。
(4)利用超声组织定征参量的瞬态变化同步构建动态变化的B模式超声图像，背向散射积分和衰减参数及其减影图像。
(5)将多参量、多模式同步动态监控图像对比分析，研究不同的损伤发展阶段下多参量、多模式同步动态监控图像。
一种针对透明仿生物组织体模的瞬态物理过程光学和声学同步动态监控成像方法。
采用上述系统，针对透明仿生物组织体模，同时使用光学图像、B模式超声图像，背向散射积分和衰减参数及其减影成像方法进行动态监控成像。其操作步骤如下：
(1)将离体样品放入监控系统中。
(2)使用全数字化超声设备同步采集目标区域超声回波射频数据。
(3)使用这些射频数据经过谱差法和回波功率谱法处理，同时对不同区域的声衰减和背向散射积分参量的同步瞬态变化进行检测。
(4)利用超声组织定征参量的瞬态变化同步构建动态变化的B模式超声图像，背向散射积分和衰减参数及其减影图像。
(5)使用光学信息采集系统可以采集聚焦透明仿生物组织体模中损伤、空化汽化微泡及其活动的动态变化过程。
(6)将光学和多参量、多模式声学动态同步监控图像与损伤动态变化进行同步对比分析，研究不同的损伤发展阶段下光学和声学监控图像的变化以及与损伤的对应关系。
本发明与现有技术相比，具有下列优点：
与传统的监控系统相比，本发明将热消融治疗系统与包括声学、光学和温度测量在内的多种监控系统集成同步，实现了瞬态物理过程中多模式、多参量同步检测成像监控。
与传统的检测方法相比，本发明同步检测多区域超声组织定征参量、空化汽化微泡及其活动和多点温度的变化，实现了瞬态物理过程的多参量同步监测与损伤动态变化同步监控。
与传统的监控成像模式相比，本发明同时使用动态的光学图像、B模式超声图像、衰减参数和背向散射积分及其减影图像对目标区域进行监控成像，实现了瞬态物理过程中多模式、多参量同步动态监控成像。
附图说明
图1是本发明提出的瞬态物理过程中多模式、多参量同步检测成像监控系统的框图；
图2是多模式、多参量同步检测成像监控系统软件的工作时序图；
图3是瞬态物理过程中多区域超声组织定征参量与多点温度同步检测方法的流程图；
图4是瞬态物理过程中多模式空化汽化微泡及其活动同步检测方法的流程图；
图5是超声组织定征参量数字减影成像方法的流程图；
图6是肿瘤热消融瞬态物理过程中多区域超声组织定征参量与多点温度同步监测及损伤动态变化同步监控的结果；
图7是肿瘤热消融瞬态物理过程中多模式空化汽化微泡及其活动同步检测与损伤动态变化同步监控的结果；
图8是肿瘤热消融瞬态物理过程中多模式、多参量同步动态监控成像图；
图中所示的标号分别为：主控电脑(PC)1、数据采集设备2、全数字化B超3、超声换能器4、摄像机5、指示灯6、样品7、采集卡8、接收器9、检测换能器10、温度采集设备11、热电偶12。
具体实施方式
为了更清楚的理解本发明，以下结合附图和发明人给出的实施方式作进一步的详细说明。
一种集空化汽化微泡及其活动检测、多区域超声组织定征参量瞬态变化检测、动态光学信息采集和多点温度监测于一体的瞬态物理过程中多模式、多参量同步检测成像监控系统。
参见图1，本发明所涉及的适合瞬态物理过程多模式、多参量同步检测成像监控系统包括：主控电脑PC1、数据采集设备2、全数字化B超3、超声换能器4、摄像机5、指示灯6、样品7、采集卡8、接收器9、检测换能器10、温度采集设备11、热电偶12。主控PC1与数据采集设备2、指示灯6、采集卡8、温度采集设备11相连；全数字化B超3与数据采集设备2、阵列换能器4相连；接收器9与采集卡8、换能器10相连；温度采集设备11与热电偶12相连。
主控PC1控制全数字化B超将其超声换能器4采集到样品7的射频数据输出到数据采集设备2，最终使其保存于主控PC1上。主控PC1触发数据采集卡8，按照预制的时序采集来自换能器10的射频数据。温度采集设备11采集来自热电偶12的温度数据。摄像机5进行拍摄，指示灯6受主控PC1的控制，而从在光学图像中表征时序。
整个系统的工作时序如图2所示。系统开始工作时，主控PC1触发采集卡8，使其按预制的时序采集来自换能器10的射频信号，并且主控PC1还控制指示灯，使其点亮。此外，温度采集设备也从此刻开始对温度进行采集，摄像机开始进行拍摄。控PC1控制数据采集设备2采集相应帧数的B超数据。采集卡8按照预制的方式进行数据采集，例如每50ms采集1ms，即：T1＝50ms，T2＝1ms。
一种适合瞬态物理过程中多区域超声组织定征参量与多点温度同步监测及损伤动态变化同步监控的方法。
本发明同时监测瞬态物理过程中不同区域的超声组织定征参量和多点温度的同步瞬态变化，以及与损伤动态变化之间的同步对应关系。参见图3，主控电脑通过时序控制同步采离体样品内部的温度变化和超声回波射频数据。使用温度监控设备监测多区域的温度变化过程，温度监控设备前端为热电偶，通过穿刺针引导插入离体样品内部，通过B超图像观察引导定位，将热电偶尖端的测温点放置在目标区域。若是透明仿体，可直接通过肉眼观察，将热电偶尖端放置在目标点。热电偶连接后端的温度采集设备，测量和显示前端热电偶尖端的温度变化，并记录保存在主控电脑中。由于后端的温度采集设备有多路通道，可同时连接多路热电偶，达到对瞬态物理过程中多区域的温度变化同步监控。
使用全数字化超声设备同步采集超声回波射频数据，并按照预制将射频数据划分为多个目标区域，对每个区域的数据单独处理，便可以同时对多个区域的超声组织定征参量的同步瞬态变化进行检测。本发明中使用对数谱差法估计衰减系数。选取组织中单个区域内的射频数据按照谱差法处理，就可以得到该区域内的平均衰减系数。本发明中使用回波信号功率谱估计背向散射积分(IBS)。选取组织中单个区域内的射频数据按照功率谱算法处理，就可以得到该区域内的平均背向散射积分。
一种适合瞬态物理过程的多模式空化汽化微泡及其活动同步监测与损伤动态变化同步监控的方法。
本发明同时使用B模式超声图，背向散射积分减影图，被动空化检测和光学图像对空化汽化微泡及其活动进行监测，并与损伤的动态变化同步对比分析。参见图4，具体方式是使用全数字化超声设备同步采集超声回波射频数据，使用这些射频数据同步构建动态变化的常规B模式超声图和背向散射积分减影图。常规B模式超声图反映了组织内部声阻抗的差异，是目前普遍使用较为成熟的超声成像方法。而为了消除组织结构不均及声束扩散效应的影响，采用与减影法相结合的成像方法对目标区域的变化进行成像。背向散射积分的减影值ΔIBS计算如下：
其中，Tb和Ta分别表示聚焦超声照射不同时刻的成像时间，P(f，Ta)和P(f，Tb)分别为不同时刻的组织背向散射信号功率谱，f0为换能器中心频率，Δf为计算时所选的半带宽。背向散射积分减影成像通过滑动窗技术来实现。选取某一窗函数并设置一定窗宽(n)，按照一定的步长，将每一扫描线的射频数据划分为若干具有一定重叠部分的数据段。利用上式，计算肿瘤热消融过程中不同时刻所采集的两帧数据中相应数据段的背向散射积分减影值，将计算所得值作为该数据段中心处的值。这样就得到一幅背向散射积分减影图像。空化汽化微泡及其活动通过增强对超声的散射、反射等，改变超声传播特性，从而引起B模式超声图和背向散射积分减影图上的变化，因此我们可以通过B模式超声图和背向散射积分减影图对空化汽化微泡及其活动进行同步动态监控。
使用声学空化检测方法对微泡及其活动同步监测。参见图4，具体方式是使用空化检测换能器同步接收空化汽化微泡及其活动散射的超声回波射频信号，并由宽带接收器放大，再经由数据采集卡采集到保存在电脑硬盘中。从接收到的射频数据中提取空化声信号特征参数，具体是将所得每段时域信号通过FFT转换到频域，提取其次谐波和宽带噪声，分别表征该段时间内的稳态和惯性空化活动的信息。所有时域信号按照上述方法处理，并按照每段信号的时间坐标返回时域，便得到了微泡活动瞬态变化曲线。因为仿体是透明的，所以可以观察到仿体内部的动态微泡活动，所以使用光学信息采集系统采集透明仿生物组织体模中空化汽化微泡及其活动的动态变化过程，
一种适合瞬态物理过程的多模式、多参量同步动态监控成像方法。
本发明将多参量、多模式同步动态监控图像对比分析，研究不同的损伤发展阶段下多参量、多模式同步动态监控图像。具体实施方式为同时使用B模式超声图像，衰减参数和背向散射积分及其减影成像方法进行动态监控成像。参见图5，使用全数字化超声设备同步采集超声回波射频数据。使用这些射频数据经过超声组织定征算法处理，同时对声衰减和背向散射积分等超声组织定征参量的同步瞬态变化进行检测。利用超声组织定征参量的瞬态变化同步构建动态变化的B模式超声图像，衰减参数和背向散射积分及其减影图像。衰减参数的减影值Δα计算如下：
Δα＝α(Ta)-α(Tb)
其中，Tb和Ta分别表示不同时刻的成像时间。衰减参数减影成像仍然可通过滑动窗技术来实现。选取某一窗函数并设置一定窗宽(n)，按照一定的步长，将每一扫描线的射频数据划分为若干具有一定重叠部分的数据段。利用上式，计算不同时刻所采集的两帧数据中相应数据段的衰减参数减影值，将计算所得值作为该数据段中心处的值。这样就得到一幅衰减参数减影图像。
一种针对透明仿生物组织体模的瞬态物理过程光学和声学同步动态监控成像方法。
本发明针对透明仿生物组织体模，同时使用光学图像、B模式超声图像，衰减参数和背向散射积分及其减影图像进行动态监控成像。使用全数字化超声设备同步采集超声回波射频数据，使用这些射频数据经过超声组织定征算法处理，同时对声衰减、背向散射积分和非线性等超声组织定征参量的同步瞬态变化进行检测。利用超声组织定征参量的瞬态变化同步构建动态变化的B模式超声图像，衰减参数和背向散射积分及其减影图像。同时使用光学信息采集系统可以采集透明仿生物组织体模中损伤、空化汽化微泡及其活动的动态变化过程。将光学和多参量、多模式声学动态同步监控图像与损伤动态变化进行同步对比分析，研究不同的损伤发展阶段下光学和声学监控图像的变化以及与损伤的对应关系。
以下是发明人给出的实施例，但并不局限于这些实施例。
实施例1：肿瘤热消融(包括微波、高强度聚焦超声以及激光)治疗瞬态物理过程中，对透明仿体内部多区域的温度和超声组织定征参量的同步监测与损伤动态变化。肿瘤热消融治疗瞬态物理过程中，通过肉眼观察，将若干根(本实施例中使用4根)热电偶尖端放置在透明仿体内部预定位置。热电偶连接后端的温度采集设备，测量和显示前端这些热电偶尖端的温度变化，并记录保存在电脑中。通过主控电脑的时序控制，温度数据从肿瘤热消融治疗开始，每间隔50ms采集一次，直至治疗结束后一段时间。
本实施例中使用聚焦超声进行热消融治疗。使用声功率为100瓦，照射30秒。热消融治疗过程中，使用全数字化超声设备同步采集肿瘤热消融治疗过程中治疗区域超声回波射频数据。按照目标区域，将射频数据划分为多个模块，每个区域内的射频数据分别经过处理，可同时对肿瘤热消融治疗过程中多个区域的声衰减和背向散射积分等超声组织定征参量的同步瞬态变化进行检测。图6是肿瘤热消融治疗瞬态物理过程中，多区域的温度与超声组织定征参量同步监测与损伤动态变化同步监控结果。肿瘤热消融治疗开始后，治疗区域的温度上升，但是不同区域的温度随损伤发展的变化不同。同时，不同治疗区域的声学参数，例如衰减参数和背向散射参数随损伤发展的变化也不同。
实施例2：肿瘤热消融(包括微波、高强度聚焦超声以及激光)治疗瞬态物理过程中，多模式空化监测与损伤的同步变化。肿瘤热消融治疗瞬态物理过程中，同时使用B模式超声图，背向散射积分减影图，被动空化检测和光学图像对肿瘤热消融治疗过程中空化汽化微泡及其活动进行监测，并与损伤的动态变化同步对照，如图7所示。
本实施例中使用聚焦超声进行热消融治疗。使用声功率为80瓦，连续照射50秒。在热消融治疗过程中，使用全数字化超声设备同步采集肿瘤热消融治疗过程中治疗区域超声回波射频数据。使用这些射频数据同步构建肿瘤热消融治疗过程中动态变化的背向散射积分减影图(图7a)和常规B模式超声图(图7b)。空化汽化微泡及其活动通过增强对超声的散射、反射等，改变超声传播特性，从而引起B模式超声图和背向散射积分减影图上的变化，因此我们可以通过B模式超声图和背向散射积分减影图对空化汽化微泡及其活动进行同步动态监控。热消融治疗过程中，同步使用声学空化检测方法对微泡及其活动同步监测。具体方式是使用空化检测换能器同步接收空化汽化微泡及其活动散射的超声回波射频信号，并由宽带接收器放大，再经由数据采集卡采集到保存在电脑硬盘中。从接收到的射频数据中提取空化声信号特征参数，具体是将所得波形通过FFT转换到频域，提取次谐波和宽带噪声，分别表征稳态和惯性空化活动的信息。热消融治疗过程中，使用光学信息采集系统采集透明仿生物组织体模中空化汽化微泡及其活动的动态变化过程。因为仿体是透明的，所以在热消融治疗过程中，我们可以使用录像机拍摄仿体内部的动态微泡活动(图7c和d)。肿瘤热消融治疗过程中B模式超声图，背向散射积分减影图，被动空化检测和光学图像对肿瘤热消融治疗过程中空化汽化微泡及其活动的监测与损伤动态变化进行同步对比分析，研究不同的损伤发展阶段下微泡的瞬态变化以及与损伤的相互影响。
实施例3：肿瘤热消融(包括微波、高强度聚焦超声以及激光)治疗过程中，多模式、多参量同步动态监控成像。在热消融治疗过程中，同时使用光学图像、B模式超声图像、背向散射积分减影成像方法对肿瘤热消融过程进行动态监控成像，如图8所示。
本实施例中使用聚焦超声进行热消融治疗。使用声功率为65瓦，占空比方式照射总时间为41秒。热消融治疗过程中，使用光学信息采集系统可以采集聚焦超声热消融治疗过程中的透明仿生物组织体模中损伤、空化汽化微泡及其活动的动态变化过程，如图8(a)所示。热消融治疗过程中，使用全数字化超声设备同步采集肿瘤热消融治疗过程中治疗区域超声回波射频数据。使用这些射频数据经过算法程序处理，同时对肿瘤热消融治疗过程中不同区域的声衰减、背向散射积分和非线性等超声组织定征参量的同步瞬态变化进行检测。利用超声组织定征参量在肿瘤热消融治疗过程的瞬态变化同步构建肿瘤热消融治疗过程中动态变化的B模式超声图像(图8b)和背向散射积分减影图像(图8c)。将肿瘤热消融治疗过程中光学和多参量、多模式声学动态同步监控图像与损伤动态变化进行同步对比分析，可研究不同的损伤发展阶段下光学和声学监控图像的变化以及与损伤的对应关系。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施方式仅限于此，对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单的推演或替换，都应当视为属于本发明由所提交的权利要求书确定专利保护范围。
Claims (6)
1、一种瞬态物理过程中多模式和多参量同步检测成像监控系统，其特征在于：包括主控PC(1)、数据采集设备(2)、全数字化B超(3)、超声换能器(4)、摄像机(5)、指示灯(6)、样品(7)、采集卡(8)、接收器(9)、检测换能器(10)、温度采集设备(11)和热电偶(12)；在主控PC(1)上连接有数据采集设备(2)、指示灯(6)、采集卡(8)和温度采集设备(11)；阵列换能器4通过全数字化B超(3)与数据采集设备(2)连接；热电偶(12)连接在温度采集设备(11)上；检测换能器(10)通过接收器(9)与采集卡(8)连接；样品(7)附近设置有超声换能器(4)、摄像机(5)、指示灯(6)、检测换能器(10)和热电偶(12)；按同步时序，超声换能器(4)采集到样品7的射频数据通过数据采集设备(2)输入到主控PC(1)并保存在主控PC(1)上；按同步时序，检测换能器(10)采集到样品7的射频数据通过接收器(9)和采集卡(8)输入到主控PC(1)并保存在主控PC(1)上；按同步时序，热电偶(12)采集到样品7的温度数据通过温度采集设备(11)输入到主控PC(1)并保存在主控PC(1)上；按同步时序，摄像机(5)进行拍摄，指示灯(6)受主控PC(1)的控制，而从在光学图像中表征同步时序。
2、根据权利要求1所述一种瞬态物理过程中多模式和多参量同步检测成像监控系统，其特征在于：所述同步时序是指主控PC(1)控制数据采集设备(2)、指示灯(6)、采集卡(8)和温度采集设备(11)同时开始工作的时间点为初始时间，以初始时间为基点按照相同的时间步长采集数据。
3、根据权利要求1所述一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法，其特征在于：
(1)将热电偶(12)插入样品(7)；按同步时序，热电偶(12)采集到样品(7)的温度数据通过温度采集设备(11)输入到主控PC(1)并保存在主控PC(1)上；
(2)将超声换能器(4)对准样品(7)；按同步时序，超声换能器(4)采集到样品(7)的射频数据通过数据采集设备(2)输入到主控PC(1)并保存在主控PC(1)上；主控PC(1)上的射频数据按照对数普差法计算得到样品(7)的平均衰减系数，主控PC(1)上的射频数据按照背向散射普法计算得到样品(7)的平均背向散射积分；
(3)将步骤(1)中得到的温度数据和步骤(2)中得到的平均衰减系数及平均背向散射积分进行比较，得到样品(7)的温度与射频数据的对应关系，由对应关系判断样品(7)的组织定征参量。
4、根据权利要求1所述一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法，其特征在于：
(1)将检测换能器(10)对准样品(7)；按同步时序，检测换能器(10)同步接收样品(7)的空化汽化微泡及其活动散射的超声回波射频信号，超声回波射频信号由宽带接收器放大，再经由采集卡(8)采集到保存在主控PC(1)中；主控PC(1)中的每段超声回波射频信号经FFT转换到频域，通过频域提取其次谐波和宽带噪声，根据每段信号对应的时间转换到时域，获得稳态空化随时间的变化参数和惯性空化随时间的变化参数；
(2)将超声换能器(4)对准样品(7)；按同步时序，超声换能器(4)采集到样品(7)的射频数据通过数据采集设备(2)输入到主控PC(1)并保存在主控PC(1)上；对主控PC(1)上的射频数据构建B超声图像和背向散射积分减影图像；
(3)将步骤(1)中得到的稳态空化随时间的变化参数与惯性空化随时间的变化参数和步骤(2)中得到的B超声图像和背向散射积分减影图像进行比较，得到样品(7)的空化与射频数据的对应关系，由对应关系判断样品(7)的组织定征参量。
5、根据权利要求1所述一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法，其特征在于：
(1)将超声换能器(4)对准样品(7)；按同步时序，超声换能器(4)采集到样品(7)的射频数据通过数据采集设备(2)输入到主控PC(1)并保存在主控PC(1)上；
主控PC(1)上的射频数据按照对数普差法计算得到样品(7)的平均衰减系数，主控PC(1)上的射频数据按照背向散射普法计算得到样品(7)的平均背向散射积分；
(2)由步骤(1)中得到的平均衰减系数构建衰减成像和衰减减影成像；由步骤(1)中得到的平均背向散射积分构建平均背向散射积分成像和平均背向散射积分减影成像；
(3)根据步骤(2)的结果判断样品(7)的组织定征参量。
6、根据权利要求1所述一种瞬态物理过程中多模式和多参量同步检测成像监控系统的监控方法，其特征在于：
(1)将检测换能器(10)对准样品(7)；按同步时序，检测换能器(10)同步接收样品(7)的空化汽化微泡及其活动散射的超声回波射频信号，超声回波射频信号由宽带接收器放大，再经由采集卡(8)采集到保存在主控PC(1)中；主控PC(1)中的每段超声回波射频信号经FFT转换到频域，通过频域提取其次谐波和宽带噪声，根据每段信号对应的时间转换到时域，获得稳态空化随时间的变化参数和惯性空化随时间的变化参数；
(2)将超声换能器(4)对准样品(7)；按同步时序，超声换能器(4)采集到样品(7)的射频数据通过数据采集设备(2)输入到主控PC(1)并保存在主控PC(1)上；对主控PC(1)上的射频数据构建B超声图像和背向散射积分减影图像；
(3)将热电偶(12)插入样品(7)；按同步时序，热电偶(12)采集到样品(7)的温度数据通过温度采集设备(11)输入到主控PC(1)并保存在主控PC(1)上；
(4)将摄像机(5)对准样品(7)、指示灯(6)设置在摄像机(5)旁边，主控PC(1)按同步时序控制摄像机(5)和指示灯(6)工作；
(5)将步骤(1)、步骤(2)、步骤(3)和步骤(4)中得到的结果进行比较，得到样品(7)的组织定征参量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100235294A CN101623203B (zh) | 2009-08-07 | 2009-08-07 | 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100235294A CN101623203B (zh) | 2009-08-07 | 2009-08-07 | 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101623203A true CN101623203A (zh) | 2010-01-13 |
| CN101623203B CN101623203B (zh) | 2011-02-09 |
Family
ID=41519302
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100235294A Expired - Fee Related CN101623203B (zh) | 2009-08-07 | 2009-08-07 | 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101623203B (zh) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102049109A (zh) * | 2011-01-06 | 2011-05-11 | 西安交通大学 | 一种聚焦超声场内流动微泡空化产生附加热量补偿流动热量损失的方法 |
| CN104156908A (zh) * | 2014-08-27 | 2014-11-19 | 广西大学 | 高时间分辨率的频域瞬态成像方法 |
| CN104190005A (zh) * | 2014-08-01 | 2014-12-10 | 深圳市普罗惠仁医学科技有限公司 | 一种基于仿生体系的超声聚焦评估系统及方法 |
| CN105380679A (zh) * | 2015-12-15 | 2016-03-09 | 飞依诺科技(苏州)有限公司 | 超声成像设备的控制方法及系统 |
| CN106037815A (zh) * | 2016-05-17 | 2016-10-26 | 西安交通大学 | 一种用于热凝固监测的超声回波统计参量成像系统与方法 |
| CN109171998A (zh) * | 2018-10-22 | 2019-01-11 | 西安交通大学 | 基于超声深度学习的热消融区域识别监测成像方法与系统 |
| CN109259822A (zh) * | 2018-11-07 | 2019-01-25 | 西安交通大学 | 基于被动空化检测成像的聚焦超声治疗三维引导与动态实时监控系统及方法 |
| CN110243454A (zh) * | 2019-07-08 | 2019-09-17 | 西安交通大学 | 粘弹媒质中基于双频叠加超声脉冲的微泡增效动力学实验系统及方法 |
| CN110313940A (zh) * | 2019-08-01 | 2019-10-11 | 无锡海斯凯尔医学技术有限公司 | 信号衰减计算方法、装置、设备及计算机可读存储介质 |
| CN111445754A (zh) * | 2020-04-30 | 2020-07-24 | 中国医学科学院生物医学工程研究所 | 一种模拟激光外科手术的辅助训练系统 |
| CN112566694A (zh) * | 2018-05-25 | 2021-03-26 | 卡尔迪亚韦弗公司 | 包括用于使空化泡成像的装置的超声处理设备 |
| CN113907792A (zh) * | 2021-08-31 | 2022-01-11 | 西安交通大学 | 非线性解相关融合超声多参量的热凝固检测与成像方法 |
| US11491352B2 (en) * | 2018-06-05 | 2022-11-08 | Korea Institute Of Science And Technology | High-low intensity focused ultrasound treatment apparatus |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4932414A (en) * | 1987-11-02 | 1990-06-12 | Cornell Research Foundation, Inc. | System of therapeutic ultrasound and real-time ultrasonic scanning |
| DE4207463C2 (de) * | 1992-03-10 | 1996-03-28 | Siemens Ag | Anordnung zur Therapie von Gewebe mit Ultraschall |
| US6533726B1 (en) * | 1999-08-09 | 2003-03-18 | Riverside Research Institute | System and method for ultrasonic harmonic imaging for therapy guidance and monitoring |
| CN100512908C (zh) * | 2003-06-11 | 2009-07-15 | 重庆海扶(Hifu)技术有限公司 | 高强度超声波治疗实时监控系统及方法 |
| CN100386057C (zh) * | 2005-10-31 | 2008-05-07 | 西安交通大学 | 基于包膜微泡的灌注成像与超声控制释放的系统和方法 |
| CN1847824B (zh) * | 2006-02-27 | 2011-01-05 | 西安交通大学 | 基于超高速摄影技术的超声场中微泡行为分析系统和方法 |
-
2009
- 2009-08-07 CN CN2009100235294A patent/CN101623203B/zh not_active Expired - Fee Related
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102049109A (zh) * | 2011-01-06 | 2011-05-11 | 西安交通大学 | 一种聚焦超声场内流动微泡空化产生附加热量补偿流动热量损失的方法 |
| CN104190005A (zh) * | 2014-08-01 | 2014-12-10 | 深圳市普罗惠仁医学科技有限公司 | 一种基于仿生体系的超声聚焦评估系统及方法 |
| CN104190005B (zh) * | 2014-08-01 | 2017-09-19 | 深圳市普罗惠仁医学科技有限公司 | 一种基于仿生体系的超声聚焦评估系统及方法 |
| CN104156908A (zh) * | 2014-08-27 | 2014-11-19 | 广西大学 | 高时间分辨率的频域瞬态成像方法 |
| CN104156908B (zh) * | 2014-08-27 | 2018-11-02 | 广西大学 | 高时间分辨率的频域瞬态成像方法 |
| CN105380679A (zh) * | 2015-12-15 | 2016-03-09 | 飞依诺科技(苏州)有限公司 | 超声成像设备的控制方法及系统 |
| CN105380679B (zh) * | 2015-12-15 | 2018-11-20 | 飞依诺科技(苏州)有限公司 | 超声成像设备的控制方法及系统 |
| CN106037815A (zh) * | 2016-05-17 | 2016-10-26 | 西安交通大学 | 一种用于热凝固监测的超声回波统计参量成像系统与方法 |
| CN112566694A (zh) * | 2018-05-25 | 2021-03-26 | 卡尔迪亚韦弗公司 | 包括用于使空化泡成像的装置的超声处理设备 |
| US11491352B2 (en) * | 2018-06-05 | 2022-11-08 | Korea Institute Of Science And Technology | High-low intensity focused ultrasound treatment apparatus |
| 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 |
| CN109259822A (zh) * | 2018-11-07 | 2019-01-25 | 西安交通大学 | 基于被动空化检测成像的聚焦超声治疗三维引导与动态实时监控系统及方法 |
| CN110243454B (zh) * | 2019-07-08 | 2020-06-19 | 西安交通大学 | 粘弹媒质中基于双频叠加超声脉冲的微泡增效动力学实验系统及方法 |
| CN110243454A (zh) * | 2019-07-08 | 2019-09-17 | 西安交通大学 | 粘弹媒质中基于双频叠加超声脉冲的微泡增效动力学实验系统及方法 |
| CN110313940A (zh) * | 2019-08-01 | 2019-10-11 | 无锡海斯凯尔医学技术有限公司 | 信号衰减计算方法、装置、设备及计算机可读存储介质 |
| CN111445754A (zh) * | 2020-04-30 | 2020-07-24 | 中国医学科学院生物医学工程研究所 | 一种模拟激光外科手术的辅助训练系统 |
| CN113907792A (zh) * | 2021-08-31 | 2022-01-11 | 西安交通大学 | 非线性解相关融合超声多参量的热凝固检测与成像方法 |
| CN113907792B (zh) * | 2021-08-31 | 2022-08-05 | 西安交通大学 | 非线性解相关融合超声多参量的热凝固检测与成像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101623203B (zh) | 2011-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101623203B (zh) | 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法 | |
| US5119819A (en) | Method and apparatus for non-invasive monitoring of blood glucose | |
| US6519485B2 (en) | Minimally invasive system for assessment of organ function | |
| DE60121179T2 (de) | Photoakustische prüfvorrichtung und bilderzeugungssystem | |
| US8226559B2 (en) | Method and apparatus for real-time temperature measuring for high-intensity focused ultrasound(HIFU) therapy system | |
| CN104135936B (zh) | 光声图像化方法及装置 | |
| JPS63122923A (ja) | 超音波測温装置 | |
| RU2007147913A (ru) | Оптоакустический мониторинг в реальном времени с помощью электрофизиологических катетеров | |
| CN101612034B (zh) | 重构混浊介质光学参数的时间分辨测量系统及方法 | |
| KR20090103408A (ko) | 초음파를 이용한 생체 조직 탄성도 및 경화 측정 시스템 | |
| CN106037818B (zh) | 一种光声和超声双模式骨质疏松症检测方法与装置 | |
| KR20080108918A (ko) | 초음파와 광의 상호 작용을 이용하여 피검체의 내부 정보를얻는 의료 기구 | |
| CN110361357B (zh) | 一种用于皮肤检测的单阵元光声谱信号获取系统及方法 | |
| CN106037815A (zh) | 一种用于热凝固监测的超声回波统计参量成像系统与方法 | |
| Kulkarni et al. | A feasibility study on noninvasive blood glucose measurement using photoacoustic method | |
| CN104706363A (zh) | 一种复合型光声无损动态血糖检测仪 | |
| CN107228904A (zh) | 一种光致超声的血糖无创检测装置及方法 | |
| EP3165167A1 (en) | Optimization and feedback control of hifu power deposition through the analysis of detected signal characteristics | |
| CN204666826U (zh) | 一种太赫兹连续安检成像装置 | |
| US8981300B2 (en) | Electromagnetic wave pulse measuring device and method, and application device using the same | |
| CN104116522B (zh) | 一种基于超声波衰减测量肌肉疲劳程度的装置和方法 | |
| CN102879335A (zh) | 一种便携式无创、实时光声粘弹检测仪 | |
| CN110960267B (zh) | 一种病理细胞检测用智能穿刺设备 | |
| IL96029A (en) | Method and instrument for characterizing returned ultrasonic pulses | |
| CN106725324A (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 | ||
| C14 | Grant of patent or utility model | ||
| 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: 20110209 Termination date: 20210807 |