实验报告六.docxVIP

南京信息工程大学数字图像处理实验报告实验课程 数字图像处理 实验名称 图像频域增强实验日期 2015-5-22 指导老师 徐旦华专业 信息与计算科学 年级 2012 姓名 杨鑫 学号 20121314059 得分 实验目的1．掌握怎样利用傅立叶变换进行频域滤波2．掌握频域滤波的概念及方法3．熟练掌握频域空间的各类滤波器4．利用MATLAB程序进行频域滤波实验内容利用MATLAB提供的低通滤波器实现图像信号的滤波运算，并与空间滤波进行比较。命令文件clcf=imread(moon.tif); subplot(2,4,1);imshow(f);title(原图像); fftf=fft2(f);%二维离散傅里叶变换sfftf=fftshift(fftf); RR=real(sfftf);ff=imag(sfftf);S=sqrt(RR.^2+ff.^2);S=(S-min(min(S)))/(max(max(S))-min(min(S)))*255;%归一化subplot(2,4,2);imshow(S);title(频谱图像) ;%显示频谱图像h=fspecial(sobel);%产生空间‘sobel’模板%freqz2(h) %查看相应频域滤波器的图像PQ= size(f);%产生滤波时所需大小的矩阵H=freqz2(h,PQ(1),PQ(2)); %产生频域中的‘sobel’滤波器H1=ifftshift(H);%重排数据序列，使得原点位于频率矩阵的左上角subplot(2,4,5);imshow(abs(H),[]) title(图形滤波) ; subplot(2,4,6);imshow(abs(H1),[])title(重排图形滤波) ; gs=imfilter(double(f),h); %用模板h进行空域滤波gf=dftfilt(f,H1);%用滤波器对图像进行频域滤波subplot(2,4,3);imshow(gs,[]) title(模板h滤波) ; subplot(2,4,7);imshow(gf,[])title(滤波器滤波) ; subplot(2,4,4);imshow(abs(gs),[])title(模板绝对滤波) ; subplot(2,4,8);imshow(abs(gf),[])title(滤波器绝对滤波) ; 函数文件function [ gs] = dftfilt( f,H1)f=fft2(f);gs=f*H1; end图像结果：结论：图像信号的滤波运算和空间滤波进行比较，其中空间滤波比较直观，可以直接显示出滤波后的图像，信号滤波效果较好，但不直观，需要进行分析后得出结论。利用MATLAB提供的高通滤波器对图像进行处理。clcI = imread(moon.tif);subplot(1,2,1),imshow(I)title(原始图像) h=fspecial(laplacian); I1=filter2(h,I);subplot(1,2,2),imshow(I1)title(laplacian高通滤波器滤波) 结论：对图像采用laplacian高通滤波器进行锐化滤波，由于高通滤波器会允许高频部分通过，低频部分被阻，图像边缘为高频部分，因此图像边界轮廓线十分明显。