[发明专利]一种血管路径确定方法、装置、终端及存储介质

权利要求书

1.一种血管路径确定方法，其特征在于，所述方法包括：

步骤a，获取电子计算机断层扫描数字X线摄影DR图像中血管路径的起始点和结束点，将所述起始点设为第一预定点；

步骤b，以所述第一预定点为球心，根据预设球半径得到感兴趣区域，在位于所述感兴趣区域边缘的像素点中，将势能函数取值最小的像素点确定为第二预定点，所述势能函数用于确定所述感兴趣区域边缘中任意一个像素点的势能；

步骤c，获取所述第二预定点在所述血管路径上的切线方向，将所述切线方向代入预设的最小活动图函数后，得到所述第一预定点至所述第二预定点的最小活动图函数，对所述第一预定点至所述第二预定点的最小活动图函数求解，得到所述第一预定点至所述第二预定点的最小路径；

步骤d，如果所述第二预定点不是所述结束点，则将所述第二预定点设置为第一预定点，转至执行步骤b；如果所述第二预定点是所述结束点，则将所有最小活动图函数求解得到的最小路径合并为所述血管路径。

2.根据权利要求1所述的方法，其特征在于，所述获取所述第二预定点在所述血管路径上的切线方向，包括：

确定与所述第二预定点之间的距离小于预设距离且灰度值小于预设阈值的第三预定点和第四预定点，根据所述第三预定点和所述第四预定点，确定出所述第二预定点在所述血管路径上的切线方向。

3.根据权利要求1所述的方法，其特征在于，所述将所述切线方向代入预设的最小活动图函数，包括：

将所述切线方向代入预设的势能函数，得到所述第一预定点至所述第二预定点的势能函数；

将所述第一预定点至所述第二预定点的势能函数代入所述预设的最小活动图函数，得到所述第一预定点至所述第二预定点的最小活动图函数；

其中，所述预设的势能函数为：

所述预设的最小活动图函数为：

所述第一预定点至所述第二预定点的最小活动图函数为：

其中，所述grey为所述感兴趣区域内所有像素点的灰度值，所述mean_grey为感兴趣区域内所有像素点的灰度平均值，所述λ为常数，所述β为常数，所述T₁为所述第二预定点在所述血管路径上的切线方向，所述x₁为所述第一预定点，所述x₂为所述第二预定点，所述x_k为第k预定点，所述n为大于1的自然数，所述Ω为所述P(x)以所述x₁为下限、以所述x₂为上限的积分区间，所述p为所述血管路径上的任意点，所述x为所述Ω内的任意一点，所述θ为所述切线方向的单位向量与所述x和所述x_k连线的单位向量的夹角，所述s为所述第一预定点与所述第二预定点之间路径的弧长参数，所述P(x)为所述预设的势能函数，所述为所述第一预定点与所述第二预定点之间的路径，所述为从所述第一预定点到所述第二预定点延路径积分得到的最小能量。

4.根据权利要求1所述的方法，其特征在于，所述对所述第一预定点至所述第二预定点的最小活动图函数求解，得到所述第一预定点至所述第二预定点的最小路径，包括：

根据最速下降方法，对所述最小活动图函数求解，得到所述第一预定点至所述第二预定点的最小路径。

5.根据权利要求1-4中任一所述的方法，其特征在于，所述获取电子计算机断层扫描数字X线摄影DR图像中血管路径的起始点和结束点之后，所述方法还包括：

以所述起始点为球心，根据预设球半径得到预定区域；在所述预定区域中确定出灰度值最低的点，将所述灰度值最低的点确定为起始点；

以所述结束点为球心，根据预设球半径得到预定区域；在所述预定区域中确定出灰度值最低的点，将所述灰度值最低的点确定为结束点。