删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

一种改进的基于奇异值分解的亚像素级图像配准算法

本站小编 Free考研考试/2021-12-25

凌程, 耿修瑞, 杨炜暾, 赵永超
中国科学院电子学研究所 中国科学院空间信息处理与应用系统技术重点实验室, 北京 100190;中国科学院大学, 北京 100049
2017年11月13日 收稿; 2018年1月9日 收修改稿
基金项目: 高分5号应用共性关键技术项目(30-Y20A28-9004-15/17)和国家重大科研仪器研制项目(41427805)资助
通信作者: 凌程, E-mail:lingcheng15@mails.ucas.ac.cn

摘要: 基于奇异值分解(SVD)的相位相关法是一种经典的具有亚像素级精度的图像配准算法,但当两幅待配准图的平移量较大或噪声较强时,该配准算法中的积分法得到的相位解缠结果往往不可靠。根据线性相位的单调变化特性,提出一种改进的相位解缠算法,通过比较相邻相位差和趋势斜率的一致性判断是否进行校正,从而得到真实相位值。针对真实光学图像的实验表明,该方法可以有效地对积分法所得到的结果进行校正,进而提高匹配精度。
关键词: 图像配准亚像素奇异值分解相位解缠积分法
An improved subpixel image registration algorithm based on singular value decomposition
LING Cheng, GENG Xiurui, YANG Weitun, ZHAO Yongchao
Key Laboratory of Technology in Geo-spatial Information Processing and Application System of CAS, Institute of Electronics, Chinese Academy of Science, Beijing 100190, China; University of Chinese Academy of Science, Beijing 100049, China


