肝脏医学图像动脉区域分割方法及系统与流程

本发明涉及一种医学影像处理系统，特别涉及一种肝脏医学图像的处理系统。

背景技术：

通常，肝病患者进行外科手术需要术前采集肝脏ct造影图像，显示出肝脏动脉的解剖结构，定位病灶供血的肝脏动脉，尽量减小手术对肝脏组织的影响。目前，动脉区域的分割是根据医生临床以及肉眼观察判断进行估计，往往由于不够准确而影响手术效果或造成手术失败。因此，医疗领域一直都在努力研发一种能够自动分割肝脏动脉区域的系统。

如中国专利申请201510313660.x号公开的一种基于ct图像的肝动脉分割方法，其包括：载入肝动脉ct图像数据和对应的肝掩模数据，并定位动脉种子点；对ct图像数据进行分割，分别得到包含脊椎骨及其相连肋骨和肝动脉的第一二值体数据以及分别得到包含脊椎骨及其相连肋骨和肝动脉的第二二值体数据；在第二二值体数据中剔除与第一二值体数据中相同的部分，获得第三二值体数据；利用干动脉种子点在第三二值体数据上进行三维区域增长，得到肝动脉分割结果。然而，该分割方法对于分割区域中存在血管交叉或与其他组织附着的情况缺乏精细处理，导致获得的肝动脉分割图像精确度较低，不利于病灶识别。

又如中国专利申请201410132571.0号公开的一种冠状动脉分割方法及装置，包括如下步骤：(1)确定原始图像中主动脉的位置，并对原始图像进行增强滤波，得到增强滤波图像；(2)根据主动脉位置、原始图像以及增强滤波图像的灰度，得到左、右冠脉的起始区域；(3)基于至少一项预设分割阈值分割增强率波图像，获取预设阈值对应的候选区域，并根据该候选区域和左、右冠脉的起始区域，确定左、右冠脉的候选区域；(4)根据主动脉和左、右动脉候选区域的位置关系，分别对左、右冠脉候选区域进行筛选，筛选后和起始区域进行连接，得到左、右冠脉。然而，该方法分割得到的冠状动脉存在血管交叉或与其他组织附着的问题。

因此，提供一种能够快速高效准确的在肝脏医学图像中分割肝脏动脉血管的肝脏医学图像动脉区域分割方法及系统法成为业内急需解决的问题

技术实现要素：

本发明的目的是提供一种肝脏医学图像动脉区域分割方法及系统，该方法快速高效，解决了血管交叉或与其他组织附着的问题，分割处理后得到的肝脏动脉分割图像精确度高。

为了实现上述目的，本发明一个方面提供了一种肝脏医学图像动脉区域分割方法，包括：(1)去除肝脏ct图像上的肌肉组织及肝组织，露出肝脏血管轮廓；(2)选取肝脏血管轮廓上的距离最远的呈对角线的两点，以选中的两点为对角线建立长方体，以长方体所包含的数据空间作为运行区域；(3)通过hessian矩阵函数计算运行区域中每个点的hessian矩阵函数数值，获得函数数值集合，在函数数值集合中选择指定步长提取种子点，其中，指定步长设定为小于5个体素；(4)计算步骤(3)中获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率，参考已设定的两个主曲率的参数范围，剔除主曲率不属于参数范围内的种子点，以剔除肝脏血管轮廓中的非管状组织，保留主曲率属于参数范围内的种子点；(5)基于拓扑细化法追踪步骤(4)保留的种子点，获得血管的中心线，比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界；以及(6)确定血管数目，根据中心线到血管边界的距离测算血管半径，根据血管半径设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模，获取肝脏动脉分割图像。

可选择地，步骤(5)中，两个主曲率的商值均为0。

可选择地，该方法进一步包括步骤(1)之前，对肝脏ct图像进行增强处理。

其中，hessian矩阵，通常称为黑塞矩阵，又译作海森矩阵、海瑟矩阵、海塞矩阵等，是一个多元函数的二阶偏导数构成的方阵，描述了函数的局部曲率。黑塞矩阵最早于19世纪由德国数学家ludwigottohesse提出，并以其名字命名。在工程实际问题的优化设计中，常常将目标函数在某点邻域展开成泰勒多项式来逼近原函数，此时函数在某点泰勒展开式的矩阵形式中会涉及到黑塞矩阵。

