遥感图像频域增强处理.doc

成都信息工程学院 1 实验名称：遥感图像频域增强处理 2 实验目的 对图像进行傅立叶变换，使用低通、高通及其它滤波器，并尝试改变滤波器窗口大小，看滤波后的图像差异，从而了解图像频率域增强的处理方法。 3实验原理： 傅里叶变换是指非周期函数的正弦或余弦和乘以加权函数的积分表示，数字图像处理中所用的傅里叶变换均属于FFT，傅里叶变换分为连续傅里叶变换和离散傅里叶变换，在数字图像处理中经常用到的是二维离散傅里叶变换。 4数据来源： 数据时通过国际科学数据服务平台下载的攀枝花地区的TM图像经过裁剪处理后，分辨率为30米，坐标系统为WGS—84 各波段波长与分辨率 波段 波长/nm 像素大小/m B1 450~520 30 B2 520~600 30 B3 630~690 30 B4 760~900 30 B5 1550~1750 30 B6 10400~12500 120 B7 2080~2530 30 5实验步骤 5.1快速傅里叶变化 5.1.1打开一个图像如下图 5.1.2在主菜单中，选择Filter—FFT Filtering—Forward FFT，在Forward FFT Input File对话框中，选择输入图像文件。 选择文件后点击OK 选择存储位置后点OK 转换过程如下图 5.1.3经过快速傅里叶变换（FFT）之后生成的图像为： 从图上看，中间很亮的部分集中了图像的低频信息；外围较暗的部分集中了图像的高频信息。 5.2 FFT滤波器的选择 5.2.1在Display窗口中显示一幅FFT图像。 5.2.2在主菜单中，选择Filters—FFT Filtering—Filter Definition。在Filter Definition选择对话框中，选择当前显示的FFT图像的Display窗口，单击OK。 5.2.3在Filter Definition对话框中，选择Filter _Type—滤波器类型。选择不同的滤波器。 Circular Pass为低通滤波器 Circular Cut为高通滤波器 Band Pass和Band Cut滤波器 5.3反向FFT变换 5.3.1在主菜单中，选择Filters—FFT Filtering—Inverse FFT，选择houfft文件，点击OK。 5.3.2出现如下对话框，选择刚才保存的滤波器类型，点OK 5.3.3选择输出路径，点OK 完成傅里叶变换。 6实验结果与分析 下图分别为低通滤波、高通滤波、Band Pass滤波的处理图像。 傅里叶变换能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。通过傅里叶变换可以将图像从空间域转到频率域进行操作处理，不同的滤波器有不同的效果，要根据实际情况选择合适的滤波器才能得到想要的结果。 7实验心得 总的来说，实验做的过程很坎坷，首先是选图的问题，我首次选的图像比较大，处理速度慢是一个问题，而且处理过程中还会出错；再有在进行滤波半径选择的时候选的都是100，所以实验效果不太明显。但是实验的整体思路还是了解的。 通过这次实验我明白了，实验数据的选择的重要性。在今后中会多加注意这方面的问题。