CN102521834A - 采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 - Google Patents
采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 Download PDFInfo
- Publication number
- CN102521834A CN102521834A CN2011104109122A CN201110410912A CN102521834A CN 102521834 A CN102521834 A CN 102521834A CN 2011104109122 A CN2011104109122 A CN 2011104109122A CN 201110410912 A CN201110410912 A CN 201110410912A CN 102521834 A CN102521834 A CN 102521834A
- Authority
- CN
- China
- Prior art keywords
- alpha
- beta
- sin
- cot
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Image Analysis (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种采用极对数坐标表示的分数阶傅里叶变换(FRFT)的时频域图像配准算法,本发明将信号在分数阶Fourier域上的表示,同时融合了信号在时域和频域的信息,采用相位相关技术,将基准图像和待配准图像作FRFT变换,确定其平移参数,并通过对数-极坐标变换得到旋转、缩放等配准参数。本发明能够全面反映信号随时间变化的频率特征。
Description
技术领域
本发明涉及一种图像分析技术,具体涉及采用分数阶傅里叶变换的图像配准方法。
背景技术
基于傅里叶(Fourier)变换的图像自动配准方法,由于不需要寻找控制点和传感器参数进行图像自动配准的优点,被广泛应用于数字图像的实时处理中。
大多数信号的频率成分是随时间变化的,即非平稳信号。而傅里叶变换只能从整体上指出信号中曾经出现过的频率成分,不能展示信号的频率是如何随时间变化的。图像的二维傅里叶变换幅值在其零频率的较小邻域上呈现有较大的波动尖峰状,仅仅采用傅里叶变换已不能准确给出低频段上的近似值,这将大大降低该配准算法的稳定性,不利于频域配准方法进一步向存在平移、旋转和尺度的多谱段、多传感器图像配准问题扩展。
作为Fourier变换的一种广义形式,分数阶傅里叶变换(FRFT)可以解释为信号在时频平面内坐标轴绕原点逆时针旋转任意角度后构成的分数阶傅里叶域上的表示方式。如果信号的Fourier变换可看成将其在时间轴上逆时针旋转π/2到频率轴上的表示,则FRFT可以看成将信号在时间轴上逆时针旋转任意角度α到u轴上的表示(u轴被称为分数阶傅里叶域)。从本质上讲,信号在分数阶傅里叶域上的表示,同时融合了信号在时域和频域的信息,因此被认为是一种时频分析方法。
发明内容
本发明的目的在于将分数阶傅里叶变换(FRFT)引入到图像配准领域,提出了基于对数极坐标分数阶傅里叶变换的时频域配准算法。非平稳信号分析的时频联合分布可以将一维的时域信号映射到二维的时频域平面,全面反映信号随时间变化的频率特征。
为达到上述目的,本发明采用如下的技术方案:
一种采用极对数坐标表示的分数阶傅里叶变换的图像配准方法,该方法包括如下具体步骤:
(1)求出尺度因子和旋转角度
针对二维的情况,若基准图像I1和待配准图像I2具有平移、旋转和尺度关系
I1[s(xcosθ0+ysinθ0)+Δx,s(-xsinθ0+ycosθ0)+Δy]=I2(x,y)
其中θ0为旋转角,s为尺度因子,(Δx,Δy)为平移参数;
对上式r坐标轴取对数得到
采用上式和相位相关技术,在(logr,θ)坐标系下,先求出尺度因子s和旋转角度θ0;
(2)对图像做校正
对待配准图像做尺度因子s和旋转角度θ0的校正得到图像
(3)求出平移参数
设基准图像I1和待配准图像I2满足
I1(x+Δx,y+Δy)=I2(x,y)
式中*代表复共轭,对互功率谱Sxy(w)逆分数阶傅里叶变换可得到I1和I2的相位相关矩阵Rxy(x,y),相关矩阵峰值所对应的位置即为平移参数(Δx,Δy)。
由于自然界的大多数信号是非平稳信号,图像的功率谱幅值在其零频率的较小邻域上呈现有较大的波动尖峰状。采用分数傅立叶变换的配准方法,在时频平面内,将配准信号的相位相关映射到任意分数阶的时频域轴上。通过分数阶时频轴在时频域平面内的任意旋转,同时考虑了信号在频域和时间域的相关性,是一种时频分析方法。该方法可以准确的给出低频段上的近似值,从而提高配准算法的稳定性。
附图说明
以下结合附图和具体实施方式来进一步说明本发明。
图1为分数傅立叶时频域坐标示意图。
图2为相位相关检测示意图。
图3为极对数坐标示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
傅里叶域方法的基本思想是相位相关法的平移检测特性。图像在对数极坐标频域中,旋转和尺度关系转化为平移关系,可通过相位相关检测。鉴于傅里叶变换和分数阶傅里叶变换的关系,本发明提出将用于检测旋转和缩放的频域表示变为分数阶傅里叶变换。
本发明的方法主要包括以下几个步骤(参见图1至图3):
1.利用FRFT求出尺度因子和旋转角度
定义:从线性积分变换的角度,u′域的函数f(u′)的p阶分数阶傅里叶变换可以定义为一个线性积分运算:
其中Kp(u,u′)≡Aαexp[jπ(u2cotα-2uu′cscα+u′2cotα)],称为分数阶傅里叶变换的核函数, p≠2n,n是整数。
式(1)中α仅出现在三角函数的参数位置上,因此,以p(或α)为参数的定义是以4(或2π)为周期的,因此只需考察区间p∈[-2,2](或α∈[-π,π])即可。定义函数的零阶变换等于该函数本身,±2阶变换等于f(-u),并且核函数Kp(u,u′)在所有的p取值上连续。
假设基准图像I1和待配准图像I2具有平移、旋转和尺度关系
I1[s(xcosθ0+ysinθ0)+Δx,s(-xsinθ0+ycosθ0)+Δy]=I2(x,y) (2)
其中θ0为旋转角,s为尺度因子,(Δx,Δy)为平移参数。
(2)式在极坐标下的FRFT变换为
对上式中r坐标轴取对数
采用(5)式和相位相关技术,相位相关技术,与平移检测的方法相同,即对互功率谱求逆分数傅立叶变换,计算相位相关矩阵,通过检测相位相关矩阵的峰值,得到s和θ0可以求出尺度因子s和旋转角度θ0。
2.对图像做校正
3.求出平移参数
假设图像I1和I2满足
I1(x+Δx,y+Δy)=I2(x,y) (6)
式中*代表复共轭,对互功率谱Sxy(w)逆分数阶傅里叶变换可得到I1和I2的相位相关矩阵Rxy(x,y),相关矩阵峰值所对应的位置即为平移参数(Δx,Δy);
本发明可以解决非平稳信号图像配准中,传统的傅里叶变换配准方法只能从整体上指出信号中曾经出现过的频率成分、无法准确给出低频段近似值的缺陷。采用分数傅里叶变换的图像配准方法,在时频平面内,将配准信号的相位相关映射到任意分数阶的时频域轴上。通过分数阶时频轴在时频域平面内的任意旋转,同时融合了信号在时域和频域的信息,从而提高低频段信号的配准稳定性。
基于上述方案,其具体实施如下:
1.求出尺度因子和旋转角度
假设图像I1和I2具有平移、旋转和尺度关系,在经过FRFT变换后的关系为
二维分数阶傅里叶变换中的角度α和β直接影响了变换的性质,相应的分数阶傅里叶变换的参数p1和p2,变换范围均为(-2,2],本发明中将两个参数设为相同的值。
2.对图像做校正
3.求出平移参数
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (1)
1.一种采用极对数坐标表示的分数阶傅里叶变换的图像配准方法,其特征在于,所述方法包括如下具体步骤:
(1)求出尺度因子和旋转角度
针对二维的情况,若基准图像I1和待配准图像I2具有平移、旋转和尺度关系:
I1[s(xcosθ0+ysinθ0)+Δx,s(-xsinθ0+ycosθ0)+Δy]=I2(x,y)
其中θ0为旋转角,s为尺度因子,(Δx,Δy)为平移参数;
对上式r坐标轴取对数得到
采用上式和相位相关技术,在(1ogr,θ)坐标系下,先求出尺度因子s和旋转角度θ0;
(2)对图像做校正
对待配准图像做尺度因子s和旋转角度θ0的校正得到图像
(3)求出平移参数
设基准图像I1和待配准图像I2满足
I1(x+Δx,y+Δy)=I2(x,y)
式中*代表复共轭,对互功率谱Sxy(w)逆分数阶傅里叶变换可得到I1和I2的相位相关矩阵Rxy(x,y),相关矩阵峰值所对应的位置即为平移参数(Δx,Δy)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104109122A CN102521834A (zh) | 2011-12-12 | 2011-12-12 | 采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104109122A CN102521834A (zh) | 2011-12-12 | 2011-12-12 | 采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102521834A true CN102521834A (zh) | 2012-06-27 |
Family
ID=46292738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011104109122A Pending CN102521834A (zh) | 2011-12-12 | 2011-12-12 | 采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102521834A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104778728A (zh) * | 2015-05-04 | 2015-07-15 | 华北电力大学(保定) | 一种基于相位相关原理的天空图像云团运动速度计算方法 |
CN105550990A (zh) * | 2015-12-10 | 2016-05-04 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的傅里叶配准方法及装置 |
CN107516322A (zh) * | 2017-08-11 | 2017-12-26 | 浙江大学 | 一种基于对数极空间的图像物体大小和旋转估计计算方法 |
CN108876826A (zh) * | 2017-05-10 | 2018-11-23 | 深圳先进技术研究院 | 一种图像匹配方法及系统 |
CN109242907A (zh) * | 2018-09-29 | 2019-01-18 | 武汉光庭信息技术股份有限公司 | 一种基于照地高速摄像头的车辆定位方法及装置 |
CN109239680A (zh) * | 2018-10-12 | 2019-01-18 | 西安工业大学 | 一种低截获概率雷达lfm信号的参数估计方法 |
CN112686933A (zh) * | 2020-12-29 | 2021-04-20 | 中国科学院长春光学精密机械与物理研究所 | 基于改进互功率谱的星上图像配准叠加增强方法及系统 |
CN113240743A (zh) * | 2021-05-18 | 2021-08-10 | 浙江大学 | 基于神经网络的异构图像位姿估计及配准方法、装置及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6373970B1 (en) * | 1998-12-29 | 2002-04-16 | General Electric Company | Image registration using fourier phase matching |
CN102054272A (zh) * | 2010-11-03 | 2011-05-11 | 浙江大学 | 一种基于联合分数傅里叶变换的位移矢量探测方法及装置 |
-
2011
- 2011-12-12 CN CN2011104109122A patent/CN102521834A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6373970B1 (en) * | 1998-12-29 | 2002-04-16 | General Electric Company | Image registration using fourier phase matching |
CN102054272A (zh) * | 2010-11-03 | 2011-05-11 | 浙江大学 | 一种基于联合分数傅里叶变换的位移矢量探测方法及装置 |
Non-Patent Citations (1)
Title |
---|
彭静等: "一种基于极对数分数阶傅里叶变换的图像配准方法", 《第七届和谐人机环境联合学术会议HHME2011》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104778728A (zh) * | 2015-05-04 | 2015-07-15 | 华北电力大学(保定) | 一种基于相位相关原理的天空图像云团运动速度计算方法 |
CN105550990A (zh) * | 2015-12-10 | 2016-05-04 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的傅里叶配准方法及装置 |
CN105550990B (zh) * | 2015-12-10 | 2017-07-28 | 北京中科紫鑫科技有限责任公司 | 一种基于傅里叶配准的dna图像处理方法及装置 |
CN108876826B (zh) * | 2017-05-10 | 2021-09-21 | 深圳先进技术研究院 | 一种图像匹配方法及系统 |
CN108876826A (zh) * | 2017-05-10 | 2018-11-23 | 深圳先进技术研究院 | 一种图像匹配方法及系统 |
CN107516322A (zh) * | 2017-08-11 | 2017-12-26 | 浙江大学 | 一种基于对数极空间的图像物体大小和旋转估计计算方法 |
CN107516322B (zh) * | 2017-08-11 | 2020-08-18 | 浙江大学 | 一种基于对数极空间的图像物体大小和旋转估计计算方法 |
CN109242907A (zh) * | 2018-09-29 | 2019-01-18 | 武汉光庭信息技术股份有限公司 | 一种基于照地高速摄像头的车辆定位方法及装置 |
CN109239680A (zh) * | 2018-10-12 | 2019-01-18 | 西安工业大学 | 一种低截获概率雷达lfm信号的参数估计方法 |
CN112686933A (zh) * | 2020-12-29 | 2021-04-20 | 中国科学院长春光学精密机械与物理研究所 | 基于改进互功率谱的星上图像配准叠加增强方法及系统 |
CN112686933B (zh) * | 2020-12-29 | 2024-03-08 | 中国科学院长春光学精密机械与物理研究所 | 基于改进互功率谱的星上图像配准叠加增强方法及系统 |
CN113240743A (zh) * | 2021-05-18 | 2021-08-10 | 浙江大学 | 基于神经网络的异构图像位姿估计及配准方法、装置及介质 |
CN113240743B (zh) * | 2021-05-18 | 2022-03-25 | 浙江大学 | 基于神经网络的异构图像位姿估计及配准方法、装置及介质 |
WO2022241877A1 (zh) * | 2021-05-18 | 2022-11-24 | 浙江大学 | 基于神经网络的异构图像位姿估计及配准方法、装置及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102521834A (zh) | 采用极对数坐标表示的分数阶傅里叶变换的图像配准方法 | |
Ye et al. | Multi-point displacement monitoring of bridges using a vision-based approach | |
CN105654423B (zh) | 基于区域的遥感图像配准方法 | |
CN103941089A (zh) | 基于dft的正弦信号频率估计方法 | |
CN104730519B (zh) | 一种采用误差迭代补偿的高精度相位解缠方法 | |
CN104007318B (zh) | 获取信号时频函数的方法 | |
CN103456015A (zh) | 基于最优分数域Gabor谱特征的SAR目标检测方法 | |
CN104103047A (zh) | 一种心电图像倾斜度校正方法 | |
CN103791899B (zh) | 一种基于成像传感器视差的卫星姿态颤振探测方法 | |
CN103839262A (zh) | 一种基于直线和fft的sar图像配准方法 | |
CN102324045A (zh) | 基于Radon变换和极谐波变换的不变矩目标识别方法 | |
CN104616276A (zh) | 数字图像篡改盲检测系统 | |
CN105491315A (zh) | 一种投影仪伽马校正方法 | |
CN104049253B (zh) | 一种前向散射雷达阴影逆合成孔径成像方法 | |
CN104142495A (zh) | 一种基于频谱的斜视sar点目标插值与剖面截取方法 | |
CN104463849A (zh) | 一种基于颜色不变量与相位相关的彩色图像配准方法 | |
CN103530654A (zh) | 一种二维图形的对称轴的检测方法 | |
CN102609943A (zh) | 一种基于线性Radon变换算法的图像处理方法 | |
CN208155328U (zh) | 一种数据序列展示装置 | |
CN102175226A (zh) | 基于显著性特征的目标识别方法 | |
US20170146370A1 (en) | Method for determining an orthogonality error between two sensor signals | |
CN107422322A (zh) | 针对金属球目标的sar图像仿真方法 | |
Gonzalez | Fourier based registration of differentially scaled images | |
Qu et al. | A high-precision registration algorithm for heterologous image based on effective sub-graph extraction and feature points bidirectional matching | |
CN102003948B (zh) | 一种高精度光学条纹位相提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120627 |