图 像 变 换 (傅 立 叶 变 换 ), 图 像 增 强 , 边 缘 检 测 , 滤 波 , 图 像 压 缩 等
实 验 工 具 :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, 为 白 色 , 其 他 区