一种超声弹性成像性能增强方法.pdf

摘要
申请专利号：	CN201210486012.0	申请日：	2012.11.26
公开号：	CN102920480A	公开日：	2013.02.13
当前法律状态：	终止	有效性：	无权
法律详情：	未缴年费专利权终止IPC(主分类):A61B 8/08申请日:20121126授权公告日:20141022终止日期:20161126\|\|\|授权\|\|\|实质审查的生效IPC(主分类):A61B 8/00申请日:20121126\|\|\|公开
IPC分类号：	A61B8/00	主分类号：	A61B8/00
申请人：	重庆理工大学
发明人：	崔少国; 彭彩碧; 黄贤英
地址：	400054 重庆市巴南区李家沱红光大道69号
优先权：
专利代理机构：	重庆博凯知识产权代理有限公司 50212	代理人：	穆祥维
PDF下载：	PDF下载

内容摘要

本发明公开了一种超声弹性成像性能增强方法，选择高质量的可靠的相邻已估计位移作为当前窗的初始迭代位移；使用二维相位根搜索位移估计改进算法进行迭代产生位移估计值并计算对应相关系数；使用二维窗位移位置估计改进算法估计位移对应的位置，然后通过线性插值法将位移估计值修正到追踪窗中心点；使用梯度法进行纵向应变估计；对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。该超声弹性成像性能增强方法可以避免错误位移的传播，减少组织横向运动及信号解相关对轴向位移估计的影响，抑制信号幅度随机波动引起的幅度调制噪声，提高成像品质，增强图像性能，使图像更加连续、平滑、清晰。

权利要求书

一种超声弹性成像性能增强方法，其特征在于，该方法包括如下步骤：
（1）、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；
（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；
（3）、第一行追踪窗初始迭代位移均设为0，即d’_m,1,0= 0，其中下标m表示列号，1表示第一行，0表示初始迭代；
（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下：
④
上式中：d_m,n,l表示第m列第n行追踪窗进行l次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；m，n表示估计窗的位置下标；S_1b和S_2b分别表示压缩前和压缩后的基带信号；₀表示中心角频率；*代表共轭；j是复信号；arg是求反正切操作；
（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；
⑤
上式中：N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式④的参数涵义相同；
（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L+1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0≤L≤2，即：
                ⑥
其中：b满足，其中0<=b<=L，m+b是上一行2L+1范围内已估计位移对应相关系数最高的列号，表示求最大值函数；
（7）、使用步骤（4）中的二维相位根搜索位移估计算法迭代产生第二行轴向位移估计值；
（8）、同第二行轴向位移估计方法一样，估计出第三行到最后一行的轴向位移；
（9）、对于步骤（3）～步骤（8）产生的位移值，使用二维窗位移位置估计改进算法估计其对应的轴向位置；该二维窗位移位置估计改进算法如下：
      ⑦
式中，表示第m列第n行追踪窗对应位移的轴向位置坐标，和分别是侧向和轴向相邻追踪窗之间的间隔，和分别是追踪窗的侧向窗宽和轴向窗长，是追踪窗的轴向位移，是权值，，其中和分别是压缩前和压缩后信号的包络；
（10）、根据位移位置估计，使用线性插值法将位移估计值修正到追踪窗中心点，方法如下：
当<时，     ⑧
当>时，    ⑨
上式中，表示第m列第n行追踪窗修正后的位移，表示该窗由步骤（4）得到修正前的位移，是位移对应的纵向坐标，是追踪窗中心点的纵坐标；
（11）、使用梯度法进行纵向应变估计；
（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。

说明书

