CN103040488A - 一种实时超声弹性成像位移估计方法和*** - Google Patents

一种实时超声弹性成像位移估计方法和*** Download PDF

Info

Publication number
CN103040488A
CN103040488A CN2012105612359A CN201210561235A CN103040488A CN 103040488 A CN103040488 A CN 103040488A CN 2012105612359 A CN2012105612359 A CN 2012105612359A CN 201210561235 A CN201210561235 A CN 201210561235A CN 103040488 A CN103040488 A CN 103040488A
Authority
CN
China
Prior art keywords
displacement
real
displacement field
field
algorithm
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
Application number
CN2012105612359A
Other languages
English (en)
Other versions
CN103040488B (zh
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.)
Shenzhen University
Original Assignee
Shenzhen 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 Shenzhen University filed Critical Shenzhen University
Priority to CN201210561235.9A priority Critical patent/CN103040488B/zh
Publication of CN103040488A publication Critical patent/CN103040488A/zh
Application granted granted Critical
Publication of CN103040488B publication Critical patent/CN103040488B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种实时超声弹性成像位移估计方法和***，其方法首先采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场
Figure 783871DEST_PATH_IMAGE001
，所述两帧射频信号之间存在位移变化；然后设定一先验估计tp，所述先验估计tp等于位移场
Figure 701011DEST_PATH_IMAGE001
的轴向位移
Figure 2012105612359100004DEST_PATH_IMAGE002
；使得两帧射频信号的位移保持1/4个波长之内最后，根据先验估计tp和位移场
Figure 543065DEST_PATH_IMAGE001
的侧向位移
Figure 938275DEST_PATH_IMAGE003
，计算相移，得到轴向位移场
Figure DEST_PATH_IMAGE004
。其解决了传统相移算法中容易出现的相位混叠问题，在兼顾运算速度的前提下，有效地提高了估计的精度。

Description

