CN107274389A - 基于ct三维序列图像的股骨和髋臼解剖参数获取方法 - Google Patents

Abstract

本发明公开一种基于CT三维序列图像的股骨和髋臼解剖参数获取方法，包括步骤：获取具有股骨和髋臼的CT三维序列图像；股骨与髋臼的分割：对CT三维序列图像依次进行粗分割、分离以及精确分割，获得股骨精确分割图像和髋臼精确分割图像；股骨解剖参数计算：根据股骨精确分割图像，构造股骨表面的三角网格表面模型，计算股骨直径参数、髓腔内径参数以及股骨头的球面拟合参数；髋臼解剖参数计算：选择髋臼精确分割图像构造髋臼表面三角网格表面模型的表面若干点，计算髋臼的前后径。本发明直接对CT三维序列图像进行处理，快速、自动地实现股骨和髋臼精确分割及解剖参数获取，辅助人工骨头假体的个性化设计与建模。

Description

基于CT三维序列图像的股骨和髋臼解剖参数获取方法

技术领域

本发明涉及股骨和髋臼的参数计算技术领域，更具体地说，本发明涉及 一种基于CT三维序列图像的股骨和髋臼解剖参数获取方法。

背景技术

在治疗股骨头坏死的半髋关节置换临床手术中，不同患者的髋臼骨骼形 态差异很大，如果使用标准件式的人工股骨头假体进行置换，术后效果不理 想，因此，人工股骨头假体的个性化设计就显得十分必要。

现多采用如下的技术方案解决这一个问题；1.人工进行分割，操作人员 手动进行分割提取；2.以大量髋关节图像数据作为训练集，训练得出髋臼的 统计学形状或者利用广义hough变换得到其平均形状，然后对于具体患者的 髋臼形状进行匹配后调整；3.以一副人工分割好的数据，作为先验形状，然 后对于具体患者的髋臼形状进行匹配后调整；4.首先得到髋关节的二值图像， 可得到股骨头和髋关节粘连在一起的初始轮廓，然后利用图像的灰度和邻域 信息将粘连部分分开，最后再利用梯度信息进一步矫正。

但是上述现有技术存在如下缺陷：首先，每个患者的髋关节形状差异很 大，得到的先验形状很难和具体患者的匹配，并且，医生测量髋臼和股骨解 剖参数一般是在二维图像上手动进行的，如X射线透视图或CT的二维断层 图，这就造成了测量存在很大的误差，因此由先验形状得到的分割结果不准 确，无法解决髋臼和股骨头个性化分割与问题；2.整个过程需要复杂的人工 交互，对操作人员的经验要求较高，不同操作人员得到的分割结果可能差别 很大；3.如果使用智能算法，需要积累大量的数据集样本，才能保证结果的 相对正确，实现难度很大，一般的非医学院附属医院很难积累足够的样本。

在一些已公开的专利文献中，也提出了使用图像处理的方式来解决上述 技术问题。但是，这些已公开的技术方案，基本思路还是通过在CT断层图 像中逐层分割股骨，以实现股骨分割和建模的目的，并未提供一个有效的髋 臼、股骨同时精确分割的方法，也不能对髋臼测的解剖参数进行自动测量。

发明内容

针对上述技术中存在的不足之处，本发明提供一种基于CT三维序列图 像的股骨和髋臼解剖参数获取方法，直接对CT三维序列图像进行建模，快 速、自动地实现股骨和髋臼的精确分割及解剖参数获取，辅助人工骨头假体 的个性化设计。

为了实现根据本发明的这些目的和其它优点，本发明通过以下技术方案 实现：

本发明提供一种基于CT三维序列图像的股骨和髋臼解剖参数获取方法， 其包括以下步骤：

获取具有股骨和髋臼的CT三维序列图像；

股骨与髋臼的分割：对所述CT三维序列图像依次进行粗分割、分离以 及精确分割，获得股骨精确分割图像和髋臼精确分割图像；

股骨解剖参数计算：根据所述股骨精确分割图像，构造股骨表面的三角 网格表面模型，计算股骨直径参数、髓腔内径参数以及股骨头的球面拟合参 数；

髋臼解剖参数计算：根据所述髋臼精确分割图像，构造髋臼表面的三角 网格表面模型，选择髋臼的三角网格表面模型表面的若干点，计算髋臼的前 后径；

