[发明专利]颅外轮廓图像通用配准方法

摘要

本发明公开了一种颅外轮廓图像通用配准方法，包括分别提取参考图像和浮动图像的外轮廓点云；分别计算参考图像和浮动图像的质心；分别将两模态的质心平移至与原点重合，得到新的参考模态和浮动模态；分别计算得出参考模态和浮动模态的第一根特征轴向量和第二根特征轴向量；根据求得的特征轴向量构造旋转算子，并完成浮动图像的旋转；平移和旋转后完成配准。本发明几何意义简单直观，经过实验和客观参考对比分析，配准精度高，并具有通用性，便于医生对病人病情进行准确地分析。

主权项

1.一种颅外轮廓图像通用配准方法，其特征是：包括下列步骤：(1)分别提取参考图像和浮动图像的外轮廓点云；(2)分别计算参考图像和浮动图像的质心；(3)分别将两模态的质心平移至与原点重合，得到新的参考模态和浮动模态；(4)分别计算得出参考模态和浮动模态的第一根特征轴向量和第二根特征轴向量；(5)根据求得的特征轴向量构造旋转算子，并完成浮动图像的旋转；(6)平移和旋转后完成配准；步骤(4)的具体方法是：(1)计算外轮廓点云中各点到质心的欧几里得距离，并找出距离值最大时所对应的点；在得到了外轮廓数据之后，在轮廓点云中找到一个点，该点满足到质心o的欧几里得距离最远；即计算轮廓上任意点到相应质心o的欧几里得距离，取欧几里得距离最远时所对应的轮廓上的点；为了便于计算，此处取欧几里得距离的平方，如式(1)：所述的式(1)中o＝xoe1+yoe2+zoe3，X_i＝x_i﹒e₁+y_i﹒e₂+z_i﹒e₃，i＝1,2,L,n；n为刚体外轮廓点云总数，x_o、y_o、z_o为常数，x_i、y_i、z_i为常数，e₁、e₂、e₃是几何代数子空间的一组基；(2)以过所述的质心和所述的最大欧几里得距离对应的外轮廓上的点的矢量作为相应模态的第一根轴向量；刚体的第一根特征轴向量记作参考模态和浮动模态的第一根特征轴向量分别是(3)以第一根轴向量作为直线方向向量，计算对应模态上所有点到此直线的距离，以映射向量的距离的平方作为距离测度，寻找到该直线距离最大的点；以作为直线方向向量，计算刚体上所有点到此直线的距离，利用映射向量的距离的平方来表示为距离测度，寻找最大值，如式(2)：所述的式(2)中o＝xoe1+yoe2+zoe3，X_i＝x_i﹒e₁+y_i﹒e₂+z_i﹒e₃，i＝1,2,L,n；，n为刚体外轮廓点云总数，x_o、y_o、z_o为常数，x_i、y_i、z_i为常数，e₁、e₂、e₃是几何代数子空间的一组基；(4)取所述的到第一根轴向量距离最远的点到所述的直线方向向量的映射向量作为第二根特征轴向量；相应模态的第二根特征轴向量记作参考模态和浮动模态的第二根特征轴向量分别是(5)单位化处理之后的第一根特征轴和第二根特征轴就是所要建立的几何相对特征不变量；参考模态和浮动模态的第一根特征轴向量分别是参考模态和浮动模态的第二根特征轴向量分别是根据计算得到的几何相对特征不变量构造旋转算子；步骤(3)中所述的质心分别为o_r和o_f，所述的新的参考模态和浮动模态分别是其中，n_r是新的参考模态外轮廓点的个数，n_f是新的浮动模态外轮廓点的个数，i＝1,2,3,...,n_r，和分别是新的参考模态和浮动模态外轮廓点云。