实验要求:

Objective:

To further understand the well-known algorithm Fast Fourier Transform (FFT) and verify its effectiveness to calculate the discrete Fourier transform (DFT).

Main requirements:

Ability of programming with C, C++, or Matlab.

Instruction manual:

The purpose of this project is to develop a 2-D FFT program “package” that will be used in several other projects that follow. Your implementation must have the capabilities to:

(a) Multiply the input image by (-1)x+y to center the transform for filtering.

(b) Multiply the resulting (complex) array by a real function (in the sense that the the real coefficients multiply both the real and imaginary parts of the transforms). Recall that multiplication of two images is done on pairs of corresponding elements.

(c) Compute the inverse Fourier transform.

(d) Multiply the result by (-1)x+y and take the real part.

(e) Compute the spectrum.

Basically, this project implements Fig. 4.5. If you are using MATLAB, then your Fourier transform program will not be limited to images whose size are integer powers of 2. If you are implementing the program yourself, then the FFT routine you are using may be limited to integer powers of 2. In this case, you may need to zoom or shrink an image to the proper size by using the program you developed in Project 02-04.

An approximation: To simplify this and the following projects (with the exception of Project 04-05), you may ignore image padding (Section 4.6.3). Although your results will not be strictly correct, significant simplifications will be gained not only in image sizes, but also in the need for cropping the final result. The principles will not be affected by this approximation.

实验要求我们编写程序实现2维图像的快速傅里叶变换,并运用于频率域滤波。

简要介绍一下频率域滤波的步骤:

  1. 给定一幅大小为M*N的输入图像f(x,y),计算得到填充参数P和Q。通常,取P=2M,Q=2N;
  2. 添加必要数量的0,形成大小为P*Q的填充后的图像fp(x,y);
  3. 用(-1)^(x+y)乘以fp(x,y),将图像移到其傅里叶变换的中心处;
  4. 计算(3)中的图像DFT,得到F(u,v);
  5. 生成一个实的、对称的滤波器函数H(u,v),其大小与fp(x,y)相同,为P*Q,中心在(P/2, Q/2)。用阵列相乘形成乘积G(u,v);
  6. 计算(5)中结果的逆DFT,取出实部,得到gp(x,y);
  7. 从gp(x,y)的左上象限提取M*N的图像,得到最终处理结果g(x,y)。

上程序:

%%
close all;
clc;
clear all; %%
% 读取图像
img = imread('gray_image.jpg');
figure(1)
imshow(img);
title('original A'); % 得到填充参数P和Q
[M, N] = size(img);
P = 2 * M;
Q = 2 * N; img = double(img);
% 添加必要数量的0
img_fp = zeros(P, Q);
img_fp(1:M, 1:N) = img(1:M, 1:N); figure(2);
imshow(img_fp, []);
title('image B'); % 用(-1)^(x+y)乘以图像的结果
for x = 1:P
for y = 1:Q
img_fp(x, y) = img_fp(x, y) .* (-1)^(x+y);
end
end % figure(3);
% imshow(img_fp, []);
% title('image C'); % 对图像做快速傅里叶变换
img_Fp = fft2(img_fp); % figure(4);
% imshow(img_Fp, []);
% title('image D'); % H = ones(P, Q);
% H(P/2, Q/2) = 0; % H = zeros(P, Q);
% H(P/2, Q/2) = 1;
% H(P/2 - 1, Q/2 - 1) = 1; r = 30;
H = ones(P, Q);
for x = 1:P
for y = 1:Q
d = sqrt((x-M)^2 + (y-N)^2);
if d > r
H(x, y) = 0;
else
H(x, y) = 1;
end % if x == P/2 && y == Q/2
% H(x, y) = 1;
% else
% H(x, y) = 1;
% end
end
end figure(5);
imshow(H, []);
title('image E'); img_G = img_Fp .* H; % figure(6);
% imshow(img_G, []);
% title('image F'); img_g = ifft2(img_G);
img_g = real(img_g); for x = 1:P
for y = 1:Q
img_g(x, y) = img_g(x, y) .* (-1)^(x+y);
end
end % figure(7);
% imshow(img_g, []);
% title('image G'); img_o = img_g(1:M, 1:N); figure(8);
imshow(img_o, []);
title('result H'); imwrite(img_o, 'result.jpg');

实验结果:



上图是原图像;



上图是步骤1和2,扩展后的图像fp(x,y);



上图是滤波器函数H(u,v),中心一个圆内为1,圆外都是0,表示一个低通滤波器;



上图是逆DFT计算得到的gp(x,y);



上图为最终的结果,从gp(x,y)提取出的g(x,y)图像;

明显得,可以发现,低通滤波器使原图像变得模糊了。