其中，所述股骨直径参数和所述髓腔内径参数分别包括最大外接圆直径 D_外、最小内接圆直径D_内以及最佳拟合圆直径D_佳；所述股骨头的球面拟合 参数包括股骨头的拟合球半径和颈干角。

优选的是，所述CT三维序列图像的粗分割，包括以下步骤：

采用阈值分割方法分割CT三维序列图像中的骨骼，得到骨骼的初分割 图像；

对骨骼初分割图像依次进行半径不同的形态闭操作，获得股骨和髋臼的 二次分割图像；

对股骨和髋臼的二次分割图像分别进行基于二值迭代空洞填充算法的空 洞填充，获得两个骨骼粗分割图像。

优选的是，所述股骨与髋臼的分离，包括以下步骤：

对两个所述骨骼粗分割图像分别依次进行阈值增长分割、半径不同的形 态学开操作，获得去髋臼的股骨图像和去除股骨的髋臼图像。

优选的是，所述股骨与髋臼的精确分割，包括以下步骤：

分别以所述股骨图像和所述髋臼图像为初始曲面，初始化水平集函数LS _股骨；

基于区域信息和边缘信息构造水平集函数，采用窄带法对水平集函数LS _股骨进行水平集演化，分别获得股骨精确分割图像和髋臼精确分割图像。

优选的是，计算股骨直径参数，包括以下步骤：

使用等值面提取算法构造股骨的表面模型，并进一步进行三角网格化， 获得股骨的三角网格表面模型；

通过指定中心或自动生成来提取股骨的中心线；

分别以股骨中心线上等间距分布的若干点为圆心，计算股骨三角网格表 面模型的断面所对应的最大外接圆直径D_外、最小内接圆直径D_内和最佳拟合 圆直径D_佳。

优选的是，若指定中心来提取股骨的中心线，即，指定股骨中心线提取 的起点、终点以及中心线经过的中间端点，则执行以下步骤：

以股骨头表面的一点为起点，以股骨最下端的一点为终点，采用最大通 量算法计算股骨三角网格表面模型中股骨的第一中心线；

以大转子表面的一点为起点，以股骨最下端的一点为终点，采用最大通 量算法计算股骨三角网格表面模型中股骨的第二中心线；

分别遍历第一中心线和第二中心线中的所有离散点，分别构造第一中心 线和第二中心线的离散点集合CL_S1、CL_S2，则，CL_S1＝{p₁，…，p_n}，CL_S2＝{q₁，…， q_m}，p_n和q_m分为为点的三维坐标；

比较任意p_n和q_m：若|p_n-q_m|＜t且|p_n-q_m|＜|p_n-q_m’|(m’＝1,..,m-1,m+1,…,M)，则使用(p_n+q_m)/2分别取代pn和q_m进行坐标更新；坐标更新后的p_n、q_m分别 记作p’_n和q’_m；并将修改后p’_n和q’_m设置为不可修改，获得将第一中心线和 第二中心线合并成的一条二分叉的第三中心线；

使用端点检测算法检测端点坐标；

或，若自动生成股骨的中心线，则执行以下步骤：

基于股骨三角网格表面模型，采用无端点约束的中心线提取算法，提取 股骨的第四中心线；

手动选择不必要的中心线分支，并将其删除。

优选的是，计算股骨头的球面拟合参数，包括以下步骤：

离散化第一中心线，获得分布于第一中心线上的离散点集；

构造第一中心线上任一离散点的法平面，法平面与股骨三角网格表面模 型产生一个闭合截面曲线；

离散化闭合截面曲线，获得分布于闭合截面曲线上的离散点集；

对闭合截面曲线上的离散点集进行圆拟合，计算拟合得到圆周的半径；

以第一中心线的起点为起点，以第一中心线的终端为终点，搜索圆周半 径的最大阶跃点；

构造法平面，最大阶跃点位于法平面上，法平面的法向量为第一中心线 在最大阶跃点的主切向量；法平面法向量的方向指向第一中心线的端点侧；

以法平面为切平面，去除法向量侧的三角面片模型，获得股骨头的表面 模型；

