基于空间距离与数据分布相似性度量的医学图像配准方法技术

本发明专利技术公开了一种基于空间距离与数据分布相似性度量的医学图像配准方法，属于医学图像配准领域。本发明专利技术提出了结合mp‑dissimilarity和lp‑norm两种相似性度量的匹配策略，显著提高了在医学多模图像配准中的准确率，对医疗工作提供了帮助。本发明专利技术只是针对于描述子的相似性度量，也可以将它运用在其它基于局部特征的图像配准工作中。

【技术实现步骤摘要】
基于空间距离与数据分布相似性度量的医学图像配准方法
本专利技术涉及医学图像配准领域，特别涉及一种基于空间距离与数据分布相似性度量的医学图像配准方法。
技术介绍
随着计算机技术以及生物医学的快速发展，医学影像学为临床诊断提供了不同模态的图像。比如计算机断层成像(CT)、核磁共振图像(MRI)、数字减影血管造影成像(DSA)、单光子辐射计算机断层成像(SPECT)以及正电子辐射断层成像(PET)。不同模态的医学图像对人体相关器官的表示信息是不同的，比如CT图像空间分辨率比较高，得到的图像更加清晰，有利于定位病灶；MRI对于软组织成像清晰，对于确定病灶的范围有优势；PET以及SPET成像虽然空间分辨率不如前面的两种成像，但是却能够提供器官的功能和代谢信息，所以对于医疗工作者迫切希望得到不同成像信息集成的结果。然而不同模态的医学图像不仅分辨率不同而且还存在空间上的形变，因此，图像配准对于处理这些不同模态的图像是必要的。图像配准的目的就是将不同模式的图像通过空间变换映射到同一坐标系中，在医学图像中，可以使相应器官的影像在空间中的位置一致。通过配准后得到的图像就可以同时反映形态和功能的信息，为医疗工作提供更可靠地依据。因为图像配准的显著优点，它在病灶定位、PACS系统、放射治疗计划、指导神经手术以及检查治疗效果许多医疗工作中有着广泛的应用。基于局部特征的图像配准方法中，最后都需要度量两幅图像中的特征点相似性。最常用的都是lp-norm可是这种方法只是简单地计算特征点之间的空间距离。另外一种度量方法，mp-dissimilarity着重于用数据之间的分布情况表示特征点相似性。但是现有的方法对医学图像的配准准确率低。
技术实现思路
为了弥补现有技术的不足，解决现有技术中医学图像配准准确率低的问题，本专利技术提供了一种基于空间距离与数据分布相似性度量的医学图像配准方法。本专利技术验证了mp-dissimilarity用于图像配准中的可行性，并且验证了两种相似性度量方法具有互补信息，从而提出了结合这两种方法的匹配策略。实验结果证明，本专利技术提出的匹配策略能够显著地提高特征点匹配的准确度。术语解释1、SIFT，尺度不变特征转化，一种经典的特征提取方法。因为具有尺度不变性、旋转不变性以及对噪声、光线的变化都不敏感，获得了良好的效果。2、mp-dissimilarity，一种将数据分布作为主要因素的相似性度量方法。这种方法的基本思想是两对点具有相同距离时，周围有更多距离的一对相似性较差。3、lp-norm是一种典型的用空间距离度量相似性的方法。一般有三种类型，由参数p决定，其中p是一个变参数。当p＝1时，就是曼哈顿距离，当p＝2时，就是欧几里得距离，当p→∞时，就是切比雪夫距离。在图像处理以及模式识别中被广泛使用的相似性判断方法。其中，mp-dissimilarity与lp-norm中的P不是同一参数(是在不同公式中代表不同参数)。在mp-dissimilarity的相关公式中，P代表控制维数影响的参数；在lp-norm的相关公式中，P表示具体的某一范式，P＝2时为欧几里得距离。本专利技术的技术方案为：一种基于空间距离与数据分布相似性度量的医学图像配准方法，包括步骤：1)输入待配准的两幅多模医学图像，其中一幅为参考图像，另一幅为目标图像；2)采用SIFT算法提取待配准图像中的特征点，给每个特征点赋予主方向，得到特征点的位置、尺度以及方向信息，使得特征点具有旋转不变性，生成128维的特征点描述子；3)目标图像与参考图像特征点匹配a)定义两个向量x,y的lp-norm为：其中，xi以及yi分别表示向量x,y中的第i个元素；p的值为2；abs表示取绝对值；d表示每个向量的维数；b)定义两个向量x,y的mp-dissimilarity为：其中，R是R(x,y)∈Rd的缩写，表示由向量x，y组成的区域；z是R中随机选择的一个点；Pi(z∈R|o(z))表示z这个点落在R中第i维的概率；d表示x,y总的维数；p作为一个参数，控制维数的影响；c)使用mp-dissimilarity度量描述子的相似性假设两幅图像中的描述子为：参考图像中共有Nq个描述子，目标图像中共有Nt个描述子，并且每一个描述子都是由d维的数组构成，当用mp-dissimilarity计算由Dq和Dt组成的描述子时，D为d×(Nq+Nt)维数组；其中，D表示Dq和Dt的并集；D中的任意两个描述子的mp-dissimilarity：e)使用mp-dissimilarity匹配描述子表示参考图像中的一个描述子，和分别表示在目标图像中Mp-dissimilarity距离与最近以及次近的描述子，匹配条件为：(MPi,Nq+j1)表示与的mp-dissimilarity距离，(MPi,Nq+j2)表示与的mp-dissimilarity距离；threshold是设定好的阈值；f)结合lp-norm与mp-dissimilarity寻找匹配的特征点。作为优选方案，步骤2)的具体为：Ⅰ)建立尺度空间，并且检测尺度空间极值：搜索所有尺度上的图像位置，通过高斯微分函数来识别潜在的对于尺度和旋转不变的兴趣点；Ⅱ)关键点定位，找到兴趣点后，在每个候选的位置上，通过拟合精细的模型确定位置和尺度；Ⅲ)方向确定，基于图像局部的梯度方向，分配给每个关键点位置一个或多个方向；所有后面的对图像数据的操作都相对于关键点的方向、尺度和位置进行变换，从而提供对于这些变换的不变性；Ⅳ)关键点描述，在每个关键点周围的邻域内，在选定的尺度上测量图像局部的梯度；这些梯度被变换成一种表示，这种表示允许局部形状的变形和光照变化；这样所有的关键点都被表示为128维的数组，也被称为描述子。作为一种优选方案，步骤f)具体为：找到参考图像中用欧式距离匹配的描述子，再求这些描述子与目标图像中所有描述子Dt之间的mp-dissimilarity；D1＝Dmq∪Dt；其中，Dmq表示与Med对应的在参考图像中的描述子；Med表示使用欧式距离得到的匹配对。作为另一种优选方案，步骤f)具体为：先找到目标图像中用欧式距离匹配的描述子，再求这些描述子与参考图像中所有描述子Dq之间的mp-dissimilarity；D2＝Dmt∪Dq；其中，Dmt表示与Med对应的在目标图像中的描述子，Med表示使用欧式距离得到的匹配对。作为再一种优选方案，步骤f)具体为：用欧式距离找到两幅图像中匹配的描述子，然后计算这些描述子之间的mp-dissimilarity；D3＝Dmq∪Dmt；其中，Dmq表示与Med对应的在参考图像中的描述子；Dmt表示与Med对应的在目标图像中的描述子；Med表示使用欧式距离得到的匹配对。作为又一种优选方案，步骤f)具体为：分别用欧氏距离以及mp-dissimilarity寻找两幅图像中匹配的描述子，然后取它们的交集；Mf＝Med∩Mmp；其中，Med表示使用欧式距离得到的匹配对；Mmp表示使用mp-dissimilarity得到的匹配对。本专利技术的有益效果为：1、本专利技术提出了结合mp-dissimilarity和lp-norm两种相似性度量的匹配策略，显著提高了在医学多模图像配准中的准确率，对医疗工作提供了帮助。2、本专利技术只是针对于描述子的本文档来自技高网...