一种实时超声弹性成像位移估计方法和***
技术领域
本发明涉及超声弹性成像技术领域，尤其涉及一种实时超声弹性成像位移估计方法和***。 
背景技术
软组织的病变常常伴随着组织硬度的改变。比如，乳腺或者***实质性肿瘤的硬度要大于正常组织，而液性囊肿的硬度小于正常组织。传统上，临床医生用触诊来定性判断肿块的软硬程度。然而，当肿块埋藏于较深的软组织内部，或者体积较小时，通过触诊就很难发现这些组织病变。超声弹性成像是近年来发展比较快的一种新的超声成像技术，可以对组织力学属性（如组织弹性、组织硬度）进行成像，为临床疾病诊断提供更加丰富的功能信息。
然而，现有技术的超声弹性成像技术存在一定的不足：
一、现有算法在速度和精确度方面存在着相互制约的问题。目前绝大部分的实时性好的算法都是一维的（即只考虑组织在受压时的轴向运动），其中包括题目为：“Real time tissue elasticity imaging using the combined autocorrelation method”，“Motion estimation in ultrasound images using time domain cross correlation with prior estimates”等文章。这些算法假设在同一根射频（Radio Frequency, RF）信号上组织位移连续，因此可以利用前一个位移的估计值作为下一次估计的先验值，从而减小互相关的搜索范围。这些一维算法在结构相对简单的体模上进行实验时可以获得很好的结果。由于对人体组织施压后，其产生位移场比较复杂，当这些算法用于活体实验时很难取得满意的结果。
二、二维位移估计算法的鲁棒性较好，但传统二维估计算法的缺点是运算量过大，无法满足临床实时成像的需求。近些年来，一些研究组提出的新算法，比如动态规划算法“Ultrasound elastography: a dynamic programming approach ”、光流法“A motion estimation refinement framework for real-time tissue axial strain estimation with freehand ultrasound” 以及快速互相关算法“A fast normalized cross-correlation method for motion estimation”等，这些算法同时考虑了组织受压产生的侧向位移，鲁棒性比较好，同时具有较好的实时性，但是其在处理组织大位移方面，容易出现误差。
有鉴于此，现有技术有待改进和提高。
发明内容
鉴于现有技术的不足，本发明目的在于提供一种实时超声弹性成像位移估计方法和***。旨在解决现有技术中针对复杂的位移场容易出现误差等问题。
本发明的技术方案如下：
一种实时超声弹性成像位移估计方法，其中，包括以下步骤：
S1、采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场                                                ，所述两帧射频信号之间存在位移变化；
S2、设定一先验估计tp，所述先验估计tp等于位移场
Figure 671409DEST_PATH_IMAGE001
的轴向位移
Figure 2012105612359100002DEST_PATH_IMAGE002
；使得经过预估值修正之后的两帧射频信号的位移保持1/4个波长之内；
S3、根据先验估计tp和位移场
Figure 556188DEST_PATH_IMAGE001
的侧向位移
Figure 2012105612359100002DEST_PATH_IMAGE003
，计算相移，得到轴向位移场
Figure 2012105612359100002DEST_PATH_IMAGE004
；
其中，u、v 分别代表轴向和侧向位移量。
优选地，所述的实时超声弹性成像位移估计方法，其中，所述步骤S1中，二维匹配算法分为二维相关检测算法和绝对差分求和算法。
优选地，所述的实时超声弹性成像位移估计方法，其中，所述二维匹配算法包括：最小均方误差函数、最小平均绝对值函数和/或最大匹配像素统计。
优选地，所述的实时超声弹性成像位移估计方法，其中，所述步骤S1具体包括以下步骤：
S11、选择两帧射频信号，以第一射频信号
Figure 2012105612359100002DEST_PATH_IMAGE005
中第一区域作为目标区域，在第二射频信号
Figure 2012105612359100002DEST_PATH_IMAGE006
寻找最佳匹配的第二区域，其中，相关系数定义为
Figure 2012105612359100002DEST_PATH_IMAGE007
；
其中，
Figure 2012105612359100002DEST_PATH_IMAGE008
，
Figure 2012105612359100002DEST_PATH_IMAGE009
、
Figure 2012105612359100002DEST_PATH_IMAGE010
为对应图像内的平均像素；
S12、计算第一区域和第二区域的绝对差分和，求出位移场
Figure 2012105612359100002DEST_PATH_IMAGE011
；
其中，
Figure 2012105612359100002DEST_PATH_IMAGE012
， u、v 分别代表轴向和侧向位移量；m、 n分别是搜索区域的轴向与横向尺寸。
优选地，所述的实时超声弹性成像位移估计方法，其中，所述步骤S3具体包括：
S31、令
Figure 2012105612359100002DEST_PATH_IMAGE013
，
Figure 2012105612359100002DEST_PATH_IMAGE014
，
代入
Figure 2012105612359100002DEST_PATH_IMAGE015
，得到相移
Figure 2012105612359100002DEST_PATH_IMAGE016
；
其中，t0 是预先设定的相关窗的长度， *表示复数共轭；
S32、将相移
Figure 573954DEST_PATH_IMAGE016
代入
Figure 2012105612359100002DEST_PATH_IMAGE017
，得到轴向位移场
Figure 2012105612359100002DEST_PATH_IMAGE018
；
其中，c 为组织内的超声传播速度，w0 是超声探头的中心角频率。
一种实时超声弹性成像位移估计***，其中，包括：
预估模块，用于采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场，所述两帧射频信号之间存在位移变化；
先验估计模块，用于设定一先验估计tp，所述先验估计tp等于位移场
Figure 348355DEST_PATH_IMAGE019
的轴向位移
Figure 2012105612359100002DEST_PATH_IMAGE020
；使得经过预估值修正之后的两帧射频信号的位移保持1/4个波长之内；
位移计算模块，用于根据先验估计tp和位移场
Figure 17234DEST_PATH_IMAGE019
的侧向位移
Figure DEST_PATH_IMAGE021
，计算相移，得到轴向位移场
Figure 72914DEST_PATH_IMAGE018
。
有益效果：
本申请的实时超声弹性成像位移估计方法和***，解决了传统相移算法中容易出现的相位混叠问题，在兼顾运算速度的前提下，有效地提高了估计的精度。 
附图说明
图1为本申请的实时超声弹性成像位移估计方法的流程图。
图2为本申请的实时超声弹性成像位移估计***的结构框图。
图 3a、图3c分别为应变0.5%和3%时，仿真的射频信号用本申请的实时超声弹性成像位移估计方法得到的位移图，图3b、图3d分别为应变为0.5%和3%时有限元软件输出的真实位移图。
图 4a、图4b分别为应变0.5%和3%时，图3b、图3d中实线标出部分的位移值。其中实线为本申请的实时超声弹性成像位移估计方法得到的位移值，虚线为有限元软件输出的真实位移值。
图 5a、图5c分别为应变0.5%和3%时，仿真的射频信号用本申请的实时超声弹性成像位移估计方法得到的应变图，图5b、图5d分别为应变为0.5%和3%时有限软软件输出的真实应变图。
图6 a为从弹性体模获得的应变分布图。
图6b为从肝癌病人经过射频消融之后超声获取的应变分布图。
具体实施方式
本发明提供一种实时超声弹性成像位移估计方法和***，为使本发明的目的、技术方案及效果更加清楚、明确，以下对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，其为本申请的实时超声弹性成像位移估计方法的流程图。如图所示，所述实时超声弹性成像位移估计方法包括以下步骤：
S1、采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场
Figure 889561DEST_PATH_IMAGE019
，所述两帧射频信号之间存在位移变化；
S2、设定一先验估计tp，所述先验估计tp等于位移场
Figure 626572DEST_PATH_IMAGE019
的轴向位移
Figure 149958DEST_PATH_IMAGE020
；使得两帧射频信号的位移保持1/4个波长之内；
S3、根据先验估计tp和位移场
Figure 579802DEST_PATH_IMAGE019
的侧向位移
Figure 883744DEST_PATH_IMAGE021
，计算相移，得到轴向位移场
Figure 158868DEST_PATH_IMAGE018
。
下面分别针对上述步骤进行详细介绍：
所述步骤S1为采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场，所述两帧射频信号之间存在位移变化（如采用压缩等方式）。其中，所述二维快匹配算法通过跟踪特定散射子在不同视频图像之间的变化来得到组织的运动信息。根据匹配算法的不同可分为二维相关检测算法和绝对差分求和算法（SAD）。
在本实施例中，所述步骤S1具体包括以下步骤：
S11、选择两帧射频信号，以第一射频信号
Figure 137505DEST_PATH_IMAGE005
中第一区域作为目标区域，在第二射频信号
Figure 2012105612359100002DEST_PATH_IMAGE022
寻找最佳匹配的第二区域，其中，相关系数定义为
Figure 430209DEST_PATH_IMAGE007
；
其中，
Figure DEST_PATH_IMAGE023
，
Figure 2012105612359100002DEST_PATH_IMAGE024
、为对应图像内的平均像素；对应于相关系数最大的位置就是第二区域。基本相关算法需要在全局范围内计算两帧图像的互相关系数，运算量很大。人们通过灵活选择匹配准则和搜索方法来降低相关算法的运算量。可供选择的其他匹配准则主要有，最小均方误差函数（MSE），最小平均绝对值函数（MAD），最大匹配像素统计（MPC）等。在检测目标速度较快或者实时性要求较高的场合，二维相关搜索算法在计算方面达不到要求。故此提出不计算计算量较大的相关函数，而是计算两个检测区的绝对差分和即SAD。即
S12、计算第一区域和第二区域的绝对差分和，求出位移场
；
其中，
Figure 571340DEST_PATH_IMAGE023
， u、v 分别代表轴向和侧向位移量；m、 n分别是搜索区域的轴向与横向尺寸。根据上面的计算公式可知，二维块匹配算法的主要缺点就是只能得到精确到采样点间隔的位移分布，而对于亚采样点间隔的位移只能通过插值得到，而不同的插值算法对于位移估计的精确程度有较大影响。
所述步骤S2为设定一先验估计tp，所述先验估计tp等于位移场
Figure 131634DEST_PATH_IMAGE019
的轴向位移
Figure 903281DEST_PATH_IMAGE020
；使得两帧射频信号的位移保持1/4个波长之内。在传统相移算法中，有一个比较明显的缺点就是位移是根据相移来计算的。而相位信息最大的范围是
Figure DEST_PATH_IMAGE027
，所以当位移过大时容易产生相位混叠。故此，本发明提出了一个比较简单的方法，即加入先验估计tp.。使得两帧射频信号的位移保持在较小的范围之内（即1/4个波长）。
所述步骤S3为根据先验估计tp和位移场
Figure 853920DEST_PATH_IMAGE019
的侧向位移
Figure 798742DEST_PATH_IMAGE021
，计算相移，得到轴向位移场
Figure 151226DEST_PATH_IMAGE018
。在本实施例中，步骤S3具体包括：
S31、令，
Figure 531709DEST_PATH_IMAGE014
，
带入
Figure 2012105612359100002DEST_PATH_IMAGE028
，得到相移
Figure 516107DEST_PATH_IMAGE016
；
其中，t0 是预先设定的相关窗的长度， *表示复数共轭；
S32、将相移
Figure 988677DEST_PATH_IMAGE016
带入
Figure 102126DEST_PATH_IMAGE017
，得到轴向位移场；
其中，c 为声速，w0 是探头的中心角频率。
进一步地，在求解tp的过程中，根据不同的处理方式也得到了不同的位移估计方法。比如“Real time tissue elasticity imaging using the combined autocorrelation method”,以及“A strain-based ultrasound elastography using phase shift with prior estimates and meshfree shape function”。
概括来说，本申请的实时超声弹性成像位移估计方法的实施例中，首先通过SAD得到精确到整数间隔的位移估计
Figure 313982DEST_PATH_IMAGE019
，第一步得到的位移场，作为第二步相移算法的先验估计。根据侧向位移
Figure 659830DEST_PATH_IMAGE021
，得到RF信号压缩后的侧向位置，避免了传统一维算法在单根线上进行追踪导致的错误估计。而轴向位移分量
Figure 134673DEST_PATH_IMAGE020
作为轴向位置的先验估计，避免了传统相移算法出现的相位混叠问题。
与此对应，本发明还提供了一种实时超声弹性成像位移估计***，如图2所示，其包括：
预估模块100，用于采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场
Figure 162672DEST_PATH_IMAGE019
，所述两帧射频信号之间存在位移变化；
先验估计模块200，用于设定一先验估计tp，所述先验估计tp等于位移场的轴向位移
Figure 596245DEST_PATH_IMAGE020
；使得两帧射频信号的位移保持1/4个波长之内；
位移计算模块300，用于根据先验估计tp和位移场
Figure 496068DEST_PATH_IMAGE019
的侧向位移
Figure 625960DEST_PATH_IMAGE021
，计算相移，得到轴向位移场
Figure 927628DEST_PATH_IMAGE018
。
为了验证该方法的有效性，我们用仿真的射频信号，弹性体模和病人数据进行了实验，结果表明，该算法对应变较大且侧向位移比较复杂的位移场也具有较好的鲁棒性。图 3a、图3c分别为应变0.5%和3%时，仿真的射频信号用本申请的实时超声弹性成像位移估计方法得到的位移图，图3b、图3d分别为应变为0.5%和3%时有限元软件输出的真实位移图。图 4a、图4b分别为应变0.5%和3%时，图3b、图3d中实线标出部分的位移值。其中实线为本申请的实时超声弹性成像位移估计方法得到的位移值，虚线为有限元软件输出的真实位移值。其中灰度表示向组织施加同一外力时，受压组织内部的位移大小，其中浅色表示位移小的区域，深色表示位移大的区域。概括来说，图3a、图3c以及图4a、图4c分别显示了当应变为0.5%，3%时，仿真的射频信号用本算法得到的位移图以及应变图。与图3b、图3d和4a、图4c的有限元仿真结果相比可知，对于应变较小（0.5%）或者较大（3%）时，该算法都可以得到比较精确的位移估计结果。即本方案的估计结果与真实值吻合度非常高：在组织大变形的情况下算法仍有很高的估算准确性。
请继续参阅图5a、图5b、图5c和图5d，其中，图 5a、图5c分别为应变0.5%和3%时，仿真的射频信号用本申请的实时超声弹性成像位移估计方法得到的应变图，图5b、图5d分别为应变为0.5%和3%时有限软软件输出的真实应变图。其中灰度表示向组织施加同一外力时，受压组织内部的应变大小，应变为组织受压后的变形与其原始的尺寸之比。其中浅色表示应变大的区域，深色表示应变小的区域。在同一外力作用下，应变大的区域硬度较小，反之，应变小的区域硬度较大。因此，超声弹性成像中用应变大小表征组织的硬度。
同时，图6a和图6b给出了该方法用于弹性体模以及病人数据时的结果。其中图6a中的白色圆圈内暗色区域表示硬度较大的区域，图6b中白色箭头所指暗色区域表示经消融后硬化的肝组织区域。对于体模数据本文提出的算法也可以计算得到比较平滑的应变场。对于位移分布极其复杂的病人数据，该方法同样可以得到理想的结果，弹性图像可以清晰显示经消融后硬化的肝组织区域。
综上所述，本发明的实时超声弹性成像位移估计方法和***，其方法首先采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场
Figure 288202DEST_PATH_IMAGE019
，所述两帧射频信号之间存在位移变化；然后设定一先验估计tp，所述先验估计tp等于位移场
Figure 675321DEST_PATH_IMAGE019
的轴向位移
Figure DEST_PATH_IMAGE029
；使得两帧射频信号的位移保持1/4个波长之内最后，根据先验估计tp和位移场的侧向位移
Figure 998035DEST_PATH_IMAGE021
，计算相移，得到轴向位移场
Figure 795090DEST_PATH_IMAGE018
。其解决了传统相移算法中容易出现的相位混叠问题，在兼顾运算速度的前提下，有效地提高了估计的精度。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。 