对骨头的表面模型依次进行面片化简、平滑和游离面去除，得到股骨头 表面模型的后处理曲面，后处理曲面的顶点记作P_i＝(x_i,y_i,z_i)，由P_i拟合球面 的公式为：Err＝Σ_i[(x_i-x₀)²+(y_i-y₀)²+(z_i-z₀)²-r₂]；通过最小化Err得到拟合球面的 参数，r为拟合球面的半径，x₀、y₀和z₀为拟合球面的球心；

使用第一中心线上的离散点拟合直线，其拟合得到的直线记作第一直线， 计算其方向向量；使用第二中心线上的离散点拟合直线，其拟合得到的直线 记作第二直线，计算其方向向量；计算第一直线和第二直线的内积，即可得 到第一直线的方向向量与第二直线的方向向量之间的夹角，即为颈干角。

优选的是，髋臼解剖参数计算，包括以下步骤：

采用等值面提取算法构造髋臼的表面模型，并进一步进行三角网格化， 获得髋臼的三角网格表面模型；

选择髋臼三角网格表面模型表面的若干个点，计算髋臼的前后径。

优选的是，在获取CT三维序列图像与所述股骨与髋臼分割之间，还包 括步骤：

对具有股骨和髋臼的CT三维序列图像进行预处理，所述预处理至少包 括采用三维高斯平滑算法对所述CT三维序列图像进行平滑滤波处理。

优选的是，对所述CT三维序列图像进行粗分割与分离之间，还包括步 骤：采用多尺度Hessian矩阵对股骨和髋臼间隙部分进行增强。

本发明至少包括以下有益效果：

本发明提供的基于CT三维序列图像的股骨和髋臼解剖参数获取方法， 通过对获取的具有股骨和髋臼的CT三维序列图像，依次进行股骨与髋臼的 粗分割、分离、精确分割的处理方法以精确分割股骨和髋臼，在分割的同时 自动测量计算股骨直径参数、髓腔内径参数、股骨头的球面拟合参数以及髋 臼前后径等股骨与髋臼的参数；整个分割和测量过程直接在CT三维序列图 像上进行，无需手动在单个CT的二维断层图像上进行分割，快速、自动地 实现股骨、髋臼以及髓腔的精确分割以及相关解剖参数获取，辅助人工骨头 假体的个性化设计，智能化水平较高。

本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将 通过对本发明的研究和实践而为本领域的技术人员所理解。

附图说明

图1为本发明基于CT三维序列图像的股骨和髋臼解剖参数获取方法流 程图；

图2为本发明所述的CT三维序列图像粗分割的方法流程图；

图3为本发明所述的CT三维序列图像中股骨与髋臼分离的方法流程图；

图4为本发明所述的CT三维序列图像中股骨与髋臼精确分割的方法流 程图。

具体实施方式

下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照 说明书文字能够据以实施。

应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不 排除一个或多个其它元件或其组合的存在或添加。

如图1所示，本发明提供一种基于CT三维序列图像的股骨和髋臼解剖 参数获取方法，其包括以下步骤：

S10，获取具有股骨和髋臼的CT三维序列图像；

S20，股骨与髋臼的分割：对CT三维序列图像依次进行粗分割、分离以 及精确分割，获得股骨精确分割图像和髋臼精确分割图像；

S30，股骨解剖参数计算：根据股骨精确分割图像，构造股骨表面的三角 网格表面模型，计算股骨直径参数、髓腔内径参数以及股骨头的球面拟合参 数；

S40，髋臼解剖参数计算：根据髋臼精确分割图像，构造髋臼表面的三角 网格表面模型，选择髋臼的三角网格表面模型表面的若干点，计算髋臼的前 后径。

上述实施方式中，步骤S10中，具有股骨和髋臼的CT三维序列图像为 CT三维序列图像。

步骤20中，如图2所示，对CT三维序列图像的粗分割包括以下步骤：

S21，采用阈值分割方法分割CT三维序列图像中的骨骼，得到骨骼的初 分割图像。阈值分割方法的低阈值，记作T_l、高阈值为CT三维序列图像的 CT值最大值，记作T_max，则阈值分割后获得的骨骼初分割图像记作I₁。

S22，对骨骼初分割图像依次进行半径不同的形态闭操作，获得股骨和髋 臼的二次分割图像。具体为，使用半径为R_e1的模板对I₁进行膨胀操作，在 膨胀操作后的图像上使用半径为R_r1的模板进行腐蚀操作，并且R_e＜R_r，形态 学闭操作后的图像记作IL₁；使用半径为R_e2的模板对I₁进行膨胀操作，在膨 胀操作后的图像上使用半径为R_r2的模板进行腐蚀操作，并且R_e2＝R_r2，形态 学闭操作后的图像记作IL₂。

