fft变换水印技术.ppt

第五讲 FFT变换技术 一、函数的傅里叶变换 傅里叶变换是光电信息处理中应用最广泛的一种变换。如整个模拟光学系统分析、变换和处理，数字图像增强、图像复原、图像编码、图像分割、图像分析与描述等，每一类处理方法和手段都用到傅里叶变换。例如：在进行图像的低通滤波、高通滤波中，可以借助傅里叶变换，把要求在空域中解决的问题转换到频率域中去解决。 1、一维连续函数的傅里叶变换 令f(x) 为实变量x的连续图像函数，f(x)的傅里叶变换定义如下： 　　　　Ｆ(u)=∫f(x)exp(-j2πux)dx (1) 积分上下限 为(-∞,+∞)． u为频率变量，Ｆ(u)为f(x)的傅里叶变换 若已知Ｆ(u)，则利用傅里叶反变换可以由下式求得： 　　　　f(x)＝ ∫Ｆ(u)exp(j2πux)du (2) (1)(2)被称为傅里叶变换对． 可以证明，如果f(x)是连续可积的，而且Ｆ(u)也是可积的，可以证明此变换对存在．实际上，这些条件几乎总是可以满足的 1、一维连续函数的傅里叶变换 用Ｆ(.)和F-1(.)表示傅里叶变换和傅里叶反变换．通常f(x)是实函数，Ｆ(u)是复数． Ｆ(u)=R(u)+jI(u) R(u),I(u)分别为Ｆ(u)的实部和虚部，也可表示成 Ｆ(u)=|Ｆ(u)|ejΦ(u) |Ｆ(u)|=[R2(u)+I2(u)]1/2 Φ(u)=arctan(I(u)/R(u)) |Ｆ(u)|称为f(x)的傅里叶谱，而Φ(u)为相位谱，傅里叶谱的平方称为f(x)的能量． ２、二维连续函数的傅里叶变换 F(u,v)=∫∫f(x,y)exp{-j2π(ux+vy)}dxdy f(x,y)= ∫∫F(u,v)exp{j2π(ux+vy)}dxdy F(u,v)是f(x,y)的二维傅里叶变换。 F(u,v)同样可以表示成傅里叶谱（幅值谱），和相位谱。也可以表示成为实部R(u,v)和虚部I(u,v)。u,v 通常称为频率变量。可以这样理解其物理意义：在傅里叶变换中exp{-j2π(ux+vy)}和exp{j2π(ux+vy)}称为变换的基函数。可以用欧拉公式将傅里叶变换的基函数表示成如下形式： exp{-j2π(ux+vy)}＝cos(2 π(ux+vy))-j sin(2 π(ux+vy)) ２、二维连续函数的傅里叶变换 exp{-j2π(ux+vy)}＝cos(2 π(ux+vy))-j sin(2 π(ux+vy)) ３、离散傅里叶变换 对一个连续函数f(x)等间隔采样，得到一个离散序列，设共采集了N个样本数据，则这个离散序列可以表示为｛f(0),f(1),…..f(N-1)}。令x为离散实变量，u为离散频率变量，可将离散傅里叶变换对定义为： 　 N-1 F(u)=(1/N)∑f(x)exp{-j2πux/N} u=0,1…N-1 x=0 N-1 f(x)= ∑F(u)exp{j2πux/N} x=0,1…N-1 u=0 同学们自己证明：将f(x)代入F(u)中证明相等。 …………….。 ４、函数傅里叶变换的性质 １、共轭性 　　f*(x,y) ４、函数傅里叶变换的性质 3、线性性 　　af１(x,y)+bf2(x,y) ４、函数傅里叶变换的性质 ５、平移性 f(x,y)exp(j2π(u0x+v0y)/N) ４、函数傅里叶变换的性质 6、连续变换性 　FT{FT{f(x,y)}}=f(-x,-y) ４、函数傅里叶变换的性质 8、旋转性 　若将原始图像旋转θ角，则其傅里叶变换也旋转θ角 f(r, θ+ θ0) ４、函数傅里叶变换的性质 10、导数定理 FT{g(m,n)(x,y)}=(j2πu)m(j2πv)nG(u,v) 5、图像傅里叶变换举例 5、图像傅里叶变换水印 6、光学傅里叶变换技术 单透镜傅里叶变换系统 双透镜傅里叶变换系统 双透镜傅里叶变换系统 * * x y 1/v 1/u 1/(u2+v2)1/2 为空间周期 相位差２π的等位线 N-1 N r=u ∑exp(j2πrx/N)exp(-j2πux/N)={ X=0