一种无人机影像滤波分频拼接方法
【专利摘要】本发明属于无人机遥感影像处理领域，涉及一种无人机影像滤波分频拼接方法，主要用来解决无人机影像拼接易出现鬼影、拼接缝问题。本发明首先利用高斯低通滤波将待拼接无人机影像对分解成低频影像对和高频影像对，然后对低频影像对采用加权平滑融合方法进行拼接处理，利用改进的动态规划算法搜索最优拼接线完成高频影像的拼接，最后将拼接好的低频影像对和高频影像对进行线性合成得到最终的拼接影像。本发明不仅可以较好的解决无人机影像拼接过程中出现的鬼影问题，最大化的避免拼接缝的出现，而且对亮度差异明显的无人机影像对也能取得良好的拼接效果。
【专利说明】
-种无人机影像滤波分频拼接方法
技术领域
[0001] 本发明属于无人机遥感影像处理领域，特别是设及一种无人机影像滤波分频拼接 方法。
【背景技术】
[0002] 无人机遥感系统具有灵活、实时、成本低的特点，在小区域和飞行困难地区快速获 取高分辨率影像方面具有显著优势，已经成为提高测绘成果现势性的有力手段，被广泛的 应用于军事、民用测绘领域中。但受飞行高度、相机视角的限制，单幅无人机影像的影像范 围很难完全包含感兴趣的区域，使其应用受到了很大的限制。为了扩大视场范围，获得更多 目标区域的信息，就需要对多幅具有一定重叠度的无人机影像进行拼接。
[0003] 由于无人机质量较轻，易受高空风力影响，飞行姿态不稳定，得到的影像重叠不规 贝1J、曝光不均匀，使得拼接影像易出现鬼影、拼接过渡不连续，存在明显的拼接缝等现象。针 对运些问题，国内外的学者进行了大量研究，解决办法主要有加权平滑融合方法、基于最优 拼接线方法。
[0004] 加权平滑融合方法，主要是根据重叠区域内的像素点在两幅影像上的距离关系， 采用加权平滑处理来消除拼接缝，简单易于实现。但该方法是在影像的空间域直接对像素 进行操作，当两幅影像存在配准误差、或者重叠区域纹理结构丰富的时候，容易产生鬼影、 模糊现象，减弱影像纹理细节特征的表现。
[0005] 基于最优拼接线方法旨在原始影像对上寻找一条灰度、几何结构差异较小的最优 路径。虽然该方法可W保证重叠区域的纹理结构一致，避免出现鬼影现象，但并不能保证在 克服鬼影的同时避免出现拼接缝。尤其当相邻影像存在明显亮度差异时，即使找到了最优 拼接线也不能完全消除拼接缝，无法实现重叠区域的自然过渡。
[0006] 通过对两类方法的分析，当待拼接的影像对纹理结构丰富，或者存在配准误差、福 射亮度差异大的时候，依靠单一的拼接方法无法取得满意的效果。
【发明内容】
[0007] 本发明针对现有技术的不足，提供一种基于滤波分频的无人机影像拼接方法。通 过该方法可W将加权平滑融合方法和基于最优拼接线方法结合起来，充分利用两类方法的 优势，做到在克服鬼影的同时较好的消除拼接缝，实现无人机影像的无缝拼接。
[000引为了达到上述目的，本发明提供的技术方案是一种无人机影像滤波分频拼接方 法，包括W下步骤：
[0009] 步骤1，利用高斯低通滤波将待拼接的无人机影像对分解成低频影像对和高频影 像对；
[0010] 步骤2，对低频影像对采用加权平滑融合方法进行拼接处理；
[0011] 步骤3,利用改进的动态规划算法捜索最优拼接线完成高频影像对的拼接，其中改 进的动态规划算法捜索最优拼接线包括W下子步骤，
[0012] 步骤3.1，定义最优拼接线的能量函数值；
[0013] 步骤3.2,定义最优拼接线的捜索方向；
[0014] 步骤3.3，定义最优拼接线选取准则，确定最优拼接线；
[0015] 步骤4,将拼接好的低频影像对和高频影像对进行线性合成得到平滑无缝的无人 机拼接影像。
[0016] 而且，所述步骤2中，对低频影像对采用加权平滑融合方法进行拼接处理的实现方 式如下，
[0017]假设两幅低频影像在拼接影像中的重叠区域为(60，非重叠区域分别为〇)1和&2， 令Ii表示区域(61中影像的灰度，12表示区域(62中影像的灰度，Ip表示拼接影像的灰度，贝U 融合公式如下，
[001引
…
[0019] 其中，（x，y)是像素点在拼接影像中的行列号，A是权重因子，通过A可W强调离重 叠中屯、近的像素的贡献，减少图像边缘处像素的影响，W减小误差；
[0020] 假设W和H分别表示重叠区域的宽度和高度，若重叠区域不规则，则利用其最小外 接矩形进行计算，A W和Ah分别为重叠区域中像素点到左下角点在水平、竖直方向上的距 离；当两幅影像仅存在水平位移时，取(W- A w)/W;仅存在竖直位移时，取A= (H- A h)/H; 当两幅影像存在双向位移时，此时通过下式计算，
[0021]
巧
[0022] 其中，iii为竖直方向的参数调整因子，化为水平方向的参数调整因子。
[0023] 而且，所述步骤3.1中，定义最优拼接线的能量函数值的实现方式如下，
[0024] 定义能量函数为，
[0025]
[0026] .(4)
[0027] (5)
[002引
[0029]
17)
[0030] 其中，dm(x，y)表示像素点（x，y)在差值影像中S通道的灰度值;Cdif(x，y)是在差 值影像中像素点(x，y)周围7X7邻域内=通道灰度均值的最大值;Gm(x，y)是在差值影像中 像素点（X, y)的8方向梯度值的最大值，Sn表示8个方向的Sobe 1算子模板，n = 1对应0°，n = 2 对应45°，n = 3对应90°，n = 4对应 135°，n = 5对应 180°，n = 6对应225°，n = 7对应270°，n = 8 对应315°，Cedge(x，y)取像素点（x，y)周围7X7邻域内Gm(x，y)均值的最大值。
[0031] 而且，所述步骤3.2中，定义重叠区域内像素对左或右各1个像素点和正上方各3个 像素点共8个像素点为捜索方向。
[0032] 而且，所述步骤3.3中，定义选取最长线为最优拼接线，拼接线的走向要适应重叠 区域的形状，且尽量靠近重叠影像的中屯、线的准则，确定最优拼接线，实现方式如下，
[0033] 假设重叠区域差值影像大小为mXn，若重叠区域不规则，可利用其最小外接矩形 将其扩展为规则区域，能量矩阵为C，能量累积矩阵为Cuml和Cum2,方向矩阵为Dirl和Dir2， 线长矩阵为Lenl、Len2,四类矩阵的大小均等于重叠区域;能量矩阵C可由步骤3.1计算得 到，方向矩阵和能量累积矩阵根据像素点的能量值和其捜索方向上像素点能量值计算并比 较得到，〇^1、〇11111、1^6111和0山2、(：皿2、1^6112分别为从左向右和从右向左进行捜索时的方向 矩阵、能量累积矩阵和线长矩阵;最优拼接线的捜索流程如下，
[0034] 1)对于第1行的像素点（l，y)，令其能量累积值均等于其能量值，方向值为0,线长 值为1，即
[0035]
(8)
[0036] 2)对于第1行W外的其他行依次处理，记当前处理第X行，假设当前像素点的坐标 为(x，y)，首先从左至右进行捜索，设M为像素点(x，y)捜索方向像素点中最小能量累积值，
[00371
片）
[0038] 则能量累积值和线长值计算如下，
[0039]
(10)
[0040] 从左向右进行捜索时，比较CumlU,y)左侧1个像素点及正上方3个像素点和对应 的Cum2(x，y)正上方3个像素点共屯个方向像素点的能量累积值，选取能量累积值最小的像 素点为连接点，即M值最小的点，修改D i r 1 (X, y)为该点的方向值，并更新Cum 1 (X, y)为当前 计算得到的和值;特殊的，对位于左边界的点，只比较Cum l(x，y)正上方2个像素点和对应 的Cum2(x，y)正上方2个像素点共四个方向像素点的能量累积值;而右边界上的点，需比较 Cum2(x，y)正上方2个像素点和对应的CumlU, y)正上方2个像素点及左侧1个像素点共五个 方向像素点的能量累积值，设1为能量累计值最小的像素点对应的拼接线长度，更新Lenl (x，y)的值为1+1;
[0041] 然后从右至左进行扫描，其过程与从左向右相似，但在逆向扫描的过程中，需要比 较化m2(x，y)右侧1个像素点及正上方3个像素点和对应的CumlU, y)正上方3个像素点共7 个方向的像素点，找到能量累积值最小的像素点后，W该点作为连接点，修改Dir2(x，y)为 该点的方向值，更新化m2 (X，y)为当前计算得到的和值，并更新Len2 (X，y)的值为1+1;特殊 的，对于右边界上的点，比较Cum2(X，y)正上方2个像素点和对应的Cuml (X，y)正上方2个像 素点共四个方向像素点的能量累积值;而位于左边界的点，只比较Cum l(x，y)正上方2个像 素点和对应的Cum2(x，y)正上方2个像素点及右侧1个像素点共五个方向像素点能量累积 值，当执行到最左侧的时候，该行结束，跳到下一行继续执行从左至右的扫描；
[0042] 3)当最后一行扫描结束后，得到了最终的能量累积矩阵、方向矩阵和线长矩阵，统 计Lenl和Len2中最后一行线长值的大小，将最大线长值对应的点设置为终点；
[0043] 4)根据方向矩阵由终止点开始进行逆向寻找，并标识找到的点，当执行到第一行 时寻找结束，最优拼接线为所有被标识的点的连线。
[0044] 与现有技术相比，本发明利用加权平滑融合方法对低频影像对进行拼接处理，可 W很好的消除影像中的拼接缝，而基于最优拼接线对高频影像对的拼接则可避免拼接影像 出现的鬼影现象，最大程度的保证拼接线两侧纹理结构的一致性。
【附图说明】
[0045] 图1为本发明实施例的流程图。
[0046] 图2为本发明实施例低频影像对的加权平滑融合规则图，其中图2a为仅具有水平 位移的示意图，图化为具有双向位移的示意图。
[0047] 图3为本发明实施例最优拼接线捜索方法示意图，其中图3a为本发明实施例最优 拼接线捜索方法从左到右捜索时Dirl的取值示意图，图3b为本发明实施例最优拼接线捜索 方法从右到左捜索时Dir2的取值示意图。
【具体实施方式】
[0048] 本发明提供一种滤波分频的无人机影像拼接方法，首先利用高斯低通滤波将待拼 接的无人机影像对分解成低频影像对和高频影像对，然后对低频影像对采用加权平滑融合 方法进行拼接，利用改进的动态规划算法捜索最优拼接线完成高频影像对的拼接，最后将 拼接好的低频影像和高频影像进行线性合成得到平滑无缝的无人机拼接影像。
[0049] 下面结合附图和实施例对本发明的技术方案作进一步说明。
[0050] 如图1所示，本发明实施例的流程包括W下步骤：
[0051] 步骤1，利用高斯低通滤波将待拼接的无人机影像对分解成低频影像对和高频影 像对。
[0052] 本发明利用高斯低通滤波对待拼接的无人机影像对进行分解，可W分别获得其 高频、低频部分影像。低频部分影像主要反映无人机影像的色调、亮度信息，利用加权平滑 融合方法可W很好的调整拼接影像的整体色调，消除拼接缝，实现重叠区域的自然过渡。高 频部分影像主要包含无人机影像的轮廓、几何结构信息，直接利用加权平滑融合方法进行 拼接容易出现鬼影、模糊现象，而通过动态规划算法捜索一条最优拼接线，使其两侧色调、 亮度、几何结构差异最小，可W有效地减少鬼影、模糊现象的出现。
[005；3] [ 1 ]Wang M, Pan J , Chen S Q , et al. A method of removing the uneven illumination phenomenon for optical remote sensing image[C]. IEEE international geoscience remote sensing symposium(IGARSS), Seoul, SOUTH Korea, 2005.
[0054] 步骤2,对低频影像对采用加权平滑融合方法进行拼接处理，具体实施如下：
[0055] 对分解得到的低频影像对，采用加权平滑融合方法进行拼接。如图2a和图化所示， 假设两幅低频影像在拼接影像中的重叠区域为(60，非重叠区域分别为(61和(62，令Ii表示 区域(61中影像的灰度，12表示区域(62中影像的灰度，Ip表示拼接影像的灰度，则融合公式 为：
[0化6] (1)
[0057] 其中，（x，y)是像素点在拼接影像中的行列号，A是权重因子，通过A可W强调离重 叠中屯、近的像素的贡献，减少图像边缘处像素的影响，W减小误差。
[0058] 假设W和H分别表示重叠区域的宽度和高度(若重叠区域不规则，则利用其最小外 接矩形进行计算），A W和A h分别为重叠区域中像素点到左下角点在水平、竖直方向上的距 离。当两幅影像仅存在水平位移时，取人=(W- A w)/W;仅存在竖直位移时，取A= (H- A h)/H; 当两幅影像存在双向位移时，需要考虑水平、竖直两个方向的位移，设m为竖直方向的参数 调整因子，为水平方向的参数调整因子，此时通过下式计算：
[0059]
但)
[0060] 步骤3,利用改进的动态规划算法捜索最优拼接线完成高频影像对的拼接；
[0061] 本发明在Duplaquet提出的拼接线捜索方法基础上，重新定义了最优拼接线的能 量准则，并增加了最优拼接线的捜索方向，提出了新的最优拼接线选取准则，选取最优拼接 线完成高频影像对的拼接。具体包括W下步骤：
[0062] 步骤3.1，定义最优拼接线的能量函数值；
[0063] 假设Imgl和Img2是两幅待拼接的原始影像，本发明的定义能量函数为：
[0064] C(x,y) =Cdif (x,y)+Cedge(x,y) (3)
[00化]其中；
[0066]
[0067]
[006引
[0069] 巧）
[0070] _^v下um、义,：yy-(X/JV|分承,林、义,：y"工圧'|且巧ク'|分下二《^^a?Ui口。八化'l且。レdlf、义,：yy是在差值 影像中像素点(x，y)周围7X7邻域内=通道灰度均值的最大值。Gm(x，y)是在差值影像中像 素点（x，y)的8方向梯度值的最大值，如式(7)所示Sn表示8个方向的Sobel算子模板(n = l对 应0°，n = 2对应45°，n = 3对应90°，n = 4对应 135°，n = 5对应 180°，n = 6对应225°，n = 7对应 270°，n = 8对应315°)点<186山7)取像素点^，7)周围7乂7邻域内6。山7)均值的最大值。在 计算CedgeU, y)时加入邻域结构的变化信息，可提前预测是否靠近结构变化大的区域，为拼 接线的下一步捜索提供转变的空间，提前避免出现穿越现象。
[0071] 步骤3.2,定义最优拼接线的捜索方向；
[0072] DupIaquet提出的方法只检测当前点下方的3个相邻方向，当重叠区域有房屋、树 木等高大建筑物时，由于房顶像点的变形与地面点的变形不一致，拼接线极有可能会穿过 房屋，出现错位现象。为避免运种情况的发生，本发明提出一种双通道捜索拼接线的思路， 检测重叠区域内像素对左或右各1个像素点和正上方各3个像素点共8个像素点为捜索方 向。如图3a和图3b所示，设点P为当前点，重新定义八个比较方向，分别为1、2、3、4、11、22、 33、66。
[0073] 步骤3.3，定义最优拼接线选取准则，确定最优拼接线；
[0074] Duplaquet提出的拼接线捜索方法可W保证得到的备选拼接线长度均相等，并选 取能量值最小的拼接线为最优拼接线。徐亚明改进了 Duplaquet的捜索方法，且认为能量值 与拼接线长度的比值，即平均强度值最小者最优。本发明认为刻意的追求拼接线的能量值 最小，或者平均强度值最小并不能作为最优拼接线的选取准则，因为二者在寻线的过程中 都会限制拼接线的长度，当遇到的障碍物面积较大或者分布密集的时候，不可避免会穿越 而过。况且拼接线的捜索本就是一个能量最小值聚合的过程，每一条线都是由邻域中能量 值最小的像素点组合而成。因此，本发明提出选择最长线为最优拼接线，此外拼接线的走向 要适应重叠区域的形状，且尽量靠近重叠影像的中屯、线。
[0075] 假设重叠区域差值影像大小为mXn(若重叠区域不规则，可利用其最小外接矩形 将其扩展为规则区域），能量矩阵为C，能量累积矩阵为Cuml和Cum2,方向矩阵为Dirl和 Dir2,线长矩阵为Lenl、Len2,四类矩阵的大小均等于重叠区域。能量矩阵C可由式(3)计算 得到，方向矩阵和能量累积矩阵根据像素点的能量值和其捜索方向上像素点能量值计算并 比较得到，方向矩阵与捜索方向对应，其元素取值只有8个数值，DirlXuml、Lenl和Dir2、 Cum2、Len2分别为从左向右和从右向左进行捜索时的方向矩阵、能量累积矩阵和线长矩阵。 本发明改进的最优拼接线的捜索流程如下：
[0076] 1)对于第1行的像素点（l，y)，令其能量累积矩阵值均等于其能量值，方向值为0， 线长值为1，即
[0077]
巧）
[007引2似第X行为例，阐述其他行的执行过程。假设当前点的坐标为(x，y)，首先从左至右进行桿去.巧'M兩條畫占 r Y V、巧去市向條畫占由爲/1、能吾受巧估.
[0079] (9)
[0080]
[0081]
[0082] 如图3a所示，从左向右进行捜索时，比较CumlU, y)左侧1个像素点及正上方3个像 素点（1、2、3、4)和对应的Cum2(x，y)正上方3个像素点（11、22、33)共屯个方向像素点的能量 累积值，选取能量累积值最小的像素点为连接点，即M值最小的点，修改DirlU,y)为该点的 方向值(运里的方向值指的是捜索方向1、2、3、4、11、22、33中的某一个值），并更新Cuml(X， y)为当前计算得到的和值。特殊的，对位于左边界的点，只比较Cum l(x，y)正上方2个像素 点(2、3)和对应的Oim2 (X，y)正上方2个像素点（22、33)共四个方向像素点能量累积值;而右 边界上的点，需比较Cum2 (X，y)正上方2个像素点（11、22)和对应的化ml (X，y)正上方2个像 素点及左侧1个像素点（1、2、4)共五个方向像素点的能量累积值。设1为能量累计值最小的 像素点对应的拼接线长度，更新Lenl (X，y)的值为1+1。
[0083] 如图3b所示，然后从右至左进行扫描，其过程与从左向右相似，但在逆向扫描的过 程中，需要比较Cum2(x，y)正上方3个像素点及右侦Ul个像素点（11、22、33、66)和对应的Cuml (x，y)正上方3个像素点（1、2、3)共屯个方向的像素点，找到能量累积值最小的像素点后，W 该点作为连接点，修改Dir2(x，y)为该点的方向值(运里的方向值指的是捜索方向1、2、3、 11、22、33、66中的某一个值），更新加1112^,7)为当前计算得到的和值，并更新1^6]12^,7)的 值为1 + 1。特殊的，对于右边界上的点，比较Cum2(x，y)正上方2个像素点（11、22)和对应的 CumlU,y)正上方2个像素点（1、2)共个四方向像素点的能量累积值;而位于左边界的点，只 比较Oim l(x，y)正上方2个像素点（2、3)和对应的Cum2(x，y)正上方2个像素点及右侧1个像 素点（22、33、66)共五个方向像素点能量累积值。当执行到最左侧的时候，该行结束，I?兆到下 一行继续执行从左至右的扫描。
[0084] 3)当最后一行扫描结束后，得到了最终的能量累积矩阵（由像素点的能量累积值 按行列构成）、方向矩阵（由像素点的方向值按行列构成)和线长矩阵（由像素点的线长值按 行列构成）。统计Lenl和Len2中最后一行线长值的大小，将最大线长值对应的点设置为终 点。
[0085] 4)根据方向矩阵由终止点开始进行逆向寻找，并标识找到的点。当执行到第一行 时寻找结束，最优拼接线为所有被标识的点的连线。
[0086] 步骤4,线性合成
[0087] 最后将拼接好的高频、低频影像对进行线性合成W得到平滑无缝的无人机拼接影 像。
[0088] [2]Ting Lu,Shutao Li,Wei Fu.Fusion Based Seamless Mosaic for Remote Sensing Images[J].Sensing and Imaging,2014,15(1):1-14.参见图 1，无人机影像Imgl 和Img2分别进行高斯低通滤波后得到相应低频影像和高频影像，分别经加权平滑融合获取 拼接影像低频部分，经改进的动态规划捜索拼接线方法获取拼接影像高频部分，最后线性 合成拼接影像。
[0089] 具体实施时，W上流程可采用计算机软件技术实现自动运行流程。
[0090] 本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领 域的技术人员可W对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替 代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
【主权项】
1. 一种基于滤波分频的无人机影像拼接方法，其特征在于，包括如下步骤： 步骤1，利用高斯低通滤波将待拼接的无人机影像对分解成低频影像对和高频影像对； 步骤2，对低频影像对采用加权平滑融合方法进行拼接处理； 步骤3,利用改进的动态规划算法搜索最优拼接线完成高频影像对的拼接，其中改进的 动态规划算法搜索最优拼接线包括以下子步骤， 步骤3.1，定义最优拼接线的能量函数值； 步骤3.2，定义最优拼接线的搜索方向； 步骤3.3，定义最优拼接线选取准则，确定最优拼接线； 步骤4,将拼接好的低频影像对和高频影像对进行线性合成得到平滑无缝的无人机拼 接影像。2. 如权利要求1所述的一种基于滤波分频的无人机影像拼接方法，其特征在于:所述步 骤2中，对低频影像对采用加权平滑融合方法进行拼接处理的实现方式如下， 假设两幅低频影像在拼接影像中的重叠区域为Φ〇,非重叠区域分别为Φ#ΡΦ2,令1:表 示区域Φ冲影像的灰度，12表示区域Φ2中影像的灰度，ΙΡ表示拼接影像的灰度，则融合公 式如下，(1) 其中，（x，y)是像素点在拼接影像中的行列号，λ是权重因子，通过λ可以强调离重叠中 心近的像素的贡献，减少图像边缘处像素的影响，以减小误差； 假设W和Η分别表示重叠区域的宽度和高度，若重叠区域不规则，则利用其最小外接矩 形进行计算，Aw和Ah分别为重叠区域中像素点到左下角点在水平、竖直方向上的距离；当 两幅影像仅存在水平位移时，取λ= (W- △ w)/W;仅存在竖直位移时，取λ= (H- △ h)/H;当两 幅影像存在双向位移时，此时通过下式计算， (2) * 1 \ 广/ * -\ / i J 其中，&为竖直方向的参数调整因子，μ2为水平方向的参数调整因子。 3 .如权利要求2所述的一种基于滤波分频的无人机影像拼接方法，其特征在于:所述步 骤3.1中，定义最优拼接线的能量函数值的实现方式如下， 定义能量函数为， C(x,y)=Cdif(x,y)+Cedge(x,y) (3) / κ、、. (4)其中，cU( x，y)表不像素点（x，y)在差值影像中二通道的灰度值;Cdif (x，y)是在差值影像 中像素点(x，y)周围7X7邻域内三通道灰度均值的最大值;Gm(x，y)是在差值影像中像素点 (X，y)的8方向梯度值的最大值，S n表示8个方向的Sobe 1算子模板，η = 1对应0°，η = 2对应 45°，η = 3对应90°，η = 4对应 135°，η = 5对应 180°，η = 6对应225°，η = 7对应270°，η = 8对应 315°，Cedge(x，y)取像素点（x，y)周围7X7邻域内Gm(x，y)均值的最大值。4. 如权利要求3所述的一种基于滤波分频的无人机影像拼接方法，其特征在于:所述步 骤3.2中，定义重叠区域内像素对左或右各1个像素点和正上方各3个像素点共8个像素点为 搜索方向。5. 如权利要求4所述的一种基于滤波分频的无人机影像拼接方法，其特征在于:所述步 骤3.3中，定义选取最长线为最优拼接线，拼接线的走向要适应重叠区域的形状，且尽量靠 近重叠影像的中心线的准则，确定最优拼接线，实现方式如下， 假设重叠区域差值影像大小为mXn，若重叠区域不规则，可利用其最小外接矩形将其 扩展为规则区域，能量矩阵为C，能量累积矩阵为Cuml和Cum2，方向矩阵为Dirl和Dir2，线长 矩阵为Lenl、Len2,四类矩阵的大小均等于重叠区域;能量矩阵C可由步骤3.1计算得到，方 向矩阵和能量累积矩阵根据像素点的能量值和其搜索方向上像素点能量值计算并比较得 至1」，0化1、〇11111、1^111和0化2、〇11112、1^112分别为从左向右和从右向左进行搜索时的方向矩 阵、能量累积矩阵和线长矩阵;最优拼接线的搜索流程如下， 1) 对于第1行的像素点（l，y)，令其能量累积值均等于其能量值，方向值为〇,线长值为 1，即觸 2) 对于第1行以外的其他行依次处理，记当前处理第X行，假设当前像素点的坐标为(X， y)，首先从左至右进行搜索，设Μ为像素点(x，y)搜索方向像素点中最小能量累积值， 则能量累积值和线长值计算如下， (9)从左向右进行搜索时，比较Cuml(x，y)左侧1个像素点及正上方3个像素点和对应的 Cum2(X，y)正上方3个像素点共七个方向像素点的能量累积值，选取能量累积值最小的像素 点为连接点，即Μ值最小的点，修改D i r 1 (X，y)为该点的方向值，并更新Cum 1 (X，y)为当前计 算得到的和值;特殊的，对位于左边界的点，只比较Cum l(x，y)正上方2个像素点和对应的 Cum2(X，y)正上方2个像素点共四个方向像素点的能量累积值；而右边界上的点，需比较 Cum2(X，y)正上方2个像素点和对应的Cuml(X，y)正上方2个像素点及左侧1个像素点共五个 方向像素点的能量累积值，设1为能量累计值最小的像素点对应的拼接线长度，更新Lenl (x，y)的值为1+1; 然后从右至左进行扫描，其过程与从左向右相似，但在逆向扫描的过程中，需要比较 Cum2(x，y)右侧1个像素点及正上方3个像素点和对应的Cuml(x，y)正上方3个像素点共7个 方向的像素点，找到能量累积值最小的像素点后，以该点作为连接点，修改Dir2(x，y)为该 点的方向值，更新Cum2 (X，y)为当前计算得到的和值，并更新Len2 (X，y)的值为1+1;特殊的， 对于右边界上的点，比较Cum2(x，y)正上方2个像素点和对应的Cuml(x，y)正上方2个像素点 共四个方向像素点的能量累积值;而位于左边界的点，只比较Cum l(x，y)正上方2个像素点 和对应的Cum2(x，y)正上方2个像素点及右侧1个像素点共五个方向像素点能量累积值，当 执行到最左侧的时候，该行结束，跳到下一行继续执行从左至右的扫描； 3) 当最后一行扫描结束后，得到了最终的能量累积矩阵、方向矩阵和线长矩阵，统计 Lenl和Len2中最后一行线长值的大小，将最大线长值对应的点设置为终点； 4) 根据方向矩阵由终止点开始进行逆向寻找，并标识找到的点，当执行到第一行时寻 找结束，最优拼接线为所有被标识的点的连线。
【文档编号】G06T3/40GK105957018SQ201610556890
【公开日】2016年9月21日
【申请日】2016年7月15日
【发明人】赵双明, 李婉
【申请人】武汉大学