一种CT图像中脊髓的提取方法发明专利.docxVIP

一种CT图像中脊髓的提取方法 技术领域 本发明涉及图像处理领域，尤其涉及一种CT图像中脊髓的提取方法。 背景技术 在放疗的准备阶段，医生通过病人的CT图像计算放疗过程中所使用的放射线路经 和剂量，同时在CT图像中勾勒(分割)出一些重要的器官，以使其在放疗中避开放射线的照 射，或者降低其所受到的放射线照射剂量。脊髓是神经系统的核心组成部分，也是放疗中必 须保护的重要区域。在制定针对头部(如脑、颈)、胸部(如乳房、肺、食道)和腹部(如肝脏、脾 脏、肾脏、胰脏)等癌症的放疗计划时，为了保护脊髓，医生都需要将脊髓从预先得到的CT图 像中分割出来。 目前，对于三维CT图像中椎管的自动分割可以分为两类。第一类方法需要人工确 定椎管的位置，比如人工选取椎管内部的一个或多个种子点，然后从种子点出发，进一步自 动分割出椎管部分，是一种半自动方法(请参考[L.G.,et al.:Computer Analysis of Images and Patterns,456-463.(2005)])。 第二类方法可以自动检测出椎管的大致位置，然后对椎管进行提取。例如，通过简 单的圆霍夫变换寻找轴切面上的圆形特征来定位椎管，然后通过模糊区域增长的方法在三 维空间中分割出椎管；使用三维广义霍夫变换直接定位出任意一节脊椎骨的空间位置，并 确定其所包围的椎管中心点作为种子点，然后通过一种自适应的算法，将小管状结构逐段 连接到种子点上(请参考[T.Klinder,et al.:Medical Image Analysis,vol.13,471. (2009)])；以及使用一套基于解剖结构图谱的决策支持系统，用于定位二维轴切面上椎管 的位置，然后使用区域增长或活动轮廓模型的方法来分割三维椎管。但是，这类方法过于复 杂，计算量比较大。 发明内容 本发明解决的问题是提供一种CT图像中脊髓的提取方法，不仅能自动、快速地提 取CT图像中的脊髓，还具有更高的稳定性和灵活性。 为了解决上述问题，本发明提供了一种CT图像中脊髓的提取方法，包括如下步骤： (1)输入所述三维CT图像，通过二维简化的椎管模型，得到所述三维CT图像中所述 椎管的中轴线位置的三维广义霍夫变换空间； (2)根据所述三维广义霍夫变换空间，确定所述椎管中轴线在所述三维CT图像中 的位置； (3)根据所述椎管中轴线的位置，通过区域增长方法，得到所述椎管的分割掩模， 即得到所述三维CT图像中的脊髓。 上述所述一种CT图像中脊髓的提取方法，其中，得到所述椎管中轴线位置的三维 广义霍夫变换空间的过程如下： 1)通过边缘检测方法，得到第k帧轴切面图像中的边缘点位置和所述边缘点位 置处的正规化梯度向量其中，k为所述三维CT图像中第k帧二维轴切面图像，j＝ 1,...,Nk为所述第k帧轴切面图像中的第j个边缘点，Nk为所述得到的第k帧轴切面图像中的 所有边缘点数； 2)根据所述边缘点位置和所述正规化梯度向量得到所述边缘点映射到广义 霍夫变换空间中的位置和并进行投票，得到所述广义霍夫空间中的点分布函 数和 3)将所述点分布函数和进行合并，得到所述第k帧轴切面的广义霍夫 变换空间Hk(c)； 4)将所有轴切面的广义霍夫变换空间在轴向进行级联，得到H(c,k)； 5)将所述H(c,k)进行平滑，得到所述椎管的中轴线的三维广义霍夫变换空间H(c, k)。 上述所述一种CT图像中脊髓的提取方法，其中，所述点分布函数的公式为： 其中，为投票函数，I[·]为指标函数；所述?r1为所述二维简化模型中椎孔的半径。 上述所述一种CT图像中脊髓的提取方法，其中，所述点分布函数的公式为： 其中，为投票函数，I[·]为指标函数，为所述 正规化梯度向量沿着人体冠状轴的分量；所述r2为所述二维简化模型 中椎体的半径；为所述二维简化模型中椎孔中轴线与轴切面的交点以及椎体中轴线与轴 切面的交点之间的相对位移。 上述所述一种CT图像中脊髓的提取方法，其中，得到所述Hk(c)的公式为或其中，w为所述所占的 权重。 上述所述一种CT图像中脊髓的提取方法，其中，所述H(c,k)的公式为： 其中，k0和ke分别为所述输入的三维CT图像中感兴趣区域的起始帧数和终止帧数。 上述所述一种CT图像中脊髓的提取方法，其中，确定所述椎管中轴线在所述三维 CT图像中的位置的过程为： 1)确定所述H(c,k)的极大值的像素点位置(cmax,kmax)； 2)根据快速行进算法和速度函数F(c,k)，得到所述三维CT图像中每个像素点到所 述像素点(cmax,kmax)的距离函数D(c,k；cmax,kmax)； 3)分别确定在所述第k0帧和第ke帧轴切面图像中距离函数D(c,k；cmax,kmax)的最小 值位置c0和ce； 4)根据最大梯度下降法，分别得到所述c0和ce到所述像素点(c