一种融合显性度量和隐性度量的多模图像配准方法
【专利摘要】为解决基于特征映射的图像配准算法初始误差过大容易引起失配的问题，本发明提出了一种融合显性度量和隐性度量的分层图像配准方法。该方法首先采用互信息作为显性度量，利用粒子群算法获得初始的图像变换参数，对待配准图像进行变换后采用基于特征映射的配准算法进一步优化获取更精确的图像配准参数。为加快图像配准速度并提高配准精度，采用了图像金字塔结构进行分层优化求解。本发明提出的方法适用于多模图像配准，能够提高图像配准的成功率，减小了图像配准误差。
【专利说明】一种融合显性度量和隐性度量的多模图像配准方法
【技术领域】
[0001] 本发明属于图像配准【技术领域】，涉及一种融合显性度量和隐性度量的多模图像配 准方法。
【背景技术】
[0002] 图像配准就是将不同时间、不同传感器或不同视角下获取的两幅或多幅图像进行 匹配、对齐、叠加的过程，是许多图像处理和计算机视觉应用的关键步骤之一，包括图像融 合、图像拼接和视觉伺服等。用于图像配准常用的变换模型包括平移变换、相似变换、仿射 变换和透视变换。
[0003] 由于多模图像成像机理的差异，待配准图像间的灰度相关性不强，采用互相关系 数等作为相似性度量很难获得满意的配准结果，如何提高多模图像配准的成功率，并减小 图像配准误差，是国内外研究学者一直致力于解决的关键问题之一。基于显性度量的图像 配准算法，采用的度量指标包括互信息，相关比，对齐度等像素分布统计特性，由于能够描 述像素间的非线性映射关系，已成功应用于多模医学图像、可见光与红外图像配准领域。基 于隐性度量的多模图像配准算法思想是参考图像中梯度幅值大的点，经过映射后在待配准 图像中对应点的梯度幅值也比较大，所以配准过程中不直接利用图像灰度间的分布关系， 而是通过待配准图像中特征显著点集的特征映射优化得到配准参数。
[0004] 由于多源图像成像机理不同，上述假设有时会失效，此外当图像间变换参数较大 时，可能造成优化算法发散，从而造成图像配准失效。通过融合多种具有不同特性的图像配 准算法，并采用高效的优化算法，能够提高多模图像配准的成功率，减小图像配准的误差。
【发明内容】
[0005] 本发明的目的是为提高多源图像配准算法的鲁棒性，提出了一种融合显性度量和 隐性度量的图像配准方法，采用图像金字塔结构逐级获得更加精确的配准参数，有效解决 了基于隐性度量的图像配准算法初始误差过大容易引起失配的问题。
[0006] 为解决上述技术问题，本发明采用的技术方案包括以下步骤：
[0007] 步骤（1).对原始图像进行一层高斯金字塔分解，获得低分辨率图像；
[0008] 所述的原始图像包括参考图像、待配准图像；
[0009] 步骤（2).在步骤（1)低分辨率图像中，首先采用基于显性度量的方法获得配准 参数初始值P〇,优化算法选用粒子群优化（PS0)算法；然后以p0作为初始值，采用基于隐 性度量的方法获得配准参数值pl，优化算法选用融合粒子群优化（PS0)算法和方向加速 (Powell)算法的优化策略；
[0010] 步骤（3).对原始图像（即更高一层分辨率图像）进行配准，首先利用步骤（2)中 获得的配准参数值pl对待配准图像进行变换，然后再采用步骤（2)中的方法基于参考图像 和经过变换后的待配准图像计算得到更精准的配准参数值P2 ;
[0011] 步骤（4).根据pl和p2计算最终的配准参数值p，并对待配准图像进行仿射变换， 采用双线性插值法得到最终的配准图像。
[0012] 本发明首先采用互信息作为显性度量，利用粒子群算法获得初始的图像变换参 数，然后采用基于隐身度量的配准算法进一步优化获取更精确的图像配准参数。为加快图 像配准速度并提高配准精度，采用了图像金字塔结构进行分层优化求解。本发明与现有技 术相比，其显著优点为：（1)基于隐性度量的图像配准方法采用融合了智能随机优化（PS0) 和共轭梯度直接优化（Powell)的优化策略，在进行参数寻优时，可有效避免陷入局部极大 值的情况。（2)融合了两种具有不同机理的图像配准算法，基于显性度量的算法利用了图像 全部像素的信息，而基于隐性度量的算法则利用了图像中显著特征点的信息，采用分层结 构来逐步获得配准参数。经过先后两次配准会提高配准算法的鲁棒性和精度，解决基于隐 性度量的图像配准算法初始误差过大容易引起失配的问题，提高了多模图像配准的成功率 和精度。
【专利附图】
【附图说明】
[0013] 图1为融合显性度量和隐性度量的多模图像配准算法；
[0014] 图2为仿射变换情况下可见光与可见光图像配准和融合结果，其中（a)为参考图 像，（b)为待配准图像，（c)为最终的配准图像，（d)为融合图像；
[0015] 图3为仿射变换情况下可见光与红外图像配准和融合结果，其中（a)为参考图像， (b)为待配准图像，（c)为最终的配准图像，（d)为融合图像；
[0016] 图4为相似变换情况下可见光与红外图像配准和融合结果，其中（a)为参考图像， (b)为待配准图像，（c)为最终的配准图像，（d)为融合图像；
[0017] 图5为Ku波段和P波段SAR图像配准和融合结果，其中（a)为参考图像，（b)为 待配准图像，（c)为最终的配准图像，（d)为融合图像；
[0018] 图6为P波段和X波段SAR图像配准和融合结果，其中（a)为参考图像，（b)为待 配准图像，（c)为最终的配准图像，（d)为融合图像。
【具体实施方式】
[0019] 以下结合附图对本发明作进一步说明。
[0020] 由图1所示，本发明方法具体实施步骤如下：
[0021] 步骤（1).对原始图像进行一层高斯金字塔分解，获得低分辨率图像；
[0022] 所述的原始图像包括参考图像、待配准图像；
[0023] 步骤（2).在步骤（1)低分辨率图像中，首先采用基于显性度量的方法获得配准 参数初始值P〇,优化算法选用粒子群优化（PS0)算法；然后以p0作为初始值，采用基于隐 性度量的方法获得配准参数值pl，优化算法选用融合粒子群优化（PS0)算法和方向加速 (Powell)算法的优化策略，具体如下：
[0024] 2. 1采用显性度量方法对步骤（1)低分辨率图像获取配准参数初始值p0,优化算 法选用粒子群优化（PS0)算法；
[0025] 所述的显性度量方法中用于图像配准的相似性度量包括互信息、相关比、对齐度 等度量指标；本发明利用互信息通过统计方法描述两幅图像像素间的非线性分布，其公式 为：
[0026] 1(1^ I2) = Hd^+Hd^-Hd!，I2) 式（1);
[0027] 其中，1(1 I2)表示互信息，H(Ii)和H(I2)分别表示图像L和I2的熵，Hd I2) 是两幅图像的联合熵，具体分别表示为：
【权利要求】
1. 一种融合显性度量和隐性度量的多模图像配准方法，其特征在于该方法包括以下步 骤： 步骤（1).对原始图像进行一层高斯金字塔分解，获得低分辨率图像； 所述的原始图像包括参考图像、待配准图像； 步骤（2).在低分辨率图像中，首先采用基于显性度量的方法获得配准参数初始值pO， 优化算法选用粒子群优化PSO算法；然后以pO作为初始值，采用基于隐性度量的方法获得 配准参数值pl，优化算法选用融合粒子群优化PSO算法和方向加速Powell算法的优化策 略； 步骤（3).对原始图像进行配准，首先利用步骤（2)中获得的配准参数值pi对待配准 图像进行变换，然后再采用步骤（2)中的方法基于参考图像和经过变换后的待配准图像计 算得到更精准的配准参数值P2 ; 步骤（4).根据pi和p2计算最终的配准参数值p3,并对待配准图像进行仿射变换，采 用双线性插值法得到最终的配准图像。
2. 根据权利要求1所述的一种融合显性度量和隐性度量的多模图像配准方法，其特征 在于步骤（2)具体包括以下步骤： 2. 1采用显性度量方法对低分辨率图像获取配准参数初始值p0,优化算法选用粒子群 优化PSO算法，具体是： 2. I. 1首选利用显性度量方法中互信息通过统计方法描述两幅图像像素间的非线性分 布，其公式为： Ui1, i2) = Ha1KH(I2)-Hd1, i2) 式⑴； 其中，I (I1, I2)表示互信息，H(I1)和H(I2)分别表示图像I1和I 2的熵，Hd1, I2)是两 幅图像的联合熵，具体分别表示为：
其中i G I1, j G 12, &(/)和分别是图像I1和I2中像素值为i和j的概率， 力是像素值为i和j的联合概率分布； 2. 1.2根据互信息作为判优依据，在仿射变换搜索空间中，利用粒子群优化PSO算法的 搜索策略获得配准参数初始值P〇,具体如下： 选用图像配准变换模型为仿射变换模型，具体变换公式为：
待求解的配准参数写成向量形式为P= [a b c d e f]T ;选取图像中心（x^y。）作为坐 标原点，变换前像素点原始坐标（Xpy1)^1 = Xtl-Xc^y1 = y^y。则变换后对应像素点坐标为 (x 2，y 2)，其中 x 2 - x2+xc，y 2 - y2+yc; 采用粒子群优化PSO算法优化配准参数的具体步骤如下： 步骤（a).初始化： 设定粒子个数为m，迭代次数为T，每个粒子的位置初始值Pi0和粒子移动速度K，均为 在一定区间范围内产生的随机数； 步骤（b).评价每一个粒子，计算粒子的适应度值： 设当前迭代次数为t，给定粒子K，对于基于显性度量的配准算法，首先基于式（5)对图 像I1进行变换，然后利用式（1)?式（4)计算得到互信息，将互信息作为该粒子的适应度； 然后更新第k个粒子移动过程中适应度值取最大的位置，记作Pbesti ;最后计算所有粒子 中适应度值最大的位置，记为Gbest ; 步骤（c).粒子更新： 采用下式（6)和式（7)对每一个粒子的速度F/和位置P/:进行更新； Vk = ' +ci' ran^O ? (Pbest, - P[ ' )+c, ? rand() ? (Gbest -P'k1) 式（6); n=pi-'+n 式⑵； 其中，C1和C2是加速常量，分别调节向全局最好粒子和个体最好粒子方向飞行的最大 步长；rand()是0到1之间的随机数； 步骤（d).检验是否符合结束条件： 判断迭代次数是否达到了预设值，若是则停止迭代，输出最优解；若否则跳转执行步骤 (b)； 2. 2以2. 1得到的配准参数pO作为初始值，采用基于隐性度量的方法获得配准参数值 pl，优化算法选用融合粒子群优化PSO算法和方向加速Powell算法的优化策略，具体步骤 如下： 步骤（a).初始点集检测： 将图像I1按行列分均匀分割成1〇*1〇个子图像，每个大小是
其中（M，N)为图 像I1的宽度和高度；然后在每个子图像中采用式（8)逐像素点计算梯度值；选择每个子图 像中梯度幅值前10%的像素点，构成特征点集S1= ，其中Vi= (Xpyi):
其中，Ix(i，j) =〇.5(I(i，j+l)-I(i，j-l))，Iy(i，j) =0.5(I(i+l，j)-I(i-l，j))， I (i，j)表示图像第i行第j列的像素值； 步骤（b).特征点集映射，计算隐性度量值： 给定配准参数初始值P〇,利用式（5)对特征点集S1里的每个坐标进行变换，得到点集 S2，采用式（9)计算该配准参数pO对应的隐性度量值；
其中，0彡Coi彡1是自适应权重； 步骤（c).配准参数优化求解，采用融合PSO和Powell算法的优化策略求取参数pi : 1) 采用PO作为初始值对粒子进行初始化，利用PSO算法对参数进行优化求解，此时采 用式（9)计算隐性度量值，将其作为粒子的适应度值； 2) 迭代结束后，采用Powell搜索法对全局最优粒子进行进一步寻优；Powell算法中参 数寻优顺序为a - b - d - e - c - f，一维搜索方法选用牛顿迭代法。
3.根据权利要求1所述的一种融合显性度量和隐性度量的多模图像配准方法，其特征 在于步骤（4)具体是： 由步骤⑵得到的第一层配准参数Pl = La1 Id1 C1 Cl1 ei fJT，以及第二层配准参数p2 =[a2 b2 c2 d2 e2 f2]T，由公式（10)得到最终配准参数 p3 = [a b c d e f]T : B - &2&1+"^2屯 b = aA+t^ei c = B2C^b2Tx+C2 式（10); d = Ci2B^e2Ci1e = Ci2C^e2G1 f = d^i+esfi+fs 利用上述得到的最终配准结果p3对待配准图像进行仿射变换映射，采用双线性插值 法计算变换后图像的像素值，从而得到最终的配准后图像。
【文档编号】G06T7/00GK104331891SQ201410613265
【公开日】2015年2月4日 申请日期:2014年11月4日 优先权日:2014年11月4日
【发明者】谷雨, 徐英, 彭冬亮, 苟书鑫 申请人:杭州电子科技大学