[发明专利]用于计算CT扫描的重建图像中的运动向量场的方法及装置

权利要求书

1.一种用于计算CT扫描的重建图像中的运动向量场的方法，包括：

依据针对CT扫描的初始重建图像中的血管中心线位置，确定需计算运动向量场的目标区域和时间控制点；

计算所述目标区域的运动水平参数；

计算各所述时间控制点对应的运动向量场数值；

确定所述目标区域的运动向量场方向；

依据所述运动水平参数、所述运动向量场方向和各所述时间控制点对应的运动向量场数值，计算所述目标区域的运动向量场。

2.根据权利要求1所述的方法，其中，

在确定所述目标区域的运动向量场方向之后，还包括：确定用于调整所述运动向量场的大小的扩展参数；以及

计算所述目标区域的运动向量场，包括：依据所述运动水平参数、所述运动向量场方向、所述扩展参数和各所述时间控制点对应的运动向量场数值，计算所述目标区域的运动向量场。

3.根据权利要求2所述的方法，其中，依据所述运动水平参数、所述运动向量场方向、所述扩展参数和各所述时间控制点对应的运动向量场数值，计算所述目标区域的运动向量场的公式为：

其中，为时间控制点t_n对应的运动向量场数值，

为时间控制点t_n对应的运动向量，为t_n时刻运动向量的X方向数值，为t_n时刻运动向量的Y方向数值，所述目标区域的运动向量场包括各个所述时间控制点对应的运动向量，

ω为所述扩展参数，

η为所述运动向量场方向，

为所述运动水平参数。

4.根据权利要求2所述的方法，其中，确定用于调整所述运动向量场的大小的扩展参数包括：

依据所述运动水平参数、所述运动向量场方向和各所述时间控制点对应的运动向量场数值，计算初始运动向量场；

设定若干间隔相等的初始扩展参数，并计算利用各所述初始扩展参数对所述初始运动向量场进行扩展后的扩展运动向量场；

利用各所述扩展运动向量场对所述目标区域进行运动补偿，以重建出各所述初始扩展参数对应的扩展图像，并计算各所述扩展图像的边界梯度标准差；

依据各所述初始扩展参数以及与各所述初始扩展参数对应的边界梯度标准差拟合出标准差变化曲线；

确定所述标准差变化曲线中的边界梯度标准差最大值，并将所述边界梯度标准差最大值对应的初始扩展参数确定为所述扩展参数。

5.根据权利要求1所述的方法，其中，计算所述目标区域的运动水平参数，包括：

计算所述目标区域的边界梯度标准差；

依据所述边界梯度标准差和预设的阈值确定所述运动水平参数。

6.根据权利要求5所述的方法，其中，所述计算所述目标区域的边界梯度标准差，包括：

根据标准差的定义，利用边界梯度绝对值计算所述目标区域的边界梯度标准差，其中，

所述边界梯度绝对值的公式为：

其中，E(a,b)为像素点(a,b)处的边界梯度绝对值，

Δx为X方向上的像素点采样间隔，Δy为Y方向上的像素点采样间隔；

a表示该像素点在X方向上的索引，b为该像素点在Y方向上的索引；

f(a,b)为像素点(a,b)在所述目标区域的重建图像上的数值。

7.根据权利要求6所述的方法，其中，依据所述边界梯度标准差和预设的阈值确定所述运动水平参数的公式为：

其中，为所述预设的阈值，StdE为所述边界梯度标准差，为所述运动水平参数。