[发明专利]基于FFT的碎纸片纵切拼接复原算法

说明书

技术领域

本发明主要运用于图像拼接（文字拼接）领域，设计了一种基于FFT的碎纸片纵切拼接复原算法。  

背景技术

快速傅里叶变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅里叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的发现，但是对于在计算机系统或者说数字系统中应用离散傅立叶变换，可以说是进了一大步。

快速傅氏变换算法可以运用到图像拼接的领域。 

发明内容

将原理图编码，编码后使之转化为矩阵。其中，为了减少算法的计算量和时间复杂度，编码采用的编码方式为，将所有白色的像素设定为0，其他均为1。编码后之后，得到矩阵，便于后续处理。

将编码后的数据进行识别然后进行FFT处理，对各个碎片进行频域变化，进行频谱分析。求其矩阵边缘相互之间的自相关性，将自相关性最大的两张碎片进行拼接。

附图说明

图1为算法流程图。

具体实施方案

为了探寻两个相邻碎片间的关系，本算法将各个碎片编码后矩阵的边缘列进行FFT快速傅里叶变换，因为最终图片的最左边全为白色，所以首先可以得到原图最左边一列，再将该碎片的最右列和其他碎片的最左列逐一进行互相关系数的计算，得出一族互相关系数，选取其中结果最大的一列的碎片进行匹配，得到最终结果。

傅里叶原理表明：任何连续测量的时序或信号，都可以表示为不同频率的正弦波信号的无限叠加。而根据该原理创立的傅里叶变换算法利用直接测量到的原始信号，以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。

对于信号的分析，很多情况下仅依靠时域是远远不够的，很多时域上不好分析的信号转入频域的分析中，立即变得便于分析。本算法依据分析的需要，分别将每个灰度矩阵的左上角和右上角看作起点，将0-1二值矩阵一列的取值看作y轴，将离散的点看作t轴，对其进行FFT处理。得到每个矩阵的边缘列的频谱。对它的特征进行提取，并分别计算它与其它碎纸片边缘的互相关系数，进行匹配。笔者认为，频谱特征接近，即互相关系数接近的两个碎纸片处在相邻的位置。

该算法的流程图如图1所示：

图片编码矩阵的FFT处理：

① FFT快速傅里叶变换与DFT离散傅里叶变换算法的比较：

对于任一离散序列                                                的点，其DFT计算公式如下所示：

设是基-2数，也就是2的整数次幂，即：，其中为正整数。这样，可以首先将序列前后对半分开，将点的  DFT写成前后两部分，如下式所示：

因为所以：

当为偶数时，为奇数时，所以可将进一步分解为偶数组和奇数组，为偶数时，

为奇数时，

令：

这样，和都是点的序列，将其分别带入下式就能清楚地看到，这两式表示的是两个点的DFT运算，如下式所示：

从上文可以发现，FFT算法对于DFT算法的改进，主要基于的取值特性，进行了一种合并，从而改进运算。DFT算法的运算量是与成正比的，,而利用的某些特性，可以将一个大点数DFT的运算分解为若干小点数的DFT运算组合，以减少运算量。该种FFT算法运算量与成正比，改进效果明显，所以本算法采用该种FFT算法进行计算。

② FFT算法对于碎片的处理与匹配：

对图片进行二值编码，首先确定出最左的碎片，之后再将该碎片矩阵的最右列进行FFT处理得到其频谱，之后选取其他各碎片对应矩阵最左列的频谱与其最接近的进行匹配。

但是，这样进行碎片的匹配需要添加较多的人工干预，为了让计算机进行自动的匹配，本专利接下来引入了相关性去对频谱的匹配程度进行定量的衡量。

归一化互相关系数模型