一种基于刚性变换的多模态图像融合方法技术

本发明专利技术公开了一种基于刚性变换的多模态图像融合方法，包括以下的步骤：S1：通过变换矩阵构建参考图像R与浮动图像F之间的刚性变换模型；S2：构建参考图像R与浮动图像F之间的相似度目标函数f；S3：利用模拟退火算法对相似度目标函数f进行全局寻优，得到相似度目标函数f全局最优解对应的变换矩阵；S4：将步骤S3得到的全局最优解对应的变换矩阵中的平移量及旋转量作为局部寻优的搜索起点，利用鲍威尔算法对相似度目标函数f进行局部寻优，得到相似度目标函数f局部最优解对应的变换矩阵，将局部最优解对应的变换矩阵代入步骤S1的刚性变换模型中，得到准确配准后的浮动图像F；S5：将参考图像R与准确配准后的浮动图像F进行加权融合。

【技术实现步骤摘要】

本专利技术涉及图像融合方法，特别是涉及一种基于刚性变换的多模态图像融合方法。
技术介绍
在临床应用中，经常需要综合多种设备采集到的影像以做出诊断，不同设备采集到的影像可以提供不同的信息，因此利用图像配准、融合的方法将多种设备上对同一部位采集到的影像进行融合有很重要的临床价值。现有技术中的图像融合方法往往存在如下的缺陷：1)很多医生利用手动去配准两幅图像，利用PS等工具做融合，手动意味着人为的因素存在，得到结果不可靠；2)融合的关键在于配准，现有的配准技术使用基于特征点的配准方法，这种针对照相机照片的配准方法不适合做医学配准，医学图像中的特征点不可靠；3)现有的融合方法配准精度不高，在医学图像中需要达到亚像素级别的精度。
技术实现思路
专利技术目的：本专利技术的目的是提供一种能够解决现有技术中存在的缺陷的基于刚性变换的多模态图像融合方法。技术方案：为达到此目的，本专利技术采用以下技术方案：本专利技术所述的基于刚性变换的多模态图像融合方法，包括以下的步骤：S1：通过变换矩阵构建参考图像R与浮动图像F之间的刚性变换模型；S2：构建参考图像R与浮动图像F之间的相似度目标函数f；S3：利用模拟退火算法对相似度目标函数f进行全局寻优，得到相似度目标函数f全局最优解对应的变换矩阵；S4：将步骤S3得到的全局最优解对应的变换矩阵中的平移量及旋转量作为局部寻优的搜索起点，利用鲍威尔算法对相似度目标函数f进行局部寻优，得到相似度目标函数f局部最优解对应的变换矩阵，将局部最优解对应的变换矩阵代入步骤S1的刚性变换模型中，得到准确配准后的浮动图像F；S5：将参考图像R与准确配准后的浮动图像F进行加权融合。进一步，所述步骤S1中的刚性变换模型如式(1)所示： x f y f 1 = c o s θ - s i n θ x s i n θ cos θ y 0 0 1 · x r y r 1 - - - ( 1 ) ]]>式(1)中，xr为浮动图像F中任意一点的横坐标，yr为浮动图像F中该点的纵坐标，xf为浮动图像F中该点变换到参考图像R中得到的相应点的横坐标，yf为浮动图像F中该点变换到参考图像R中得到的相应点的纵坐标，x为浮动图像F的水平平移量，y为浮动图像F的垂直平移量，θ为浮动图像F的逆时针旋转角度，为变换矩阵。进一步，所述步骤S2中的相似度目标函数f如式(2)所示： f ( w ) = N M I ( R , F ) = H ( R ) + H ( F ) H ( R , F ) - - - ( 2 ) ]]>式(2)中，w＝(x,y,θ)，H(R)为参考图像R的熵，H(F)为浮动图像F的熵，H(R,F)为参考图像R和浮动图像F的联合熵。进一步，所述步骤S3包括以下的步骤：S3.1：将(x,y,θ)＝(0,0,0)作为初始解w0，计算相应的相似度目标函数f(w0)；S3.2：利用扰动函数对w0进行扰动，得到新解w，并计算相应的相似度目标函数f(w)；S3.3：判断f(w)和f(w0)的大小：如果f(w)>f(w0)，则将w赋给w0，然后继续进行步骤S3.4；如果f(w)<f(w0)，则根据Metropolis准则决定是否接受w，如果接受，则将w赋给w0，然后继续进行步骤S3.4，如果不接受，则进行步骤S3.5；S3.4：判断是否达到最大迭代次数：如果已达到最大迭代次数，则跳至步骤S3.5；如果未达到最大迭代次数，则返回步骤S3.2；S3.5：判断此时的相似度目标函数f是否满足要求：如果满足要求，则结束；如果不满足要求，则重新设定最大迭代次数，然后返回步骤S3.2。进一步，所述步骤S4包括以下的步骤：S4.1：将所述步骤S3中得到的全局最优解对应的变换矩阵中本文档来自技高网...

【技术保护点】
一种基于刚性变换的多模态图像融合方法，其特征在于：包括以下的步骤：S1：通过变换矩阵构建参考图像R与浮动图像F之间的刚性变换模型；S2：构建参考图像R与浮动图像F之间的相似度目标函数f；S3：利用模拟退火算法对相似度目标函数f进行全局寻优，得到相似度目标函数f全局最优解对应的变换矩阵；S4：将步骤S3得到的全局最优解对应的变换矩阵中的平移量及旋转量作为局部寻优的搜索起点，利用鲍威尔算法对相似度目标函数f进行局部寻优，得到相似度目标函数f局部最优解对应的变换矩阵，将局部最优解对应的变换矩阵代入步骤S1的刚性变换模型中，得到准确配准后的浮动图像F；S5：将参考图像R与准确配准后的浮动图像F进行加权融合。

【技术特征摘要】
1.一种基于刚性变换的多模态图像融合方法，其特征在于：包括以下的步骤：S1：通过变换矩阵构建参考图像R与浮动图像F之间的刚性变换模型；S2：构建参考图像R与浮动图像F之间的相似度目标函数f；S3：利用模拟退火算法对相似度目标函数f进行全局寻优，得到相似度目标函数f全局最优解对应的变换矩阵；S4：将步骤S3得到的全局最优解对应的变换矩阵中的平移量及旋转量作为局部寻优的搜索起点，利用鲍威尔算法对相似度目标函数f进行局部寻优，得到相似度目标函数f局部最优解对应的变换矩阵，将局部最优解对应的变换矩阵代入步骤S1的刚性变换模型中，得到准确配准后的浮动图像F；S5：将参考图像R与准确配准后的浮动图像F进行加权融合。2.根据权利要求1所述的基于刚性变换的多模态图像融合方法，其特征在于：所述步骤S1中的刚性变换模型如式(1)所示： x f y f 1 = c o s θ - s i n θ x sin θ cos θ y 0 0 1 · x r y r 1 - - - ( 1 ) ]]>式(1)中，xr为浮动图像F中任意一点的横坐标，yr为浮动图像F中该点的纵坐标，xf为浮动图像F中该点变换到参考图像R中得到的相应点的横坐标，yf为浮动图像F中该点变换到参考图像R中得到的相应点的纵坐标，x为浮动图像F的水平平移量，y为浮动图像F的垂直平移量，θ为浮动图像F的逆时针旋转角度，为变换矩阵。3.根据权利要求1所述的基于刚性变换的多模...

【专利技术属性】
技术研发人员：刘金虎，王东辉，加冯刚，
申请(专利权)人：南京普爱医疗设备股份有限公司，
类型：发明
国别省市：江苏;32