一种基于智慧医疗的医学图像分割与显示方法与流程

1.本发明属于图像处理领域，具体涉及一种基于智慧医疗的医学图像分割与显示方法。


背景技术：

2.据最新的医疗数据显示，血管类疾病已严重影响当代人的生命健康，成为致死率较高的疾病之一。比如动脉粥样硬化、炎症性血管疾病、血管真性肿瘤性疾病等等。血管类疾病中常见的诱因是血管狭窄、堵塞、破裂，以及斑块等等。目前临床上对于血管病变程度与血管狭窄化程度的评估，通常使用基于管腔成像的方法，如数字减影血管造影术(digital subtraction angiography,dsa)、ct血管成像(computed tomography angiography,cta)、磁共振血管成像(magnetic resonance angiography,mra)，以及高分辨率磁共振血管成像(high-resolution magnetic resonance angiography,hrmra)等。
3.其中，mra和hrmra作为一种对患者无创的成像方法，可以清晰地检测到血管壁结构并进行分析，扫描得到的磁共振图像对于软组织的分辨率高，没有骨伪影，图像质量好，且能够使用多种序列扫描得到具有不同成像特点的组织结构，在血管的显示上具有明显的优越性。
4.但由于mra和hrmra技术得到的亮血序列、黑血序列对应的图像均为二维图像，在临床上，医生需要凭借经验结合两种图像的信息，来获得血管的综合情况，以进行血管病变分析。但二维图像具有局限性，不利于简便快速地获得血管的真实信息。并且，从上述成像方法所得到的图像中，无法直观、快速地得到关于血管狭窄程度的分析数据，不利于临床上对血管病灶区域的定位分析。


技术实现要素：