【技术保护点】
一种基于空间距离与数据分布相似性度量的医学图像配准方法，其特征在于，包括步骤：1)输入待配准的两幅多模医学图像，其中一幅为参考图像，另一幅为目标图像；2)采用SIFT算法提取待配准图像中的特征点，给每个特征点赋予主方向，得到特征点的位置、尺度以及方向信息，使得特征点具有旋转不变性，生成128维的特征点描述子；3)目标图像与参考图像特征点匹配a)定义两个向量x,y的

【技术特征摘要】
1.一种基于空间距离与数据分布相似性度量的医学图像配准方法，其特征在于，包括步骤：1)输入待配准的两幅多模医学图像，其中一幅为参考图像，另一幅为目标图像；2)采用SIFT算法提取待配准图像中的特征点，给每个特征点赋予主方向，得到特征点的位置、尺度以及方向信息，使得特征点具有旋转不变性，生成128维的特征点描述子；3)目标图像与参考图像特征点匹配a)定义两个向量x,y的为：其中，xi以及yi分别表示向量x,y中的第i个元素；p的值为2；abs表示取绝对值；d表示每个向量的维数；b)定义两个向量x,y的mp-dissimilarity为：其中，R是R(x,y)∈Rd的缩写，表示由向量x，y组成的区域；z是R中随机选择的一个点；Pi(z∈R|o(z))表示z这个点落在R中第i维的概率；d表示x,y总的维数；p作为一个参数，控制维数的影响；c)使用mp-dissimilarity度量描述子的相似性假设两幅图像中的描述子为：参考图像中共有Nq个描述子，目标图像中共有Nt个描述子，并且每一个描述子都是由d维的数组构成，当用mp-dissimilarity计算由Dq和Dt组成的描述子时，D为d×(Nq+Nt)维数组；其中，D表示Dq和Dt的并集；D中的任意两个描述子的mp-dissimilarity：e)使用mp-dissimilarity匹配描述子表示参考图像中的一个描述子，和分别表示在目标图像中Mp-dissimilarity距离与最近以及次近的描述子，匹配条件为：(MPi,Nq+j1)表示与的mp-dissimilarity距离，(MPi,Nq+j2)表示与的mp-dissimilari...

【专利技术属性】
技术研发人员：吕国华，颜晗，任晓强，董祥军，
申请(专利权)人：齐鲁工业大学，
类型：发明
国别省市：山东,37