1、问题背景
此问题来源于冈萨雷斯的数字图像处理(第三版)5.19题,自己做的时候,没发现哪里有问题,一个同学的提问使我重新审视了这个问题,并感觉对傅里叶变换有了新的理解。
根据下面一部分的基础知识,我们可以计算出线性运动对应的模糊滤波器 H ( u , v ) H(u,v) H(u,v),但是对于图像的旋转,这种非线性的图像变换关系,并没有给出直接的计算方法分理出 H ( u , v ) H(u,v) H(u,v)。比较自然地,我们考虑使用极坐标变换 ( r , θ ) (r,\theta) (r,θ),使得图像的旋转随着 θ \theta θ线性变化,进而求出恢复后的图像。
2、基础知识
在图像获取的过程中,图像与传感器之间的线性运动会使最终得到的图像出现模糊现象。对于这种线性模糊,往往可以通过积分来处理:假设图像 f ( x , y ) f(x,y) f(x,y)进行平面运动, x 0 ( t ) x_0(t) x0(t)和 y 0 ( t ) y_0(t) y0(t)分别是在 x x x和 y y y方向上随时间变化的分量,那么在记录介质任意点的曝光总数是通过对时间间隔内瞬时曝光量的积分得到的,设 T T T为曝光时间,则有
g ( x , y ) = ∫ 0 T f [ x − x 0 ( t ) , y − y 0 ( t ) ] d t (2.1) g(x, y)=\int_{0}^{T} f\left[x-x_{0}(t), y-y_{0}(t)\right] \mathrm{d} t\tag{2.1} g(x,y)=∫0Tf[x−x0(t),y−y0(t)]dt(2.1)
其中, g ( x , y ) g(x,y) g(x,y)为模糊后的图像。
根据傅里叶变换的定义,可得
G ( u , v ) = ∫ − ∞ ∞ ∫ − ∞ ∞ g ( x , y ) e − j 2 π ( u x + v y ) d x d y = ∫ − ∞ ∞ ∫ − ∞ ∞ [ ∫ 0 T f [ x − x 0 ( t ) , y − y 0 ( t ) ] d t ] e − j 2 π ( u x + v y ) d x d y (2.2) \begin{aligned} G(u, v) &=\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} g(x, y) \mathrm{e}^{-\mathrm{j} 2 \pi(u x+v y)} \mathrm{d} x \mathrm{d} y \\ &=\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}\left[\int_{0}^{T} f\left[x-x_{0}(t), y-y_{0}(t)\right] \mathrm{d} t\right] \mathrm{e}^{-\mathrm{j} 2 \pi(u x+v y)} \mathrm{d} x \mathrm{d} y \end{aligned}\tag{2.2} G(u,v)=∫−∞∞∫−∞∞g(x,y)e−j2π(ux+vy)dxdy=∫−∞∞∫−∞∞[∫0Tf[x−x0(t),y−y0(t)]dt]e−j2π(ux+vy)dxdy(2.2)
改变积分顺序,可得
G ( u , v ) = ∫ 0 T [ ∫ − ∞ ∞ ∫ − ∞ ∞ f [ x − x 0 ( t ) , y − y 0 ( t ) ] e − j 2 π ( u x + v y ) d x d y ] d t (2.3) G(u, v)=\int_{0}^{T}\left[\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} f\left[x-x_{0}(t), y-y_{0}(t)\right] \mathrm{e}^{-\mathrm{j} 2 \pi(u x+v y)} \mathrm{d} x \mathrm{d} y\right] \mathrm{d} t\tag{2.3} G(u,v)=∫0T[∫−∞∞∫−∞∞f[x