5.为了在临床应用上，简便、快速、直观地获得血管的真实信息和关于血管狭窄程度的分析数据，以进行血管病变分析。本发明实施例提供了一种基于智慧医疗的医学图像分割与显示方法。包括：
6.获取血管部位的亮血图像组、黑血图像组和增强黑血图像组；其中，所述亮血图像组、所述黑血图像组、所述增强黑血图像组分别包括k个亮血图像、黑血图像和增强黑血图像；所述亮血图像组、所述黑血图像组、所述增强黑血图像组中的图像一一对应；k为大于2的自然数；
7.针对所述亮血图像组中每一个亮血图像，以所述增强黑血图像组中对应的增强黑血图像为基准，利用基于互信息和图像金字塔的配准方法进行图像配准，得到配准后亮血图像组；
8.利用所述配准后亮血图像组，对所述增强黑血图像组中的增强黑血图像进行流空伪影消除操作，得到伪影消除增强黑血图像组；
9.利用所述配准后亮血图像组建立血液三维模型；
10.利用所述配准后亮血图像组建立血液边界扩展的血管三维模型；
11.基于所述伪影消除增强黑血图像组和所述黑血图像组，得到造影增强三维模型；
12.基于所述血液三维模型、所述血管三维模型、所述造影增强三维模型，得到血管造影增强三维模型；
13.获取所述血管造影增强三维模型中各段血管的表征血管狭窄程度的目标参数的数值，并利用各段血管的所述目标参数的数值对所述血管造影增强三维模型进行标记，得到血管造影增强三维狭窄化分析模型；
14.本发明的方案能够在临床应用上，简便、快速、直观地获得颅内血管的真实信息和关于颅内血管狭窄程度的分析数据，辅助医生更加精确、直观地进行病灶分析和判断。
15.当然，实施本发明的任一产品或方法并不一定需要同时达到以上所述的所有优点。
16.以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
17.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
18.图1为本发明实施例提供的一种基于智慧医疗的医学图像分割与显示方法的流程示意图；
19.图2为本发明实施例所提供的灰度线性变换及参数设置示意图
20.图3为本发明实施例的颅内血管磁共振图像的待配准区域示意图；
21.图4(a)为本发明实施例的颅内血管磁共振图像的亮血高斯金字塔和黑血高斯金字塔；图4(b)为本发明实施例的颅内血管磁共振图像的亮血拉普拉斯金字塔和黑血拉普拉斯金字塔；
22.图5为本发明实施例的颅内血管磁共振图像的拉普拉斯金字塔图像的配准结果；
23.图6为本发明实施例的颅内血管磁共振图像基于互信息的高斯金字塔图像配准步骤示意图；
24.图7为本发明实施例的不同迭代次数下的归一化互信息；
25.图8为多种配准方法的颅内血管磁共振图像的配准结果；
26.图9为本发明实施例的灰度线性变换结果图；
27.图10为本发明实施例的图像二值化结果图；
28.图11为本发明实施例针对颅内血管的流空伪影消除结果；
29.图12为本发明实施例针对颅内血管的血液三维模型效果图、血管三维模型效果图和造影增强三维模型效果图；
30.图13为本发明实施例针对颅内血管的血管造影增强三维模型效果图；
31.图14为本发明实施例的颅内血管的血管造影增强三维狭窄化分析模型效果图；
32.图15为本发明实施例的颅内血管的血管造影增强三维狭窄化分析模型和切面图显示效果图；
33.图16为本发明实施例提供的颅内血管的血管造影增强三维狭窄化分析模型裸眼3d全息可视化图；
34.图17为本发明实施例提供的对颅内血管的血管造影增强三维狭窄化分析模型的裸眼3d全息显示结果进行手势识别的示意图；
35.图18为本发明实施例提供的颅内血管的血管造影增强三维狭窄化分析模型的3d打印结果图。
具体实施方式
36.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行描述。
37.为了在临床应用上，简便、快速、直观地获得血管的真实信息和关于血管狭窄程度的分析数据，以进行血管病变分析。本发明实施例提供了一种基于智慧医疗的医学图像分割与显示方法。
38.如图1所示，本发明实施例所提供的一种基于智慧医疗的医学图像分割与显示方法，可以包括如下步骤：
39.s1，获取血管部位的亮血图像组、黑血图像组和增强黑血图像组；
40.其中，亮血图像组、黑血图像组、增强黑血图像组分别包括k个亮血图像、黑血图像和增强黑血图像；亮血图像组、黑血图像组、增强黑血图像组中的图像一一对应；k为大于2的自然数；
41.本发明实施例中，血管可以是颅内血管、心血管、眼底血管等组织部位的血管，在此，并不对本发明实施例的血管部位进行限制。
42.亮血图像组、黑血图像组、增强黑血图像组中的k个图像是一一对应的，其中对应方式是按照扫描时间形成的图像次序相同。
43.s2，针对亮血图像组中每一个亮血图像，以增强黑血图像组中对应的增强黑血图像为基准，利用基于互信息和图像金字塔的配准方法进行图像配准，得到配准后亮血图像组；
44.该步骤实际上是完成每一个亮血图像的图像配准，即将待配准的亮血图像作为浮动图像，将该亮血图像对应的增强黑血图像作为参考图像，利用基于互信息的相似性度量，并引入图像金字塔方法，进行图像配准。
45.可选的一种实施方式中，s2可以包括s21～s27：
46.s21，对每一个亮血图像和对应的增强黑血图像进行预处理，得到第一亮血图像和第一黑血图像；
47.可选的另一种实施方式中，s21可以包括s211和s212：
48.s211，针对每一个亮血图像，以对应的增强黑血图像为基准，将该亮血图像进行坐标变换和图像插值，并使用基于互信息的相似性度量及预定搜索策略，得到预配准后的第一亮血图像；
49.s211步骤实际上是对亮血图像以增强黑血图像为基准进行的图像预配准。
50.具体的，增强黑血图像和亮血图像为待配准的图像，可以根据亮血图像的dicom文件中的方位标签信息，以增强黑血图像的坐标系为基准，将增强黑血图像作为参考图像，将亮血图像作为浮动图像，将亮血图像进行坐标变换，实现旋转亮血图像至与增强黑血图像
相同坐标系的目的，旋转后亮血图像的扫描方向也变为呈冠状面。
51.为了便于理解本发明实施例的方法，以下结合图像配准过程进行简要介绍，具体过程可以借鉴相关现有技术进行了解。
52.对于两幅图像a和b的配准，实际上就是将图像a中的每个坐标位置通过一个映射关系，对应到图像b中。具体的坐标变换方法可以包括刚体变换、仿射变换、投影变换和非线性变换等。但是在坐标变换过程中，浮动图像的坐标系会发生拉伸或形变，经过坐标变换后的图像像素坐标与原图的采样网格并不会完全重合，即原先为整数的像素坐标点经过坐标变换后可能不再是整数，导致图像的有些区域丢失掉部分像素，因此在图像坐标变换的过程中，需要同时对图像进行重采样插值，来确定经过坐标变换后图像像素坐标点的灰度值，便于后续处理。图像插值方法可以采用最近邻插值、双线性插值和双三次插值等方法中的任意一种。
53.使用图像插值方法将缺失的像素点进行图像恢复之后，还需要使用某种相似性度量来计算参考图像和变化后的浮动图像的相似性，接着利用搜索策略找到最优的相似性度量，循环反复迭代求优，直到两幅图像的相似性度量达到最优时，迭代停止，最后根据确定的空间变换矩阵(旋转矩阵)对浮动图像进行坐标转换，实现图像完全配准。待配准的图像在经历一个迭代算法优化后，能够计算得到两幅图像的空间位置配准关系以及配准图像，使得配准后的浮动图像与参考图像相似度达到最高。
54.在本发明实施例中，采用基于互信息的相似性度量，具体可以为互信息或归一化互信息。
55.图像配准实质上是一个多参数的优化问题，即通过使用某种搜索策略对图像进行空间坐标变化，最终使得两幅图像的相似性度量达到最优，其中搜索策略与空间坐标变化在实际计算过程中是彼此交叉进行的。算法思想是在每次迭代中计算两幅图像之间的相似性度量，并通过平移或旋转等坐标变换的操作调整浮动图像，同时对图像进行插值，一直到两幅图像的相似性度量最大为止。目前常用的搜索策略包括梯度下降优化器、基于进化策略(evolution strategy,es)的(1 1)-es等等，本发明实施例中的预定搜索策略可以根据需要选择。
56.通过该步骤的预配准，可以初步实现在同一坐标系下对比相同扫描层面的磁共振图像，但是由于亮血序列和黑血序列扫描的时间不同，且患者可能在扫描前后发生了轻微的运动，所以上述操作只是一个粗糙的坐标变换，仅通过预配准并不能实现多模态磁共振图像的完全配准，但是该步骤可以为后续的精确配准环节省略不必要的处理过程，提高处理速度。
57.s212，从对应的增强黑血图像中，提取与第一亮血图像的扫描范围相同的区域内容，形成第一黑血图像。
58.由于血管成像在不同磁共振序列中的扫描范围不同，当亮血图像经过图像坐标变换后，其冠状面的信息并没有增强黑血图像的信息丰富，因此为能够更快速、准确地配准好两种序列下的同一区域，可以根据第一亮血图像的扫描区域，在增强黑血图像中提取出相同的扫描区域，缩小后续图像的配准范围。
59.可选的，s212可以包括以下步骤：
60.1.获得第一亮血图像中血管的边缘轮廓信息；
61.具体可以使用sobel边缘检测方法等方法获取边缘轮廓信息。边缘轮廓信息包含各个边缘点的坐标值。
62.2.在边缘轮廓信息中，提取横坐标、纵坐标的最小值和最大值，基于获得的四个坐标值确定初始提取边框；
63.也就是在边缘轮廓信息中，提取最小横坐标值、最大横坐标值、最小纵坐标值、最大纵坐标值，利用这四个坐标值确定方形边框的四个顶点，从而得到初始提取边框；
64.3.在第一亮血图像的尺寸范围内，将初始提取边框沿四个方向分别扩大预设数量个像素大小，得到最终提取边框；
65.其中，四个方向分别是横纵坐标的正负方向；预设数量根据血管图像的类型合理选择，目的是保证扩大后的最终提取边框不超过第一亮血图像的尺寸范围，比如预设数量可以为20等。
66.4.在对应的增强黑血图像中提取最终提取边框中对应的区域内容，形成第一黑血图像。
67.依据最终提取边框划定的坐标范围，提取对应的增强黑血图像中对应区域的内容，将提取出的内容形成第一黑血图像。该步骤通过提取待配准区域来获取两种模态下磁共振图像的共同扫描范围，有利于后续的快速配准。
68.本发明实施例的上述两步预处理过程有着非常重要的作用，经过该预处理后的图像可以更多地关注有用信息，排除无关信息，在实际使用中，使用该图像预处理可以提高图像配准、识别的可靠性。
69.本发明实施例中，为提高图像配准的准确度，避免图像在配准过程中收敛到局部最大值，选择使用多分辨率策略来解决局部极值的问题，同时利用多分辨率策略在满足图像配准精度的条件下，提高算法执行速度，增加鲁棒性。因此采用图像金字塔方法。构建图像金字塔就是一种通过增加模型复杂度来提高配准精度和速度的有效方式，即在配准过程中，按从粗配准到精配准的顺序执行，首先对低分辨率的图像进行配准，接着在低分辨率图像配准完成的基础上，对分辨率高的图像进行配准。可选的，可以采用以下步骤：
70.s22，基于下采样处理，由第一亮血图像得到亮血高斯金字塔，由第一黑血图像得到黑血高斯金字塔；其中，亮血高斯金字塔和黑血高斯金字塔中包括从下至上分辨率依次变小的m个图像；m为大于3的自然数；
71.可选的一种实施方式中，s22可以包括以下步骤：
72.获取第i层的输入图像，对第i层的输入图像利用高斯内核进行滤波，并删除滤波后图像的偶数行和偶数列，得到高斯金字塔的第i层图像gi，并将第i层图像gi作为第i 1层的输入图像，得到高斯金字塔的第i 1层图像g
i 1
；
73.其中，i＝1、2，
…
，m-1；当高斯金字塔为亮血高斯金字塔时，第1层的输入图像为第一亮血图像，当高斯金字塔为黑血高斯金字塔时，第1层的输入图像为第一黑血图像。
74.具体的，高斯金字塔中的多个图像是同一原图像不同分辨率所对应的图像。高斯金字塔通过高斯滤波和下采样来获取图像，它的每一层构建步骤可以分为两步：首先使用高斯滤波对图像进行平滑滤波，即采用高斯内核进行滤波；接着删除滤波后图像的偶数行和偶数列，即将低一层图像的宽和高缩小一半，得到当前层图像，因此当前层图像为低一层图像大小的四分之一，通过不断地迭代以上步骤，最终可以得到高斯金字塔。
75.虽然使用二维的高斯模板可以实现模糊图像的效果，但是当一个点在边界，周围没有足够的点时，会由于权重矩阵的关系导致边缘图像缺失，因此本发明实施例对二维高斯模板进行优化。可以将二维高斯滤波拆分成两个独立的一维高斯滤波，分别在横纵两个方向上进行图像滤波。将高斯函数分离，不仅能够消除二维高斯模板产生的边缘，还可以极大地加快程序的运行速度。相较于其他模糊滤波器，高斯滤波既可以实现图像的模糊效果，又能更好地保留边际效果。
76.本步骤中通过对经过预处理后的第一亮血图像和第一黑血图像进行上述处理，可以得到亮血高斯金字塔和黑血高斯金字塔。其中，图像层数m可以为4。
77.由于高斯金字塔是向下采样，即缩小图像，因此会丢失掉图像的一部分数据。因此，本发明实施例为避免图像在缩放过程中的数据缺失，恢复细节数据，在后续步骤使用拉普拉斯金字塔，配合高斯金字塔一起实现图像重建，在高斯金字塔图像的基础上突出细节。
78.s23，基于上采样处理，利用亮血高斯金字塔得到亮血拉普拉斯金字塔，利用黑血高斯金字塔得到黑血拉普拉斯金字塔；其中，亮血拉普拉斯金字塔和黑血拉普拉斯金字塔中包括从下至上分辨率依次变小的m-1个图像；
79.可选的一种实施方式中，s23可以包括以下步骤：
80.对高斯金字塔的第i 1层图像g
i 1
进行上采样，并用数据0填充新增的行和列，得到填充图像；
81.对填充图像利用高斯内核进行卷积，获得填充像素的近似值，得到放大图像；
82.将高斯金字塔的第i层图像gi与放大图像相减，得到拉普拉斯金字塔的第i层图像li；
83.其中，当高斯金字塔为亮血高斯金字塔时，拉普拉斯金字塔为亮血拉普拉斯金字塔，当高斯金字塔为黑血高斯金字塔时，拉普拉斯金字塔为黑血拉普拉斯金字塔。
84.由于拉普拉斯金字塔是图像经过下采样操作后与原图的残差，因此从下至上对比，拉普拉斯金字塔比高斯金字塔结构少一层高层图像。
85.具体的，生成拉普拉斯金字塔结构的数学公式如(1)所示，其中li表示第i层拉普拉斯金字塔(亮血拉普拉斯金字塔或黑血拉普拉斯金字塔)，gi表示第i层高斯金字塔(亮血高斯金字塔或黑血高斯金字塔)，而up操作为向上采样放大图像，符号是卷积符号，是在构建高斯金字塔中使用的高斯内核。此公式表明了拉普拉斯金字塔实质上是使用原图像减去先缩小、再放大的图像的残差数据构成的，是一种残差预测金字塔，核心思想是用来存储图像经过下采样操作后与原图的差异，保留图像的高频信息，目的是为能够完整地恢复出每一层级进行下采样操作前的图像。由于之前下采样操作丢失的一部分信息并不能通过上采样来完全恢复，即下采样是不可逆的，所以图像先经过下采样，再进行上采样后的显示效果比原图模糊。通过存储图像经过下采样操作后与原图的残差，能够在高斯金字塔图像的基础上给不同频率层的图像增加细节，对细节等进行突出。
[0086][0087]
对应于4层的高斯金字塔，本步骤可以得到图像层数为3的亮血拉普拉斯金字塔和黑血拉普拉斯金字塔。
[0088]
s24，对亮血拉普拉斯金字塔和黑血拉普拉斯金字塔中对应层的图像进行配准，得
到配准的亮血拉普拉斯金字塔；
[0089]
可选的一种实施方式中，s24可以包括以下步骤：
[0090]
针对亮血拉普拉斯金字塔和黑血拉普拉斯金字塔中的每一层，将该层对应的黑血拉普拉斯图像作为参考图像，将该层对应的亮血拉普拉斯图像作为浮动图像，使用基于互信息的相似性度量及预定搜索策略，实现图像配准，得到配准后的该层亮血拉普拉斯图像；
[0091]
由配准后的多层亮血拉普拉斯图像，依据分辨率依次减小的顺序，从下至上构成配准的亮血拉普拉斯金字塔；
[0092]
其中，黑血拉普拉斯图像为黑血拉普拉斯金字塔中的图像，亮血拉普拉斯图像为亮血拉普拉斯金字塔中的图像。
[0093]
该步骤中的配准过程类似于前述的预配准过程，通过对亮血拉普拉斯图像进行坐标变换、图像插值，使用基于互信息的相似性度量及预定搜索策略，实现图像配准，可以得到配准后的亮血拉普拉斯图像。
[0094]
s25，利用配准的亮血拉普拉斯金字塔作为叠加信息，对亮血高斯金字塔和黑血高斯金字塔中各层图像进行自上而下的配准，得到配准的亮血高斯金字塔；
[0095]
针对s25，利用配准的亮血拉普拉斯金字塔作为叠加信息，对亮血高斯金字塔和黑血高斯金字塔中各层图像进行自上而下的配准，需要对高斯金字塔中不同分辨率的图像进行配准，由于低分辨率图像的配准可以更容易把握住图像的本质特征，因此本发明实施例在低分辨率图像配准的基础上配准高分辨率图像，即将高斯金字塔图像从上到下配准，将上一层图像的配准结果作为下一层图像配准的输入。
[0096]
可选的一种实施方式中，s25可以包括以下步骤：
[0097]
对亮血高斯金字塔和黑血高斯金字塔中自上而下的第j层，将该层对应的黑血高斯图像作为参考图像，将该层对应的亮血高斯图像作为浮动图像，使用基于互信息的相似性度量及预定搜索策略，实现图像配准，得到配准后的第j层亮血高斯图像；
[0098]
将配准后的第j层亮血高斯图像进行上采样操作，并与配准后的对应层亮血拉普拉斯图像相加，利用相加后的图像替换亮血高斯金字塔中第j 1层的亮血高斯图像；
[0099]
将第j 1层的黑血高斯图像作为参考图像，将替换后的第j 1层的亮血高斯图像作为浮动图像，使用预定相似性度量及预定搜索策略，实现图像配准，得到配准后的第j 1层亮血高斯图像；其中j＝1，2，
…
，m-1，黑血高斯图像为黑血高斯金字塔中的图像，亮血高斯图像为亮血高斯金字塔中的图像。
[0100]
重复以上操作，直至完成底层高斯金字塔图像的高分辨率配准，得到配准的亮血高斯金字塔。其中的亮血图像的坐标系与黑血图像的坐标系一致，且图像具备较高的相似性。配准过程类似于前述的预配准过程，不再赘述。
[0101]
s26，基于配准的亮血高斯金字塔得到该亮血图像对应的配准后亮血图像；
[0102]
该步骤中，获取配准的亮血高斯金字塔中的底层图像作为配准后亮血图像。
[0103]
s27，由k个亮血图像分别对应的配准后亮血图像，得到配准后亮血图像组。
[0104]
所有亮血图像均配准完成后，由k个配准后亮血图像可以得到配准后亮血图像组。每一个配准后亮血图像和对应的增强黑血图像可以作为配准后图像对。
[0105]
通过以上步骤，可以实现亮血图像和增强黑血图像的图像配准，在本发明实施例提供的配准方案中，基于互信息作为相似性度量，可以提高配准精度；同时引入了图像金字
塔算法，它是一种通过增加模型复杂度来提高配准精度和速度的有效方式，即首先对分辨率较低的图像进行粗配准，接着在粗配准的基础上，对分辨率高的图像进行精配准，通过运用高斯金字塔和拉普拉斯金字塔来对血管图像进行分解和重建，模拟了人眼在不同距离下观测一幅图像的效果，更容易获取血管图像的本质特征。使用金字塔算法对血管部位的磁共振亮血图像与黑血图像进行配准，可以提高配准效率，使图像从低分辨率到高分辨率逐层提高配准精度。通过上述图像配准可以将亮血图像和增强黑血图像统一在相同坐标系下，可以方便医生理解黑血序列和亮血序列对应的血管图像，简便、快速地得到诊断所需的综合信息，为后续的医疗诊断、制定手术计划、放射治疗计划等提供准确可靠的参考信息。本发明实施例提供的配准方案，可以给其余医学图像配准提供一种较好的参考方式，具有很大的临床应用价值。同时，本发明实施例的图像配准过程是后续消除流空伪影的重要基础。
[0106]
在图像配准后，可以对配准后增强黑血图像中的流空伪影进行消除，其中流空伪影出现的原因是在血管壁成像过程中，由于血管太细小，走向迂曲处血液流速较慢，以及周围血液和组织液可能有信号污染等问题，导致在黑血序列扫描得到的图像中，本应该为黑色的血液信息反而表现为亮色，从而模拟正常个体的壁增厚或斑块外观，夸大血管狭窄程度。本发明实施例考虑利用配准后亮血图像中的血液信息，对配准后增强黑血图像中信号显示不正确的血液信息进行修正，将配准后亮血图像中的血液信息嵌入配准后增强黑血图像中，以达到图像融合的效果。具体可以通过以下步骤实现：
[0107]
s3，利用配准后亮血图像组，对增强黑血图像组中的增强黑血图像进行流空伪影消除操作，得到伪影消除增强黑血图像组；
[0108]
可选的一种实施方式中，s3可以包括s31～s34：
[0109]
s31，针对每一个配准后亮血图像，提高该配准后亮血图像的对比度，得到对比增强亮血图像；
[0110]
步骤可选的一种实施方式中，可以根据在亮血图像中血液呈高信号，而周围组织呈低信号的特点，对配准后亮血图像进行灰度线性变换，调整图像灰度范围，实现提高图像对比度的目的。
[0111]
比如，一种配准后亮血图像所使用的灰度线性变换及参数设置如图2所示，图2为本发明实施例所提供的灰度线性变换及参数设置示意图。利用图2所示的灰度线性变换，可以将原来的配准后亮血图像f中较小的灰度值变化区间扩展为新的配准后亮血图像f1(增强亮血图像)中的较大灰度值变化区间，调整图像灰度范围，实现提高配准后亮血图像的对比度的目的。通过该步骤，可以得到对比增强亮血图像，实现针对k个配准后亮血图像得到对应的k个对比增强亮血图像。由于医学图像像素范围较大，可能是-1000～ 1000，通过该步骤，可以将像素范围归一化为0～255，成为符合一般图像处理的像素范围，可便于后续处理。
[0112]
关于灰度线性变换的具体过程可以参见相关现有技术，在此不再赘述。
[0113]
s32，从对比增强亮血图像中提取出血液信息，得到亮血特征图；
[0114]
可选的一种实施方式中，s32可以包括以下步骤：
[0115]
s321，利用预设图像二值化方法确定第一阈值；
[0116]
s322，利用第一阈值，从对比增强亮血图像中提取出血液信息；
[0117]
该步骤使用的方法称为阈值分割。
[0118]
s323，由提取出的血液信息得到亮血特征图。
[0119]
本发明实施例通过预设图像二值化方法可以将对比增强亮血图像中的血液信息突出显示为白色，将无关信息显示为黑色，以便于提取出血液信息对应的亮血特征图。本发明实施例中的预设图像二值化方法可以包括最大类间方差法otsu、kittle等等。
[0120]
血液信息的提取公式如(2)所示，其中t(x,y)为对比增强亮血图像灰度值，f(x,y)为亮血特征图灰度值，t为第一阈值。
[0121][0122]
s33，将亮血特征图与该配准后亮血图像对应的增强黑血图像，依据预设融合公式进行图像融合，得到该增强黑血图像对应的流空伪影消除的目标增强黑血图像；
[0123]
在该步骤中，首先建立亮血特征图与对应的增强黑血图像之间的空间映射关系，将亮血特征图映射到对应的增强黑血图像中，依据预设融合公式进行图像融合，其中预设融合公式为：
[0124][0125]
其中，f(x,y)为亮血特征图的灰度值，r(x,y)为对应的增强黑血图像的灰度值，g(x,y)为融合后的目标增强黑血图像的灰度值。
[0126]
经过以上操作，可以将对应的增强黑血图像中本应该为黑色，却表现为亮色的流空伪影灰度值更改为黑色，从而实现消除流空伪影的目的。
[0127]
s34，由k个增强黑血图像对应的目标增强黑血图像，得到伪影消除增强黑血图像组。
[0128]
所有增强黑血图像均完成流空伪影消除后，可以得到伪影消除增强黑血图像组。
[0129]
s4，利用配准后亮血图像组建立血液三维模型；
[0130]
目前的三维重建技术包括移动立方体(marching cubes,mc)方法、最大密度投影(maximum intensity projection,mip)方法、表面阴影遮盖方法(ssd)、容积漫游技术(vrt)、曲面重建方法(cpr)、虚拟内镜技术(ve)等等。本发明实施例可以采用任意一种三维重建方法建立血液三维模型。该血液三维模型可以初步模拟三维血管，直观地显示血管的走向及病灶区域等。
[0131]
可选的一种实施方式中，s4可以包括s41～s43：
[0132]
s41，获取k个对比增强亮血图像构成的第一三维体数据；
[0133]
可以获取s31步骤得到的k个对比增强亮血图像。本领域技术人员可以理解的是，k个对比增强亮血图像其实是层叠为一个三维立方体数据的。为了方便区分，本发明实施例中将其命名为第一三维体数据。
[0134]
s42，利用最大类间方差法计算第一三维体数据中，居中的第二三维体数据对应的第二阈值；
[0135]
该步骤仍是利用最大类间方差法otsu确定阈值，但和s321中利用该方法确定第一阈值有所不同，本步骤中是使用最大类间方差法otsu求出位于第一三维体数据这个大的三
维立方体中，靠近中间部分的一个小立方体(称为第二三维体数据)中的多个对比增强亮血图像所对应的一个阈值作为第二阈值。因为在对比增强图像中，血液信息基本集中于图像的中部，那么，针对第一三维体数据中，选取居中的小的立方体数据(第二三维体数据)确定第二阈值，可以减小阈值计算量，提高计算速度，且该第二阈值准确适用于第一三维体数据中所有血液信息。
[0136]
对于第二三维体数据的大小，可以首先确定第一三维体数据的中心点，然后以预设边长在立方体对应的六个方向延伸，从而确定第二三维体数据的大小，其中，预设边长可以根据经验值确定，比如为第一三维体数据这个立方体的边长的1/4等。
[0137]
s43，将第二阈值作为移动立方体方法的输入阈值，利用移动立方体方法对第一三维体数据进行处理，得到血液三维模型。
[0138]
如前，移动立方体方法(简称mc)是一种三维重建方法，可以依据给定的输入阈值，对第一三维体数据进行处理，直接得到血液三维模型。
[0139]
移动立方体方法相较于其他的面绘制算法，具有生成网格质量好的优点。关于移动立方体方法对第一三维体数据的具体处理过程，请参见相关的现有技术，在此不做赘述。
[0140]
s5，利用配准后亮血图像组建立血液边界扩展的血管三维模型；
[0141]
s4步骤中得到的是血液三维模型，其表征的其实是血液的流向和区域分布，由于实际中血液外围存在有血管壁，因此血液三维模型其实并不能完全代表真实的血管情况。
[0142]
因此，在s5步骤中，可以对配准后亮血图像中的血液边界进行扩展，使得扩展后的血液边界能够涵盖血管壁的范围，形成一个中空管的效果，再对扩展血液边界后的二维图像利用三维重建方法生成三维模型，进而得到比s4中的血液三维模型更接近真实血管情况的血管三维模型。
[0143]
关于血液边界的扩展可以通过检测配准后亮血图像中血液边界像素点，将检测到的像素点向预设方向扩展预设个像素点实现，预设个像素点可以根据大量血管直径及血管壁厚度数据所得到的经验值来选取。当然，本发明实施例中扩展血液边界的方式不限于此。
[0144]
可选的一种实施方式中，s5可以包括s51～s55：
[0145]
s51，获取k个亮血特征图；
[0146]
即获取s32步骤得到的k个亮血特征图。
[0147]
s52，针对每个亮血特征图，利用膨胀操作扩大该亮血特征图中血液的边界，得到该亮血特征图对应的扩展亮血特征图；
[0148]
一种可选的实施方式中，可以利用半径为1的圆形内核对亮血特征图进行多步膨胀，直到达到最大梯度位置停止膨胀，从而确定血管外壁边界，实现血管壁的分割，得到亮血特征图对应的扩展亮血特征图。由于血管壁紧贴血液，且管壁极薄，假设膨胀后的范围就是血管壁的所在范围，这步操作即可将血液附近的血管壁所在区域包括进来，作为血管壁造影增强特性的搜索范围。
[0149]
关于膨胀操作的具体实施过程可以参见相关现有技术，在此不再赘述。
[0150]
s53，将该亮血特征图对应的扩展亮血特征图与该亮血特征图求差，得到该亮血特征图对应的差值特征图；
[0151]
该步骤针对每个亮血特征图得到的差值特征图是一个类似于中空血管的二维平面图。同样的，该差值特征图的像素值也仅有0和255。
[0152]
s54，确定第三阈值；
[0153]
该步骤可以根据经验值为所有差值特征图选定一个像素值作为第三阈值，比如可以选取100～200之间的任意一个值，如128作为第三阈值。
[0154]
s55，将第三阈值作为移动立方体方法的输入阈值，利用移动立方体方法对k个差值特征图进行处理，得到血液边界扩展的血管三维模型。
[0155]
移动立方体方法利用第三阈值作为输入阈值，可以由k个差值特征图得到血液边界扩展的血管三维模型。关于移动立方体方法的具体实施过程在此不再赘述。
[0156]
s6，基于伪影消除增强黑血图像组和黑血图像组，得到造影增强三维模型；
[0157]
可选的一种实施方式中，s6包括s61和s62：
[0158]
s61，将伪影消除增强黑血图像组和黑血图像组中对应图像相减，得到k个造影增强图；
[0159]
每个目标增强黑血图像和对应的黑血图像相减，可以得到具有造影增强效果的造影增强图，当所有目标增强黑血图像和对应的黑血图像均相减后，可以得到k个造影增强图，可以理解的是，这k个造影增强图均是二维图。
[0160]
s62，利用k个造影增强图建立造影增强三维模型。
[0161]
该步骤可以利用移动立方体方法实现，具体参见s4和s5，在此不再赘述。
[0162]
s7，基于血液三维模型、血管三维模型、造影增强三维模型，得到血管造影增强三维模型；
[0163]
可选的一种实施方式中，s7可以包括s71和s72：
[0164]
s71，保留造影增强三维模型中与血管三维模型的重叠部分，得到保留后造影增强三维模型；
[0165]
由于s6得到的造影增强三维模型并非只包含了血管的造影增强，需要排除无关组织的增强特性，因此使用s5得到的血管三维模型中血管壁造影增强特性的搜索范围，来判断s6得到的造影增强三维模型是否位于血液附近的血管壁区域，即判断造影增强三维模型中是否有与血管三维模型的重叠部分，如果是，则表明重叠部分位于搜索范围之内，则需要保留该重叠部分，因此得到保留后造影增强三维模型。
[0166]
s72，将保留后造影增强三维模型与血液三维模型融合，得到血管造影增强三维模型。
[0167]
将表征血管造影增强的保留后造影增强三维模型，与表征血液信息的血液三维模型进行融合，可以直观显示有明显造影增强的血管壁，可以清晰地看到血管的哪个部位范围内造影增强效果最为明显，那么该区域可能出现粥样硬化或易损性斑块。
[0168]
可选的一种实施方式中，血管造影增强三维模型中可以获得造影增强定量分析，具体的，可以针对血管造影增强三维模型中血管壁上任意一个点，得到斑块强化指数ce，其中ce定义为：
[0169][0170]
其中，s
prebbmr
和s
postbbmr
分别为黑血图像和造影增强黑血图像中的信号强度。
[0171]
本领域技术人员可以理解的是，s
prebbmr
和s
postbbmr
分别是拍摄黑血图像和造影增强黑血图像后，图像中携带的信息。本发明实施例利用上述信息得到血管壁边沿各个点的斑
块强化指数ce，并将其体现在血管造影增强三维模型中，可以方便医生获取更为详细的血管信息，具体的，当ce大于一个斑块阈值，比如0.5时，表示该处血管壁上出现了斑块，因此，通过测量血管壁区域的斑块强化指数，有助于鉴别责任动脉斑块等，可以提供有价值的诊断辅助信息。
[0172]
关于两个三维模型的融合技术可以采用现有技术实现，在此不做赘述。
[0173]
s8，获取血管造影增强三维模型中各段血管的表征血管狭窄程度的目标参数的数值，并利用各段血管的目标参数的数值对血管造影增强三维模型进行标记，得到血管造影增强三维狭窄化分析模型。
[0174]
可选的一种实施方式中，s8可以包括s81～s84：
[0175]
s81，针对血管造影增强三维模型中的每一段血管，从预设的三个方位进行切分，获得各方位的二维切面图；
[0176]
该步骤中，可以先将血管造影增强三维模型中的血管进行划分，针对每一段血管，从预设的三个方位进行切分，获得各方位的二维切面图。
[0177]
其中，预设的三个方位包括：轴位、冠状位和矢状位。
[0178]
关于对血管造影增强三维模型进行某一方位切分，得到该方位的二维切面图，可以采用现有技术实现，在此不再赘述。
[0179]
s82，将每个方位的二维切面图中的血管进行腐蚀操作，记录血管腐蚀至单个像素时的目标腐蚀次数；
[0180]
腐蚀操作是形态学运算的一种，腐蚀操作能够消去物体的边缘数据，腐蚀后的物体比原先面积小，甚至会完全消失，腐蚀还可以断开一些细小狭长的连通区域。
[0181]
当血管较粗时，可以进行多次腐蚀操作，而血管较细小时，却只能进行少数几次腐蚀操作，本发明实施例是根据腐蚀操作达到血管对应部分达到单个像素时的次数估计血管的粗细程度。
[0182]
在s82步骤中，针对轴位的二维切面图中的血管进行腐蚀操作，记录该方位的二维切面图中的血管腐蚀至单个像素时对应的目标腐蚀次数n1；针对冠状位的二维切面图中的血管进行腐蚀操作，记录该方位的二维切面图中的血管腐蚀至单个像素时对应的目标腐蚀次数n2；针对矢状位的二维切面图中的血管进行腐蚀操作，记录该方位的二维切面图中的血管腐蚀至单个像素时对应的目标腐蚀次数n3。
[0183]
s83，根据该段血管在三个方位分别对应的目标腐蚀次数，得到表征该段血管狭窄程度的目标参数的数值；
[0184]
可选的一种实施方式中，目标参数包括狭窄率和/或扁平度；本领域技术人员可以理解，这两个参数均可以表征血管狭窄程度。
[0185]
当目标参数包括狭窄率时，s83可以包括：
[0186]
根据n1、n2、n3，利用血管的狭窄率公式得到该段血管的狭窄率的数值；其中，狭窄率公式为：
[0187][0188]
其中，分辨率为各方位二维切面图的分辨率(三个方位的二维切面图的分辨率相同)，狭窄率的数值越小表明血管越狭窄。
[0189]
当目标参数包括扁平度时，s83可以包括：
[0190]
根据n1、n2、n3，利用血管的扁平度公式得到该段血管的扁平度的数值；其中，扁平度公式为：
[0191][0192]
扁平度的数值越大表明血管越狭窄。
[0193]
s84，利用各段血管的目标参数的数值对血管造影增强三维模型进行标记，得到血管造影增强三维狭窄化分析模型。
[0194]
通过上述步骤，可以得到每一段血管的目标参数的数值，那么可以将各段血管的这些数值在血管造影增强三维模型上进行标记，得到血管造影增强三维狭窄化分析模型。即将各个点的目标参数的数值嵌入在该血管造影增强三维狭窄化分析模型中，以便于在需要时，可以提取出并显示每个点的目标参数的数值，便于医生在观察整体三维血管状态时，及时获取各个位置的血管狭窄程度的数据。
[0195]
可选的一种实施方式中，s84可以包括：
[0196]
利用各段血管的目标参数的数值，采用与各数值对应的颜色对血管造影增强三维模型进行标记，得到血管造影增强三维狭窄化分析模型。
[0197]
为了便于直观显示，可以将不同的数值利用不同的颜色标记在血管造影增强三维模型上得到血管造影增强三维狭窄化分析模型，比如针对狭窄率数值由小至大可以用由浅至深的多个颜色对应标记，针对扁平度数值，由于数值较少，可能只有2种数值，那么可以采用与狭窄率区分的两种颜色对应标记。采用不同色调的颜色显示可以更加直观地表现出血管的狭窄化程度，便于引起医生的注意。
[0198]
其中，优选的实施方式中，可以在一个血管造影增强三维狭窄化分析模型上用不同数值对应的颜色标记狭窄率数值，在另一个血管造影增强三维狭窄化分析模型上用不同数值对应的颜色标记扁平度数值，以便于医生分别观察狭窄率情况和扁平度情况。
[0199]
并且进一步的，由于医生习惯于观察切平面的二维医学图像，本发明实施例在血管狭窄化分析的同时，还可以提供各方位对应的二维切面图，即将该血管造影增强三维狭窄化分析模型中的每一个点对应的当前点所在冠状面、矢状面、轴状面图像一并显示出来。在显示该血管造影增强三维狭窄化分析模型时，还可以使用诸如三种颜色的点，实现两个点测距，三个点测量角度的功能，并在显示屏左下方显示出来，在显示屏右下方显示当前选取模型的体积大小。以便于医生获取颅内血管更为详细的数据。
[0200]
s9，将血管造影增强三维狭窄化分析模型进行显示。
[0201]
上述步骤得到对血管狭窄进行标记的血管造影增强三维狭窄化分析模型，可以直接通过软件在计算机显示屏上显示，当然，也可以采用其它更为直观的方法进行显示。
[0202]
作为本发明的一种实施方式，将血管造影增强三维狭窄化分析模型进行显示，具体可以采用裸眼3d全息显示系统进行显示。本发明的方案并不需要采取任何穿戴式设备，如vr或mr眼镜，而是采用裸眼3d全息显示系统，将前后左右四个角度的图像通过软件分别投射到金字塔全息玻璃上，便于多位医生围在一起，同时看清楚颅内血管三维结构以及病变位置；且具有成像空间大、分辨率高、安静、方便讨论，及成本较低等优点。
[0203]
为了进一步增强得到血管模型的立体感觉，以及增加医生的观看代入感，在裸眼3d全息显示的基础上，也可以进一步采用手势识别对所述裸眼3d全息显示的血管造影增强三维狭窄化分析模型进行操作，比如，手势识别可以采用leap motion体感控制器，以进行手动缩放、旋转、切割、虚拟手术等操作。本发明的方案采用的手势识别技术，可以利用红外led 灰阶照相机的方式采集双手数据；前者利用双目视觉原理测量深度，后者提取关键点，从而重建出手掌在真实三维世界中的信息。
[0204]
当然，采用手势识别也可以针对计算机显示屏上输出的血管造影增强三维狭窄化分析模型直接进行手动缩放、旋转、切割、虚拟手术等操作。采用手势识别的优点是体积小，识别精度高，不受环境光源的限制，还能进行测距。作为本发明的另一种实施方式，将血管造影增强三维狭窄化分析模型进行显示，具体可以将所述血管造影增强三维狭窄化分析模型导出为stl文件并通过3d打印显示。通过将最终得到的血管模型进行3d打印显示，并对比正常的血管三维模型，可以直观的看出哪个位置血管出现狭窄，发生了病变。
[0205]
需要说明的是，上述用于显示的裸眼3d全息显示、手势识别和3d打印均可采用现有技术中相应的技术，在此不再赘述。
[0206]
本发明实施例所提供的方案中，首先对磁共振技术扫描得到的亮血图像和增强黑血图像采用基于互信息和图像金字塔的配准方法进行图像配准，可以提高配准效率，使图像从低分辨率到高分辨率逐层提高配准精度。通过上述图像配准可以将亮血图像和增强黑血图像统一在相同坐标系下。其次利用配准后亮血图像对增强黑血图像进行流空伪影消除操作，可以显示更准确、全面的血管信息。本发明实施例所提供的方案是从图像后处理的角度对流空伪影进行消除，无需使用新的成像技术、成像模式或脉冲序列，因此可以简便、准确、快速地消除流空伪影，并可以在临床应用实现较好的推广。再次，利用配准后亮血图像建立血液三维模型，和血液边界扩展的血管三维模型，并基于伪影消除增强黑血图像组和黑血图像组得到具有造影增强效果的造影增强三维模型；之后基于血液三维模型、血管三维模型、造影增强三维模型，得到具有造影增强效果的血管壁对应的血管造影增强三维模型。最后利用血管造影增强三维模型中表征血管狭窄程度的目标参数的数值进行标记，得到血管造影增强三维狭窄化分析模型。该血管造影增强三维狭窄化分析模型实现了血管的三维可视化，无需医生通过想象力还原血管组织结构及病症特征等，可以提供形象的血管三维空间信息，便于直观观测，便于定位与显示狭窄的病灶区域。能够在临床应用上，简便、快速、直观地获得血管的真实信息和关于血管狭窄程度的分析数据。
[0207]
以下以颅内血管为例，说明本发明实施例所提供的基于智慧医疗的医学图像分割与显示方法的实施过程和实施效果。实施过程可以包括以下步骤：
[0208]
步骤一，获取血管部位的亮血图像组、黑血图像组和增强黑血图像组；
[0209]
步骤二，针对亮血图像组中每一个亮血图像，以增强黑血图像组中对应的增强黑血图像为基准，利用基于互信息和图像金字塔的配准方法进行图像配准，得到配准后亮血图像组；
[0210]
在该步骤可以参见前文s2，包括：
[0211]
(一)，对每一个亮血图像和对应的增强黑血图像进行预处理，得到第一亮血图像和第一黑血图像；该预处理过程可以分为两个主要步骤：
[0212]
(1)预配准：
[0213]
由于颅内血管可以看作是一个刚体，因此，该步骤选用刚体变换作为坐标变换方法。具体预配准过程参见s211步骤。
[0214]
本发明实施例对亮血图像的图像插值方法进行仿真实验，先将原图缩小50％，然后使用不同插值算法得到与原图大小相同的效果图，并与原图像做对比。并利用均方根误差rmse、峰值信噪比psnr、归一化互相关系数ncc、归一化互信息nmi以及耗时time进行评判。从实验数据得知，双三次插值的精度明显优于最近邻插值和双线性插值。因此，本发明实施例采用双三次插值。
[0215]
实验中，使用梯度下降优化器和(1 1)-es这两种搜索策略分别对160幅亮血图像与相应扫描层面的160幅增强黑血图像进行配准，结果显示，未使用优化器进行配准的图像中，增强黑血图像与亮血图像并未重合，且阴影较多；当使用梯度下降优化器配准图像时，虽然较未使用前配准效果较好，但在脑灰质处仍出现了明显的不重合现象；而使用(1 1)-es优化器的图像中，配准结果精确，图像中不重合的阴影部分完全消失。从归一化互信息nmi、归一化互相关系数ncc与算法耗时time进行评价，(1 1)-es的配准图像效果显示更清晰，优于梯度下降优化器；因此本发明实施例使用(1 1)-es作为预定搜索策略。
[0216]
(2)统一扫描区域：
[0217]
从增强黑血图像中，提取与第一亮血图像的扫描范围相同的区域内容，形成第一黑血图像。具体过程参见s212步骤，在此不再赘述。
[0218]
参见图3，图3为本发明实施例的颅内血管磁共振图像的待配准区域示意图，其中左图为预配准后的第一亮血图像，右图为增强黑血图像，方框为增强黑血图像中待提取区域。这个区域包含了在颅内血管磁共振图像中，亮血序列以及黑血序列的共同扫描范围，通过确定待提取区域，能够更快速地关注到有用信息。
[0219]
(二)，在预处理之后，可以采用基于互信息和图像金字塔的配准方法对第一亮血图像和第一黑血图像进行图像配准，具体参见前文步骤s22～s27相关内容，在此不再赘述。具体包括：
[0220]
①
基于下采样处理，由第一亮血图像得到亮血高斯金字塔，由第一黑血图像得到黑血高斯金字塔；
[0221]
其中，亮血高斯金字塔和黑血高斯金字塔中包括从下至上分辨率依次变小的4个图像；如图4(a)所示，图4(a)为本发明实施例的颅内血管磁共振图像的亮血高斯金字塔和黑血高斯金字塔。
[0222]
②
基于上采样处理，利用亮血高斯金字塔得到亮血拉普拉斯金字塔，利用黑血高斯金字塔得到黑血拉普拉斯金字塔；
[0223]
其中，亮血拉普拉斯金字塔和黑血拉普拉斯金字塔中包括从下至上分辨率依次变小的3个图像；亮血拉普拉斯金字塔和黑血拉普拉斯金字塔的生成过程参见前文s23，在此不再赘述。如图4(b)所示，图4(b)为本发明实施例的颅内血管磁共振图像的亮血拉普拉斯金字塔和黑血拉普拉斯金字塔。图像显示使用了伽马矫正实现更清晰的效果，伽马值为0.5。
[0224]
③
对亮血拉普拉斯金字塔和黑血拉普拉斯金字塔中对应层的图像进行配准，得到配准的亮血拉普拉斯金字塔；
[0225]
在该步骤中，使用互信息作为两幅图的相似性度量，选用(1 1)-es作为搜索策略，
在每一次图像配准进行坐标变换后，循环迭代计算两幅图的互信息，直至互信息达到最大，图像配准完成。具体过程参见前文s24，在此不再赘述。
[0226]
结果如图5所示，图5为本发明实施例的颅内血管磁共振图像的拉普拉斯金字塔图像的配准结果，左图为黑血拉普拉斯金字塔中的参考图像，中图为亮血拉普拉斯金字塔中已配准好的图像，右图为左、中两幅图像直接叠加后的效果图，叠加图像显示采用蒙太奇效果，使用伪彩色透明处理增强黑血图像与亮血图像，其中紫色为增强黑血拉普拉斯金字塔图像，绿色为亮血拉普拉斯金字塔图像(附图图像是原图经过灰度处理的图像，颜色未示出)。
[0227]
④
利用配准的亮血拉普拉斯金字塔作为叠加信息，对亮血高斯金字塔和黑血高斯金字塔中各层图像进行自上而下的配准，得到配准的亮血高斯金字塔；
[0228]
该步骤参见前文s25，基于互信息的高斯金字塔图像配准具体步骤如图6所示，图6为本发明实施例的颅内血管磁共振图像基于互信息的高斯金字塔图像配准步骤示意图。首先对顶层的低分辨率的黑血高斯图像和顶层的低分辨率亮血高斯图像进行基于互信息的配准；接着将已配准好的亮血高斯图像进行上采样操作，并与保留了高频信息，且根据上述操作已配准好的对应层的亮血拉普拉斯图像相加，作为下一层亮血高斯图像；接着将上述操作得到的亮血高斯图像作为输入图像，再与对应层的黑血高斯图像进行配准，重复以上操作，直至完成底层高斯金字塔图像的高分辨率配准。
[0229]
在基于互信息的高斯金字塔图像配准中，需要对每一层亮血高斯图像和黑血高斯图像以归一化互信息作为相似性度量进行配准，通过循环迭代计算两幅图的nmi，直至nmi达到最大。其中当迭代次数太小时，不能完成图像的准确配准，但是当迭代次数太大时，计算量会急剧增加，图7为本发明实施例的不同迭代次数下的归一化互信息，当高斯金字塔中分辨率最高的底层图像配准达到最大nmi值且数据稳定时，停止迭代。
[0230]
为验证本发明实施例的基于互信息和图像金字塔的图像配准方法的有效性与实用性进行了对比实验，使用五位患者的颅内血管磁共振图像；同时选取仅使用dicom图像方位标签信息进行配准的算法，以及基于互信息度量的配准算法，与本发明实施例的基于互信息和图像金字塔的配准方法进行对比，其中基于互信息度量的算法是通过多参数优化方法寻找参考图像和浮动图像间的最佳变换，使得两幅图像的互信息值最大，并没有使用到图像金字塔算法。实验平台是matlab r2016b。针对实验的图像配准结果，采用定性分析与定量分析相结合。在定性分析方面，由于多模态医学图像间存在较大的灰度差异，将配准图像与参考图像相减得到的差值图像无法有效地反映出多模态医学图像的配准结果，因此本发明实施例通过将配准图像与参考图像进行重叠，获得可以反映出配准图像和参考图像对齐程度的彩色重叠图像，通过彩色重叠图像对多模态配准算法的配准效果进行定性分析，图8显示了多模态颅内血管磁共振图像的配准结果，图8为多种配准方法的颅内血管磁共振图像的配准结果。其中，(a)为参考图像；(b)为浮动图像；(c)为基于图像方位标签信息的重叠图像；(d)为基于互信息度量的重叠图像；(e)为本发明实施例的基于互信息和图像金字塔的图像配准方法的重叠图像。其中附图均为原图的灰度图，彩色未示出。在定量分析方面，采用归一化互相关系数ncc，归一化互信息nmi作为评价指标，当归一化互相关系数ncc和归一化互信息nmi的值越大时，表示图像配准精度越高，表1显示了不同配准算法的评价指标结果分析，受篇幅所限仅示出患者a数据。
[0231]
表1不同配准方法的结果分析
[0232][0233]a中的值是基于患者多幅图像配准的评价指标平均值
±
均方误差
[0234]
从图8的重叠图像可以明显看出，基于互信息度量的方法出现了较大的配准偏移，可能是因为仅使用基于互信息度量的方法容易陷入局部最优值，而非全局最优值；基于图像方位标签信息的配准效果表现同样欠佳，图像出现了部分不重叠情况；而基于互信息和图像金字塔的配准方法的图像效果表现良好，图像显示更加清晰，图像几乎完全重叠在一起。由表1可知，从ncc与nmi这两个评价指标来看，相较于只使用dicom图像的方位标签信息的配准算法，以及基于互信息度量的配准算法，本发明实施例提出的基于互信息和图像金字塔的配准方法在配准精度上有所提高，可以良好地处理多模态颅内血管磁共振图像的配准。
[0235]
⑤
基于配准的亮血高斯金字塔得到该亮血图像对应的配准后亮血图像；
[0236]
获取配准的亮血高斯金字塔中的底层图像作为配准后亮血图像，将配准后亮血图像和对应的增强黑血图像作为配准后图像对。
[0237]
⑥
由k个亮血图像分别对应的配准后亮血图像，得到配准后亮血图像组。
[0238]
本发明实施例中使用基于互信息与图像金字塔的图像配准方法，不仅考虑灰度信息的相关性，而且利用高斯金字塔提高配准效率，使图像从低分辨率到高分辨率，逐层提高配准精度。
[0239]
步骤三，利用配准后亮血图像组，对增强黑血图像组中的增强黑血图像进行流空伪影消除操作，得到伪影消除增强黑血图像组；具体参见前文s3步骤。结果如下：
[0240]
如图9所示，图9为本发明实施例的灰度线性变换结果图。图9中，左图为配准后亮血图像，右图为其经过灰度线性变换后的结果图，可以看到右图中血液部分与周围像素相比对比度明显增强。采用最大类间方差法otsu得到亮血特征图，结果如图10所示，图10为本发明实施例的图像二值化结果图，图10中左图为对比增强亮血图像，右图为其经过阈值提取后的血液信息。可以看到，右图中显示为亮色的部分仅为与血液相关的信息。图11为本发明实施例针对颅内血管的流空伪影消除结果。图11中左图为增强黑血图像原图，右图为流空伪影消除后的增强黑血图像，箭头所示处出现了流空伪影，对比可见流空伪影的消除效果较为明显。
[0241]
步骤四，利用配准后亮血图像组建立血液三维模型；
[0242]
步骤五，利用配准后亮血图像组建立血液边界扩展的血管三维模型；
[0243]
步骤六，基于伪影消除增强黑血图像组和黑血图像组，得到造影增强三维模型；
[0244]
步骤七，基于血液三维模型、血管三维模型、造影增强三维模型，得到血管造影增强三维模型。
[0245]
步骤八，获取血管造影增强三维模型中各段血管的表征血管狭窄程度的目标参数
的数值，并利用各段血管的目标参数的数值对血管造影增强三维模型进行标记，得到血管造影增强三维狭窄化分析模型。
[0246]
步骤九，将血管造影增强三维狭窄化分析模型进行显示。
[0247]
关于步骤四至步骤九的具体过程不再赘述。具体结果如下，参见图12，图12为本发明实施例针对颅内血管的血液三维模型效果图、血管三维模型效果图和造影增强三维模型效果图；图12中从左至右分别为血液三维模型效果图、血管三维模型效果图和造影增强三维模型效果图。其中各图均进行了灰度处理，在实际中，可以以不同颜色进行显示，比如血管三维模型效果图和造影增强三维模型效果图分别以蓝色和红色进行显示。
[0248]
图13为本发明实施例针对颅内血管的血管造影增强三维模型效果图。其中白色线圈内的亮色部分为出现了造影增强的颅内血管部位，即该处可能出现了颅内动脉粥样硬化的病症或者易损性斑块，其余部分为没有出现造影增强的血管部位，在实际中，图13中可以用不同颜色进行区分，比如蓝色为没有出现造影增强的血管部位，红色处为出现了造影增强的血管部位。并且该血管造影增强三维模型可以实现旋转、放大缩小等基本功能，从而辅助医生定位病灶区域，作出更精准的判断。
[0249]
图14为本发明实施例的颅内血管的血管造影增强三维狭窄化分析模型效果图。其中左图为狭窄率标记效果，右图为扁平度标记效果。在实际中，模型上显示有不同的颜色，可以区分狭窄化程度，比如血管较细的部位为暖色调，最狭窄的部位为红色，血管较粗的部位为冷色调，最粗的部位显示为绿色等，白色箭头所示处为颅内血管的突变狭窄，不同色调的颜色显示可以更加直观地表现出血管的狭窄化。附图中是经过灰度处理的效果，颜色未示出。
[0250]
由于医生习惯于观察切平面的二维医学图像，本发明实施例可以在提供模拟化三维血管狭窄分析模型的同时提供三个方位的二维切面图，即可以将该模拟化三维血管狭窄分析模型中的每一个点对应的当前点所在冠状面、矢状面、轴状面图像显示出来。参见图15，图15为本发明实施例的颅内血管的血管造影增强三维狭窄化分析模型和切面图显示效果图。图15中血管的暖色调处可能出现了血管狭窄化，冷色调处没有明显的血管狭窄化现象，且图像右边的三幅二维图像从上至下分别为当前点所在的轴状面、矢状面以及冠状面成像；在显示该模拟化三维血管狭窄分析模型时，还可以使用诸如红绿蓝三种颜色的点，实现两个点测距，三个点测量角度的功能，并在显示屏左下方显示出来，在显示屏右下方显示当前选取模型的体积大小。以便于医生获取颅内血管更为详细的数据。
[0251]
参见图16，图16是本发明实施例提供的颅内血管的血管造影增强三维狭窄化分析模型裸眼3d全息可视化图，图中，前视、后视、左视、右视四个视图组合在一起，实现裸眼3d全息可视化。参见图17，图17是本发明实施例提供的对颅内血管的血管造影增强三维狭窄化分析模型的裸眼3d全息显示结果进行手势识别的示意图。请参见图18，图18是本发明实施例提供的颅内血管的血管造影增强三维狭窄化分析模型的3d打印结果图。图16-图18提供的显示方法，均是为了进一步将得到的颅内血管的血管造影增强三维狭窄化分析模型更直观的进行展示，以便于医生对颅内病灶判断时代入感更强。
[0252]
本发明实施例所提供的方案中得到的血管造影增强三维狭窄化分析模型，实现了颅内血管的三维可视化，无需医生通过想象力还原血管组织结构及病症特征等，可以提供形象的血管三维空间信息，便于直观观测，便于定位与显示狭窄的病灶区域。能够在临床应
用上，简便、快速、直观地获得血管的真实信息和关于颅内血管狭窄程度的分析数据。
[0253]
注：本发明实施例中的患者实验数据均来源于陕西省人民医院，图像可用作一般的科学研究。
[0254]
以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。