图 像 变 换 (傅 立 叶 变 换 ), 图 像 增 强 , 边 缘 检 测 , 滤 波 , 图 像 压 缩 等 . 实 验 工 具 :MATLAB 软 件 课 程 设 计 时 间 :2008 年 12 月 实 验 部 分 1. 图 像 变 换 程 序 代 码 及 说 明 clear all N=100; f=zeros(50,50); %产 生 一 个 50*50 的 全 0 数 组 f(15:35,23:28)=1; %定 义 图 像 数 组 ， 从 15 行 到 35 行 ， 23 列 到 28 列 附 值 为 1， 为 白 色 ， 其 他 区 域 为 黑 色 figure(1) %创 建 窗 口 的 图 形 对 象 ， 句 柄 为 1 imshow(f,'notruesize') %显 示 图 像 f F=fft2(f,N,N); %在 二 维 傅 立 叶 变 换 前 把 f 截 断 或 者 添 加 0， 使 其 成 为 N*N 的 数 组 F2=fftshift(abs(F)); %把 傅 立 叶 变 换 的 零 频 率 部 分 移 到 频 谱 的 中 间 figure(2) x=1:N;y=1:N; %定 义 x 和 y 的 范 围 mesh(x,y,F2(x,y));colormap(gray);colorbar %绘 制 立 体 网 状 图 ， 将 图 形 对 象 的 色 度 改 为 灰 度 图 像 ， colorbar 给 坐 标 轴 添 加 色 彩 条 %构 建 一 个 类 似 于 figure(1)的 矩 形 函 数 N=200; f=zeros(100,100); f(30:70,45:55)=1; %定 义 图 像 数 组 ， 从 30 行 到 70 行 ， 45 列 到 55 列 附 值 为 1， 为 白 色 ， 其 他 区 域 为 黑 色 imshow(f,'notruesize'); %然 后 对 f 进 行 二 维 快 速 傅 立 叶 变 换 ： 以 下 列 出 你自己编写的 代 码 … N=200; f=zeros(100,100); f(30:70,45:55)=1; imshow(f,'notruesize'); title('原始图 像 '); F=fft2(f,N,N); %对 图 像 f 进 行 二 维 快 速 傅 立 叶 变 换 grid on %打开网 格线 axis on %打开坐 标 轴 imshow(F,[-1,5],'notruesize'); %显 示 傅 立 叶 变 换 后 的 图 像 ， 图 像 数 据的 值 域 为 [-1,5] x=1:N;y=1:N; title('二 维 快 速 傅 立 叶 变 换 后 的 图 像 '); mesh(abs(F)); %绘 制 F 的 频 谱 图 title('傅 立 叶 变 换 后 的 频 谱 图 '); %然 后 对 上述二 维 快 速 傅 立 叶 变 ...