一种超声弹性成像性能增强方法
技术领域
本发明属于医学超声弹性成像技术领域，具体涉及一种超声弹性成像性能增强方法，主要包括位移估计及位移修正技术。
背景技术
超声（Ultrasound）弹性成像(Strain Imaging)是医学超声成像的新模式，主要是利用超声对组织力学特性参数进行检测并进行成像。当人体组织发生病变时，通常表现为其力学参数发生改变。如乳腺癌、前列腺癌和肝纤维化等疾病发生时，其对应的弹性模量（硬度）发生改变。弹性成像为临床上多种疾病尤其肿瘤的检测与诊断提供重要手段。
准静态压缩法是超声弹性成像常用的成像方法，该方法通过对组织施加静态或准静态激励，然后通过对组织压缩前后回波信号（/或图像）分析估计出组织的位移(Displacement)分布；再然后通过对位移分布进行梯度操作产生组织内部的应变(Strain)分布。在相同应力下，组织的不同应变反映了组织的弹性（硬度）信息，将应变分布映射成灰度图像就是反映组织弹性特征的弹性图像。
位移估计和应变估计是弹性成像过程的两大核心环节，算法的特性直接影响成像质量。相位根搜索（Phase Root Seeking）法是弹性成像位移估计常采用的算法之一，该算法简单高效且具有较高的实时性。其基本原理：对压缩前的一维估计窗，在压缩后信号中使用牛顿迭代法搜索与之互相关函数的相位过零点的信号窗，两信号窗的时移即为所求位移。第一行估计窗初始迭代位移设为0，其它各行的初始迭代位移使用同一列上一个已估计位移。
该算法在实际应用中存在以下缺陷：(1)会出现错误位移传播。该算法将扫描线上上一个已估计位移作为当前追踪窗的初始迭代位移，如果上一个位移恰好是错误的，这个错误位移将会继续向下传播，尤其当对复杂组织成像或帧间应变较大时该现象较明显；(2)当组织存在横向运动时，图像出现横向不连续性。算法使用一维信号窗进行轴向位移估计，如果在压缩时组织产生横向运动（这在准静态压缩超声弹性成像中常会发生），压缩前后信号解相关增强，一维窗进行轴向位移估计时常出现图像不连续性，严重时会产生条纹状弹性图像。
弹性成像中应变估计通常采用梯度法，通过对位移场进行轴向梯度操作产生轴向应变场。该方法的原理是基于位移估计值是追踪窗中心点位移的假设。实际上相位根搜索位移估计法由于受信号幅度随机波动的影响，该算法产生的位移估计值并不是窗中心点的真正位移。此时仍以传统梯度法进行应变估计将产生显著的幅度调制噪声。
发明内容
针对现有技术中存在的上述不足，本发明提供了一种主要为了避免错误位移传播和产生更高精度位移估计和应变估计值，从而提高弹性质量，获得高品质弹性图像的超声弹性成像性能增强方法。
为了解决上述技术问题，本发明采用了如下技术方案：
一种超声弹性成像性能增强方法，该方法包括如下步骤：
（1）、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；
（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；
（3）、第一行追踪窗初始迭代位移均设为0，即d’_m,1,0= 0，其中下标m表示列号，1表示第一行，0表示初始迭代（第0次）；
（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下：
④
上式中：d_m,n,l表示第m列第n行追踪窗进行l次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；m，n表示估计窗的位置下标；S_1b和S_2b分别表示压缩前和压缩后的基带信号；₀表示中心角频率；*代表共轭；j是复信号；arg是求反正切操作；
（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；
⑤
上式中：N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式④的参数涵义相同；
（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L+1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0≤L≤2，即：
                ⑥
其中：b满足，其中0<=b<=L，m+b是上一行2L+1范围内已估计位移对应相关系数最高的列号，表示求最大值函数；
（7）、使用步骤（4）中的二维相位根搜索位移估计算法迭代产生第二行轴向位移估计值；
（8）、同第二行轴向位移估计方法一样，估计出第三行到最后一行的轴向位移；
（9）、对于步骤（3）～步骤（8）产生的位移值，使用二维窗位移位置估计改进算法估计其对应的轴向位置；该二维窗位移位置估计改进算法如下：
      ⑦
