资源描述：
《1022019855-李盼盼-基于傅里叶变换的图像拼接》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、基于傅里叶变换的图像拼接李盼盼(河南大学数学与信息科学学院开封475004)摘要：图像拼接在图像处理中应用广泛。本文主要研究了傅里叶变换，将图像变换到频域进行分析。利用时域的相对位移等效于频域互功率谱的相位，可以计算出图片的相对偏移量。然后，利用渐进渐出的加权平均法方法消除拼接缝隙，简单且能达到较好的结果。关键词：图像拼接；傅里叶变换；对数极坐标变换1引言图像拼接（mosaic）技术就是将一组重叠图像的集合拼接成一幅大的分辨率更高的图片[1]。图像拼接在我们的生活应用广泛。由于全景图的硬件设备（如广角镜头）比较昂贵，普遍应用不太可能。图像拼接技术就可以让我们不
2、依赖全景图的硬件，而得到大场景的照片。2图像拼接步骤图像拼接框架如图2-1所示，图像拼接的关键技术是图像配准与图像融合两个步骤[2]。图像采集图像融合图像配准图像预处理图1图像拼接流程示意图2.1图像采集在图像拼接之前，需要采集需要的图片。对于采集的图片，相邻的图片之间必须有足够的重叠部分。为了在后续的图像配准中获得较好的图像匹配效果，通常在相同环境下采集，以减少环境因素带来的图像相同区域之间的差异。图像采集一般通过旋转，平移相机等方式获取。本文通过相机平移采集照片。2.2图像预处理9图像预处理的主要目的是为了减少环境及其他因素对图像配准的影响，提高配准精度和
3、拼接处理速度。图像预处理主要承担两项任务[3]：（1）对使用的相机进行校准，获取其内参数；（2）对图像失真进行校正，获取相机的校正参数。2.3图像配准数字图像拼接的质量主要取决于图像的配准精度。本文利用傅里叶变换后，时域的平移量等效于互功率谱的相位，其傅里叶反变换是平移运动坐标上的脉冲，搜索最大值的位置（x0，y0）就是两幅图像的相对平移量。2.4图像融合图像融合技术主要分为三个层次：像素级，特征级，决策级[4]。具体参考文献[5]。常用的图像融合算法有：（1）直接平均法。直接平均法是将两幅图片的重叠区域像素的灰度值叠加，然后求其平均值，用其作为融合后图像的像
4、素灰度值。公式为：.（2）加权平均法加权平均法是把图像重叠区域的像素值进行加权相加，然后计算平均值。公式为：,f,f1,f2意义同上，w1+w2=1，w1>0，w2>0.本文融合采用渐进渐出法（w1=d1,w2=d2，在重叠区域，d1由1渐变为0，d2由0渐变成1）。3相位相关法二维傅里叶变换定义9定义二维傅里叶变换为：,其中，u=1,2……，m;v=1,2……，n.傅里叶逆变换定义为：,其中，x=1，…….m;y=1,……，n.二维离散傅里叶变换的幅度谱定义为：.相位谱定义为：.能量谱定义为：.设f1和f2分别表示两幅待拼接的数字图像，第二副图像在第一副图像
5、的相对位置上有如下的坐标转换关系：.其中（x0，y0）表示相对平移量。对上式两边做傅里叶变换，则有下式：.由上式可知，它们之间的相位差可以等效为互功率谱的相位：.由于.对上式做傅里叶逆变换，相位差的傅里叶逆变换是平移坐标上的脉冲函数，只需找到最大值得位置，对应的坐标为两幅图的相对位移。4基于相位相关法的图像拼接算法94.1算法流程(1)导入图象001，,002,003，并进行傅里叶变换图片lib1图片lib29图片schooldoor1图片schooldoor2图2图像拼接示例(4幅源图)9图片相对偏移量计算函数为：function[x0,y0]=pyjs(p
6、ic1,pic2)[pic1xpic1y]=size(pic1);[pic2xpic2y]=size(pic2);x=max(pic1x,pic2x);y=max(pic1y,pic2y);p1=[zeros(pic1x,y-pic1y),pic1;zeros(x-pic1x,y)];p2=[pic2,zeros(pic2x,y-pic2y);zeros(x-pic2x,y)];f1=fft2(p1);f2=fft2(p2);glp=zeros(x,y);fb=f1.*conj(f2);fc=abs(fb)+eps;p3=ifft2(fb./fc);[d,e]
7、=max(p3);[f,g]=max(d);x0=e(g);y0=g;导入图像计算相对平移量程序如下:clear;clc;cp='E:matlabprogram';lib1=imread([cp,'lib1.png']);lib2=imread([cp,'lib2.png']);scd1=imread([cp,'schooldoor1.png']);scd2=imread([cp,'schooldoor2.png']);[x0,y0]=pyjs(lib1,lib2);9（3）图像融合，加权平均法的matlab融合函数如下：functionpicture=
8、ronghe(pic1,pic2,a,
