1.一种基于颅外轮廓特征几何不变量多信息3D医学图像配准方法,其特征是:包括下列步骤:(1)分别提取参考图像和浮动图像的外轮廓点云;
(2)分别计算参考图像和浮动图像的质心;
(3)分别将两模态的质心平移至与原点重合,得到新的参考模态和浮动模态;
(4)分别计算得出参考模态和浮动模态的第一根特征轴向量和第二根特征轴向量;
(5)根据求得的特征轴向量构造旋转算子,并完成浮动图像的旋转;
(6)平移和旋转后完成配准;
步骤(4)的具体方法是:
(1)计算外轮廓点云中各点到质心的欧几里得距离,并找出距离值最大时所对应的点;
在得到了外轮廓数据之后,在轮廓点云中找到一个点,该点满足到质心o的欧几里得距离最远;即计算轮廓上任意点到相应质心o的欧几里得距离,取欧几里得距离最远时所对应的轮廓上的点;为了便于计算,此处取欧几里得距离的平方,如式(1):所述的式(1)中o=xoe1+yoe2+zoe3,
Xi=xi﹒e1+yi﹒e2+zi﹒e3,i=1,2,…,n;n为刚体外轮廓点云总数,xo、yo、zo为常数,xi、yi、zi为常数,e1、e2、e3是几何代数 子空间的一组基;
(2)以过所述的质心和所述的最大欧几里得距离对应的外轮廓上的点的矢量作为相应模态的第一根轴向量;刚体的第一根特征轴向量记作 参考模态和浮动模态的第一根特征轴向量分别是(3)以第一根轴向量作为直线方向向量,计算对应模态上所有点到此直线的距离,以映射向量的距离的平方作为距离测度,寻找到该直线距离最大的点;以 作为直线方向向量,计算刚体上所有点到此直线的距离,利用映射向量的距离的平方来表示为距离测度,寻找最大值,如式(2):所述的式(2)中o=xoe1+yoe2+zoe3,
Xi=xi﹒e1+yi﹒e2+zi﹒e3,i=1,2,…,n;n为刚体外轮廓点云总数,xo、yo、zo为常数,xi、yi、zi为常数,e1、e2、e3是几何代数 子空间的一组基;
(4)取所述的直线距离最大的点到所述的直线方向向量的映射向量作为第二根特征轴向量;相应模态的第二根特征轴向量记作 参考模态和浮动模态的第二根特征轴向量分别是(5)单位化处理之后的第一根特征轴和第二根特征轴就是所要建立的几何相对特征不变量;参考模态和浮动模态的第一根特征轴向量分别是 参考模态和浮动模态的第二根特征轴向量分别是 根据计算得到的几何相对特征不变量构造旋转算子。
2.根据权利要求1所述的基于颅外轮廓特征几何不变量多信息3D医学图像配准方法,其特征是:步骤(1)的具体方法:根据参考图像与浮动图像的3D模态数据构成的离散点云,提取参考图像和浮动图像的外轮廓点云,分别记作 和 其中,nr是参考模态外轮廓点的个数,nf是浮动模态外轮廓点的个数,i=1,2,3,...,n, 和 分别是参考模态和浮动模态外轮廓点云。