式中，表示第m列第n行追踪窗对应位移的轴向位置坐标，和分别是侧向和轴向相邻追踪窗之间的间隔，和分别是追踪窗的侧向窗宽和轴向窗长，是追踪窗的轴向位移，是权值，，其中和分别是压缩前和压缩后信号的包络；
（10）、根据位移位置估计，使用线性插值法将位移估计值修正到追踪窗中心点，方法如下：
当<时，     ⑧
当>时，    ⑨
上式中，表示第m列第n行追踪窗修正后的位移，表示该窗修正前的位移（步骤（4）得到的位移），是位移对应的纵向坐标，是追踪窗中心点的纵坐标；
（11）、使用梯度法进行纵向应变估计；
（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。
本发明的有益效果是：初始迭代位移选择方法的基本思想是选择上一行与当前位置相邻的已估计位移中质量最高的位移作为其初始迭代位移，此方法可以有效避免错误位移传播。基于二维信号窗的轴向位移估计方法可以减少组织横向运动及信号解相关对位移估计的影响，使图像更加连续平滑。基于位移位置估计的位移修正方法的基本思想是使用位移位置估计将已估计的位移使用插值法修正到追踪窗中心点，可以抑制信号幅度随机波动引起的弹性图像噪声，提高成像品质，增强超声弹性成像性能。
附图说明
图1为初始迭代位移选择方法示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步说明。
一种超声弹性成像性能增强方法，该方法包括如下步骤：
（1）、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；
（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；
（3）、第一行追踪窗初始迭代位移均设为0，即d’_m,1,0= 0，其中下标m表示列号，1表示第一行，0表示初始迭代（第0次）；
（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下：
④
上式中：d_m,n,l表示第m列第n行追踪窗进行l次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；m，n表示估计窗的位置下标；S_1b和S_2b分别表示压缩前和压缩后的基带信号；₀表示中心角频率；*代表共轭；j是复信号；arg是求反正切操作；
（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；
⑤
上式中：N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式④的参数涵义相同；
（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L+1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0≤L≤2，即：
                ⑥
其中：b满足，其中0<=b<=L，m+b是上一行2L+1范围内已估计位移对应相关系数最高的列号，表示求最大值函数；
（7）、使用步骤（4）中的二维相位根搜索位移估计算法迭代产生第二行轴向位移估计值；
（8）、同第二行轴向位移估计方法一样，估计出第三行到最后一行的轴向位移；
（9）、对于步骤（3）～步骤（8）产生的位移值，使用二维窗位移位置估计改进算法估计其对应的轴向位置；该二维窗位移位置估计改进算法如下：
      ⑦
式中，表示第m列第n行追踪窗对应位移的轴向位置坐标，和分别是侧向和轴向相邻追踪窗之间的间隔，和分别是追踪窗的侧向窗宽和轴向窗长，是追踪窗的轴向位移，是权值，，其中和分别是压缩前和压缩后信号的包络；
（10）、根据位移位置估计，使用线性插值法将位移估计值修正到追踪窗中心点，方法如下：
当<时，     ⑧
当>时，    ⑨
上式中，表示第m列第n行追踪窗修正后的位移，表示该窗修正前的位移（步骤（4）得到的位移），是位移对应的纵向坐标，是追踪窗中心点的纵坐标；
（11）、使用梯度法进行纵向应变估计；
（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。
图1为初始迭代位移选择方法示意图，图中d’表示初始迭代位移，d表示使用d’作为初始迭代位移估计出来的位移（根据步骤（4）），下标表示追踪窗所在位置坐标，d后面括号里的数值表示该位移对应的相关系数（根据步骤（5））。空心箭头表示初始位移的来源，实心箭头表示位移估计方向。以图1中第4行第2列方格为例，当前预求第22列23行追踪窗的位移，而上一行已估计的3个相邻位移（这里L=1）d_21,22、d_22,22 和d_23,22中，右上侧窗（第23列22行）对应位移d_23,22有最大相关系数0.99，因此使用该位移作为当前窗（第22列23行）的初始迭代位移，即：d’_22,23= d_23,22。然后使用步骤（4）进行迭代产生最终的位移d_22,23(0.98)。括号里0.98表示该位移对应的相关系数为0.98，使用步骤（5）计算得到。
通过本发明的方法产生的弹性图像，可以有效避免错误位移传播，减少横向运动对位移估计的影响以及抑制弹性成像幅度调制噪声，可显著提高成像品质，增强超声弹性成像性能。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

资源描述

《一种超声弹性成像性能增强方法.pdf》由会员分享，可在线阅读，更多相关《一种超声弹性成像性能增强方法.pdf（9页珍藏版）》请在专利查询网上搜索。

1、10申请公布号CN102920480A43申请公布日20130213CN102920480ACN102920480A21申请号201210486012022申请日20121126A61B8/0020060171申请人重庆理工大学地址400054重庆市巴南区李家沱红光大道69号72发明人崔少国彭彩碧黄贤英74专利代理机构重庆博凯知识产权代理有限公司50212代理人穆祥维54发明名称一种超声弹性成像性能增强方法57摘要本发明公开了一种超声弹性成像性能增强方法，选择高质量的可靠的相邻已估计位移作为当前窗的初始迭代位移；使用二维相位根搜索位移估计改进算法进行迭代产生位移估计值并计算对应相关系数；使用二。

