[发明专利]超声多普勒血流成像的处理方法及处理系统

权利要求书

1.一种超声多普勒血流成像的处理方法，其特征在于，所述方法包括以下步骤：

S1、分别获取每个扫查包下、每条扫查线上、每个采样点对应的回波RF信号；分别对获取的回波RF信号做希尔伯特变换以形成正交IQ复信号；

S2、在慢时间方向上，分别对获得的正交IQ复信号做滤波处理，以获得对应的滤波变换信号；

S3、在快时间方向，分别对获得的滤波变换信号做短时傅里叶变换，使其结果形成频域矩阵；

S4、获取当前所述频域矩阵中，每个快时间方向上的频率分量；

S5、根据一维自相关算法，估计当前频域矩阵中每列数据在其对应频率分量上的平均速度和能量；

S6、对当前频域矩阵中对应频率分量的平均速度和能量进行线性拟合，获取用于最终的频谱显示的最终能量和最终速度；

其中，所述步骤S3具体包括：P1、在同一深度下，对应每个扫查包、在每条扫查线的快时间方向上分别获取L个采样点对应的滤波变换信号，L小于等于N，N表示每条扫查线对应的原始采样点数量；P2、对每条扫查线对应的L个滤波变换信号序列进行加窗补零后，做L_fft阶短时傅里叶变换，使其结果形成对应一个扫查包的M*L_fft频域矩阵，M表示每个扫查包设置的扫查线数量；L_fft表示短时傅里叶变换的阶数；

所述步骤S4具体包括：M1、获取当前深度下，M*L_fft频域矩阵对应的中心频率位置；M2、以获取的中心频率位置为中心，获取当前深度下，感兴趣带宽的截止频率位置；M3、根据获取的感兴趣带宽的截止频率位置获取其对应的M*K频域矩阵，K表示感兴趣带宽的截止频率所对应的列的总和；M4、获取M*K频域矩阵中，每列数据对应的频率分量。

2.根据权利要求1所述的超声多普勒血流成像的处理方法，其特征在于，

所述步骤S5包括：根据一维自相关算法，估计M*K频域矩阵中每列数据在频率分量f(k)上的平均速度和能量，并将其集合标记为第一平均速度序列，和第一能量序列；

所述步骤S5和步骤S6之间，所述方法还包括：

N1、根据预设的高速度阈值，低速度阈值，噪声能量阈值和组织能量阈值对M*K频域矩阵对应的第一平均速度序列和第一能量序列进行第一次预处理；

若相同频率分量对应的第一平均速度序列中的平均速度值大于所述高速度阈值，对应能量序列中的能量值小于噪声能量阈值；

或相同频率分量对应的第一平均速度序列中的平均速度值小于所述低速度阈值，对应能量序列中的能量值大于组织能量阈值；

则剔除相同频率分量对应第一平均速度序列和第一能量序列对应的中的平均速度值和能量值，形成第二平均速度序列和第二能量序列。

3.根据权利要求2所述的超声多普勒血流成像的处理方法，其特征在于，所述步骤N1后，所述方法还包括：

N2、获取第二平均速度序列对应的方差值及速度平均值，并根据其对第二平均速度序列和第二能量序列进行第二次预处理；

若第二平均速度序列中的任一平均速度值与第二平均速度序列的速度平均值差值的平方大于第二平均速度的方差值与第一预设常数值的乘积，则剔除相同频率分量对应第二平均速度序列和第二能量序列中的平均速度值和能量值，形成第三平均速度序列和第三能量序列。