Abstract: Phase correlation method based on singular value decomposition (SVD) is a known subpixel image registration algorithm. However, when the translation between two images to be registered is large or there is high noise, the phase unwrapping results obtained using the integral method in the registration algorithm are often unreliable. In this work, an improved phase unwrapping algorithm is proposed based on the monotonic variation of linear phase. By comparing the consistency of the adjacent phase differences and the trend slopes, the algorithm determines whether to make corrections, so as to get the real phase. Experiments on real optical images show that the proposed method effectively corrects the results obtained using the integral method, and hence increases matching accuracy.
Keywords: image registrationsubpixelsingular value decompositionphase unwrappingintegral method
图像配准是一项非常重要的图像预处理技术。随着数字图像处理的发展,图像配准技术已经被广泛应用于高分辨率图像生成[1]、医学图像对齐[2]、高精度3D重建[3]、遥感图像变化检测、多源图像融合等领域。这些关键领域中,需要高精度的图像配准,即亚像素级的图像配准。现有的亚像素级图像配准算法主要分为3类[4-5]:基于插值的方法[6-8]、解最优化问题法[9-10]和扩展的相位相关法[2, 11-13]。其中,基于插值的方法采用插值技术对相似性函数或者图像进行插值重采样,以插值后的最大值位置作为配准位置,以此达到亚像素的精度。Guizar-Sicairos等[7]对相似形函数内插法进行快速改进;Yousef等[8]在此基础上进一步提高速度,在保证配准精度的同时大大减少计算复杂度。本文在此将Yousef等[8]提出的快速相似性函数内插法称为fast-SSDFT方法。解最优化问题法的目标函数为衡量两幅待配准图像之间相似程度的函数,支持各种变换模型,具有参数选择灵活,计算精度高的特点,但是其计算量非常大,收敛概率和得到全局最优解的概率也需要提高。基于相位相关的方法利用傅里叶变换的特性,将时域上的平移量转化为求解频域上的相位相关,具有精度高、对噪声鲁棒和计算复杂度低等优点。
目前所有扩展的相位相关法都是基于傅里叶变换的平移特性,即两幅待配准的图像在空间域上存在平移,转换到频域上则表现为线性相位,通过对频域上线性相位的求解来获得空间域上的平移量。传统的相位相关法直接将归一化互功率谱矩阵进行傅里叶反变换,通过判断冲激函数的峰值位置得到平移参数的值,但是这种方法只能得到整数像素值。Hoge[11]提出,通过对归一化互功率谱矩阵进行奇异值分解,用最大奇异值所对应的奇异向量估算出平移量,可以达到亚像素的精度。该方法通过选取最大奇异值,抛弃小的奇异值,起到对噪声过滤的作用,因此该方法对平移参数的估计具有较高的精度。另外为了求得亚像素精度的平移量,需要对得到的奇异向量的相位值进行一维相位解缠,然后用最小二乘法估计出平移量。Tong等[12]利用Ransac算法[14]替代最小二乘法对解缠后的相位值进行估计,在噪声较大的情况下,Ransac算法能够滤除噪声带来的异常值,使得估计出的平移参数精度更高。Feng等[13]针对低信噪比的图像提出一种鲁棒的亚像素配准算法,该算法通过计算两幅图像的互相关系数[15]的质心得到精确的平移参数,从而达到亚像素的配准精度,本文在此将其称为CC_centriod方法。
在理想无噪的情况下,缠绕相位具有一定的周期性,可以通过积分的方法进行相位解缠。通过积分的方法求得真实相位具有一个前提,即相邻两个像元间的真实相位差的绝对值要小于π[16],但是如果两幅图像之间的平移量超过某个阈值或者受噪声影响时,所得到的一维缠绕相位向量相邻的真实相位差的绝对值会超出π,因此传统的积分法进行相位解缠就会变得不再适用。本文提出一种改进的相位解缠方法对传统积分法得到的结果进行校正,改进的算法可以避免积分法受到此前提和噪声的影响,使得估计出的平移参数更精确。本文把基于奇异值分解(SVD)的相位相关法结合改进的相位解缠方法称为SVD-IIM(SVD-improved integral method)。
1 基于SVD分解的相位相关法1.1 方法原理和流程令x0y0表示两幅大小为M×N的图像之间在行方向和列方向上的平移量,可以表示成
$g\left( {x,y} \right) = f\left( {x - {x_0},y - {y_0}} \right).$ (1)
式中:f(x, y)和g(x, y)为二维离散函数,其中变量x=1, 2, 3, …, M和变量y=1, 2, 3, …, N。根据傅里叶变换的性质可以得到
$G\left( {u,v} \right) = F\left( {u,v} \right)\exp \left\{ { - {\rm{j}}2{\rm{ \mathsf{ π} }}\left( {u{x_0}/M + v{y_0}/N} \right)} \right\}.$ (2)
式中:F(u, v)和G(u, v)分别是f(x, y)和g(x, y)的二维离散傅里叶变换函数,其中变量u=1, 2, 3, …, M和变量v=1, 2, 3, …, N。因此可以得到归一化互功率谱函数Q(u, v)
$\begin{array}{l}Q\left( {u,v} \right) = \frac{{G\left( {u,v} \right)F{{\left( {u,v} \right)}^ * }}}{{\left| {G\left( {u,v} \right)F{{\left( {u,v} \right)}^ * }} \right|}}\\\;\;\;\;\;\;\;\;\;\;\; = \exp \left\{ { - {\rm{j}}2{\rm{ \mathsf{ π} }}\left( {u{x_0}/M + v{y_0}/N} \right)} \right\}.\end{array}$ (3)
式中:*代表共轭,传统的相位相关法直接计算归一化互功率谱函数Q(u, v)的傅里叶反变换:
$Q\left( {x,y} \right) = \delta \left( {x - {x_0},y - {y_0}} \right).$ (4)
在理想情况下,归一化互功率谱矩阵Q在理论上是一个秩为1的矩阵,因此归一化互功率谱函数Q(u, v)可以进行如下分解:
$\begin{array}{l}Q\left( {u,v} \right) = \exp \left\{ { - {\rm{j}}2{\rm{ \mathsf{ π} }}\left( {u{x_0}/M + v{y_0}/N} \right)} \right\}\\\;\;\;\;\;\;\;\;\;\;\; = \exp \left\{ { - {\rm{j}}2{\rm{ \mathsf{ π} }}u{x_0}/M} \right\}\exp \left\{ { - {\rm{j}}2{\rm{ \mathsf{ π} }}v{y_0}/N} \right\}\\\;\;\;\;\;\;\;\;\;\;\; = {\mathit{\boldsymbol{q}}_{{x_0}}}\left( u \right)\mathit{\boldsymbol{q}}_{{y_0}}^{\rm{H}}\left( v \right).\end{array}$ (5)
式中:H代表共轭转置,令向量qx0qy0为:
${\mathit{\boldsymbol{q}}_{{x_0}}} = {\left( {{{\rm{e}}^{ - {\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/M} \right){x_0}}},{{\rm{e}}^{ - {\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/M} \right)2{x_0}}}, \cdots ,{{\rm{e}}^{ - {\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/M} \right)M{x_0}}}} \right)^{\rm{T}}},$ (6)
${\mathit{\boldsymbol{q}}_{{y_0}}} = {\left( {{{\rm{e}}^{{\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/N} \right){y_0}}},{{\rm{e}}^{{\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/N} \right)2{y_0}}}, \cdots ,{{\rm{e}}^{{\rm{j}}\left( {2{\rm{ \mathsf{ π} }}/N} \right)N{y_0}}}} \right)^{\rm{T}}}.$ (7)
因此归一化互功率谱矩阵Q可以表示为
$\mathit{\boldsymbol{Q}} = {\mathit{\boldsymbol{q}}_{{x_0}}} * \mathit{\boldsymbol{q}}_{{y_0}}^{\rm{H}}.$ (8)
但是,在一般情况下,Q并不是一个秩为1的矩阵,因此文献[11]提出可以通过对Q矩阵进行SVD分解,用最大奇异值和相对应的奇异向量对Q进行秩1估计,由于奇异值最大的分量受噪声影响较小,因而秩1估计可以起到抑制噪声的作用。对Q进行秩1估计:
$\mathit{\boldsymbol{Q}} \approx {\mathit{\boldsymbol{u}}_1} * {\sigma _1} * \mathit{\boldsymbol{v}}_1^{\rm{H}}.$ (9)
式中:σ1为最大的奇异值,u1v1σ1所对应的奇异向量,因此有u1qx0v1qy0。对u1v1进行相位解缠,得到解缠后的相位pxpy,并用最小二乘法拟合pxpy,得到pxpy的斜率kxky,通过得到的斜率最终求得平移参数:
${x_0} = - \frac{{{k_x}M}}{{2{\rm{ \mathsf{ π} }}}},{y_0} = - \frac{{{k_x}N}}{{2{\rm{ \mathsf{ π} }}}}.$ (10)
图 1所示的流程图,为减弱边缘效应的影响,对待配准的图像f(x, y)和g(x, y)分别进行Black窗处理后,再计算其相位相关矩阵Q。为进一步减弱频谱混叠和边缘效应的影响,需要滤除Q矩阵中的高频分量和幅度值小于特定阈值的频率分量。以直流分量为中心,以0.3N(N代表图像短边长度)为半径,滤除距离直流分量中心在此半径之外的高频分量;计算直流分量周围5×5邻域的幅度平均值,并乘以因子0.01~0.05,作为阈值,滤除小于此阈值的频率分量[12]。然后对Q进行SVD分解并按照式(9)对其进行秩1估计,得到最大奇异值所对应的奇异向量u1v1,用积分法对u1v1的相位值进行相位解缠,并用最小二乘法估计出解缠后相位的斜率kxky,最后利用式(10)计算得到平移参数x0y0
Fig. 1
Download: JPG
larger image
图 1 基于SVD分解的相位相关法流程图 Fig. 1 Flow chart of phase correlation method based on SVD
图 1 基于SVD分解的相位相关法流程图