本发明的另一个方面在于提供一种肝脏医学图像动脉区域分割系统，包括：数据读取模块，其用于读取肝脏ct图像；运行区域划分模块，其与数据读取模块通信连接，用于在肝脏ct图像上选择距离最远的呈对角线的两点，以选中的两点为对角线建立长方体，以长方体中所包含的数据空间作为运行区域；函数计算模块，其与运行区域划分模块通信连接，用于通过hessian矩阵函数计算运行区域中每个点的hessian矩阵函数数值，获得函数数值集合，在函数数值集合中选择指定步长提取种子点，其中，指定步长设定为小于5个体素；种子点确认模块，其与函数计算模块通信连接，用于计算函数计算模块获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率，以这两个主曲率为上下限设定参数范围，剔除主曲率不属于参数范围内的种子点，保留主曲率属于参数范围内的种子点；中心线提取模块，其与种子点确认模块通信连接，用于基于拓扑细化法追踪种子点确认模块保留的种子点，获得血管的中心线，比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界；以及三维绘制模块，其与中心线提取模块通信连接，用于在确定血管数目后，根据中心线到血管边界的距离测算血管半径，根据血管半径设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模。

可选择地，进一步包括与数据读取模块及运行区域划分模块通信连接的图像预处理模块，图像预处理模块包括通信连接的图像平滑单元和图像增强单元，其中，图像平滑单元用于通过空域法在空间域中对每幅肝脏ct图像的像素灰度值直接进行运算处理，以滤除每幅肝脏ct图像中的噪声；图像增强单元用于尖锐化增强处理经过图像平滑处理单元处理后的每幅肝脏ct图像。

可选择地，种子点确认模块还用于计算动脉血管的形状指数。

可选择地，种子点确认模块包括通信连接的法曲率计算单元及种子点保留单元，其中，法曲率计算单元用于计算函数计算模块获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率；种子点保留单元用于参考已设定的两个主曲率的参数范围，剔除主曲率不属于参数范围内的种子点，保留主曲率属于参数范围内的种子点。

可选择地，中心线提取模块包括通信连接的中心点探测单元以及血管边界探测单元，其中，中心点探测单元用于基于拓扑细化法追踪种子点确认模块保留的种子点，获得血管的中心线；血管边界探测单元用于比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界。

可选择地，三维绘制模块包括：通信连接的血管数目计算单元、血管半径测量单元、以及区域合并绘制单元，其中，血管数目计算单元用于对肝脏ct图像进行区域裁剪确定血管数目；血管半径测量单元用于根据中心线到血管边界的距离测算血管半径；区域合并绘制单元用于根据血管半径，设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模。

可选择地，还包括与三维绘制模块通信连接的肝脏动脉数据库，其用于存储肝脏来源信息、对应肝脏的肝脏ct图像、以及对应的肝脏动脉分割图像。

可选择地，肝脏来源信息至少包含上传的对应肝脏ct图像的患病类型、患者性别、患者年龄、生活地区以及就诊医院等来源信息。

可选择地，肝脏来源信息还可以包含：患者生活习惯、生化检查信息、患者典型症状、体征、医学图像、影像诊断结果、治疗方案、不良反应及主治医师等。

可选择地，还包括与数据库通信连接的检索模块，其用于在肝脏动脉数据库中检索目标肝脏的肝脏ct图像及对应的肝脏动脉分割图像。

可选择地，函数计算模块设定的参数还可以包括：动脉血管高斯曲率、平均曲率以及弯曲度等相关参数。

可选择地，主曲率和形状指数都可以描述各类物体的不同形状，理想管状结构的两个主曲率的商值约为0，形状指数约等于0.75。血管的形状指数范围是0.65～0.8。

可选择地，主曲率的参数范围可为0.5～0.7。

可选择地，血管数目计算单元设定肝脏ct图像上每根血管有且只有一个亮区，即每根血管通过各切片只有一次，血管平稳的改变方向，所有血管在第一次切片开始，即没有新的神经血管出现在其余的切片。

本发明的有益效果是：(1)、本发明以中心线、长轴和短轴建立椭圆血管模型，由于各个血管不是均匀的，不同血管的不同部位半径不同，但是对于每根血管来说，半径的变化是连续的，避免将血管设定为圆形造成误差使获得的分割图像更接近真实血管；(2)、本发明除去肝脏ct图像中的非管状组织，三维绘制模块进一步确定动脉血管数目，避免因三维空间中血管发生附着或交叉导致提取到的中心线丢失影响分割准确度；(3)本发明的分割系统通过函数计算模块设置血管参数，例如主曲率和形状指数等，可以刻画动脉血管曲面的波动、凹凸等微小形状变化。

附图说明

图1为本发明的肝脏医学图像动脉区域分割系统的构造示意图。

图2为本发明的肝脏医学图像动脉区域分割方法的流程示意图。

具体实施方式

