本发明涉及医学影像处理，具体涉及一种基于核磁共振图像的脑血管硬化检测方法。
背景技术：
1、脑卒中已成为世界范围内第二大致残或致命的疾病，其主要成因为脑动脉瘤、动静脉畸形等严重的脑血管疾病。根据世界卫生组织相关统计，全世界每年有1500万人中风，其中500万人死亡，另有500万人终身残疾。由此可见，针对脑血管疾病的早期诊断、安全治疗以及术后康复检查等方面的工作已刻不容缓。
2、近些年，随着医学影像技术与计算机辅助诊断技术的飞速发展，磁共振血管造影(magnetic resonance angiography，mra)等技术在脑血管疾病的临床诊断和治疗中被普遍采用。在临床诊断中，医生从图像中获取的脑血管中心线信息，可以在脑血管疾病的治疗上起到良好的协助作用。因此面向mra图像的脑血管中心线提取技术成为医学图像处理领域的核心问题之一。由于人们常常是对心脏冠脉血管和视网膜血管进行研究，然而对于脑血管的研究少之又少。
技术实现思路
1、本发明的目的在于提供一种基于核磁共振图像的脑血管硬化检测方法，通过分析处理三维时间飞跃磁共振血管图像，为医务人员诊断脑血管硬化提供参考。
2、为了实现上述目的，本发明采用的技术方案如下：
3、一种基于核磁共振图像的脑血管硬化检测方法，包括以下步骤：
4、s1：获取检测对象的脑血管图像，并对脑血管图像进行预处理；
5、s2：基于fmm-mrf脑血管分割模型，提取脑血管特征分割图；
6、s3：获取脑血管中心线；
7、s4：根据血管中心线计算得到血管直径分布图；
8、s5：计算血管狭窄程度，并根据脑血管狭窄测量标准评估脑血管硬化程度。
9、进一步的，所述s1中脑血管图像通过三维时间飞跃磁共振血管成像(3d-tof-mra)技术获取。
10、进一步的，所述s1中预处理具体包括以下步骤：
11、s11：对三维时间飞跃磁共振血管图像进行偏置场修正；
12、s12：去除脑血管图像中的颅骨结构；
13、s13：将去除颅骨的血管图像输入基于hessian矩阵特征值的多尺度滤波器，得到增强后的血管特征图。
14、进一步的，所述s2具体包括以下子步骤：
15、s21：采用k-means聚类方法获取有限混合模型的初始参数；
16、s22：采用最大期望算法对有限混合模型的初始参数进行更新优化，进而建立有限混合模型；
17、s23：将脑血管图像输入有限混合模型得到血管分割结果；
18、s24：血管分割结果作为马尔可夫随机场统计过程中的初始标记场，结合血管特征图进一步优化分割结果。
19、进一步的，所述s3具体包括以下子步骤：
20、s31：通过细化算法获取脑血管骨架线；
21、s32：提取脑血管骨架线血管分支点集；
22、s33：根据局部边缘距离场进行多分支点加权平均，将多个血管分支点合并为一个分支点，进而优化脑血管骨架线；
23、s34：通过优化后的脑血管骨架线转化为脑血管中心线。
24、进一步的，所述s4具体包括以下子步骤：
25、s41：提取脑血管中心线的每一点；
26、s42：采用基于mom评价模型的跟踪算法测量每一点的直径值。
27、进一步的，所述s42具体包括以下子步骤：
28、s421：计算感兴趣区域的每一中心点的血管方向；
29、s422：获取垂直于每一中心点的血管方向的左右边缘点；
30、s423：计算左右边缘点距中心点的距离，即中心点所在血管半径；
31、s424：根据血管半径构建脑血管直径分布图。
32、进一步的，所述s5具体包括：
33、s51：根据脑血管直径分布图获取颅内动脉管径和残存血管最小管径；
34、s52：按照下列公式计算狭窄程度：
35、脑血管狭窄程度＝(颅内动脉管径-残存血管最小管径)/颅内动脉管径×100％；
36、s53：参照北美症状性动脉内膜切除术研究协助组制定的脑血管狭窄测量标准，评估脑血管硬化程度。
37、本发明的有益效果：
38、(1)本发明通过对三维时间飞跃磁共振血管图像进行预处理，并且采用多尺度滤波器增强血管目标，能够能加精确的提取目标血管的中心线，同时抑制非血管结构的干扰，有效地去除部分噪声干扰；
39、(2)本发明通过颅骨剔除技术和脑血管特征分割来精简感兴趣区域；由于血管的颅内容积小，传统em算法对fmm参数估计过程是有偏的，由此提出基于血管知识的em算法；将多尺度血管增强结果与fmm中血管容积估计权值相结合，进一步得到具有概率意义的血管特征图，并将其嵌入到mrf过程，构建了双能量约束的mrf血管优化过程，在分割精度和鲁棒性上相对传统统计模型明显提高。
40、(3)本发明针对中心线提取与介入路径规划：在细化算法得到骨架线的基础上，结合边缘距离场解决了分支点不唯一问题，得到精确的脑血管中心线；
41、(4)本发明通过图像处理得到的脑血管中心线，进一步计算血管直径，为医务人员诊断脑血管硬化提供可靠数据支持。
1.一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s1中脑血管图像通过三维时间飞跃磁共振血管成像技术获取。
3.根据权利要求2所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s1中预处理具体包括以下步骤：
4.根据权利要求3所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s2具体包括以下子步骤：
5.根据权利要求1所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s3具体包括以下子步骤：
6.根据权利要求5所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s4具体包括以下子步骤：
7.根据权利要求6所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s42具体包括以下子步骤：
8.根据权利要求1所述的一种基于核磁共振图像的脑血管硬化检测方法，其特征在于，所述s5具体包括：