数字图像处理实验(5):PROJECT 04-01 [Multiple Uses],Two-Dimensional Fast Fourier Transform 标签: 图像处理MATLAB数字图像处理的更多相关文章

  1. 数字图像处理实验(总计23个)汇总 标签: 图像处理MATLAB 2017-05-31 10:30 175人阅读 评论(0)

    以下这些实验中的代码全部是我自己编写调试通过的,到此,最后进行一下汇总. 数字图像处理实验(1):PROJECT 02-01, Image Printing Program Based on Half ...

  2. 数字图像处理实验(12):PROJECT 05-03,Periodic Noise Reduction Using a Notch Filter 标签: 图像处理MATLAB 2017-0

    实验要求: Objective: To understand the principle of the notch filter and its periodic noise reducing abi ...

  3. 移植u-boot.2012.04.01

    /*************************************************** *u-boot版本:u-boot2012.04.01 *gcc版本:arm-linux-gcc ...

  4. MATLAB数字图像处理(二)图像增强

    1         图像增强 1.1            直方图均衡化 对于灰度图像,可以使用直方图均衡化的方法使得原图像的灰度直方图修正为均匀的直方图. 代码如下: I2=histeq(I1); ...

  5. ok6410 u-boot-2012.04.01移植六完善MLC NAND支持

    继ok6410 u-boot-2012.04.01移植四.五后,开发板基本已支持MLC NAND,支持DM9000.但是通过NAND命令更新u-boot到NAND,还存在问题,需要根据u-boot的n ...

  6. ok6410 u-boot-2012.04.01移植七完善u-boot移植(u-boot移植结束)

    继ok6410 u-boot-2012.04.01移植六后,开发板已支持MLC NAND.DM9000等.但还需要完善比如环境变量.mtdpart分区.裁剪.制作补丁等.下面的工作就是完善移植的u-b ...

  7. ok6410 u-boot-2012.04.01移植五支持DM9000

    继ok6410 u-boot-2012.04.01移植四后,开发板基本已支持MLC NAND,但还有一些细节地方修改,这节增加DM9000支持,通过网卡tftp程序到内存,接着通过NAND命令写到NA ...

  8. ok6410 u-boot-2012.04.01移植二修改源码支持单板

    继ok6410 u-boot-2012.04.01移植一后修改代码,对ok6410单板初始化,主要包括时钟.串口.NAND.DDR等初始化.这些工作在以前的裸板程序都写了,直接拿来用.我觉得先写裸板程 ...

  9. 移植u-boot-2012.04.01到JZ2440

    开发环境:Ubuntu 12.04 开发板:JZ2440  256M NandFlash  64M SDRAM 交叉编译器:arm-linux-gcc-4.3.2 u-boot:u-boot-2012 ...

随机推荐

  1. WPF 变量转换的实现

    有时候,我们传入的是一个值,但是真正显示的需要是另一个值,这时候就需要转换.比如我们传入一个枚举值,而不同的枚举值对于的图片是不一样的. 这时候就需要一个转换规则.WPF里面给我们提供了一个接口IVa ...

  2. Codeforces554B:Ohana Cleans Up

    B. Ohana Cleans Up Time Limit: 2000ms Memory Limit: 262144KB 64-bit integer IO format: %I64d      Ja ...

  3. 自己封装的json工具类

    package com.develop.util; import java.util.ArrayList; import java.util.HashMap; import java.util.Ite ...

  4. Visual Studio 2015里面汇编工具Asm Dude的配置!

    最近开始学习汇编,也开始使用Visual Studio 2015写汇编程序,比较了半天,最后觉得Asm Dude应该是目前Visual Studio 2015非常好的汇编开发插件了,但是如果默认安装上 ...

  5. hadoop集群默认配置和常用配置【转】

    转自http://www.cnblogs.com/ggjucheng/archive/2012/04/17/2454590.html 获取默认配置 配置hadoop,主要是配置core-site.xm ...

  6. iText 文本

    iText中用文本块(Chunk).短语(Phrase)和段落(paragraph)处理文本. 文本块(Chunk)是处理文本的最小单位,有一串带格式(包括字体.颜色.大小)的字符串组成.如以下代码就 ...

  7. Installshield自动安装IIS组件

    原文:Installshield自动安装IIS组件 一. 工程类型:IS2010 installscript 二.原理: 1. XP和 Server2003系统:由于系统默认没有自带IIS组件,一般情 ...

  8. linux jdk环境变量

    export JAVA_HOME=/usr/share/jdk8 export PATH=$JAVA_HOME/bin:$PATH export CLASSPATH=.:$JAVA_HOME/lib/ ...

  9. vueSSR全栈(项目实战 mac)

    1.准备安装及指定版本 参考安装类中的 安装部分(node,npm,webpack) nuxt 官网下载nuxt脚手架(可以自定义版本) 需要下载MongoDB  redis 以及数据库可视化工具 具 ...

  10. 并发编程(十五)——定时器 ScheduledThreadPoolExecutor 实现原理与源码深度解析

    在上一篇线程池的文章<并发编程(十一)—— Java 线程池 实现原理与源码深度解析(一)>中从ThreadPoolExecutor源码分析了其运行机制.限于篇幅,留下了Scheduled ...