一种基于高斯曲率流的脊椎特征点自动识别方法技术

一种基于高斯曲率流的脊椎特征点自动识别方法，首先，对人体脊椎的CT图像进行三维重建，得到所需标注特征点的椎体三维模型，对步骤一得到的椎体模型进行网格划分；其次，手动选取模型上某点，计算在手动选取点周围极小半径为R的球形空间内每个顶点的平均曲率值；最后，选取手动选取点周围高斯曲率最大的n个点，对高斯曲率最大的n个点，分别与选取点作内积，夹角最小的即为所求的点，能够提高特征点标注的精确度，从而提高数据测量精确度，具有精确度高、易于实现的特点。

An automatic recognition method of spine feature points based on Gauss curvature flow

An automatic identification method, the feature points of Gauss spine curvature flow based on the first, CT images of the human spine were reconstructed to obtain the required labeling of vertebral 3D model feature points, mesh model step by the vertebral body; secondly, manually select a point on the model, calculate the average curvature for each vertex sphere of the R value in the manual selection of points around the minimum radius; finally, select the manual selection of n points around the point of maximum Gauss curvature, n points of maximum curvature and Gauss, respectively selected as the inner product, minimum angle is required, can improve the accuracy of feature point and annotation. To improve the accuracy of measurement data, has the characteristics of high accuracy, easy to implement.

【技术实现步骤摘要】
一种基于高斯曲率流的脊椎特征点自动识别方法
本专利技术属于医学图像处理
，特别涉及一种基于高斯曲率流的脊椎特征点自动识别方法。
技术介绍
脊椎位于背部正中，是人体的支撑骨骼，对人体的内脏有着重要的保护作用。由于近年来人们的生活和工作方式的改变，人们长期保持同一种姿势，致使脊椎的负担加重，脊椎疾病的发病率越来越高，但是因为脊椎复杂的结构以及脊椎手术本身的特点导致脊椎手术的难度很高，而且伴随着很大的风险。脊椎模型的精确定位和配准问题仍然是还没有解决的诊疗难题。在脊椎模型的定位和配准中，主要任务是特征点的定位。为了标注特征点，一般是采用纯手动的标注的方法。这种方法标注出特征点不准确，导致数据的精度不高，模型匹配不准确。
技术实现思路
为了克服上述现有技术的不足，本专利技术的目的在于提出一种基于高斯曲率流的脊椎特征点自动识别方法，能够有效提高特征点标注的准确性，从而提高数据测量精确度，具有精确度高、易于实现的特点。为了实现上述目的，本专利技术采用的技术方案是：一种基于高斯曲率流的脊椎特征点自动识别方法，步骤如下：步骤一，对人体脊椎的CT图像进行三维重建，得到所需标注特征点的椎体三维模型；步骤二，对步骤一得到的椎体模型进行网格划分；步骤三，手动选取模型上某点，计算在手动选取点周围极小半径为R的球形空间内每个顶点的高斯曲率值的绝对值；步骤四，选取手动选取点周围高斯曲率最大的n个点；步骤五，对高斯曲率最大的n个点，分别与选取点作内积，夹角最小的即为所求的点。所述步骤四和步骤五中的n≤5。所述步骤三的具体过程如下：步骤3-1，按照如下步骤计算出椎体模型上每个顶点的泰森多边形区域面积AM：a.遍历脊椎三维图像表面的三角面片，获得其三个顶点的坐标v0、v1和v2，并得出边向量e0＝v1-v0，e1＝v2-v1，e2＝v0-v2；b.根据向量的点积公式，将边向量(e1,e2)、(e2,e0)和(e0,e1)分别带入，得到三角面片各个顶角的度数，再根据cell的拓扑关系，计算脊椎图像三维表面上每个点邻域三角形对应的角度和其中#f为三角形面片的个数；c.计算由这三个顶点组成的三角形的面积A，同理根据cell的拓扑关系可得到每个点附近三角形的面积AM；步骤3-2，利用以下公式遍历计算范围内每个点的高斯曲率值：其中，#f表示以点xi为顶点的三角形的数目，AM为包含点xi的泰森多边形区域的面积。所述步骤五的具体过程如下：步骤5-1，利用以下公式计算高斯曲率值最大的n个点与手动选取点的两个向量内积：＜pi-O,p0-O＞＝xix0+yiy0+ziz0其中，O为坐标系原点，p0为手动选取点，pi为平均曲率值最大的n个点之一；步骤5-2，利用以下公式计算两向量的夹角θ：与手动选取点夹角最小的点即为经过修正后的特征点。与现有技术相比，本专利技术具有以下有益效果：经过对手动选取点的高斯曲率修正后获得的特征点的精度比较高，且容易实现；算法只对手动选取的特征点周围r范围内n个顶点进行运算，效率较高；使用高斯曲率值自适应计算得到的特征点，经过对手动选取点的修正，较之前通过计算高斯曲率值手动选取特征点，更加趋近于真实特征点；较之通过平均曲率选取特征点，因为某个特征点的高斯曲率值远远大于其平均曲率值，固本改进算法能很好地避免由于人为因素产生的偏差，大大提高了特征点标注的准确性，对脊柱高发疾病的诊断和治疗具有一定的参考价值。附图说明图1为本专利技术方法的流程图。图2(a)为三维重建得到人体脊椎模型正面图，2(b)为三维重建得到人体脊椎模型背面图。图3(a)为模型区域生长分割原始图像，3(b)为模型区域生长分割结果图。图4(a)为腰段L1-L3节模型建立效果图，4(b)为腰段L1节腰段L1节。图5(a)为腰椎L1节模型俯视图，5(b)为腰椎L1节模型主视图。图6为将椎体模型导入到mimics中的显示。图7(a)为腰椎L1-L3节体网格模型实体单元表示，7(b)为腰椎L1-L3节体网格模型节点表示。图8为脊椎腰骶段的右侧面观特征点标记图。图9为脊椎腰骶段的上面观特征点标记图。图10(a)为采用高斯曲率手动选择特征点结果的仰视图，10(b)为采用高斯曲率手动选择特征点结果的左视图，(c)为采用高斯曲率手动选择特征点结果的右视图。图11(a)为采用文献“三维脊椎模型特征点标注的研究与实现”选择特征点结果的仰视图，11(b)为采用文献“三维脊椎模型特征点标注的研究与实现”选择特征点结果的左视图，11(c)为采用文献“三维脊椎模型特征点标注的研究与实现”选择特征点结果的右视图。图12(a)为采用本方法自适应选择特征点结果的仰视图，12(b)为采用本方法自适应选择特征点结果的左视图，12(c)为采用本方法自适应选择特征点结果的右视图。图13为泰森多边形法示意图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。参见图1，一种基于高斯曲率流的脊椎特征点自动识别方法，包括以下步骤：步骤一，对人体脊椎的CT图像进行三维重建，得到所需标注特征点的椎体三维模型。参见图2、图3、图4，对脊椎CT数据进行增强和阈值提取，分割出目标椎体的范围。对椎体的分割可以使用mimics软件实现，如图5、图6所示。在模型构建过程中，可以采用图像增强，图像阈值分割和区域生长算法，去噪点等步骤保证模型的准确。步骤二，对步骤一得到的椎体模型进行网格划分。设置对椎体进行网格划分的参数，利用mimics软件生成模型的三角面片格式，如图7所示。利用该方式生成的三角面片能够自动优化局部复杂结构的网格划分，在保证整体网格质量的同时，降低网格数量，简化操作流程。步骤三，在针对算法开发的特征点标注软件中，使用鼠标选取特征点，脊椎采样特征点的定义如图8、图9所示。计算在手动选取点周围极小半径为R的球形空间内每个顶点的高斯曲率值kG。步骤3-1，按照如下步骤计算出椎体模型上每个顶点的泰森多边形区域面积AM:遍历脊椎三维图像表面的三角面片，获得其三个顶点的坐标v0、v1和v2，并得出边向量e0＝v1-v0，e1＝v2-v1，e2＝v0-v2。根据向量的点积公式，将边向量(e1,e2)、(e2,e0)和(e0,e1)分别带入，可以得到三角面片各个顶角的度数θj，再根据cell的拓扑关系，计算脊椎图像三维表面上每个点邻域三角形对应的角度和计算由这三个顶点组成的三角形的面积A，同理根据cell的拓扑关系可得到每个点附近三角形的面积AM，计算公式如下：步骤3-2，利用以下公式遍历计算范围内每个点的高斯曲率值：其中，#f表示以点xi为顶点的三角形的数目，AM为包含点xi的泰森多边形区域的面积，θj所表示的角度如图13所示，表示点xi邻域三角形对应的角度和，近似为点xi邻域法线映射到单位球上的面积，再除以相应区域的面积就刻画了该点附近曲面的弯曲程度；步骤四，在手动选取点周围半径为R的球形空间内，选取高斯曲率最大的5个点作为候选特征点。以手动选取点为中心，首先确定一个半径R，对半径R内的球形空间内所有顶点的高斯曲率值，进行降序排序，取得高斯曲率值最大的5个点作为候选特征点。步骤五，对5个候选特征点，分别与手动选取点作内积，夹角最小的为所求的点。利用以下公式分别计算高斯曲率值最大的5个点与手动选取点的内积。＜pi-O,p0-O＞＝x本文档来自技高网...