2、维窗位移位置估计改进算法估计位移对应的位置，然后通过线性插值法将位移估计值修正到追踪窗中心点；使用梯度法进行纵向应变估计；对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。该超声弹性成像性能增强方法可以避免错误位移的传播，减少组织横向运动及信号解相关对轴向位移估计的影响，抑制信号幅度随机波动引起的幅度调制噪声，提高成像品质，增强图像性能，使图像更加连续、平滑、清晰。51INTCL权利要求书2页说明书5页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书5页附图1页1/2页21一种超声弹性成像性能增强方法，其特征在于，该方法包括如下步骤（1）、在准静态压缩组织过。

3、程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；（3）、第一行追踪窗初始迭代位移均设为0，即DM,1,00，其中下标M表示列号，1表示第一行，0表示初始迭代；（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下上式中DM,N,L表示第M列第N行追踪窗进行L次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；M，N表示估计窗的位置下标；S1B和S2B分别表示压缩前和压缩后的基带信号；0表示中心角频率；代表共轭；J是复信号；ARG是求反正切操作；（5）、计算该。

4、位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；上式中N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式的参数涵义相同；（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0L2，即其中B满足，其中0时，上式中，表示第M列第N行追踪窗修正后的位移，表示该窗由步骤（4）得到修正前的位移，是位移对应的纵向坐标，是追踪窗中心点的纵坐标；（11）、使用梯度法进行纵向应变估计；（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。权利要求书CN102920480A1。

5、/5页4一种超声弹性成像性能增强方法技术领域0001本发明属于医学超声弹性成像技术领域，具体涉及一种超声弹性成像性能增强方法，主要包括位移估计及位移修正技术。背景技术0002超声（ULTRASOUND）弹性成像STRAINIMAGING是医学超声成像的新模式，主要是利用超声对组织力学特性参数进行检测并进行成像。当人体组织发生病变时，通常表现为其力学参数发生改变。如乳腺癌、前列腺癌和肝纤维化等疾病发生时，其对应的弹性模量（硬度）发生改变。弹性成像为临床上多种疾病尤其肿瘤的检测与诊断提供重要手段。0003准静态压缩法是超声弹性成像常用的成像方法，该方法通过对组织施加静态或准静态激励，然后通过对组织。

6、压缩前后回波信号（/或图像）分析估计出组织的位移DISPLACEMENT分布；再然后通过对位移分布进行梯度操作产生组织内部的应变STRAIN分布。在相同应力下，组织的不同应变反映了组织的弹性（硬度）信息，将应变分布映射成灰度图像就是反映组织弹性特征的弹性图像。0004位移估计和应变估计是弹性成像过程的两大核心环节，算法的特性直接影响成像质量。相位根搜索（PHASEROOTSEEKING）法是弹性成像位移估计常采用的算法之一，该算法简单高效且具有较高的实时性。其基本原理对压缩前的一维估计窗，在压缩后信号中使用牛顿迭代法搜索与之互相关函数的相位过零点的信号窗，两信号窗的时移即为所求位移。第一行估计。

7、窗初始迭代位移设为0，其它各行的初始迭代位移使用同一列上一个已估计位移。0005该算法在实际应用中存在以下缺陷1会出现错误位移传播。该算法将扫描线上上一个已估计位移作为当前追踪窗的初始迭代位移，如果上一个位移恰好是错误的，这个错误位移将会继续向下传播，尤其当对复杂组织成像或帧间应变较大时该现象较明显；2当组织存在横向运动时，图像出现横向不连续性。算法使用一维信号窗进行轴向位移估计，如果在压缩时组织产生横向运动（这在准静态压缩超声弹性成像中常会发生），压缩前后信号解相关增强，一维窗进行轴向位移估计时常出现图像不连续性，严重时会产生条纹状弹性图像。0006弹性成像中应变估计通常采用梯度法，通过对位。

