第9章图像增强扩展阅读.pdfVIP

图像增强补充案例 （配第9 章）  测试数据 数据目录：附带光盘下的..\chapter09\data\ 文件名 说明 Landsat8_OLI_b 1.dat 某地Landsat 8 OLI 图像的某波段，灰度值被拉伸到0 ~ 100 Landsat8_OLI_b2.dat 某地Landsat 8 OLI 图像的某波段，灰度值被拉伸到0 ~ 255 Landsat8_OLI_b2.dat 图像被裁剪后数据，行列数为 394 ×394 ，用 Landsat8_OLI_FFT_test.dat 于傅里叶变换  IDL 程序 程序目录：附带光盘下的..\chapter09\programs\ 文件名 说明 histomatch.sav IDL 代码，用于直方图匹配，具体参见电子补充材料 9.0 启动ENVI+IDL 单击开始 所有程序 ENVI 5.2 ENVI+IDL ，稍等片刻，计算机将依次 启动IDL 8.4.1 和ENVI 5.2. 1 软件。 9.1 直方图匹配 直方图匹配就是将原图像的直方图以参考图像的直方图为标准作变换，使两 幅图像的直方图相同或近似，从而使原图像具有与参考图像类似的色调和方差。 由于ENVI5.2.1 的窗口界面还没有直方图匹配功能，本书提供封装好的IDL 程序 histomatch.sav 可实现该功能。本实验以Landsat8_OLI_b1.dat 图像为原始图像， 以Landsat8_OLI_b2.dat 图像为参考图像。IDL 程序实现操作步骤如下： （1）打开图像。在菜单栏中，单击File Open ，选择Landsat8_OLI_b1.dat 和Landsat8_OLI_b2.dat 文件，则成功加载了两幅图像。 （2 ）将ENVI 软件中的图像导出为IDL 中的变量。在工具箱中，选择Raster Management IDL Export to IDL Variable ，在Export to IDL Input File 中选中 1 Landsat8_OLI_b1.dat 文件。单击OK 。然后，在Export Variable Name 对话框中， New Variable Name 文本框输入变量名为“data_o ”（原始图像）。然后，按照上述 方法把Landsat8_OLI_b2.dat 文件导入到IDL 中，变量名为“data_r ”（参考图像）。 （3 ）导入直方图匹配函数。在IDL 控制台中输入如下代码，注意需根据实 际情况替换histomatch.sav文件的路径。 RESTORE,D:\RSDIP\chapter09\programs\histomatch.sav;导入直方图匹配 函数 result=HISTOMATCH(data_o,data_r);data_o为原始图像，data_r为参考图像 4 IDL ENVI ENVI Raster （）将 中的变量导入为 中的图像。在 工具箱中，选择 Management IDL Import from IDL Variable ，在弹出的Import IDL Variables 对 result ENVI 9.1 话框中选择变量 导入到 中保存的变量。匹配结果如图 所示，左 上角图像为原图，右上角图像为参考图像，左下角图像为匹配后的图像。在工具 栏点击按钮 ，可查看图像直方图，效果如图 9.2 所示，可以看出直方图匹配