一种基于医学图像的器官骨架提取方法
【专利摘要】本发明提供一种基于医学图像的器官骨架提取方法，在利用DICOM格式医学图像的基础上，以提取一个完整的器官骨架为目标，围绕基于结构张量的局部结构提取、基于K-means模糊聚类的点云模型分割、基于点云的骨架提取方法展开研究，首先从DICOM格式医学图像中读取三维数据，并根据三维数据的局部结构特征对三维数据进行特征提取得到点云，然后对得到的点云通过K-means聚类进行模型分割，最终对得到的分割点云进行骨架提取。
【专利说明】一种基于医学图像的器官骨架提取方法
【技术领域】
[0001] 本发明涉及一种基于医学图像的器官骨架提取方法。
【背景技术】
[0002] 随着计算机在各行业的广泛应用，人们开始不满足于仅能显示二维的图像，更希 望计算机能表现出真实感强三维世界。计算机图形学、虚拟现实等技术的发展，三维数据成 为继图像、音频、视频以来又一重要的多媒体数据形式。三维数据就是使用各种采集仪采集 得到的，它记录了有限体表面在离散点上的各种物理参量。它包括最基本信息是物体的各 离散点三维坐标，其它的可以包括物体表面颜色、透明度、纹理特征等。由于三维模型的信 息量很大，这也使得在许多应用中出现占用内存空间过大、运行计算负载过重，所以需要一 种"紧凑的"方式来尽可能完整、全面地表示描述三维模型的结构特征信息。其中最常用一 种简化表示方式就是使用一维曲线，一般称为骨架。利用物体的骨架来描述对象是一种既 能强调物体的结构特征，也能提高内存使用效率与数据压缩的好方法。
[0003] 为了找到"紧凑的"一维曲线（骨架），本发明提出了一种基于医学图像的器官骨 架提取方法，以提取完整的器官骨架结构为目标，围绕三维数据局部结构、模糊分类、骨架 提取展开研究。首先通过对医学文件读取得到CT的原始扫描数据，并对三维数据计算其二 阶黑塞矩阵的特征值，选取特征值满足一定条件的三维数据，然后结合K-means模糊聚类， 对得到的三维数据进行分类，最后通过Ll-中值定义和最小化能量泛函，并加入正则化得 到一个完整的器官骨架结构。
【发明内容】
[0004] 本发明解决的技术问题是：提供了一种基于医学图像的器官骨架提取方法，且可 以直接对原始扫描数据进行骨架提取，适用度广。
[0005] 本发明采用的技术方案为：一种基于医学图像的器官骨架提取方法，主要包括以 下三个步骤：
[0006] 步骤（1)、读取三维医学图像数据：通过对DICOM格式的医学CT数据读取，得到医 学器官的三维数据，作为整个算法原始输入数据。
[0007] 步骤（2)、提取三维数据的局部特征结构：通过对得到三维数据固有的结构特征 进行分析来提取局部结构；利用三维数据的黑塞矩阵来刻画局部几何结构，并通过计黑塞 矩阵的特征值进行组合来提取合适的局部几何结构，最大特征值所对应的特征向量表征了 从一种材质过度到另一种材质时最为迅速的变化方向，而最小特征值所对应的方向则指明 了某种材质在图像内部流动的方向；因此可以对三维图像内部的材质属性变化进行定量刻 画；得到的三维离散数据（也就是点云模型）作为整个骨架提取的输入数据；
[0008] 步骤（3)、点云模型分割：基于步骤（1)中的数据对数据进一步进行优化，利用 K-means对步骤（1)中数据进行分类，得到一个基于灰度值的分类结果数据，以此作为骨架 提取的最终输入数据；
[0009]步骤（4)、骨架结构提取：步骤（3)中得到的分类结果数据作为骨架提取的输入数 据，基于步骤（3)中的分类结果数据，利用Ll-中值理论，并定义一个能量泛函，通过求解能 量最小的中心点来代表骨架分支；通过结合正则化，得到一个一维连续的骨架结构；
[0010] 以上的所有操作都基于交互的模式，用户通过调整参数，提取合适的局部特征结 构，而且对于数据中含有噪声和孤立点，甚至对于数据部分缺失都可以提取一个完整的骨 架借结构。
[0011] 以上所有操作都是用Qt和OpenGL来实现的，具有很好的交互性，用户可以通过选 定合适参数对提取不同形状，具有很高的适用性。
[0012] 本发明的原理在于：
[0013] (1)通过微分几何理论来对图像特征进行分析，源自图像微分分析的结构张量可 紧凑的刻画图像局部材质变化的方向以及材质间的过度变化关系，所以可以利用黑塞矩阵 来表达图像特征，通过选取黑塞矩阵特征值满足一定条件的点，可以实现对三维数据的局 部几何结构进行提取。
[0014] (2)为了解决复杂多分支模型情况下对骨架提取的影响，本发明利用K-means聚 类方法，对从（1)中获取的数据进行聚类，该方法可以对点云模型进行有意义的分割。
[0015] (3)为了对点云模型进行骨架提取，本发明依据Ll-中值理论，定义一个能量泛 函，并求解使能量泛函在局部领域中最小值的点作为骨架节点，通过加入正则化，最终得到 一个完整的一维连续的骨架结构。
[0016] 本发明与现有技术相比的优点在于：
[0017] 1、本发明可以直接使用原始CT扫描数据，对于复杂的医学CT扫描数据可以直接 应用而不用进行噪声等的预处理。本发明对于含有噪声、孤立点，甚至是数据部分缺失都能 提取完整的骨架结构。
[0018] 2、良好的交互性，本发明基于Qt来进行交互框架的开发，用户可以通过手动调整 参数来提取图像的不同局部结构，适用度高。
【专利附图】
【附图说明】
[0019] 图1为基于医学图像的器官骨架提取方法流程图；
[0020] 图2为正方体盒子点云模型图；
[0021] 图3为动脉血管的点云模型图；
[0022] 图4k= 3正方体盒子聚类结果图；
[0023] 图5k= 4正方体盒子聚类结果图；
[0024] 图6k= 3动脉血管聚类结果图；
[0025] 图7k= 4动脉血管聚类结果图；
[0026] 图8正方体盒子点云模型的骨架提取结果图；
[0027] 图9动脉血管点云模型的骨架提取结果图。
【具体实施方式】
[0028] 图1给出了基于医学图像的器官骨架提取方法的总体处理流程图，下面结合其他 附图及【具体实施方式】进一步说明本发明。
[0029] 本发明提供一种基于医学图像的器官骨架提取方法，主要步骤介绍如下
[0030] 1、提取三维数据局部特征结构方法
[0031] 本发明采用微分几何理论来对图像特征进行分析，利用图像微分分析的结构张量 可紧凑的刻画图像局部材质变化的方向以及材质间的过渡变化关系，由图像的二阶导数构 成的具有对称性质的黑塞矩阵H能够很好地表示曲面法线的变化情况，进而可以用图像的 二阶导数构成的黑塞矩阵来刻画三维数据的局部几何结构。
[0032] 本发明重点在于利用黑塞矩阵的不同组合来提取合适的局部几何结构，从某种意 义上说，最大特征值所对应的特征向量表征了从一种材质过度到另一种材质时最为迅速的 变化方向，而最小特征值所对应的方向则指明了某种材质在图像内部流动的方向。并且，最 大和最小特征值恰恰蕴含了图像的主曲率信息。因此，从各向异性热传递的角度来说，同时 考虑黑塞矩阵的特征值和特征向量无疑可以克服各向同性滤波器的一些局限性。同时，提 取三维图像的局部结构并对三维图像内部的材质属性变化进行定量刻画也正是我们的目 的所在。
[0033] 通过选择黑塞矩阵H特征值满足不同条件的三维数据来提取不同局部几何结构， 通过将黑塞矩阵的特征值进行不同形式的组合运算，分析了它们与图像局部结构（平面、 块和噪声）之间的关系：
[0034] Rsheet =(I入 3丨 + 1 入 2丨)/1 入ll
[0035] Rblob = (2| 入」_| 入2|_| 入 3|)/2| 入」
[0036] Rtube = (|Aj-|A2|-|A3I)/!A1
[0037]
【权利要求】
1. 一种基于医学图像的器官骨架提取方法，其特征在于包括以下三个步骤： 步骤（1)、读取三维医学图像数据：通过对DICOM格式的医学CT数据读取，得到医学器 官的三维数据，作为整个算法原始输入数据。 步骤（2)、提取三维数据的局部特征结构：通过对得到三维数据固有的结构特征进行 分析来提取局部结构；利用三维数据的黑塞矩阵来刻画局部几何结构，并通过计黑塞矩阵 的特征值进行组合来提取合适的局部几何结构，最大特征值所对应的特征向量表征了从一 种材质过度到另一种材质时最为迅速的变化方向，而最小特征值所对应的方向则指明了某 种材质在图像内部流动的方向；因此可以对三维图像内部的材质属性变化进行定量刻画； 得到的三维离散数据点云模型作为整个骨架提取的输入数据； 步骤（3)、点云模型分割：基于步骤（1)中的数据对数据进一步进行优化，利用K-means 对步骤（1)中数据进行分类，得到一个基于灰度值的分类结果数据，以此作为骨架提取的 最终输入数据； 步骤（4)、骨架结构提取：步骤（3)中得到的分类结果数据作为骨架提取的输入数据， 基于步骤（3)中的分类结果数据，利用Ll-中值理论，并定义一个能量泛函，通过求解能量 最小的中心点来代表骨架分支；通过结合正则化，得到一个一维连续的骨架结构； 以上的所有操作都基于交互的模式，用户通过调整参数，提取合适的局部特征结构， 而且对于数据中含有噪声和孤立点，甚至对于数据部分缺失都能提取一个完整的骨架借结 构。
2. 根据权利要求1所述的基于医学图像的器官骨架提取方法，其特征在于：步骤（2) 中所述的提取三维数据的局部特征结构，该方法使用一组从CT上读取的原始扫描数据，利 用三维数据的局部几何结构，通过求解每个三维数据的黑塞矩阵的特征值，将其提取局部 特征结构转化为选取满足一定特征值条件的三维数据。
3. 根据权利要求1所述的基于医学图像的器官骨架提取方法，其特征在于：步骤（3) 中所述的点云模型分割方法，利用K-means聚类对步骤（2)中得到具有几何结构的数据进 行分类，减少分支间在提取骨架时的互相干扰。
4. 根据权利要求1所述的基于医学图像的器官骨架提取方法，其特征在于：步骤（4) 中所述的骨架结构提取利用Ll-中值理论，找出一组局部邻域最小值点作为骨架分支点、 定义能量泛函来刻画局部邻域最小值，并通过最小化能量泛函求出骨架的分支点；利用 Ll-中值理论求局部邻域最小值点对噪声不敏感，且当点的维数大于二时，求得的局部领域 最小值点唯一；由于单独应用Ll-中值理论，会造成中心点聚积到一起从而产生一个稀疏 的分布，所以需要用正则化理论控制中心点的分布，来实现骨架结构提取。
5. 根据权利要求1所述的基于医学图像的器官骨架提取方法，其特征在于：步骤（2) (3) (4)中所述所有操作都是基于交互的模式；步骤（2)中可以选择点的黑塞矩阵特征值满 足不同条件的三维离散点来实现提取感兴趣的局部结构；步骤（3)中用户可以选择类的个 数，不同的类可以产生不同的效果；步骤（4)用户可以选择采样的数量等利于优化骨架结 构的参数；所有操作和交互都是利用OpenGL和Qt技术进行处理。
【文档编号】G06K9/62GK104318563SQ201410571271
【公开日】2015年1月28日 申请日期:2014年10月22日 优先权日:2014年10月22日
【发明者】李帅, 张东波, 郝爱民, 秦洪 申请人:北京航空航天大学