基于CT影像的血管三维分割方法技术

基于CT影像的血管三维分割方法，包括：1)读取CT造影图像数据，把血管开始出现的层记为第一层，并生成第一层种子点，方法为：获取第一层灰度值的平均值；把本层图像的中心点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值的点生长为种子点；2)从第二层往下逐层生长，方法为：分别计算上一层种子点和非种子点的灰度值的平均值，对它们加权计算得到本层种子点灰度最小值Mi；从上一层的种子点向本层生长；按上述方法逐层往下生长，直到没有新的种子点生成为止。所有种子点的集合组成了血管三维分割图像。其目的在于提供一种减少对操作人员经验的依赖、提高分割准确度和完整性的基于CT影像的血管三维分割方法。

【技术实现步骤摘要】
基于CT影像的血管三维分割方法
本专利技术涉及医学
图像处理方法，具体的说是一种基于CT影像的血管三维分割方法。
技术介绍
过去的近三十年，医学成像技术得到了飞速发展。新的成像技术能够得到人体结构、组织功能在空间和时间域上的图像。医学成像让研究人员和医生在非接触的条件下了解病人潜在的疾病信息，已经成为手术导航和模拟，以及跟踪疾病变化的重要手段。随着医学成像在疾病诊断和治疗中的广泛应用，如何通过各种成像手段提取人体解剖结构中对临床有用的信息成为一个重要的问题。这些成像手段包括X-ray、CT(ComputedTomography)、MR(MagneticResonance)和PET(PositronEmissionTomographyimaging)等。在成像技术发展的同时，不断发展的图像分析工具提供了大量从图像上提取用于图像理解的重要特征的方法。然而，目前利用计算机准确有效的量化和分析人体内部结构和功能信息的技术仍然处在十分有限的水平。医学图像特殊的成像过程和图像性质使其需要特定的图像处理方法。在对血管进行分割时，比较常见的是区域生长法，区域生长法根据图像中血管所占领区域内的性质,如灰度、纹理、色彩等,把血管所占领的区域找出来。它有两个重要的准则：性质的相似性和空间的相邻性,假定空间上相邻、性质上相似的像素属于同一物体。以灰度为例,在血管图像内部用某一准则找一个初选点P,从P点出发按八邻域搜寻,设定一个值H,凡是与P点的灰度差不超过H的点,都认为在物体上。用此原则向四周搜索,直到找不到这样的点为止。区域生长法是一种半自动的方法,它需要预先给出一个或多个种子点和判断准则(比如灰度值范围)，参见图1，区域生长法包括如下步骤：(1)选择种子点，并设置生长区域内相邻点的灰度差范围，把所有的种子点放在一个待处理队列中，比如堆栈。(2)从队列中取出种子点，得到种子点相邻三维空间中的8个点，比较它们与种子点的灰度差，如果在灰度范围内，则该点视为生长区域内的点，并将该点放入队列；否则视为边缘。(3)重复步骤(2)直到遍历完队列中所有的点。所有生长出来的点的集合即是分割出来的图像。该方法需要人工辅助，种子点的选择和灰度范围设定都完全依赖操作人员的经验，而且在CT造影图像中，血管不同部位的像素灰度值是不同的，不同层之间存在灰度的渐变，需要按层进行灰度范围调整。传统区域生长法不能按层进行调整，难以适应真实血管的灰度变化，可能把一些别的组织生长进来，而真正的血管却不能生长进来。
技术实现思路
本专利技术的目的在于提供一种减少对操作人员经验的依赖、提高分割准确度和完整性，让所获得的血管三维更加完整、清晰、准确，对图像的处理效率更高的基于CT影像的血管三维分割方法。本专利技术基于CT影像的血管三维分割方法，包括：1)读取血管CT造影图像数据，把血管开始出现的层命名为第一层，并在该层生成种子点，所述种子点生成方法如下：101)获取所述第一层灰度值的平均值G1；102)从所述第一层中选取一个图像点，比较所述图像点的灰度值与平均值G1的大小，如果大于等于平均值G1，把所述图像点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；如果所述图像点的灰度值小于平均值G1，把本层距离图像点最近且灰度值大于或等于平均值的点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；2)与所述第一层相邻且位于第一层下方的层命名为第二层，从所述第二层往下逐层进行生长，方法为：201)获取上一层种子点的灰度值的平均值Ai-1，和上一层非种子点的灰度值的平均值Bi-1，对平均值Ai-1和平均值Bi-1加权计算得到本层的种子点灰度最小值Mi；202)从上一层的种子点向本层进行生长，与种子点相邻且灰度值大于种子点灰度最小值Mi的点都生长为种子点；按步骤201和202所述方法逐层往下进行种子点的生长，直到没有新的种子点生成为止。所有生长出来的种子点的集合即组成了血管三维分割图像。优选地，所述步骤202还包括：获取本层种子点数量Ni，并与上一层的种子点数量Ni-1进行比较，如果Ni和Ni-1接近，即|Ni-1-Ni|/Ni-1≤10％，则本层种子点生长完成；否则，用最大类间方差法计算本层的阈值T，比较种子点灰度最小值Mi和阈值T，如果|Mi-T|≤S，则本层种子点生长完成，其中所述S为取值范围为[10,30]的正整数；否则，对Mi进行调整，并重新执行步骤202，所述对Mi进行调整的方法如下：如果Ni大于Ni-1，如果Ni小于Ni-1，按调整后的Mi重新执行步骤202。优选地，本专利技术的方法还包括：3)设含有生成种子点的最下面那层为倒数第一层，在从上往下逐层种子点生长完成后，从倒数第二层开始逐层往上进行种子点回溯生长，方法为：301)获取下一层种子点的灰度值的平均值Ci+1，和下一层非种子点的灰度值的平均值Di+1，对平均值Ci+1和平均值Di+1加权计算得到本层的种子点灰度最小值Mi；302)从下一层的种子点向本层进行生长，与种子点相邻且灰度值大于种子点灰度最小值Mi的点都生长为种子点；按上述方法逐层往上进行种子点的生长，直到所述第一层的种子点生长完成。优选地，所述步骤302还包括：获取本层种子点数量Ni，并与下一层的种子点数量Ni+1进行比较，如果Ni和Ni+1接近，即|Ni+1-Ni|/Ni+1≤10％，则本层种子点生长完成；否则，用最大类间方差法计算本层的阈值T，比较种子点灰度最小值Mi和阈值T，如果|Mi-T|≤S，则本层种子点生长完成，其中所述S为取值范围为[10,30]的正整数；否则，对Mi进行调整，并重新执行步骤302，所述对Mi进行调整的方法如下：如果Ni大于Ni+1，如果Ni小于Ni+1，按调整后的Mi重新执行步骤302。优选地，设从上往下种子点生长和从下往上种子点回溯生长为一轮种子点生长，每一轮种子点生长完成后统计所有生成的种子点数量，并与上一轮种子点生长完成后统计生成的所有种子点数量比较，若种子点数量增多了，则继续进行下一轮种子点生长；若种子点数量没有增多，则所有生长出来的种子点的集合组成了血管三维分割图像。优选地，步骤102中所述图像点位于该层的中心区域。优选地，步骤102中如果所述图像点的灰度值小于平均值G1，还可提示人工介入指定第一个种子点。优选地，所述血管CT造影图像数据为DICOM格式。本专利技术基于CT影像的血管三维分割方法，自动生成种子点，并自动计算种子点灰度范围，减少了对操作人员经验的依赖；另外，根据CT影像中血管不同部位的像素灰度值不同、不同层之间存在灰度的渐变的特点，实现了按临层的灰度统计动态调整生成种子点的灰度范围，从而提高血管分割的准确度和完整性，让所获得的血管三维更加完整、清晰、准确，对图像的处理效率更高。下面结合附图对本专利技术基于CT影像的血管三维分割方法作进一步说明。附图说明图1为
技术介绍
中传统区域生长法的流程图；图2为本专利技术的基于CT影像的血管三维分割方法的一种实施方式流程图；图3为一种血管空间结构示意图；图4为一种复杂的血管空间结构示意图。具体实施方式本专利技术基于CT影像的血管三维分割方法，包括：1)读取血管CT造影图像数据，图像数据格式为DI本文档来自技高网...

