1.一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,该方法包括以下步骤:获取每个时刻下采集的脑血管造影图像;
根据脑血管造影图像中每个像素点与相邻的脑血管造影图像中对应的像素点之间的灰度差异,获取脑血管造影图像中每个像素点的正常程度;根据像素点周围的像素点,获取像素点的局部区域不同方向的灰度共生矩阵,结合像素点以及周围像素点的正常程度,获取灰度共生矩阵中每个非零元素的聚类权重;根据灰度共生矩阵中每个非零元素的位置坐标与聚类权重,对所有非零元素进行聚类,获取若干簇类;根据每个簇类在灰度共生矩阵中的分布位置,获取像素点的边缘特征;
将每个脑血管造影图像分成若干子图像;根据子图像中每个像素点的边缘特征,获取子图像中脑血管大致延伸方向;根据子图像中脑血管大致延伸方向,结合子图像中每个像素点的梯度值与梯度方向,获取子图像中每个像素点为脑血管边缘像素点的可能性;
根据像素点为脑血管边缘像素点的可能性,赋予每个像素点平滑参数并对脑血管造影图像进行去噪。
2.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取脑血管造影图像中每个像素点的正常程度,包括的具体方法为:对于任意脑血管造影图像,将与所述脑血管造影图像的采集时刻的相邻时刻下所采集的脑血管造影图像,记为所述脑血管造影图像的对比图像;
将所述脑血管造影图像中任意像素点记为目标像素点,并将所述脑血管造影图像的对比图像与目标像素点坐标位置相同的像素点,作为相关像素点;
以目标像素点为中心,构建一个局部窗口;将局部窗口内的像素点记为目标像素点的局部像素点;
获取相关像素点的局部像素点;
根据目标像素点的局部像素点与相关像素点的局部像素点之间的灰度差异,得到目标像素点的正常程度,所述目标像素点的正常程度与所述灰度差异呈负相关关系。
3.根据权利要求2所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取像素点的局部区域不同方向的灰度共生矩阵,包括的具体方法为:对于任意脑血管造影图像中任意像素点,将所述像素点的局部像素点组成的区域记为所述像素点的局部区域,获取所述像素点的局部区域0°、45°、90°、135°方向的灰度共生矩阵,并将每个灰度共生矩阵中135°方向的对角线,分别记为每个灰度共生矩阵的主对角线。
4.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取灰度共生矩阵中每个非零元素的聚类权重,包括的具体方法为:对于像素点的局部区域0°方向的灰度共生矩阵中任意非零元素,标记所述非零元素在像素点的局部区域对应的像素点记为标记像素点,将所有标记像素点正常程度的均值与所述非零元素的积作为所述非零元素的聚类权重;
获取像素点的局部区域90°方向的灰度共生矩阵中任意非零元素的聚类权重。
5.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述根据灰度共生矩阵中每个非零元素的位置坐标与聚类权重,对所有非零元素进行聚类,获取若干簇类,包括的具体方法为:在像素点的局部区域0°方向的灰度共生矩阵中随机设置3个初始聚类中心,对于任意非零元素,获取所述非零元素与每个初始聚类中心之间的加权距离,其权重为所述非零元素的聚类权重,将最小的加权距离对应的初始聚类中心,作为所述非零元素所归属的聚类中心;
得到每个非零元素所归属的聚类中心,将归属于相同聚类中心的非零元素归为同一初始簇类;
对于任意初始簇类,根据初始簇类中所有非零元素的在灰度共生矩阵中的位置坐标以及聚类权重,获取初始簇类的加权中心;
得到每个初始簇类中的加权中心,将每个初始簇类中的加权中心作为一个新的聚类中心,根据新的聚类中心与非零元素之间的加权距离,迭代获取加权中心,直至加权中心的位置不发生改变,得到像素点的局部区域0°方向的灰度共生矩阵中的3个簇类;
获取像素点的局部区域90°方向的灰度共生矩阵中的3个簇类。
6.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取像素点的边缘特征,包括的具体方法为:对于任意像素点的局部区域0°方向的灰度共生矩阵中的3个簇类,获取每个簇类的加权中心与灰度共生矩阵的主对角线之间的距离,将距离最大簇类记为边缘簇类,将其余簇类记为内部簇类;将两个内部簇类的加权中心到灰度共生矩阵的主对角线之间的距离和记为内部簇类距离,对边缘簇类的加权中心到灰度共生矩阵的主对角线之间的距离与内部簇类距离的比值进行归一化,将归一化结果作为像素点在0°方向上的边缘特征;
获取像素点在45°、90°、135°方向上的边缘特征。
7.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述将每个脑血管造影图像分成若干子图像,包括的具体方法为:预设一个 大小的子图像规格,将每个脑血管造影图像分成若干个规格为 的子图像。
8.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取子图像中脑血管大致延伸方向,包括的具体方法为:对于任意子图像,根据子图像中所有像素点在0°方向上的边缘特征,获取0°方向为子图像中脑血管大致延伸方向的可能性,其具体的计算公式为:式中, 表示0°方向为子图像中脑血管大致延伸方向的可能性;表示子图像中像素点的数量;表示子图像中所有像素点在0°方向上的边缘特征的均值; 表示子图像中第 个像素点在0°方向上的边缘特征;表示取绝对值函数; 表示以自然常数为底数的指数函数;
获取45°、90°、135°方向为子图像中脑血管大致延伸方向的可能性,将0°、45°、90°、
135°方向为子图像中脑血管大致延伸方向的可能性中,最大值对应的方向作为子图像中脑血管大致延伸方向。
9.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述获取子图像中每个像素点为脑血管边缘像素点的可能性,包括的具体方法为:对于任意子图像中的任意像素点,根据所述子图像中脑血管大致延伸方向以及所述像素点的梯度值与梯度方向,获取像素点为脑血管边缘像素点的可能性,其具体的计算公式为:式中, 表示所述子图像中第 个像素点为脑血管边缘像素点的可能性; 表示所述子图像中第 个像素点的梯度值; 表示脑血管造影图像中最大的梯度值; 表示所述子图像中第 个像素点的梯度方向与所述子图像中脑血管大致延伸方向之间的夹角;
表示取正弦函数。
10.根据权利要求1所述一种面向脑血管造影x射线检查确定血液动力的方法,其特征在于,所述赋予每个像素点平滑参数并对脑血管造影图像进行去噪,包括的具体方法为:预设一个初始平滑参数,对于任意像素点,将所述像素点为脑血管边缘像素点的可能性的反比例归一化结果,与初始平滑参数的积,作为所述像素点的平滑参数;
获取每个脑血管造影图像中所有像素点的平滑系数,利用非局部均值滤波算法,结合每个脑血管造影图像中所有像素点的平滑参数,对每个脑血管造影图像进行去噪。