[发明专利]声辐射力脉冲成像估算方法和系统

权利要求书

1.一种声辐射力脉冲成像估算方法，包括：

数据划分步骤，读取通过声辐射力脉冲成像所采集的输入数据，并将所述输入数据分成预定数量数据矩阵；

数据滤波变换步骤，分别对每个数据矩阵进行前向滤波处理和逆向滤波处理，并对滤波处理后的数据矩阵进行希尔伯特变换得到解析数据；

位移求取步骤，对所述解析数据进行互相关运算，求取得到预定数量位移矩阵；

位移合并步骤，对每个位移矩形进行时间累加，取累积和最大的列位移，再将预定数量位移矩阵抽取的列位移合并，生成新的位移矩阵；

插值步骤，对得到的新的位移矩阵进行时间方向上的三次样条插值；

拟合步骤，对插值结果进行拉东变换，得到拉东变换矩阵；

估算步骤，获取所述声辐射力脉冲成像的参数，根据所述参数及拉东变换矩阵求取剪切波速率。

2.根据权利要求1所述的声辐射力脉冲成像估算方法，其特征在于，所述数据滤波变换步骤中，采用多个流按顺序依次对预定数量数据矩阵进行前向滤波处理、逆向滤波处理和希尔伯特变换处理。

3.根据权利要求1所述的声辐射力脉冲成像估算方法，其特征在于，所述数据滤波变换步骤中前向滤波处理的公式为：

当i≤iParalen时：

tOutPut(k,i)=Σj=0i(param[j]*input(k,i-j))+Σj=i+1iParalen-2(param[j]*input(k,0))]]>

当i＞iParalen时：

tOutPut(k,i)=Σj=0iParalen-2(param[j]*input(k,i-j))]]>

其中，param是滤波器，iParalen是滤波器长度，input为待滤波矩阵，k、i分别表示待滤波矩阵的第k行和第i列，tOutPut是滤波结果矩阵；

所述数据滤波变换步骤中逆向滤波处理的公式为：

当i≥iParalen时：

output(k,Width-1-i)=Σj=0iParalen-1(Param[j]*input(k,Width-1+j-i))]]>

其中，input为待滤波矩阵，k、i分别表示待滤波矩阵的第k行和第i列，output为逆向滤波的滤波结果矩阵，Width为滤波器宽度；

所述前向滤波处理的内核函数和逆向滤波处理均采用图形处理器并行处理，所述图形处理器调用的内核函数的blocks数量和threads数量相同，所述blocks数量为待滤波矩阵的行数，所述threads数量为待滤波矩阵的列数。