S23，对股骨和髋臼的二次分割图像分别进行基于二值迭代空洞填充算法 的空洞填充，获得两个骨骼粗分割图像。具体为，对IL₁和IL₂使用模板半径 为R₁的模板进行基于二值迭代空洞填充算法的空洞填充，以填补步骤S22中 二次分割后图像中的空洞，获得两个粗分割图像，分别记作IL_1-1和IL_2-1；以 IL_1-1为模板，得到分割后的骨骼粗分割图像，记作I_1-1；以IL_2-1为模板，得到 分割后的骨骼粗分割图像，记作I_2-1。

步骤20中，如图3所示，对CT三维序列图像中股骨与髋臼的分离，包 括以下步骤：S24，对两个骨骼粗分割图像分别依次进行阈值增长分割、半径 不同的形态学开操作，获得去髋臼的股骨图像和去除股骨的髋臼图像。股骨 与髋臼的分离，采用区域增长方法，具体为，在骨骼粗分割图像I_1-1和骨骼粗 分割图像I_2-1中，分别选择一个或若干个位于股骨侧的种子点；设置低阈值为 T₂且T₂≥T₁、阈值增长的高阈值T_l＝T_max来对两幅粗分割图像分别进行区域增长 分割；对区域增长后的分割结果使用形态学开操作，即先进行腐蚀操作再进行膨胀操作，腐蚀和膨胀的模板半径都为R₂，得到股骨初步分割后去除髋臼 的结果，记作I_股骨-1，以及髋臼初步分割后去除股骨的结果，记作I_髋臼-1。

步骤20中，如图4所示，对CT三维序列图像中股骨与髋臼的精确分割， 包括以下步骤：S25，分别以股骨图像和髋臼图像为初始曲面，初始化水平集 函数LS_股骨；S26，基于区域信息和边缘信息构造水平集函数，采用窄带法对 水平集函数LS_股骨进行水平集演化，分别获得股骨精确分割图像和髋臼精确分 割图像。具体地，分别以I_股骨-1和I_髋臼-1为初始曲面，初始化水平集函数LS_股骨；采用区域信息和边缘信息构造水平集函数，在水平集演化过程中使用窄带 法，分别得到股骨和髋臼的精确分割图像，分别记作I_股骨-2和I_髋臼-2。

步骤30中，股骨直径参数和髓腔内径参数分别包括最大外接圆直径D_外、 最小内接圆直径D_内以及最佳拟合圆直径D_佳；股骨头的球面拟合参数包括股 骨头的拟合球半径和颈干角。

具体地，计算股骨直径参数，包括以下步骤：

S31，使用等值面提取算法构造股骨的表面模型，并进一步进行三角网格 化，获得股骨的三角网格表面模型，记作S_股骨。

S32，通过指定中心或自动生成来提取股骨的中心线，其中，

1.若指定中心来提取股骨的中心线，即指定中心线提取的起点、终点以 及中心线经过的中间端点，则执行以下步骤；

以股骨头表面的一点P_s1为起点，以股骨最下端的一点P_t为终点，采用最 大通量算法计算股骨三角网格表面模型中股骨的第一中心线，记作CL₁；

以大转子表面的一点P_s2为起点，以股骨最下端的一点P_t为终点，采用最 大通量算法计算股骨三角网格表面模型中股骨的第二中心线，记作CL₂；

分别遍历第一中心线CL₁和第二中心线CL₂中的所有离散点，分别构造 第一中心线CL₁和第二中心线CL₂的离散点集合CL_S1、CL_S2，则，CL_S1＝{p₁，…， p_n}，CL_S2＝{q₁，…，q_m}，p_n和q_m分为为点的三维坐标；

比较任意p_n和q_m：若|p_n-q_m|＜t且|p_n-q_m|＜|p_n-q_m’|(m’＝1,..,m-1,m+1,…,M)，则使用(p_n+q_m)/2分别取代pn和q_m进行坐标更新；坐标更新后的p_n、q_m分别 记作p’_n和q’_m；并将修改后p’_n和q’_m设置为不可修改，获得将第一中心线 CL₁和第二中心线CL₂合并成的一条二分叉的第三中心线CL₃；