Claims (6)

1.一种实时超声弹性成像位移估计方法，其特征在于，包括以下步骤：
S1、采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场                                                
Figure 2012105612359100001DEST_PATH_IMAGE001
，所述两帧射频信号之间存在位移变化；
S2、设定一先验估计tp，所述先验估计tp等于位移场
Figure 2012105612359100001DEST_PATH_IMAGE002
的轴向位移
Figure 2012105612359100001DEST_PATH_IMAGE003
；使得经过预估值修正之后的两帧射频信号的位移保持1/4个波长之内；
S3、根据先验估计tp和位移场
Figure 947566DEST_PATH_IMAGE002
的侧向位移
Figure 2012105612359100001DEST_PATH_IMAGE004
，计算相移，得到轴向位移场
Figure 2012105612359100001DEST_PATH_IMAGE005
；
其中，u、v 分别代表轴向和侧向位移量。
2.根据权利要求1所述的实时超声弹性成像位移估计方法，其特征在于，所述步骤S1中，二维匹配算法分为二维相关检测算法和绝对差分求和算法。
3.根据权利要求2所述的实时超声弹性成像位移估计方法，其特征在于，所述二维匹配算法包括：最小均方误差函数、最小平均绝对值函数和/或最大匹配像素统计。
4.根据权利要求1所述的实时超声弹性成像位移估计方法，其特征在于，所述步骤S1具体包括以下步骤：
S11、选择两帧射频信号，以第一射频信号中第一区域作为目标区域，在第二射频信号
Figure DEST_PATH_IMAGE007
寻找最佳匹配的第二区域，其中，相关系数定义为
Figure 2012105612359100001DEST_PATH_IMAGE008
；
其中，，
Figure 2012105612359100001DEST_PATH_IMAGE010
，
Figure DEST_PATH_IMAGE011
为对应图像内的平均像素；
S12、计算第一区域和第二区域的绝对差分和，求出位移场
；
其中，
Figure 308009DEST_PATH_IMAGE009
， u、v 分别代表轴向和侧向位移量；m、 n分别是搜索区域的轴向与横向尺寸。
5.根据权利要求1所述的实时超声弹性成像位移估计方法，其特征在于，所述步骤S3具体包括：
S31、令
Figure DEST_PATH_IMAGE013
，
Figure 2012105612359100001DEST_PATH_IMAGE014
，
代入
Figure DEST_PATH_IMAGE015
，得到相移
Figure 2012105612359100001DEST_PATH_IMAGE016
；
其中，t0 是预先设定的相关窗的长度， *表示复数共轭；
S32、将相移
Figure 882472DEST_PATH_IMAGE016
代入
Figure DEST_PATH_IMAGE017
，得到轴向位移场；
其中，c 为组织内的超声传播速度，w0 是超声探头的中心角频率。
6.一种实时超声弹性成像位移估计***，其特征在于，包括：
预估模块，用于采用二维匹配算法对待处理的两帧射频信号进行预估，得到位移场
Figure 8877DEST_PATH_IMAGE002
，所述两帧射频信号之间存在位移变化；
先验估计模块，用于设定一先验估计tp，所述先验估计tp等于位移场
Figure 720481DEST_PATH_IMAGE002
的轴向位移
Figure 647986DEST_PATH_IMAGE003
；使得经过预估值修正之后的两帧射频信号的位移保持1/4个波长之内；
位移计算模块，用于根据先验估计tp和位移场
Figure 786843DEST_PATH_IMAGE002
的侧向位移，计算相移，得到轴向位移场
Figure 619987DEST_PATH_IMAGE005
。
CN201210561235.9A 2012-12-21 2012-12-21 一种实时超声弹性成像位移估计方法和*** Expired - Fee Related CN103040488B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210561235.9A CN103040488B (zh) 2012-12-21 2012-12-21 一种实时超声弹性成像位移估计方法和***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210561235.9A CN103040488B (zh) 2012-12-21 2012-12-21 一种实时超声弹性成像位移估计方法和***