【技术保护点】
基于CT影像的血管三维分割方法，其特征在于，包括：1)读取血管CT造影图像数据，把血管开始出现的层命名为第一层，并在该层生成种子点，所述种子点生成方法如下：101)获取所述第一层灰度值的平均值G1；102)从所述第一层中选取一个图像点，比较所述图像点的灰度值与平均值G1的大小，如果大于等于平均值G1，把所述图像点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；如果所述图像点的灰度值小于平均值G1，把本层距离图像点最近且灰度值大于或等于平均值的点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；2)与所述第一层相邻且位于第一层下方的层命名为第二层，从所述第二层往下逐层进行生长，方法为：201)获取上一层种子点的灰度值的平均值Ai‑1，和上一层非种子点的灰度值的平均值Bi‑1，对平均值Ai‑1和平均值Bi‑1加权计算得到本层的种子点灰度最小值Mi；202)从上一层的种子点向本层进行生长，与种子点相邻且灰度值大于种子点灰度最小值Mi的点都生长为种子点；按步骤201和202所述方法逐层往下进行种子点的生长，直到没有新的种子点生成为止。所有生长出来的种子点的集合即组成了血管三维分割图像。...

【技术特征摘要】
1.基于CT影像的血管三维分割方法，其特征在于，包括：1)读取血管CT造影图像数据，把血管开始出现的层命名为第一层，并在该层生成种子点，所述种子点生成方法如下：101)获取所述第一层灰度值的平均值G1；102)从所述第一层中选取一个图像点，比较所述图像点的灰度值与平均值G1的大小，如果大于等于平均值G1，把所述图像点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；如果所述图像点的灰度值小于平均值G1，把本层距离图像点最近且灰度值大于或等于平均值的点记为第一个种子点，并在本层进行生长，与种子点相邻且灰度值大于平均值G1的点都生长为种子点；2)与所述第一层相邻且位于第一层下方的层命名为第二层，从所述第二层往下逐层进行生长，方法为：201)获取上一层种子点的灰度值的平均值Ai-1，和上一层非种子点的灰度值的平均值Bi-1，对平均值Ai-1和平均值Bi-1加权计算得到本层的种子点灰度最小值Mi；202)从上一层的种子点向本层进行生长，与种子点相邻且灰度值大于种子点灰度最小值Mi的点都生长为种子点；按步骤201和202所述方法逐层往下进行种子点的生长，直到没有新的种子点生成为止；所有生长出来的种子点的集合即组成了血管三维分割图像。2.根据权利要求1所述的基于CT影像的血管三维分割方法，其特征在于：所述步骤202还包括：获取本层种子点数量Ni，并与上一层的种子点数量Ni-1进行比较，如果Ni和Ni-1接近，即|Ni-1-Ni|/Ni-1≤10％，则本层种子点生长完成；否则，用最大类间方差法计算本层的阈值T，比较种子点灰度最小值Mi和阈值T，如果|Mi-T|≤S，则本层种子点生长完成，其中所述S为取值范围为[10,30]的正整数；否则，对Mi进行调整，并重新执行步骤202，所述对Mi进行调整的方法如下：如果Ni大于Ni-1，如果Ni小于Ni-1，按调整后的Mi重新执行步骤202。3.根据权利要求2所述的基于CT影像的血管三维分割方法，其特征在于，还包括：3)设含有生成种子点的最下面那...

【专利技术属性】
技术研发人员：唐铁兵，刘芳丽，程悦，
申请(专利权)人：北京欣方悦医疗科技有限公司，
类型：发明
国别省市：北京;11