8、移场进行轴向梯度操作产生轴向应变场。该方法的原理是基于位移估计值是追踪窗中心点位移的假设。实际上相位根搜索位移估计法由于受信号幅度随机波动的影响，该算法产生的位移估计值并不是窗中心点的真正位移。此时仍以传统梯度法进行应变估计将产生显著的幅度调制噪声。发明内容0007针对现有技术中存在的上述不足，本发明提供了一种主要为了避免错误位移传播和产生更高精度位移估计和应变估计值，从而提高弹性质量，获得高品质弹性图像的超声弹性成像性能增强方法。说明书CN102920480A2/5页50008为了解决上述技术问题，本发明采用了如下技术方案一种超声弹性成像性能增强方法，该方法包括如下步骤（1）、在准静态压缩组。

9、织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；（3）、第一行追踪窗初始迭代位移均设为0，即DM,1,00，其中下标M表示列号，1表示第一行，0表示初始迭代（第0次）；（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下上式中DM,N,L表示第M列第N行追踪窗进行L次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；M，N表示估计窗的位置下标；S1B和S2B分别表示压缩前和压缩后的基带信号；0表示中心角频率；代表共轭；J是复信号；ARG是求反正切操作；。

10、（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；上式中N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式的参数涵义相同；（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0L2，即其中B满足，其中0时，上式中，表示第M列第N行追踪窗修正后的位移，表示该窗修正前的位移（步骤（4）得到的位移），是位移对应的纵向坐标，是追踪窗中心点的纵坐标；（11）、使用梯度法进行纵向应变估计；（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。0009本发明。

11、的有益效果是初始迭代位移选择方法的基本思想是选择上一行与当前位置相邻的已估计位移中质量最高的位移作为其初始迭代位移，此方法可以有效避免错误位移传播。基于二维信号窗的轴向位移估计方法可以减少组织横向运动及信号解相关对位移估计的影响，使图像更加连续平滑。基于位移位置估计的位移修正方法的基本思想是使用位移位置估计将已估计的位移使用插值法修正到追踪窗中心点，可以抑制信号幅度随机波动引起的弹性图像噪声，提高成像品质，增强超声弹性成像性能。附图说明0010图1为初始迭代位移选择方法示意图。具体实施方式0011下面结合附图和具体实施方式对本发明作进一步说明。0012一种超声弹性成像性能增强方法，该方法包括如。

12、下步骤（1）、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；说明书CN102920480A4/5页7（3）、第一行追踪窗初始迭代位移均设为0，即DM,1,00，其中下标M表示列号，1表示第一行，0表示初始迭代（第0次）；（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下上式中DM,N,L表示第M列第N行追踪窗进行L次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；M，N表示估计窗的位置下标；S1B和S2B分别表示压缩前和压缩后的基带。

13、信号；0表示中心角频率；代表共轭；J是复信号；ARG是求反正切操作；（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；上式中N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式的参数涵义相同；（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0L2，即其中B满足，其中0时，上式中，表示第M列第N行追踪窗修正后的位移，表示该窗修正前的位移（步骤（4）得到的位移），是位移对应的纵向坐标，是追踪窗中心点的纵坐标；（11）、使用梯度法进行纵向应变估计；（12）、对应。

14、变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。0013图1为初始迭代位移选择方法示意图，图中D表示初始迭代位移，D表示使用D作为初始迭代位移估计出来的位移（根据步骤（4），下标表示追踪窗所在位置坐标，D后面括号里的数值表示该位移对应的相关系数（根据步骤（5）。空心箭头表示初始位移的来源，实心箭头表示位移估计方向。以图1中第4行第2列方格为例，当前预求第22列23行追踪窗的位移，而上一行已估计的3个相邻位移（这里L1）D21,22、D22,22和D23,22中，右上侧窗（第23列22行）对应位移D23,22有最大相关系数099，因此使用该位移作为当前窗（第22列23行）的初始迭代位移，即。

15、D22,23D23,22。然后使用步骤（4）进行迭代产生最终的位移D22,23098。括号里098表示该位移对应的相关系数为098，使用步骤（5）计算得到。0014通过本发明的方法产生的弹性图像，可以有效避免错误位移传播，减少横向运动对位移估计的影响以及抑制弹性成像幅度调制噪声，可显著提高成像品质，增强超声弹性成像性能。0015最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。说明书CN102920480A1/1页9图1说明书附图CN102920480A。

展开阅读全文