Fig. 1 Flow chart of phase correlation method based on SVD -->

1.2 积分法相位解缠原理在理想情况下,图像的采样率满足Nyquist采样定理,采样频率必须大于信号最高频率的2倍,解缠相位中相邻像素点之间的相位差值不可能超过半个周期,当满足此条件时必然能通过积分法,由缠绕相位解缠出正确的真实相位值。令?(m)为周期缠绕前的真实相位值,$φ$(m)为相应的缠绕相位,相位缠绕算子可以如下定义:
$\omega \left( {\phi \left( m \right)} \right) = \varphi \left( m \right) = \phi \left( m \right) = 2{\rm{ \mathsf{ π} }}k\left( m \right).$ (11)
其中-π < ω(?(m)) < π,得到属于(-π, π)的缠绕相位$φ$(m),定义差分算子Δ,因此对于真实相位值有
$\Delta \phi \left( m \right) = \phi \left( {m + 1} \right) - \phi \left( m \right).$ (12)
由于真实相位值中相邻像素点之间的相位差值不可能超过半个周期,所以有
$ - {\rm{ \mathsf{ π} }} < \Delta \phi \left( m \right) < {\rm{ \mathsf{ π} }}.$ (13)
对相邻缠绕相位也进行差分运算得
$\begin{array}{l}\Delta \varphi \left( m \right) = \varphi \left( {m + 1} \right) - \varphi \left( m \right) = \Delta \phi \left( m \right) + \\\;\;\;\;\;\;\;\;\;\;\;\;\;\;2{\rm{ \mathsf{ π} }}\Delta k\left( m \right).\end{array}$ (14)
对该相位差使用缠绕算子得
${\rm{ \mathsf{ π} }}\left( {\Delta \varphi \left( m \right)} \right) = \Delta \varphi \left( m \right) + 2{\rm{ \mathsf{ π} }}\Delta k\left( m \right) + 2{\rm{ \mathsf{ π} }}k'\left( m \right).$ (15)
根据缠绕算子ω的定义,ω$φ$(m))必然在区间(-π, π)内,由式(13)可知,Δ?(m)也属于区间(-π, π),因此有
$\Delta k\left( m \right) + k'\left( m \right) = 0.$ (16)
所以可将式(15)简化为
$\omega \left( {\Delta \varphi \left( m \right)} \right) = \Delta \phi \left( m \right).$ (17)
因此最终的真实相位可以表示成
$\phi \left( m \right) = \phi \left( 0 \right) + \sum\limits_{n = 0}^{m - 1} {\omega \left( {\Delta \varphi \left( n \right)} \right)} .$ (18)
从以上推导可以看出,使用积分法进行相位解缠依赖于一个重要的前提,即真实相位中相邻像素点之间的相位差值必须在(-π, π)范围内。而在实际遥感图像配准应用中,由于有地形突变、相位混叠、平移量较大等各种因素的存在,相位信息中将会引入噪声,在计算解缠相位的过程中,并不确定式(13)是否成立,所以通过积分法即式(18)计算所得到的相位值只是相位的估计值,其准确性并不能保证[17]
1.3 改进的相位解缠方法在使用积分法进行相位解缠时,式(18)中的缠绕算子ω会将缠绕相位差折叠到区间(-π, π)里,当真实相位差的绝对值大于π时,则会出现2kπ(k≠0)的误差(此误差在大多数情况下为2π),因此本文提出一种改进算法,能够针对这个误差进行校正。以估计平移参数x0为例,首先利用积分法对u1进行相位解缠得到解缠相位的初始结果px,根据式(6)可知,在理想状态下px是线性递增或递减的直线,如果px不满足这样一个递增或递减的趋势,那么本文就认为有必要对其进行校正处理。利用最小二乘法估计初始结果px的斜率,得到初始斜率kx1,虽然初始斜率kx1和真实斜率存在一定误差,但是已经足够描述递增或递减的趋势,因此令kx1为趋势斜率,对px进行校正。如果px中相邻相位之差通过加上或者减去2π的误差校正之后更符合趋势斜率,那么执行此校正操作,否则不进行校正。完成一次px的校正之后,需要用最小二乘法重新估计px的斜率,得到新的趋势斜率,直到前后趋势斜率不发生变化时,完成校正。最后按照式(10)利用稳定的趋势斜率精确计算出平移参数x0
图 2所示,如果当真实相位差的绝对值大于π时,使用积分法进行相位解缠会出现图 2(a)所示的现象,如果对?(i+1)进行加上或减去2π的校正操作后更符合趋势斜率kx,那么执行此校正操作,使得校正后的?(i+1)-?(i)更符合趋势斜率,如图 2(b)所示。
Fig. 2
Download: JPG
larger image
图 2 相位校正示意图 Fig. 2 Schematic of phase correlation
图 2 相位校正示意图

