遥感数字图像处理-matlab-频域处理.docVIP

北京师范大学 地理学与遥感科学学院 蓝建航 作者邮箱： jianhanglan@ 欢迎志同道合都联系 目录 TOC \o 1-4 \h \u 27401 实验五：频域图像增强 2 26354 一、 实验目的 2 3263 二、 实验内容 2 3628 1. 傅立叶变换 2 13081 ①实验代码 2 939 ②实验图像 3 25365 ③算法分析 3 25571 2. 图像填充 4 3338 ①实验代码 4 13522 ②实验图像 5 8688 ③算法分析 6 674 3. 频域中滤波 6 10266 ①实验代码 6 25087 ②实验图像 6 29336 ③算法分析 7 3687 4. 滤波比较 7 30939 ①实验代码 8 19399 ②实验图像 8 30808 ③算法分析 10 11800 5. 低通滤波 10 6188 ①实验代码 10 16098 ②实验图像 11 18204 ③算法分析 14 30664 6. 高通滤波 15 3688 ①实验代码 15 653 ②实验图像 17 10303 ③算法分析 19 29239 三、 实验总结 20  实验五：频域图像增强 2024-蓝建航-地理科学 实验目的 1、 理解傅立叶变换的基本原理，掌握频域的频谱，相位和能量谱等基本概念； 2、 掌握傅立叶变换的主要性质，掌握傅立叶变换进行图像填充的必要性和方法，对比填充 前后进行图像傅立叶变换的差别（折叠误差） ； 3、 掌握不同滤波器的特征和频域进行图像滤波的基本步骤； 4、 掌握用 MATLAB 实现低通滤波，高通滤波和高频增强滤波，比较和分析滤波结果（模 糊、纹理和振铃等）； 实验内容 傅立叶变换 对一幅图像进行傅立叶变换，求其频谱和相位谱，并通过傅立叶反变换重构图像。本例 中主要熟悉 MATLAB 的几个傅立叶变换函数: fft2, fftshift, ifft2, ifftshift ①实验代码 { close all; f=imread(C:\实验五\DIPLab5_Data\pea_red.jpg); figure,imshow(f,[]); F=fft2(double(f)); % 傅立叶变换到频域 S=log(1+abs(F)); % 对输出的傅立叶频谱的原点不调整 figure,imshow(S,[]); Fc=fftshift(F); % 对输出的傅立叶频谱的原点进行调整(平移) Sc=log(1+abs(Fc)); % 求傅立叶谱并进行对数拉伸变换 figure,imshow(Sc,[]); ff=real(ifft2(F)); % 傅立叶反变换重构图像,取实部 figure,imshow(ff,[]); % figure,imshow(angle(F),[]) %求相位谱 figure, imshow(gscale(uint8(abs(double(f)-double(ff))))); % 变换前后两图像的差异 } ②实验图像 原孔雀图像 傅立叶频谱的原点不调整 傅立叶频谱的原点进行调整 傅立叶反变换重构图像,取实部 相位谱 变换前后两图像的差异 ③算法分析 图像通过傅立叶变换后的图像原点不在中心，须通过原点调整，其实部不能直接看不出与原图的关系。观察可知反变换后的图像与原图像无差别（相减为0）， 2. 图像填充 掌握折叠误差产生的原因及处理方法，以矩形函数图像为例，对比图像填充前后，进行傅立叶变换的差异。 ①实验代码 { close all; %没填充的滤波: f=imread(C:\实验五\DIPLab5_Data\square.tif); [M,N]=size(f); imshow(f,[]) impixelinfo; F=fft2(double(f)); Fc= log(1+abs(fftshift(F))); figure, imshow(Fc); sig=5; H=lpfilter(gaussian,M,N,sig); G=H.*F; g=real(ifft2(G)); g=gscale(g); figure,imshow(g,[]); %填充后的频域滤波 PQ=paddedsize(size(f)); Fp=fft2(double(f),PQ(1),PQ(2)); %计算 FFT Hp=lpfilter(