matlab图像处理各种程序代码要点.pdfVIP

图像变换 (傅立叶变换 ), 图像增强 , 边缘检测 , 滤波 , 图像压缩等 . 实验工具 :MA TLAB 软件 课程设计时间 :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( 傅立叶变换后的频谱图 ); %然后对上述二维快速傅立叶变换提高分辨率： 要提高二维快速傅立叶变换的分辨率，在采样率一定的情况下，增大采 样点数 N 即可。对应的频谱图见测试结果。 N=300; f=zeros(100,100); f(30:70,45:55)=1; figure(1) imshow(f,notruesize); title( 原始图像 ); F=fft2(f,N,N); axis on figure(2) imshow(F,[-1,5],notruesize); title( 二维快速傅立叶变换后的图像 ); x=1:N;y=1:N; figure(3) mesh(abs(F)); title( 傅立叶变换后的频谱图 ); N=400; f=zeros(100,100); f(30:70,45:55)=1; figure(1) imshow(f,notruesize); title( 原始图像 ); F=fft2(f,N,N); axis on figure(2) imshow(F,[-1,5],notruesize); title( 二维快速傅立叶变换后的图像 ); x=1:N