Fig. 2 Schematic of phase correlation -->

2 实验结果及分析本文选取的实验数据为高分2号(GF-2)卫星和资源3号(ZY-3)卫星的高分辨率多光谱图像数据(部分图像数据如图 3所示),分别进行校正效果实验和配准精度对比实验。其中校正效果实验主要对比传统积分法和改进的相位解缠方法对缠绕相位的解缠效果以及相应的配准精度,配准精度对比则实验对比在不同的平移量、加入不同大小的高斯白噪声、不同混叠情况下,本文所提SVD-IIM方法与已有的SVD方法[11]、SVD-RANSAC方法[12],CC-centriod[13]方法和fast-SSDFT[8]方法这4种方法的配准精度。
Fig. 3
Download: JPG
larger image
图 3 GF-2数据(a,b)和ZY-3数据(c,d) Fig. 3 Imagery of GF-2(a, b) and ZY-3(c, d)
图 3 GF-2数据(a,b)和ZY-3数据(c,d)

Fig. 3 Imagery of GF-2(a, b) and ZY-3(c, d) -->

2.1 校正效果实验对上文所述的高分2号数据在光谱维上取均值得到其灰度图像,在其中选取一幅大小为1 280像素×1 280像素的原始参考图像,如图 4(a)所示。将其在原灰度图像上分别在行方向平移541个像素和在列方向平移548个像素,得到另一幅平移后的原始待配准图像,如图 4(b)所示。那么这两幅灰度图像的平移关系为(541, 548),分别对两幅图像进行10倍降采样得到两幅大小为128像素×128像素的图像,如图 4(c)4(d)所示,因此其平移关系变为(54.1, 54.8)[18]
Fig. 4
Download: JPG
larger image
图 4 原始配准图和10倍降采样后的配准图 Fig. 4 Original images and images after 10 times downsampling used for registration
图 4 原始配准图和10倍降采样后的配准图

