目录
序言:对图像而言,空间频率是指单位长度内亮度作为周期性变化的次数!
一、傅里叶变换基础
三、频率域滤波
四、二维离散变换的性质
五、实现
六、相关资源
序言:对图像而言,空间频率是指单位长度内亮度作为周期性变化的次数!
一、傅里叶变换基础
一维傅里叶变换数学推导
首先,我们知道傅里叶级数,形如:
是由三角函数系
经过组合所产生的三角级数。a0,an,bn 称为傅里叶系数。为了求傅里叶系数,我们先引出下面的预备知识。
三角函数系(2)中任何两个不相同的函数的乘积在 [?T/2,T/2]上的积分等于0。即:
而三角函数系(2)中任何两个相同函数的乘积在 [?T/2,T/2]上的积分都不等于0。
对(1)式在 [?T/2,T/2] 逐项积分,
由上面的预备知识,我们知道上式右边括号内的积分等于0,所以
故
现以 cos(nωx)乘(1)式两端,再同上逐项积分,我们得到:
等式右边除了
外,其他积分项都等于0,即
故
同理,用 sin(nωx)乘以(1) 两端并在 [?T/2,T/2]逐项积分,可得
三个系数求完后,我们尝试把它化为指数形式。这时我们需要用到欧拉公式
联立两式得到
将上式代入(1)中,替换 cos(nωx)和 sin(nωx),得到
我们令
整理一下得到
结合上面三个式子,我们可以写成如下:
故,将上式代入(3)式,我们得到傅里叶变换的指数公式:
二、傅里叶变换的理解
不了解傅里叶变换的可以先看这个链接
如何直观形象、生动有趣地给文科学生介绍傅里叶变换?
傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工(例如去噪)。最后还可以利用傅立叶反变换将这些频域信号转换成时域信号(得到去噪后的图片)。
在处理图像时,图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。如:大面积的沙漠在图像中是一片灰度变化缓慢的区域(图像比较柔和),对应的频率值很低;而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域,对应的频率值较高(有可能是噪声)。傅立叶变换在实际中有非常明显的物理意义,设 f 是一个能量有限的模拟信号,则其傅立叶变换就表示 f 的谱。从纯粹的数学意义上看,傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看,傅立叶变换是将图像从空间域转换到频率域。其逆变换是将图像从频率域转换到空间域。换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数。(这段摘自他人)
三、频率域滤波
空间域和频率域中的线性滤波的基础都是卷积定理
很简单,频率域滤波就是通过傅里叶变换后在频率域中进行滤波,再通过傅里叶反变换得到滤波后的图像。
频率域滤波的步骤:
7.通过的左上限提取M*N,得到最终的处理结果。
低通滤波器
低通滤波器:使低频通过而使高频衰减的滤波器
1. 被低通滤波的图像比原始图像少尖锐的细节部分而突出平滑过渡部分
2. 对比空间域滤波的平滑处理,如均值滤波器
其中满足 的点的轨迹为一个圆。
高通滤波器
高通滤波器:使高频通过而使低频衰减的滤波器
1.被高通滤波的图像比原始图像少灰度级的平滑过渡而突出边缘等细节部分
2.对比空间域的梯度算子、拉普拉斯算子
同样,我们可以知道,经过高通滤波后图像中的边缘和其他灰度急剧转变得到了增强,图像失去了大部分原来图像平缓的部分。
理想高通滤波器(IHPF):
布特沃斯高通滤波器(BHPF):
高斯高通滤波器(GHPF):
频率域拉普拉斯算子:对应时域使用二阶微分对图像进行锐化操作,由时域的拉普拉斯算子得频率域的拉普拉斯算子
选择性滤波:带阻滤波器、带通滤波器、陷滤波器等。
四、二维离散变换的性质
五、实现
- 二维DFT可分性:f(x,y)的二维DFT可通过计算f(x,y)的每一行的一维变换,然后沿计算结果的每一列计算一维变换得到。
- 使用DFT算法计算IDFT。
- 快速傅里叶变换FFT。
六、相关资源
为什么要进行傅里叶变换,其物理意义是什么?
傅里叶级数和傅里叶变换是什么关系?
如果看了此文你还不懂傅里叶变换,那就过来掐死我吧【完整版】
(精心整理)图像的傅里叶变换