数字图像处理第四章——频率与滤波

本文链接：https://blog.csdn.net/qq_67725022/article/details/140053842

一、背景

傅里叶级数：任何周期函数都可以表示为不同频率的正弦和/或余弦之和的形式，每个正弦项和/或余弦项乘以不同的系数。

傅里叶变换：非周期函数(但该曲线下的面积是有限的)也可以用正弦和/或余弦乘以加权函数的积分来表示。

二、基本概念

1、复数

复数 $C$ ：

$C=R+jI$

( R和Ⅰ是实数, $j=\sqrt{-1}$ )。

复数 $C$ 的共轭 $C^{*}$ ：

$C^{*}=R-jI$ 。

极坐标下：

$C=\left | C \right |\left ( \cos \theta +j \sin \theta \right )$ ( $\left | C \right |=\sqrt{R^{2}+I^{2}}$ ).

欧拉公式：

$e^{j\theta }=\cos\theta +j\sin\theta$ 。

2、傅里叶级数

$f(t)=\sum_{n=-\infty }^{\infty}c_{n}e^{j\frac{2\pi n}{T}t}$

其中： $c_{n}=\frac{1}{T}\int_{-T/2}^{T/2}f(t)e^{-j\frac{2\pi n}{T}t}dt,n=0,\pm 1,\pm 2,...$ 是系数

3、冲激及其取样特性

连续变量t在t =0处的单位冲激表示为 $\delta (t)$ ，

$\delta (t)=\left\{\begin{matrix} \infty , &t=0 \\0, &t\neq 0 \end{matrix}\right.$

它还被限制为满足等式

$\int_{-\infty }^{\infty }\delta (t)dt=1$

物理上，如果我们把t解释为时间，那么一个冲激可看成是幅度无限、持续时间为0、具有单位面积的尖峰信号。一个冲激具有关于如下积分的所谓取样( sifting)特性:

$\int_{-\infty }^{\infty }f(t)\delta (t)dt=f(0)$

的取样特性的一种更为一般的说明涉及位于任意点t的冲激，表示为 $\delta (t-t_{0})$ 。在这种情况下，取样特性变为

$\int_{-\infty }^{\infty }f(t)\delta (t-t_{0})dt=f(t_{0})$

令x表示一个离散变量。单位离散冲激8(x)在离散系统中的作用与处理连续变量时冲激8(t)的作用相同。与其连续形式不同的是，离散冲激是一个普通函数。

冲激串 $s_{\Delta T}(t)$ ,它定义为无限多个分离的周期冲激单元 $\Delta T$ 之和:

$s_{\Delta T}(t)=\sum_{n=-\infty }^{\infty }\delta (t-n\Delta T )$

冲激可以是连续的或离散的。

4、连续变量函数的傅里叶变换

由 $\Im \left \{ f(t) \right \}$ 表示的连续变量t的连续函数f(t)的傅里叶变换由下式定义:

$\Im \left \{ f(t) \right \}=\int_{-\infty }^{\infty }f(t)e^{-j2\pi \mu t}dt$

傅里叶变换对：

$F(\mu )=\int_{-\infty }^{\infty }f(t)e^{-j2\pi \mu t}dt$

$f(t )=\int_{-\infty }^{\infty }F(\mu )e^{j2\pi \mu t}du$

用欧拉公式;

$F(\mu )=\int_{-\infty }^{\infty }f(t)[\cos(2\pi \mu t)-j\sin(2\pi \mu t)]dt$

5、卷积

f(t)★h(t)和 H(u)F(u)是傅里叶变换对，卷积定理：

三、取样和取样函数的傅里叶变换

1、取样

模拟取样的一种方法是用一个AT单位间隔的冲激串作为取样函数去乘以f(t)，即

2、取样函数的傅里叶变换

F(u)和S(u)的卷积:

3、取样定理

如果以超过函数最高频率的两倍的取样率来获得样本，连续的带限函数可以完全地从它的样本集来恢复。一个连续的带限函数用取样率大于函数最高频率的两倍的取样来表示不会有信息损失。

4、混淆

由函数欠取样导致的被来自邻近周期的频率破坏这种效果就是周知的频率混淆，或简称为混淆。

5、由取样后的数据重建(复原)函数

完美重建的函数是用取样值加权的sinc函数的无限和，并且有重要的特性，即重建的函数恒等于在多个△T的整数增量处的样本值。

四、单变量的离散傅里叶变换(DFT)

1、由取样后的函数的连续变换得到DFT

离散傅里叶变换

傅里叶反变换(IDFT)

2、取样和频率间隔间的关系

五、两个变量的函数的扩展

1、二维冲激及其取样特性

两个连续变量t和z的冲激 $\delta (t,z)$ ：

二维冲激在积分下呈现了取样特性:

更一般地，对位于坐标 $(t_{0},z_{0})$ 处的冲激，有

取样特性在冲激所在的位置产生函数f(t,z)的值。

2、二维连续傅里叶变换对

3、二维取样和二维取样定理

二维取样可用取样函数(二维冲激串)建模:

4、图像中的混淆

一维混淆

图像内插和重取样

莫尔(波纹)模式

5、二维离散傅里叶变换及其反变换

二维傅里叶变换(DFT) :

傅里叶反变换(IDFT)：

六、二维离散傅里叶变换的一些性质

1、空间和频率间隔的关系

2、平移和旋转

3、周期性

4、对称性

5、傅里叶谱和相角

傅里叶谱：

相角：

6、二维卷积定理

将一维的扩展至两个变量可得到下面的二维循环卷积的表达式：

七、频率域滤波基础

1、频率域的其他特性

变化最慢的频率分量(u = v=0)与图像的平均灰度成正比。当我们远离变换的原点时，低频对应于图像中变化缓慢的灰度分量（例如，图像的平滑部分）。当我们从原点移开更远一些时，较高的频率开始对应于图像中越来越快的灰度变化（例如，边缘或噪声等尖锐部分）。

2、频率域滤波基础

频率域滤波由修改一幅图像的傅里叶变换然后计算其反变换得到处理后的结果组成。若给定一幅大小为M×N的数字图像f(x, y)，基本滤波公式有如下形式:

3、频率域滤波的基本步骤

①用 $(-1)^{x+y}$ 乘以输入图像f(x,y)，使其原点中心化
②对步骤①的结果执行DFT，得到关于中心对称的频谱F(uv)③生成一个实的、中心对称的频域滤波器H(u,v)
④对滤波器H(u,v)、频谱F(u,v)执行阵列相乘（对应元素逐个进行相乘)，形成乘积G(u,v)=H(u,v)F(u,v)，其中G(m,n)=H(m,n)F(m,n)且0≤m≤M-1,0≤n≤N-1
⑤对步骤④的结果G(u,v)执行反DFT，并取其结果的实部
⑥用 $(-1)^{x+y}$ 乘以步骤⑤的反DFT结果的实部，得到滤波结果g(x,y)