[数学]傅里叶、小波变换.ppt

4.4 图像傅里叶变换 傅里叶变换是数学上，特别是工程数学上常用 的变换方法。Matlab中的二维快速傅里叶变换函 数是fft2，该函数对应的逆傅里叶变换函数是 ifft2。 图像傅里叶变换函数 在这一节中，还是通过Matlab中的傅里叶变换函数直观上理解分析傅里叶变换。 图像傅里叶变换的数学描述 例4-21 A0= imread(D:\0.bmp); A=double(A0); s=size(A); M=s(1); N=s(2); for u=0:M-1 for x=0:M-1 P(u+1,x+1)= (1/M^0.5)*exp(-j*2*pi*(u*x/M)); end end for v=0:N-1 for y=0:N-1 Q(v+1,y+1)= (1/N^0.5)*exp(-j*2*pi*(v*y/N)); end end C=fft2(A); F=P*A*Q; subplot(1,3,1); imshow(A) subplot(1,3,2); imshow(C) subplot(1,3,3); imshow(F) 4.5 图像小波变换 Matlab中有很多关于小波变换的函数，其中二维小波变换函数都可以用在图像处理、图像压缩等领域。另外Matlab也提供了许多小波函数专门用于图像处理。 二维小波变换 1.二维小波变换 母小波函数是指满足式 的函数，二维母小波函数的构造可由一维母小波函数的张量积形成，也可以采用非张量积的方法构造。 Matlab图像小波变换函数 假定函数 为二维母小波函数。 令 （4-11） 由（4-11）式定义的，产生于二维母小波函数，依赖于参数a、b、c的一组函数就是小波基函数。 二维连续小波变换定义为： （4-12） 二维离散小波变换定义为： （4-13） 令式（4-13）中的，即得到常用的一种离散小波变换： （4-14） 【注】母小波函数与小波基函数都是小波函数。 2.小波函数 Haar小波就是一个常用的小波函数。 （1）Haar小波 Haar小波是最早使用的、简单的紧支撑小波，Haar小波实际上是Daubechies小波族中的一个特例。 （4-15） 利用（4-15）可以构造出二维小波函数。 （2）Daubechies小波 Daubechies小波也有相应的数学描述。其中在Matlab中称为’db1’ 的Daubechies小波与Haar小波是相同的。下面不研究其数学描述，只通过图形显示一组二维Daubechies小波函数。 * * * * 1.图像傅里叶变换函数fft2 【例4-16】利用傅里叶变换函数变换图像，观察分析变换结果。 A= imread(D:\0371.bmp); B= imread(D:\00.bmp); A1=fft2(A); B1=fft2(B); subplot(1,4,1); imshow(A) subplot(1,4,2); imshow(A1) subplot(1,4,3); imshow(B) subplot(1,4,4); imshow(B1) 设计左面程序进行傅里叶变换。 程序的运行结果是图4-12所示。 (a) (b) (c) (d) (a) 原图像0371，(b) 图像0371变换结果，(c) 图像00，(b) 图像00变换结果 图4-12 图像傅里叶变换结果 由于图4-12中的原图像简单，所以能够看出变换后结果的一些特征。如果把上面程序中的原图像换为图4-13中的(a)与(c)，那么得到的变换数据显示出来，视觉上觉得杂乱无章。 (a) (b) (c) (d) (a) 原图像0010，(b) 图像0010变换结果，(c) 图像0041，(b) 图像0041变换结果 图4-13 复杂图像傅里叶变换结果 2.图像逆傅里叶变换函数ifft2 【例4-17】利用傅里叶变换函数变换图像，然后用逆傅里叶变换函数复原图像，观察分析比较。 设计如下程序： A= imread(D:\00