1022019855-李盼盼-基于傅里叶变换图像拼接.docVIP

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