Fig. 4 Original images and images after 10 times downsampling used for registration -->

按照图 1所示的流程计算图 4(c)图 4(d)的归一化互功率谱矩阵Q(u, v),并滤除其中的高频分量和小于特定阈值的频率分量,再对Q矩阵进行SVD分解得到最大奇异值对应的奇异向量uv,分别利用传统积分法和改进的相位解缠方法对uv的相位进行相位解缠。利用最小二乘法分别对图 5(b)5(d)所示的解缠相位值的斜率进行估计,其配准结果如表 1所示,显然传统积分法的误差远远大于改进的相位解缠方法的误差。由图 5(b)5(d)可以看出,改进的相位解缠算法能对传统积分法得到的结果进行很好的校正,因此当参考图像和待配准图像之间平移量较大时,即当相邻像素点之间相位差的绝对值接近或超过p时,利用传统的积分法相位解缠得到的结果并不可靠,而经过改进的相位解缠方法可以精确地估计出平移参数。
Fig. 5
Download: JPG
larger image
图 5 相位解缠结果示意图 Fig. 5 Schematic of the phase unwrapping results
图 5 相位解缠结果示意图

Fig. 5 Schematic of the phase unwrapping results -->


Table 1
表 1 传统积分法和改进相位解缠方法配准精度比较Table 1 Precision comparison between the traditional integral method and the improved phase unwrapping method
改进的相位解缠方法 传统积分法
真实平移量 (54.1, 54.8) (54.1, 54.8)
估计值 (54.088 5, 54.961 2) (25.071 1, 20.880 1)
绝对误差 (0.011 5, 0.161 2) (29.028 9, 33.919 9)

表 1 传统积分法和改进相位解缠方法配准精度比较Table 1 Precision comparison between the traditional integral method and the improved phase unwrapping method