使用端点检测算法检测端点坐标，则第三中心线CL_3的起点记为P_d、终点 记为P_t；

2.若自动生成股骨的中心线，则执行以下步骤：

基于股骨三角网格表面模型，采用无端点约束的中心线提取算法，提取 股骨的第四中心线CL₄；

手动选择不必要的中心线分支，并将其删除。

S33，分别以股骨中心线上等间距分布的若干点为圆心，计算股骨三角网 格表面模型的断面所对应的最大外接圆直径D_外、最小内接圆直径D_内和最佳 拟合圆直径D_佳。

步骤S30中，髓腔内径参数的计算，包括步骤：以S_股骨为ROI区域，将 S_股骨外侧标记为背景；采用阈值分割，分割得到髓腔；采用等值面提取算法 的方法得到髓腔表面并进行三角面片化处理，髓腔的三角面片表面记作S_髓腔； 使用自动的中心线提取算法计算S_髓腔的中心线；以中心线上等间距分布若干 点为圆心，计算S_股骨断面所对应的最大外接圆直径D_外，最小内接圆直径D_内 和最佳拟合圆直径D_佳及其最小值，并记录该点坐标。

步骤S30中，计算股骨头的球面拟合参数，包括以下步骤：

1)、离散化第一中心线CL₁，获得分布于第一中心线CL₁上的离散点集， 记作PCL₁；

2)、构造第一中心线CL₁上任一离散点PCL₁的法平面PN，法平面PN 与S_股骨三角网格表面模型产生一个闭合截面曲线，记作C；

3)、离散化闭合截面曲线C，获得分布于闭合截面曲线C上的离散点集， 记作P_c，则P_c＝{p_c1,…,p_cN’}；

对闭合截面曲线上的离散点集P_c进行圆拟合，计算拟合得到圆周的半径 RPC；

4)、以第一中心线CL₁的起点P_s1为起点，以第一中心线P_d的终端为终 点，搜索圆周半径的最大阶跃点P_J；

5)、构造法平面PN，最大阶跃点P_J位于法平面上，法平面PN的法向量 为第一中心线CL₁在最大阶跃点P_J的主切向量；法平面PN法向量的方向指 向第一中心线CL₁的端点侧P_d；

6)、以法平面PN为切平面，去除法向量侧的三角面片模型，获得股骨 头的表面模型，记作S_股骨；

7)、对骨头的表面模型依次进行面片化简、平滑和游离面去除，得到股 骨头表面模型S_股骨的后处理曲面，后处理曲面的顶点记作P_i＝(x_i,y_i,z_i)，由P_i拟合球面的公式为：Err＝Σ_i[(x_i-x₀)²+(y_i-y₀)²+(z_i-z₀)²-r₂]；通过最小化Err得到拟 合球面的参数，r为拟合球面的半径，x₀、y₀和z₀为拟合球面的球心；

8)、使用第一中心线CL₁上的离散点拟合直线，其拟合得到的直线记作 第一直线CL_1L，计算其方向向量CL_1D；使用第二中心线CL₂上的离散点拟合 直线，其拟合得到的直线记作第二直线CL_2L，计算其方向向量CL_2D；计算第 一直线CL_1L和第二直线CL_2L的内积，即可得到第一直线CL_2L的方向向量 CL_1D与第二直线CL_2L的方向向量CL_2D之间的夹角，即为颈干角。

步骤S30中，髋臼解剖参数计算，包括以下步骤：

采用等值面提取算法构造髋臼的表面模型，并进一步进行三角网格化， 获得髋臼的三角网格表面模型；

选择髋臼三角网格表面模型表面的若干个点，计算髋臼的前后径。

在步骤S10获取CT三维序列图像与步骤S20股骨与髋臼分割之间，还 包括步骤：

S50，对具有股骨和髋臼的CT三维序列图像进行预处理，预处理至少包 括采用三维高斯平滑算法对CT三维序列图像进行平滑滤波处理。

步骤S20中，对CT三维序列图像进行粗分割与分离之间，作为优选， 还包括步骤：采用多尺度Hessian矩阵对股骨和髋臼间隙部分进行增强。具体 地，包括以下步骤：

1)在原始CT三维序列图像中选择包含股骨和髋臼整个间隙的ROI区域， 对该ROI区域中体素计算其Hessian矩阵的特征值，分别记作λ₁,λ₂和λ₃；

