电子信息工程专业综合课程设计报告.docVIP

电子信息工程专业综合课程设计报告        课程设计内容: 图像变换(傅立叶变换), 图像增强, 边缘检测, 滤波, 图像压缩等. 实验工具: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(傅立叶变换后的频谱图); 　　%DC系数移动（使用函数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 %绘制立体网状图，将图形对象的色度改为灰度