CN104605891A - 检测剪切波在生物组织中传输速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 - Google Patents
检测剪切波在生物组织中传输速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 Download PDFInfo
- Publication number
- CN104605891A CN104605891A CN201410852544.0A CN201410852544A CN104605891A CN 104605891 A CN104605891 A CN 104605891A CN 201410852544 A CN201410852544 A CN 201410852544A CN 104605891 A CN104605891 A CN 104605891A
- Authority
- CN
- China
- Prior art keywords
- tissue
- shearing wave
- signal
- echo
- detected
- 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
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
- 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/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了检测剪切波在生物组织中传输速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法。本发明所述的检测剪切波在生物组织中传输速度的方法,其特征在于,利用待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在该两位置上实时回波信号与原始回波信号的MAC值产生变化的时间差,计算该两位置间剪切波传播的速度。本发明中采用模态置信准则对剪切波进行跟踪,可以实现组织在剪切波作用下的微形变估计算,提高图像的信噪比、分辨率、对比度;灵敏度高。另一方面,模态置信准则是对超声回波射频信号的时域计算方法,计算流程一致性较好,可以允许多项连续的模态置信因子计算任务在相同时间的并发操作即并行计算使整幅弹性图像的复杂计算操作获得较高的加速比。因此,该方法可以有效的提高弹性成像算法的计算速度并实现剪切波的传播过程的实时跟踪。
Description
技术领域
本发明涉及检测剪切波在生物组织中传输速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法。
背景技术
在医学诊断中,生物组织的弹性信息对于诊断疾病具有非常重要的参考价值。一般认为,软组织病理的改变往往与组织的硬度相关。如组织的硬度或弹性会随着肿瘤在组织中的生长或扩散而改变。甚至在肿瘤在出现形态学改变前,其硬度已经发生变化,导致正常组织与病变组织之间存在较大的弹性差异。这种利用生物组织的力学特性的差异,通过区分外力作用下不同生物组织的响应(应变、位移等)而进行成像的方法即为弹性成像(elastography)。
对目前的超声弹性成像技术来说,根据加载方式不同,可将其大致分为两类:静态/准静态超声弹性成像和动态超声弹性成像。目前发展较为成熟的超声弹性成像技术多属于静态/准静态超声弹性成像技术,其基本原理是利用采用手动加压的方法使待检目标发生变形,通过向变形前后的特检目标发射超声波并接收回波信号,经分析处理后获得位移、应变、应变率和速度等待检目标的响应参数,最终估计出待检目标的杨氏模量、剪切模量、泊松比和拉梅常数等材料力学属性的相对值,进而定性地了解目标不同组分的硬度差异,并根据这一差异诊断待检目标的病理状况。
但类似的超声弹性成像方法对成像的深度和位置都会有限制。手动加压的加载方式决定了这种超声弹性成像方法不能对深层组织进行成像,相应产生的应变与位移也会受到施加压力的大小、加载时间的长短等人为因素的影响。所以类似方法只能定性的显示浅表组织的弹性信息,重复性差。为了上述缺陷,可行的方法之一是采用聚焦超声束激励的方法引起组织运动,再利用超声成像技术探测组织的运动,通过组织运动特性重建组织体弹性参数分布,即动态超声弹性成像技术。作为一种具有代表性的动态超声弹性成像技术,剪切波弹性成像技术(shear wave elastography,SWE)是超声弹性成像技术的一种替代方案。其基本思想是通过检测声辐射力激发组织而产生的剪切波的传播特性进行成像的方法。一般是通过超声探头发射的强度在安全阈值内的声辐射脉冲在待检组织内形成声辐射力(Acoustic Radiation Force)对下层组织进行施压,压力在声波传播方向上推动组织,组织相应地产生复原力,该复原力会产生机械波,尤其是沿组织横向传播的剪切波。而机械波的传播性质与材料的力学特征紧密相关,剪切波作为组织复原力产生的一种机械横波,其传播性质必然包含着组织内的诸多力学信息。简单地说,剪切波速度越快,说明组织杨氏模量值越大,即组织的硬度越大。利用这一性质,即可通过组织内传播的剪切波速度反推得到组织的杨氏模量绝对值。所以对剪切波弹性成像方法来说,其核心技术是跟踪剪切波在组织内的传播过程并据此计算其传播速度。剪切波在生物体内传播速度缓慢,一般为1~10m/s,故可利用超声成像方法捕捉剪切波的传播过程。
对一般的动态超声弹性成像方法而言,其工作原理是通过检测声辐射力作用前后的超声回波信号,应用相应的变形估计算法对这些信号进行分析处理,进而根据胡克定理确定出待检对象弹性信息的空间分布。这类技术急待改进的方面主要体现在:改善***的性能水平,有效提高超声换能器的电声转换效率,输出高能量的超声波信号产生更大的超声辐射力,实现待检材料弹性参数较大动态范围的检测,以进一步拓展该技术的应用空间;开发利用超声回波实现组织微小形变检测的计算方法,提高图像的信噪比、分辨率、对比度及帧频进而实现组织力学参数空间分布的实时检测,以提高该技术的临床应用价值。
发明内容
本发明的目的之一是为了克服现有技术中的不足,提供一种可对深层生物组织检测的检测剪切波在生物组织中传输速度的方法。
为实现以上目的,本发明通过以下技术方案实现:
检测剪切波在生物组织中传输速度的方法,其特征在于,利用待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在该两位置上实时回波信号与原始回波信号的MAC值产生变化的时间差,计算该两位置间剪切波传播的速度。
优选地是,按照以下公式计算实时回波信号与原始回波信号之间的模态置信度MAC值:
其中,ωk(n)代表窗口长度为n的第k个窗函数,ψij代表在第i次发射波束后第j个通道接收到的回波信号,代表t=0时刻下,第j个通道采集的原始回波信号,MACijk代表通过ωk(n)截取的实时回波信号与ψij及的模态置信度。
优选地是,根据对待检测组织各位置对超声回波信号的MAC值,组成矩阵;并记录在不同时刻下待检测组织各位置对超声回波信号的MAC值,分别组成矩阵;依据待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在两位置上实时回波信号与原始回波信号的MAC值矩阵产生变化的时间差,计算该距离间剪切波传播的平均速度。
优选地是,所述的实时回波信号为实时超声回波信号;所述原始回波信号为原始超声回波信号。
优选地是,包括步骤:
①.待检测组织发射超声波束,并接收超声回波信息作为参考的原始组织回波信号,记为
②.对待检测组织发射聚焦波束,使待检测组织产生剪切波;
③.在待检测组织产生剪切波后,对待检测组织的待检测区域发射超声波束并接收反射回来的回波信号;
④.重复步骤③,对待检测区域多次发射超声波束,记录超声回波信号,记为ψij;直到剪切波传播至待检测区域之外;
⑤.通过窗函数ωk对步骤④中实时接收的超声回波信号ψij和步骤①中记录原始超声射频回波信号进行多次截取,计算MAC值。
优选地是,所述的步骤②中,对待检测组织发射声波,在组织内部产生声辐射力,利用声辐射力在声波传播方向上推动待检测组织,使待检测组织产生剪切波。
优选地是,所述的步骤①中,通过超声仪器的发射端(12)控制阵列超声换能器(13)对待检测组织(14)发射超声波束,并通过接收端(15)和数据采集模块(16)记录不同通 道接收到的超声回波信息作为参考的原始组织回波信号,记为并传回至数据存储器(17);
所述的步骤②中,通过发射端(12)控制阵列式超声换能器(13)对待检测组织(14)发射聚焦波束,在组织内部产生声辐射力,声辐射力在声波传播方向上推动待检测组织,使待检测组织产生剪切波;
所述的步骤③中,在声辐射力对待检测组织进行激励后,控制发射端(12)对待检测组织的待检测区域发射超声波束并接收反射回来的超声回波信号(15);
所述的步骤④中,重复步骤③,对待检测区域多次发射超声波束,通过数据采集模块(16)实时记录各通道接收到的超声射频回波信号并传回至数据存储器(17),记为ψij,直到剪切波传播至待检测区域之外;
所述的步骤⑤中,通过窗函数ωk(21)对步骤④中实时接收的超声射频回波信号ψij(23)和步骤①中记录原始超声射频回波信号(22)进行多次截取,并计算截取后的超声射频回波信号与原始回波信号之间的模态置信度MAC值。
优选地是,所述生物组织为人体组织。
优选地是,所述人体组织为深层组织。
本发明的目的之二是为了克服现有技术中的不足,提供一种可对深层生物组织检测的检测剪生物组织弹性的方法。
为实现以上目的,本发明通过以下技术方案实现:
检测生物组织弹性的方法,其特征在于,采用前述的方法检测剪切波在生物组织中的传播速度,再利用剪切波在生物组织中的传播速度与弹性的线性关系计算得到组织的弹性特性数据。
本发明的目的之三是为了克服现有技术中的不足,提供一种可对深层生物组织检测的生物组织弹性成像方法。
为实现以上目的,本发明通过以下技术方案实现:
生物组织弹性成像方法,其特征在于,利用前述的方法计算生物组织的弹性,并对弹性数据生成图像。优选为灰度图像。
本专利提供的基于模态置信准则获得剪切波弹性成像的方法,能对人体深层组织的弹性信息进行测量,具有实时、定量、快速、高分辨率等优点。其与现有超声弹性成像技术相比最大的不同在于本方法采用模态置信准则对受到剪切波影响的组织进行分析,以跟踪剪切波的传播过程。进而据此计算组织的弹性信息为目的进行组织弹性信息定量分析,通过比较正常组织与病变组织的弹性信息差异来实现组织病变的早期诊断。
本发明中涉及的剪切波弹性成像技术有别于静态弹性成像技术的接触式激励方法,利用超声辐射力作为驱动力源,声波穿透组织并在深层位置聚焦产生作用力,该作用力非接触的作用于深层组织上使其产生在深层组织内传播的剪切波,进而通过跟踪该剪切波的传播过程对深层组织的弹性信息进行测量。因此,现有的静态弹性成像技术多应用于人体浅表组织的弹性信息测量(如***等),而本发明涉及的剪切波弹性成像方法可对深层组织的弹性信息进行测量,尤其是距离表皮5cm、甚至是10cm以上的深层组织,如肝脏等。
本发明中采用模态置信准则对剪切波进行跟踪,可以实现组织在剪切波作用下的微形变估计算,提高图像的信噪比、分辨率、对比度;灵敏度高。另一方面,模态置信准则是对超声回波射频信号的时域计算方法,计算流程一致性较好,可以允许多项连续的模态置信因子计算任务在相同时间的并发操作即并行计算。我们能够利用数据并行性使整幅弹性图像的复杂计算操作获得较高的加速比。因此,该方法可以有效的提高弹性成像算法的计算速度并实现剪切波的传播过程的实时跟踪,计算速度快。
附图说明
图1为本发明使用的超声***结构示意图。
图2为本发明步骤①、②实施方法示意图。
图3为本发明步骤⑤实施方法示意图。
图4为本发明步骤⑥实施方法示意图。
图5为本发明中的MAC值灰度成像示意图。
具体实施方式
模态置信准则(Modal Assurance Criteria,MAC)一般是用来计算机械振动中的待测振型与参考振型的线性独立性。其定义可以表示如下:
其中,和分别是待辨识的波形信号特性向量和参考波形信号特性向量。
在利用这种方法进行任意波形信号的匹配分析中,首先需要对波形信号进行离散化,而后对离散后的波形信号构造其特性向量,再根据参考波形信号的特性向量和待检测的波形信号的特性向量构造总体模态置信度矩阵。模态置信度矩阵非对角元越小,则说明两者波形独立性越好,反之则意味着两者波形相关性越大。或者也可以根据参考波形信号的特性向量和待检测的波形信号的特性向量直接计算两者之间的模态置信度。计算得到的模态置信度的数值介于0至1之间,越接近0代表两者独立性越大,越接近1则代表两者相关性越大。
根据上述原理,本发明专利描述了一种利用模态置信准则(Modal Assurance Criteria,MAC)对待检组织内由声辐射力激励产生的剪切波进行捕捉并定量计算组织弹性信息的剪切波弹性成像方法,使用户可以得到定量显示组织弹性信息的图像。其基本原理如下:首先记录待检组织初始状态下超声回波信号,而后通过阵列超声换能器对待检测组织发射一个或多个推动脉冲波束,使组织受到激励产生剪切波;随后,对组织的待检测区域发射超声波束并实时接收反射回来的回波信号直到剪切波传播至待检测区域之外;通过比较待检测组织内各位置的原始回波信号与实时获得的回波信号在不同时刻下的模态置信度矩阵,可分析得到不同时刻下待检测组织各位置上的由剪切波引起的组织位移情况,追踪剪切波在组织中的传播过程,进而得到待测组织内各位置上的杨氏弹性模量绝对值,并映射得到弹性成像。其与现有方法的不同在于,本发明专利首次采用模态置信准则对待检组织内剪切波传播过程进行捕捉,进而获得组织各部分的定量弹性信息。再通过比较正常软组织与病***组织的弹性信息 差异进行诊断。该方法可有利地促进剪切波弹性成像在实际临床应用中的推广和应用,加强人们对组织病变过程及特性的认识,尤其是对早期(未扩散Ⅶ)肿瘤的预防和诊治具有很高的参考和实用价值。
以下结合附图对本发明进一步说明。
参考图1显示了根据本专利的原理构建的用于测量剪切波的超声***。
生物组织弹性成像方法,包括步骤:
①.首先通过发射端12控制阵列超声换能器13对待检测组织14发射超声波束,并通过接收端15和数据采集模块16记录不同通道接收到的超声回波信息作为参考的原始组织回波信号,记为并传回至数据存储器17。
具体地,首先需将阵列式超声换能器放置于待检组织上方,各换能器单元由多通道的电子线路控制发射超声脉冲波束和接收相应的回波信息。然后采用计算机控制信号发射器输出一个采控制信号至现场可编程门阵列(匹配电路)并分频得到两个时钟触发信号:一路时钟信号触发信号发生器产生一个短暂的正弦波电压信号,换能器受到电压激励,产生超声信号;在超声信号发射后,换能器阵列转为接收模式,对回波信号进行滤波、放大、时间增益补偿后,通过多路模数转换器以高于信号带宽四倍的采样频率进行同步采样,从而获得无混叠失真的高信噪比RF数据。另一路时钟触发信号控制信号发生器产生一个时延电压。其目的是为了预留时间让成像超声换能器采集施加声辐射力之前的射频回波信号,以便与施加声辐射力之后的射频信号进行比较并通过算法重构出仿体(或生物组织)的弹性特征。
在上述实验基础上,按照图2所示的工作序列进行实施,首先在***发射端,对待检组织发射非聚焦的平面波(plane wave,PW)超声波束41,同时向***接收端发射一个5-10V的脉冲触发信号,***接收端在收到触发信号后需经历短暂的硬件延时45,而后开始采集组织反射回来的原始超声回波信号48(采样深度5cm),并传回至数据采集卡内存(16GB DDR2 Memory Buffer)以此作为原始的参考超声射频回波信号,记为
②.通过发射端12控制阵列式超声换能器13对待检测组织14发射聚焦超声波束,在待检测组织内部产生声辐射力,声辐射力在声波传播方向上推动待检测组织,使待检测组织产生剪切波;
③.在声辐射力对待检测组织进行激励后,控制发射端12对待检测组织的待检测区域发射超声波束并接收反射回来的超声回波信号15。
④.重复步骤③,对待检测区域多次发射超声波束,通过数据收集模块16实时记录各通道接收到的超声射频回波信号并传回至数据存储器17、记为ψij,直到剪切波传播至待检测区域之外。
具体地,步骤③、④、⑤的实施方法:如图2所示,对超声换能器(或外部激励源如高强度聚焦超声换能器(High Intensity Focused Ultrasound,HIFU))施加10个周期的工作频率与换能器中心频率相同的正弦波电压信号42,并通过设定一定的延时序列控制超声换能器聚焦深度(如2.5cm),以此产生脉冲聚焦声辐射力;与此同时,通过两路钟信号触发信号发射器控制成像探头发射实时超声信号43,并由接收器接收实时超声射频信号44,脉冲时间间隔46为114微秒,采样深度为5cm,直至剪切波传播至组织待检测区域之外的区域47,并将接收到的实时超声射频回波信号44反馈至数字采集卡内存。
⑤.通过窗函数ωk 21对步骤④中实时接收的超声射频回波信号ψij 23和步骤①中记录原始超声射频回波信号22依次进行k次截取,k为自然数,每次截取沿横轴的长度为n(即窗函数取样长度为n),n为大于零的整数。n的数值,根据取样频率和图像分辨率确定,比如取样频率为1000M,分辨率为128×128,则如依序取样,则n值最小应为1000/128(取整)。如每次取样范围之间有交叉,则n值可比依据以上方法计算值大。优选取样范围有交叉,如图3所示,第一次取样A1到B1之间与第二次取样A2至B2之间有部分交叉。多次取样之间有交叉,可通过选择的窗函数(如汉宁窗等)来降低交叉部分的信号对结果造成的影响,以提高计算精度。如图3所示,第1次截取,利用汉宁窗函数ω1同时对原始超声射频回波信号22和实时超声射频回波信号23中A1至A2段的超声射频回波信号进行截取,按照取样频率,获得A1到B1之间的一个或多个超声回波信号(单位:电压V);第2次截取时,利用汉宁窗函数ω2对A2至B2段的超声射频回波信号进行截取;依次类推,直至第N次截取时,对AN至BN段的超声射频回波信号进行截取。每次截取沿横轴的长度相同,即A1到B1之间的长度与A2至B2之间的长度相同。而后分别计算每次截取后的多个超声射频回波信号与原始回波信号之间的模态置信度MAC值;计算公式如下;
其中,ωk(n)代表窗口长度为n的第k个窗函数,ψij代表在第i次发射超声波束后第j个通道接收到的超声射频回波信号(单位:电压V),代表t=0时刻下,第j个通道采集 的原始超声射频回波信号(单位:电压V),MACijk代表通过ωk(n)截取的超声信号与ψij及的模态置信度。
为了说明方便,示意图中只显示了以汉宁窗对超声射频回波信号进行截取的配置方案,本发明专利适用的窗函数也可以为矩形窗、海明窗、高斯窗、切比雪夫窗等,除非其他技术人员对相关方案的变化和改型偏离了本发明的范围,否则都应包括在本发明的保护范围内。
⑥.利用待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在该两位置上实时回波信号与原始回波信号的MAC值产生变化的时间差,计算该距离间剪切波传播的速度。
具体地,根据对待检测组织各位置对超声回波信号的MAC值,组成矩阵;并记录在不同时刻下待检测组织各位置对超声回波信号的MAC值,分别组成矩阵;依据待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在两位置上实时回波信号与原始回波信号的MAC值矩阵产生变化的时间差,计算该距离间剪切波传播的平均速度。
通过比较待检测组织内各位置的原始回波信号与实时获得的回波信号在不同时刻下的MAC值矩阵,分析得到不同时刻下待检测组织各位置上的由剪切波引起的组织位移情况,追踪剪切波在组织中的传播过程。参考图4,依据待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在两位置上实时回波信号与原始回波信号31的MAC值矩阵32产生变化的时间差,可以计算该距离间剪切波传播的平均速度。若这两个位置间的距离足够小,则可以近似地认为该平均速度为剪切波的瞬时传播速度,并以此构建组织的剪切波波速映射图33(图4所示)。
⑦.根据捕捉、记录沿待检测组织各方向传播的剪切波在不同位置上的瞬时波速,可以对待测组织的剪切波波速映射成像,进一步反推得到待测组织内各位置上的杨氏弹性模量绝对值,并在组织剪切波波速图像的基础上对待检测组织进行弹性成像,最终输出至显示器18。
利用计算机存储设备保存数据采集卡内存中的超声射频回波信号,并利用采集的原始超声射频回波信号48与采集的实时超声射频回波信号44计算得到了不同时刻下的模态置信度矩阵灰度映射图,该图像可以对剪切波在组织内的传播过程进行跟踪并计算剪切波传播速度。进而通过剪切波波速与组织的弹性信息的线性关系可以计算得到组织的弹性特性。重复上述步骤,最终可以达到对待检组织进行实时、定量弹性成像的技术目的。
以上所述的生物组织弹性成像方法,通过检测剪切波在生物组织中传输速度。利用剪切波与生物组织弹性的线性关系计算生物组织弹性,尤其适用于人体组织,特别是深层生物组 织。
本发明提供的基于模态置信准则获得剪切波弹性成像的方法,能对人体深层组织的弹性信息进行测量,具有实时、定量、快速、高分辨率等优点。其与现有超声弹性成像技术相比最大的不同在于本方法采用模态置信准则对受到剪切波影响的组织进行分析,以跟踪剪切波的传播过程。进而据此计算组织的弹性信息为目的进行组织弹性信息定量分析,通过比较正常组织与病变组织的弹性信息差异来实现组织病变的早期诊断。在具体实施例中,我们对各项同性的均匀超声体模进行了试验,检测深度为5cm至6cm,按照本专利阐述的方法对实际采集到的超声射频回波信号计算得到了不同时刻下的模态置信度矩阵映射图(图5),试验结果清晰地展示了由声辐射力激励产生的剪切波在仿体中的传播过程:声辐射力在0时刻下对组织进行激励,在2.645ms时刻下,待检测区域的模态置信度矩阵基本未产生变化,这表示待检测区域内组织并未发现剪切波。而后随着时间的推移,剪切波逐渐传播至检测区域内,第一边界线51、第二边界线52、第三边界线53和第四边界线54分别代表了在15.870ms、18.515ms、21.160ms、23.805ms时在剪切波作用下的待检测组织的形变边界线,时间间隔均为2.645ms。通过各边界线之间的距离及检测到各边界线的时间间隔可以计算得到不同时刻下剪切波在仿体内的平均传播速度。如18.515ms和21.160ms时刻下对应的剪切波传播位置为第三边界线52和第四边界线53,沿对角线方向的传播距离约为3.5mm,传播时间为2.645ms,进而得出期间的剪切波平均传播速度约为1.323m/s,符合理论阈值(1至10m/s)。本具体实施例旨在说明模态置信准则在剪切波弹性成像技术中的可实施性以及基于模态置信映射图计算剪切波传播速度的基本方法。若提高图像采集帧频,可以获得更加精确的数据。比如将图像采集帧频提高至3780KHz,即最小可以识别时间间隔为264.5μs的剪切波传播情况,小于实施例中展示的间隔时间2.645ms。根据以上的方法,可以对待检测组织采集大量的数据,并对大量数据进行分析、计算,能够获得更加精确的结果。
获得剪切波在该检测区域的传播速度后,再根据剪切波的传播速度与组织弹性的线性对应关系,可获知该位置的组织弹性。根据该组织的弹性数据,可对其健康状况作出判断。
本发明尤其适合于检测深层组织的弹性数据,距离表皮5cm、甚至10cm以上的组织,无法通过触诊方式接触到的人体组织器官,如肝脏等。
综上所述,本技术方法可有利地促进超声弹性成像,尤其是剪切波弹性成像的进一步推广和应用,加强人们对组织癌变的早期诊断能力。
本专利提供的基于模态置信准则获得剪切波弹性成像的方法,能对人体深层组织的弹性信息进行测量,具有实时、定量、快速、高分辨率等优点。其与现有超声弹性成像技术相比 最大的不同在于本方法采用模态置信准则对受到剪切波影响的组织进行分析,以跟踪剪切波的传播过程。进而据此计算组织的弹性信息为目的进行组织弹性信息定量分析,通过比较正常组织与病变组织的弹性信息差异来实现组织病变的早期诊断。
本发明中涉及的剪切波弹性成像技术有别于静态弹性成像技术的接触式激励方法,利用超声辐射力作为驱动力源,声波穿透组织并在深层位置聚焦产生作用力,该作用力非接触的作用于深层组织上使其产生在深层组织内传播的剪切波,进而通过跟踪该剪切波的传播过程对深层组织的弹性信息进行测量。因此,现有的静态弹性成像技术多应用于人体浅表组织的弹性信息测量(如***等),而本发明涉及的剪切波弹性成像方法可对深层组织的弹性信息进行测量(如肝脏等)。本发明中采用模态置信准则对剪切波进行跟踪,可以实现组织在剪切波作用下的微小形变检测,提高图像的信噪比、分辨率、对比度;灵敏度高。另一方面,模态置信准则是对超声回波射频信号的时域计算方法,计算流程一致性较好,可以允许多项连续的模态置信因子计算任务在相同时间的并发操作即并行计算。我们能够利用数据并行性使整幅弹性图像的复杂计算操作获得较高的加速比。因此,该方法可以有效的提高弹性成像算法的计算速度并实现剪切波的传播过程的实时跟踪,计算速度快。
本发明中的实施例仅用于对本发明进行说明,并不构成对权利要求范围的限制,本领域内技术人员可以想到的其他实质上等同的替代,均在本发明保护范围内。
Claims (11)
1.检测剪切波在生物组织中传输速度的方法,其特征在于,利用待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在该两位置上实时回波信号与原始回波信号的MAC值产生变化的时间差,计算该两位置间剪切波传播的速度。
2.根据权利要求1所述的检测剪切波在生物组织中传输速度的方法,其特征在于,
按照以下公式计算实时回波信号与原始回波信号之间的模态置信度MAC值:
其中,ωk(n)代表窗口长度为n的第k个取样窗,ψij代表在第i次发射波束后第j个通道接收到的回波信号,代表t=0时刻下,第j个通道采集的原始回波信号,MACijk代表通过ωk(n)截取的实时回波信号与ψij及的模态置信度。
3.根据权利要求1或2所述的检测剪切波在生物组织中传输速度的方法,其特征在于,根据对待检测组织各位置对超声回波信号的MAC值,组成矩阵;并记录在不同时刻下待检测组织各位置对超声回波信号的MAC值,分别组成矩阵;依据待检测组织内沿剪切波传播方向上的任意两位置的距离和分别在两位置上实时回波信号与原始回波信号的MAC值矩阵产生变化的时间差,计算该距离间剪切波传播的平均速度。
4.根据权利要求1所述的检测剪切波在生物组织中传输速度的方法,其特征在于,所述的实时回波信号为实时超声回波信号;所述原始回波信号为原始超声回波信号。
5.根据权利要求1所述的检测剪切波在生物组织中传输速度的方法,其特征在于,包括步骤:
①.待检测组织发射超声波束,并接收超声回波信息作为参考的原始组织回波信号,记为
②.对待检测组织发射聚焦波束,使待检测组织产生剪切波;
③.在待检测组织产生剪切波后,对待检测组织的待检测区域发射超声波束并接收反射回来的回波信号;
④.重复步骤③,对待检测区域多次发射超声波束,记录超声回波信号,记为ψij;直到剪切波传播至待检测区域之外;
⑤.通过窗函数ωk对步骤④中实时接收的超声回波信号ψij和步骤①中记录原始超声射频回波信号进行多次截取,计算MAC值。
6.根据权利要求5所述的检测剪切波在生物组织中传输速度的方法,其特征在于,所述的步骤②中,对待检测组织发射声波,在组织内部产生声辐射力,利用声辐射力在声波传播方向上推动待检测组织,使待检测组织产生剪切波。
7.根据权利要求5所述的检测剪切波在生物组织中传输速度的方法,其特征在于,
所述的步骤①中,通过超声仪器的发射端(12)控制阵列超声换能器(13)对待检测组织(14)发射超声波束,并通过接收端(15)和数据采集模块(16)记录不同通道接收到的超声回波信息作为参考的原始组织回波信号,记为并传回至数据存储器(17);
所述的步骤②中,通过发射端(12)控制阵列式超声换能器(13)对待检测组织(14)发射聚焦波束,在组织内部产生声辐射力,声辐射力在声波传播方向上推动待检测组织,使待检测组织产生剪切波;
所述的步骤③中,在声辐射力对待检测组织进行激励后,控制发射端(12)对待检测组织的待检测区域发射超声波束并接收反射回来的超声回波信号(15);
所述的步骤④中,重复步骤③,对待检测区域多次发射超声波束,通过数据采集模块(16)实时记录各通道接收到的超声射频回波信号并传回至数据存储器(17),记为ψij,直到剪切波传播至待检测区域之外;
所述的步骤⑤中,通过窗函数ωk(21)对步骤④中实时接收的超声射频回波信号ψij(23)和步骤①中记录原始超声射频回波信号(22)进行多次截取,并计算截取后的超声射频回波信号与原始回波信号之间的模态置信度MAC值。
8.根据权利要求1所述的检测剪切波在生物组织中传输速度的方法,其特征在于,所述生物组织为人体组织。
9.根据权利要求8所述的检测剪切波在生物组织中传输速度的方法,其特征在于,所述人体组织为深层组织。
10.检测生物组织弹性的方法,其特征在于,采用权利要求1至7任一权利要求所述的方法检测剪切波在生物组织中的传播速度,再利用剪切波在生物组织中的传播速度与弹性的线性关系计算得到组织的弹性特性数据。
11.生物组织弹性成像方法,其特征在于,利用权利要求1至7任一权利要求所述的方法计算生物组织的弹性,并对弹性数据生成映射图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410852544.0A CN104605891B (zh) | 2014-12-31 | 2014-12-31 | 检测剪切波在生物组织中传播速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410852544.0A CN104605891B (zh) | 2014-12-31 | 2014-12-31 | 检测剪切波在生物组织中传播速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104605891A true CN104605891A (zh) | 2015-05-13 |
CN104605891B CN104605891B (zh) | 2017-05-31 |
Family
ID=53140743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410852544.0A Active CN104605891B (zh) | 2014-12-31 | 2014-12-31 | 检测剪切波在生物组织中传播速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104605891B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107095692A (zh) * | 2016-02-19 | 2017-08-29 | 乐普(北京)医疗器械股份有限公司 | 超声波成像***、超声波成像方法及一维位移扫描方法 |
CN107530060A (zh) * | 2015-05-14 | 2018-01-02 | 安德拉有限公司 | 用于对生物组织结构进行成像的***和方法 |
WO2018082458A1 (zh) * | 2016-11-04 | 2018-05-11 | 声泰特(成都)科技有限公司 | 一种定量剪切波弹性成像方法及*** |
CN108065964A (zh) * | 2018-01-16 | 2018-05-25 | 中国科学院苏州生物医学工程技术研究所 | 一种超声成像方法、装置、设备及超声成像探头 |
CN109828030A (zh) * | 2019-03-28 | 2019-05-31 | 深圳中凯剑无损检测设备科技有限公司 | 一种基于声场特征的反射体形貌提取***和提取方法 |
CN110151270A (zh) * | 2019-05-10 | 2019-08-23 | 中国科学院苏州生物医学工程技术研究所 | 振动穿刺装置、频率调整方法及超声成像*** |
CN110418609A (zh) * | 2017-10-19 | 2019-11-05 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及弹性对比测量方法 |
CN110946618A (zh) * | 2019-11-24 | 2020-04-03 | 华南理工大学 | 一种基于经验流形的弹性成像方法 |
CN111050661A (zh) * | 2018-05-15 | 2020-04-21 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性测量方法及剪切波弹性成像*** |
CN111281435A (zh) * | 2019-12-24 | 2020-06-16 | 徐浩 | 基于水平剪切波的材料弹性测量方法及材料弹性测量*** |
WO2020140917A1 (zh) * | 2019-01-02 | 2020-07-09 | 无锡海斯凯尔医学技术有限公司 | 用于检测组织硬度的方法、设备及*** |
CN111449629A (zh) * | 2020-04-28 | 2020-07-28 | 北京信息科技大学 | 一种光学相干弹性成像方法及装置 |
CN112702953A (zh) * | 2018-09-06 | 2021-04-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性成像方法和装置、及计算机存储介质 |
CN112731359A (zh) * | 2020-12-31 | 2021-04-30 | 无锡祥生医疗科技股份有限公司 | 超声探头的速度确定方法、装置及存储介质 |
CN113116381A (zh) * | 2019-12-30 | 2021-07-16 | 重庆融海超声医学工程研究中心有限公司 | 生物组织的性质变化检测方法、变性控制方法及对应装置 |
CN114287968A (zh) * | 2019-09-27 | 2022-04-08 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法、***及计算机可读存储介质 |
US11304678B2 (en) * | 2017-01-10 | 2022-04-19 | Koninklijke Philips N.V. | Systems, methods, and apparatuses for confidence mapping of shear wave imaging |
CN114748037A (zh) * | 2022-03-14 | 2022-07-15 | 江勇 | 一种人工植入肝脏弹性检测装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
CN1586410A (zh) * | 2004-08-06 | 2005-03-02 | 清华大学 | 一种变尺度的生物组织位移估计方法 |
CN1586408A (zh) * | 2004-08-20 | 2005-03-02 | 清华大学 | 一种多尺度的生物组织位移估计方法 |
US20100016718A1 (en) * | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
CN101675356A (zh) * | 2007-05-08 | 2010-03-17 | 柏林夏瑞蒂医学院 | 用于身体组织的弹性成像检查的方法和装置 |
WO2013026141A1 (en) * | 2011-08-19 | 2013-02-28 | The University Of British Columbia | Elastography using ultrasound imaging of a thin volume |
-
2014
- 2014-12-31 CN CN201410852544.0A patent/CN104605891B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
CN1586410A (zh) * | 2004-08-06 | 2005-03-02 | 清华大学 | 一种变尺度的生物组织位移估计方法 |
CN1586408A (zh) * | 2004-08-20 | 2005-03-02 | 清华大学 | 一种多尺度的生物组织位移估计方法 |
CN101675356A (zh) * | 2007-05-08 | 2010-03-17 | 柏林夏瑞蒂医学院 | 用于身体组织的弹性成像检查的方法和装置 |
US20100016718A1 (en) * | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
WO2013026141A1 (en) * | 2011-08-19 | 2013-02-28 | The University Of British Columbia | Elastography using ultrasound imaging of a thin volume |
Non-Patent Citations (1)
Title |
---|
MIROSLAV PASTOR,MICHAL BINDA,TOMAS HARCARIK: "Modal Assurance Criterion", 《PROCEDIA ENGINEERING》 * |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107530060A (zh) * | 2015-05-14 | 2018-01-02 | 安德拉有限公司 | 用于对生物组织结构进行成像的***和方法 |
CN107095692A (zh) * | 2016-02-19 | 2017-08-29 | 乐普(北京)医疗器械股份有限公司 | 超声波成像***、超声波成像方法及一维位移扫描方法 |
CN107095692B (zh) * | 2016-02-19 | 2024-02-09 | 乐普(北京)医疗器械股份有限公司 | 超声波成像***、超声波成像方法及一维位移扫描方法 |
WO2018082458A1 (zh) * | 2016-11-04 | 2018-05-11 | 声泰特(成都)科技有限公司 | 一种定量剪切波弹性成像方法及*** |
US11357480B2 (en) | 2016-11-04 | 2022-06-14 | Saset Chengdu Technology Ltd. | Quantitative shear wave elasticity imaging method and system |
US11304678B2 (en) * | 2017-01-10 | 2022-04-19 | Koninklijke Philips N.V. | Systems, methods, and apparatuses for confidence mapping of shear wave imaging |
CN110418609B (zh) * | 2017-10-19 | 2021-04-20 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及弹性对比测量方法 |
CN110418609A (zh) * | 2017-10-19 | 2019-11-05 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及弹性对比测量方法 |
CN108065964A (zh) * | 2018-01-16 | 2018-05-25 | 中国科学院苏州生物医学工程技术研究所 | 一种超声成像方法、装置、设备及超声成像探头 |
CN108065964B (zh) * | 2018-01-16 | 2021-04-20 | 中国科学院苏州生物医学工程技术研究所 | 一种超声成像方法、装置、设备及超声成像探头 |
CN111050661A (zh) * | 2018-05-15 | 2020-04-21 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性测量方法及剪切波弹性成像*** |
CN111050661B (zh) * | 2018-05-15 | 2023-08-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性测量方法及剪切波弹性成像*** |
US11998392B2 (en) | 2018-05-15 | 2024-06-04 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Shear wave elasticity measurement method and shear wave elasticity imaging system |
CN112702953A (zh) * | 2018-09-06 | 2021-04-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性成像方法和装置、及计算机存储介质 |
CN111388012A (zh) * | 2019-01-02 | 2020-07-10 | 无锡海斯凯尔医学技术有限公司 | 用于检测组织硬度的方法、设备及*** |
WO2020140917A1 (zh) * | 2019-01-02 | 2020-07-09 | 无锡海斯凯尔医学技术有限公司 | 用于检测组织硬度的方法、设备及*** |
CN111388012B (zh) * | 2019-01-02 | 2021-10-08 | 无锡海斯凯尔医学技术有限公司 | 用于检测组织硬度的方法、设备及*** |
CN109828030A (zh) * | 2019-03-28 | 2019-05-31 | 深圳中凯剑无损检测设备科技有限公司 | 一种基于声场特征的反射体形貌提取***和提取方法 |
CN110151270A (zh) * | 2019-05-10 | 2019-08-23 | 中国科学院苏州生物医学工程技术研究所 | 振动穿刺装置、频率调整方法及超声成像*** |
CN110151270B (zh) * | 2019-05-10 | 2024-01-30 | 中国科学院苏州生物医学工程技术研究所 | 振动穿刺装置、频率调整方法及超声成像*** |
CN114287968A (zh) * | 2019-09-27 | 2022-04-08 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法、***及计算机可读存储介质 |
CN110946618A (zh) * | 2019-11-24 | 2020-04-03 | 华南理工大学 | 一种基于经验流形的弹性成像方法 |
CN110946618B (zh) * | 2019-11-24 | 2022-03-29 | 华南理工大学 | 一种基于经验流形的弹性成像方法 |
CN111281435A (zh) * | 2019-12-24 | 2020-06-16 | 徐浩 | 基于水平剪切波的材料弹性测量方法及材料弹性测量*** |
CN113116381A (zh) * | 2019-12-30 | 2021-07-16 | 重庆融海超声医学工程研究中心有限公司 | 生物组织的性质变化检测方法、变性控制方法及对应装置 |
CN111449629A (zh) * | 2020-04-28 | 2020-07-28 | 北京信息科技大学 | 一种光学相干弹性成像方法及装置 |
CN112731359A (zh) * | 2020-12-31 | 2021-04-30 | 无锡祥生医疗科技股份有限公司 | 超声探头的速度确定方法、装置及存储介质 |
CN112731359B (zh) * | 2020-12-31 | 2024-04-09 | 无锡祥生医疗科技股份有限公司 | 超声探头的速度确定方法、装置及存储介质 |
CN114748037B (zh) * | 2022-03-14 | 2022-11-29 | 江勇 | 一种人工植入肝脏弹性检测装置 |
CN114748037A (zh) * | 2022-03-14 | 2022-07-15 | 江勇 | 一种人工植入肝脏弹性检测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104605891B (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104605891A (zh) | 检测剪切波在生物组织中传输速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法 | |
CN102283679B (zh) | 弹性测量的超声成像***及测量生物组织弹性的方法 | |
CA2068740C (en) | Method and apparatus for measurement and imaging of tissue compressibility or compliance | |
CN103269639B (zh) | 利用质心估计剪切波速度 | |
US9492139B2 (en) | Non-imaging low frequency ultrasonic testing and diagnostic evaluation system | |
JP6140953B2 (ja) | 医用超音波剪断波イメージングにおける分類前処理 | |
US10795007B2 (en) | Method and a device for imaging a visco-elastic medium | |
CN105816205B (zh) | 声辐射力脉冲成像中的稀疏跟踪 | |
US9585631B2 (en) | Devices, methods, and systems for measuring elastic properties of biological tissues using acoustic force | |
CN107613878A (zh) | 用于检测物体运动的超声成像***和方法 | |
JP2010526626A (ja) | 関心領域の粘弾性の平均値を測定するための方法および装置 | |
RU2641067C1 (ru) | Ультразвуковая система и способ измерения с использованием сдвиговой волны | |
CN109717899A (zh) | 超声医学成像中根据剪切速率的组织粘弹性估计 | |
EP2073713B1 (en) | Method and apparatus for acoustoelastic extraction of strain and material properties | |
Van Pamel et al. | Numerical simulations of ultrasonic array imaging of highly scattering materials | |
JP2014512923A (ja) | オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 | |
JP2017217147A (ja) | 超音波診断装置及び超音波信号処理方法 | |
KR20140086626A (ko) | 전단파의 변위 산출 방법, 전단파를 이용한 피검체의 기계적 계수 산출 방법 및 이를 포함하는 시스템 | |
CN108852416B (zh) | 一种剪切波传播速度的确定方法及装置 | |
Bochud et al. | Probabilistic inverse problem to characterize tissue-equivalent material mechanical properties | |
JP7354632B2 (ja) | 超音波診断装置、および、超音波診断装置の制御方法 | |
CN108784744A (zh) | 用于剪切波成像的可变焦点 | |
CN102046092A (zh) | 组织应变分析 | |
Sultan et al. | Cortical Bone Thickness Assessment from Multi-frequency Ultrasound RF Data using a Convolutional Architecture with Multi-head Attention | |
JP7302651B2 (ja) | 超音波信号処理装置、超音波診断装置、超音波信号処理方法、およびプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |