三维血管图像处理方法、装置和计算机程序产品制造方法及图纸

本申请涉及一种三维血管图像处理方法、装置和计算机程序产品，其中，该三维血管图像处理方法包括：获取待处理的初始三维血管图像，计算初始三维血管图像中的每个体素的目标血管概率，得到第一血管概率图像；采用预设的分割阈值，对第一血管概率图像进行二值分割处理，得到初始分割图像；在初始分割图像的第一值的表示区域中，将尺寸小于预设值的连通区域以第二值表示，得到目标分割图像；基于目标分割图像，对初始三维血管图像进行伪影抑制处理，得到目标三维血管图像。通过本申请，解决了现有三维血管图像中，强度较大的伪影与小血管具有接近的灰度值，导致用户无法通过设置合适的窗宽窗位，在保留小血管的前提下不显示伪影的问题。的问题。的问题。

【技术实现步骤摘要】
三维血管图像处理方法、装置和计算机程序产品


[0001]本申请涉及医疗图像处理领域，特别是涉及一种三维血管图像处理方法、装置和计算机程序产品。

技术介绍

[0002]三维血管图像可以展示血管的三维空间形态，但是由于各种原因，在三维血管图像中存在各种伪影。示例性地，三维血管造影(3D DSA)可以重建被造影剂填充的血管的三维图像，允许用户以任意角度观察血管树的空间结构。在三维血管造影图像中，通常包括：运动伪影、不一致性伪影和欠采样伪影等。
[0003]通常在显示三维血管图像时，通过设置合适的窗宽窗位可以过滤掉部分伪影。其原理在于，伪影相较于血管，其灰度值较低。因此，可以通过设置合适的窗宽窗位，不显示图像中灰度值较低的部分。但是，伪影的灰度值与其自身强度相关，强度越大其灰度值越大；同时，小血管相较于大血管，其灰度值较低。
[0004]在三维血管图像的成像过程中，若是成像条件比较好，在患者无运动、造影剂注射均匀时，三维血管图像中伪影的强度会较小，则可以较为有效地通过设置合适的窗宽窗位去除伪影；若是成像条件比较差，出现扫描时血管未被造影剂充盈、扫描末期造影剂已停止注射等情况时，三维血管图像中伪影的强度会较大。
[0005]当三维血管图像中伪影的强度较大时，其灰度值与小血管的灰度值比较接近，甚至会高于小血管的灰度值。此时若想通过设置较高的窗宽窗位，使得图像中的伪影不显示，则必然会导致图像中的小血管无法显示，从而极大地牺牲血管细节结构。
[0006]针对三维血管图像中，强度较大的伪影与小血管具有接近的灰度值，导致用户无法通过设置合适的窗宽窗位，在保留小血管的前提下不显示伪影的问题，目前还需要一些有效的解决方案。

技术实现思路

