[发明专利]脑血管形态特征定量分析方法

说明书

本发明属于医学影像技术领域，具体为一种脑血管形态特征定量分析方法。本发明利用血管网络的结构特点，通过构建血管网络邻接矩阵实现各项定量特征的计算，具有很强的普适性，可用于X射线、CT、MRI等不同采集方式获取的血管影像。分析结果包含血管分支数目、分支长度、血管半径、血管弯曲度、血管网络复杂度等多种特征，能够有效提高医生分析血管定量参数的效率；本发明具有旋转不变性，血管影像经过旋转后不影响分析结果，可实现二维血管影像以及三维血管影像的自动分析。本发明在心血管疾病以及未来其他重大血管疾病的智能诊治中发挥巨大的作用，极具市场潜力和经济社会效益。

技术领域

本发明属于医学影像技术领域，具体涉及一种脑血管形态特征定量分析方法。

背景技术

众多疾病的产生与发展往往伴随着组织器官内血管的异常和病变。随着医学影像技术的飞速发展，利用各种成像方法无创地观察组织血管，分析和研究血管的变化情况，在临床中发挥重要作用。

血管异常与病变是心脑血管疾病、肝硬化、各色肿瘤等多种疾病主要特征。长期以来，组织血管定量分析的方法是利用组织切片，测量血管密度、血管长度和管径等相关参数。此类方法属有创的方法，在实施过程中，需重复进行血管染色，借助光学显微镜进行观测。测量结果受染色结果的影响较大，同时光学显微镜的观测范围十分有限，影响测量结果的准确性。随着成像技术的不断发展，采用X射线，CT，PET，MR等成像技术，能够实现无创的获取组织血管影像，进一步可利用算法进行血管定量分析。

目前临床医生利用医学影像分析血管分布和形态变化，研究相关疾病的预防、发现和诊治方法，往往只能凭借主观的经验观察，对疾病的严重程度、治疗效果等的分析过分依赖医生本身的主观判断。同时，伴随着成像技术的进步和完善，影像所包含的信息更为丰富与复杂，加重了医生观测与分析数据的负担，进而影响诊断的准确性和治疗效果。因此，对血管进行自动的定量分析具有十分重要的实用价值。现有的血管定量分析方法能够实现简单的分支、半径计算，但缺少其他定量方法，难以用于肿瘤鉴别。

本发明利用CT、MR等影像，自动提取脑血管的定量特征，包括血管密度类特征，血管弯曲度特征和血管网络特征等。利用这些血管定量特征可用于识别、评估脑血管疾病，如脑卒中、脑动脉瘤、脑动脉狭窄等。

发明内容

针对现有技术存在的不足，本发明旨在提供一种自动脑血管形态特征定量分析方法。该方法适用于不同采集方式获取的血管图像，能够自动计算血管的分支、半径、弯曲度等特征，可便于临床分析与诊治。

本发明提供的脑血管形态特征定量分析方法，具体步骤如下：

步骤S1，从不同的采集方式中，获取组织影像的血管结构；

步骤S2，利用形态学细化算法(Kollmannsberger Pet al.The small world ofosteocytes:connectomics of the lacuno-canalicular network in bone.New Journalof Physics.2017)，对步骤S1中的血管结构进行细化，获取血管结构的中心线，中心线所有点的集合记为C；

步骤S3，判断血管中心线上点之间的邻接关系，构建中心线上所有点的邻接矩阵，记为A；

步骤S4，利用邻接矩阵计算中心线上所有点的度，根据度的大小将所有点分为结点与边结点：度为2的点是边结点，记为E，其他点为结点，记为N；

步骤S5，利用遍历搜索方法，构建结点的邻接矩阵，记为A_node，并得到中心线各分支的长度与分支个数；

步骤S6，计算所有边结点所处位置的血管半径；

步骤S7，计算所有边结点的弯曲程度，求出每条分支的平均弯曲程度；

步骤S8，利用邻接矩阵A，计算血管网络的聚类系数、网络结构熵。