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

图像变换(傅立叶变换), 图像增强, 边缘检测, 滤波, 图像压缩等. 实验工具:MATLAB 软件 课程设计时间:2008 年 12 月实 验 部 分 图像变换程序代码及说明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;y=1:N; figure(3) mesh(abs(F)); title(傅立叶变换后的频谱图); 系数移动（使用函数 fftshift） N=200; 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(fftshift(abs(F))); %把零频率部分移到频谱的中间title(傅立叶变换后的频谱图); %滤波器频率响应x=1:N;y=1:N; mesh(x,y,F2(x,y));colormap(gray);colorbar %绘制立体网状图，将图形对象的色度改为灰度图像，colorbar 给坐标轴添加色彩条测试结果图： 1》对 f 进行二维快速傅立叶变换 （请自己运行查看） 2》对上述二维快速傅立叶变换提高分辨率N=300 时： N=400 时 ： N=800 时： N=1000 时： （请自己运行查看） 从傅立叶变换的频谱图中可以看出，提高分辨率以后，其边缘更加平滑，锯齿状明显减弱。但其傅立叶变换后的图像没有明显改变。 3》DC 系数移