2)由于股骨和髋臼之间间隙是一个近似于具有一定厚度的曲面，其局部 结构可以看做是盘状结构，因此其间隙的局部Hessian矩阵的特征具有以下特 征：λ₁<<λ₂≤λ₃；

3)构造盘状目标函数：

其中，Ra＝λ₂/λ₃，α、β和c是敏感性参 数；一般情况下，α和β取值为0.5，c根据图像灰度范围调整；当体素位于 盘状目标中时，v输出一个较大的响应。

4)由于股骨和髋臼之间间隙厚度存在一定差异，使用多尺度Hessian滤 波器对图像进行滤波：

其中，δ是空间尺度因子，当δ与血管实际宽度接近或匹配时，血管相似 性函数将输出一个较大的值。在本发明中，δ的取值范围介于[d_min/4,d_max/4]， 即δ_min＝d_min/4，δ_max＝d_max/4。d_min和d_max表示图像中血管的最小宽度和最大宽 度。在δ的取值范围内，将其等分为N个尺度，逐级增大δ，计算v的输出 相应，并计算其最大响应：

5)设置盘状目标检测阈值t，若v₀≥t，则认为该像素位于盘状目标内， 将该像素置为背景值，否则认为该像素位于盘状目标之外，保持原灰度值不 变。

本发明提供的基于CT三维序列图像的股骨和髋臼解剖参数获取方法， 通过对获取的具有股骨和髋臼的CT三维序列图像，依次进行股骨与髋臼的 粗分割、分离、精确分割的处理方法以精确分割股骨和髋臼，在分割的同时 自动测量计算股骨直径参数、髓腔内径参数、股骨头的球面拟合参数以及髋 臼前后径等股骨与髋臼的参数；整个分割和测量过程直接在CT三维序列图 像上进行，无需手动在单个CT的二维断层图像上进行分割，快速、自动地 实现股骨、髋臼以及髓腔的精确分割以及相关参数获取，辅助人工骨头假体 的个性化设计，智能化水平较高。

尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方 式中所列运用。它完全可以被适用于各种适合本发明的领域。对于熟悉本领 域的人员而言可容易地实现另外的修改。因此在不背离权利要求及等同范围 所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (10)

1.一种基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，其包括以下步骤：

获取具有股骨和髋臼的CT三维序列图像；

股骨与髋臼的分割：对所述CT三维序列图像依次进行粗分割、分离以及精确分割，获得股骨精确分割图像和髋臼精确分割图像；

股骨解剖参数计算：根据所述股骨精确分割图像，构造股骨表面的三角网格表面模型，计算股骨直径参数、髓腔内径参数以及股骨头的球面拟合参数；

髋臼解剖参数计算：根据所述髋臼精确分割图像，构造髋臼表面的三角网格表面模型，选择髋臼的三角网格表面模型表面的若干点，计算髋臼的前后径；

其中，所述股骨直径参数和所述髓腔内径参数分别包括最大外接圆直径D_外、最小内接圆直径D_内以及最佳拟合圆直径D_佳；所述股骨头的球面拟合参数包括股骨头的拟合球半径和颈干角。

2.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，所述CT三维序列图像的粗分割，包括以下步骤：

采用阈值分割方法分割CT三维序列图像中的骨骼，得到骨骼的初分割图像；

对骨骼初分割图像依次进行半径不同的形态闭操作，获得股骨和髋臼的二次分割图像；

对股骨和髋臼的二次分割图像分别进行基于二值迭代空洞填充算法的空洞填充，获得两个骨骼粗分割图像。

3.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，所述股骨与髋臼的分离，包括以下步骤：

对两个所述骨骼粗分割图像分别依次进行阈值增长分割、半径不同的形态学开操作，获得去髋臼的股骨图像和去除股骨的髋臼图像。

4.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，所述股骨与髋臼的精确分割，包括以下步骤：

分别以所述股骨图像和所述髋臼图像为初始曲面，初始化水平集函数LS_股骨；

基于区域信息和边缘信息构造水平集函数，采用窄带法对水平集函数LS_股骨进行水平集演化，分别获得股骨精确分割图像和髋臼精确分割图像。

5.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，计算股骨直径参数，包括以下步骤：

