[发明专利]CT图像的运动向量场的计算方法及装置

权利要求书

1.一种CT图像的运动向量场的计算方法，其特征在于，该方法包括：

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

计算所述目标运动区域的运动水平参数，计算各个时间控制点对应的运动向量场的数值，以及，确定运动向量场的方向；

依据所述运动水平参数、运动向量场的数值和方向，计算所述初始CT重建图像对应的目标运动向量场。

2.根据权利要求1所述的方法，其特征在于，在确定所述初始CT重建图像对应的目标运动向量场之前，还包括：计算所述运动向量场的扩展参数，包括：

确定若干间隔相等的初始扩展参数；

依据各个初始扩展参数计算对应的各个边界梯度标准差；

依据所述各个边界梯度标准差拟合出标准差变换曲线；

选择所述标准差变换曲线中最大的边界梯度标准差；

将所述选择的最大的边界梯度标准差对应的初始扩展参数确定为目标扩展参数；

则在确定所述初始CT重建图像对应的目标运动向量场之后，还包括：

将所述目标运动向量场与所述目标扩展参数的乘积确定为最终的补偿运动向量场。

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

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

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

4.根据权利要求1所述的方法，其特征在于，所述计算各个时间控制点对应的运动向量场的数值，包括：

针对各个时间控制点中的任一当前时间控制点，计算该当前时间控制点对应的投影数据相对于运动向量场的导数；

将所述运动水平参数与所述导数的乘积作为当前时间控制点对应的运动向量场的当前数值；

判断所述当前时间控制点是否为最后一个时间控制点，如果否，则选取下一个时间控制点作为当前时间控制点，并执行所述计算该当前时间控制点对应的投影数据相对于运动向量场的导数的步骤，直至所有的时间控制点对应的运动向量场的数值都计算完毕。

5.根据权利要求1所述的方法，其特征在于，所述确定运动向量场的方向，包括：

确定所述运动向量场的初始运动方向参数包含的若干初始方向；

选取若干初始方向中的任一方向作为当前方向，计算所述当前方向下的当前运动向量场；

依据所述当前运动向量场对所述目标运动区域进行补偿以重建出当前目标运动图像；

计算所述当前目标运动图像的当前边界梯度标准差；

判断所述当前方向是否为最后一个方向，如果是，则确定若干初始方向分别对应的若干边界梯度标准差中最大的边界梯度标准差，并将所述最大的边界梯度标准差对应的方向确定为所述运动向量场的目标方向；

如果否，则选择下一个方向作为当前方向，执行所述计算所述当前方向下的当前运动向量场的步骤。

6.一种CT图像的运动向量场的计算装置，其特征在于，包括：

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

计算运动水平参数模块，用于计算所述目标运动区域的运动水平参数；

计算数值模块，用于计算各个时间控制点对应的运动向量场的数值；

确定方向模块，用于确定运动向量场的方向；

计算向量场模块，用于依据所述运动水平参数、运动向量场的数值和方向，确定所述初始CT重建图像对应的目标运动向量场。