一种二维迭代的超声弹性成像应变估计方法技术

本发明专利技术公开了一种二维迭代的超声弹性成像应变估计方法，运用到超声弹性成像系统中，具有以下优点：重复利用了轴向应变和横向应变的相关性，使用轴向或横向应变估计信息，采用时间拉伸方法对齐了横向或轴向RF回波信号，增强了横向或轴向RF回波信号的相关性，提高了横向或轴向应变估计的准确度；采用子带技术和互相关方法对横向应变进行估计，提高了横向应变估计的准确度；采用三次样条插值方法对横向RF回波信号进行插值，增加了横向RF回波信号的采样率，提高了横向应变估计的准确度；反复多次迭代求解，不断提高轴向应变和横向应变估计的准确度和精度，得到了更准确的泊松比分布图，有助于医生对早期肿瘤等疾病的诊断。

【技术实现步骤摘要】
一种二维迭代的超声弹性成像应变估计方法
本专利技术属于医生超声成像
，具体涉及一种二维迭代的超声弹性成像应变估计方法。
技术介绍
人体组织的硬度(或称之为弹性)是人体组织很重要的一个机械属性。临床医生一直以来对生物组织的机械属性都有研究，很多病变，特别是肿瘤和癌症，和人体器官硬度变化密切相关，临床医生可以根据人体器官硬度的变化推断出该器官是否有疾病。比如，乳腺肿瘤在病变过程中，其在对应的乳房部位会变硬。超声弹性成像是一种新的医学超声成像技术，可以定性或定量测出人体组织的各器官的软硬度信息、应变和弹性模量，医生通过超声弹性成像获得的人体组织的应变和弹性模量信息对癌症、肿瘤等疾病的早期诊断具有重要的意义。超声弹性成像过程的主要步骤包括：(1)用探头轻轻沿扫描平面挤压人体表面，使人体组织内部产生应力场，人体组织产生轻微的形变；(2)根据压缩前、后的超声回波射频RF回波信号，估计每条A线的每个相关窗口的位移，得到一帧位移分布图；(3)对位移进行求导，得到一帧应变分布图，并将它以图像的方式显示在计算机屏幕上。在现有的超声弹性成像系统中，一般只考虑对轴向应变(沿声束传播方向)进行估计，这是因为在弹性成像过程中，探头是沿轴向对人体组织进行施压，人体组织形变主要沿轴向方向。但在一个实际的弹性成像系统中，人体组织内应力场是三维的，有轴向应力、横向应力(在扫描平面垂直声束传播方向)和正向应力(垂直扫描平面)。由于现有的超声探头扫描得到都是二维的RF回波信号，很难从二维RF回波信号求解正向应力，所以目前的弹性成像系统中都忽略对正向应力的估计。横向位移和应变估计可以提供额外的人体机械属性信息，利用横向位移和应变估计技术，泊松比可以被估计出来，泊松比是一个很重要的生物组织机械属性,泊松比有助于医生对早期肿瘤等疾病的诊断。对于一个二维的应变场(同时存在轴向和横向应变)，横向应变会导致轴向回波信号的解相关噪声，轴向应变也会导致横向回波信号的解相关噪声，并且通常超声系统横向采样率比轴向采样率低很多(通常低几倍或十几倍)，这些因素给横向位移估计带来了很大难度，导致横向应变估计的准确度和精度不高。在实际的超声弹性成像系统中，由于人体组织的横向位移比较弱小，并且探头的横向采样率比较低，轴向应变也会导致横向回波信号的解相关噪声，这些因素导致对横向应变估计非常困难，所以现有的大部分超声弹性成像方法都只是轴向应变(沿声束传播方向)进行估计，没有对横向应变进行估计，这样缺失了横向应变信息，造成无法得到人体组织的泊松比信息。医生由于缺乏横向应变和泊松比的信息，可能导致对肿瘤等疾病的诊断不够准确。对于一个二维的应变场，横向应变和轴向应变是相关的，现有的超声弹性成像应变估计方法几乎都没有考虑到这种相关性，其实可以利用横向应变的信息提高轴向应变的估计准确度。还有现有的横向应变估计方法得到的横向应变的准确度不高，从而得到的泊松比信息不够准确，影响了超声弹性成像的普及和应用。
技术实现思路
针对现有技术中的上述不足，本专利技术提供的二维迭代的超声弹性成像应变估计方法解决了现有的超声弹性成像方法中没有对横向应变进行估计，造成横向应变信息缺失，导致无法得到人体组织的泊松比信息，进而对肿瘤等疾病的诊断不够准确的问题。为了达到上述专利技术目的，本专利技术采用的技术方案为：一种二维迭代的超声弹性成像应变估计方法，包括以下步骤：S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像Si(z,x)；S4、利用轴向应变估计图像Si(z,x)的应变信息对压缩后的一帧RF回波信号数据进行轴向的时间处理；S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一帧横向应变估计图像Li(z,x)；S7、利用横向应变估计图像Li(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF回波信号数据进行横向的时间处理；S8、判断迭代次数i与最大迭代次数m的关系；若迭代次数i<m，则i＝i+1，并返回步骤S3；否则得到最终的轴向应变图像和横向应变图像并进入步骤S9；S9、将最终的轴向应变图像除以最终的横向应变图像得到泊松比分布图，得到准确的轴向应变和横向应变估计。进一步地，所述步骤S1中：压缩前、后的两帧射频RF回波信号是二维的，沿探头声波发射方向是轴向，垂直于轴向方向是横向；压缩前、后的每一条射频RF回波信号是一维的，所述一维的压缩前、后的RF回波信号为：r1(x)＝s(x)p(x)+n1(x)(1)式中，r1(x)、r2(x)分别表示压缩前和压缩后RF回波信号；s(x)表示人体组织散射子函数，p(x)是点扩散函数；n1(x)和n2(x)表示系统白噪声；α是压缩系数，且α＝1-ε，ε是外加应变。进一步地，所述步骤S2中的最大迭代次数m为4。进一步地，所述步骤S3为：采用互相关方法根据压缩前、后的两帧数据沿轴向进行位移估计，得到一帧轴向应变估计图像；S31、将压缩前、后相对应的两条轴向RF回波信号划分成若干个部分重叠的窗口；上述若干个重叠的窗口的重叠率为75％；S32、求每个窗口的位移估计；S33、求每个窗口的应变估计；上述每个窗口的应变估计是对应窗口位移估计的导数。进一步地，所述步骤S4中的时间处理为时间拉伸处理或时间压缩处理；所述步骤S4具体为：S41、将一帧轴向应变估计图像的一条轴向应变估计数据分成若干个不重叠的窗口；S42、求出每个窗口的平均应变估计值；S43、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的轴向RF回波信号；所述步骤S43中的，轴向的时间处理后RF回波信号为：r3(x)＝r2(αx)＝s(x)*p(αx)+n3(x)(3)式中，r3(x)表示轴向的时间处理后的RF回波信号，α为拉伸或压缩系数，n3(x)表示系统白噪声。进一步地，所述步骤S5中的插值处理方法为三次样条插值方法。进一步地，所述步骤S6具体为使用子带划分方法和互相关方法对重建的横向RF回波信号中的压缩前、后的两帧RF回波信号数据进行横向位移估计；其具体方法为：S61、将重建的压缩前、后的每条轴向RF回波信号划分成三个子带信号；S62、对每一对压缩前、后的子带信号，采用互相关方法沿横向进行位移估计得到三个子带横向应变估计图像；S63、对三个子带横向应变估计图像进行加权平均，得到一幅横向应变估计图像。进一步地，所述步骤S7中的时间处理为时间拉伸处理或时间压缩处理；具体为：S71、将一帧横向应变估计图像的一条横向应变估计数据分成若干个不重叠的窗口；S72、求出每个窗口的平均应变估计值；S73、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的横向RF回波信号；所述步骤S73中的，横向的时间处理后RF回波信号为：r4(x)＝r2(αx)＝s(x)*p(αx)+n4(x)(4)式中，r4(x)表示横向的时间处理后的RF回波信号，α为拉伸或压缩系数，n4(x)表示系统白噪声。本文档来自技高网...

【技术保护点】
1.一种二维迭代的超声弹性成像应变估计方法，其特征在于，包括以下步骤：S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像Si(z,x)；S4、利用轴向应变估计图像Si(z,x)的应变信息对压缩后的一帧RF回波信号数据进行轴向的时间处理；S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一帧横向应变估计图像Li(z,x)；S7、利用横向应变估计图像Li(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF回波信号数据进行横向的时间处理；S8、判断迭代次数i与最大迭代次数m的关系；若迭代次数i

【技术特征摘要】
1.一种二维迭代的超声弹性成像应变估计方法，其特征在于，包括以下步骤：S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像Si(z,x)；S4、利用轴向应变估计图像Si(z,x)的应变信息对压缩后的一帧RF回波信号数据进行轴向的时间处理；S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一帧横向应变估计图像Li(z,x)；S7、利用横向应变估计图像Li(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF回波信号数据进行横向的时间处理；S8、判断迭代次数i与最大迭代次数m的关系；若迭代次数i<m，则i＝i+1，并返回步骤S3；否则得到最终的轴向应变图像和横向应变图像并进入步骤S9；S9、将最终的轴向应变图像除以最终的横向应变图像得到泊松比分布图，得到准确的轴向应变和横向应变估计。2.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S1中：压缩前、后的两帧射频RF回波信号是二维的，沿探头声波发射方向是轴向，垂直于轴向方向是横向；压缩前、后的每一条射频RF回波信号是一维的，所述一维的压缩前、后的RF回波信号为：r1(x)＝s(x)p(x)+n1(x)(1)式中，r1(x)、r2(x)分别表示压缩前和压缩后RF回波信号；s(x)表示人体组织散射子函数，p(x)是点扩散函数；n1(x)和n2(x)表示系统白噪声；α是压缩系数，且α＝1-ε，ε是外加应变。3.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S2中的最大迭代次数m为4。4.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S3为：采用互相关方法根据压缩前、后的两帧数据沿轴向进行位移估计，得到一帧轴向应变估计图像；其方法具体为：S31、将压缩前、后相对应的两条轴向RF回波信号划分成若干个部分重叠的窗口；上述若干个重叠的窗口的重叠率为75％；S...

【专利技术属性】
技术研发人员：彭辉，铁菊红，
申请(专利权)人：成都信息工程大学，
类型：发明
国别省市：四川,51