2.2 配准精度对比实验对上文所述的高分2号数据和资源3号数据在光谱维上取均值,利用得到的灰度图像进行配准精度对比试验,选取大小为1 280像素×1 280像素的子图为原始参考图像,以其为基准进行不同程度的平移,得到原始待配准图像。分别对原始参考图像和原始待配准图像进行10倍降采样得到实验所用的大小为128像素×128像素的配准图像。图 6图 7为本文提出的SVD-IIM方法与其他4种方法分别在高分2号数据和资源3号数据上进行的对比结果。其中图 6(a)6(d)图 7(a)7(d)为对于不同平移量的绝对误差对比示意图,其中图 6(a)图 7(a)平移量为57~62个像素,图 6(d)图 7(d)平移量为55~60个像素,步长均为0.1个像素;图 6(b)6(e)图 7(b)7(e)为在加入不同高斯白噪声情况下的绝对误差对比示意图,将配准图像归一化至0~256,加入标准差为6~10的高斯白噪声,真实平移量为(59.4, 7.8);对待配准图减采样之前的原图进行不同程度的高斯模糊可以模拟不同情况下的混叠,而高斯模糊等价于低通滤波器[2]图 6(c)6(f)图 7(c)7(f)为在不同混叠程度下的绝对误差对比示意图,高斯模糊窗大小为25×25,模糊因子为高斯函数的标准差, 真实平移量为(59.4, 7.8)。图中所有绝对误差为100幅随机选取的待配准图像配准误差的平均值。
Fig. 6
Download: JPG
larger image
图 6 GF-2数据配准精度对比图 Fig. 6 The contrast diagram of GF-2 data registration accuracy
图 6 GF-2数据配准精度对比图

Fig. 6 The contrast diagram of GF-2 data registration accuracy -->


Fig. 7
Download: JPG
larger image
图 7 ZY-3数据配准精度对比图 Fig. 7 The contrast diagram of ZY-3 data registration accuracy
图 7 ZY-3数据配准精度对比图

Fig. 7 The contrast diagram of ZY-3 data registration accuracy -->

