实验4.DFT及频域滤波资料.docxVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 实验名称：数字图像DFT及频域滤波 实验目的 了解数字图像各种正交变换的概念、原理和用途。 深熟练掌握数字图像的 DFT/DCT 的原理、方法和实现流程，熟悉两种变换的性质，并能对数字图像 DFT 及 DCT 的结果进行必要解释。 熟悉熟悉和掌握利用 MATLAB 工具进行数字图像 FFT 及 DCT 的基本步骤、MATLAB 函数使用及具体变换的处理流程，并能根据需要进行必要的频谱分析和可视化显示。 熟悉利用空域滤波器构建对应的频域滤波器的方法和关键步骤。 熟悉和掌握几种典型的频域低通滤波器及高通滤波器的原理、特性和作用。 搞清空域图像处理与频域图像处理的异同，包括处理流程、各自的优势等。掌握 频域滤波的基本原理和基本流程，并能编写出相应的程序代码。 实验原理 FFT实验 傅立叶变换：非周期函数表示为正弦和/或余弦乘以加权函数的积分。 一维连续Fourier变换 对函数f(x)进行傅立叶变换得到F(u) MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 逆变换，即将F(u)变换到f(x)为 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 一维离散Fourier变换 正变换(DFT) MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 逆变换(IDFT) MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) 用幅值和相位表示傅立叶变换 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5) 用途：频域滤波，图像增强，频谱分析。 实验步骤 模型图像的 FFT 实验 图 SEQ 图 \* ARABIC 1 图 SEQ 图 \* ARABIC 1 （a） （b） 对变换结果做频谱中心化处理，并分别显示出其 2D 频谱图以及对应的 3D频谱图。 对以上两幅原始图像 FFT 后的频谱图进行分析，可以得出什么样的结论或验证了DFT的什么性质。 图 SEQ 图 \* ARABIC 2 实际图像的 FFT 实验 读任意读取一幅灰度图像，对其进行 FFT 变换，分别画出变换前的原始图像及其 FFT 后频谱的 2D 平面图，要求变换结果做频谱中心化处理。 画出对应的3D中心化频谱。 图 SEQ 图 \* ARABIC 3 数字图像的频域滤波处理 设定截止频率 D0 =100，试分别构建 256×256 的频域理想低通滤波器（ILPF）和频域理想高通滤波器（IHPF）。编写 MATLAB 代码，分别画出它们的频域滤波器响应 3D 图及其对应的 2D 投影平面图。 任意读取一幅数字图像。编写 MATLAB 代码，分别利用理想低通滤波器（ILPF） 和理想高通滤波器（IHPF）对其进行频域滤波处理。具体可参考如下步骤： 消除折叠现象的填充：P=2M，Q=2N，其中 M、N 为原图像的尺寸； 原图 f (x, y) 的傅立叶变换，产生 F(u,v) ； 频谱中心化处理（当然也可以不变换，则 H(u,v) 要改变）； 用滤波器函数 H(u,v) 乘以 F(u,v) ； 傅立叶反变换； 取实数部分，绝对值很小的虚数部分是浮点运算存在误差造成的； 空间域中心还原变换（反中心化） 。 分别显示 ILPF 和 IHPF 两种滤波器滤波前的原始图像、滤波器频域响应（2D）及滤波结果图像。 对滤波结果进行必要分析，如振铃现象、图像模糊、图像变暗等