[0007]在本专利技术中提供了一种三维血管图像处理方法、装置和计算机程序产品，以解决现有三维血管图像中，强度较大的伪影与小血管具有接近的灰度值，导致用户无法通过设置合适的窗宽窗位，在保留小血管的前提下不显示伪影的问题。
[0008]第一个方面，在本专利技术中提供了一种三维血管图像处理方法，所述方法包括：
[0009]获取待处理的初始三维血管图像，计算所述初始三维血管图像中的每个体素的目标血管概率，得到第一血管概率图像；
[0010]采用预设的分割阈值，对所述第一血管概率图像进行二值分割处理，以第一值表示高于所述分割阈值的图像部分，以第二值表示低于所述分割阈值的图像部分，得到初始分割图像；
[0011]在所述初始分割图像的所述第一值的表示区域中，将尺寸小于预设值的连通区域以所述第二值表示，得到目标分割图像；
[0012]基于所述目标分割图像，对所述初始三维血管图像进行伪影抑制处理，得到目标三维血管图像。
[0013]在其中的一些实施例中，所述基于所述目标分割图像，对所述初始三维血管图像进行伪影抑制处理，得到目标三维血管图像，包括：
[0014]根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，得到所述目标三维血管图像。
[0015]在其中的一些实施例中，所述根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，包括：
[0016]将所述第一血管概率图像进行概率反转，得到所述第一血管概率图像的反转图像；
[0017]根据所述目标分割图像，在所述反转图像中，确定对应所述第一值的血管区域和对应所述第二值的非血管区域，并降低所述非血管区域的血管概率，使得所述非血管区域的血管概率低于所述血管区域的血管概率，得到第二血管概率图像；
[0018]根据所述第一血管概率图像和所述第二血管概率图像，调整所述初始三维血管图像中各体素的灰度值。
[0019]在其中的一些实施例中，所述根据所述第一血管概率图像和所述第二血管概率图像，调整所述初始三维血管图像中各体素的灰度值，包括：
[0020]按照第一权重叠加所述第一血管概率图像和所述第二血管概率图像，得到第三血管概率图像；
[0021]根据所述第三血管概率图像，调整所述初始三维血管图像中各体素的灰度值。
[0022]在其中的一些实施例中，所述降低所述非血管区域的血管概率，包括：
[0023]将所述非血管区域的血管概率置零。
[0024]在其中的一些实施例中，所述根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，包括：
[0025]根据所述目标分割图像，在所述第一血管概率图像中确定对应所述第二值的非血管区域，并将所述非血管区域的血管概率置零，得到第四血管概率图像；
[0026]根据所述第四血管概率图像，调整所述初始三维血管图像中各体素的灰度值。
[0027]在其中的一些实施例中，所述根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，得到所述目标三维血管图像，包括：
[0028]根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，得到中间三维血管图像；
[0029]按照第二权重叠加所述中间三维血管图像和所述初始三维血管图像，得到所述目标三维血管图像。
[0030]在其中的一些实施例中，所述计算所述初始三维血管图像中的每个体素的目标血管概率，包括：
[0031]分别采用若干不同标准差对所述初始三维血管图像进行高斯滤波处理，并计算高斯滤波处理后的所述初始三维血管图像中每个所述体素的血管概率；
[0032]针对每个所述体素，确定所述体素对应于不同标准差的多个血管概率中的最大血管概率，根据所述最大血管概率确定所述目标血管概率。
[0033]第二个方面，在本专利技术中提供了一种三维血管图像处理装置，所述装置包括：
[0034]图像获取模块，用于获取待处理的初始三维血管图像，计算所述初始三维血管图像中的每个体素的目标血管概率，得到第一血管概率图像；
[0035]图像分割模块，用于采用预设的分割阈值，对所述第一血管概率图像进行二值分割处理，以第一值表示高于所述分割阈值的图像部分，并以第二值表示低于所述分割阈值的图像部分，得到初始分割图像；
[0036]区域移除模块，用于在所述初始分割图像的所述第一值的表示区域中，将尺寸小于预设值的连通区域以所述第二值表示，得到目标分割图像；
[0037]伪影抑制模块，用于基于所述目标分割图像，对所述初始三维血管图像进行伪影抑制处理，得到目标三维血管图像。
[0038]第三个方面，在本专利技术中提供了一种电子装置，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一个方面所述的三维血管图像处理方法。
[0039]第四个方面，在本专利技术中提供了一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述第一个方面所述的三维血管图像处理方法。
[0040]第五个方面，在本专利技术中提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现上述第一个方面所述的三维血管图像处理本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种三维血管图像处理方法，其特征在于，所述方法包括：获取待处理的初始三维血管图像，计算所述初始三维血管图像中的每个体素的目标血管概率，得到第一血管概率图像；采用预设的分割阈值，对所述第一血管概率图像进行二值分割处理，以第一值表示高于所述分割阈值的图像部分，以第二值表示低于所述分割阈值的图像部分，得到初始分割图像；在所述初始分割图像的所述第一值的表示区域中，将尺寸小于预设值的连通区域以所述第二值表示，得到目标分割图像；基于所述目标分割图像，对所述初始三维血管图像进行伪影抑制处理，得到目标三维血管图像。2.根据权利要求1所述的三维血管图像处理方法，其特征在于，所述基于所述目标分割图像，对所述初始三维血管图像进行伪影抑制处理，得到目标三维血管图像，包括：根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，得到所述目标三维血管图像。3.根据权利要求2所述的三维血管图像处理方法，其特征在于，所述根据所述目标分割图像，调整所述初始三维血管图像中各体素的灰度值，包括：将所述第一血管概率图像进行概率反转，得到所述第一血管概率图像的反转图像；根据所述目标分割图像，在所述反转图像中，确定对应所述第一值的血管区域和对应所述第二值的非血管区域，并降低所述非血管区域的血管概率，使得所述非血管区域的血管概率低于所述血管区域的血管概率，得到第二血管概率图像；根据所述第一血管概率图像和所述第二血管概率图像，调整所述初始三维血管图像中各体素的灰度值。4.根据权利要求3所述的三维血管图像处理方法，其特征在于，所述根据所述第一血管概率图像和所述第二血管概率图像，调整所述初始三维血管图像中各体素的灰度值，包括：按照第一权重叠加所述第一血管概率图像和所述第二血管概率图像，得到第三血管概率图像；根据所述第三血管概率图像，调整所述初始三维血管图像中各体素的灰度值。5.根据权利要求3所述的三维血管图像处理方法，其特征在于，所述降低所述非血管区域的血管概率，包括：将所述非血管区域的血管概率置零。6.根据权利要求2所述的三维血管图像处...

【专利技术属性】
技术研发人员：宋聿晗，郭立祺，
申请(专利权)人：深圳市联影高端医疗装备创新研究院，
类型：发明
国别省市：