一种基于TS201的归一化互相关快速算法

0 引言

　　图像匹配指在已知目标基准图的子图集合中，寻找与实时图像最相似的子图，以达到目标识别与定位目的的图像处理技术。由于归一化互相关算法的实现方案简单，对灰度值的线性变化具有适应性、抗白噪声能力强，以及匹配性能稳定的特点，已成为图像匹配技术中应用最为广泛的匹配算法。该算法不仅可以应用于民用图像处理技术领域，也可以应用于SAR成像和红外成像制导系统等军事领域当中。虽然归一化互相关算法的匹配能力强，性能稳定，但该算法的计算量较大，不易应用于实时系统当中。针对这一现象，讨论了一种基于ADSP-TS201处理器，在保证全图遍历条件下，使用递推与多模板思想构建的归一化互相关快速算法。实验证实，该快速算法在保证算法性能的同时，有效地提高了算法的执行效率。

　　l 归一化互相关算法实现过程分析

　　归一化互相关算法是基于灰度图像匹配技术中最为经典的算法，具有抗白噪声能力强，匹配准确性高的突出特点，但该算法的计算量较大，严重影响了算法应用的实时性。

　　设F(x，y)为一幅M×N的基准图像，B(x，y)是一幅m×n的实时图像、Sx,y是一幅与实时图像做互相关运算的匹配子图，且m<M，n<N。

　　标准归一化互相关算法计算公式如下：

　　观察上述公式可知，对应数据的减法运算、乘法运算和累加运算占据了标准算法中绝大部分的计算量。针对TS20l处理器，可以借助其宽总线、大内存及并行指令的特点实现图像数据的高速存取和并发计算，提升标准算法的执行效率，但该处理器仅支持32位取值，计算单元不支持8位数据乘法的特性，也给提升标准算法的执行效率带来了一定的限制。实验发现，导致TS201执行归一化互相关算法效率下降的主要原因有两个：一个是标准算法本身计算量庞大；另一个是处理器32位取值的特性使匹配子图数据的获取困难。

　　2 递推多模板方案快速算法

　　2．1 递推方案

　　递推多模板方案快速算法首先要解决的问题就是归一化互相关算法本身计算量庞大的问题。观察式(1)可知，式(2)为匹配子图与实时图协方差计算过程，式(3)为匹配子图方差的计算过程，式(4)为实时图的方差计算过程。对于实时图的方差计算部分，在整个匹配过程中只需要计算一次，而对于匹配子图的方差计算和两张图像的协方差计算则需要多次重复执行，共需要约(M-m)x(N-n)×3×m×n次加乘计算，计算量庞大。本文采用匹配子图递推方案来减少计算量。

　　结合实际图像，最大互相关算法实际是将实时图在匹配基准图像上做逐点逐行平移，并与其覆盖的匹配子图做相关运算，求取相关系数并比对，以获取最佳匹配位置。在实时图平移的过程中，其覆盖的新匹配子图数据遵循如下递推原则：

　　(1)当实时图进行水平平行移动时，因Sx,y+1是Sx,y在基准图中右移1列的位置所对应的新匹配子图，即Sx,y+1的前n-1列是Sx,y的后n-1列，有：

　　(2)当实时图进行垂直平行移动时，因Sx+1,y是Sx,y在基准图中下移1行的位置所对应的新匹配子图，即Sx+1,y的前m-1行是Sx,y的后m-1行，有：

　　利用此递推原则，可以将归一化互相关式中的式(3)进行简化。观察式(3)中数据计算结果是以各个数据平方和的累加形式出现的，所以先将式(5)和式(6)进行平方运算，得到对应的式(7)和式(8)。

　　使用式(7)和式(8)对式(3)进行变换，得到全图遍历过程中使用的匹配子图递推方差计算：

　　在全图遍历的过程中，式(9)每次需要得到的计算结果只是在上次累加结果的基础上，减去移出数据的累加值和增加移人数据的累加值，中间重复数据的累加结果不再重复计算，从而减少了计算量。匹配子图递推方案将会使单次子图方差计算的计算量由2×m×n次加乘减少至2×m或者2×n次加乘。

　　2．2 多模板方案

　　递推多模板方案的快速算法需要解决的另一个重要问题是匹配子图的数据准备问题。有两个不利因素导致产生此问题，一个是匹配子图数据准备过程重复进行(M-m)×(N-n)次，占据了大量处理时钟；另一个由式(2)的分析可知，算法运行过程中，每次相关匹配运算需要将实时图数据块与匹配子图数据块中的像素一一对应进行乘法运算。由于TS201自身只支持32位取值(字取值)，所以无法将直接获取数据的方法获得的起始像素不位于字首的匹配子图数据块用于计算，也即不能利用DMA使用双缓冲的常规方式准备数据。

　　本文采用多模板计算方案来解决匹配子图数据的准备问题。

多模板计算方案首先将处理器定位于4字取值工作模式，加载并行指令，使处理器的计算能力最大化；其次将实时图予以适当的横向放大，构建多个与匹配子图等高的大实时图块，并分别在大实时图块范围内将真实实时