【数字图像处理】Two-Dimensional Fourier Transform and Filtering

实验要求：

该实验的目的是开发一个 2-D FFT 程序包。要求程序能完成下面的功能:
(1.a) 用因子 (-1)x+y 乘以输入图像，以实现滤波的中心化变换；
(1.b) 计算付立叶变换；
(1.c) 用一个实矩阵乘以一个复数矩阵，即用实矩阵中的元素同时乘以复数矩阵对应位置上的复数的实部与虚部。
可以通过调用两个图像的乘法程序来实现对应元素的相乘；
(1.d) 计算反付立叶变换；
(1.e) 结果乘以 (-1)x+y，并取其实部；
(1.f) 计算频谱。
(1.g) 计算图 4.41(a)的中心化后的付立叶谱；
(1.f) 使用二维图像和三维图形两种形式显示计算的幅度谱和相位谱；
该实验主要实现 4.7.3 节的算法。若试验中使用 MATLAB 语言实现, 则在(1.b)和(1.d)中二维付立叶变换
和反变换程序可以使用 Matlab 的函数 FFT2 和 IFFT2。
(2.a) 编写一个实现公式(4.8-6)所示的高斯低通滤波器的程序。
(2.b) 对图 4.41(a)进行高斯低通滤波以获得图 4.48(c)。
(3.a) 如公式(4.9-1)所示，原始图像减去实验(2.b)中的图像得到锐化后的图像。观察结果图像与图 4.56 中得
到的高斯高通滤波结果图的差异， 并解释其中的原因。
(3.b) 调整高斯低通滤波器变量，使(3.a)中的减法锐化图与图 4.56(c)接近为止，并解释你的结果。

技术论述：