数字图像频域变换.ppt
《数字图像频域变换.ppt》由会员分享,可在线阅读,更多相关《数字图像频域变换.ppt(79页珍藏版)》请在三一办公上搜索。
1、数字图象处理,北京大学计算机研究所 陈晓鸥,第三节 频域变换,傅立叶变换导言理论基础、连续与离散的傅立叶变换二维傅立叶变换特性可分离性、周期与共轭对称、平移性、旋转特性、线性与相似性、均值性、拉普拉斯、卷积与相关快速傅立叶变换FFT算法、逆向FFT算法、算法实现,第三节 频域变换,2.2.1 傅立叶变换导言理论基础连续与离散的傅立叶变换,第三节 频域变换:理论基础,理论基础线性系统卷积与相关,第三节 频域变换:理论基础,线性系统系统的定义:接受一个输入,并产生相应输出的任何实体。系统的输入是一个或两个变量的函数,输出是相同变量的另一个函数。,系统,x(t)输入,y(t)输出,第三节 频域变换:
2、理论基础,线性系统线性系统的定义:对于某特定系统,有:x1(t)y1(t)x2(t)y2(t)该系统是线性的当且仅当:x1(t)+x2(t)y1(t)+y2(t)从而有:a*x1(t)a*y1(t),第三节 频域变换:理论基础,线性系统线性系统移不变性的定义:对于某线性系统,有:x(t)y(t)当输入信号沿时间轴平移T,有:x(t-T)y(t-T)则称该线性系统具有移不变性,第三节 频域变换:理论基础,卷积卷积的定义离散一维卷积二维卷积的定义离散二维卷积相关的定义,第三节 频域变换:理论基础,卷积的定义对于一个线性系统的输入f(t)和输出h(t),如果有一个一般表达式,来说明他们的关系,对线性
3、系统的分析,将大有帮助卷积积分就是这样的一般表达式 h(t)=g(t-)f()d 记为:h=g*f-g(t)称为冲激响应函数,第三节 频域变换:理论基础,离散一维卷积 h(i)=f(i)*g(i)=f(j)g(i-j)j二维卷积的定义 h(x,y)=f*g=f(u,v)g(x u,y v)dudv-,第三节 频域变换:理论基础,离散二维卷积h(x,y)=f*g=f(m,n)g(x m,y n)m n相关的定义 h(t)=g(t+)f()d 记为:y=g x-,第三节 频域变换,连续与离散的傅立叶变换一维连续傅立叶变换二维连续傅立叶变换离散傅立叶变换离散傅立叶变换的计算与显示,第三节 频域变换:
4、傅立叶变换,一维连续傅立叶变换:定义设 f(x)为实变量x的连续函数,f(x)的傅立叶变换表示为Ff(x),定义为:Ff(x)=F(u)=f(x)exp(-j2ux)dx 其中 j2=-1-如果给定F(u),f(x)可以由傅立叶逆变换得到:FF(u)=f(x)=F(u)exp(j2ux)du,第三节 频域变换:傅立叶变换,一维连续傅立叶变换:几个概念 假设函数f(x)为实函数。但一个实函数的傅立叶变换可能为复函数:F(u)=R(u)+jI(u)(1)f(x)的傅立叶模记为:|F(u)|F(u)|=R2(u)+I2(u)1/2(2)f(x)的傅立叶模平方记为:P(u)P(u)=|F(u)|2=R
5、2(u)+I2(u),第三节 频域变换:傅立叶变换,一维连续傅立叶变换:几个概念(3)f(x)的傅立叶相位记为:(u)(u)=tan-1(I(u)/R(u)(4)傅立叶变换中的变量u通常称为频率变量 这个名称源于尤拉公式中的指数项 exp-j2ux=cos2ux-jsin2ux 如果把傅立叶变换的积分解释为离散项的和,则易推出F(u)是一组sin和cos函数项的无限和,其中u的每个值决定了其相应cos,sin函数对的频率。,第三节 频域变换:傅立叶变换,二维连续傅立叶变换 如果f(x,y)连续可积,并且F(u,v)可积,则存在以下傅立叶变换对,其中u,v为频率变量:Ff(x,y)=F(u,v)
6、=f(x,y)exp-j2(ux+vy)dxdy-FF(u,v)=f(x,y)=F(u,v)expj2(ux+vy)dudv-,第三节 频域变换:傅立叶变换,二维连续傅立叶变换 二维傅立叶模、相位和模平方分别为:模:|F(u,v)|=R2(u,v)+I2(u,v)1/2 相位:(u,v)=tan-1(I(u,v)/R(u,v)模平方:P(u,v)=|F(u,v)|2=R2(u,v)+I2(u,v),第三节 频域变换:傅立叶变换,离散傅立叶变换 假设连续函数f(x),通过取N个x单位的采样点,被离散化为一个序列:f(x0),f(x0+x),f(x0+2x),f(x0+N1 x)无论将x作为离散的
7、或连续的变量,在子序列中来研究都将是方便的,仅仅依赖于讨论的上下文。为作到此要求定义:f(x)=f(x0+xx),第三节 频域变换:傅立叶变换,离散傅立叶变换 其中假设x现在的离散值是:0,1,2,N-1。f(x0),f(x0+x),f(x0+2x),.,f(x0+N1x)表示相对与连续函数的任意N个统一的空间采样,第三节 频域变换:傅立叶变换,离散傅立叶变换函数f(x0+xx)的离散傅立叶变换对有:N-1F(u)=1/N f(x)exp-j2ux/N x=0u=0,1,2,.N-1 N-1 f(x)=F(u)expj2ux/N u=0 x=0,1,2,.N-1,第三节 频域变换:傅立叶变换,
8、离散傅立叶变换:二维 M-1 N-1F(u,v)=1/MNf(x,y)exp-j2(ux/M+vy/N)x=0 y=0u=0,1,2,M-1;v=0,1,2,.N-1 M-1 N-1 f(x,y)=F(u,v)expj2(ux/M+vy/N)u=0 v=0 x=0,1,2,.N-1;y=0,1,2,.N-1,第三节 频域变换:傅立叶变换,离散傅立叶变换的计算与显示离散傅立叶变换的计算举例离散傅立叶变换的显示,第三节 频域变换:傅立叶变换,离散傅立叶变换的计算举例,x,f(x0)=f(x0+x),0,1,2,3,1,2,3,4,第三节 频域变换:傅立叶变换,F(0)=1/4f(x)exp0=1/
9、4f(0)+f1(1)+f(2)+f(3)=1/4(2+3+4+4)=3.25F(1)=1/4f(x)exp-j2x/4)=1/4(2e0+3e j2/4+4e j22/4+4e j23/4)=1/4(-2+j)F(2)=-1/4(1+j0)F(3)=-1/4(2+j),第三节 频域变换:傅立叶变换,离散傅立叶变换的计算举例因为,函数f(x,y)的傅立叶变换是f(x,y)积分的函数,因此计算每一个傅立叶变换值,原函数f(x,y)的每一个点都需要参予,第三节 频域变换:傅立叶变换,离散傅立叶变换的显示 通过对傅立叶变换模,来显示傅立叶变换图象。由于模的值域大于显示的值域,因此要进行动态值域的压缩
10、D(u,v)=c log(1+|F(u,v)|)其中:c=255/k;k=max(log(1+|F(u,v)|)值域0,k的上限(最大值),第三节 频域变换:傅立叶变换,离散傅立叶变换的显示,第三节 频域变换:傅立叶变换,离散傅立叶变换的显示,第三节 频域变换:二维傅立叶变换特性,2.2.2 二维傅立叶变换特性可分离性周期与共轭对称平移性旋转特性,线性与相似性 均值性 拉普拉斯 卷积与相关,第三节 频域变换:二维傅立叶变换特性,可分离性二维离散傅立叶变换DFT可分离性的基本思想是:二维DFT可分离为两次一维DFT应用:二维快速傅立叶算法FFT,是通过计算两次一维FFT实现的,第三节 频域变换:
11、二维傅立叶变换特性,可分离性的定义 M-1 N-1 F(u,v)=1/MN f(x,y)e(-j2vy/N)e(-j2ux/M)x=0 y=0u=0,1,2,M-1;v=0,1,2,.N-1 M-1 N-1 f(x,y)=F(u,v)e(j2vy/N)e(j2ux/M)u=0 v=0 x=0,1,2,.N-1;y=0,1,2,.N-1,第三节 频域变换:二维傅立叶变换特性,可分离性成立的推导先对行(y变量)做变换:N-1F(x,v)=1/Nf(x,y)exp(-j2vy/N)y=0然后对列(x变量)进行变换:M-1F(u,v)=1/MF(x,v)exp(-j2ux/M)x=0,第三节 频域变换
12、:二维傅立叶变换特性,先对行做变换:,然后对列进行变换:,f(x,y),(0,0),(N-1,M-1),x,y,F(x,v),(0,0),(N-1,M-1),x,v,F(x,v),(0,0),(N-1,M-1),x,v,F(u,v),(0,0),(N-1,M-1),u,v,第三节 频域变换:二维傅立叶变换特性,平移性定理平移性的描述函数自变量的位移的傅立叶变换产生一个复系数 Ff(x-a,y-b)=exp-j2(au+bv)F(u,v),第三节 频域变换:二维傅立叶变换特性,平移性成立的证明用一维函数为例进行证明:设位移为a,f(x-a)的傅立叶变换为:Ff(x-a)=F(u)=f(x-a)e
13、xp(-j2ux)dx-将积分乘以 1=exp(-j2au)exp(j2au)且设:v=x-a,dv=dx,第三节 频域变换:二维傅立叶变换特性,平移性成立的证明Ff(x-a)=F(u)=f(x-a)exp(-j2ux)exp(-j2au)exp(j2au)dx-=exp(-j2au)f(x-a)exp(-j2ux)exp(j2ua)dx-,第三节 频域变换:二维傅立叶变换特性,=exp(-j2au)f(x-a)exp-j2u(x-a)dx-=exp(-j2au)f(v)exp-j2uvdv-=exp(-j2au)F(u),第三节 频域变换:二维傅立叶变换特性,周期与共轭对称周期性的描述:离散
14、傅立叶变换DFT和它的逆变换是以N为周期的对于一维傅立叶变换有:F(u)=F(u+N)对于二维傅立叶变换有:F(u,v)=F(u+M,v+N),第三节 频域变换:二维傅立叶变换特性,周期与共轭对称共轭对称性的描述:傅立叶变换结果是以原点为中心的共轭对称函数对于一维傅立叶变换有:F(u)=F*(-u)对于二维傅立叶变换有:F(u,v)=F*(-u,-v),第三节 频域变换:二维傅立叶变换特性,共轭对称性证明以一维傅立叶变换为例证明:F(u)=f(x)exp-j2uxdx=f(x)expj2(-u)xdx=f(x)exp-j2(-u)x*dx(取共轭复数)=F*(-u),第三节 频域变换:二维傅立
15、叶变换特性,旋转特性旋转特性描述:如果f(x,y)旋转了一个角度,那么f(x,y)旋转后的图象的傅立叶变换也旋转了相同的角度。设:a(x,y)=x cos()-y sin()b(x,y)=x sin()+y cos()Ff(a(x,y),b(x,y)F(a(u,v),b(u,v),第三节 频域变换:二维傅立叶变换特性,旋转特性结论:对图象的旋转变换和傅立叶变换的顺序是可交换的FRf(x,y)RFf(x,y),第三节 频域变换:二维傅立叶变换特性,线性与相似性线性的描述:傅立叶变换是线性系统、函数和的傅立叶变换是可分离的设:f(x,y)的傅立叶变换为Ff(x,y)g(x,y)的傅立叶变换为Fg(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字图像 变换
链接地址:https://www.31ppt.com/p-6165365.html