图 6(a)~6(c)图 7(a)~7(c)所示,对于不同的平移量、加入不同大小的高斯白噪声、不同混叠情况下,本文提出的SVD-IIM方法的误差都在0附近,而原有的SVD方法和SVD-RANSAC方法的误差都在10个像素以上,因此可以进一步确定在平移量较大时传统积分法得到的相位解缠结果是错误的,而改进的相位解缠方法所得到的结果是可靠的,并且精度较高。
图 6(d)~6(f)图 7(d)~7(f)对比在不同平移量,加入不同高斯白噪声,不同混叠情况下的3种方法的配准精度,可以看出本文提出的SVD-IIM方法的误差相比于CC-centriod方法和fast-SSDFT方法的误差更小,并且更稳定。对于较大平移量的配准问题,参考图像和待配准图像间存在较大的非匹配区域,这些非匹配区域的信息在估计平移参数时如同噪声,因此在存在大量噪声的情况下,本文提出的SVD-IIM方法可以将配准精度控制在0.1个像素范围内。
3 结束语本文针对基于奇异值分解(SVD)的相位相关法中传统积分法进行相位解缠的可靠性问题,利用线性相位单调变化的特性,提出一种改进的相位解缠方法,对传统积分法得到的结果进行校正,可以精确得到真实相位值,使得对于两幅平移量较大的待配准图像可以稳定并准确地估计出其平移参数。本文用真实光学图像的实验验证了本文提出的SVD-IIM方法在配准问题上相比于原有的SVD方法、SVD-RANSAC方法、CC-centriod方法和fast-SSDFT方法更具鲁棒性和准确性,具有重要的应用价值。
参考文献
[1] Vandewalle P, Süsstrunk S, Vetterli M. A frequency domain approach to registration of aliased images with application to supper-resolution[J]. Eurasip Journal on Advances in Signal Processing, 2006(1): 1-14.
[2] Stone H S, Orchard M T, Chang E, et al. A fast direct Fourier-based algorithm for subpixel registration of images[J]. IEEE Trans on Geoscience and Remote Sensing, 2001, 39(10): 2235-2243. DOI:10.1109/36.957286
[3] Muquit M A, Shibahara T. A high-accuracy passive 3D measurement system using phase-based image matching[J]. Ieice Transactions on Fundamentals of Electronics Communications and Computer Sciences, 2006, 89(3): 686-697.
[4] 黎俊, 彭启民, 范植华. 亚像素级图像配准算法研究[J]. 中国图像图像学报, 2008, 13(11): 2070-2075.
[5] 卢浩, 刘团结, 尤红建. 亚像素级的图像配准方法[J]. 国外电子测量技术, 2012, 31(4): 45-49. DOI:10.3969/j.issn.1002-8978.2012.04.013
[6] Tian Q, Huhns M N. Algorithms for subpixel registration[J]. Computer Vision, Graphics, and Image Processing, 1986, 35(2): 220-233. DOI:10.1016/0734-189X(86)90028-9
[7] Guizar-Sicairos M, Thurman S T, Fienup J R. Efficient subpixel image registration algorithms[J]. Optics Letters, 2008, 33(2): 156-158.
[8] Yousef A, Li J, Karim M. High-speed image registration algorithm with subpixel accuracy[J]. IEEE Signal Processing Letters, 2015, 22(10): 1796-1800. DOI:10.1109/LSP.2015.2437881
[9] Thevenaz P, Unser M. Optimization of mutual information for multiresolution image registration[J]. IEEE Transaction on Image Processing, 2000, 9(12): 2083-2099. DOI:10.1109/83.887976
[10] 徐宝昌, 陈哲, 赵龙. 一种改进的最小二乘景象匹配算法[J]. 北京航空航天大学学报, 2005, 31(8): 848-852. DOI:10.3969/j.issn.1001-5965.2005.08.005
[11] Hoge W S. A subspace identification extension to the phase correlation method[J]. IEEE Trans on Medical Imaging, 2003, 22(2): 277-280.
[12] Tong X, Ye Z, Xu Y, et al. A noval subpixel phase correlation method using singular value decomposition and unified random sample consensus[J]. IEEE Trans on Geoscience and Remote Sensing, 2015, 63(3): 207-224.
[13] Feng S, Deng L, Shu G, et al. A subpixel registration algorithm for low PNSR images[C]//2012 IEEE fifth International Conference on Advanced Computational Intelligence, Nanjing, 2012: 626-630.
[14] Fischler M A, Bolles R C. Random sample consensus:a paradigm for model fitting with applications to image analysis and automated cartography[M]. United States of America: ACM, 1981: 726-740.
[15] Zitová B, Flusser J. Image registration methods:a survey[J]. Image and Vision Computing, 2003, 21(11): 977-1000. DOI:10.1016/S0262-8856(03)00137-9
[16] Itoh K. Analysis of the phase unwrapping algorithm[J]. Applied Optics, 1982, 21(14): 2470. DOI:10.1364/AO.21.002470
[17] 邓晓龙. InSAR相位解缠算法研究[D].哈尔滨: 哈尔滨工业大学, 2013.
[18] Foroosh H, Zerubia J B, Berthod M. Extension of phase correlation to subpixel registration[J]. IEEE Transactions on Image Processing, 2002, 11(3): 188-200.


