基于傅里叶变换图像配准

最新推荐文章于 2023-10-03 17:47:05 发布

Adagrad

最新推荐文章于 2023-10-03 17:47:05 发布

阅读量3.3w

点赞数 17

本文链接：https://blog.csdn.net/u013049912/article/details/85986855

版权

1、理论

首先给出傅里叶变换的旋转、平移、比例放大特性，及如何利用这些特性进行图像配准。

1.1平移变化

如果图像 $f_{2}(x,y)$ 是图像 $f_{1}(x,y)$ 经平移 $(x_{0},y_{0})$ 后的图像，即 $f_{2}(x,y)=f_{1}(x-x_{0},y-y_{0})$ ，则对应的傅里叶变化F1和F2的关系为：

$F_{2}(\xi ,\eta )=e^{-j2\pi (\xi x_{0}+\eta y_{0})}F_{1}(\xi,\eta)$ (1)

且对应频域中两个图像的互能量谱为：

$\frac{F_{1}(\xi,\eta)F_{2}^*}{|F_{1}(\xi,\eta)F_{2}^*|}=e^{j2\pi (\xi x_{0}+\eta y_{0})}$ (2)

式中 $F_{2}^*$ 为 $F_{2}$ 的复共轭。平移理论表明，互能量谱的相位等于图像间的相位差。通过对互能量谱进行反变换，就可得到一个冲击函数 $\delta (x-x_{0},y-y_{0})$ 。此函数在偏僻位置处有明显的尖锐峰值。其他位置的值接近于零，所以据此就能找到两图像间的偏移量。

1.2 没有尺度变化的旋转特性

如果 $f_{2}(x,y)$ 是 $f_{1}(x,y)$ 经平移 $(x_{0},y_{0})$ 、旋转 $\theta _{0}$ 得到的图像，即：

$f_{2}(x,y)=f_{1}(xcos \theta _{0}+ysin\theta_{0}-x_{0},-xsin\theta_{0}+ycos\theta_{0}-y_{0})$ (3)

根据傅里叶的旋转和平移特性，变换后两图像间的关系为：

$F_{2}(\xi,\eta)=e^{-j2\pi (\xi x_{0}+\eta y_{0})}F_{1}(\xi cos\theta_{0}+\eta sin\theta_{0},-\xi sin\theta_{0}+\eta cos\theta_{0})$ (4)

假定M1，M2为F1和F2的能量，则：

$M_{2}(\xi,\eta)=M_{1}(\xi cos\theta_{0}+\eta sin\theta_{0},-\xi sin\theta_{0}+\eta cos\theta_{0})$ (5)

由公式（5）可以看出，F1和F2的能量是相同的，不过其中一个是另一个旋转后的副本。直角坐标中的旋转对应着极坐标角度的平移。因此将公式（5）进行极坐标描述为：

$M_{1}(\rho ,\theta)=M_{2}(\rho,\theta-\theta_{0})$ (6)

进而利用相位相关理论，可得到 $\theta_{0}$ 。

1.3 带比例放大的变换特性

如果f1为f2分别在水平和垂直方向上进行比例缩放后的图像，缩放因子为(a,b)，根据傅里叶尺度变换特性，有：

$F_{2}(\xi,\eta)=\frac{1}{ab}F_{1}(\xi/a,\eta/b)$ (7)

通过对数轴变换，比例变换可转换为平移变换（忽略乘积因子1/ab）：

$F_{2}(log\xi,log\eta)=F_{1}(log\xi -loga,log\eta-logb)$ (8)

上式可写成：

$F_{2}(x,y)=F_{1}(x-c,y-d)$ (9)

式中 $x=log\xi,y=log\eta,c=loga,d=logb$ 。则平移(c,d)可通过相位相关技术得到，尺度因子(a,b)由(c,d)得到：

$a=e^{c},b=e^{d}$ (10)

如果(x,y)尺度变化为(x/a,y/a)，则其极坐标可描述为：

$\rho_{1}=(x^2+y^2)^{1/2}\quad \theta_{1}=arctan(y/x)$ (11)

$\rho_{2}=((x/a)^2+(y/a)^2)^{1/2}=\rho_{1}/a$

$\theta_{2}=arctan((y/a)/(x/a))=\theta_{1}$ (12)

进而，如果f1为f2经平移、旋转和比例缩放后的图像，则它们的极坐标描述的能量谱间的关系为：

$M_{1}(\rho,\theta)=M_{2}(\rho/a,\theta-\theta_{0})$ (13)

$M_{1}(log\rho,\theta)=M_{2}(log\rho-loga,\theta-\theta_{0})$ (14)

$M_{1}(\xi,\theta)=M_{2}(\xi-d,\theta-\theta_{0})$ (15)

式中， $\xi=log\rho;d=loga$ 。利用上式和相位相关技术可得到比例因子a和旋转角 $\theta$ ，分别对要配准的图像进行比例变换和旋转后，再利用相位相关技术可求出图像间的偏移量。

2、总结

由上讨论可以看出，首先求出比例因子及旋转角，按此值对欲配准图像变换后，求出平移量，再进行变换可得到配准好的图像。具体步骤如下：

（1）对原始图像进行傅里叶变换，并求出各自的能量。

（2）高通滤波。

$H(\xi,\eta)=(1.0-X(\xi,\eta))(2.0-X(\xi,\eta))$ (16)

$X(\xi,\eta)=[cos(\pi \xi)cos(\pi \eta)]\quad-0.5\leqslant \xi,\eta\leqslant 0.5$ (17)

（3）将滤波后的各图像的能量转换成对数-极坐标形式，并求其互能量谱，得到比例系数及旋转角。

（4）将欲配准的图像经旋转、比例放大后再与参考图像一起计算互能量谱，从而得到平移量。

Adagrad

关注

17
点赞
踩
62

收藏

觉得还不错? 一键收藏
17
评论
基于傅里叶变换图像配准

1、理论首先给出傅里叶变换的旋转、平移、比例放大特性，及如何利用这些特性进行图像配准。1.1平移变化如果图像是图像经平移后的图像，即，则对应的傅里叶变化F1和F2的关系为：(1)且对应频域中两个图像的互能量谱为：(2)式中为的复共轭。平移理论表明，互能量谱的相位等于图像间的相位差。通过对互能量谱进行反变换，就可得到一个冲击函数。此函数在偏僻位置处有明显的尖锐峰值。其他位...
复制链接

扫一扫