数字图像处理04.pptVIP

数字图像处理;

本章主要内容;

?空间域是由f(x,y)所张成的坐标系.

?频域系统是由F(u,v)所张成的坐标系.

?傅里叶逆变换可表示为;

在频域原点处变换的值称为傅里叶变换的直流(DC)分量

频域原点

数学上,F(0,0)

MATLAB坐标体系为F(1,1)

对于图像,常定为频率矩形中心

傅里叶变换

空间域(x,y)频域(u,v)傅里叶逆变换;

4.2在MATLAB中计算并可视化二维DFT

?格式1：F=fft2(f)

?格式2：F=fft2(f,P,Q)P,Q的含义？

快速傅里叶逆变换(函数ifft2)

?格式：f=ifft2(F)

将变换的原点移动到频率矩形中心(函数fftshift)

?格式：Fc=fftshift(F)

将频率矩形中心移动到变换的原点(函数ifftshift)

?格式：F=ifftshift(Fc)

■上述函数使用演示：图4.3（f43.m） ;4.3频域滤波

4.3.1基本概念

?空间域和频域线性滤波的基础都是卷积定理。

f(x,y)*h(x,y)H(u,v)F(u,v)f(x,y)h(x,y)H(u,v)*F(u,v)

?滤波传递函数H(u,v)

?低通滤波器—衰减高频分量而保持低频分量不变;

4.3频域滤波;4.3频域滤波

4.3.2DFT滤波的基本步骤

1.使用函数paddedsize获得填充参数

PQ=paddedsize(size(f))

2.得到使用填充的傅里叶变换

F=fft2(f,PQ(1),PQ(2))

3.生成一个滤波函数

4.将变换乘以滤波函数:G=H.*F

5.获得G逆傅里叶变换的实部

g=real(ifft2(G))

6.将g左上部的矩形修剪为原始大小

g=g(1:size(f,1),1:siae(f,2));

4.3频域滤波

4.3.2DFT滤波的基本步骤;

4.3频域滤波

4.3.3用于频域滤波的M函数

Functiong=dftfilt(f,H)

F=fft2(f,size(H,1),size(H,2))

G=H.*F

g=real(ifft(G))

g=g(1:size(f,1),1:size(f,2));4.4从空间滤波器获得频域滤波器

IPT提供了生成空间滤波器的函数（3.5）

fspecial

对给定的空间滤波器h，对应的频域滤波器H为:

H=fft2(h,PQ(1),PQ(2))

函数freqz2:计算FIR滤波器的频率响应

H=freqz2(h,R,C)

■演示;4.5在频域中直接生成滤波器

4.5.1建立用于实现频域滤波器的网格数组

函数dftuv

函数meshgrid;

4.5在频域中直接生成滤波器

4.5.2低通频域滤波器

理想低通滤波器（ILPF）;

4.5在频域中直接生成滤波器

4.5.3线框图与表面图

函数mesh

函数gridon/off

函数axison/off

函数view;4.6锐化频域滤波器

4.6.1基本的高通滤波器

Hhp(u,v)=1-Hlp(u,v);4.6.2高频强调滤波器

高通滤波器偏离了直流项,从而把图像的平均值降低到了零。

Hhfe(u,v)=a+bHhp(u,v);

●

17