matlab图像处理各种程序代码.doc

图像变换(傅立叶变换), 图像增强, 边缘检测, 滤波, 图像压缩等. 实验工具: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(傅立叶变换后的频谱图); %然后对上述二维快速傅立叶变换提高分辨率： ??? 要提高二维快速傅立叶变换的分辨率，在采样率一定的情况下，增大采 样点数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》