本发明专利技术公开了一种超声弹性成像性能增强方法，选择高质量的可靠的相邻已估计位移作为当前窗的初始迭代位移；使用二维相位根搜索位移估计改进算法进行迭代产生位移估计值并计算对应相关系数；使用二维窗位移位置估计改进算法估计位移对应的位置，然后通过线性插值法将位移估计值修正到追踪窗中心点；使用梯度法进行纵向应变估计；对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。该超声弹性成像性能增强方法可以避免错误位移的传播，减少组织横向运动及信号解相关对轴向位移估计的影响，抑制信号幅度随机波动引起的幅度调制噪声，提高成像品质，增强图像性能，使图像更加连续、平滑、清晰。
【技术实现步骤摘要】
本专利技术属于医学超声弹性成像
，具体涉及,主要包括位移估计及位移修正技术。
技术介绍
超声(Ultrasound)弹性成像(Strain Imaging)是医学超声成像的新模式,主要是利用超声对组织力学特性参数进行检测并进行成像。当人体组织发生病变时，通常表现为其力学参数发生改变。如乳腺癌、前列腺癌和肝纤维化等疾病发生时，其对应的弹性模量(硬度)发生改变。弹性成像为临床上多种疾病尤其肿瘤的检测与诊断提供重要手段。准静态压缩法是超声弹性成像常用的成像方法，该方法通过对组织施加静态或准静态激励，然后通过对组织压缩前后回波信号(/或图像)分析估计出组织的位移(Displacement)分布；再然后通过对位移分布进行梯度操作产生组织内部的应变(Strain)分布。在相同应力下，组织的不同应变反映了组织的弹性(硬度)信息，将应变分布映射成灰度图像就是反映组织弹性特征的弹性图像。位移估计和应变估计是弹性成像过程的两大核心环节，算法的特性直接影响成像质量。相位根搜索(Phase Root Seeking)法是弹性成像位移估计常采用的算法之一,该算法简单高效且具有较高的实时性。其基本原理对压缩前的一维估计窗，在压缩后信号中使用牛顿迭代法搜索与之互相关函数的相位过零点的信号窗，两信号窗的时移即为所求位移。第一行估计窗初始迭代位移设为0，其它各行的初始迭代位移使用同一列上一个已估计位移。该算法在实际应用中存在以下缺陷(I)会出现错误位移传播。该算法将扫描线上上一个已估计位移作为当前追踪窗的初始迭代位移，如果上一个位移恰好是错误的，这个错误位移将会继续向下传播，尤其当对复杂组织成像或帧间应变较大时该现象较明显；(2)当组织存在横向运动时，图像出现横向不连续性。算法使用一维信号窗进行轴向位移估计，如果在压缩时组织产生横向运动(这在准静态压缩超声弹性成像中常会发生)，压缩前后信号解相关增强，一维窗进行轴向位移估计时常出现图像不连续性，严重时会产生条纹状弹性图像。弹性成像中应变估计通常采用梯度法，通过对位移场进行轴向梯度操作产生轴向应变场。该方法的原理是基于位移估计值是追踪窗中心点位移的假设。实际上相位根搜索位移估计法由于受信号幅度随机波动的影响，该算法产生的位移估计值并不是窗中心点的真正位移。此时仍以传统梯度法进行应变估计将产生显著的幅度调制噪声。
技术实现思路
针对现有技术中存在的上述不足，本专利技术提供了一种主要为了避免错误位移传播和广生更闻精度位移估计和应变估计值，从而提闻弹性质量，获得闻品质弹性图像的超声弹性成像性能增强方法。为了解决上述技术问题，本专利技术米用了如下技术方案，该方法包括如下步骤(1)、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换 成基带IQ信号；(2)、对压缩前信号，以位移估计点为中心取二维运动追踪窗；(3)、第一行追踪窗初始迭代位移均设为0，即d’ffl；1；0= 0，其中下标m表示列号，I表示 第一行，O表示初始迭代(第O次);(4)、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下权利要求1.，其特征在于，该方法包括如下步骤 (1)、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号； (2)、对压缩前信号，以位移估计点为中心取二维运动追踪窗； (3)、第一行追踪窗初始迭代位移均设为O，即d’ffl；1；0= O，其中下标m表示列号，I表示第一行，O表示初始迭代； (4)、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下全文摘要本专利技术公开了，选择高质量的可靠的相邻已估计位移作为当前窗的初始迭代位移；使用二维相位根搜索位移估计改进算法进行迭代产生位移估计值并计算对应相关系数；使用二维窗位移位置估计改进算法估计位移对应的位置，然后通过线性插值法将位移估计值修正到追踪窗中心点；使用梯度法进行纵向应变估计；对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。该超声弹性成像性能增强方法可以避免错误位移的传播，减少组织横向运动及信号解相关对轴向位移估计的影响，抑制信号幅度随机波动引起的幅度调制噪声，提高成像品质，增强图像性能，使图像更加连续、平滑、清晰。文档编号A61B8/00GK102920480SQ20121048601公开日2013年2月13日 申请日期2012年11月26日 优先权日2012年11月26日专利技术者崔少国, 彭彩碧, 黄贤英 申请人:重庆理工大学本文档来自技高网...
【技术保护点】
一种超声弹性成像性能增强方法，其特征在于，该方法包括如下步骤：（1）、在准静态压缩组织过程中，超声探头发射超声波，取压缩前后两帧RF回波，转换成基带IQ信号；（2）、对压缩前信号，以位移估计点为中心取二维运动追踪窗；（3）、第一行追踪窗初始迭代位移均设为0，即d’m,1,0=？0，其中下标m表示列号，1表示第一行，0表示初始迭代；（4）、使用二维相位根搜索位移估计算法，迭代产生第一行轴向位移估计值，方法如下：？④上式中：dm,n,l表示第m列第n行追踪窗进行l次迭代后的位移估计值；为侧向窗宽；为轴向窗长；和分别是相邻估计窗之间侧向和轴向的间隔；m，n表示估计窗的位置下标；S1b和S2b分别表示压缩前和压缩后的基带信号；0？表示中心角频率；*代表共轭；j是复信号；arg是求反正切操作；（5）、计算该位移对应的相关系，为下一行初始迭代位移选择做准备；使用复相关归一化相关系数法；⑤上式中：N是常数，是设定的位移估计迭代次数；其它参数涵义与步骤（4）中公式④的参数涵义相同；（6）、对于第二行位移估计，先从它相邻上一行中以它为轴左右各取L列的2L+1个已估计位移中选择相关系数最高对应的位移作为它的初始位移，0≤L≤2，即：？？？？？？？？？？？？？？？？？⑥其中：b满足，其中0时，？？？？⑨上式中，表示第m列第n行追踪窗修正后的位移，表示该窗由步骤（4）得到修正前的位移，是位移对应的纵向坐标，是追踪窗中心点的纵坐标；（11）、使用梯度法进行纵向应变估计；（12）、对应变估计矩阵进行下采样和灰度映射，扫描转换生成弹性图像。412001dest_path_image001.jpg,2012104860120100001dest_path_image002.jpg,321795dest_path_image003.jpg,2012104860120100001dest_path_image004.jpg,320975dest_path_image005.jpg,2012104860120100001dest_path_image006.jpg,729959dest_path_image007.jpg,2012104860120100001dest_path_image008.jpg,830901dest_path_image009.jpg,2012104860120100001dest_path_image010.jpg,796583dest_path_image011.jpg,2012104860120100001dest_path_image012.jpg,2012104860120100001dest_path_image013.jpg,102800dest_path_image004.jpg,977959dest_path_image005.jpg,815465dest_path_image002.jpg,568526dest_path_image003.jpg,870194dest_path_image014.jpg,2012104860120100001dest_path_image015.jpg,856867dest_path_image016.jpg,2012104860120100001dest_path_image017.jpg,181669dest_path_image018.jpg,738421dest_path_image013.jpg,2012104860120100001dest_path_image019.jpg,566700dest_path_image020.jpg,49240dest_path_image013.jpg,658076dest_path_image019.jpg,2012104860120100001dest_path_image021.jpg,956203dest_path_image022.jpg,2012104860120100001dest_path_image023.jpg,904567dest_path_image013.jpg,560938dest_path_image024.jpg,657070dest_path_image019.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员：崔少国，彭彩碧，黄贤英，
申请(专利权)人：重庆理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。