CN108784736A - 一种二维迭代的超声弹性成像应变估计方法 - Google Patents

Abstract

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

Description

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

技术领域

本发明属于医生超声成像技术领域，具体涉及一种二维迭代的超声弹性成像应变估计方法。

背景技术

人体组织的硬度(或称之为弹性)是人体组织很重要的一个机械属性。临床医生一直以来对生物组织的机械属性都有研究，很多病变，特别是肿瘤和癌症，和人体器官硬度变化密切相关，临床医生可以根据人体器官硬度的变化推断出该器官是否有疾病。比如，乳腺肿瘤在病变过程中，其在对应的乳房部位会变硬。超声弹性成像是一种新的医学超声成像技术，可以定性或定量测出人体组织的各器官的软硬度信息、应变和弹性模量，医生通过超声弹性成像获得的人体组织的应变和弹性模量信息对癌症、肿瘤等疾病的早期诊断具有重要的意义。

超声弹性成像过程的主要步骤包括：

(1)用探头轻轻沿扫描平面挤压人体表面，使人体组织内部产生应力场，人体组织产生轻微的形变；

(2)根据压缩前、后的超声回波射频RF回波信号，估计每条A线的每个相关窗口的位移，得到一帧位移分布图；

(3)对位移进行求导，得到一帧应变分布图，并将它以图像的方式显示在计算机屏幕上。

在现有的超声弹性成像系统中，一般只考虑对轴向应变(沿声束传播方向)进行估计，这是因为在弹性成像过程中，探头是沿轴向对人体组织进行施压，人体组织形变主要沿轴向方向。但在一个实际的弹性成像系统中，人体组织内应力场是三维的，有轴向应力、横向应力(在扫描平面垂直声束传播方向)和正向应力(垂直扫描平面)。由于现有的超声探头扫描得到都是二维的RF回波信号，很难从二维RF回波信号求解正向应力，所以目前的弹性成像系统中都忽略对正向应力的估计。横向位移和应变估计可以提供额外的人体机械属性信息，利用横向位移和应变估计技术，泊松比可以被估计出来，泊松比是一个很重要的生物组织机械属性,泊松比有助于医生对早期肿瘤等疾病的诊断。对于一个二维的应变场(同时存在轴向和横向应变)，横向应变会导致轴向回波信号的解相关噪声，轴向应变也会导致横向回波信号的解相关噪声，并且通常超声系统横向采样率比轴向采样率低很多(通常低几倍或十几倍)，这些因素给横向位移估计带来了很大难度，导致横向应变估计的准确度和精度不高。

在实际的超声弹性成像系统中，由于人体组织的横向位移比较弱小，并且探头的横向采样率比较低，轴向应变也会导致横向回波信号的解相关噪声，这些因素导致对横向应变估计非常困难，所以现有的大部分超声弹性成像方法都只是轴向应变(沿声束传播方向)进行估计，没有对横向应变进行估计，这样缺失了横向应变信息，造成无法得到人体组织的泊松比信息。医生由于缺乏横向应变和泊松比的信息，可能导致对肿瘤等疾病的诊断不够准确。

对于一个二维的应变场，横向应变和轴向应变是相关的，现有的超声弹性成像应变估计方法几乎都没有考虑到这种相关性，其实可以利用横向应变的信息提高轴向应变的估计准确度。还有现有的横向应变估计方法得到的横向应变的准确度不高，从而得到的泊松比信息不够准确，影响了超声弹性成像的普及和应用。

发明内容

针对现有技术中的上述不足，本发明提供的二维迭代的超声弹性成像应变估计方法解决了现有的超声弹性成像方法中没有对横向应变进行估计，造成横向应变信息缺失，导致无法得到人体组织的泊松比信息，进而对肿瘤等疾病的诊断不够准确的问题。

为了达到上述发明目的，本发明采用的技术方案为：一种二维迭代的超声弹性成像应变估计方法，包括以下步骤：

S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；

S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；

S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像S_i(z,x)；

S4、利用轴向应变估计图像S_i(z,x)的应变信息对压缩后的一帧RF回波信号数据进行轴向的时间处理；

S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；

S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一帧横向应变估计图像L_i(z,x)；

