【技术实现步骤摘要】
本专利技术涉及图像处理和图像配准，尤其是一种基于相位谱与幅度谱重组的一维相位相关图像配准方法。
技术介绍
1、图像配准就是将同一场景的不同图像转换到同样的坐标系统中的过程；常常用于不同视场、不同成像模式或不同拍摄时间，但又有重叠或相似的两幅或多幅图像之间的匹配或对齐。
2、傅里叶变换将一个函数表示为不同频率的正弦和余弦函数的叠加；它将一个函数从时域(时间域)转换到频域(频率域)，并展示了该函数所包含的不同频率成分的振幅和相位信息；相位谱描述了信号在不同频率上的相位信息，而幅度谱描述了信号在不同频率上的振幅大小或能量分布；通过分析这两个谱，可以获得关于信号的频率特性和结构信息；相位相关是一种通过比较相位信息来推断图像或信号之间位移量的方法；它通过计算两个图像或信号的相位谱，并利用互功率谱计算它们之间的相位相关函数并得到相位相关性的结果，该结果表示两个频域图像之间的相似程度；该值越接近1，表示两个图像的相位谱越相似；与传统的图像配准方法相比，相位相关算法具有良好的精度和可靠性等优势，因此被广泛地研究和应用在实际产品。
3、然而，在现有技术中，由于图像的相位相关求偏移量需要对图像的相位谱进行二维计算，相较于一维计算更为耗时。当两幅图像的特征信息较少或背景噪声过多时，直接将二维图像转为一维信号并使用一维相位相关算法会引入二维图像幅度谱与噪声的干扰，导致计算得到的偏移量产生较大误差。因此，该方法的使用受到严格限制，不适合大范围推广应用。
技术实现思路
1、为了克服现有
2、本专利技术解决其技术问题所采用的技术方案是：
3、基于相位谱与幅度谱重组的一维相位相关图像配准方法，包括如下步骤：
4、步骤一：对二维图像进行平滑处理，再根据先验知识，对处理过的图像a(x,y)与图像b(x,y)选取重叠区域，a(x,y)上的重叠区域记为f(x,y)，b(x,y)上的重叠区域记为g(x,y)，记重叠区域在水平方向的长度为m，在竖直方向的长度为n，即0≤x＜m，0≤y＜n，从而减少冗余信息；
5、步骤二：分别对图像f(x,y)与图像g(x,y)进行傅里叶变换，分离其相位谱与振幅谱；傅里叶变换的数学表达式如下：
6、
7、其中，f(u,v)是频域中的复数函数，表示f(x,y)在频率域上的表示；e(-i2π(ux+vy))是二维复指数函数，u和v分别表示频域中的横向和纵向频率；
8、振幅谱af(u,v)可以通过对f(k)进行求模得到：
9、af(u,v)＝|f(u,v)|
10、相位谱φf(u,v)可以通过对f(k)进行求辐角得到：
11、φf(u,v)＝arg(f(u,v))
12、同样地，求得图像g(x,y)对应的频域中的复数函数g(u,v)与振幅谱ag(u,v)、相位谱φg(u,v)；
13、步骤三：将相位谱φf(u,v)、φg(u,v)分别与振幅谱af(u,v)或ag(u,v)融合；
14、作为一种举例说明，取振幅谱af(u,v)，融合得到的频域中的复合函数为ff(u,v)、gf(u,v),公式为：
15、
16、
17、将频域中的复合函数ff(u,v)、gf(u,v)转换回时域，公式为：
18、
19、
20、其中，m为重叠区域在水平方向的长度，n为在竖直方向的长度，0≤x＜m，0≤y＜n。
21、步骤四：分别对图像的x方向和y方向求导，以得到x方向与y方向的边缘信息，即计算出图像ff(x,y)的水平梯度fdx(x,y)和垂直梯度fdy(x,y)；同理，计算出图像gf(x,y)的水平梯度gdx(x,y)和垂直梯度gdy(x,y)；
22、步骤五：通过按行/列累加、求行/列平均值的计算方法，将图像fdx(x,y)投影为x轴方向的一维信号，记为fx(y)，将图像fdy(x,y)投影为y轴方向的一维信号，记为fy(x)；同理，将图像gdx(x,y)投影到x轴上记为gx(y)，以及将图像gdy(x,y)投影到y轴上记为gy(x)；
23、步骤六：对fx(y)、fy(x)、gx(y)、gy(x)可进行预处理从而减小噪声等因素对配准的影响；fx(y)、fy(x)、gx(y)、gy(x)进行一维傅里叶变换，将一维信号转为频域空间fx(k)、fy(k)、gx(k)、gy(k)；根据fx(k)和gx(k)求取水平方向的互功率谱同理求得垂直方向的互功率谱对求得的分别做反傅里叶变换得到狄拉克函数通过寻找峰值的坐标即可找到偏移量。当fx(y)与gx(y)、fy(x)与gy(x)相似时，将出现明显的峰值ax、ay，峰值的位置δxδy，即为fx(y)与gx(y)、fy(x)与gy(x)的相对偏移。
24、本专利技术的有益效果：
25、本专利技术克服了两幅图像因为特征点过少和图像存在噪音、暗角导致计算出的偏移量误差较大的问题，并提升了计算速度；
26、为了避免冗余信息的干扰，只选取重叠区域进行处理；
27、对重叠区域分别进行傅里叶变换并分离它们的相位谱与幅度谱；将两个相位谱分别与其中一个幅度谱结合，并通过逆傅里叶变换得到它们在时域的表示，相对于直接使用相位谱的方法减小了相位谱存在噪声相位信息而对结果产生的误差；
28、对重新结合相位谱与幅度谱的图像计算出x方向和y方向的梯度图像，通过按行/列相加、求行/列平均值等计算方法，将梯度图像投影到相关方向轴转为一维信号，相对于直接求一维信号的方法消除了因引入幅度谱产生的误差；
29、对一维信号做去噪声、去背景等处理；利用相位相关算法应用在处理过的一维信号上，计算得出两幅图像的偏移量，解决了由于图像特征信息过少、光照不均、存在冗余信息与噪声等因素找不准图像偏移量的问题，并提升了计算速度。
本文档来自技高网...【技术保护点】
1.基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，对图像求其频域表示采用：二维离散傅里叶变换DFT或快速傅里叶变换FFT方法。
3.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，对FF(u，v)、GF(u，v)求其时域表示采用：逆傅里叶变换IDFT或逆快速傅里叶变换IFFT。
4.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，所述水平梯度fdx(x，y)、gdx(x，y)和垂直梯度fdy(x，y)、gdy(x，y)是指某像素在x和y方向上的变化率。
5.根据权利要求4所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，计算图像的水平梯度和垂直梯度的方法包括：直接计算法和算子求解法。
6.根据权利要求5所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，所述直接计算法采用：当前像素相邻的像素比较；
...【技术特征摘要】
1.基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，对图像求其频域表示采用：二维离散傅里叶变换dft或快速傅里叶变换fft方法。
3.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，对ff(u，v)、gf(u，v)求其时域表示采用：逆傅里叶变换idft或逆快速傅里叶变换ifft。
4.根据权利要求1所述的基于相位谱与幅度谱重组的一维相位相关图像配准方法，其特征在于，所述水平梯度fdx(x，y)、gdx(x，y)和垂直梯度fdy(x，y)、gdy(x，y)是指某像素在x和y方向上的变化率。
5.根据权利要求4所述的基于相位...
【专利技术属性】
技术研发人员：窦晓慧，洪喜，王楠，
申请(专利权)人：长光辰英杭州科学仪器有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。