Publications (2)

Publication Number Publication Date
CN103040488A true CN103040488A (zh) 2013-04-17
CN103040488B CN103040488B (zh) 2014-06-04

Family

ID=48053507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210561235.9A Expired - Fee Related CN103040488B (zh) 2012-12-21 2012-12-21 一种实时超声弹性成像位移估计方法和***

Country Status (1)

Country Link
CN (1) CN103040488B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103549978A (zh) * 2013-11-08 2014-02-05 中国科学院声学研究所 一种瞬时弹性成像快速时移估计方法
CN103654865A (zh) * 2013-12-26 2014-03-26 华南理工大学 基于最大互信息的超声弹性成像组织位移估计方法
CN103735287A (zh) * 2013-12-05 2014-04-23 中国科学院苏州生物医学工程技术研究所 一种血管内超声弹性成像二维多级混合位移估计方法
CN105232087A (zh) * 2015-11-05 2016-01-13 无锡祥生医学影像有限责任公司 超声弹性成像实时处理***
CN105266849A (zh) * 2014-07-09 2016-01-27 无锡祥生医学影像有限责任公司 实时超声弹性成像方法和***
CN109640830A (zh) * 2016-07-14 2019-04-16 医视特有限公司 基于先例的超声聚焦
CN110215233A (zh) * 2019-04-30 2019-09-10 深圳大学 一种基于超声平面波扫描的分段式脉搏波成像方法
CN113925598A (zh) * 2020-07-13 2022-01-14 通用电气精准医疗有限责任公司 利用超声成像***对组织进行热监测的方法和***
CN115063599A (zh) * 2022-06-16 2022-09-16 昆明理工大学 一种应用于中小型水库大坝监测的小波光流估计与图像相关变形识别方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0312140A (ja) * 1989-06-08 1991-01-21 Yokogawa Medical Syst Ltd 超音波診断装置
JPH04164440A (ja) * 1990-10-29 1992-06-10 Noritoshi Nakabachi 超音波による心臓各部位の振動の計測方法
CN1586411A (zh) * 2004-08-06 2005-03-02 清华大学 一种二维综合互相关的生物组织位移估计方法
CN101569543A (zh) * 2008-04-29 2009-11-04 香港理工大学 弹性成像的二维位移估计方法
CN102048560A (zh) * 2010-12-14 2011-05-11 深圳市蓝韵实业有限公司 一种采用双尺度的生物组织位移估计方法
CN102764141A (zh) * 2012-07-20 2012-11-07 中国科学院深圳先进技术研究院 弹性成像方法和***及其中的生物组织位移估计方法和***
CN102824193A (zh) * 2011-06-14 2012-12-19 深圳迈瑞生物医疗电子股份有限公司 一种弹性成像中的位移检测方法、装置及***

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0312140A (ja) * 1989-06-08 1991-01-21 Yokogawa Medical Syst Ltd 超音波診断装置
JPH04164440A (ja) * 1990-10-29 1992-06-10 Noritoshi Nakabachi 超音波による心臓各部位の振動の計測方法
CN1586411A (zh) * 2004-08-06 2005-03-02 清华大学 一种二维综合互相关的生物组织位移估计方法
CN101569543A (zh) * 2008-04-29 2009-11-04 香港理工大学 弹性成像的二维位移估计方法
CN102048560A (zh) * 2010-12-14 2011-05-11 深圳市蓝韵实业有限公司 一种采用双尺度的生物组织位移估计方法
CN102824193A (zh) * 2011-06-14 2012-12-19 深圳迈瑞生物医疗电子股份有限公司 一种弹性成像中的位移检测方法、装置及***
CN102764141A (zh) * 2012-07-20 2012-11-07 中国科学院深圳先进技术研究院 弹性成像方法和***及其中的生物组织位移估计方法和***

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103549978A (zh) * 2013-11-08 2014-02-05 中国科学院声学研究所 一种瞬时弹性成像快速时移估计方法
CN103735287B (zh) * 2013-12-05 2015-11-18 中国科学院苏州生物医学工程技术研究所 一种血管内超声弹性成像二维多级混合位移估计方法
CN103735287A (zh) * 2013-12-05 2014-04-23 中国科学院苏州生物医学工程技术研究所 一种血管内超声弹性成像二维多级混合位移估计方法
CN103654865B (zh) * 2013-12-26 2016-06-22 华南理工大学 基于最大互信息的超声弹性成像组织位移估计方法
CN103654865A (zh) * 2013-12-26 2014-03-26 华南理工大学 基于最大互信息的超声弹性成像组织位移估计方法
CN105266849A (zh) * 2014-07-09 2016-01-27 无锡祥生医学影像有限责任公司 实时超声弹性成像方法和***
CN105232087A (zh) * 2015-11-05 2016-01-13 无锡祥生医学影像有限责任公司 超声弹性成像实时处理***
CN105232087B (zh) * 2015-11-05 2018-01-09 无锡祥生医疗科技股份有限公司 超声弹性成像实时处理***
CN109640830A (zh) * 2016-07-14 2019-04-16 医视特有限公司 基于先例的超声聚焦
CN110215233A (zh) * 2019-04-30 2019-09-10 深圳大学 一种基于超声平面波扫描的分段式脉搏波成像方法
CN113925598A (zh) * 2020-07-13 2022-01-14 通用电气精准医疗有限责任公司 利用超声成像***对组织进行热监测的方法和***
CN113925598B (zh) * 2020-07-13 2024-06-11 通用电气精准医疗有限责任公司 利用超声成像***对组织进行热监测的方法和***
CN115063599A (zh) * 2022-06-16 2022-09-16 昆明理工大学 一种应用于中小型水库大坝监测的小波光流估计与图像相关变形识别方法