S7、利用横向应变估计图像L_i(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF回波信号数据进行横向的时间处理；

S8、判断迭代次数i与最大迭代次数m的关系；

若迭代次数i<m，则i＝i+1，并返回步骤S3；

否则得到最终的轴向应变图像和横向应变图像并进入步骤S9；

S9、将最终的轴向应变图像除以最终的横向应变图像得到泊松比分布图，得到准确的轴向应变和横向应变估计。

进一步地，所述步骤S1中：

压缩前、后的两帧射频RF回波信号是二维的，沿探头声波发射方向是轴向，垂直于轴向方向是横向；

压缩前、后的每一条射频RF回波信号是一维的，所述一维的压缩前、后的RF回波信号为：

r₁(x)＝s(x)p(x)+n₁(x) (1)

式中，r₁(x)、r₂(x)分别表示压缩前和压缩后RF回波信号；s(x)表示人体组织散射子函数，p(x)是点扩散函数；n₁(x)和n₂(x)表示系统白噪声；α是压缩系数，且α＝1-ε，ε是外加应变。

进一步地，所述步骤S2中的最大迭代次数m为4。

进一步地，所述步骤S3为：采用互相关方法根据压缩前、后的两帧数据沿轴向进行位移估计，得到一帧轴向应变估计图像；

S31、将压缩前、后相对应的两条轴向RF回波信号划分成若干个部分重叠的窗口；

上述若干个重叠的窗口的重叠率为75％；

S32、求每个窗口的位移估计；

S33、求每个窗口的应变估计；

上述每个窗口的应变估计是对应窗口位移估计的导数。

进一步地，所述步骤S4中的时间处理为时间拉伸处理或时间压缩处理；

所述步骤S4具体为：

S41、将一帧轴向应变估计图像的一条轴向应变估计数据分成若干个不重叠的窗口；

S42、求出每个窗口的平均应变估计值；

S43、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的轴向RF回波信号；

所述步骤S43中的，轴向的时间处理后RF回波信号为：

r₃(x)＝r₂(αx)＝s(x)*p(αx)+n₃(x) (3)

式中，r₃(x)表示轴向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₃(x)表示系统白噪声。

进一步地，所述步骤S5中的插值处理方法为三次样条插值方法。

进一步地，所述步骤S6具体为使用子带划分方法和互相关方法对重建的横向RF回波信号中的压缩前、后的两帧RF回波信号数据进行横向位移估计；

其具体方法为：

S61、将重建的压缩前、后的每条轴向RF回波信号划分成三个子带信号；

S62、对每一对压缩前、后的子带信号，采用互相关方法沿横向进行位移估计得到三个子带横向应变估计图像；

S63、对三个子带横向应变估计图像进行加权平均，得到一幅横向应变估计图像。

进一步地，所述步骤S7中的时间处理为时间拉伸处理或时间压缩处理；

具体为：

S71、将一帧横向应变估计图像的一条横向应变估计数据分成若干个不重叠的窗口；

S72、求出每个窗口的平均应变估计值；

S73、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的横向RF回波信号；

所述步骤S73中的，横向的时间处理后RF回波信号为：

r4(x)＝r₂(αx)＝s(x)*p(αx)+n₄(x) (4)

式中，r₄(x)表示横向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₄(x)表示系统白噪声。

进一步地，所述步骤S9中泊松比为人体弹性组织的重要机械属性，泊松比ν为：

上式中，ε_l是横向应变，ε_a是轴向应变。

本发明提供的二维迭代的超声弹性成像应变估计方法运用到超声弹性成像系统中，具有以下优点：

(1)重复利用了轴向应变和横向应变的相关性，使用轴向或横向应变估计信息，采用时间拉伸方法对齐了横向或轴向RF回波信号，增强了横向或轴向RF回波信号的相关性，提高了横向或轴向应变估计的准确度；

(2)采用子带技术和互相关方法对横向应变进行估计，提高了横向应变估计的准确度；

(3)采用三次样条插值方法对横向RF回波信号进行插值，增加了横向RF回波信号的采样率，提高了横向应变估计的准确度；

(4)反复多次迭代求解，不断提高轴向应变和横向应变估计的准确度和精度，得到了更准确的泊松比分布图，有助于医生对早期肿瘤等疾病的诊断。

附图说明

图1为本发明提供的实施例中二维迭代的超声弹性成像应变估计方法实现流程图。

图2为本发明提供的实施例中采用互相关方法得到一帧轴向应变估计图像方法实现流程图。

图3为本发明提供的实施例中轴向的时间处理方法实现流程图。

图4为本发明提供的实施例中使用子带划分方法和互相关方法进行横向位移估计实现流程图。

图5为本发明提供的实施例中横向的时间处理方法实现流程图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

如图1所示，一种二维迭代的超声弹性成像应变估计方法，包括以下步骤：

S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；

上述步骤S1中：

压缩前、后的两帧射频RF回波信号是二维的，沿探头声波发射方向是轴向，垂直于轴向方向是横向；一帧二维离散化的RF回波信号可以看做是一个二维信号矩阵，轴向就是矩阵的列方向，横向就是矩阵的行方向，轴向条数(列数)取决于探头的宽度和扫描间隔，横向条数(行数)取决于沿轴向的信号采样率。所以，一个二维的RF回波信号可以看做由若干条一维的轴向信号构成，或者可以看做是由若干条一维的横向信号构成。

压缩前、后的每一条射频RF回波信号是一维的，一维的压缩后的RF回波信号为：

r₁(x)＝s(x)p(x)+n₁(x) (1)

式中，r₁(x)、r₂(x)分别表示压缩前、后的一条RF回波信号；s(x)表示人体组织散射子函数，p(x)是点扩散函数；n₁(x)和n₂(x)表示系统白噪声；α是压缩系数，且α＝1-ε，ε是外加应变。

S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；

上述最大迭代次数m一般为4。

S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像S_i(z,x)；

其中i表示迭代次数，z和x分别表示轴向和横向坐标；

上述步骤S3具体为：采用互相关方法根据压缩前、后的两帧数据沿轴向进行位移估计，得到一帧轴向应变估计图像；

其方法具体为：

S31、将压缩前、后相对应的两条轴向RF回波信号划分成若干个部分重叠的窗口；

上述若干个重叠的窗口的重叠率为75％；

S32、求每个窗口的位移估计；

其中，第n个窗口的位移估计的计算公式为：

式中，dn是第n个窗口的位移估计；

f_n(τ)是第n个窗口的互相关函数；

argmax(f_n(τ))表示求f_n(τ)函数最大值所对应的参数τ；

ΔT表示相邻两个窗口的间距；

T是窗口的长度；

r₁(x)、r₂(x)分别表示压缩前和压缩后的RF回波信号。

S33、求每个窗口的应变估计；

上述每个窗口的应变估计是对应窗口位移估计的导数，第n个窗口的应变估计的计算公式为：

式中，s_n表示第n个窗口的应变估计；

d_n、d_n+1分别表示第n和n+1窗口的位移估计；

_ΔT表示相邻两个窗口的间距。

S4、利用轴向应变估计图像S_i(z,x)的应变信息对压缩后的一帧RF数据进行轴向的时间处理；

其中，时间处理为时间拉伸处理或时间压缩处理，其目的为将压缩前的RF数据与压缩后的RF数据进行轴向的对齐，降低散射子横向位移对轴向信号的解相关噪声；

上述步骤S4具体为：

S41、将一帧轴向应变估计图像的一条轴向应变估计数据分成若干个不重叠的窗口；

S42、求出每个窗口的平均应变估计值；

S43、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的轴向RF回波信号；

所述步骤S43中的，轴向的时间处理后RF回波信号为：

r₃(x)＝r₂(αx)＝s(x)*p(αx)+n₃(x) (5)

式中，r₃(x)表示轴向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₃(x)表示系统白噪声。

上述公式(5)中，假设应变在一条RF回波信号上是均匀的、一致的，但在实际过程中，应变在整个弹性成像区域是不均匀的，在局部相对小的区域应变变化比较小，可以看作是均匀的。

S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；

上述步骤S5的目的是增加横向采样频率，提高横向应变估计的准确度和精度；插值方法可采用线性插值和三次样条插值等，由于三次样条插值与线性插值相比能够增加相关系数，所以本发明方法采用三次样条插值法进行插值，其具体方法为：

对压缩前、后的每条横向RF回波信号做N:1三次样条采样处理：设一条横向RF回波信号由L个采样点构成，设相邻两个采样点之间的连线符合三次多项式，在相邻两个采样点之间按三次多项式插入N-1个点，插值后的RF回波信号后有N*L个采样点。

S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一幅横向应变估计图像L_i(z,x)；

所述步骤S6具体为使用子带划分方法和互相关方法对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，其具体方法为：

S61、将重建的压缩前、后的每条轴向RF回波信号划分成三个子带信号；

其子带划分公式为：

s_i(x)＝BPF{r(x),f_i1,f_i2} (6)

式中，

S_i(x)是第i个子带信号；

r(x)是一条轴向的RF信号；

f_i1和f_i2分别是第i个子带的开始频率和截止频率；

BPF是带通滤波器。

S62、对每一对压缩前、后的子带信号，采用互相关方法沿横向进行位移估计得到三个子带横向应变估计图像；

S63、对三个子带横向应变估计图像进行加权平均，得到一帧横向应变估计图像。

S7、利用横向应变估计图像L_i(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF数据进行横向的时间处理；

所述步骤S7中的时间处理为时间拉伸处理或时间压缩处理；其目的为将压缩前的RF数据与压缩后的RF数据进行轴向的对齐，降低散射子横向位移对轴向信号的解相关噪声。

具体为：

S71、将一帧横向应变估计图像的一条横向应变估计数据分成若干个不重叠的窗口；

S72、求出每个窗口的平均应变估计值；

S73、通过每个窗口的平均估计值拉伸或压缩该窗口内的对应的横向RF回波信号。

步骤S73中的，横向的时间处理后RF回波信号为：

r₄(x)＝r₂(αx)＝s(x)*p(αx)+n₄(x)

r₄(x)表示横向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₄(x)表示系统白噪声。

S8、判断迭代次数i与最大迭代次数m的关系；

若迭代次数i<m，则i＝i+1，并返回步骤S3；

否则得到最终的轴向应变图像和横向应变图像并进入步骤S9；

上述步骤S4和步骤S7中的横向或轴向应变估计方法是多次迭代的过程，通过在最大迭代次数范围内进行多次迭代，每次迭代得到的轴向应变图像或横向应变图像质量都比上一次要好，弹性信噪比更低，所以每次迭代时，时间拉伸(或压缩)的窗口大小应该不一样，从大到小递减变化，这样时间拉伸或压缩窗口的平均应变估计值更加接近真实值，降低了由于时间拉伸或压缩引起的人为噪声。

S9、将最终的轴向应变图像除以最终的横向应变图像得到泊松比分布图，得到准确的轴向应变和横向应变估计。

所述泊松比为人体弹性组织的重要机械属性，泊松比ν为：

上式中，ε_l是横向应变，ε_a是轴向应变，其中的横向应变和轴向应变分别为轴向应变估计和横向应变图像的应变信息。

本发明提供的二维迭代的超声弹性成像应变估计方法运用到超声弹性成像系统中，具有以下优点：

(1)重复利用了轴向应变和横向应变的相关性，使用轴向或横向应变估计信息，采用时间拉伸方法对齐了横向或轴向RF回波信号，增强了横向或轴向RF回波信号的相关性，提高了横向或轴向应变估计的准确度；

(2)采用子带技术和互相关方法对横向应变进行估计，提高了横向应变估计的准确度；

(3)采用三次样条插值方法对横向RF回波信号进行插值，增加了横向RF回波信号的采样率，提高了横向应变估计的准确度；

(4)反复多次迭代求解，不断提高轴向应变和横向应变估计的准确度和精度，得到了更准确的泊松比分布图，有助于医生对早期肿瘤等疾病的诊断。

Claims (9)

1.一种二维迭代的超声弹性成像应变估计方法，其特征在于，包括以下步骤：

S1、使用超声探头挤压人体组织，实时获取压缩前、后的两帧射频RF回波信号数据；

S2、设置迭代次数i，最大迭代次数为m，且初始迭代次数i＝0；

S3、根据当前迭代的RF回波信号数据，沿轴向方向进行位移估计，得到一帧轴向应变估计图像S_i(z,x)；

S4、利用轴向应变估计图像S_i(z,x)的应变信息对压缩后的一帧RF回波信号数据进行轴向的时间处理；

S5、将时间处理过的压缩后的一帧RF回波信号数据和压缩前的一帧RF回波信号数据进行横向的N：1的插值处理，重建横向RF回波信号；

S6、对重建的横向RF回波信号中的压缩前、后两帧数据进行横向位移估计，得到一帧横向应变估计图像L_i(z,x)；

S7、利用横向应变估计图像L_i(z,x)的应变信息将重建的RF回波信号中的压缩后的一帧RF回波信号数据进行横向的时间处理；

S8、判断迭代次数i与最大迭代次数m的关系；

若迭代次数i<m，则i＝i+1，并返回步骤S3；

否则得到最终的轴向应变图像和横向应变图像并进入步骤S9；

S9、将最终的轴向应变图像除以最终的横向应变图像得到泊松比分布图，得到准确的轴向应变和横向应变估计。

2.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S1中：

压缩前、后的两帧射频RF回波信号是二维的，沿探头声波发射方向是轴向，垂直于轴向方向是横向；

压缩前、后的每一条射频RF回波信号是一维的，所述一维的压缩前、后的RF回波信号为：

r₁(x)＝s(x)p(x)+n₁(x) (1)

式中，r₁(x)、r₂(x)分别表示压缩前和压缩后RF回波信号；s(x)表示人体组织散射子函数，p(x)是点扩散函数；n₁(x)和n₂(x)表示系统白噪声；α是压缩系数，且α＝1-ε，ε是外加应变。

3.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S2中的最大迭代次数m为4。

4.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S3为：采用互相关方法根据压缩前、后的两帧数据沿轴向进行位移估计，得到一帧轴向应变估计图像；

其方法具体为：

S31、将压缩前、后相对应的两条轴向RF回波信号划分成若干个部分重叠的窗口；

上述若干个重叠的窗口的重叠率为75％；

S32、求每个窗口的位移估计；

S33、求每个窗口的应变估计；

上述每个窗口的应变估计是对应窗口位移估计的导数。

5.根据权利要求2所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S4中的时间处理为时间拉伸处理或时间压缩处理；

所述步骤S4具体为：

S41、将一帧轴向应变估计图像的一条轴向应变估计数据分成若干个不重叠的窗口；

S42、求出每个窗口的平均应变估计值；

S43、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的轴向RF回波信号；

所述步骤S43中的，轴向的时间处理后RF回波信号为：

r₃(x)＝r₂(αx)＝s(x)*p(αx)+n₃(x) (3)

式中，r₃(x)表示轴向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₃(x)表示系统白噪声。

6.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S5中的插值处理方法为三次样条插值方法。

7.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S6具体为使用子带划分方法和互相关方法对重建的横向RF回波信号中的压缩前、后的两帧RF回波信号数据进行横向位移估计；

其具体方法为：

S61、将重建的压缩前、后的每条轴向RF回波信号划分成三个子带信号；

S62、对每一对压缩前、后的子带信号，采用互相关方法沿横向进行位移估计得到三个子带横向应变估计图像；

S63、对三个子带横向应变估计图像进行加权平均，得到一幅横向应变估计图像。

8.根据权利要求2所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S7中的时间处理为时间拉伸处理或时间压缩处理，具体为：

S71、将一帧横向应变估计图像的一条横向应变估计数据分成若干个不重叠的窗口；

S72、求出每个窗口的平均应变估计值；

S73、通过每个窗口的平均应变估计值拉伸或压缩该窗口内的对应的横向RF回波信号；

所述步骤S73中的，横向的时间处理后RF回波信号为：

r4(x)＝r₂(αx)＝s(x)*p(αx)+n₄(x) (4)

式中，r₄(x)表示横向的时间处理后的RF回波信号，α为拉伸或压缩系数，n₄(x)表示系统白噪声。

9.根据权利要求1所述的二维迭代的超声弹性成像应变估计方法，其特征在于，所述步骤S9中泊松比为人体弹性组织的重要机械属性，泊松比ν为：

上式中，ε_l是横向应变，ε_a是轴向应变。