傅里叶变换及其在图像处理中的应用
- 格式:doc
- 大小:368.00 KB
- 文档页数:5
傅里叶变换在图像去噪中的应用优化探讨图像去噪是数字图像处理领域中的一个重要问题,目的是通过消除图像中的噪声,恢复图像的清晰度和细节。
傅里叶变换作为一种有效的信号处理工具,在图像去噪中被广泛应用。
本文将探讨傅里叶变换在图像去噪中的应用优化方法。
一、傅里叶变换的基本原理傅里叶变换是将一个时域函数转化为其频域表示的一种数学变换方法。
在图像处理中,傅里叶变换可以将图像分解为一系列频率成分。
其基本公式如下:F(u, v) = ∬f(x, y)e^(-i2π(ux+vy))dxdy其中F(u, v)表示频域中的图像,f(x, y)表示时域中的图像。
傅里叶变换将图像从空间域转换到频域,使得频域中不同频率成分的信息可以更清晰地被提取和处理。
二、傅里叶变换在图像去噪中的应用图像去噪是通过去除图像中的噪声来提高图像质量的过程。
传统的图像去噪方法包括均值滤波、中值滤波等。
然而,这些方法往往会模糊图像细节,因此需要一种更加有效的方法来保持图像的清晰度。
傅里叶变换在图像去噪中的应用主要体现在频域滤波上。
通过将图像从空间域转换到频域,可以很容易地对图像进行频域滤波操作。
常见的频域滤波方法包括低通滤波和高通滤波。
低通滤波可以滤除图像中高频成分,从而去除图像中的噪声;高通滤波可以强调图像中的高频成分,使得图像的细节更加清晰。
三、傅里叶变换在图像去噪中的优化方法尽管傅里叶变换在图像去噪中具有广泛应用,但是它也存在一些问题,例如频谱泄漏、边缘模糊等。
为了优化傅里叶变换在图像去噪中的效果,研究人员提出了一些改进方法。
1. 加窗函数加窗函数可以有效缓解频谱泄漏问题。
常见的窗函数包括汉宁窗、汉明窗等。
通过在时域中对图像进行窗函数处理,可以减小傅里叶变换中的泄漏现象,从而提高去噪效果。
2. 频域滤波器设计传统的频域滤波器设计方法主要包括理想滤波器和巴特沃斯滤波器。
然而,这些方法会引入一些额外的问题,如振铃和削波等。
为了解决这些问题,研究人员提出了更加复杂的滤波器设计方法,如维纳滤波器和自适应滤波器。
傅里叶变换的五种应用场景傅里叶变换是一种重要的数学工具,在信号处理、图像处理、通信系统、物理学等领域都有广泛的应用。
本文将深入探讨傅里叶变换的五种应用场景,并分享对这些应用的观点和理解。
一、信号处理傅里叶变换在信号处理领域中扮演着不可或缺的角色。
信号可以是时间域中的连续信号也可以是离散信号,通过傅里叶变换可以将这些信号从时间域转化为频率域。
在频率域中,我们可以更清晰地观察信号的周期性和频谱特征。
这对于音频处理、图像处理、视频处理等都非常有用。
傅里叶变换的应用使得我们能够分析信号的频率成分、滤波去噪,甚至进行信号的压缩与解压缩。
二、图像处理图像处理是另一个广泛应用傅里叶变换的领域。
通过将图像进行傅里叶变换,我们可以将图像从空间域转换到频率域。
在频率域中,我们可以观察到图像中不同频率的成分,并对图像进行频率滤波、图像增强以及减少噪声的操作。
傅里叶变换的应用还包括图像压缩和图像恢复等方面。
例如,在JPEG图像压缩中,傅里叶变换被用来将图像编码成频域数据,从而实现图像的压缩。
三、通信系统在通信系统中,傅里叶变换起着至关重要的作用。
通过将信号进行傅里叶变换,我们可以将信号转换到频率域,进而对信号进行调制、解调、频谱分析等。
例如,正交频分多路复用技术(OFDM)是一种常用于现代通信系统中的调制技术。
OFDM基于傅里叶变换将高速数据流分成多个低速子流,并在不同频率上进行传输。
傅里叶变换的应用使得OFDM技术能够高效地利用频谱资源和抵御多径干扰。
四、物理学在物理学中,傅里叶变换也是一种应用广泛的数学工具。
不同物理现象可以通过傅里叶变换转换到频率域进行分析。
例如,在声学领域中,通过对声音信号进行傅里叶变换,我们可以观察到声音的频谱成分,从而对声音进行分析和处理。
在量子力学领域,傅里叶变换也被广泛应用于波函数的分析和计算。
五、其他领域除了上述提到的领域,傅里叶变换还在其他各个科学领域有着重要的应用。
例如,在生物医学领域中,傅里叶变换被用于对生物信号(如心电图、脑电图)进行频谱分析与滤波处理,以便提取有价值的信息。
图像处理之傅⾥叶变换图像处理之傅⾥叶变换⼀、傅⾥叶变换傅⾥叶变换的作⽤:⾼频:变化剧烈的灰度分量,例如边界低频:变化缓慢的灰度分量,例如⼀⽚⼤海滤波:低通滤波器:只保留低频,会使得图像模糊⾼通滤波器:只保留⾼频,会使得图像细节增强OpenCV:opencv中主要就是cv2.dft()和cv2.idft(),输⼊图像需要先转换成np.float32 格式。
得到的结果中频率为0的部分会在左上⾓,通常要转换到中⼼位置,可以通过shift变换来实现。
cv2.dft()返回的结果是双通道的(实部,虚部),通常还需要转换成图像格式才能展⽰(0,255)。
import numpy as npimport cv2from matplotlib import pyplot as pltimg = cv2.imread('lena.jpg',0)img_float32 = np.float32(img)dft = cv2.dft(img_float32, flags = cv2.DFT_COMPLEX_OUTPUT)dft_shift = np.fft.fftshift(dft)# 得到灰度图能表⽰的形式magnitude_spectrum = 20*np.log(cv2.magnitude(dft_shift[:,:,0],dft_shift[:,:,1]))plt.subplot(121),plt.imshow(img, cmap = 'gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122),plt.imshow(magnitude_spectrum, cmap = 'gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.show()import numpy as npimport cv2from matplotlib import pyplot as pltimg = cv2.imread('lena.jpg',0)img_float32 = np.float32(img)dft = cv2.dft(img_float32, flags = cv2.DFT_COMPLEX_OUTPUT) #时域转换到频域dft_shift = np.fft.fftshift(dft) #将低频部分拉到中⼼处rows, cols = img.shapecrow, ccol = int(rows/2) , int(cols/2) #确定掩膜的中⼼位置坐标# 低通滤波mask = np.zeros((rows, cols, 2), np.uint8)mask[crow-30:crow+30, ccol-30:ccol+30] = 1# IDFTfshift = dft_shift*mask #去掉⾼频部分,只显⽰低频部分f_ishift = np.fft.ifftshift(fshift) #将低频部分从中⼼点处还原img_back = cv2.idft(f_ishift) #从频域逆变换到时域img_back = cv2.magnitude(img_back[:,:,0],img_back[:,:,1]) #该函数通过实部和虚部⽤来计算⼆维⽮量的幅值plt.subplot(121),plt.imshow(img, cmap = 'gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122),plt.imshow(img_back, cmap = 'gray')plt.title('Result'), plt.xticks([]), plt.yticks([])plt.show()img = cv2.imread('lena.jpg',0)img_float32 = np.float32(img)dft = cv2.dft(img_float32, flags = cv2.DFT_COMPLEX_OUTPUT) dft_shift = np.fft.fftshift(dft)rows, cols = img.shapecrow, ccol = int(rows/2) , int(cols/2) # 中⼼位置# ⾼通滤波mask = np.ones((rows, cols, 2), np.uint8)mask[crow-30:crow+30, ccol-30:ccol+30] = 0# IDFTfshift = dft_shift*maskf_ishift = np.fft.ifftshift(fshift)img_back = cv2.idft(f_ishift)img_back = cv2.magnitude(img_back[:,:,0],img_back[:,:,1]) plt.subplot(121),plt.imshow(img, cmap = 'gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122),plt.imshow(img_back, cmap = 'gray')plt.title('Result'), plt.xticks([]), plt.yticks([])plt.show()。
【数字图像处理】傅⾥叶变换在图像处理中的应⽤1.理解⼆维傅⾥叶变换的定义1.1⼆维傅⾥叶变换1.2⼆维离散傅⾥叶变换1.3⽤FFT计算⼆维离散傅⾥叶变换1.3图像傅⾥叶变换的物理意义2.⼆维傅⾥叶变换有哪些性质?2.1⼆维离散傅⾥叶变换的性质2.2⼆维离散傅⾥叶变换图像性质3.任给⼀幅图像,对其进⾏⼆维傅⾥叶变换和逆变换4.附录 94.1matlab代码4.2参考⽂献⽬录1.理解⼆维傅⾥叶变换的定义1.1⼆维傅⾥叶变换⼆维Fourier变换:逆变换:1.2⼆维离散傅⾥叶变换⼀个图像尺⼨为M×N的函数的离散傅⾥叶变换由以下等式给出:其中和。
其中变量u和v⽤于确定它们的频率,频域系统是由所张成的坐标系,其中和⽤做(频率)变量。
空间域是由f(x,y)所张成的坐标系。
可以得到频谱系统在频谱图四⾓处沿和⽅向的频谱分量均为0。
离散傅⾥叶逆变换由下式给出:令R和I分别表⽰F的实部和需部,则傅⾥叶频谱,相位⾓,功率谱(幅度)定义如下:1.3⽤FFT计算⼆维离散傅⾥叶变换⼆维离散傅⾥叶变换的定义为:⼆维离散傅⾥叶变换可通过两次⼀维离散傅⾥叶变换来实现:1)作⼀维N点DFT(对每个m做⼀次,共M次)2)作M点的DFT(对每个k做⼀次,共N次)这两次离散傅⾥叶变换都可以⽤快速算法求得,若M和N都是2的幂,则可使⽤基⼆FFT算法,所需要乘法次数为⽽直接计算⼆维离散傅⾥叶变换所需的乘法次数为(M+N)MN,当M和N⽐较⼤时⽤⽤FFT运算,可节约很多运算量。
1.3图像傅⾥叶变换的物理意义图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平⾯空间上的梯度。
如:⼤⾯积的沙漠在图像中是⼀⽚灰度变化缓慢的区域,对应的频率值很低;⽽对于地表属性变换剧烈的边缘区域在图像中是⼀⽚灰度变化剧烈的区域,对应的频率值较⾼。
傅⾥叶变换在实际中有⾮常明显的物理意义,设f是⼀个能量有限的模拟信号,则其傅⾥叶变换就表⽰f的频谱。
从纯粹的数学意义上看,傅⾥叶变换是将⼀个函数转换为⼀系列周期函数来处理的。
傅里叶定律应用实例傅里叶定律是一种将任意周期性函数分解为一组正弦或余弦函数的方法。
它有非常广泛的应用,例如在信号处理、图像处理、量子力学、声音波谱分析等领域。
1. 信号处理和音频压缩傅里叶定律可以用于压缩音频或其他信号。
通过将信号分解为一组正弦或余弦函数的和,可以找到一个足够小的子集来代表原始信号。
这使得信号的存储空间更小,并且可以更快地传输。
现代音频压缩算法如MP3就使用了傅里叶变换来分解音频信号。
2. 图像处理在图像处理中,傅里叶变换可以用来分析和处理图像。
通过将图像分解为其频率成分,可以实现许多图像处理操作,例如去噪、过滤、锐化和边缘检测。
傅里叶变换还可以用于图像压缩,通常与离散余弦变换(DCT)结合使用。
3. 量子力学傅里叶变换在量子力学中也有广泛的应用。
傅里叶变换可以用于将一个波函数从空间域转换为能量域,这对于解决一些量子力学问题非常有用。
傅里叶变换还可以用于分析和处理量子力学中的能级和自旋。
4. 声音波谱分析傅里叶变换可以用于分析声音波形成分的频率。
在声音波形中,每个频率成分可以表示为正弦或余弦波的组合。
通过使用傅里叶变换,可以将波形转换为频域,以便更好地理解声音的波形结构。
除了上述应用,傅里叶定律还有其他一些重要的作用。
下面进一步探讨一下它在不同领域的应用:5. 数字信号处理傅里叶变换在数字信号处理中扮演着非常重要的角色。
通过将信号从时域转换为频域,可以更好地理解信号的性质和特征。
可以使用傅里叶变换来从一个信号中分离出特定的频率成分,以便更好地对信号进行分析。
6. 机器学习在机器学习中,傅里叶变换可以用来处理图像和声音等数据。
可以使用傅里叶变换将图像从空间域转换为频域,以便更好地识别图像中的模式和特征。
同样地,傅里叶变换也可以用来处理声音数据,以便更好地识别声音信号中的模式和特征。
7. 通信系统在通信系统中,傅里叶变换可以用于信号传输和处理。
通过分析信号频率成分,可以更好地理解信号的性质,并且可以更好地设计和优化通信系统。
傅里叶变换在图像处理中的应用摘要傅里叶变换是一种重要的信号分析工具,在平稳信号的分析方面具有十分重要的地位,线性系统中,常利用傅里叶变换进行分析和处理。
本文对傅里叶变换和数字图像处理的相关概念进行了介绍,并主要针对傅里叶变换在数字图像处理中的应用进行分析和研究,对图像处理领域的学习很有帮助。
关键词傅里叶变换;信号分析;平稳信号;数字图像处理前言随着信号处理领域的不断发展,越来越多信号分析工具得到了相关学者的研究。
傅里叶变换于19世纪就已经被研究人员提出,在之后的研究和应用中,傅里叶变换也一直是重要的信号处理工具[1-2]。
信息时代的到来使数字图像处理技术也开始飞速进步,它与信号处理等技术息息相关,因此傅里叶变换在图像处理中也得到了重要的应用[3]。
传统的处理方式往往适合在时域对图像进行处理分析,而与傅里叶变换相结合便使图像处理技术得以在频域进行,傅里叶变换常用于线性系统中的处理,因此,可以很好地和图像处理领域相联系,有效提高数字图像处理的效率和精度[4]。
1 傅里叶变换的概述最早在1807年,法国工程师傅里叶首先提出了有关傅里叶级数这一理论,首次提到可以將一个周期性的信号展开成多个复正弦信号相加的形式,这一理论引起了学者们的注意。
十几年之后,傅里叶正式提出了傅里叶变换的概念。
通过傅里叶变换,我们可以将一个信号由时域转换到频域进行信号处理和分析,并且通过傅里叶变换的提出才加深了人们对于频率这个概念的理解。
因此,在傅里叶变换被提出之后,在信号分析领域提出了从频域进行分析这个新思路,使人们对信号的特性进行了一些新的方面的研究。
很多对信号的处理问题以往通过时域分析很难真的得到充分的解释,傅里叶变换这个思路使很多问题变得显而易见。
对于傅里叶变换之后的研究中,出现了关于傅里叶变换的快速算法,使得傅里叶变换更加具有实际应用价值,也对处理离散的数字信号起了重要的作用。
2 基于傅里叶变换的图像处理在对图像进行处理的过程中,图像中包含许多线性变化的元素,而其中的频率便是十分重要的物理量,而这种包含频率信息的元素正适合应用傅里叶变换进行处理,因此,傅里叶变换在图像处理领域得到了广泛的应用。
图像处理与傅里叶变换1背景傅里叶变换是一个非常复杂的理论,我们在图像处理中集中关注于其傅里叶离散变换离散傅立叶变换(Discre t e Fourie r Transf o rm) 。
1.1离散傅立叶变换图象是由灰度(R GB )组成的二维离散数据矩阵,则对它进行傅立叶变换是离散的傅立叶变换。
对图像数据f (x,y)(x=0,1,… ,M-1; y=0,1,… ,N-1)。
则其离散傅立叶变换定义可表示为:式中,u=0,1,…, M-1;v= 0,1,…, N-1其逆变换为式中,x=0,1,…, M-1;y= 0,1,…, N-1在图象处理中,一般总是选择方形数据,即M=N影像f(x,y)的振幅谱或傅立叶频谱: 相位谱: 能量谱(功率谱) )1(2exp ),(1),(1010∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+-=M x N y N vy M ux i y x f MN v u F π)2(2exp ),(1),(1010∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+=M u N v N vy M ux i v u F MN y x f π),(),(),(22v u I v u R v u F +=[]),(/),(),(v u R v u I arctg v u =ϕ),(),(),(),(222v u I v u R v u F v u E +==1.2快速傅里叶变化可分离性的优点是二维的傅立叶变换或逆变换由两个连续的一维傅立叶变换变换来实现,对于一个影像f (x,y),可以先沿着其每一列求一维傅立叶变换,再对其每一行再求一维变换正变化逆变换 由于二维的傅立叶变换具有可分离性,故只讨论一维快速傅立叶变换。
正变换逆变换由于计算机进行运算的时间主要取决于所用的乘法的次数。
按照上式进行一维离散由空间域向频率域傅立叶变换时,对于N 个F ∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡⨯⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+=10101010)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N v N u N u N v N vy i v u F N N ux i v u F N N vy ux i v u F NN y x f πππ∑-=⎥⎦⎤⎢⎣⎡-=102exp )(1)(N x N ux i x f N u F π∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡-⨯⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡+-=10101010)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N y N x N x N y N vy i y x f N N ux i y x f NN vy ux i y x f NN v u F πππ∑-=⎥⎦⎤⎢⎣⎡=102exp )(1)(N u N ux i u F N x f π(u)值,中的每一个都要进行N 次运算,运算时间与N 2成正比。
fft快速傅里叶变换应用场景一、引言傅里叶变换是信号处理中常用的基本工具之一,它可以将时域信号转化为频域信号,从而对信号进行频谱分析。
但是,传统的傅里叶变换算法计算复杂度较高,对于实时性要求较高的应用场景不太适合。
因此,快速傅里叶变换(FFT)应运而生。
本文将介绍FFT快速傅里叶变换在各种应用场景中的具体应用。
二、图像处理1. 图像压缩图像压缩是指通过某种算法将图像数据压缩到更小的存储空间中,以减少存储空间和传输带宽。
FFT快速傅里叶变换可以将图像从时域转化为频域,然后对频域信息进行压缩。
这样做的好处是可以去除一些高频成分和低频成分,从而减少冗余数据。
2. 图像滤波图像滤波是指通过某种算法对图像进行降噪或增强处理。
FFT快速傅里叶变换可以将图像从时域转化为频域,在频域中进行滤波操作。
例如,在高通滤波器中,可以将低频成分滤除,从而增强图像的高频细节。
三、音频处理1. 音频压缩音频压缩是指通过某种算法将音频数据压缩到更小的存储空间中,以减少存储空间和传输带宽。
FFT快速傅里叶变换可以将音频从时域转化为频域,然后对频域信息进行压缩。
这样做的好处是可以去除一些高频成分和低频成分,从而减少冗余数据。
2. 音乐合成音乐合成是指通过某种算法将多个声音信号合并为一个复合声音信号。
FFT快速傅里叶变换可以将多个声音信号从时域转化为频域,在频域中进行加和操作。
这样做的好处是可以避免在时域中信号相加时出现相位问题。
四、通信领域1. 无线电通信在无线电通信中,FFT快速傅里叶变换被广泛应用于OFDM(正交分组多路复用)调制技术中。
OFDM技术利用FFT技术将高速数据流分割成多个低速子载波,在每个子载波上进行调制和解调,从而提高了无线电信号的传输速率和抗干扰能力。
2. 有线通信在有线通信中,FFT快速傅里叶变换被广泛应用于数字信号处理中。
例如,在数字电视中,FFT技术可以将视频和音频数据分离出来,从而实现高清晰度的视频和清晰的声音。
图像处理与傅里叶变换1背景傅里叶变换是一个非常复杂的理论,我们在图像处理中集中关注于其傅里叶离散变换离散傅立叶变换(Discrete Fourier Transform) 。
1.1离散傅立叶变换图象是由灰度(RGB )组成的二维离散数据矩阵,则对它进行傅立叶变换是离散的傅立叶变换。
对图像数据f(x,y)(x=0,1,… ,M-1; y=0,1,… ,N-1)。
则其离散傅立叶变换定义可表示为:式中,u=0,1,…, M-1;v= 0,1,…, N-1 其逆变换为式中,x=0,1,…, M-1;y= 0,1,…, N-1在图象处理中,一般总是选择方形数据,即M=N影像f(x,y)的振幅谱或傅立叶频谱: 相位谱:能量谱(功率谱) )1(2exp ),(1),(101∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+-=M x N y N vy M uxi y x f MNv u F π)2(2exp ),(1),(101∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+=M u N v N vy M uxi v u F MNy x f π),(),(),(22v u I v u R v u F +=[]),(/),(),(v u R v u I arctg v u =ϕ),(),(),(),(222v u I v u R v u F v u E +==1.2快速傅里叶变化可分离性的优点是二维的傅立叶变换或逆变换由两个连续的一维傅立叶变换变换来实现,对于一个影像f(x,y),可以先沿着其每一列求一维傅立叶变换,再对其每一行再求一维变换正变化逆变换由于二维的傅立叶变换具有可分离性,故只讨论一维快速傅立叶变换。
正变换 逆变换由于计算机进行运算的时间主要取决于所用的乘法的次数。
按照上式进行一维离散由空间域向频率域傅立叶变换时,对于N 个F∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡⨯⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+=110101)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N v N u N u N v N vy i v u F NN ux i v u F N N vy ux i v u F NNy x f πππ∑-=⎥⎦⎤⎢⎣⎡-=12exp )(1)(N x N ux i x f Nu F π∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡-⨯⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡+-=11101)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N y N x N x N y N vy i y x f NN ux i y x f NN vy ux i y x f NNv u F πππ∑-=⎥⎦⎤⎢⎣⎡=12exp )(1)(N u N ux i u F Nx f π(u)值,中的每一个都要进行N 次运算,运算时间与N 2成正比。
傅里叶变换及其在数字图像处理中的应用
王家硕 学号:1252015
一、 Fourier 变换
1. 一维连续傅里叶变换
设 f (x)为x 的实变函数,如果f (x)满足下面的狄里赫莱条件: (1)具有有限个间隔点。
(2)具有有限个极点。
(3)绝对可积。
则 f (x )的傅里叶变换(Fourier Transformation ,FT )定义为: Fourier 正变换:dt e t f t f f F t j ⎰
+∞
∞
--=
=ωω)()]([)(;
Fourier 逆变换:ωωπ
ωd e f t F f t f t j ⎰
∞
+∞
---=
=)(21)]([)(1
,
式中:1-=
j ,ω 为频域变量。
f (x )与F (w )构成傅里叶变换对,可以证明傅里叶变换对总是存在的。
由于f (x )为实函数,则它的傅里叶变换F (w )通常是复函数,于是F (w )可写成
F (w ) = R (w ) + j I (w ) (1)
式中:R (w )和I (w )分别是F (w )的实部和虚部。
公式1可表示为指数形式:
式中:
F (w ) 为f (x )的傅里叶幅度谱,f (w )为f (x )的相位谱。
2. 二维连续傅里叶变换
如果二维函数f (x , y )是连续可积的,即∞<⎰⎰
+∞
∞
-dxdy y x f |),(,且F (u , v )是可积的,
则二维连续傅里叶变换对可表示为:
dt e y x f v u F t j ⎰
⎰+∞
∞
--+∞∞
-=
ω),(),(
dt e v u F y x F t j ⎰
⎰
∞
+∞
-∞
+∞
-=
ω),(),(
对于图像 f (x, y),F(u, v)是它的频谱。
变量u 是对应于x 轴的空间频率,变量v 是
对应于y 轴的空间频率,与在一维的情况类似,可定义二维傅里叶变换的幅度谱和相位谱为:
3.一维离散傅里叶变换
对一个连续函数f (x)等间隔采样可得到一个离散序列。
设共采样N个,则这个离散序列可表示为{ f (0), f (1), f (2), , f (N -1)}。
则其离散傅里叶变换F(u)为:
离散傅里叶反变换(IDFT)为:
式中:x = 0, 1, 2, , N -1。
令,则上述公式变成:
4.二维离散傅里叶变换
二维离散傅里叶变换:
二、傅里叶变换在图像处理中的应用
傅立叶变换在图像处理中有非常重要的作用,被广泛应用于图像增强与图像去噪、图像分割之边缘检测、图像特征提取(形状、纹理)、图像压缩等方面。
1.基于傅里叶变换的图像增强
在图像处理中,图像高频分量:图像突变部分;在某些情况下指图像边缘信息,某些情况下指噪声,更多是两者的混合;低频分量:图像变化平缓的部分,也就是图像轮廓信息。
高通滤波器:让图像使低频分量抑制,高频分量通过。
低通滤波器:与高通相反,让图像使高频分量抑制,低频分量通过。
图像频域滤波增强技术是在频率域空间对图像进行滤波,因此需要将图像从空间域通过傅里叶变换频率域,具体操作如下:
假定原图像f(x,y),经傅里叶变换为F(u,v),频率域增强就是选择合适的滤波器函数H(u,v)对F(u,v)的频谱成分进行调整,然后经逆傅里叶变换得到增强的图像g(x,y)。
该过程可以通过下面的流程描述:
(1)对原始图像f(x,y)进行傅里叶变换得到F(u,v);
(2)将F(u,v)与传递函数H(u,v)进行卷积运算得到g(u,v);
(3)将g(u,v)进行傅里叶逆变换得到增强图像g(x,y)。
频域滤波的核心在于如何确定传递函数,即H(u,v)。
常用的频率域低通滤波器H(u,v)有4种。
理想低通滤波器、巴特沃斯(Butterworth)低通滤波器、指数低通滤波器、梯形低通滤波器。
2.图像增强MatLab实现
clear all;
I1=imread('cameraman.tif'); %转载图像
fftI1=fft2(I1); %二维离散傅立叶变换
sfftI1=fftshift(fftI1); %直流分量移到频谱中心
RR1=real(sfftI1); %取傅立叶变换的实部
II1=imag(sfftI1); %取傅立叶变换的虚部
A1=sqrt(RR1.^2+II1.^2); %计算频谱幅值
A1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*225;%归一化
subplot(3,2,2);imshow(A1); %显示原图像的频谱
figure,imshow(I1,[]); %把图像显示出来
% I1=imnoise(I1,'salt & pepper');
% figure,imshow(I1);
snoise=0.1*randn(size(I1));
I1=imadd(I1,im2uint8(snoise)); % 受随机噪声干扰
figure,imshow(I1); %显示噪声图像
f=double(I1); %图像存储类型转换
g=fft2(f); %傅立叶变换
g=fftshift(g); %转换数据矩阵
[N1,N2]=size(g); %测量图像尺寸参数
n=2;
d0=50;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
for j=1:N2
d=sqrt((i-n1)^+(j-n2)^2)
h=1/(1+0.414*(d/d0)^(2*n)); %计算Butterworth低通转换函数 result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
X3=uint8(real(X2));
figure,imshow(X3) % 显示频域增强后的图像
3.运行结果
原始图像原始图像频谱
加入噪声图像 增强图像。