[发明专利]一种基于多角度序列化图像空间特征点集的脑血管图像分割方法

权利要求书

1.一种基于多角度序列化图像空间特征点集的脑血管图像分割方法，其特征在于：包括以下步骤：

S1：对脑血管的活片图像和蒙片图像进行配准：在活片图像和蒙片图像上确定几何特征点，找出匹配的特征点对，通过匹配的特征点对之间的关系对蒙片图像进行处理获得序列化减影图像；

S2：提取序列化减影图像的几何特征点集；

S3：对序列化减影图像中血管边缘的特征点进行局部位置调整：通过灰度值梯度来调整特征点的位置，将血管边缘的特征点移动到血管内；

S4：采用空间旋转坐标系模型剔除减影图像中非血管位置的特征点：建立一个空间旋转坐标系模型，利用相邻减影图像上特征点的位置信息，求出特征点在其它减影图像中的理论位置，根据其他减影图像中特征点的位置，采取投票的方式剔除非血管位置的特征点；

S5：确定序列化减影图像的自适应分割阈值；

S6：基于区域生长和自适应阈值的血管图像分割：将S4中获得的特征点作为种子点，将S5中获得的自适应分割阈值作为生长法则的标准，采用区域生长算法分割序列化减影图像获得纯净的脑血管图像；

S4中具体采用如下方式：

S41：从当前的减影图像中选取一特征点P₁(x₁,y₁)，在当前的减影图像中的相邻减影图像中计算与P₁满足如下关系的特征点P₂(x₂,y₂)：

|x₂-x₁|≤4

|y₂-y₁|≤2

S42：根据P₁的横坐标x₁、P₂的横坐标x₂以及P₁和P₂之间的角度α＝2°计算特征点P₁到旋转中轴的距离r：

r=x22-2x1x2cosα+x12sinα]]>

S43：计算特征点P₁在其他减影图像中投影的理论位置P_n：设当前处理的减影图像A的序号为i₁，其他减影图像B的序号为i_n，通过序号的差值求出两幅减影图像之间的旋转的角度α：

α＝|2×(i_n-i₁)|；

S44：根据上述计算出的r和旋转角度α求出特征点P₁在当前减影图像B中的投影的位置p_n(x_n,y_n)；

x_n＝x₁±r×cosα

y_n＝y₁

S45：从减影图像B中查找是否存在满足以下要求的特征点(x_m,y_m)，如果存在则给特征点P₁投一票，否则不投票；

|x_m-x_n|≤4

|y_m-y_n|≤2

S46：设定阈值，剔除得票数低于该阈值的特征点后获得特征点集。