【技术保护点】
一种基于高斯曲率流的脊椎特征点自动识别方法，其特征在于，步骤如下：步骤一，对人体脊椎的CT图像进行三维重建，得到所需标注特征点的椎体三维模型；步骤二，对步骤一得到的椎体模型进行网格划分；步骤三，手动选取模型上某点，计算在手动选取点周围极小半径为R的球形空间内每个顶点的高斯曲率值的绝对值；步骤四，选取手动选取点周围高斯曲率最大的n个点；步骤五，对高斯曲率最大的n个点，分别与选取点作内积，夹角最小的即为所求的点。

【技术特征摘要】
1.一种基于高斯曲率流的脊椎特征点自动识别方法，其特征在于，步骤如下：步骤一，对人体脊椎的CT图像进行三维重建，得到所需标注特征点的椎体三维模型；步骤二，对步骤一得到的椎体模型进行网格划分；步骤三，手动选取模型上某点，计算在手动选取点周围极小半径为R的球形空间内每个顶点的高斯曲率值的绝对值；步骤四，选取手动选取点周围高斯曲率最大的n个点；步骤五，对高斯曲率最大的n个点，分别与选取点作内积，夹角最小的即为所求的点。2.根据权利要求1所述的一种基于高斯曲率流的脊椎特征点自动识别方法，其特征在于，所述步骤四和步骤五中的n≤5。3.根据权利要求1所述的一种基于高斯曲率流的脊椎特征点自动识别方法，其特征在于，所述步骤三的具体过程如下：步骤3-1，按照如下步骤计算出椎体模型上每个顶点的泰森多边形区域面积AM:a.遍历脊椎三维图像表面的三角面片，获得其三个顶点的坐标v0、v1和v2，并得出边向量e0＝v1-v0，e1＝v2-v1，e2＝v0-v2；b.根据向量的点积公式，将边向量(e1,e2)、(e2,e0)和(e0,e1)分别带入，可以得到三角面片各个顶角的度数，再根据cel...

【专利技术属性】
技术研发人员：惠宇，武君胜，鱼滨，贺伟，杨柳，杜静，李航，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61