傅里叶变换
傅里叶变换的定义
傅里叶变换是将图像表示为不同幅度、频率和相位的复指数的和。傅里叶变换在广泛的图像处理应用中起着关键作用,包括增强、分析、恢复和压缩。
如果f(米,n)是两个离散空间变量的函数吗米而且n,则二维傅里叶变换的f(米,n)是由关系定义的吗
的变量ω1而且ω2为频率变量;单位是弧度/样本。F(ω1,ω2)常被称为频域的代表f(米,n).F(ω1,ω2)复值函数在内和内都是周期函数吗ω1而且ω2,带周期 .由于周期性,通常只有范围 会显示出来。请注意,F(0,0)是所有值的和吗f(米,n).出于这个原因,F(0,0)常被称为恒定的组件或直流分量傅里叶变换。(DC代表直流电;它是一个电气工程术语,指的是恒压电源,而不是电压呈正弦变化的电源。)
变换的逆是一种操作,当对变换后的图像执行时,会产生原始图像。二维傅里叶反变换由
粗略地说,这个方程意味着f(米,n)可以表示为无限个具有不同频率的复指数(正弦信号)的和。频率上贡献的幅度和相位(ω1,ω2)是由F(ω1,ω2).
可视化傅里叶变换
为了说明这一点,考虑一个函数f(米,n)它在矩形区域内等于1,其他地方都是0。为了简化图表,f(米,n)是一个连续函数,即使变量米而且n是离散的。
矩形函数
下图作为网格图显示了傅里叶变换的幅值,
上图所示的矩形函数。幅值的网格图是傅里叶变换可视化的常用方法。
矩形函数的模像
在情节中心的顶点是F(0,0)中所有值的和f(米,n).图表也表明了这一点F(ω1,ω2)在高水平频率比在高垂直频率有更多的能量。这反映了水平横截面的事实f(米,n)为窄脉冲,垂直截面为宽脉冲。窄脉冲的高频含量比宽脉冲高。
另一种常见的可视化傅里叶变换的方法是显示
如图所示。
矩形函数的傅里叶变换的对数
使用对数有助于得到傅里叶变换的细节F(ω1,ω2)非常接近0。
其他简单形状的傅里叶变换示例如下所示。
一些简单形状的傅里叶变换
离散傅里叶变换
在计算机上处理傅里叶变换通常涉及一种称为离散傅里叶变换(DFT)的变换形式。离散变换是一种输入和输出值为离散样本的变换,便于计算机操作。使用这种形式的转换有两个主要原因:
DFT的输入和输出都是离散的,便于计算机操作。
有一种计算DFT的快速算法称为快速傅里叶变换(FFT)。
DFT通常定义为一个离散函数f(米,n)它只在有限区域内是非零的 而且 .二维米——- - - - - -NDFT和逆米——- - - - - -NDFT关系由
而且
的值F(p,问)的DFT系数是f(米,n).零频率系数,F(0,0)通常称为“直流分量”。直流是电气工程术语,代表直流电。(注意MATLAB中的矩阵指标®总是从1开始,而不是0;因此,矩阵元素f(1,1)而且F(1,1)对应于数学量f(0,0)而且F(0,0)分别)。
MATLAB函数fft
,fft2
,fftn
实现了快速傅里叶变换算法,分别计算一维DFT、二维DFT和n维DFT。的函数传输线
,ifft2
,ifftn
计算DFT逆。
和傅里叶变换的关系
DFT系数F(p,问)是傅里叶变换的样本F(ω1,ω2).
可视化离散傅里叶变换
构造一个矩阵
f
这和函数很相似f(米,n)在傅里叶变换的定义.记住,f(米,n)在矩形区域内等于1,其他地方等于0。使用二值图像来表示f(米,n).F = 0 (30,30);F (5:24,13:17) = 1;imshow (f,‘InitialMagnification’,‘适合’)
计算和可视化的30 × 30 DFT
f
通过这些命令。F = fft2(F);F2 = log(abs(F));imshow (F2, [1 - 5], ' InitialMagnification ', '适应');colormap(飞机);colorbar
无填充的离散傅里叶变换计算
图中显示的傅里叶变换与图中的不同可视化傅里叶变换.首先,傅里叶变换的采样要粗糙得多。其次,零频系数显示在左上角,而不是传统的中心位置。
为了获得更精细的傅里叶变换采样,将零填充添加到
f
当计算其DFT时。零填充和DFT计算可以用这个命令在一个步骤中执行。F = fft2(F,256,256);
命令zero-pads
f
为256 × 256,然后再计算DFT。imshow(日志(abs (F)), [1 5]);colormap(飞机);colorbar
用填充计算的离散傅里叶变换
然而,零频率系数仍然显示在左上角,而不是中心。您可以通过使用函数来解决这个问题
fftshift
的象限互换F
零频率系数在中心。F = fft2(F,256,256);F2 = fftshift(F);imshow(日志(abs (F2)), [1 5]);colormap(飞机);colorbar
所得到的图与图中所示的图相同可视化傅里叶变换.
傅里叶变换的应用
本节介绍傅里叶变换的许多图像处理相关应用中的几个。
线性滤波器的频率响应
线性滤波器脉冲响应的傅里叶变换给出了滤波器的频率响应。这个函数freqz2
计算并显示滤波器的频率响应。高斯卷积核的频率响应表明,该滤波器通过低频,衰减高频。
H = fspecial('高斯');freqz2 (h)
高斯滤波器的频率响应
看到频域线性滤波器的设计有关线性滤波、滤波器设计和频率响应的更多信息。
使用傅里叶变换进行快速卷积
这个例子展示了如何使用傅里叶变换对两个矩阵进行快速卷积。傅里叶变换的一个关键性质是两个傅里叶变换的乘积对应于相关空间函数的卷积。这个性质,连同快速傅里叶变换,形成了快速卷积算法的基础。
注意:基于fft的卷积方法最常用于大输入。对于小的输入,通常使用imfilter
函数。
创建两个简单的矩阵,一个
而且B
.一个
是m × n矩阵B
是p × q矩阵。
A =魔术(3);B = ones(3);
在一个
而且B
使它们至少是(M+P-1)-by-(N+Q-1)。(通常是一个
而且B
0被填充为2的幂,因为fft2
对于这些尺寸是最快的。)该示例将矩阵填充为8 × 8。
A(8,8) = 0;B(8,8) = 0;
计算的二维DFT一个
而且B
使用fft2
函数。将两个DFT乘在一起,并使用ifft2
函数。
C = fft2(fft2(A).*fft2(B));
提取结果的非零部分,并去除舍入错误引起的虚部。
C = C(1:5,1:5);C =实(C)
C =5×58.0000 9.0000 15.0000 7.0000 6.0000 11.0000 17.0000 30.000 19.0000 13.0000 15.0000 30.000 45.0000 30.000 15.0000 7.0000 21.0000 30.000 23.0000 9.0000 4.0000 13.0000 15.0000 11.0000 2.0000
执行基于fft的相关性来定位图像特征
这个例子展示了如何使用傅里叶变换来进行相关,这与卷积密切相关。相关性可以用来定位图像中的特征。在这种情况下,通常称为相关性模板匹配.
将一个示例图像读入工作区。
Bw = imread(“text.png”);
通过从图像中提取字母“a”来创建一个匹配模板。类的交互语法也可以创建模板imcrop
函数。
A = bw(32:45,88:98);
通过将模板图像旋转180度,然后使用基于fft的卷积技术,计算模板图像与原始图像的相关性。(如果将卷积核旋转180度,卷积就等价于相关性。)要将模板与图像匹配,请使用fft2
而且ifft2
功能。在生成的图像中,明亮的峰值对应字母的出现。
C =实际(ifft2 (fft2 (bw)。* fft2 (rot90 (a, 2), 256256)));图imshow (C, [])将图像缩放到适当的显示范围。
若要查看模板在图像中的位置,请找到最大像素值,然后定义小于此最大值的阈值。阈值图像将这些峰值的位置显示为阈值相关图像中的白点。(为了使位置在图中更容易看到,示例将阈值图像放大以扩大点的大小。)
马克斯(C (:))
Ans = 68
脱粒= 60;使用略小于max的阈值。D = C >脱粒;Se = strel(“磁盘”5);E = imdilate(D,se);图imshow (E)显示值超过阈值的像素。