请参照图1，根据本发明的一种非限制性实施方式，提供了一种肝脏医学图像动脉区域分割系统，包括：数据读取模块10、运行区域划分模块20、函数计算模块30、种子点确认模块40、中心线提取模块50、以及三维绘制模块60。

其中，数据读取模块10用于读取肝脏ct图像。

运行区域划分模块20用于在肝脏ct图像上选择距离最远的呈对角线的两点，以选中的两点为对角线建立长方体，以长方体中所包含的数据空间作为运行区域。

函数计算模块30用于通过hessian矩阵函数计算运行区域中每个点的hessian矩阵函数数值，获得函数数值集合，在函数数值集合中选择指定步长提取种子点，其中，指定步长设定为小于5个体素。

种子点确认模块40用于计算函数计算模块获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率，参考已设定的两个主曲率的参数范围，剔除主曲率不属于参数范围内的种子点，保留主曲率属于参数范围内的种子点。

中心线提取模块50用于基于拓扑细化法追踪种子点确认模块保留的种子点，获得血管的中心线，比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界。

三维绘制模块60用于在确定血管数目后，根据中心线到血管边界的距离测算血管半径，根据血管半径设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模。

作为一种可替代实施方式，还包括与数据读取模块10及运行区域划分模块20通信连接的图像预处理模块70，图像预处理模块70包括图像平滑单元710和图像增强单元720。其中，图像平滑单元710用于通过空域法在空间域中对每幅肝脏ct图像的像素灰度值直接进行运算处理，从而滤除每幅肝脏ct图像中的噪声。图像增强单元720用于尖锐化增强处理经过图像平滑处理单元处理后的每幅肝脏ct图像。

作为另一种可替代实施方式，种子点确认模块40包括通信连接的法曲率计算单元410及种子点保留单元420。其中，法曲率计算单元410用于计算函数计算模块获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率。种子点保留单元420用于参考已设定的两个主曲率的参数范围，剔除主曲率不属于参数范围内的种子点，保留主曲率属于参数范围内的种子点。

中心线提取模块50包括中心点探测单元510以及血管边界探测单元520。其中，中心点探测单元510用于基于拓扑细化法追踪种子点确认模块保留的种子点，获得血管的中心线。血管边界探测单元520用于比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界。

三维绘制模块60包括血管数目计算单元610、血管半径测量单元620、以及区域合并绘制单元630。其中，血管数目计算单元610用于对肝脏ct图像进行区域裁剪确定血管数目。血管半径测量单元620用于根据中心线到血管边界的距离测算血管半径。区域合并绘制单元630用于根据血管半径，设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模。

作为又一种可替代实施方式，肝脏医学图像动脉区域分割系统还包括肝脏动脉数据库80及检索模块90。其中，肝脏动脉数据库80用于存储肝脏来源信息、对应肝脏的肝脏ct图像、以及对应的肝脏动脉分割图像。检索模块90用于在肝脏动脉数据库80中检索目标肝脏的肝脏ct图像及对应的肝脏动脉分割图像。

下面，结合图2，详细叙述本发明的肝脏医学图像动脉区域分割方法。如图2所示，首先，在步骤s1中，去除肝脏ct图像上的肌肉组织及肝组织，露出肝脏血管轮廓。接着，在步骤s2中，选取肝脏血管轮廓上的距离最远的呈对角线的两点，以选中的两点为对角线建立长方体，以长方体所包含的数据空间作为运行区域。然后，在步骤s3中，通过hessian矩阵函数计算运行区域中每个点的hessian矩阵函数数值，获得函数数值集合，在函数数值集合中选择指定步长提取种子点，其中，指定步长设定为小于5个体素。随后，在步骤s4中，计算步骤s3中获得的每个种子点的法曲率，将法曲率的最大值和最小值设定为种子点的两个主曲率，参考已设定的两个主曲率的参数范围，剔除主曲率不属于参数范围内的种子点，以剔除肝脏血管轮廓中的非管状组织，保留主曲率属于参数范围内的种子点；在步骤s5中，基于拓扑细化法追踪步骤s4保留的种子点，获得血管的中心线，比较垂直于拓扑细化法的追踪方向的像素的灰度值，设定血管边界的灰度值得到血管边界。最后，在步骤s6中，确定血管数目，根据中心线到血管边界的距离测算血管半径，根据血管半径设定血管的长轴和短轴，以椭圆来描绘血管截面形状建模。

尽管在此已详细描述本发明的优选实施方式，但要理解的是本发明并不局限于这里详细描述和示出的具体结构和步骤，在不偏离本发明的实质和范围的情况下可由本领域的技术人员实现其它的变型和变体。