相关话题/图像 数据 计算 实验 相位

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 建筑物的高分辨率SAR图像仿真方法
    任苗苗1,2,潘卓1,徐向辉1,苏海涛31.中国科学院电子学研究所,北京100190;2.中国科学院大学,北京100049;3.中国人民解放军62026部队,西安7100322017年8月12日收稿;2018年1月24日收修改稿基金项目:国家自然科学基金(61471340)和中国科学院电子学研究所" ...
    本站小编 Free考研考试 2021-12-25
  • 基于车载激光扫描数据的城区道路自动提取
    赵海鹏1,2,习晓环1,王成1,雷钊1,21.中国科学院遥感与数字地球研究所中国科学院数字地球重点实验室,北京100094;2.中国科学院大学资源与环境学院,北京1000492017年7月24日收稿;2017年9月30日收修改稿基金项目:国家科技部重大科学仪器设备开发专项(2013YQ120343) ...
    本站小编 Free考研考试 2021-12-25
  • 沉积环境下涡轮叶片前缘气膜冷却的实验研究*
    在燃气涡轮发动机运行的过程中,涡轮部件会受到污染物的沉积、腐蚀和侵蚀。当飞机在起降或低空飞行时,燃气涡轮发动机往往会吸入灰尘、砂砾或其他颗粒物,特别是在沙漠或火山活跃地区。此外,航空燃料在高温环境下生成的金属化合物的副产物,也是极易沉积的颗粒物。其中,钒燃烧的产物V2O3会在涡轮叶片上沉积,在过剩热 ...
    本站小编 Free考研考试 2021-12-25
  • 分离式飞机应急数据记录跟踪系统设计与试验*
    在法航AF447和马航MH370等海上空难中,由于飞行数据记录器(FlightDataRecorder,FDR),也称“黑匣子”,固定在机身中,空难发生时记录器随飞机残骸沉入海底,给打捞和定位造成困难,甚至无法打捞,严重制约了及时救援和事故调查。特别是马航MH370事故,数国历时几年竭尽全力的救援和 ...
    本站小编 Free考研考试 2021-12-25
  • 气冷涡轮导叶流热耦合计算及机理*
    涡轮内部的流场与温度场是高度耦合的,一方面涡轮内部流场存在着如激波、二次涡系等复杂流动结构,会影响到部件温度场分布及其换热特性,另一方面,固体域的温度分布及冷却射流与主流的相互作用也将反过来影响流场的流动特征[1]。因此,理解和掌握涡轮内部流场温度场相互作用机制,对于提高冷却效率和准确预测热负荷有着 ...
    本站小编 Free考研考试 2021-12-25
  • 轴向通流旋转盘腔流动换热稳态实验与数值模拟*
    当前,提高涡轮前温度是提升航空发动机综合性能的主要途径之一。为使涡轮等热端部件能够在愈发恶劣的环境下正常、持久地工作,需要从压气机级间抽取冷气对其进行冷却。冷却气体用量的增加可以提高冷却效果,但提高了主流损失,对整机性能影响巨大;冷气由压气机级间进入盘腔后会与压气机转盘和盘罩换热,形成复杂的流动结构 ...
    本站小编 Free考研考试 2021-12-25
  • 基于残缺图像样本的生成对抗网络图像修复方法*
    图像修复是将图像或视频中丢失或损坏部分重建的过程。常见的图像修复方法主要有2类:传统的修复方法[1-10]和基于神经网络或深度神经网络学习的修复方法。传统的修复方法根据图像的边缘、纹理和结构,利用图像中的残留信息推断被破坏的内容,完成图像修复。例如,[1-2]提出了流体动力学修复方法,使用边界梯度匹 ...
    本站小编 Free考研考试 2021-12-25
  • 无人直升机系留气动载荷CFD计算分析*
    无人直升机停放在地面和舰船甲板等环境中时有可能被风吹动或倾覆,因此在地面或舰船甲板等停放时需要进行系留。在选取系留点和设计系留绳索时需先得到大风条件下无人直升机的系留载荷[1-5]。系留载荷的计算需要考虑无人直升机的自身重力、轮胎与地面的摩擦力、系留绳索的预紧力和直升机受到的气动载荷等。其中,系留气 ...
    本站小编 Free考研考试 2021-12-25
  • 优化Landweber迭代快速电磁层析成像图像重建算法*
    电学层析成像(ElectricalTomography,ET)技术作为一种高新检测技术,包括电阻层析成像(ElectricalResistanceTomography,ERT)[1-2]、电磁层析成像(ElectromagneticTomography,EMT)[3-4]、电容层析成像(Electr ...
    本站小编 Free考研考试 2021-12-25
  • 基于改进加权响应面的结构可靠度计算方法*
    在现有的结构可靠度分析方法中,一次二阶矩法[1]、二次二阶矩法[2-3]的精度较低,并且在非线性程度较高的情况下还会遇到无法收敛的问题。蒙特卡罗法[4-5]虽然能够得到精确解,但需要大量的抽样和计算时间,限制了其实际应用。响应面法[6]采用多项式函数来近似极限状态函数,原理简单、易于操作且计算效率较 ...
    本站小编 Free考研考试 2021-12-25