Also Published As

Publication number Publication date
CN103040488B (zh) 2014-06-04

Similar Documents

Publication Publication Date Title
CN103040488B (zh) 一种实时超声弹性成像位移估计方法和***
Crosby et al. 3-D speckle tracking for assessment of regional left ventricular function
CN101901481B (zh) 一种图像拼接方法
Rabben et al. Ultrasound-based vessel wall tracking: an auto-correlation technique with RF center frequency estimation
CN101569543B (zh) 弹性成像的二维位移估计方法
CN104688266A (zh) 三维弹性超声成像中的运动修正
US9607405B2 (en) Method and device for detecting displacement in elastography
Basarab et al. A method for vector displacement estimation with ultrasound imaging and its application for thyroid nodular disease
Behar et al. Improving motion estimation by accounting for local image distortion
CN104546014A (zh) 一种用于生物组织弹性测量的剪切波速度估计方法
CN103735287B (zh) 一种血管内超声弹性成像二维多级混合位移估计方法
CN105266849A (zh) 实时超声弹性成像方法和***
Deprez et al. 3D estimation of soft biological tissue deformation from radio-frequency ultrasound volume acquisitions
Ara et al. Phase-based direct average strain estimation for elastography
Omidyeganeh et al. Estimation of strain elastography from ultrasound radio-frequency data by utilizing analytic gradient of the similarity metric
JP5559464B2 (ja) 超音波弾性映像を形成するためのシステム及び方法
CN102217953A (zh) 基于多邻域辅助二维超声形变组织图像跟踪方法及装置
US20160063742A1 (en) Method and system for enhanced frame rate upconversion in ultrasound imaging
CN102217952A (zh) 基于心肌运动轨迹的向量环图生成方法和装置
CN107205725B (zh) 通过实时超声应变成像对心肌梗塞的评估
WO2016039100A1 (ja) 超音波診断装置
Yuan et al. A strain-based ultrasound elastography using phase shift with prior estimates and meshfree shape function
US20100130862A1 (en) Providing Volume Information On A Periodically Moving Target Object In An Ultrasound System
Manuchehri et al. A robust time delay estimation method for ultrasonic echo signals and elastography
Foroughi et al. Tracked ultrasound elastography (TrUE)

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: 20140604

Termination date: 20171221