使用等值面提取算法构造股骨的表面模型，并进一步进行三角网格化，获得股骨的三角网格表面模型；

通过指定中心或自动生成来提取股骨的中心线；

分别以股骨中心线上等间距分布的若干点为圆心，计算股骨三角网格表面模型的断面所对应的最大外接圆直径D_外、最小内接圆直径D_内和最佳拟合圆直径D_佳。

6.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，若指定中心来提取股骨的中心线，即，指定股骨中心线提取的起点、终点以及中心线经过的中间端点，则执行以下步骤：

以股骨头表面的一点为起点，以股骨最下端的一点为终点，采用最大通量算法计算股骨三角网格表面模型中股骨的第一中心线；

以大转子表面的一点为起点，以股骨最下端的一点为终点，采用最大通量算法计算股骨三角网格表面模型中股骨的第二中心线；

分别遍历第一中心线和第二中心线中的所有离散点，分别构造第一中心线和第二中心线的离散点集合CL_S1、CL_S2，则，CL_S1＝{p₁，…，p_n}，CL_S2＝{q₁，…，q_m}，p_n和q_m分为为点的三维坐标；

比较任意p_n和q_m：若|p_n-q_m|＜t且|p_n-q_m|＜|p_n-q_m’|(m’＝1,..,m-1,m+1,…,M)，则使用(p_n+q_m)/2分别取代pn和q_m进行坐标更新；坐标更新后的p_n、q_m分别记作p’_n和q’_m；并将修改后p’_n和q’_m设置为不可修改，获得将第一中心线和第二中心线合并成的一条二分叉的第三中心线；

使用端点检测算法检测端点坐标；

或，若自动生成股骨的中心线，则执行以下步骤：

基于股骨三角网格表面模型，采用无端点约束的中心线提取算法，提取股骨的第四中心线；

手动选择不必要的中心线分支，并将其删除。

7.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，计算股骨头的球面拟合参数，包括以下步骤：

离散化第一中心线，获得分布于第一中心线上的离散点集；

构造第一中心线上任一离散点的法平面，法平面与股骨三角网格表面模型产生一个闭合截面曲线；

离散化闭合截面曲线，获得分布于闭合截面曲线上的离散点集；

对闭合截面曲线上的离散点集进行圆拟合，计算拟合得到圆周的半径；

以第一中心线的起点为起点，以第一中心线的终端为终点，搜索圆周半径的最大阶跃点；

构造法平面，最大阶跃点位于法平面上，法平面的法向量为第一中心线在最大阶跃点的主切向量；法平面法向量的方向指向第一中心线的端点侧；

以法平面为切平面，去除法向量侧的三角面片模型，获得股骨头的表面模型；

对骨头的表面模型依次进行面片化简、平滑和游离面去除，得到股骨头表面模型的后处理曲面，后处理曲面的顶点记作P_i＝(xi,y_i,z_i)，由P_i拟合球面的公式为：Err＝Σ_i[(x_i-x₀)²+(y_i-y₀)²+(z_i-z₀)²-r₂]；通过最小化Err得到拟合球面的参数，r为拟合球面的半径，x₀、y₀和z₀为拟合球面的球心；

使用第一中心线上的离散点拟合直线，其拟合得到的直线记作第一直线，计算其方向向量；使用第二中心线上的离散点拟合直线，其拟合得到的直线记作第二直线，计算其方向向量；计算第一直线和第二直线的内积，即可得到第一直线的方向向量与第二直线的方向向量之间的夹角，即为颈干角。

8.如权利要求1所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，髋臼解剖参数计算，包括以下步骤：

采用等值面提取算法构造髋臼的表面模型，并进一步进行三角网格化，获得髋臼的三角网格表面模型；

选择髋臼三角网格表面模型表面的若干个点，计算髋臼的前后径。

9.如权利要求1-8中任一项所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，在获取CT三维序列图像与所述股骨与髋臼分割之间，还包括步骤：

对具有股骨和髋臼的CT三维序列图像进行预处理，所述预处理至少包括采用三维高斯平滑算法对所述CT三维序列图像进行平滑滤波处理。

10.如权利要求1-8中任一项所述的基于CT三维序列图像的股骨和髋臼解剖参数获取方法，其特征在于，对所述CT三维序列图像进行粗分割与分离之间，还包括步骤：采用多尺度Hessian矩阵对股骨和髋臼间隙部分进行增强。