[发明专利]一种弹性减影成像方法

权利要求书

1.一种弹性减影成像方法，其特征在于，该方法包括以下步骤：

1)对损伤引入前的组织施加一定的压缩，并采集压缩前后的超声射频数据；

2)对损伤引入后的组织施加一定的压缩，并采集压缩前后的超声射频数据；

3)根据步骤1)所采集的超声射频数据计算损伤引入前的组织的位移场；

4)根据步骤2)所采集的超声射频数据计算损伤引入后的组织的位移场；

5)根据步骤3)、4)得到的结果，计算弹性减影图像。

2.如权利要求1所述方法，其特征在于，所述步骤5)计算弹性减影图像，具体步骤如下：

5-1)通过对损伤引入前的组织的位移在深度方向求微分来计算第k条扫描线上第i个数据段对应组织的应变ε1_i，k，如下所示：

ϵ1i,k=ϵ1(z,k)|z=(M+i)·Δz,k=∂d1(z,k)∂z|z=(M+i)·Δz,k]]>

其中ε1(z，k)、d₁(z，k)分别是损伤引入前的组织第k条扫描线上中心深度为z的数据段对应组织的应变和位移，Δz为相邻两点数据对应的实际距离，M为与扫描线上数据段长度相关的常数(2M+1为数据段长度)；

以此获得损伤引入前的每条扫描线上每个数据段对应组织的应变所构成的组织应变分布S1：

S₁＝{ε1_i，k}|_{i＝1、2、...、(L-2M)，k＝1、2、...、N}

其中L为每条扫描线包含的数据点数，N为采集的超声射频数据包含的扫描线数目；

5-2)用上述步骤5-1)计算应变的方法，计算损伤引入后的每条扫描线上每个数据段对应组织的应变所构成的组织应变分布S2：

S₂＝{ε2_i，k}|_{i＝1、2、...、(L-2M)，k＝1、2、...、N}

其中ε2_i，k是损伤引入后的组织第k条扫描线上第i个数据段对应组织的应变；

5-3)根据步骤5-1)和5-2)的结果，将损伤引入后的组织应变分布S2减去损伤引入前的组织应变分布S1即得到弹性剪影成像S_sub，即：

S_sub＝S₂-S₁＝{ε2_i，k-ε1_i，k}|_{i＝1、2、...、(L-2M)，k＝1、2、...、N}。

3.如权利要求1所述方法；其特征在于，所述步骤5)计算弹性减影图像，具体步骤如下：

5-1)对步骤4)计算的损伤引入后的组织的位移场进行如下规整化处理：

d2′(z,k)=d2(z,k)·d1((L-M)·Δz,k)d2((L-M)·Δz,k)|z=(M+1)·Δz,(M+2)·Δz,···,(L-M)·Δz,k=1,2,...,N]]>

其中d₂(z，k)、d′₂(z，k)分别是损伤引入后的组织第k条扫描线上中心深度为z的数据段对应组织的位移及对该位移规整化处理后的位移，d₁((L-M)●Δz，k)、d₂((L-M)●Δz，k)分别为损伤引入前后第k条扫描线上最后一个数据段对应组织的位移，L为每条扫描线包含的数据点数，Δz为相邻两点数据对应的实际距离，M为与扫描线上数据段长度相关的常数(2M+1为数据段长度)，N为采集的超声射频数据包含的扫描线数目；从而得到新的规整化后的位移场：

{d′₂(z，k)}|_{z＝(M+1)●Δz、(M+2)●Δz、…、(L M)●Δz，k＝1、2、...、N}；

5-2)利用规整化后的位移场和步骤3)计算的损伤引入前的位移场进行减影操作，得到减影位移场D_sub：

D_sub＝{d′₂(z，k)-d₁(z，k)}|_{z＝(M+1)●Δz、(M+2)●Δz、…、(L M)●Δz，k＝1、2、...、N}其中d₁(z，k)是损伤引入前的组织第k条扫描线上中心深度为z的数据段对应组织的位移；

5-3)对得到的减影位移场D_sub沿着组织深度的方向微分得到弹性减影成像S_sub：

Ssub=∂Dsub∂z]]>

={∂(d2′(z,k)-d1(z,k))∂z}|z=(M+1)·Δz,(M+2)·Δz,···,(L-M)·Δz,k=1,2,...,N.]]>